pith. sign in

arxiv: 2603.06778 · v3 · submitted 2026-03-06 · 🧬 q-bio.MN · math.DS

A cocktail of chemical reaction networks and mathematical epidemiology tools for positive ODE stability problems

Florin Avram , Rim Adenane , Andrei-Dan Halanay This is my paper

Pith reviewed 2026-05-15 15:33 UTC · model grok-4.3

classification 🧬 q-bio.MN math.DS
keywords chemical reaction networksnext generation matrixmathematical epidemiologypositive ODEsstability analysisstoichiometric matrixbifurcation analysis
0
0 comments X

The pith

A chemical reaction network generalization of the next generation matrix theorem provides stability criteria for positive ordinary differential equations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper combines tools from chemical reaction network theory and mathematical epidemiology to address stability problems in positive systems of ordinary differential equations. It presents a generalization of the well-known next generation matrix theorem that incorporates concepts from chemical reaction networks. This approach is then applied to review a symbolic-numeric method for analyzing bifurcations by examining the characteristic polynomial through minors of the stoichiometric matrix. Applications in both fields are discussed using specific software tools.

Core claim

The authors provide a CRN-flavored generalization of the Next Generation Matrix theorem, which is the most cited result in mathematical epidemiology, allowing for the analysis of stability in positive ODE models that can be represented as chemical reaction networks.

What carries the argument

The generalized Next Generation Matrix criterion based on the stoichiometric matrix and reaction rates satisfying specific structural hypotheses.

If this is right

  • Stability analysis can be performed by checking conditions on the stoichiometric matrix similar to the NGM approach.
  • The method applies to models in both mathematical epidemiology and chemical reaction networks.
  • The symbolic-numeric approach identifies coefficients of the characteristic polynomial as Child Selection minors.
  • Applications can be implemented using the Epid-CRN tools package in Mathematica.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This generalization may simplify stability proofs for complex biological systems by leveraging existing CRN results.
  • It could extend to other areas like systems biology where positive ODEs appear.
  • Testing the criterion on known stable and unstable models would validate its applicability.

Load-bearing premise

The epidemiological or chemical models can be faithfully represented as chemical reaction networks whose stoichiometric matrix and reaction rates satisfy the structural hypotheses needed for the generalized NGM criterion.

What would settle it

A counterexample would be a positive ODE system modeled as a CRN where the generalized NGM criterion predicts stability but the system is actually unstable, or vice versa.

Figures

Figures reproduced from arXiv: 2603.06778 by Andrei-Dan Halanay, Florin Avram, Rim Adenane.

Figure 1
Figure 1. Figure 1: Contraction and enlargement between the red dashed SIRWS network on left NSIRWS and the Boros–Rost witness NBR on right. The projection π collapses the immune pathway i → R → W into a direct interaction i → W, which is represented however as a consumption, since νiW consumes W, but dotted, to indicate catalysis by i. The enlargement E (in the sense of Banaji–Boros–Hofbauer) reverses this, inserting R as an… view at source ↗
Figure 2
Figure 2. Figure 2: Flow diagram of model (9). Λ, µs, µi, µr are the inflow rate of S and the per capita death rates of S, I, R, respectively.The linear conversion reactions S −→γss R, R −→γrr S, I −→iri R, I −→isi S model vaccinations, loss (waning) of immunity, recovery, and very fast loss of immunity. The infection and treatment rates S + I −→ βsiF (i) 2I, I −→T [i] R are functional. Remark 18. Note that the fact that the … view at source ↗
read the original abstract

We continue recent attempts to put together concepts and results of Chemical Reaction Networks theory (CRNT) and Mathematical Epidemiology (ME), for solving problems of stability of positive ODEs. We provide first an elegant CRN-flavored generalization of the most cited result in ME, the Next Generation Matrix (NGM) theorem. We review next the "symbolic-numeric approach of Vassena and Stadler, which tackles bifurcation problems by viewing the characteristic polynomial of the Jacobian at fixed points as a formal polynomial in the "symbolic reactivities", and identifies its coefficients as "Child Selection minors of the stoichiometric matrix". We also review two applications of this approach using the Mathematica package Epid-CRN tools from both CRNT and ME.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a CRN-flavored generalization of the Next Generation Matrix (NGM) theorem for stability analysis of positive ODEs, reviews the symbolic-numeric approach of Vassena and Stadler that expresses the Jacobian characteristic polynomial coefficients as Child Selection minors of the stoichiometric matrix, and presents two applications of this framework using the Epid-CRN tools Mathematica package.

Significance. If the generalization is rigorously derived and the required structural hypotheses hold for standard compartmental models, the work would provide a useful bridge between chemical reaction network theory and mathematical epidemiology, enabling unified stability and bifurcation criteria for positive systems. The software tools and explicit minor-based coefficient identification add practical reproducibility value.

