Characterizing the conditions for indefinite growth in open chemical reaction networks.

The thermodynamic and dynamical conditions necessary to observe indefinite growth in homogeneous open chemical reaction networks (CRNs) satisfying mass action kinetics are presented in Srinivas et al. [Phys. Rev. Lett. 132, 268001 (2024)10.1103/PhysRevLett.132.268001]. Unimolecular CRNs can accumulate only equilibrium concentrations of species while multimolecular CRNs are needed to produce indefinite growth with nonequilibrium concentrations. Within multimolecular CRNs, pseudo-unimolecular CRNs produce nonequilibrium concentrations with zero efficiencies. Nonequilibrium growth with efficiencies greater than zero requires dynamically nonlinear CRNs. In this paper, we provide a detailed analysis supporting these results. Mathematical proofs are provided for growth in unimolecular and pseudo-unimolecular CRNs. For multimolecular CRNs, four models displaying very distinctive topological properties are extensively studied, both numerically and partly analytically.

Thermodynamics of Growth in Open Chemical Reaction Networks.

We identify the thermodynamic conditions necessary to observe indefinite growth in homogeneous open chemical reaction networks (CRNs) satisfying mass action kinetics. We also characterize the thermodynamic efficiency of growth by considering the fraction of the chemical work supplied from the surroundings that is converted into CRN free energy. We find that indefinite growth cannot arise in CRNs chemostatted by fixing the concentration of some species at constant values, or in continuous-flow stirred tank reactors. Indefinite growth requires a constant net influx from the surroundings of at least one species. In this case, unimolecular CRNs always generate equilibrium linear growth, i.e., a continuous linear accumulation of species with equilibrium concentrations and efficiency one. Multimolecular CRNs are necessary to generate nonequilibrium growth, i.e., the continuous accumulation of species with nonequilibrium concentrations. Pseudounimolecular CRNs-a subclass of multimolecular CRNs-always generate asymptotic linear growth with zero efficiency. Our findings demonstrate the importance of the CRN topology and the chemostatting procedure in determining the dynamics and thermodynamics of growth.

Deficiency, kinetic invertibility, and catalysis in stochastic chemical reaction networks.

Stochastic chemical processes are described by the chemical master equation satisfying the law of mass-action. We first ask whether the dual master equation, which has the same steady state as the chemical master equation, but with inverted reaction currents, satisfies the law of mass-action and, hence, still describes a chemical process. We prove that the answer depends on the topological property of the underlying chemical reaction network known as deficiency. The answer is yes only for deficiency-zero networks. It is no for all other networks, implying that their steady-state currents cannot be inverted by controlling the kinetic constants of the reactions. Hence, the network deficiency imposes a form of non-invertibility to the chemical dynamics. We then ask whether catalytic chemical networks are deficiency-zero. We prove that the answer is no when they are driven out of equilibrium due to the exchange of some species with the environment.