{"id":"764f5fed-3e97-4590-9562-833625517c13","arxiv_id":"2411.11209","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A tutorial deriving and numerically illustrating the relaxation oscillations, Hopf and pitchfork bifurcations, homoclinic orbit, and canard locus of the FitzHugh-Nagumo fast-slow system.","lead":"This paper gives a step-by-step dynamical analysis of the classic FitzHugh-Nagumo neuron model, a fast-slow system with one fast and one slow variable. It computes the period of relaxation oscillations, locates Hopf and pitchfork bifurcations, and applies canard theory to predict where duck-shaped cycles appear.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Example 6.1 applies Theorem 4 with λ=¯c and l5=-1, which forces ¯c_c=+3/8ε; the printed -3/8ε has the wrong sign, so the paper's only quantitative canard-location prediction is internally inconsistent.","rationale":"After checking the singular-flow analysis, the period integral (14), the Hopf computations, and the normal-form setup, the only place where a central quantitative claim breaks down is Example 6.1. The period estimate and the b,c bifurcation diagrams are consistent with known FH-N dynamics; the homoclinic orbit is indeed only numerically observed, and the paper's own wording there is observational ('we have observed numerically'), so the abstract's 'rigorously check all conditions' is an overstatement but not the main quantitative failure. The canard sign error is internal: substituting the stated li values and λ = ¯c into the paper's own Theorem 4 gives +3/8ε. Therefore the printed formula cannot be used as is; the correction is a straightforward sign change if numerics confirm the positive sign, so the paper should be revised rather than rejected. This supports the reader's CONDITIONAL verdict, so no change to the verdict is recommended.","tokens_in":20230,"tokens_out":11566,"duration_ms":109950,"concrete_test":"For ε = 10^{-3} and ε = 10^{-4}, use MatCont (or an independent continuation code) on system (23) to locate the canard-explosion parameter c where a headless canard has an intermediate amplitude, e.g. maximum |x| close to the right fold or orbit length A ≈ 10; compute r_ε = (c_H - c)/ε. If r_ε → +3/8 as ε → 0, the printed -3/8ε is a sign typo and the correct formula is +3/8ε; if r_ε → -3/8, the normal-form identification (the convention for λ/l5) is wrong and Example 6.1 must be rewritten. An immediate algebraic cross-check is to substitute λ = ¯c and l5 = -1 into (27): the paper's own Theorem 4 then gives +3/8ε.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (28) is x' = -y + x^2(-2√3 - x), y' = ε(x + ¯c). In Example 6.1 the authors set l1=1, l2=-2√3-x, l3=0, l4=1, l5=-1, l6=0. Because the normal form (27) has y' = ε(±x l4 - λ l5 + y l6), this substitution gives -λ l5 = +λ, so λ must equal ¯c. Then A = (0 + 3(-1) - 0 + 0)/8 = -3/8 and B = (0 + 0)/2 = 0, so Theorem 4 yields λ_c = -(B + A)ε = +3/8ε, i.e. ¯c_c = +3/8ε + O(ε^{3/2}), not -3/8ε. The sign cannot be rescued by setting λ = -¯c without also changing l5 to +1; the paper states both 'canards occur for ¯c_c' and l5 = -1. The numerics in Fig. 10 are consistent with positive ¯c (c ≈ 1.150077 for ε = 0.5 and c ≈ 1.153794 for ε = 0.1 are both below c_H, so ¯c > 0), although at those ε values the O(ε^{3/2}) correction is too large to make the comparison quantitative. This formula is the paper's central quantitative canard claim and is also used to identify the Hopf bifurcation as subcritical for c. A sign error here does not invalidate the tutorial's qualitative picture, but it makes the headline canard prediction unreliable as printed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper is a tutorial-style study of the two-dimensional FitzHugh-Nagumo fast-slow system (5). It develops the critical-manifold geometry and slow flow, proves the existence of a relaxation-type periodic orbit in the singular limit, estimates its period in Eq. (14), applies Fenichel theory to obtain a first-order slow-manifold expansion, computes Hopf and pitchfork bifurcations in the two parameter cases b=0 and c=0, and uses the Krupa-Szmolyan normal form to identify canard and singular Hopf bifurcations. All theoretical statements are accompanied by numerical simulations. The paper advertises an analytic proof of relaxation oscillations, Hopf-induced periodic solutions, and a rigorous check of the conditions for homoclinic orbits and canards in the FH-N system.","tokens_in":20633,"tokens_out":18725,"duration_ms":163353,"significance":"The paper is a useful pedagogical resource: the critical-manifold analysis, the period integral (14), and the linearized Hopf/pitchfork computations are explicit, parameter-free, and independently checkable. The numerical illustrations are extensive and well matched to the qualitative geometric picture. However, the quantitative canard-location formula in Example 6.1 is affected by a sign error in the normal-form identification, and the claimed equivalence in Example 6.2 is not exact; both issues concern the advertised rigorous check of canard conditions. The underlying geometric singular-perturbation narrative appears sound, and the errors are local and fixable.","major_comments":[{"comment":"The text sets l1=1, l2=-2*sqrt(3)-x, l3=0, l4=1, l5=-1, l6=0 and lambda=cbar. In the normal form (27), y'=epsilon(+/-x l4 - lambda l5 + y l6); substituting l5=-1 gives y'=epsilon(x+lambda), so lambda must equal cbar. Then A=(-0+3(-1)-0+0)/8=-3/8 and B=0, and Theorem 4 gives lambda_c=-(B+A)epsilon=+3/8 epsilon+O(epsilon^{3/2}). The printed value -3/8 epsilon has the wrong sign. It cannot be repaired by setting lambda=-cbar without changing l5 to +1. Moreover, Figure 10 reports canard explosion near c=1.150077 for epsilon=0.5 and c=1.153794 for epsilon=0.1, both below c_H=2/sqrt(3), so cbar>0 there; this is consistent with +3/8 epsilon, not -3/8 epsilon. Since this is the paper's headline quantitative canard prediction, the sign error must be corrected.","section":"Section 6, Example 6.1, Eq. (28) and Theorem 4"},{"comment":"The stated change of variables (xbar,ybar,lambda)=(x-x+, y-8x+/3, b-3/8) does not transform system (24) into the printed system (29). Direct substitution gives xbar'=-ybar-xbar^3-3x+ xbar^2 and ybar'=epsilon[xbar-lambda(ybar+8x+/3)-3/8 ybar]; the term -epsilon lambda(8x+/3) is omitted in (29). The system identified in the normal form is therefore not equivalent to the original system. The omission is consequential: with the printed l6=-3/8, Theorem 4 gives lambda_H=-B epsilon=+3/16 epsilon, whereas the exact Hopf value in Section 5, b_H=(-4+sqrt(16+3epsilon))/epsilon, expands as b_H=3/8-9epsilon/512+O(epsilon^2), so lambda_H=-9epsilon/512+O(epsilon^2). The two asymptotics have opposite signs. The application of Theorem 4 in Example 6.2 needs to be redone with the correct transformed system and an appropriately rescaled parameter.","section":"Section 6, Example 6.2, Eq. (29)"},{"comment":"The statement that a supercritical Hopf bifurcation at E+ 'generates an unstable periodic solution' is internally contradictory under the standard terminology used in the same section for the c-Hopf case. A supercritical Hopf bifurcation produces a stable periodic orbit, while an unstable periodic orbit corresponds to a subcritical bifurcation. The text and Figure 8 locate the unstable periodic solutions for b>b_H, where the equilibrium is stable; this description should be reconciled with the normal-form application in Example 6.2.","section":"Section 5(ii) and Example 6.2"}],"minor_comments":[{"comment":"The value x+=4/3 is incorrect: for b=3/8, x+=sqrt(4-1/b)=sqrt(4-8/3)=2/sqrt(3). The correct value also gives l2=-2sqrt(3)-xbar rather than l2=-4-xbar.","section":"Example 6.2"},{"comment":"There is a recurring typo 'Hodgin-Huxley' where 'Hodgkin-Huxley' is intended; this occurs in the abstract and in Section 1.","section":"Abstract and Introduction"},{"comment":"The caption contains a typo: 'with §epsilon = 0.1' should read 'with epsilon = 0.1'.","section":"Figure 10 caption"},{"comment":"The sign convention for lambda and l5 in the normal form (27) deserves an explicit remark, since the confusion between lambda=cbar and lambda=-cbar is at the root of the sign error in Example 6.1.","section":"Theorem 4 and Example 6.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a tutorial and its main value is pedagogical. The central singular perturbation and bifurcation analysis is mostly sound, but the canard-section errors are load-bearing for the advertised quantitative claims. The sign error in Example 6.1 and the missing constant term in Example 6.2 should be fixed before publication; after that, the tutorial would be a suitable contribution. The self-citation [13] appears only in the future-work discussion and does not raise a circularity concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuinely useful tutorial on GSPT for the FitzHugh-Nagumo system, but the headline quantitative canard claim, the formula in Example 6.1, has the wrong sign, and the paper overstates what it proves about the homoclinic orbit.\n\nWhat is actually new: not the phenomena—relaxation oscillations and canards in FH-N are well established—but the explicit period integral, the first-order slow manifold expansion, the normal-form computation for the singular Hopf/canard point, and a clean bifurcation analysis of both parameter cases. These are worked out carefully and would save a student a lot of time. The singular analysis in Section 3 and the Hopf/pitchfork calculations in Section 5 checked out.\n\nThe soft spot is Example 6.1. The stress-test is right: with l5 = -1 and λ = ¯c, Theorem 4 gives λ_c = -(B + A)ε = +3/8ε, not -3/8ε. The paper's own numerics (c ≈ 1.150 and 1.154, both below cH) support the corrected sign—canards at positive ¯c. So the qualitative picture is fine, but the only quantitative canard-location prediction in the paper is printed with the wrong sign and needs correction. Note also the O(ε^{3/2}) term is large at the simulated ε, so the numerical evidence doesn't independently confirm the coefficient, but the sign definitely agrees with the corrected formula.\n\nSecond, the introduction says all conditions for a homoclinic orbit are 'rigorously checked', but the homoclinic bifurcation in Section 5 is only observed numerically. The local Hopf and pitchfork results are proved, and the canard conditions are checked via Theorem 4, but the homoclinic claim is not a proof. That should be reworded.\n\nThird, minor: the tutorial occasionally falls into 'we leave the proof to references' for Fenichel and Hopf, which is fine, but it should not be described as fully self-contained. The self-citation [13] in the discussion is benign.\n\nBottom line: the paper is a solid pedagogical contribution, not a research advance. With the sign fix and a more accurate claim about the homoclinic orbit, it would be a good fit for a journal that publishes expository/tutorial work. I would send it to a referee and expect a revise-and-resubmit.","headline":"Useful tutorial with a real sign error in the canard-location formula and an unsupported claim about rigorous homoclinic verification.","tokens_in":21148,"tokens_out":2992,"would_cite":false,"duration_ms":25883,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34E15","34C26","37G15","37C29","92C20"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper aims to prove that the fast-slow FitzHugh-Nagumo system combines relaxation oscillations, Hopf-induced periodic solutions, homoclinic orbits, and canards, with explicit formulas for the oscillation period and the canard…","keywords":["FitzHugh-Nagumo system","fast-slow systems","canards","canard explosion","relaxation oscillations","Hopf bifurcation","homoclinic bifurcation","singular perturbation theory"],"falsifier":"For $\\varepsilon=0.5$, the paper's own numerical figure places rapid canard-length growth near $c \\approx 1.150077$, whereas the formula $\\bar c_c = -3\\varepsilon/8$ with $\\bar c = 2/\\sqrt{3} - c$ puts the canard at $c \\approx 1.342$; continuing the periodic orbits of system (23) numerically and locating the actual headless-canard value of $c$ would settle which prediction is correct.","tokens_in":20031,"feed_emoji":"🧠","tokens_out":10797,"duration_ms":99119,"temperature":0.7,"pith_summary":"This tutorial-style paper seeks to prove, rather than merely observe, that the two-variable FitzHugh-Nagumo fast-slow system (derived from the classic four-variable nerve-impulse model) exhibits three intertwined dynamical mechanisms: relaxation oscillations, periodic orbits born at Hopf bifurcations, and canard orbits that connect them. The authors work out the singular limit, apply the standard persistence theory for slow manifolds, and then check bifurcation conditions by hand. They also provide quantitative formulas, including a closed-form estimate for the relaxation-oscillation period and a canard parameter value in terms of the small time-scale ratio. If the proofs hold, the paper converts a qualitative picture of excitable dynamics into concrete parameter predictions for a classical mathematical-neuroscience model.","feed_headline":"Proof links relaxation cycles, Hopf orbits, and canards in one model","feed_subtitle":"Closed-form period estimate and canard value make quantitative predictions for the classic nerve-impulse model.","key_machinery":"The central objects are the cubic critical manifold $C_0 = \\{(x,y): y = 4x - x^3\\}$, the explicit slow flow in the fast variable $x$ given by $\\dot x = (b x^3 + (1-4b)x + c)/(4-3x^2)$, and the canard normal form (27) around a regular singular fold. A regular singular fold is a fold point of the critical manifold that is simultaneously an equilibrium of the slow flow, with both the $x$-dependence and the parameter-dependence of the slow equation nonzero there. The normal form expresses the system near such a point as $x' = -y l_1 + x^2 l_2 + \\varepsilon l_3$ and $y' = \\varepsilon(\\pm x l_4 - \\lambda l_5 + y l_6)$, and its coefficients $l_1,\\dots,l_6$ determine the constants $A$ and $B$ that give the Hopf value $\\lambda_H = -B\\varepsilon$ and the canard value $\\lambda_c = -(B+A)\\varepsilon$. The paper's computations fix these coefficients in both parameter slices, obtaining $A=-3/8$ in Example 6.1 and $A<0$ in Example 6.2, which is what fixes the direction of the Hopf bifurcation and the location of the canard explosion.","core_discovery":"The central claim is that the FitzHugh-Nagumo fast-slow system (5), with fast equation $x' = -y + 4x - x^3$ and slow equation $y' = \\varepsilon(x - by - c)$, is a fully worked case study of geometric singular perturbation theory. In the singular limit $\\varepsilon = 0$, the paper proves the existence of a relaxation cycle whose period is estimated by $T_{\\gamma_0} \\approx 2\\int_{4/\\sqrt{3}}^{2/\\sqrt{3}} \\frac{4-3x^2}{bx^3+(1-4b)x+c}\\,dx$. For $0<\\varepsilon\\ll1$, it argues that the regular system retains this cycle with $O(\\varepsilon)$ period error. In the parameter slice $b=0$, it identifies a subcritical Hopf bifurcation at $c_H = 2/\\sqrt{3}$ and, via the canard normal form, derives the canard location $\\bar c_c = -3\\varepsilon/8$ in the translated coordinate $\\bar c = c_H - c$. In the slice $c=0$, it finds a pitchfork at $b=1/4$, a Hopf point at $b_H^\\varepsilon = (-4+\\sqrt{16+3\\varepsilon})/\\varepsilon$, and a homoclinic bifurcation that destroys the unstable periodic orbit; numerically, for $\\varepsilon=0.5$, this occurs at $b \\approx 0.36932$.","pith_inferences":["The paper leaves implicit that its quoted $\\bar c_c = -3\\varepsilon/8$ conflicts in sign with the normal-form formula $\\lambda_c = -(B+A)\\varepsilon$ when $B=0$ and $A=-3/8$; resolving this sign or rescaling issue would change the quantitative canard prediction.","The same normal-form computation could be applied to the two coupled FitzHugh-Nagumo systems mentioned in Section 7, potentially predicting which coupling strengths produce mixed-mode oscillations without full numerical simulation.","The period integral (14) can be evaluated in closed form for other $(b,c)$ parameter pairs, turning the singular limit into testable frequency-amplitude relations for excitable cells.","A natural continuation is to track the homoclinic value $b_{\\rm Hom}^\\varepsilon$ as $\\varepsilon \\to 0$ and check whether it approaches $3/8$ with a scaling consistent with canard explosion theory."],"forward_implications":["For $0<\\varepsilon\\ll1$, the regular relaxation cycle persists with period $T_{\\gamma_0} + O(\\varepsilon)$, so formula (14) gives a quantitative prediction for the spike interval in the nerve-impulse model.","In the slice $b=0$, the subcritical Hopf bifurcation sits exactly at $c_H = 2/\\sqrt{3}$, and stable small-amplitude cycles exist for $c$ just below this value, with the canard value indicating where the cycle expands into a large relaxation loop.","In the slice $c=0$, the pitchfork at $b=1/4$, the supercritical Hopf at $b_H^\\varepsilon$, and the homoclinic bifurcation at $b \\approx 0.36932$ for $\\varepsilon=0.5$ organize the birth, growth, and destruction of the unstable periodic orbit.","The normal-form coefficient $A=-3/8$ in Example 6.1 fixes the bifurcation as supercritical in the translated parameter and subcritical in the original parameter $c$, confirming the shape of the numerically observed branch of periodic solutions.","Since all conditions of the persistence theorem and the canard normal form are checked explicitly, the paper positions the FitzHugh-Nagumo system as a template in which every step from singular limit to canard explosion is made explicit."],"supporting_citations":[{"why":"Supplies the original Hodgkin-Huxley nerve-impulse model from which the two-variable FitzHugh-Nagumo system is reduced.","marker":"[18]"},{"why":"States the physiological reduction to two variables, providing the biological grounding for the cubic and affine vector fields used in system (5).","marker":"[20]"},{"why":"Gives the specific fast-slow form of the FitzHugh-Nagumo system studied here and the observation of two distinct dynamics near the Hopf bifurcation.","marker":"[6]"},{"why":"Provides the persistence theorem for normally hyperbolic slow manifolds that links the singular dynamics to the regular case.","marker":"[7]"},{"why":"Supplies the textbook treatment of slow manifolds, fold points, and the canard normal form used throughout the paper.","marker":"[25]"},{"why":"Establishes the fold-and-canard normal form in two dimensions that Example 6.1 and Example 6.2 use to compute A and B.","marker":"[22]"},{"why":"Provides the relaxation-oscillation and canard-explosion results that the paper imports to describe the growth of periodic orbits.","marker":"[23]"},{"why":"Supplies the Hopf and homoclinic bifurcation conditions and the coordinate-change methodology used in the normal-form derivations.","marker":"[15]"},{"why":"The Hopf theorem used to prove the existence of periodic solutions at the bifurcation values.","marker":"[19]"}],"fun_headline_variants":["Nerve model gets exact canard location and cycle period","FitzHugh-Nagumo tutorial nails canards and relaxation cycles","Closed-form canard value and period for the nerve impulse model","One model ties relaxation cycles, Hopf, and canards"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The canard-location formula rests on the assumption that a single translation of coordinates and parameter puts the system into the exact normal form whose canard expansion is known, with no additional rescaling; if that identification is off by a scaling or a sign, the formula $\\bar c_c = -3\\varepsilon/8$ does not apply to the original parameter $c$.","fun_headline_variants_meta":{"raw":{"variants":["Nerve model gets exact canard location and cycle period","FitzHugh-Nagumo tutorial nails canards and relaxation cycles","Closed-form canard value and period for the nerve impulse model","One model ties relaxation cycles, Hopf, and canards"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000215,"raw_usage":{"total_tokens":1434,"prompt_tokens":954,"completion_tokens":480,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":570,"completion_tokens_details":{"reasoning_tokens":408}},"tokens_in":570,"tokens_out":480,"duration_ms":5315,"temperature":1.0,"reasoning_tokens":408,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:49:56.202138+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For $\\varepsilon=0.5$, the paper's own numerical figure places rapid canard-length growth near $c \\approx 1.150077$, whereas the formula $\\bar c_c = -3\\varepsilon/8$ with $\\bar c = 2/\\sqrt{3} - c$ puts the canard at $c \\approx 1.342$; continuing the periodic orbits of system (23) numerically and locating the actual headless-canard value of $c$ would settle which prediction is correct.","supporting_citations":[{"cited_title":"Hopf, Abzweigung einer periodischen Losung von einer stationaren Losung eines Differ- entialsystems, Ber","cited_arxiv_id":null,"evidence_quote":"The Hopf theorem used to prove the existence of periodic solutions at the bifurcation values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the original Hodgkin-Huxley nerve-impulse model from which the two-variable FitzHugh-Nagumo system is reduced."},{"cited_title":"Keener and J","cited_arxiv_id":null,"evidence_quote":"States the physiological reduction to two variables, providing the biological grounding for the cubic and affine vector fields used in system (5)."},{"cited_title":"Doss-Bachelet, J.-P","cited_arxiv_id":null,"evidence_quote":"Gives the specific fast-slow form of the FitzHugh-Nagumo system studied here and the observation of two distinct dynamics near the Hopf bifurcation."},{"cited_title":"Fenichel, Geometric singular perturbation theory for ordinary differential equations, Jour- nal of Differential Equations 31 (1979), no","cited_arxiv_id":null,"evidence_quote":"Provides the persistence theorem for normally hyperbolic slow manifolds that links the singular dynamics to the regular case."},{"cited_title":"Kuehn, Multiple time scale dynamics , Springer–Verlag, 2015.l","cited_arxiv_id":null,"evidence_quote":"Supplies the textbook treatment of slow manifolds, fold points, and the canard normal form used throughout the paper."},{"cited_title":"Krupa and P","cited_arxiv_id":null,"evidence_quote":"Establishes the fold-and-canard normal form in two dimensions that Example 6.1 and Example 6.2 use to compute A and B."},{"cited_title":"Krupa and P","cited_arxiv_id":null,"evidence_quote":"Provides the relaxation-oscillation and canard-explosion results that the paper imports to describe the growth of periodic orbits."},{"cited_title":"Guckenheimer and P","cited_arxiv_id":null,"evidence_quote":"Supplies the Hopf and homoclinic bifurcation conditions and the coordinate-change methodology used in the normal-form derivations."}],"review_version":1}