Back to search
2104.07304

Understanding Calcium Dynamics via Process-Oriented Geometric Singular Perturbation Theory

Samuel Jelbart, Nathan Pages, Vivien Kirk, James Sneyd, Martin Wechselberger

correctmedium confidence
Category
Not specified
Journal tier
Strong Field
Processed
Sep 28, 2025, 12:56 AM

Audit review

The paper states and proves Theorem 5.9: for system (5.1) there is, for small ε, a locally unique relaxation cycle Γε that is O(ε^{1/3})-close to the singular orbit Γ = Γc ∪ Γh ∪ Γl and is exponentially attracting with Floquet exponent bounded above by −κ/ε^2. The statement and proof outline appear verbatim in the main text and Appendix A, including the multi-regime analysis (R1/R2), Krupa–Szmolyan fold blow-up giving the O(ε^{1/3}) closeness, and a Poincaré map with contractions of order e^{−a/ε^2} along the infra-slow manifold, which yields the −κ/ε^2 Floquet bound (Theorem 5.9 and Lemma 5.8; see also the explicit e^{−a/ε^2} estimates in Appendix A) . By contrast, the candidate solution reproduces the existence argument and the O(ε^{1/3}) proximity near the fold, but explicitly leaves the sharp Floquet bound −κ/ε^2 “open,” only obtaining a weaker e^{−κ1/ε} estimate; it also assumes (rather than proves) the existence/uniqueness of the fast ‘drop’ Γl and relies on heuristic monotone-wedge arguments near the degenerate axes instead of the paper’s cylindrical-plus-spherical blow-ups and foliations. Since the sharp Floquet rate and the rigorous construction across the degenerate line Sh are central parts of the theorem as stated and proved in the paper, the model solution is incomplete relative to the problem, while the paper’s analysis is correct and complete on these points .

Referee report (LaTeX)

\textbf{Recommendation:} minor revisions

\textbf{Journal Tier:} strong field

\textbf{Justification:}

This is a careful application of beyond-standard-form GSPT and blow-up methods to a biologically motivated calcium model exhibiting nonuniform time-scale separation. The existence of a three-timescale relaxation cycle, the refined O(ε\^{1/3}) closeness near the fold, and the sharp −κ/ε\^2 Floquet bound are convincingly established via a global Poincaré map built on cylindrical and successive spherical blow-ups. The analysis is technically sound and contributes a transferable blueprint for similar problems, though some Appendix material could be summarized more in the main text.