Pith. sign in

REVIEW 3 major objections 5 minor 41 references

Microscopic Variability Alters Macroscopic Rotation Speed in Stochastic Spiral Waves

T0 review · 3 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper derives a general second-order formula for how noise alters the rotation speed of spiral waves, splitting the correction into an always-slowing instantaneous part and an orbital-drift part that can go either way.

desk verdict Genuinely new second-order decomposition of noise-induced spiral wave slowing, with a quantitative validation for multiplicative temporal noise on the Barkley model; additive and attenuating noise matches are weaker than the abstract suggests. read the letter →

arxiv 2511.21710 v3 pith:WOKDW5Q6 submitted 2025-11-18 nlin.PS math-phmath.MPphysics.bio-ph

Jolien Kamphuis , Desmond Kabus , Hermen Jan Hupkes , Tim De Coster This is my paper · ORCID
classification nlin.PSmath-phmath.MPphysics.bio-ph MSC 35K5760H1537L10
keywords spiralwavesstochasticreaction-diffusionequationsangularvelocitynoise-inducedslowdownorbitaldriftBarkleymodelexcitablemediacenter-manifoldreduction
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to establish that the average rotation speed of a stochastic spiral wave is the deterministic speed plus a correction of order σ² (the noise variance), and that this correction splits into two physically distinct parts: an instantaneous term coming from how noise distorts the rotating-wave profile, which always slows rotation, and an orbital-drift term coming from fluctuations around that profile, which can either slow or speed rotation. For the Barkley excitable-medium model with multiplicative temporal noise, the theory predicts Ωσ = 1.82 − 1.00σ² and computer simulations give 1.82 − 0.97σ², a close quantitative match. If the formula holds generally, it means microscopic noise is a generic mechanism for macroscopic slowing of spiral rotation in any excitable medium — a point the authors connect to cardiac arrhythmias and neural waves.

What carries the argument

The argument is carried by three objects: (i) the exponentially localized adjoint eigenfunctions ψ0, ψ± of the deterministic spiral's linearization, which project perturbations onto the rotational and translational symmetry directions; (ii) the decomposition of the noisy solution into an instantaneous stochastic rotating wave (Φσ, ωσ) plus a perturbation φ, giving the split of the correction; and (iii) three bilinear maps B1, B2, Bod — traces of products of the noise covariance operator with these projections — that turn the abstract SPDE into numbers, so that ω^(2) = TrB1[h, R] + TrB2[R, R] and ω_od^(2) = ∫₀^∞ Tr Bod[e^{Ls}S0, e^{Ls}S0] ds. The Itô lemma and the mild Itô formula are the too

What would settle it

Recompute ω^(2) and ω_od^(2) on a 384×200 polar grid with dt=0.0005 and check that the totals in Eq. (20) (and the additive/attenuating totals) shift by less than the O(σ³) error; alternatively, measure Ωσ for additive temporal noise at σ=0.06–0.12 with 64 runs and compare the fitted slope against the predicted −11.08 and the reported −8.01.

Watch

Extended reading notes

Core claim

The central discovery is equation (17): the total average angular velocity of the stochastic spiral is Ωσ = ω0 + σ²[ω^(2) + ω_od^(2)] + O(σ³). Here ω^(2) is the second-order shift of the 'instantaneous stochastic wave' (Φσ, ωσ) — the rotating solution that would be seen if all fluctuations were frozen at their mean — and it always has the sign of slowing. The second term ω_od^(2) is the orbital-drift correction, the long-time averaged contribution of the fluctuation field φ perturbing that instantaneous wave; its sign depends on the noise type and on which part of the wave (front or back) the noise acts. The paper shows the odd-order terms vanish by symmetry and evaluates both coefficients u

Load-bearing premise

The whole quantitative match rests on the numerical evaluation of the spectral coefficients in Appendix C — a 192×100 polar grid for the adjoint eigenfunctions, trapezium integration of the semigroup with dt=0.001, and truncated Q-traces — being accurate; the paper provides no convergence or discretization-error analysis, and the reported fits for additive and attenuating temporal noise deviate noticeably from the theoretical totals.

Editorial extensions

