REVIEW 3 major objections 6 minor 42 references
Singular bifurcations in a modified Leslie-Gower model
T0 review · 3 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read 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$.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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).
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 (3)
- [Sections 3.1 and 3.2, Eqs. (11) and (15)] 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 6.5, Theorem 4(2)] 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 5, Proposition 3] 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.
minor comments (6)
- [Appendix A] 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.
- [Theorems 4 and 5] 'Haussdorf' should be 'Hausdorff' in both occurrences.
- [Figure 14 and Theorem 4] 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 3.1, after Theorem 2] The phrase 'the (complex) eigenvalues of J have negative real parts negative' contains a duplicated word and should be corrected.
- [Appendix B, Eqs. (65)-(66)] The notation switches between S and ε within the same displayed formulas; since S=ε is set only later, the notation should be harmonized.
- [Figure 17 caption] 'codimension2bifurcation' should be 'codimension-2 bifurcation'.
Circularity Check
No significant circularity: the derivation is self-contained; the suspected QH algebra error is a correctness issue, not a self-referential reduction.
full rationale
The paper's central claims are derived from the model rather than from fitted parameters or self-citations. QH in (15) is obtained by setting tr(J)=0 in (9)/(14), i.e. from the equilibrium's Jacobian, and is then used as the slow-fast Hopf / canard parameter; this is a consistency check with the model, not a prediction forced by an input. The criticality calculation in Section 5 applies the external intrinsic formula of De Maesschalck et al. [4] and verifies conditions (30) at P; the sign of sigma is explored numerically but is not used as an input to define the Hopf value. The canard-cycle Theorem 4 rests on standard GSPT transition maps from Krupa-Szmolyan [21] and the blow-up construction, whose details are standard and cited to [24] and (for the survey) to the co-authored [19]; [19] is not load-bearing for the theorem. The Takens-Bogdanov appendix solves tr=det=0 from J(E1) and then compares with the MatCont value, which is a verification, not a circular prediction. The skeptic's algebraic objection to (15) — the sign of the U1^2 term in the denominator — is a mathematical correctness concern: if (15) is wrong, Theorem 4(2) is evaluated at the wrong Q-value, but the formula still comes from the paper's own trace expression and does not reduce to the conclusion it purports to prove. Thus no circular step is exhibited.
Assumptions & free parameters
assumptions (4)
- domain assumption The prey growth rate is much larger than the predator growth rate, so S = ε is a small parameter.
- domain assumption The fold point is a generic jump point or a canard point, and the equilibrium lies on the repelling branch with 0 < U1 < up in Theorem 4.
- ad hoc to paper Non-degeneracy A-M+AM > 0 and A-M+AM > 1/Q for the degenerate blow-up analysis.
- standard math The intrinsic criticality formula of De Maesschalck et al. [4] is valid for the system.
Cite this review
Pith. "Pith review of Singular bifurcations in a modified Leslie-Gower model." pith.science (2026). https://pith.science/paper/CHCZXXP4
@misc{pith2026241118059,
author = {Pith},
title = {Pith review of: Singular bifurcations in a modified Leslie-Gower model},
year = {2026},
howpublished = {\url{https://pith.science/paper/CHCZXXP4}},
note = {Machine review of arXiv:2411.18059}
}
abstract
We study a predator-prey system with a generalist Leslie-Gower predator, a functional Holling type II response, and a weak Allee effect on the prey. The prey's population often grows much faster than its predator, allowing us to introduce a small time scale parameter $\varepsilon$ that relates the growth rates of both species, giving rise to a slow-fast system. Zhu and Liu (2022) show that, in the case of the weak Allee effect, Hopf singular bifurcation, slow-fast canard cycles, relaxation oscillations, etc., exist. Our main contribution lies in the rigorous analysis of a degenerate scenario organized by a (degenerate) transcritical bifurcation. The key tool employed is the blow-up method that desingularizes the degenerate singularity. In addition, we determine the criticality of the singular Hopf bifurcation using recent intrinsic techniques that do not require a local normal form. The theoretical analysis is complemented by a numerical bifurcation analysis, in which we numerically identify and analytically confirm the existence of a nearby Takens-Bogdanov point.
Figures
Figures from the paper (15 more)
Reference graph
Works this paper leans on
-
[4]
P. De Maesschalck, T. S. Doan, and J. Wynen. Intrinsic determination of the criticality of a slow-fast Hopf bifurcation. J. Dynam. Differential Equations, 33(4):2253–2269, 2021
work page 2021
-
[1]
C. Arancibia-Ibarra and J. Flores. Dynamics of a Leslie-Gower predator-prey model with Holling type II func- tional response, Allee effect and a generalist predator.Math. Comput. Simul., 188(2):1–22, 2021
work page 2021
-
[2]
G. Birkhoff and G.-C. Rota.Ordinary differential equations. John Wiley & Sons, Inc., New York, fourth edition, 1989
work page 1989
-
[3]
F. Courchamp, L. Berec, and J. Gascoigne.Allee effects in ecology and conservation. Oxford University Press, Oxford, 2008
work page 2008
-
[5]
P. De Maesschalck, F. Dumortier, and R. Roussarie.Canard cycles—from birth to transition, volume 73 of Ergebnisse der Mathematik und ihrer Grenzgebiete. 3. Folge. A Series of Modern Surveys in Mathematics [Results in Mathematics and Related Areas. 3rd Series. A Series of Modern Surveys in Mathematics]. Springer, Cham, [2021] ©2021
work page 2021
-
[6]
P. De Maesschalck and S. Schecter. The entry-exit function and geometric singular perturbation theory.J. Differential Equations, 260(8):6697–6715, 2016
work page 2016
- [7]
-
[8]
F. Dumortier. Techniques in the theory of local bifurcations: blow-up, normal forms, nilpotent bifurcations, singular perturbations. InBifurcations and periodic orbits of vector fields (Montreal, PQ, 1992), volume 408 of NATO Adv. Sci. Inst. Ser. C: Math. Phys. Sci., pages 19–73. Kluwer Acad. Publ., Dordrecht, 1993
work page 1992
Show all 42 references
-
[9]
Dumortier and R
F. Dumortier and R. Roussarie. Canard cycles and center manifolds.Mem. Amer. Math. Soc., 121(577), 1996
1996
-
[10]
K. Fang, Z. Zhu, F. Chen, and Z. Li. Qualitative and bifurcation analysis in a Leslie-Gower model with Allee effect. Qual. Theory Dyn. Syst., 21(3):Paper No. 86, 19, 2022
2022
-
[11]
Fenichel
N. Fenichel. Geometric singular perturbation theory for ordinary differential equations.J. Differential Equations, 31(1):53–98, 1979. SINGULAR BIFURCATIONS IN A MODIFIED LESLIE-GOWER MODEL 39
1979
-
[12]
H. I. Freedman.Deterministic mathematical models in population ecology, volume 57 ofMonographs and Text- books in Pure and Applied Mathematics. Marcel Dekker, Inc., New York, 1980
1980
-
[13]
Gao and S
Y. Gao and S. Yao. Dynamical analysis of a modified Leslie-Gower Holling-type II predator-prey stochastic model in polluted environments with interspecific competition and impulsive toxicant input.J. Biol. Dyn., 16(1):840–858, 2022
2022
-
[14]
Gonzalez-Olivares and A
E. Gonzalez-Olivares and A. Rojas-Palma. Global stability in a modified Leslie-Gower type predation model assuming mutual interference among generalist predators.Math. Biosci. Eng., 17(6):7708–7731, 2020
2020
-
[15]
Guckenheimer and P
J. Guckenheimer and P. Holmes.Nonlinear oscillations, dynamical systems, and bifurcations of vector fields, volume 42 ofApplied Mathematical Sciences. Springer-Verlag, New York, 1990
1990
-
[16]
G. Hek. Geometric singular perturbation theory in biological practice.J. Math. Biol., 60:347–86, 05 2009
2009
-
[17]
C. Holling. The components of predation as revealed by a study of small-mammal predation of the european sawfly.Can. Entomol., 91:293–320, 1959
1959
-
[18]
C. Holling. Some characteristics of simple types of predation and parasitism.Can. Entomol., 91:385–398, 1959
1959
-
[19]
Jardón-Kojakhmetov and C
H. Jardón-Kojakhmetov and C. Kuehn. A survey on the blow-up method for fast-slow systems. InMexican mathematicians in the world—trends and recent contributions, volume 775 ofContemp. Math., pages 115–160. Amer. Math. Soc., [Providence], RI, 2021
2021
-
[20]
Q. J. A. Khan, E. Balakrishnan, and G. C. Wake. Analysis of a predator-prey system with predator switching. Bull. Math. Biol., 66(1):109–123, 2004
2004
-
[21]
Krupa and P
M. Krupa and P. Szmolyan. Extending geometric singular perturbation theory to nonhyperbolic points—fold and canard points in two dimensions.SIAM J. Math. Anal., 33(2):286–314, 2001
2001
-
[22]
Krupa and P
M. Krupa and P. Szmolyan. Extending slow manifolds near transcritical and pitchfork singularities.Nonlinearity, 14(6):1473–1491, 2001
2001
-
[23]
Krupa and P
M. Krupa and P. Szmolyan. Geometric analysis of the singularly perturbed planar fold. InMultiple-time-scale dynamical systems (Minneapolis, MN, 1997), volume 122 ofIMA Vol. Math. Appl., pages 89–116. Springer, New York, 2001
1997
-
[24]
Kuehn.Multiple time scale dynamics, volume 191
C. Kuehn.Multiple time scale dynamics, volume 191. Springer, 2015
2015
-
[25]
P. H. Leslie and J. C. Gower. The properties of a stochastic model for the predator-prey type of interaction between two species.Biometrika, 47:219–234, 1960
1960
-
[26]
J. Li, S. Li, and X. Wang. Canard, homoclinic loop, and relaxation oscillations in a Lotka-Volterra system with Allee effect in predator population.Chaos, 33(7):Paper No. 073130, 14, 2023
2023
-
[27]
A. J. Lotka.Elements of Physical Biology. Williams and Wilkins, Baltimore, 1925
1925
-
[28]
Ludwig, D
D. Ludwig, D. D. Jones, and C. S. Holling.Qualitative analysis of insect outbreak systems: the spruce budworm and forest, pages 547–564. University of Chicago Press, Chicago, 2022
2022
-
[29]
Muratori and S
S. Muratori and S. Rinaldi. Remarks on competitive coexistence.SIAM J. Appl. Math., 49(5):1462–1472, 1989
1989
-
[30]
Murray.Mathematical biology
J. Murray.Mathematical biology. Springer-Verlag, NY, 2nd edition edition, 2002
2002
-
[31]
Qiu and S
H. Qiu and S. Guo. Bifurcation structures of a Leslie-Gower model with diffusion and advection.Appl. Math. Lett., 135:Paper No. 108391, 6, 2023
2023
-
[32]
Sahoo and G
D. Sahoo and G. Samanta. Oscillatory and transient dynamics of a slow-fast predator-prey system with fear and its carry-over effect.Nonlinear Anal. Real World Appl., 73:Paper No. 103888, 39, 2023
2023
-
[33]
S. J. Schreiber. Allee effects, extinctions, and chaotic transients in simple population models.Theoretical Pop- ulation Biology, 64(2):201–209, 2003
2003
-
[34]
J. Shen. Canard limit cycles and global dynamics in a singularly perturbed predator-prey system with non- monotonic functional response.Nonlinear Anal. Real World Appl., 31:146–165, 2016
2016
-
[35]
N. C. Stenseth, W. Falck, O. Bjornstad, and C. Krebs. Population regulation in snowshoe hare and canadian lynx: Asymmetric food web configurations between hare and lynx.Proceedings of the National Academy of Sciences of the United States of America, 94:5147–52, 06 1997
1997
-
[36]
F. Takens. Constrained equations; a study of implicit differential equations and their discontinuous solutions. In P. Hilton, editor,Structural Stability, the Theory of Catastrophes, and Applications in the Sciences, pages 143–234, Berlin, Heidelberg, 1976. Springer Berlin Heidelberg
1976
-
[37]
Volterra
V. Volterra. Variazioni e fluttuazioni del numero d’individui in specie animali conviventi.Mem. Acad. Lincei, 6(2):31–113, 1926
1926
-
[38]
Wen and T
Z. Wen and T. Shi. Existence and Uniqueness of a Canard Cycle with Cyclicity at Most Two in a Singularly Per- turbed Leslie–Gower Predator–Prey Model with Prey Harvesting.Internat. J. Bifur. Chaos Appl. Sci. Engrg., 34(3):Paper No. 2450036, 2024. 40 SINGULAR BIFURCATIONS IN A ...
2024
-
[39]
J. Yao, J. Huang, and R. Huzak. Cyclicity of slow–fast cycles with two canard mechanisms.Chaos, 34(5):Paper No. 053112, 11, 2024
2024
-
[40]
Yao and R
J. Yao and R. Huzak. Cyclicity of the limit periodic sets for a singularly perturbed leslie–gower predator–prey model with prey harvesting.Journal of Dynamics and Differential Equation, 2022
2022
-
[41]
Zhong and J
L. Zhong and J. Shen. Degenerate transcritical bifurcation point can be an attractor: A case study in a slow–fast modified leslie–gower model.Qual. Theory Dyn. Syst., 21(76), 2022
2022
-
[42]
Zhu and X
Z. Zhu and X. Liu. Canard cycles and relaxation oscillations in a singularly perturbed Leslie-Gower predator-prey model with Allee effect.Internat. J. Bifur. Chaos Appl. Sci. Engrg., 32(5):Paper No. 2250071, 23, 2022. Departamento de Matemáticas, UAM–Iztapalapa, 09310 Iztapala...
2022
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.