{"id":"396d7477-c552-4965-b8ce-8404d6ab0cac","arxiv_id":"2411.16373","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Inside shrimp-shaped regions of a discrete-time slow-fast membrane model, the paper reports a complete devil's staircase of periodic attractors, with early afterdepolarizations mostly periodic and delayed afterdepolarizations chaotic.","lead":"The paper studies a simple three-variable map model of nerve and heart cell electrical activity and finds that, inside a special shrimp-shaped region of its parameter space, many periodic patterns are arranged like a staircase, with infinitely many steps. The result links abstract slow-fast dynamics to early and delayed afterdepolarizations connected to cardiac arrhythmias, so it could guide how period changes in heart cell models are interpreted.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'complete devil's staircase' claim is not established: finite-time period detection and Df≈1 cannot distinguish complete from harmless staircases, and the <ISI>-1 shift questions the winding-number labels.","rationale":"As a second-pass reviewer, I read the paper as a numerical study of a slow-fast map claiming a new phenomenon: shrimp-shaped regions containing a complete devil's staircase of periodic attractors, with EAD-like attractors mostly periodic and DAD-like mostly chaotic. The phase-diagram observations are credible and the figures support striped internal structure. My concern targets the strongest quantitative claim, completeness of the staircase. The evidence for completeness is a finite number of Farey-labeled steps plus a box-counting dimension near 1, neither of which rules out a harmless staircase with finitely many steps; the dimension of a finite-jump graph is also 1. The paper asserts that gaps are finite-time artifacts but provides no convergence study as the detection horizon grows. There is also an unexplained <ISI>-1 shift in Fig. 7, which either means the labels depend on an ad hoc correction or the paper has not specified how Q and P are measured independently. A convergence test with longer runs and direct period detection would settle both issues. This does not change the reader's conditional verdict; it sharpens the reason for conditionality.","tokens_in":19005,"tokens_out":7077,"duration_ms":69696,"concrete_test":"At fixed T=0.2343864, take a gap between two labeled steps in Fig. 7A and increase the detection horizon N from 2×10^5 to 2×10^6 and 2×10^7, identifying periodic orbits by exact recurrence of x(t) (not by ISI statistics). Record the winding numbers P/Q of all resolved steps and the total xR measure classified as periodic. If the number of resolved steps does not grow (or the periodic measure does not approach the full interval), the 'complete' claim fails. Also recompute the Fig. 7 comparison using raw <ISI> instead of <ISI>-1; if the Farey-labeled 1/w curves no longer align, the winding-number labels are not independently validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"To make the central claim true, two things must hold: (i) the winding numbers P/Q labeling the staircase steps are correct, and (ii) the staircase is actually complete, i.e. steps are dense and the quasiperiodic complement has measure zero. Both are weakly supported. For (i), Eq. (6) gives <ISI>=Q/P, yet Fig. 7 plots <ISI>-1 to match 1/w with no derivation of the -1 offset. If the shift is not a harmless average of the ±1 timestamp fluctuations (in the cosine example of Sec. III A the fluctuations do average to zero), the P/Q labels are not independently confirmed by the ISI data. For (ii), the Farey mediant is exhibited for only two or three adjacent steps per shrimp, and Df≈0.95–0.98 is compatible with a harmless staircase of finitely many jumps (graph dimension 1). The claim that gaps are solely a finite-time artifact is not demonstrated; some labeled periods Q exceed or approach the 2×10^5 iteration horizon, so longer runs could resolve new steps or change labels. Thus the headline claim of a complete devil's staircase overreaches the numerical evidence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a discrete-time three-variable map-based model of cardiac action potentials, interpreted as a single myocyte. It investigates the transition from cardiac spiking (plateau spikes) to bursting, where early and delayed afterdepolarizations (EADs and DADs) appear. The central claims are: (i) inside shrimp-shaped regions of the (T, xR) parameter plane, multiple periodic attractors coexist and form an internal striped structure that, along a one-parameter cut, constitutes a complete devil's staircase of winding numbers w = P/Q; (ii) EADs are mainly periodic (inside shrimps) while DADs are chaotic (between shrimps); and (iii) near quasiperiodic orbits the maximum negative Lyapunov exponent scales as -L ~ 1/P, where P is the period of the interspike-interval sequence. The paper derives the identity <ISI> = Q/P = 1/w from the definitions of spike timestamps and uses it to label staircase steps, exhibits Farey-mediant relations among a few labeled steps, and estimates fractal dimensions Df ≈ 0.95–0.98 for the staircases.","tokens_in":19217,"tokens_out":3311,"duration_ms":32270,"significance":"If the central claim of a complete devil's staircase inside shrimps is correct, it extends the known phenomenology of shrimp-shaped regions beyond isoperiodic and quasiperiodic shrimps, providing a novel link between slow-fast dynamics, cardiac arrhythmia models, and devil's staircases. The paper has notable strengths: the definitions in Sec. III are precise, the identity in Eq. (6) is derived from first principles (though it is an internal consistency check), the numerical methods (Eckmann-Ruelle Lyapunov exponents, box-counting dimension) are standard, and a data availability link is provided. The biological interpretation is speculative but clearly labeled as such. However, the headline claim of a complete (as opposed to harmless or incomplete) devil's staircase is not supported by the finite numerical evidence presented, and the ad hoc <ISI>-1 correction raises questions about the reliability of the winding-number labels. These issues are load-bearing for the paper's main contribution.","major_comments":[{"comment":"The claim of a 'complete devil's staircase' is not established. Completeness means infinitely many steps, dense in the parameter interval, with a zero-measure quasiperiodic complement. The paper exhibits only two or three Farey-mediant relations per shrimp (e.g., w1 = 6/2016, w2 = 3/1007, w3 = 9/3023 in Fig. 7A), and the box-counting dimension Df ≈ 0.95–0.98 is compatible with a harmless staircase with finitely many jumps, whose graph dimension is also 1. The assertion that 'the gaps in the staircases are due only to the finite simulation time' (Section IV D) is not demonstrated; since some labeled periods Q (e.g., 47524 in Fig. 7C) approach or exceed the 2×10^5 iteration horizon, longer simulations could resolve new steps or change the labels. The authors should either provide a quantitative argument for completeness (e.g., scaling of the number of steps with resolution, or an analytic mechanism such as a monotone circle map with a devil's staircase) or weaken the claim to 'devil's-staircase-like structure' with finitely many resolved steps.","section":"Section IV D, Fig. 7"},{"comment":"The winding-number labels are not independently confirmed because Eq. (6) is an identity derived from the definitions of ISI and w, so plotting <ISI> against 1/w is a consistency check, not a validation. The paper then plots <ISI>−1 instead of <ISI> to match 1/w, stating only that 'we obtained a better match' without deriving the correction. The ±1 timestamp fluctuation in the cosine example in Sec. III A averages to zero, but for the actual map the shift is not shown to be uniform across the staircase; if the correction varies from step to step, the P/Q labels become unreliable and the Farey-mediant construction loses its quantitative anchor. The authors should derive the <ISI>−1 offset from the spike-timestamp conditions or show numerically that the correction is constant over each staircase.","section":"Section IV D, Eq. (6), Fig. 7"},{"comment":"The relation −L ∼ 1/P near quasiperiodic orbits is supported only by Pearson correlations (R = 0.74 to 0.91) between rescaled curves of −L and 1/P, and the paper acknowledges it is not derived from first principles. Because both quantities vary smoothly and systematically along the staircase, a high correlation does not provide strong evidence for an asymptotic power-law relation; moreover, the claim that P diverges at step boundaries is handled numerically by arbitrarily replacing diverging P with 10^6, which is not a test of divergence. The qualitative correlation is plausible, but the scaling claim as stated overreaches the evidence. The authors should either present a more direct test (e.g., collapse of data over multiple resolutions, or a derivation in a limit of the map) or rephrase Eq. (8) as a heuristic observation without the '∼' notation.","section":"Section IV E, Eq. (8)"}],"minor_comments":[{"comment":"There are several typos and inconsistencies: 'analogou' (Section IV D), 'Lyuapnov' in the Fig. 4 caption, 'Ex. (6)' instead of 'Eq. (6)' in Fig. 7, and inconsistent accenting of 'Hénon'/'Henón' in the text and references.","section":"Throughout"},{"comment":"The definition of the Farey mediant w3 = P3/Q3 with mP3 = P1 + P2 and mQ3 = Q1 + Q2 is confusing because it suggests P3 and Q3 are not reduced. Please state explicitly that w3 = (P1+P2)/(Q1+Q2) is reduced by a common factor m, and clarify that m is the greatest common divisor or otherwise specify the reduction.","section":"Eq. (7) and surrounding text"},{"comment":"The terms 'sliding' and 'locked' incommensurate phases are introduced without definition or reference. Since these terms are not standard in the nonlinear dynamics literature cited, please define them or give a precise reference.","section":"Section III B"},{"comment":"The caption refers to a 'harmless staircase going from one shrimp to another', but the text in Section IV D says the <ISI> in between shrimps is 'spurious' due to aperiodic behavior. This inconsistency should be resolved.","section":"Fig. 6G,H caption"},{"comment":"The sentence 'The <ISI> is much easier to measure' should be qualified: Eq. (5) defines <ISI> only for periodic or quasiperiodic orbits, so the staircase analysis applies only to non-chaotic regions. Please state this restriction explicitly where the comparison is made.","section":"Section IV D"},{"comment":"Reference [6] is cited as a Research Square preprint while the paper is being published; if a peer-reviewed version exists, it should be cited instead. Also, reference [41] (Rulkov 2007) is listed as arXiv:0708.1173v1, but the paper appears to be a later published version; please check the citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's central novelty—a complete devil's staircase inside shrimps—depends on extrapolating from finite numerical staircases with only a few labeled steps and a box-counting dimension close to 1. The ad hoc <ISI>-1 adjustment in Fig. 7 is a red flag because it is central to the winding-number labels. These issues are fixable in principle, but they require substantial additional analysis (e.g., systematic resolution scaling, a derivation of the timestamp correction, or a more careful statement of the claims). The manuscript is a good candidate for a journal like Chaos, but in its current form the headline claim outruns the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe paper is worth a look for one clear reason: it shows that shrimp-shaped regions in parameter space need not be isoperiodic. Inside these shrimps the authors find a sequence of stripes, each with a different period, and the sequence organizes like a devil's staircase. That is a genuine addition to the shrimp taxonomy and it is documented carefully with phase diagrams, Lyapunov exponents, and ISI period maps.\n\nThe strongest part is the empirical relation between the maximum Lyapunov exponent and the number of cycles, -L ~ 1/P, near quasiperiodic boundaries. The correlations (R = 0.74-0.91) are convincing as a numerical trend, even if the relation is not derived.\n\nNow the soft spots, in proportion.\n\nThe central claim of a 'complete' devil's staircase - infinitely many steps, dense in parameter, zero-measure quasiperiodic complement - is not established. The authors show a few labeled steps and a box-counting dimension near 1, but a harmless staircase with finitely many jumps also has graph dimension 1. The Farey mediant is exhibited for only two or three adjacent steps per shrimp, and some of the labeled periods Q approach the simulation horizon, so the staircase might be finite or might continue to merge into unresolved small steps. The paper asserts that 'gaps are due only to the finite simulation time' but does not test that by longer runs or a scaling analysis. This is the load-bearing claim, and it remains a conjecture.\n\nSecond, the winding-number labels are less secure than they look. Equation (6) gives w = P/Q = 1/<ISI>, but the paper then plots <ISI>-1 to match 1/w, citing intrinsic ±1 fluctuations. The cosine example in Sec. III A shows the fluctuations average to zero, so the -1 shift is not explained by that mechanism. If the shift is not understood, the P/Q labels - and hence the nonstandard Farey construction - are not independently confirmed by the ISI data. This is a quantitative anchor that needs a derivation.\n\nThe cardiac framing is speculative: parameters are arbitrary units, no experimental data, and the connection to arrhythmias is indirect. That said, the authors are careful to present it as prediction, so the biological part is a minor issue.\n\nThe paper is a legitimate numerical study with a novel observation. It deserves a serious referee, but the referee should push for a more honest statement of what is proven: a staircase-like structure with many steps, not a proven complete devil's staircase. I would send it to peer review, with the expectation of revisions.\n\nCheers.","headline":"A solid numerical observation of striped shrimps, but the 'complete devil's staircase' claim outruns the evidence.","tokens_in":19842,"tokens_out":3320,"would_cite":true,"duration_ms":30705,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37D45","37N25","37C27"],"pacs":["05.45.-a","87.19.Hh"],"model":"deepseek-v4-flash","headline":"The paper establishes that shrimp-shaped regions of parameter space in a discrete-time cardiac model contain infinitely many periodic attractors ordered as a complete devil's staircase, with early afterdepolarizations mainly periodic and…","keywords":["devil's staircase","shrimp-shaped regions","interspike interval","afterdepolarization","Lyapunov exponent","slow-fast dynamics","cardiac action potential","Farey sequence"],"falsifier":"Iterate the map long enough to find the exact period $Q$ of an attractor inside a shrimp and compare the Farey label $w = P/Q$ from the ISI sequence; if the labels do not match the actual rotation numbers, the staircase identification fails.","tokens_in":18767,"feed_emoji":"🫀","tokens_out":10232,"duration_ms":79070,"temperature":0.7,"pith_summary":"The paper studies a three-variable discrete-time map that models cardiac action potentials with slow-fast dynamics. It claims that the shrimp-shaped regions of its parameter space, traditionally associated with a single periodic orbit, actually contain infinitely many periodic attractors arranged as a complete devil's staircase as one parameter is varied. It also claims that early afterdepolarizations (EADs) are mostly periodic, while delayed afterdepolarizations (DADs) are chaotic, and that the maximum negative Lyapunov exponent scales roughly as the inverse of the interspike-interval period near quasiperiodic orbits. If true, these results give a quantitative handle on how periodicity changes as a cardiac plateau loses stability, which is relevant to arrhythmia mechanisms.","feed_headline":"Inside shrimps, infinite periodicities form a devil's staircase","feed_subtitle":"A cardiac map model shows arrhythmia-linked afterdepolarizations hide a Farey-organized ladder of spike periods.","key_machinery":"The object carrying the argument is the winding number $w = P/Q$, computed from the interspike interval (ISI) sequence of the map. Each ISI is the time between successive upcrossings of $x = 0$, and for a periodic attractor of period $Q$ with $P$ spikes per period, the average ISI equals $Q/P$, so $w = 1/\\langle\\text{ISI}\\rangle$. Because time is discrete, the timestamps fluctuate by $\\pm 1$, and the paper compensates by plotting $\\langle\\text{ISI}\\rangle - 1$ against $1/w$ to match the staircase steps. The maximum period $P$ of the ISI sequence, the maximum Lyapunov exponent $L$ (computed by the Eckmann\\textendash Ruelle method), and the box-counting fractal dimension $D_f$ of the staircase in $(x_R, w)$ space are the supporting measurements that identify the staircase as complete and Farey-organized.","core_discovery":"The central discovery is that shrimp-shaped regions in the $(T, x_R)$ parameter plane of the map are not isoperiodic: each shrimp is filled with stripes of constant interspike-interval period $P$, and traversing $x_R$ at fixed $T$ crosses these stripes as steps of a complete devil's staircase. The steps are labeled by rational winding numbers $w = P/Q$, where $Q$ is the period of the attractor and $P$ is the number of cycles (spikes) per period, and the labels follow a nonstandard Farey tree so that between any two steps another step exists. Quasiperiodic orbits, with maximum Lyapunov exponent $L = 0$ and diverging $P$, form the zero-measure boundaries between steps. Along the staircase the relation $-L \\sim 1/P$ holds, with the plateaus of $1/P$ coinciding with peaks of $-L$. EAD attractors are found inside shrimps and are mostly periodic, while DAD attractors occupy chaotic regions between shrimps.","pith_inferences":["One could test the staircase labeling in a continuous-time cardiac model: if the rotation number computed from the actual voltage waveform differs systematically from $w = P/Q$, the discrete-time timestamp correction is not universal.","The relation $-L \\sim 1/P$ may hold for complete devil's staircases in other dissipative maps, not only this model; checking it in a different system with a known staircase would indicate whether the scaling is generic.","The fractal dimension near 0.95\\textendash 0.98 suggests the staircase is nearly space-filling; this could imply that in a real myocyte, noise would wash out the finest steps, so only the largest steps would be physiologically observable.","If chaotic DADs are indeed more arrhythmogenic than periodic EADs, then the position of a parameter region within a shrimp (periodic interior vs. chaotic border) may correlate with clinical risk; this is not tested here."],"forward_implications":["Shrimp-shaped regions in slow-fast systems can host infinitely many periodic attractors, not just one, so isoperiodicity is not a defining feature of shrimps.","The complete devil's staircase implies that along the chosen parameter axis the quasiperiodic orbits form a set of zero measure, with periodic steps dense, so the transition from cardiac spiking to bursting passes through a dense ladder of periodicities.","The scaling $-L \\sim 1/P$ offers a practical proxy: measuring the ISI period near a quasiperiodic boundary estimates the Lyapunov exponent without computing Jacobians.","The model predicts that shifting the slow-current reversal potential $x_R$ toward EAD behavior produces a Farey-organized sequence of period changes, a signature that could be looked for in biophysically detailed cardiac myocyte models.","Since EADs are mostly periodic and DADs are chaotic, the periodicity of the ISI sequence could serve as a distinguishing marker between these two classes of afterdepolarizations."],"supporting_citations":[{"why":"Defines the complete devil's staircase and its commensurate–incommensurate transition properties, the framework used to identify the staircase inside shrimps.","marker":"7"},{"why":"Introduced shrimp-shaped regions in parameter space as isoperiodic attractors, the notion this paper overturns by showing internal periodicity structure.","marker":"4"},{"why":"Recently found quasiperiodic orbits inside shrimp-shaped regions in a predator-prey system, the result this paper extends by finding periodic staircases instead.","marker":"3"},{"why":"Provides the map-based neuron model with the logistic transfer function $F(u)=u/(1+|u|)$ and the phase-diagram classification of oscillation modes used here.","marker":"15"},{"why":"Supplies the Eckmann–Ruelle method for computing Lyapunov exponents, used to measure $L$ and establish the $-L \\sim 1/P$ scaling.","marker":"29"},{"why":"Documents nonstandard Farey sequences in a physical map, the pattern used to label the staircase steps inside the shrimps.","marker":"27"},{"why":"Gives the slow-fast analysis of this model, including the delayed Neimark–Sacker bifurcation and the identification of EAD and DAD behaviors.","marker":"6"}],"fun_headline_variants":["Shrimp-shaped regions hide a devil's staircase of spike periods","Cardiac spikes: infinite period steps inside shrimp regions","Devil's staircase found in shrimp-shaped parameter zones","Afterdepolarizations reveal Farey steps in cardiac map","Shrimp interiors: complete devil's staircase of periodicities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The winding number $w = P/Q$, measured from the discrete-time upcrossing timestamps of the ISI sequence, is assumed to correctly label the rotation number of each periodic orbit, relying on a uniform $\\pm 1$ timestamp fluctuation correction so that the Farey labels are accurate.","fun_headline_variants_meta":{"raw":{"variants":["Shrimp-shaped regions hide a devil's staircase of spike periods","Cardiac spikes: infinite period steps inside shrimp regions","Devil's staircase found in shrimp-shaped parameter zones","Afterdepolarizations reveal Farey steps in cardiac map","Shrimp interiors: complete devil's staircase of periodicities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000172,"raw_usage":{"total_tokens":1293,"prompt_tokens":979,"completion_tokens":314,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":248}},"tokens_in":595,"tokens_out":314,"duration_ms":3847,"temperature":1.0,"reasoning_tokens":248,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:12:12.079064+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Iterate the map long enough to find the exact period $Q$ of an attractor inside a shrimp and compare the Farey label $w = P/Q$ from the ISI sequence; if the labels do not match the actual rotation numbers, the staircase identification fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Recently found quasiperiodic orbits inside shrimp-shaped regions in a predator-prey system, the result this paper extends by finding periodic staircases instead."}],"review_version":1}