Boundedness and persistence of weakly reversible mass-action systemsResearch Paper
Motivation: can a weakly reversible reaction network lose a species or blow up?
Chemical reaction network theory studies the polynomial differential equations that describe concentrations of interacting species under mass-action kinetics, where each reaction proceeds at a rate proportional to the product of the concentrations of its reactants. A central theme is which properties of the dynamics are forced by the structure of the reaction graph alone, independently of the (usually unknown) rate constants. Weak reversibility — every reaction can be undone by some directed path of reactions — is the main structural condition in this theory. Two long-standing conjectures ask whether weak reversibility alone guarantees that no species dies out (persistence) and that no concentration grows without bound (boundedness). These questions matter for systems biology and chemical engineering, where extinction of a species or unbounded growth would contradict the modeled behavior, and they are closely tied to the global attractor conjecture for complex-balanced systems.
Timeline
- 1972 — Horn and Jackson develop the equilibrium and stability theory of complex-balanced mass-action systems (ARMA 1972); Horn (ARMA 1972) and Feinberg (ARMA 1972) give complex-balancing criteria and the deficiency-zero analysis.
- 1987 — Feinberg formulates the expectation that a positive trajectory of a weakly reversible system cannot converge to a boundary point (Chem. Eng. Sci. 1987, Remark 6.1.E).
- 2010 — August and Barahona state the general boundedness and persistence conclusion (IFAC 2010); the source identifies a gap in one step of their proof (p. 2).
- 2011 — Anderson separates the boundedness conjecture from the persistence conjecture and proves boundedness for a single linkage class (J. Math. Chem. 2011); he also proves the global attractor conjecture in the single-linkage case (SIAM J. Appl. Math. 2011).
- 2012–2013 — Pantea proves persistence of bounded trajectories when the stoichiometric subspace has dimension two (SIAM J. Math. Anal. 2012); Craciun, Nazarov and Pantea prove boundedness, persistence and permanence for two-species systems (SIAM J. Appl. Math. 2013).
- 2014 — Gopalkrishnan, Miller and Shiu prove permanence for strongly endotactic systems (SIAM J. Appl. Dyn. Syst. 2014).
- 2019–2020 — Boros proves every positive stoichiometric class of a weakly reversible system contains a positive equilibrium (SIAM J. Math. Anal. 2019); Boros and Hofbauer reprove single-linkage permanence (SIAM J. Appl. Dyn. Syst. 2020).
- 2026 — Craciun's revised toric-differential-inclusion approach gives persistence for bounded trajectories (arXiv:1501.02860v3).
- 2026 — The OpenAI preprint Boundedness and persistence of weakly reversible mass-action systems (OpenAI Math Release, September 25, 2026) claims both conjectures for constant positive rates in every dimension. It has not been peer reviewed and its theorem is not formally verified.
Setting
Fix species. A reaction network consists of a finite set of complexes and a set of reactions with , . It is weakly reversible if for every reaction there is a directed path of reactions from back to . Given positive rate constants , the mass-action system is
A positive trajectory is persistent if for every .
In Lean (namespace OAI.Problem326), a ReactionNetwork d has a Finset of complexes in Fin d → ℕ and a Finset of reactions between them with distinct source and target; WeaklyReversible is Relation.TransGen reachability from target back to source; massAction N κ is the vector field above; and IsGlobalForwardSolution N κ x0 x means and has derivative at every .
Formalization targets
Goal: global existence, boundedness and persistence (Theorem 1.1)
Let , let the network be weakly reversible and all rates positive. For every the solution with exists for all , and there is , depending on the network, the rates and , such that
The goal is published on the platform with status Open.
Significance
The result itself. The theorem proves the boundedness conjecture and the persistence conjecture together, for arbitrary constant positive rates, any number of linkage classes and any number of species; earlier results needed a single linkage class, two species, or a two-dimensional stoichiometric subspace. It applies whether or not the stoichiometric compatibility class is bounded. Combined with complex balance, it yields convergence to the unique positive equilibrium in each class (the global attractor conclusion; Remark 3.5, p. 15). The proof also produces, for each positive initial point, a compact convex forward-invariant polytope in the open orthant (Proposition 3.3, p. 13), which gives a short proof of Boros's existence theorem for positive equilibria (Corollary 3.4, p. 14).
Formalizing it. Statements about reaction networks quantify over all networks and all rates, so they are a natural fit for formal verification; the history includes a published proof with a gap in one step. A machine-checked proof would settle the status of the general theorem independently of that history.
Difficulty
The natural approach is to find a Lyapunov function or an invariant region whose boundary each reaction crosses inward. Individual reactions need not point inward across any fixed face, and near the orthant boundary and at infinity different monomials dominate in different regions, with ties between monomials that an arbitrary choice of face normal can reverse. The argument must control the combined flux across cuts of the reaction graph using the return paths given by weak reversibility, and must do so simultaneously for small and large concentrations; a naive comparison of monomials by total degree fails, as the source's example with shows.
Formalization scope
- Complexes are vectors in
Fin d → ℕ, so exponents are nonnegative integers and is a polynomial vector field on ; the reaction set is aFinsetwith source ≠ target. Empty reaction sets and isolated complexes are allowed, as in the source. - Rates are a function on reactions with
∀ e, 0 < κ e; is a hypothesis. - The conclusion asserts existence of some global forward solution and bounds for every global forward solution from (uniqueness is therefore not assumed). Solutions are functions
ℝ → (Fin d → ℝ)withHasDerivAtat every , including a two-sided derivative at . - The same works for all and all species; it may depend on . No uniformity over initial points is claimed.
- Needed infrastructure: Picard–Lindelöf for polynomial fields, forward invariance of convex polytopes (Nagumo-type tangency conditions), and graph-cut flux estimates. ODE invariance tools are reusable beyond this mission.
Selected references
- F. Horn and R. Jackson, General mass action kinetics, Arch. Rational Mech. Anal. 47 (1972). https://doi.org/10.1007/BF00251225
- M. Feinberg, Complex balancing in general kinetic systems, Arch. Rational Mech. Anal. 49 (1972). https://doi.org/10.1007/BF00255665
- M. Feinberg, Chemical reaction network structure and the stability of complex isothermal reactors — I, Chem. Eng. Sci. 42 (1987). https://doi.org/10.1016/0009-2509(87)80099-4
- E. August and M. Barahona, Solutions of weakly reversible chemical reaction networks are bounded and persistent, IFAC Proc. Vol. 43 (2010). https://doi.org/10.3182/20100707-3-BE-2012.0018
- D. F. Anderson, Boundedness of trajectories for weakly reversible, single linkage class reaction systems, J. Math. Chem. 49 (2011). https://doi.org/10.1007/s10910-011-9886-4
- C. Pantea, On the persistence and global stability of mass-action systems, SIAM J. Math. Anal. 44 (2012). https://doi.org/10.1137/110840509
- G. Craciun, F. Nazarov and C. Pantea, Persistence and permanence of mass-action and power-law dynamical systems, SIAM J. Appl. Math. 73 (2013). https://doi.org/10.1137/100812355
- M. Gopalkrishnan, E. Miller and A. Shiu, A geometric approach to the global attractor conjecture, SIAM J. Appl. Dyn. Syst. 13 (2014). https://doi.org/10.1137/130928170
- B. Boros, Existence of positive steady states for weakly reversible mass-action systems, SIAM J. Math. Anal. 51 (2019). https://doi.org/10.1137/17M115534X
- G. Craciun, Toric differential inclusions and a proof of the global attractor conjecture, arXiv:1501.02860v3, 2026. https://arxiv.org/abs/1501.02860v3
- OpenAI, Boundedness and persistence of weakly reversible mass-action systems, OpenAI Math Release preprint, September 25, 2026 (Theorem 1.1, p. 1). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Boundedness-and-persistence-of-weakly-reversible-mass-action-systems-September-25-2026/paper.pdf