Recognition: unknown
From the Volterra type Lyapunov functions of Rahman-Zou towards a competitive exclusion partition property for rank one models
Pith reviewed 2026-05-09 16:58 UTC · model grok-4.3
The pith
Rank-one next-generation matrices let a Perron-Volterra Lyapunov function divide parameter space into four regions each with a unique globally stable equilibrium.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For multi-strain models with rank-one next-generation matrices, a Perron-Volterra Lyapunov function built from resident Volterra entropy plus an invader linear functional certifies global asymptotic stability once the balance identity holds; when the incidence functions are additionally concave and increasing, the two-strain case yields the competitive exclusion partition property: the parameter space decomposes into four open regions, each containing a unique globally asymptotically stable equilibrium (disease-free, resident-only, invader-only, or coexistence).
What carries the argument
Perron-Volterra Lyapunov function: Volterra entropy on resident variables plus Perron-weighted linear functional on invaders, derived from the left Perron eigenvector of the transversal Jacobian, with a balance identity that cancels coupling terms.
If this is right
- Global stability of every boundary equilibrium reduces to recursive computation of invasion numbers on the siphon lattice.
- The same Lyapunov construction applies to any finite collection of singleton strains.
- Models with one scalar strain and one irreducible rank-one block inherit the same partition property.
- An algorithmic procedure in the EpidCRN package constructs the candidate Lyapunov functions, verifies the balance identity, and partitions the parameter space.
- A local Lyapunov theorem holds for stability on each siphon face.
Where Pith is reading between the lines
- The obstruction that appears for two rank-one blocks suggests that cross-equilibrium terms must be added when the next-generation matrix is no longer rank one.
- The recursive siphon-lattice construction may extend to models whose incidence functions satisfy weaker monotonicity conditions if the balance identity can be recovered by other means.
- The framework supplies an explicit certificate that could be checked numerically for higher-dimensional rank-one blocks before attempting full global analysis.
Load-bearing premise
The incidence functions must be concave and increasing and the next-generation matrix must have exact rank one, otherwise the balance identity used to cancel coupling terms may fail.
What would settle it
A concrete two-strain model with concave increasing incidence whose next-generation matrix has rank one, yet whose parameter space contains a region with two distinct stable equilibria or a region with no stable equilibrium.
read the original abstract
This paper presents a Perron-Volterra framework that unifies explicit Lyapunov constructions for multi-strain epidemic models with rank-one next-generation matrices. At each boundary equilibrium on a siphon face, the Lyapunov function consists of a Volterra entropy on resident variables plus a Perron-weighted linear functional on invaders, derived from the left Perron eigenvector of the transversal Jacobian. A balance identity cancels coupling terms, reducing global stability to recursive computation of invasion numbers on the siphon lattice. For two-strain models with concave, increasing incidence, we prove the competitive exclusion partition property (CEPP): the parameter space splits into four open regions, each possessing a unique globally asymptotically stable equilibrium (disease-free, single-strain, or coexistence) certified by an explicit Lyapunov function. The same mechanism extends to an arbitrary number of singleton strains and to models with one scalar strain and one irreducible rank-one block. We implement the algorithmic approach in the Mathematica package EpidCRN, which constructs candidate Lyapunov functions, verifies the balance identity, and partitions the parameter space recursively. For two rank-one matrix blocks, the standard ansatz fails; we characterize the obstruction and propose an augmented cross-equilibrium Lyapunov function. A local Lyapunov theorem for siphon faces is also provided. The framework offers a systematic stability analysis of rank-one models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a Perron-Volterra framework for explicit Lyapunov constructions in multi-strain epidemic models with rank-one next-generation matrices. At boundary equilibria on siphon faces, the Lyapunov function combines Volterra entropy on resident variables with a Perron-weighted linear functional on invaders, derived from the left Perron eigenvector of the transversal Jacobian. A balance identity is invoked to cancel coupling terms, reducing global stability certification to recursive computation of invasion numbers. For two-strain models with concave, increasing incidence functions, the authors establish the competitive exclusion partition property (CEPP): the parameter space partitions into four open regions, each containing a unique globally asymptotically stable equilibrium (disease-free, single-strain, or coexistence). The approach is implemented in the Mathematica package EpidCRN for automated construction, verification, and partitioning; extensions to arbitrary singleton strains, one scalar strain plus one irreducible rank-one block, and a local Lyapunov theorem for siphon faces are also presented.
Significance. If the central derivations hold, the work supplies a systematic, first-principles method for certifying global stability via explicit Lyapunov functions in rank-one models, unifying earlier Volterra-type constructions and reducing the problem to invasion-number calculations on the siphon lattice. The EpidCRN implementation provides a concrete, reproducible tool for algorithmic verification, and the characterization of the obstruction for two rank-one blocks plus the local Lyapunov theorem for siphon faces are useful additions to the literature on dynamical systems in epidemiology.
major comments (2)
- [Section deriving the balance identity and the CEPP theorem] The balance identity used to cancel all coupling terms in the Lyapunov derivative (reducing dV/dt to a linear combination of invasion numbers) is asserted to hold for general concave increasing incidence under the rank-one assumption. However, the explicit derivation showing that state-dependent Jacobian entries for strictly concave forms (e.g., saturated incidence) produce no residual term of indefinite sign is not supplied in sufficient detail; this step is load-bearing for the sign control in the CEPP theorem and the global-stability claims in all four parameter regions.
- [CEPP theorem and EpidCRN implementation section] The recursive partitioning of parameter space via invasion numbers on the siphon lattice is central to the CEPP. While the abstract states that proofs exist, the manuscript should include at least one fully worked numerical example with a specific strictly concave incidence function to verify that the identity cancels identically and that the resulting dV/dt sign is controlled solely by the invasion numbers.
minor comments (2)
- Notation for the siphon lattice and transversal Jacobian could be accompanied by a small diagram or table summarizing the resident/invader decomposition at each boundary equilibrium.
- The description of the augmented cross-equilibrium Lyapunov function for two rank-one blocks would benefit from an explicit formula or small example illustrating the augmentation.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions planned for the updated version.
read point-by-point responses
-
Referee: The balance identity used to cancel all coupling terms in the Lyapunov derivative (reducing dV/dt to a linear combination of invasion numbers) is asserted to hold for general concave increasing incidence under the rank-one assumption. However, the explicit derivation showing that state-dependent Jacobian entries for strictly concave forms (e.g., saturated incidence) produce no residual term of indefinite sign is not supplied in sufficient detail; this step is load-bearing for the sign control in the CEPP theorem and the global-stability claims in all four parameter regions.
Authors: We appreciate the referee highlighting the need for greater explicitness in this derivation. The balance identity follows from the rank-one structure of the next-generation matrix and the monotonicity/concavity of the incidence functions: the left Perron eigenvector of the transversal Jacobian ensures that the state-dependent partial derivatives of the incidence terms, when contracted against the weighted variables, cancel identically along trajectories due to the homogeneity properties and the fact that the incidence is increasing and concave. While the manuscript presents the identity and its consequences, we agree that the intermediate steps for strictly concave cases were not expanded sufficiently. In the revision we will add a detailed subsection (or appendix) that computes the relevant Jacobian entries explicitly for a general concave increasing incidence, shows the cancellation of all state-dependent residuals, and confirms that the derivative reduces precisely to the linear combination of invasion numbers with no indefinite remainder. revision: yes
-
Referee: The recursive partitioning of parameter space via invasion numbers on the siphon lattice is central to the CEPP. While the abstract states that proofs exist, the manuscript should include at least one fully worked numerical example with a specific strictly concave incidence function to verify that the identity cancels identically and that the resulting dV/dt sign is controlled solely by the invasion numbers.
Authors: We agree that a concrete numerical verification would strengthen the exposition and allow readers to check the cancellation directly. In the revised manuscript we will insert a fully worked example using a specific strictly concave incidence (e.g., the saturated form βS I / (1 + α I) with chosen parameter values). The example will (i) compute the invasion numbers on the siphon lattice, (ii) construct the Perron-Volterra Lyapunov function, (iii) differentiate it explicitly to verify that the balance identity holds with no residual terms, and (iv) confirm that the sign of dV/dt is governed solely by the invasion numbers, thereby partitioning the parameter space into the four open regions each containing a unique globally asymptotically stable equilibrium. The EpidCRN package will be used to automate the construction and verification steps within the example. revision: yes
Circularity Check
Derivation self-contained via explicit Perron-Volterra construction and balance identity
full rationale
The paper constructs the Lyapunov function explicitly from the left Perron eigenvector of the transversal Jacobian at each boundary equilibrium, then applies a balance identity (derived from the rank-one next-generation matrix and the eigenvector property) to cancel cross terms in the derivative. Global stability then reduces to sign conditions on invasion numbers computed recursively on the siphon lattice. For the two-strain concave-increasing case this yields the CEPP partition into four regions, each with a unique GAS equilibrium. The construction is algorithmic (implemented in EpidCRN) and does not rely on fitting parameters to data or on self-citations whose content is presupposed; the balance identity is asserted to hold under the stated assumptions without reducing the final stability claim to a tautology. No load-bearing step collapses to a renaming, ansatz smuggling, or fitted-input prediction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Incidence functions are concave and strictly increasing.
- domain assumption The next-generation matrix is exactly rank one.
Reference graph
Works this paper leans on
-
[1]
and Vassena, Nicola , title =
Golnik, Richard and Gatter, Thomas and Stadler, Peter F. and Vassena, Nicola , title =. 2025 , eprint =
2025
-
[2]
Nonlinearity , volume=
The lifting bifurcation problem on feed-forward networks , author=. Nonlinearity , volume=. 2018 , publisher=
2018
-
[3]
Journal of Mathematical Biology , year =
Sanz-Lorenzo, Luis and Bravo de la Parra, Rafael and Poggiale, Jean-Christophe and Auger, Pierre , title =. Journal of Mathematical Biology , year =
-
[4]
SIAM Journal on Applied Mathematics , volume=
Mass action systems: two criteria for Hopf bifurcation without Hurwitz , author=. SIAM Journal on Applied Mathematics , volume=. 2025 , publisher=
2025
-
[5]
, title =
Angeli, David and Sontag, Eduardo D. , title =. IEEE Transactions on Automatic Control , volume =
-
[6]
preprint , year=
Codimension-three Andronov-Hopf bifurcations in infectious diseases affecting juveniles , author=. preprint , year=
-
[7]
Proceedings of the Royal Society B: Biological Sciences , volume=
Immunological serotype interactions and their effect on the epidemiological pattern of dengue , author=. Proceedings of the Royal Society B: Biological Sciences , volume=. 2009 , publisher=
2009
-
[8]
Communications in Nonlinear Science and Numerical Simulation , volume=
Prey herding and predators’ feeding satiation induce multiple stability , author=. Communications in Nonlinear Science and Numerical Simulation , volume=. 2023 , publisher=
2023
-
[9]
bioRxiv , pages=
A mathematical model for cancer dynamics with treatment and saboteur bacteria , author=. bioRxiv , pages=. 2025 , publisher=
2025
-
[10]
IEEE Transactions on Control of Network Systems , year=
On metzler positive systems on hypergraphs , author=. IEEE Transactions on Control of Network Systems , year=
-
[11]
IEEE Transactions on Automatic Control , year=
Analysis of higher-order Lotka-Volterra models: Application of S-tensors and the polynomial complementarity problem , author=. IEEE Transactions on Automatic Control , year=
-
[12]
2023 62nd IEEE Conference on Decision and Control (CDC) , pages=
Species coexistence and extinction resulting from higher-order Lotka-Volterra two-faction competition , author=. 2023 62nd IEEE Conference on Decision and Control (CDC) , pages=. 2023 , organization=
2023
-
[13]
SIAM Journal on Applied Dynamical Systems , volume=
Symmetry groupoids and patterns of synchrony in coupled cell networks , author=. SIAM Journal on Applied Dynamical Systems , volume=. 2003 , publisher=
2003
-
[14]
SIAM Journal on Control and Optimization , volume=
Convergence and Equilibria Analysis of a Networked Bivirus Epidemic Model , author=. SIAM Journal on Control and Optimization , volume=. 2022 , publisher=
2022
-
[15]
1996 , publisher=
Global dynamical properties of Lotka-Volterra systems , author=. 1996 , publisher=
1996
-
[16]
, author=
Lotka-Volterra Dynamics — An Introduction. , author=. 2010 , publisher=
2010
-
[17]
arXiv preprint arXiv:2202.08707 , year=
Synchronization induced by directed higher-order interactions , author=. arXiv preprint arXiv:2202.08707 , year=
-
[18]
Theoretical Ecology , volume=
Higher order interactions and species coexistence , author=. Theoretical Ecology , volume=. 2021 , publisher=
2021
-
[19]
Proceedings of the National Academy of Sciences , volume=
Analytical note on certain rhythmic relations in organic systems , author=. Proceedings of the National Academy of Sciences , volume=. 1920 , publisher=
1920
-
[20]
ICES Journal of Marine Science , volume=
Variations and fluctuations of the number of individuals in animal species living together , author=. ICES Journal of Marine Science , volume=. 1928 , publisher=
1928
-
[21]
Journal of Mathematical Biology , volume=
Global stability in two species interactions , author=. Journal of Mathematical Biology , volume=. 1976 , publisher=
1976
-
[22]
The American Naturalist , volume=
Stability in models of mutualism , author=. The American Naturalist , volume=. 1979 , publisher=
1979
-
[23]
Journal of Mathematical Analysis and Applications , volume=
The stability of generalized Volterra equations , author=. Journal of Mathematical Analysis and Applications , volume=. 1978 , publisher=
1978
-
[24]
The American Naturalist , volume=
Arguments in favor of higher order interactions , author=. The American Naturalist , volume=. 1983 , publisher=
1983
-
[25]
Nature ecology & evolution , volume=
Higher-order interactions capture unexplained complexity in diverse communities , author=. Nature ecology & evolution , volume=. 2017 , publisher=
2017
-
[26]
Ecology letters , volume=
The mechanistic basis for higher-order interactions and non-additivity in competitive communities , author=. Ecology letters , volume=. 2019 , publisher=
2019
-
[27]
Handbook of differential equations: ordinary differential equations , volume=
Monotone dynamical systems , author=. Handbook of differential equations: ordinary differential equations , volume=. 2006 , publisher=
2006
-
[28]
A survey of results , author=
Systems of ordinary differential equations which generate an order preserving flow. A survey of results , author=. SIAM review , volume=. 1988 , publisher=
1988
-
[29]
On the Endemic Behavior of a Competitive Tri-Virus
Gracy, Sebin and Ye, Mengbin and Anderson, Brian and Uribe, Cesar A , journal=. On the Endemic Behavior of a Competitive Tri-Virus
-
[30]
Chaos, Solitons & Fractals , volume=
Bifurcation analysis in the diffusive Lotka--Volterra system: An application to market economy , author=. Chaos, Solitons & Fractals , volume=. 2009 , publisher=
2009
-
[31]
International Journal of Educational Research and Social Sciences (IJERSC) , volume=
Stability Analysis of Lotka-Volterra Model in The Case of Interaction of Local Religion and Official Religion , author=. International Journal of Educational Research and Social Sciences (IJERSC) , volume=
-
[32]
Discrete applied mathematics , volume=
Directed hypergraphs and applications , author=. Discrete applied mathematics , volume=. 1993 , publisher=
1993
-
[33]
Systems & Control Letters , volume=
Contraction analysis of virtually positive systems , author=. Systems & Control Letters , volume=. 2022 , publisher=
2022
-
[34]
SIAM Journal on Applied Mathematics , volume=
Competing subcommunities of mutualists and a generalized Kamke theorem , author=. SIAM Journal on Applied Mathematics , volume=. 1986 , publisher=
1986
-
[35]
Proceedings of the National Academy of Sciences , volume=
Coexistence in diverse communities with higher-order interactions , author=. Proceedings of the National Academy of Sciences , volume=. 2022 , publisher=
2022
-
[36]
and Spong, Mark W
Mohammadi, A. and Spong, Mark W. , journal=. Chetaev Instability Framework for Kinetostatic Compliance-Based Protein Unfolding , year=
-
[37]
IEEE Transactions on Automatic Control , year=
Explicit Solutions and Stability Properties of Homogeneous Polynomial Dynamical Systems , author=. IEEE Transactions on Automatic Control , year=
-
[38]
SIAM Journal on Matrix Analysis and Applications , volume=
M-tensors and some applications , author=. SIAM Journal on Matrix Analysis and Applications , volume=. 2014 , publisher=
2014
-
[39]
Linear Algebra and Its Applications , volume=
M-tensors and nonsingular M-tensors , author=. Linear Algebra and Its Applications , volume=. 2013 , publisher=
2013
-
[40]
Journal of Scientific Computing , volume=
Solving multi-linear systems with M-tensors , author=. Journal of Scientific Computing , volume=. 2016 , publisher=
2016
-
[41]
Multigroup
Cisneros-Velarde, Pedro and Bullo, Francesco , journal=. Multigroup. 2021 , publisher=
2021
-
[42]
Proyecciones (Antofagasta) , volume=
The Chetaev theorem for ordinary difference equations , author=. Proyecciones (Antofagasta) , volume=. 2012 , publisher=
2012
-
[43]
Numerical Linear Algebra with Applications , volume=
A survey on the spectral theory of nonnegative tensors , author=. Numerical Linear Algebra with Applications , volume=. 2013 , publisher=
2013
-
[44]
SIAM Journal on Matrix Analysis and Applications , volume=
Further results for Perron--Frobenius theorem for nonnegative tensors II , author=. SIAM Journal on Matrix Analysis and Applications , volume=. 2011 , publisher=
2011
-
[45]
An invariance principle for discontinuous dynamic systems with application to a coulomb friction oscillator , author=. J. Dyn. Sys., Meas., Control , volume=
-
[46]
Proceedings of the National Academy of Sciences , volume=
Understanding bistability in complex enzyme-driven reaction networks , author=. Proceedings of the National Academy of Sciences , volume=. 2006 , publisher=
2006
-
[47]
SIAM Journal on Applied Dynamical Systems , volume=
Local and global stability of equilibria for a class of chemical reaction networks , author=. SIAM Journal on Applied Dynamical Systems , volume=. 2013 , publisher=
2013
-
[48]
2009 European Control Conference (ECC) , pages=
A tutorial on Chemical Reaction Networks dynamics , author=. 2009 European Control Conference (ECC) , pages=. 2009 , organization=
2009
-
[49]
The Journal of Physical Chemistry A , volume=
Autonomous discovery of unknown reaction pathways from data by chemical reaction neural network , author=. The Journal of Physical Chemistry A , volume=. 2021 , publisher=
2021
-
[50]
IFAC-PapersOnLine , volume=
On the stability of the endemic equilibrium of a discrete-time networked epidemic model , author=. IFAC-PapersOnLine , volume=. 2020 , publisher=
2020
-
[51]
arXiv preprint arXiv:2206.07425 , year=
Discrete-time layered-network epidemics model with time-varying transition rates and multiple resources , author=. arXiv preprint arXiv:2206.07425 , year=
-
[52]
IEEE Transactions on Automatic Control , volume=
Analysis and control of a continuous-time bi-virus model , author=. IEEE Transactions on Automatic Control , volume=. 2019 , publisher=
2019
-
[53]
IEEE Transactions on Control Systems Technology , volume=
Analysis, estimation, and validation of discrete-time epidemic processes , author=. IEEE Transactions on Control Systems Technology , volume=. 2018 , publisher=
2018
-
[54]
Analysis, Control, and State Estimation for the Networked Competitive Multi-Virus
Zhang, Ciyuan and Gracy, Sebin and Basar, Tamer and Pare, Philip E , journal=. Analysis, Control, and State Estimation for the Networked Competitive Multi-Virus
-
[55]
IEEE Transactions on Control of Network Systems , year=
Estimation and distributed eradication of SIR epidemics on networks , author=. IEEE Transactions on Control of Network Systems , year=
-
[56]
Node-based
Liu, Fangzhou and Buss, Martin , booktitle=. Node-based. 2016 , organization=
2016
-
[57]
at-Automatisierungstechnik , volume=
Optimal filtering and control of network information epidemics , author=. at-Automatisierungstechnik , volume=. 2021 , publisher=
2021
-
[58]
Nature communications , volume=
Simplicial models of social contagion , author=. Nature communications , volume=. 2019 , publisher=
2019
- [59]
-
[60]
Physical review letters , volume=
Universal nonlinear infection kernel from heterogeneous exposure on higher-order networks , author=. Physical review letters , volume=. 2021 , publisher=
2021
-
[61]
Applied Mathematics and Computation , volume=
Competing spreading dynamics in simplicial complex , author=. Applied Mathematics and Computation , volume=. 2022 , publisher=
2022
-
[62]
ESAIM: Proceedings and surveys , volume=
An introduction to the basic reproduction number in mathematical epidemiology , author=. ESAIM: Proceedings and surveys , volume=. 2018 , publisher=
2018
-
[63]
SIAM Review , volume=
What are higher-order networks? , author=. SIAM Review , volume=. 2023 , publisher=
2023
-
[64]
Communications in Mathematical Sciences , volume=
Perron-Frobenius theorem for nonnegative tensors , author=. Communications in Mathematical Sciences , volume=. 2008 , publisher=
2008
-
[65]
SIAM Journal on Matrix Analysis and Applications , volume=
Further results for Perron--Frobenius theorem for nonnegative tensors , author=. SIAM Journal on Matrix Analysis and Applications , volume=. 2010 , publisher=
2010
-
[66]
2011 50th IEEE conference on decision and control and European control conference , pages=
Distributed control of positive systems , author=. 2011 50th IEEE conference on decision and control and European control conference , pages=. 2011 , organization=
2011
-
[67]
Journal of Symbolic Computation , volume=
Eigenvalues of a real supersymmetric tensor , author=. Journal of Symbolic Computation , volume=. 2005 , publisher=
2005
-
[68]
Computational Optimization and Applications , volume=
Iterative methods for computing U-eigenvalues of non-symmetric complex tensors with application in quantum entanglement , author=. Computational Optimization and Applications , volume=. 2020 , publisher=
2020
-
[69]
and Nedić, Angelia and Tang, Choon Yik and Beck, Carolyn L
Liu, Ji and Paré, Philip E. and Nedić, Angelia and Tang, Choon Yik and Beck, Carolyn L. and Başar, Tamer , journal=. Analysis and Control of a Continuous-Time Bi-Virus Model , year=
-
[70]
IEEE Transactions on Automatic Control , volume=
Hurwitz stability of Metzler matrices , author=. IEEE Transactions on Automatic Control , volume=. 2010 , publisher=
2010
-
[71]
TenEig - Tensor Eigenpairs Solver , author=
-
[72]
SIAM Journal on Matrix Analysis and Applications , volume=
Computing tensor eigenvalues via homotopy methods , author=. SIAM Journal on Matrix Analysis and Applications , volume=. 2016 , publisher=
2016
-
[73]
2005 , publisher=
H-infinity control and estimation of state-multiplicative linear systems , author=. 2005 , publisher=
2005
-
[74]
The Journal of Engineering , volume=
LQR controller design for quad-rotor helicopters , author=. The Journal of Engineering , volume=. 2019 , publisher=
2019
-
[75]
IEEE Transactions on Network Science and Engineering , volume=
Controllability of hypergraphs , author=. IEEE Transactions on Network Science and Engineering , volume=. 2021 , publisher=
2021
-
[76]
Linear and Multilinear Algebra , volume=
Spectral directed hypergraph theory via tensors , author=. Linear and Multilinear Algebra , volume=. 2016 , publisher=
2016
-
[77]
SIAM Journal on Applied Algebra and Geometry , volume=
Polynomial dynamical systems, reaction networks, and toric differential inclusions , author=. SIAM Journal on Applied Algebra and Geometry , volume=. 2019 , publisher=
2019
-
[78]
2008 , publisher=
Distributed control and analysis of coupled cell systems , author=. 2008 , publisher=
2008
-
[79]
Geometry, mechanics, and dynamics , pages=
Patterns of oscillation in coupled cell systems , author=. Geometry, mechanics, and dynamics , pages=. 2002 , publisher=
2002
-
[80]
SIAM Journal on Mathematics of Data Science , volume=
Three hypergraph eigenvector centralities , author=. SIAM Journal on Mathematics of Data Science , volume=. 2019 , publisher=
2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.