If this is right

  • If equation (17) is right, noise-induced slowing of spiral rotation is a generic second-order effect that holds for any noise type and any reaction-diffusion model, not a quirk of the Barkley equations.
  • Because the orbital-drift term can have either sign, parameter regimes should exist where net rotation accelerates, a reversal the authors note but do not exhibit.
  • Extrapolating measured Ωσ linearly in σ² back to σ=0 recovers the noise-free angular velocity ω0 to about 98% accuracy — a direct recipe for estimating intrinsic rotor speed from noisy cardiac or neural data.
  • The framework separates noise effects on the wavefront from those on the waveback, offering a mechanistic way to predict which ion-channel-like perturbations will speed or slow rotors.
  • Slower spirals are less prone to breakup but more prone to anchoring, so the predicted slowdown has opposite clinical consequences in fibrillation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Extension the paper leaves implicit: the same σ² decomposition should survive for three-dimensional scroll waves, where the translation modes are three-dimensional; the rotational symmetry argument carries over directly.
  • Testable extension: because additive temporal noise acts on both wavefront and waveback while multiplicative noise acts mainly on the waveback, the sign of ω_od^(2) should track the width of the waveback — a prediction one could check by varying the Barkley parameter b.
  • The discrepancy between fitted coefficients for additive (-8.01) and attenuating (-2.13) temporal noise and the theoretical totals (-11.08 and -5.28) suggests the Appendix C spectral computation may not be converged for those noise types; a grid-refinement study is the natural next check.
  • The near-zero slowdown for spatio-temporal noise raises a possibly general principle: temporal averaging over fast fluctuations partially cancels the instantaneous term, meaning the noise autocorrelation time is itself a control parameter for rotation speed.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops a perturbation theory for the average angular velocity of stochastic spiral waves. Starting from an SPDE model and assuming a deterministic rotating spiral with angular velocity ω0, the authors derive Eq. (17), which expresses the noise-corrected angular velocity as Ωσ = ω0 + σ²[ω^(2) + ω_od^(2)] + O(σ³), splitting the second-order correction into an 'instantaneous' term ω^(2) and an 'orbital-drift' term ω_od^(2). The coefficients are computed for the Barkley model with additive, multiplicative, and attenuating temporal noise, and compared with finite-difference simulations. The paper reports that multiplicative noise gives a theoretical coefficient of -1.00 versus a simulated -0.97, and it further reports simulation results showing that spatial and spatio-temporal noise also slow the spiral. The discussion connects these findings to cardiac arrhythmias and neural wave propagation.

Significance. The conceptual decomposition of stochastic frequency shifts into instantaneous and orbital-drift contributions is appealing and, if rigorously validated, would provide a general framework for predicting noise-induced rotation slowdown in excitable media. The paper's strengths include an independent numerical computation of the theoretical coefficients, open data/code, and a clear statement of the O(σ³) limitation. However, the quantitative validation is currently demonstrated for only one of three temporal noise types, and the numerical evaluation of the coefficients in Appendix C lacks a convergence analysis. These gaps prevent the central quantitative claim from being fully supported as stated.

major comments (3)
  1. [§II.B, §II.C, Eqs. (20)-(21)] The abstract states that 'Analytical predictions match Barkley-model simulations with temporal noise', but the only quantitative match shown is for multiplicative temporal noise: theory gives -1.00 and the simulation fit gives -0.97. For additive and attenuating temporal noise, the theoretical total coefficients are -11.08 and -5.28, while the fitted values reported in Section II.C are -8.01 and -2.13, i.e. deviations of 28% and 60%. The text says these noise types 'agree closely ... when we consider small values of σ', but no small-σ fits or their uncertainties are reported. To support the stated claim, the authors should report the restricted-σ² fits with confidence intervals, or modify the claim to apply only to multiplicative temporal noise.
  2. [Appendix C] The numerical evaluation of ω^(2) and ω_od^(2) uses a 192×100 polar grid, a trapezium time-stepping with dt=0.001, and a truncated semigroup/Q-trace, but no convergence, grid-refinement, or discretization-error study is provided. Because the match for multiplicative noise and the mismatches for additive/attenuating noise depend on these computed coefficients, the accuracy of this numerical evaluation is load-bearing. Please add convergence checks in grid size and time step, and report the truncation error of the semigroup integral in Eq. (16).
  3. [Appendix A] Appendix A says it 'summarizes the computations in [38]'. However, Eqs. (7), (15), and (16) are the central theoretical results of the paper, and the derivation from the Itô lemma and orthogonality conditions to the explicit bilinear-trace expressions is not self-contained. This makes it difficult for a reader to verify the theory without accessing the thesis [38]. Please provide the full derivation in the paper or in a supplementary file, or explicitly state which theorem/equation in [38] corresponds to each step.
