Solving Moment Hierarchies for Chemical Reaction Networks

Abstract

The study of Chemical Reaction Networks (CRN's) is a very active field. Earlier well-known results \cite{Feinberg:def_01:87, Anderson:product_dist:10} identify a topological quantity called deficiency, for any CRN, which, when exactly equal to zero, leads to a unique factorized steady-state for these networks. No results exist however for the steady states of non-zero-deficiency networks. In this paper, we show how to write the full moment-hierarchy for any non-zero-deficiency CRN obeying mass-action kinetics, in terms of equations for the factorial moments (FM). Using these, we can recursively predict values for lower moments from higher moments, reversing the procedure usually used to solve moment hierarchies. We show, for non-trivial examples, that in this manner we can predict to high accuracy, any moment of interest, for CRN's with non-zero deficiency and non-factorizable steady states.

version appearing in journal

References in corpus (3)

• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Comparison of different moment-closure approximations for stochastic chemical kinetics
• Dissipation in noisy chemical networks: The role of deficiency

Cited by in corpus (5)

• Intrinsic and extrinsic thermodynamics for stochastic population processes with multi-level large-deviation structure
• Flows, scaling, and the control of moment hierarchies for stochastic chemical reaction networks
• Action Functional Gradient Descent algorithm for estimating escape paths in Stochastic Chemical Reaction Networks
• Superintegrability and Deformed Oscillator Realizations of Quantum TTW Hamiltonians on Constant-Curvature Manifolds and with Reflections in a Plane
• Eikonal solutions for moment hierarchies of Chemical Reaction Networks in the limits of large particle number