{"id":"c974cf45-6770-4e1c-8bf0-a04d38723e0e","arxiv_id":"2411.18059","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper claims existence and stability of transitory canard and relaxation oscillation cycles through a degenerate transcritical point in a slow-fast Leslie-Gower predator-prey model with weak Allee effect.","lead":"This paper analyzes a predator-prey model where prey reproduce much faster than predators, and identifies new types of oscillation cycles that pass through a degenerate singular point. The result matters because it extends geometric singular perturbation theory to degenerate transcritical bifurcations in ecological models, though several key formulas contain algebraic errors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The degenerate Hopf value QH in (15) is algebraically wrong: the U1^2 denominator term has the wrong sign, giving QH ≈ 11.45 instead of ≈ 2.10 for A = 1/2, M = -1/10, S = 1/20, so Theorem 4(2) is stated at a parameter value that is not the canard point.","rationale":"The reader's weakest assumption is that QH in (11) and (15) is algebraically incorrect, and my independent derivation confirms this. The trace expression derived in Appendix A, when combined with the equilibrium equation (6) and the degenerate condition C = -AMQ, gives a denominator term (M+1-A)U1^2, not (M-1-A)U1^2 as printed in (15). The numerical discrepancy is not marginal: for the paper's own simulation parameter values A = 1/2, M = -1/10, S = 1/20, the correct degenerate canard value is about 2.10, while (15) yields about 11.45. This is load-bearing because Theorem 4(2) explicitly invokes QH from (15) to define the canard point, and Figure 15 uses Q = QH - delta to display the canard explosion. At Q ≈ 11.45, the equilibrium is not near the fold point P, so the singular canard orbit underlying Theorem 4(2) is not the one being perturbed. The error also propagates into the numerical validation, although the MatCont bifurcation diagram is computed independently and may still be correct. I do not see evidence of a fatal flaw in the blow-up strategy itself; the main theorem appears repairable by replacing (15) with the correct QH and rerunning the simulations. For this reason I recommend CONDITIONAL acceptance rather than outright rejection: the paper should not be accepted without correcting the QH formula and re-verifying the numerical claims. A secondary concern, also noted by the reader, is that Proposition 3 asserts existence of an open parameter region based on numerical slices of sigma rather than an analytical proof; this is another gap but is secondary to the QH error because the sign of sigma is evaluated at the wrong canard point in the current text.","tokens_in":29464,"tokens_out":23507,"duration_ms":189403,"concrete_test":"Independently recompute QH for the degenerate case directly from tr(J) = 0 and the equilibrium equation (6): solve U h'(U) = S(U+A) for the positive root near up with A = 1/2, M = -1/10, S = 1/20, then set Q = U/(h(U)+AM) and compare with formula (15) at the same U. If the two values differ by more than 10% (they differ by roughly a factor of five), formula (15) is wrong; rerun Figure 15 with Q chosen near the corrected QH to check whether the canard explosion and the cycles through TC are actually located at the claimed parameter values.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The main canard-cycle theorem, Theorem 4(2), explicitly identifies the canard point with Q = QH given in (15). That formula is incorrect. Using the paper's own trace expression (9) and equilibrium equation (6), as in Appendix A, the degenerate case C = -AMQ gives tr(J) = (U1+C)[U1((A-M-1)U1 + 2(1/Q - A + M - AM) + 1/Q) - S(U1+A)]. Setting tr = 0 and eliminating C via (6) yields the Hopf value QH = 3U1 / [S(U1+A) + (M+1-A)U1^2 + 2(A-M+AM)U1], equivalently Q = U1/(h(U1)+AM) at the root of U1 h'(U1) = S(U1+A). Equation (15), however, has (M-A-1)U1^2 in the denominator instead of (M+1-A)U1^2, an algebraic error of 2U1^2. For A = 1/2, M = -1/10, S = 1/20, the relevant positive root is U1 ≈ 0.545, for which the correct QH is about 2.10, while (15) gives about 11.45. Consequently, at the Q-value named in Theorem 4(2), the equilibrium is far from the fold point P, so the hypothesized canard point does not exist there, and the canard explosion shown in Figure 15 is parametrized around the wrong center. The GSPT/blow-up construction may be repairable by substituting the correct QH, but the central quantitative claim as written is not supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a rescaled Leslie–Gower predator–prey model with a weak Allee effect and a fast prey/slow predator structure. The main claims are: stability thresholds and Hopf values for the unique positive equilibrium in a generic and in a degenerate transcritical case; an intrinsic determination of the criticality of the singular Hopf bifurcation following De Maesschalck et al.; a blow-up desingularization of the degenerate point TC; existence and stability of relaxation oscillations and transitory canard cycles near the fold and TC; and a numerical bifurcation analysis that locates a nearby Takens–Bogdanov point. The central theoretical object is Theorem 4(2), which asserts that for C = -AMQ and A-M+AM > 1/Q, a stable canard cycle exists for Q = Qc(ε) near QH, with QH given by Eq. (15).","tokens_in":29795,"tokens_out":26701,"duration_ms":235986,"significance":"If the main claims were correct, the paper would make a useful contribution to the slow-fast predator-prey literature: the blow-up treatment of the degenerate transcritical point TC goes beyond the generic analysis in Zhu and Liu, the use of the intrinsic criticality formula of [4] is appropriate and avoids normal-form computations, and the numerical Takens-Bogdanov check in Appendix B is a genuine consistency test rather than a fitted result. The analysis is self-contained and transparently derives the trace and determinant expressions. However, the central quantitative claim is not supported as written because the Hopf value QH used throughout is algebraically incorrect, and the stability statement in Theorem 4(2) omits a hypothesis that the paper itself shows is necessary. The conceptual blow-up framework appears salvageable, but the current manuscript cannot be accepted in its present form.","major_comments":[{"comment":"The Hopf values in Eqs. (11) and (15) contain an algebraic sign error. Setting tr(J)=0 in Eq. (14) and solving for Q gives, in the generic case, QH = 3(U1+C) / [S(U1+A) + 2(A-M+AM)U1 - 3AM + (M+1-A)U1^2]; in the degenerate case C=-AMQ this reduces to QH = 3U1 / [S(U1+A) + 2(A-M+AM)U1 + (M+1-A)U1^2]. Equation (15) instead has (M-A-1)U1^2 in the denominator, and Eq. (11) has (M-1-A)U1^2; both differ from the correct coefficient (M+1-A)U1^2 by -2U1^2. For A=1/2, M=-1/10, S=1/20, the correct degenerate Hopf value near the fold is approximately 2.1, whereas Eq. (15) evaluated at the fold coordinate U1≈0.582 gives approximately 11.45. At Q≈11.45 the positive equilibrium is at U1≈0.91, far from the fold point P, so the identification of Q=QH with a canard point in Theorem 4(2) is false as stated, and the canard explosion shown in Figure 15 is parametrized around the wrong value of Q. This error is load-bearing for the thresholds in Theorems 2 and 3 and for the main canard-cycle theorem.","section":"Sections 3.1 and 3.2, Eqs. (11) and (15)"},{"comment":"Theorem 4(2) asserts the existence of a locally stable canard cycle for every parameter tuple satisfying C=-AMQ and A-M+AM>1/Q, without imposing supercriticality of the singular Hopf bifurcation. The paper's own Section 5 shows that the criticality coefficient σ in Eq. (31) changes sign (Figure 8), so the canard cycle is not stable in the subcritical regime. The proof of Proposition 12(2) explicitly depends on an 'appropriate choice of parameters' that makes the singular Hopf bifurcation supercritical. The theorem statement should include this hypothesis (for example, σ<0, or parameters in the open set of Proposition 3), and the stability conclusion should be restricted accordingly.","section":"Section 6.5, Theorem 4(2)"},{"comment":"Proposition 3 asserts the existence of an open set of parameters for which the Hopf bifurcation is supercritical, but the supporting evidence is numerical evaluation of σ in Figure 8; no analytic computation or proof is supplied. Since the stability assertion in Theorem 4(2) depends on this existence, the manuscript should either provide a proof (for example, by evaluating σ and its derivatives at a concrete parameter point) or state the claim as a numerical observation rather than as a proposition. As written, the formal statement is not established.","section":"Section 5, Proposition 3"}],"minor_comments":[{"comment":"The phrase 'solving for wJ11' should read 'solving for U1J11'; the derivation text also promises an expression equivalent to 1/Q - J11 but displays U1J11, which is confusing.","section":"Appendix A"},{"comment":"'Haussdorf' should be 'Hausdorff' in both occurrences.","section":"Theorems 4 and 5"},{"comment":"The panel references in Theorem 4 appear inconsistent with the caption: if the left panel is the transitory canard and the right panel is the relaxation oscillation, then Theorem 4(1) should refer to the right panel and Theorem 4(2) to the left panel.","section":"Figure 14 and Theorem 4"},{"comment":"The phrase 'the (complex) eigenvalues of J have negative real parts negative' contains a duplicated word and should be corrected.","section":"Section 3.1, after Theorem 2"},{"comment":"The notation switches between S and ε within the same displayed formulas; since S=ε is set only later, the notation should be harmonized.","section":"Appendix B, Eqs. (65)-(66)"},{"comment":"'codimension2bifurcation' should be 'codimension-2 bifurcation'.","section":"Figure 17 caption"}],"recommendation":"major_revision","confidential_remarks":"The algebraic error in the Hopf value is central and affects Theorems 2, 3, Theorem 4(2), and the numerical validation in Figure 15, but it appears to be a mechanical sign mistake that can be corrected by substituting the correct denominator. The blow-up construction and the use of the intrinsic criticality criterion are conceptually sound and likely repairable. I would recommend requiring the authors to correct the formulas, regenerate the numerical figures with the corrected QH, and add the missing supercriticality hypothesis before the paper can be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has a genuinely interesting core: it blow-ups a degenerate transcritical point TC in the slow-fast Leslie–Gower model and shows that, when the degenerate condition C = -AMQ holds, the point organizes two distinct regimes depending on the sign of A-M+AM - 1/Q. In one regime you get a transitory canard and a relaxation oscillation passing through TC. That is a real extension of Zhong and Shen and of Zhu and Liu, and the blow-up analysis itself looks careful. The use of the De Maesschalck intrinsic criticality formula is appropriate, and the numerical bifurcation work with the Takens–Bogdanov check is a nice complement.\n\nThe problem is the algebra behind the canard point. Theorem 3 and Theorem 4(2) use QH as given in (15). That formula is wrong. I checked it against the paper's own trace expression (14) and the equilibrium equation (6). For C = -AMQ, setting tr(J)=0 gives a denominator with (1+M-A)U1^2, not (M-A-1)U1^2. The difference is 2U1^2, so the sign is off. For the paper's own numerical values, A = 1/2, M = -1/10, S = 1/20, the correct degenerate Hopf value is about 2.10, while (15) gives about 11.45. That is not a typo; it changes the location of the canard point by a factor of five. Consequently, Theorem 4(2) is stated at a parameter value where the equilibrium is not at the fold, and the canard explosion in Figure 15 is centered around the wrong Q. The central quantitative claim is not supported as written.\n\nTwo smaller issues compound this. Proposition 3 calls numerical slices of the criticality coefficient a proof; as written it proves the existence of an open parameter set only at the level of the numerics shown. And Theorem 4(2) states the existence of a locally stable cycle without explicitly stating the supercriticality condition on the first Lyapunov coefficient, which the paper does establish only numerically for a region. These are fixable, but they are part of the same need-for-revision pattern.\n\nThe blow-up construction and the general idea are probably repairable: substituting the correct QH into the statements and re-running the simulations would likely restore the claims. But this is a major revision, not a minor edit.\n\nThis paper deserves a serious referee, but only because the underlying scenario is worth someone's time to get right. If I were the editor, I would send it to review, and I would expect the referee to send it back for major revision.","headline":"Interesting blow-up analysis of a degenerate transcritical point, but the central QH formula is algebraically wrong, so the main canard-cycle theorem is stated at the wrong parameter value.","tokens_in":30418,"tokens_out":7247,"would_cite":false,"duration_ms":52214,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34E15","34C23","37G10","92D25"],"pacs":[],"model":"deepseek-v4-flash","headline":"In the degenerate case $C=-AMQ$, the paper proves that a degenerate transcritical point plus a fold/canard point organizes stable relaxation and canard cycles that converge to singular cycles as $\\varepsilon\\to0$.","keywords":["slow-fast system","Leslie-Gower model","weak Allee effect","canard cycles","relaxation oscillations","singular Hopf bifurcation","blow-up method","degenerate transcritical bifurcation"],"falsifier":"Take the paper's own parameter values $A=1/2$, $M=-1/10$, $S=1/20$, set $U_1=u_p$ from (21), and solve the trace equation (14) for $Q$; compare with $Q_H$ from (15). A disagreement of the size reported by a direct check (about 2.10 versus 11.45) would show the quantitative canard location is not as claimed, though it would not by itself disprove the qualitative existence of a cycle.","tokens_in":29162,"feed_emoji":"🔄","tokens_out":10476,"duration_ms":89490,"temperature":0.7,"pith_summary":"This paper studies a Leslie-Gower predator-prey model with a Holling type II response and a weak Allee effect on the prey, in the regime where the prey evolves much faster than the predator, so the system is slow-fast. Earlier work found canard cycles and relaxation oscillations in the generic case; this paper's contribution is the degenerate case in which the slow and fast nullclines meet on the predator axis, creating a degenerate transcritical point $TC$. It claims that, when a certain parameter inequality holds, $TC$ acts not as an attractor but as a saddle-like organizing center, and together with the fold point $P$ it produces a stable relaxation cycle and a stable transitory canard cycle for small $\\varepsilon$. The proof uses blow-up desingularization, transition maps, and an intrinsic criticality formula for the singular Hopf bifurcation; the authors also numerically locate and analytically confirm a nearby codimension-two bifurcation point. If correct, the paper shows that degenerate transcritical singularities can drive sustained predator-prey oscillations rather than simple extinction.","feed_headline":"Degenerate predator-prey point still organizes stable cycles","feed_subtitle":"A slow-fast Leslie-Gower model is shown to host stable relaxation and canard cycles through a degenerate transcritical point.","key_machinery":"The load-bearing object is the degenerate transcritical point $TC=(0,-AM)$ that appears when $C=-AMQ$, where the slow nullcline meets the fast nullcline on the $v$-axis; at this point the vector field is nilpotent, so standard slow-fast theory does not apply directly. The machinery is the blow-up transformation $(u,v,\\varepsilon)=(r\\bar u,r\\bar v,r\\bar\\varepsilon)$, which resolves the nilpotent point into charts (entry, central, exit, bottom) containing only semi-hyperbolic equilibria and center manifolds. In these charts the transitions $\\Sigma_0\\to\\Sigma_2$ and $\\Sigma_1\\to\\Sigma_2$ are contractions, and together with contraction away from the singularities this yields stable cycles via a contraction-map argument. A second piece of machinery is the intrinsic formula from [4] for the criticality of a slow-fast Hopf point, used to show the Hopf bifurcation at $Q=Q_H$ is supercritical without computing a normal form.","core_discovery":"The central claim is Theorem 4: for $C=-AMQ$ and $A-M+AM>1/Q$, if the fold point $P$ is a generic jump point there is a locally stable relaxation cycle $\\gamma_\\varepsilon$ converging in Hausdorff distance to the singular cycle $\\gamma_0$ as $\\varepsilon\\to0$; if instead $P$ is a canard point at $Q=Q_H$, and $Q=Q_c(\\varepsilon)\\approx Q_H$ is chosen for a maximal canard, there is a locally stable transitory canard cycle $\\tilde\\gamma_\\varepsilon$ converging to $\\tilde\\gamma_0$. The singular cycles pass through the degenerate transcritical point $TC$, so the paper's slogan is that $TC$ organizes the oscillatory dynamics rather than being an absorbing extinction point. The paper further claims, via a normal-form-free intrinsic calculation, that the singular Hopf bifurcation at $Q_H$ is supercritical in an open parameter region, so the small cycles are attracting. It also proves existence and uniqueness of the generic relaxation oscillation when $C<-AMQ$ using an entry-exit function, and identifies the location of the codimension-two nilpotent bifurcation point analytically.","pith_inferences":["Because the proof of Theorem 4 is parameterized by $Q_H$, a corrected trace-zero identity would shift the quantitative location of the canard cycle; the qualitative existence argument via blow-up could still survive.","The same four-chart blow-up proof could be adapted to other predator-prey models whose nullclines meet degenerately on an invariant axis, potentially turning 'attractor' transcritical points into organizing centers for oscillations.","Ecologically, the saddle-like behavior of $TC$ implies long bouts of near-extinction of prey alternating with recovery; measuring residence time near $TC$ in simulations could serve as a direct check.","The intrinsic criticality formula could be applied to locate supercritical versus subcritical Hopf boundaries in higher codimension settings where normal-form computations are impractical."],"forward_implications":["For $C=-AMQ$ and $A-M+AM>1/Q$, a stable relaxation oscillation exists through $TC$ when the fold is a generic jump point (Theorem 4(1)).","When the fold is a canard point, a stable transitory canard cycle exists for $\\varepsilon$ small and $Q$ near $Q_H$, converging to the singular canard as $\\varepsilon\\to0$ (Theorem 4(2)).","The singular Hopf bifurcation at $Q_H$ is supercritical for an open set of parameters, so the bifurcating small-amplitude cycles are stable.","In the complementary regime $A-M+AM<1/Q$, the degenerate point $TC$ is attracting and no oscillatory behavior is organized by it.","A codimension-two nilpotent bifurcation point exists near $TC$, and the numerical bifurcation diagram shows canard-explosion cycles ending in a homoclinic orbit when multiple equilibria are present."],"supporting_citations":[{"why":"Establishes the baseline slow-fast analysis of this model, including canard cycles and relaxation oscillations in the non-degenerate case, which this paper extends.","marker":"[42]"},{"why":"Supplies the intrinsic formula used to determine the criticality of the singular Hopf bifurcation without a normal form.","marker":"[4]"},{"why":"Provides the fold and canard point theory and transition maps used in the proof of Theorem 4.","marker":"[21]"},{"why":"Identifies the degenerate transcritical point as an attractor in a related model, the scenario this paper shows can instead organize cycles.","marker":"[41]"},{"why":"Introduces the model and its equilibria and bifurcations, the starting point for the parameter setting and Theorem 1.","marker":"[1]"},{"why":"Gives the normal hyperbolicity theory that justifies the slow and fast manifolds away from $P$ and $TC$.","marker":"[11]"},{"why":"Provides the entry-exit function used to locate the exit point in the generic relaxation oscillation proof.","marker":"[6]"},{"why":"Numerical bifurcation software used for the codimension-two diagram that locates and confirms the nilpotent bifurcation point.","marker":"[7]"}],"fun_headline_variants":["Degenerate transcritical point spawns stable relaxation cycles","Singular point organizes canard cycles in predator-prey system","Stable cycles traced to degenerate blow-up singularity","Slow-fast model: degenerate point dictates oscillation stability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central calculation that pins down the exact value of the bifurcation parameter $Q_H$ at which the equilibrium sits on the fold/canard point is algebraically correct; if that value is wrong, the claimed canard cycle and its numerical simulations occur at the wrong parameter values.","fun_headline_variants_meta":{"raw":{"variants":["Degenerate transcritical point spawns stable relaxation cycles","Singular point organizes canard cycles in predator-prey system","Stable cycles traced to degenerate blow-up singularity","Slow-fast model: degenerate point dictates oscillation stability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000767,"raw_usage":{"total_tokens":3417,"prompt_tokens":979,"completion_tokens":2438,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":2384}},"tokens_in":595,"tokens_out":2438,"duration_ms":16983,"temperature":1.0,"reasoning_tokens":2384,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:37:35.097481+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the paper's own parameter values $A=1/2$, $M=-1/10$, $S=1/20$, set $U_1=u_p$ from (21), and solve the trace equation (14) for $Q$; compare with $Q_H$ from (15). A disagreement of the size reported by a direct check (about 2.10 versus 11.45) would show the quantitative canard location is not as claimed, though it would not by itself disprove the qualitative existence of a cycle.","supporting_citations":[{"cited_title":"Zhu and X","cited_arxiv_id":null,"evidence_quote":"Establishes the baseline slow-fast analysis of this model, including canard cycles and relaxation oscillations in the non-degenerate case, which this paper extends."},{"cited_title":"De Maesschalck, T","cited_arxiv_id":null,"evidence_quote":"Supplies the intrinsic formula used to determine the criticality of the singular Hopf bifurcation without a normal form."},{"cited_title":"Zhong and J","cited_arxiv_id":null,"evidence_quote":"Identifies the degenerate transcritical point as an attractor in a related model, the scenario this paper shows can instead organize cycles."},{"cited_title":"Arancibia-Ibarra and J","cited_arxiv_id":null,"evidence_quote":"Introduces the model and its equilibria and bifurcations, the starting point for the parameter setting and Theorem 1."},{"cited_title":"Fenichel","cited_arxiv_id":null,"evidence_quote":"Gives the normal hyperbolicity theory that justifies the slow and fast manifolds away from $P$ and $TC$."},{"cited_title":"De Maesschalck and S","cited_arxiv_id":null,"evidence_quote":"Provides the entry-exit function used to locate the exit point in the generic relaxation oscillation proof."},{"cited_title":"Dhooge, W","cited_arxiv_id":null,"evidence_quote":"Numerical bifurcation software used for the codimension-two diagram that locates and confirms the nilpotent bifurcation point."}],"review_version":1}