minor comments (5)
  1. [§II.B] Typo: 'cerate' should be 'create' in the sentence listing the noise types.
  2. [Fig. 1 and §II.B] The text refers to 'translations (Fig. 1D)' and 'rotations (Fig. 1E)', but the Fig. 1 caption identifies panel D as the rotation angle and does not mention panel E. Please align the figure panels and text references.
  3. [§II.C] The phrase 'besides temporal for the term ξ' is awkward; it should be 'besides temporal noise for the term ξ'.
  4. [Appendix A] Typo: 'The indexialways ranges' should read 'The index i always ranges'.
  5. [Eq. (21)] Equation (21) is presented as a fitted relation, but it is written with the same form as the theoretical expansion (20). Clarify that (21) is an empirical fit over the sampled σ range, not a truncated Taylor expansion verified to all orders.

Circularity Check

0 steps flagged · score 1.0 of 10

No material circularity: the theoretical coefficient in Eq. (17) is computed from the deterministic spiral linearization and noise model, while the Barkley simulations are an independent benchmark; the only notable self-citation is Appendix A's reliance on the first author's thesis [38].

full rationale

The derivation chain for Eq. (17) is: define the stochastic decomposition (2); derive coordinate/perturbation SDEs (3)-(4); expand the instantaneous stochastic wave (5)-(6) to obtain omega^(2) via Eq. (7); expand the orbital drift (9)-(11), solve the linearized mild equation (13)-(14), and obtain omega_od^(2) via Eq. (16); then combine into Eq. (17). The coefficients are functionals of the deterministic spiral Phi0, its adjoint eigenfunctions, and the noise function h, evaluated independently in Appendix C without using the simulated Omega_sigma values. The Barkley validation (Appendix D) is a separate finite-difference GPU simulation with a different discretization and domain radius; its fitted coefficient -0.97 is not inserted into the theory. The only self-citation of note is Appendix A, which summarizes the derivation of the coordinate SDEs from the first author's thesis [38]; this is prior work rather than the target prediction, and the formulas are stated in the paper, so the central result is not a self-citation loop. The numerical evaluation in Appendix C lacks convergence analysis, and the full-range fitted coefficients for additive/attenuating temporal noise deviate from the theoretical totals, but these are accuracy/falsifiability concerns, not circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central prediction (Eq 17) requires no fitted parameters: the coefficients ω^(2) and ω_od^(2) are computed from the deterministic spiral and its adjoint eigenfunctions, then compared with direct simulations. The assumptions are the standard existence/localization properties of spiral waves and the small-noise perturbation expansion.

assumptions (5)
  • domain assumption Existence of a deterministic rotating spiral wave Φ0 with angular velocity ω0 for the Barkley model at σ=0
    Assumed throughout; standard for Barkley model and supported by simulation (Section II.B).
  • domain assumption Adjoint eigenfunctions ψ0, ψ± are exponentially localized, so low-dimensional center-manifold projection describes perturbation dynamics
    Invoked in Section II.A.a; cited to [15].
  • ad hoc to paper For small σ, the stochastic solution admits a power series expansion in σ and odd terms vanish by symmetry
    Used to obtain Eqs (5)-(6), (10)-(11); no rigorous small-noise theorem is proven here.
  • ad hoc to paper Finiteness and truncation of the Q-trace and semigroup integral in Eq (16) to a numerical grid
    Appendix C computes coefficients on a 192×100 polar grid with dt=0.001; no convergence analysis.
  • standard math Itô lemma and mild Itô formula apply to the infinite-dimensional SPDE
    Used to derive Eqs (3)-(4) and (16).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Microscopic Variability Alters Macroscopic Rotation Speed in Stochastic Spiral Waves." pith.science (2026). https://pith.science/paper/WOKDW5Q6

