Autocatalysis and Boundary Stability in Chemical Reaction Systems
Abstract
The next-generation matrix method is a powerful tool for computing the basic reproduction number in compartmental mathematical models of infectious diseases. The method has been recently extended to mathematical biochemistry, where it can be used to establish parameter regions of stability and instability for boundary steady states. Several significant challenges in its application remain, however, particularly around establishing conditions under which the method is mathematically valid, computationally tractable, and biologically meaningful in the biochemical setting. In this paper, we address these challenges by shifting the interpretation from new infections in the epidemiological setting to autocatalysis in the biochemical setting. We introduce a graph, called the ALT-graph (autocatalysis-leak-transition graph), which decomposes the contribution of each reaction to the Jacobian as autocatalytic, leak, or transition edges. We then present a systematic algorithm for splitting the ALT-graph, which is guaranteed to produce a valid reproduction number, ρ(FV−1), while also decreasing computational complexity by lowering the rank of FV−1. We apply the method to models of both biochemical reaction networks and infectious disease spread.
Related
- Classification of strongly non-hyperbolic critical points of planar second-order chemical reaction systems
- Large-deviations theory for growing chemical reaction networks
- Identifying unknown time delay and spatially varying coefficients in a reaction-diffusion equation from boundary measurements
- Reaction-aware hypergraph learning of chemical representations with LARK
- Generative discovery of chemical reaction networks
BibTeX
@misc{johnston2026autocatalysis,
title = {{Autocatalysis and Boundary Stability in Chemical Reaction Systems}},
author = {Matthew D. Johnston and Pol Jarne Cupons},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.01636},
eprint = {2610.01636},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.01636},
}
RIS
TY - UNPB
TI - Autocatalysis and Boundary Stability in Chemical Reaction Systems
AU - Matthew D. Johnston
AU - Pol Jarne Cupons
PY - 2026
DA - 2026/10/01
PB - arXiv
DO - 10.48550/arxiv.2610.01636
UR - https://arxiv.org/abs/2610.01636
AB - The next-generation matrix method is a powerful tool for computing the basic reproduction number in compartmental mathematical models of infectious diseases. The method has been recently extended to mathematical biochemistry, where it can be used to establish parameter regions of stability and instability for boundary steady states. Several significant challenges in its application remain, however, particularly around establishing conditions under which the method is mathematically valid, computationally tractable, and biologically meaningful in the biochemical setting. In this paper, we address these challenges by shifting the interpretation from new infections in the epidemiological setting to autocatalysis in the biochemical setting. We introduce a graph, called the ALT-graph (autocatalysis-leak-transition graph), which decomposes the contribution of each reaction to the Jacobian as autocatalytic, leak, or transition edges. We then present a systematic algorithm for splitting the ALT-graph, which is guaranteed to produce a valid reproduction number, ρ(FV^−1), while also decreasing computational complexity by lowering the rank of FV^−1. We apply the method to models of both biochemical reaction networks and infectious disease spread.
ER -