major comments (2)
  1. [Abstract] Abstract and introduction: the claimed 'elegant CRN-flavored generalization' of the NGM theorem is load-bearing for the central contribution, yet the text provides no explicit statement or verification that common epidemiological models admit faithful CRN representations whose stoichiometric matrix and rate functions satisfy the precise conditions (deficiency zero, complex balance, or required sign patterns on Jacobian minors) under which the generalized criterion applies.
  2. [Applications] Applications section: the two reviewed applications using Epid-CRN tools must include explicit checks confirming that the target models meet the structural hypotheses of the generalized NGM result; absent such checks, the claim that the approach solves positive ODE stability problems for ME models cannot be assessed.
minor comments (2)
  1. Define 'Child Selection minors' and 'symbolic reactivities' at first use with a brief reminder of their relation to the stoichiometric matrix.
  2. Add a short table or paragraph contrasting the new CRN-NGM criterion with the classical NGM statement to clarify the precise extension.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. The comments identify important points where additional clarity and explicit verification will strengthen the presentation of the CRN-flavored generalization of the NGM theorem and its applications. We address each major comment below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract] Abstract and introduction: the claimed 'elegant CRN-flavored generalization' of the NGM theorem is load-bearing for the central contribution, yet the text provides no explicit statement or verification that common epidemiological models admit faithful CRN representations whose stoichiometric matrix and rate functions satisfy the precise conditions (deficiency zero, complex balance, or required sign patterns on Jacobian minors) under which the generalized criterion applies.

    Authors: We agree that the manuscript would benefit from an explicit statement of the precise conditions under which the generalized NGM result holds. The original text presents the generalization but does not include the requested verifications for standard models. In the revised manuscript we will add a dedicated paragraph in the introduction that states the required conditions on the stoichiometric matrix, rate functions, and sign patterns of the relevant minors, followed by explicit checks confirming that common compartmental models (e.g., SIR and SEIR) admit faithful CRN representations satisfying these hypotheses. revision: yes

  2. Referee: [Applications] Applications section: the two reviewed applications using Epid-CRN tools must include explicit checks confirming that the target models meet the structural hypotheses of the generalized NGM result; absent such checks, the claim that the approach solves positive ODE stability problems for ME models cannot be assessed.

    Authors: We concur that explicit checks are necessary to substantiate the applicability claims. The current applications section reviews the use of the Epid-CRN tools but does not document the required structural verifications. In the revised version we will insert, for each of the two applications, a short subsection that uses the package to confirm the models satisfy the hypotheses of the generalized NGM result (including deficiency-zero status where applicable and the sign patterns on the child-selection minors). These checks will be presented with the corresponding output from the symbolic-numeric procedure. revision: yes

Circularity Check

0 steps flagged

No circularity: CRN generalization of NGM theorem is structurally independent

full rationale

The paper states a CRN-flavored generalization of the standard Next Generation Matrix theorem as its central result, framed under explicit structural hypotheses on stoichiometric matrices and reaction rates drawn from CRNT. No derivation step reduces by construction to a fitted parameter, self-citation chain, or renamed input; the NGM generalization is presented as a new structural criterion rather than a re-expression of prior fitted quantities. The review of the Vassena-Stadler symbolic-numeric method is cited as external work. The derivation chain remains self-contained against the stated assumptions and does not collapse to tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Based on the abstract alone the paper relies on standard structural assumptions from CRNT (mass-action kinetics, stoichiometric matrix properties) and ME (next-generation matrix framework) without introducing new free parameters or invented entities.

axioms (2)
  • domain assumption Models can be represented as chemical reaction networks with non-negative stoichiometric coefficients and mass-action kinetics.
    Invoked to allow the CRN-flavored generalization of the NGM theorem.
  • domain assumption The Jacobian at the disease-free equilibrium admits a block structure compatible with the next-generation matrix construction.
    Standard assumption in ME that is carried over to the generalized setting.