@misc{pith2026251121710,
  author       = {Pith},
  title        = {Pith review of: Microscopic Variability Alters Macroscopic Rotation Speed in Stochastic Spiral Waves},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WOKDW5Q6}},
  note         = {Machine review of arXiv:2511.21710}
}
read the original abstract

We present a general theory for noise-induced corrections to the angular velocity of spiral waves. Stochasticity produces two second-order effects: an instantaneous term from heterogeneity that always slows rotation, and an orbital-drift term from temporal fluctuations that can either accelerate or decelerate it. For our parameters, orbital drift is weaker, producing a net slowdown. Analytical predictions match Barkley-model simulations with temporal noise. Examination of additional noise types in silico confirms angular velocity slowing. This mechanism provides a robust route by which stochasticity reshapes spiral dynamics in excitable media, with direct implications for arrhythmias and neural wave propagation.

Figures

Figures reproduced from arXiv: 2511.21710 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 3 linked inside Pith

  1. [38]

    Kamphuis,Stochastic effects on the angular velocity and translational displacements of spiral waves, Master’s thesis, Leiden University (2022)

    J. Kamphuis,Stochastic effects on the angular velocity and translational displacements of spiral waves, Master’s thesis, Leiden University (2022)

    2022

  2. [1]

    J. D. Murray,Mathematical biology: I. An introduction, Vol. 17 (Springer Science & Business Media, 2007)

    2007

  3. [2]

    A. V. Panfilov, inSpirals and Vortices: In Culture, Na- ture, and Science(Springer, 2019) pp. 209–215

    2019

  4. [3]

    A. T. Winfree, Science175, 634 (1972)

    1972

  5. [4]

    L. S. Tsimring, Reports on Progress in Physics77, 026601 (2014)

    2014

  6. [5]

    Alonso and F

    S. Alonso and F. Sagu´ es, Physical Review E63, 046205 (2001)

    2001

  7. [6]

    Song, Journal of Biological Physics49, 521 (2023)

    E. Song, Journal of Biological Physics49, 521 (2023)

    2023

  8. [7]

    Steinbock, V

    O. Steinbock, V. Zykov, and S. C. M¨ uller, Nature366, 322 (1993)

    1993

