Abstract:Living systems are usually maintained out of equilibrium and exhibit complex dynamical behaviors. The external energy supply often comes from chemical fluxes that can keep some species concentrations constant. Furthermore, the properties of the underlying chemical reaction networks (CRNs) are also instrumental in establishing robust biological functioning. Hence, capturing the emergent complexity of living systems and the role of their non-equilibrium nature is fundamental to uncover constraints and properties of the CRNs underpinning their functions. In particular, while kinetics plays a key role in shaping detailed dynamical phenomena, the range of operations of any CRN must be fundamentally constrained by thermodynamics, as they necessarily operate with a given energy budget. Here, we derive universal thermodynamic upper and lower bounds for the accessible space of species concentrations in a generic CRN. The resulting region determines the "thermodynamic space" of the CRN, a concept we introduce in this work. Moreover, we obtain similar bounds also for the affinities, shedding light on how global thermodynamic properties can limit local non-equilibrium quantities. We illustrate our results in two paradigmatic examples, the Schlögl model for bistability and a minimal self-assembly process, demonstrating how the onset of complex behaviors is intimately tangled with the presence of non-equilibrium driving. In summary, our work unveils the exact form of the accessible space in which a CRN must work as a function of its energy budget, shedding light on the non-equilibrium origin of a variety of phenomena, from amplification to pattern formation. Ultimately, by providing a general tool for analyzing CRNs, the presented framework constitutes a stepping stone to deepen our ability to predict complex out-of-equilibrium behaviors and design artificial chemical reaction systems.

On Persistence of Some $\mathcal {W}_{i}$-Endotactic Chemical Reaction Networks.

On Persistence of Some <i>W<sub>I</sub></i>-Endotactic Chemical Reaction Networks

Revisiting persistence of Chemical Reaction Network through Lyapunov PDE

Revisiting Persistence of Chemical Reaction Networks Through Lyapunov Function Partial Differential Equations

A computational approach to persistence, permanence, and endotacticity of biochemical reaction systems

Persistence of Delayed Complex Balanced Chemical Reaction Networks

Intermediates, Catalysts, Persistence, and Boundary Steady States

Autocatalysis in Reaction Networks

Emergence of Heterogeneous Structures in Chemical Reaction-Diffusion Networks.

Endotactic and strongly endotactic networks with infinitely many positive steady states

Global stability of first order endotactic reaction systems

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

On Stability of Two Kinds of Delayed Chemical Reaction Networks*

Capturing Persistence of High-Dimensional Delayed Complex Balanced Chemical Reaction Systems Via Decomposition of Semilocking Sets

Capturing persistence of delayed complex balanced chemical reaction systems via decomposition of semilocking sets

On the steady states of weakly reversible chemical reaction networks

Non-explosivity of endotactic stochastic reaction systems

Thermodynamic Space of Chemical Reaction Networks

Source-only realizations, weakly reversible deficiency one networks and dynamical equivalence

Robust Persistence and Permanence of Polynomial and Power Law Dynamical Systems