Stochastic approach to entropy production in chemical chaos

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.

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.

Slow-fast oscillatory dynamics and phantom attractors in stochastic modeling of biochemical reactions

A problem of the probabilistic analysis of stochastic phenomena in slow-fast dynamical systems modeling biochemical reactions is considered. We study how multiplicative noise induces systematic shifts of probabilistic distributions and forms "phantom" attractors in nonlinear enzymatic models. The mathematical analysis of the underlying probabilistic mechanism of such stochastic transformations is performed by the "freeze-and-average" method. Our theoretical results are supported by direct numerical simulation.

Nonequilibrium thermodynamics of non-ideal chemical reaction networks

All current formulations of nonequilibrium thermodynamics of open chemical reaction networks rely on the assumption of non-interacting species. We develop a general theory that accounts for interactions between chemical species within a mean-field approach using activity coefficients. Thermodynamic consistency requires that rate equations do not obey standard mass-action kinetics but account for the interactions with concentration dependent kinetic constants. Many features of the ideal formulations are recovered. Crucially, the thermodynamic potential and the forces driving non-ideal chemical systems out of equilibrium are identified. Our theory is general and holds for any mean-field expression of the interactions leading to lower bounded free energies.