Abstract:In this paper, we provide a detailed theoretical analysis of the numerical scheme introduced in J. Comput. Phys. 436 (2021) 110253 for the reaction kinetics of a class of chemical reaction networks that satisfies detailed balance condition. In contrast to conventional numerical approximations, which are typically constructed based on ordinary differential equations (ODEs) for the concentrations of all involved species, the scheme is developed using the equations of reaction trajectories, which can be viewed as a generalized gradient flow of physically relevant free energy. The unique solvability, positivity-preserving, and energy-stable properties are proved for the general case involving multiple reactions, under a mild condition on the stoichiometric matrix.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop and analyze a numerical scheme for simulating the kinetics of chemical reaction networks that satisfy the detailed balance condition. Specifically, the author focuses on ensuring that the numerical scheme can preserve several key properties of the physical system: 1. **Positivity - preserving**: Ensure that all species concentrations remain non - negative throughout the numerical simulation process. 2. **Energy stability**: Ensure that the total energy (or free energy) of the system decreases monotonically with time, in accordance with the energy dissipation law of the physical system. 3. **Unique solvability**: Ensure that the numerical scheme can find a unique solution at each iteration. ### Background and Problem Description Traditional numerical methods are usually based on ordinary differential equations (ODEs) to describe the concentration changes of all participating species. However, these methods face challenges when dealing with complex chemical reaction networks, especially when multiple reactions and detailed balance conditions are involved. To overcome these problems, the author proposes a new numerical scheme, which is based on the reaction trajectory equation rather than directly solving the ODEs for concentration changes. ### Reaction Kinetic Equations Consider a chemical reaction network containing \( N \) species \( X_1, X_2,\ldots, X_N \) and \( M \) reversible reactions: \[ \alpha_{l1}X_1+\alpha_{l2}X_2+\cdots+\alpha_{lN}X_N\rightleftharpoons\beta_{l1}X_1+\beta_{l2}X_2+\cdots+\beta_{lN}X_N,\quad l = 1,\ldots,M \] where \( \alpha_{li} \) and \( \beta_{li} \) are the stoichiometric coefficients of the \( l\) - th reaction. The reaction rate \( r_l(c)\) is usually expressed as the difference between the forward and reverse reaction rates: \[ r_l(c)=r_l^+(c)-r_l^-(c) \] According to the Law of Mass Action, the reaction rate can be expressed as: \[ r_l^+(c)=k_l^+c^{\alpha_l},\quad r_l^-(c)=k_l^-c^{\beta_l} \] where \( c^{\alpha_l}=\prod_{i = 1}^N c_i^{\alpha_{li}}\), \( c^{\beta_l}=\prod_{i = 1}^N c_i^{\beta_{li}}\) ### Key Features of the Numerical Scheme The numerical scheme proposed by the author has the following key features: - **Reaction Trajectory Equation**: By introducing the reaction trajectory \( R(t)\in\mathbb{R}^M\), the concentration change is expressed as \( c(t)=c_0+SR(t)\), where \( S \) is the stoichiometric matrix. - **Generalized Gradient Flow**: The reaction kinetics can be regarded as the generalized gradient flow of the free energy \( F[c]\), satisfying the energy dissipation law: \[ \frac{d}{dt}F[c(R)]=-D(R,\dot{R}) \] - **Variational Structure**: Through the variational principle, the reaction trajectory equation is discretized into an optimization problem, ensuring the stability and uniqueness of the numerical scheme. ### Research Objectives The main objective of the paper is to prove that the proposed numerical scheme can still maintain positivity, energy stability, and unique solvability in the case of multiple reactions (i.e., \( M>1\)). This provides a theoretical basis for the efficient numerical simulation of complex chemical reaction networks. Through rigorous mathematical analysis, the author proves that, in the case of satisfying the detailed balance condition, this numerical scheme can maintain the above - mentioned key properties during multiple iterations, thus providing a reliable method for the numerical simulation of chemical reaction kinetics.

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

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

A Second-Order Accurate, Operator Splitting Scheme for Reaction-Diffusion Systems in an Energetic Variational Formulation

Existence of Positive Steady States for Mass Conserving and Mass-Action Chemical Reaction Networks with a Single Terminal-Linkage Class

Efficient bound preserving and asymptotic preserving semi-implicit schemes for the fast reaction-diffusion system

Convergence analysis of a positivity-preserving numerical scheme for the Cahn-Hilliard-Stokes system with Flory-Huggins energy potential

A second-order accurate, original energy dissipative numerical scheme for chemotaxis and its convergence analysis

Persistence of Delayed Complex Balanced Chemical Reaction Networks

A Decomposition-Based Approach for Deriving Positive Steady States of a Class of Chemical Reaction Networks with Non-Mass-Action Kinetics

A Second Order Accurate, Positivity Preserving Numerical Method for the Poisson–Nernst–Planck System and Its Convergence Analysis

Global stability of first order endotactic reaction systems

Second Order Unconditional Positive Preserving Schemes for Non-equilibrium Reactive Flows with Mass and Mole Balance

Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations

A Graphic Formulation of Nonisothermal Chemical Reaction Systems and the Analysis of Detailed Balanced Networks.

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case

An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling

A Kinetic Scheme Based On Positivity Preservation For Multi-component Euler Equations

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

A Balanced Shaping Method for Stabilizing Chemical Reaction Networks

Sufficient conditions for linear stability of complex-balanced equilibria in generalized mass-action systems