Abstract:SIAM Journal on Applied Dynamical Systems, Volume 23, Issue 1, Page 256-296, March 2024. Model reduction of fast-slow chemical reaction networks based on the quasi-steady state approximation fails when the fast subsystem has first integrals. We call these first integrals approximate conservation laws. In order to define fast subsystems and identify approximate conservation laws, we use ideas from tropical geometry. We prove that any approximate conservation law evolves more slowly than all the species involved in it and therefore represents a supplementary slow variable in an extended system. By elimination of some variables of the extended system, we obtain networks without approximate conservation laws, which can be reduced by standard singular perturbation methods. The field of applications of approximate conservation laws covers the quasi-equilibrium approximation, which is well known in biochemistry. We discuss reductions of slow-fast as well as multiple timescale systems. Networks with multiple timescales have hierarchical relaxation. At a given timescale, our multiple timescale reduction method defines three subsystems composed of (i) slaved fast variables satisfying algebraic equations, (ii) slow driving variables satisfying reduced ordinary differential equations, and (iii) quenched much slower variables that are constant. The algebraic equations satisfied by fast variables define chains of nested normally hyperbolic invariant manifolds. In such chains, faster manifolds are of higher dimension and contain the slower manifolds. Our reduction methods are introduced algorithmically for networks with monomial reaction rates and linear, monomial, or polynomial approximate conservation laws. We propose symbolic algorithms to reshape and rescale the networks such that geometric singular perturbation theory can be applied to them, test the applicability of the theory, and finally reduce the networks. As a proof of concept, we apply this method to a model of the TGF-beta signaling pathway.

Conservation-dissipation Structure of Chemical Reaction Systems

Emergence of Heterogeneous Structures in Chemical Reaction-Diffusion Networks.

Universal limiting behaviour of reaction-diffusion systems with conservation laws

Trend to equilibrium for reaction-diffusion systems arising from complex balanced chemical reaction networks

Nonequilibrium Thermodynamic Formalism of Nonlinear Chemical Reaction Systems with Waage–Guldberg’s Law of Mass Action

Fine-tuning the basic forces of nature through the triple-alpha process in red giant stars

Explicit exponential convergence to equilibrium for nonlinear reaction-diffusion systems with detailed balance condition

Modeling of chemical reaction systems with detailed balance using gradient structures

Convergence to equilibrium of renormalised solutions to nonlinear chemical reaction-diffusion systems

Field Theory of Reaction-Diffusion: Mass Action with an Energetic Variational Approach

Mathematical Analysis of Chemical Reaction Systems

Reduction of Chemical Reaction Networks with Approximate Conservation Laws

On a positive-preserving, energy-stable numerical scheme to mass-action kinetics with detailed balance

Conservation-dissipation Formalism of Irreversible Thermodynamics

Discrete models for chemically reacting systems

Persistence of Delayed Complex Balanced Chemical Reaction Networks

On data and dimension in chemistry I -- irreversibility, concealment and emergent conservation laws

Universal slow dynamics of chemical reaction networks

Mathematical Formalism of Nonequilibrium Thermodynamics for Nonlinear Chemical Reaction Systems with General Rate Law

A network dynamics approach to chemical reaction networks