pith-pipeline@v0.9.0 · 5427 in / 1354 out tokens · 38576 ms · 2026-05-15T15:33:36.335508+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

  1. [1]

    [AAFG23] Rim Adenane, Florin Avram, Mohamed El Fatini, and RP Gupta,Dynamics of an sir epidemic model with limited medical resources, revisited and corrected, arXiv preprint arXiv:2310.19696 (2023). [AAHJ24] Florin Avram, Rim Adenane, Andrei D Halanay, and Matthew D Johnston,Stability in reaction network models via an extension of the next generation matr...

    work page arXiv 2023

  2. [2]

    2, 598–618

    [ADLS07] David Angeli, Patrick De Leenheer, and Eduardo D Sontag,A petri net approach to the study of persistence in chemical reaction networks, Mathematical biosciences210(2007), no. 2, 598–618. [And08] David F Anderson,Global asymptotic stability for a class of nonlinear chemical equations, SIAM Journal on Applied Mathematics68(2008), no. 5, 1464–1476. ...

    work page 2007

  3. [3]

    DOI 10.48550/arXiv.2508.15273, 2508.15273

    [BC10] Murad Banaji and Gheorghe Craciun,Graph-theoretic criteria for injectivity and unique equilibria in general chemical reaction systems, Advances in Applied Mathematics44(2010), no. 2, 168–184. [BLN20] Alex Blokhuis, David Lacoste, and Philippe Nghe,Universal motifs and the diversity of autocatalytic systems, Proceedings of the National Academy of Sc...

    work page arXiv 2010

  4. [4]

    [FW13] Elisenda Feliu and Carsten Wiuf,Simplifying biochemical models with intermediate species, J

    [FS25] Elisenda Feliu and Anne Shiu,From chemical reaction networks to algebraic and polyhedral geometry–and back again, arXiv preprint arXiv:2501.06354 (2025). [FW13] Elisenda Feliu and Carsten Wiuf,Simplifying biochemical models with intermediate species, J. R. Soc. Interface10 (2013), no. 87, 20130484. [GGSV25] Richard Golnik, Thomas Gatter, Peter F. S...

    work page arXiv 2025

  5. [5]

    19 [HC11] Wassim M Haddad and VijaySekhar Chellaboina,Nonlinear dynamical systems and control: a lyapunov-based approach, Princeton University Press,

    [GLRS+23] Awildo Gutierrez, Elijah Leake, Caelyn Rivas-Sobie, Jordy Lopez Garcia, and Anne Shiu,Deficiency of chemical reaction networks: the effect of operations that preserve multistationarity and periodic orbits, arXiv preprint arXiv:2305.19410 (2023). 19 [HC11] Wassim M Haddad and VijaySekhar Chellaboina,Nonlinear dynamical systems and control: a lyap...

    work page arXiv 2023

  6. [6]

    limit sets, SIAM Journal on Mathematical Analysis13(1982), no

    [Hir82] Morris W Hirsch,Systems of differential equations which are competitive or cooperative: I. limit sets, SIAM Journal on Mathematical Analysis13(1982), no. 2, 167–179. [HSD13] Morris W. Hirsch, Stephen Smale, and Robert L. Devaney,Differential equations, dynamical systems, and an introduction to chaos, third ed., Elsevier/Academic Press, Amsterdam,

    work page 1982

  7. [7]

    6, 1271–1277

    MR 3293130 [IT79] AN Ivanova and BL Tarnopolskii,One approach to the determination of a number of qualitative features in the behavior of kinetic systems, and realization of this approach in a computer (critical conditions, autooscillations), Kinetics and Catalysis20(1979), no. 6, 1271–1277. [JA25] MD Johnston and F Avram,The boundary reproduction number ...

    work page 1979

  8. [8]

    12, 5500–5535

    [Soa18] Pedro Soares,The lifting bifurcation problem on feed-forward networks, Nonlinearity31(2018), no. 12, 5500–5535. [SS89] Lee A Segel and Marshall Slemrod,The quasi-steady-state assumption: a case study in perturbation, SIAM review 31(1989), no. 3, 446–477. [V AA24] Nicola Vassena, Florin Avram, and Rim Adenane,Finding bifurcations in mathematical ep...

    work page 2018

  9. [9]

    [Vas23] ,Symbolic hunt of instabilities and bifurcations in reaction networks, Discrete and Continuous Dynamical Systems-B (2023), 0–0

    [Vas19] Nicola Vassena,Good and bad children in metabolic networks, arXiv preprint arXiv:1905.12272 (2019). [Vas23] ,Symbolic hunt of instabilities and bifurcations in reaction networks, Discrete and Continuous Dynamical Systems-B (2023), 0–0. [Vas25] ,Mass action systems: two criteria for hopf bifurcation without hurwitz, SIAM Journal on Applied Mathe- m...

    work page arXiv 1905

  10. [10]

    2, 775–793

    MR 2004534 [WR04] Wendi Wang and Shigui Ruan,Bifurcations in an epidemic model with constant removal rate of the infectives, Journal of Mathematical Analysis and Applications291(2004), no. 2, 775–793. [XR07] Dongmei Xiao and Shigui Ruan,Global analysis of an epidemic model with nonmonotone incidence rate, Mathe- matical biosciences208(2007), no. 2, 419–42...

    work page 2004

  11. [11]

    1, 312–324

    [ZF12] Linhua Zhou and Meng Fan,Dynamics of an sir epidemic model with limited medical resources revisited, Nonlinear Analysis: Real World Applications13(2012), no. 1, 312–324. [ZL08] Xu Zhang and Xianning Liu,Backward bifurcation of an epidemic model with saturated treatment function, Journal of mathematical analysis and applications348(2008), no. 1, 433...

    work page 2012

  12. [12]

    5, 1903–1915

    [ZXL07] Yugui Zhou, Dongmei Xiao, and Yilong Li,Bifurcations of an epidemic model with non-monotonic incidence rate of saturated mass-action, Chaos, Solitons & Fractals32(2007), no. 5, 1903–1915. 21 22

    work page 2007