Show all 41 references
  1. [8]

    Kuehn, J

    C. Kuehn, J. MacLaurin, and G. Zucal, Stochastics and Dynamics22, 2240029 (2022)

    2022

  2. [9]

    Hamster and H

    C. Hamster and H. Hupkes, Physica D Nonlinear Phe- nomena401(2020)

    2020

  3. [10]

    van den Bosch and H

    M. van den Bosch and H. J. Hupkes, Studies in Applied Mathematics155, e70114 (2025)

    2025

  4. [11]

    M. v. d. Bosch and H. J. Hupkes, arXiv preprint arXiv:2504.09350 (2025)

    2025 arXiv

  5. [12]

    S. V. Lototsky, B. L. Rozovsky,et al.,Stochastic partial differential equations, Vol. 11 (Springer, 2017)

    2017

  6. [13]

    Karczewska, arXiv preprint math/0511512 (2005)

    A. Karczewska, arXiv preprint math/0511512 (2005)

    2005 arXiv

  7. [14]

    Da Prato and J

    G. Da Prato and J. Zabczyk,Stochastic equations in in- finite dimensions, Vol. 152 (Cambridge university press, 2014)

    2014

  8. [15]

    Sandstede and A

    B. Sandstede and A. Scheel,Spiral waves: linear and nonlinear theory, Vol. 285 (American Mathematical So- ciety, 2023)

    2023

  9. [16]

    Barkley, M

    D. Barkley, M. Kness, and L. S. Tuckerman, Physical Review A42, 2489 (1990)

    1990

  10. [17]

    Da Prato, A

    G. Da Prato, A. Jentzen, and M. R¨ ockner, Transactions of the American Mathematical Society372, 3755 (2019). 6

    2019

  11. [18]

    M¨ orters and Y

    P. M¨ orters and Y. Peres,Brownian motion, Vol. 30 (Cam- bridge University Press, 2010)

    2010

  12. [19]

    C. H. Hamster and H. J. Hupkes, SIAM Journal on Math- ematical Analysis52, 1386 (2020)

    2020

  13. [20]

    M. v. d. Bosch, C. H. Hamster, and H. J. Hupkes, arXiv preprint arXiv:2502.02164 (2025)

    2025 arXiv

  14. [21]

    Shardlow, Journal of computational and applied mathematics175, 429 (2005)

    T. Shardlow, Journal of computational and applied mathematics175, 429 (2005)

    2005

  15. [22]

    R. A. Gray, A. M. Pertsov, and J. Jalife, Nature392, 75 (1998)

    1998

  16. [23]

    R. A. Gray, J. Jalife, A. V. Panfilov, W. T. Baxter, C. Cabo, J. M. Davidenko, and A. M. Pertsov, Science 270, 1222 (1995)

    1995

  17. [24]

    Ten Tusscher and A

    K. Ten Tusscher and A. Panfilov, Physical Review E68, 062902 (2003)

    2003

  18. [25]

    Z. Qu, J. N. Weiss, and A. Garfinkel, American Jour- nal of Physiology-Heart and Circulatory Physiology276, H269 (1999)

    1999

  19. [26]

    F. H. Fenton, E. M. Cherry, H. M. Hastings, and S. J. Evans, Chaos: An Interdisciplinary Journal of Nonlinear Science12, 852 (2002)

    2002

  20. [27]

    J. M. Davidenko, A. V. Pertsov, R. Salomonsz, W. Bax- ter, and J. Jalife, Nature355, 349 (1992)

    1992

  21. [28]

    Defauw, N

    A. Defauw, N. Vandersickel, P. Dawyndt, and A. V. Pan- filov, American Journal of Physiology-Heart and Circu- latory Physiology307, H1456 (2014)

    2014

  22. [29]

    Takagi, A

    S. Takagi, A. Pumir, D. Pazo, I. Efimov, V. Nikolski, and V. Krinsky, Physical review letters93, 058101 (2004)

    2004

  23. [30]

    Z. Qu, F. Xie, A. Garfinkel, and J. N. Weiss, Annals of biomedical engineering28, 755 (2000)

    2000

  24. [31]

    Dhein and A

    S. Dhein and A. Salameh, Cells10, 2422 (2021)

    2021

  25. [32]

    Panfilov, Physical review letters88, 118101 (2002)

    A. Panfilov, Physical review letters88, 118101 (2002)

    2002

  26. [33]

    Nezlobinsky, O

    T. Nezlobinsky, O. Solovyova, and A. Panfilov, Scientific Reports10, 764 (2020)

    2020

  27. [34]

    De Coster, P

    T. De Coster, P. Claus, I. V. Kazbanov, P. Haemers, R. Willems, K. R. Sipido, and A. V. Panfilov, Scientific reports8, 2050 (2018)

  28. [35]

    Roxin, N

    A. Roxin, N. Brunel, and D. Hansel, Physical review letters94, 238103 (2005)

    2005

  29. [36]

    G. Tse, E. T. Lai, A. P. Lee, B. P. Yan, and S. H. Wong, Frontiers in physiology7, 230 (2016)

    2016

  30. [37]

    Code for

    J. Kamphuis, D. Kabus, H. J. Hupkes, and T. De Coster, “Code for ”microscopic variability alters macroscopic ro- tation speed in stochastic spiral waves”,” (2025)

    2025

  31. [39]

    Kabus, H

    D. Kabus, H. Dierckx, and T. De Coster, Computer Physics Communications (2025), under Review

    2025

  32. [40]

    Kabus, L

    D. Kabus, L. Arno, L. Leenknegt, A. V. Panfilov, and H. Dierckx, PLOS ONE17, 1 (2022)

    2022

  33. [41]

    Tomii, M

    N. Tomii, M. Yamazaki, T. Ashihara, K. Nakazawa, N. Shibata, H. Honjo, and I. Sakuma, Computers in Biology and Medicine130, 104217 (2021). 7 Appendix A: Mathematical definitions Here we summarize the computations in [38] to pro- vide explicit definitions for the functions appe...

    2021

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.