Nonequilibrium thermodynamic formalism of nonlinear chemical reaction systems with Waage–Guldberg’s law of mass action

Thermodynamic bounds on time-reversal asymmetry

Quantifying irreversibility of a system using finite information constitutes a major challenge in stochastic thermodynamics. We introduce an observable that measures the time-reversal asymmetry between two states after a given time lag. Our central result is a bound on the time-reversal asymmetry in terms of the total cycle affinity driving the system out of equilibrium. This result leads to further thermodynamic bounds on the asymmetry of directed fluxes, on the asymmetry of finite-time cross-correlations, and on the cycle affinity of coarse-grained dynamics.

Information-geometric structure for chemical thermodynamics: An explicit construction of dual affine coordinates

We construct an information-geometric structure for chemical thermodynamics, applicable to a wide range of chemical reaction systems including nonideal and open systems. For this purpose, we explicitly construct dual affine coordinate systems, which completely designate an information-geometric structure, using the extent of reactions and the affinities of reactions as coordinates on a linearly constrained space of amounts of substances. The resulting structure induces a metric and a divergence (a function of two distributions of amounts), both expressed with chemical potentials. These quantities have been partially known for ideal-dilute solutions, but their extensions for nonideal solutions and the complete underlying structure are novel. The constructed geometry is a generalization of dual affine coordinates for stochastic thermodynamics. For example, the metric and the divergence are generalizations of the Fisher information and the Kullback-Leibler divergence. As an application, we identify the chemical-thermodynamic analog of the Hatano-Sasa excess entropy production using our divergence.

Thermodynamics of concentration vs flux control in chemical reaction networks

We investigate the thermodynamic implications of two control mechanisms of open chemical reaction networks. The first controls the concentrations of the species that are exchanged with the surroundings, while the other controls the exchange fluxes. We show that the two mechanisms can be mapped one into the other and that the thermodynamic theories usually developed in the framework of concentration control can be applied to flux control as well. This implies that the thermodynamic potential and the fundamental forces driving chemical reaction networks out of equilibrium can be identified in the same way for both mechanisms. By analyzing the dynamics and thermodynamics of a simple enzymatic model, we also show that while the two mechanisms are equivalent at steady state, the flux control may lead to fundamentally different regimes where systems achieve stationary growth.

Recent Advances in Conservation-Dissipation Formalism for Irreversible Processes

The main purpose of this review is to summarize the recent advances of the Conservation–Dissipation Formalism (CDF), a new way for constructing both thermodynamically compatible and mathematically stable and well-posed models for irreversible processes. The contents include but are not restricted to the CDF’s physical motivations, mathematical foundations, formulations of several classical models in mathematical physics from master equations and Fokker–Planck equations to Boltzmann equations and quasi-linear Maxwell equations, as well as novel applications in the fields of non-Fourier heat conduction, non-Newtonian viscoelastic fluids, wave propagation/transportation in geophysics and neural science, soft matter physics, etc. Connections with other popular theories in the field of non-equilibrium thermodynamics are examined too.

Thermodynamic Uncertainty Relation and Thermodynamic Speed Limit in Deterministic Chemical Reaction Networks

We generalize the thermodynamic uncertainty relation (TUR) and thermodynamic speed limit (TSL) for deterministic chemical reaction networks (CRNs). The scaled diffusion coefficient derived by considering the connection between macro- and mesoscopic CRNs plays an essential role in our results. The TUR shows that the product of the entropy production rate and the ratio of the scaled diffusion coefficient to the square of the rate of concentration change is bounded below by two. The TSL states a trade-off relation between speed and thermodynamic quantities, the entropy production, and the time-averaged scaled diffusion coefficient. The results are proved under the general setting of open and nonideal CRNs.

Nonequilibrium thermodynamics of light-induced reactions

Current formulations of nonequilibrium thermodynamics of open chemical reaction networks only consider chemostats as free-energy sources sustaining nonequilibrium behaviors. Here, we extend the theory to include incoherent light as a source of free energy. We do so by relying on a local equilibrium assumption to derive the chemical potential of photons relative to the system they interact with. This allows us to identify the thermodynamic potential and the thermodynamic forces driving light-reacting chemical systems out-of-equilibrium. We use this framework to treat two paradigmatic photochemical mechanisms describing light-induced unimolecular reactions-namely, the adiabatic and diabatic mechanisms-and highlight the different thermodynamics they lead to. Furthermore, using a thermodynamic coarse-graining procedure, we express our findings in terms of commonly measured experimental quantities, such as quantum yields.

Enzyme activities predicted by metabolite concentrations and solvent capacity in the cell

Experimental measurements or computational model predictions of the post-translational regulation of enzymes needed in a metabolic pathway is a difficult problem. Consequently, regulation is mostly known only for well-studied reactions of central metabolism in various model organisms. In this study, we use two approaches to predict enzyme regulation policies and investigate the hypothesis that regulation is driven by the need to maintain the solvent capacity in the cell. The first predictive method uses a statistical thermodynamics and metabolic control theory framework while the second method is performed using a hybrid optimization-reinforcement learning approach. Efficient regulation schemes were learned from experimental data that either agree with theoretical calculations or result in a higher cell fitness using maximum useful work as a metric. As previously hypothesized, regulation is herein shown to control the concentrations of both immediate and downstream product concentrations at physiological levels. Model predictions provide the following two novel general principles: (1) the regulation itself causes the reactions to be much further from equilibrium instead of the common assumption that highly non-equilibrium reactions are the targets for regulation; and (2) the minimal regulation needed to maintain metabolite levels at physiological concentrations maximizes the free energy dissipation rate instead of preserving a specific energy charge. The resulting energy dissipation rate is an emergent property of regulation which may be represented by a high value of the adenylate energy charge. In addition, the predictions demonstrate that the amount of regulation needed can be minimized if it is applied at the beginning or branch point of a pathway, in agreement with common notions. The approach is demonstrated for three pathways in the central metabolism of E. coli (gluconeogenesis, glycolysis-tricarboxylic acid (TCA) and pentose phosphate-TCA) that each require different regulation schemes. It is shown quantitatively that hexokinase, glucose 6-phosphate dehydrogenase and glyceraldehyde phosphate dehydrogenase, all branch points of pathways, play the largest roles in regulating central metabolism.

Stochastic thermodynamics and entropy production of chemical reaction systems

We investigate the nonequilibrium stationary states of systems consisting of chemical reactions among molecules of several chemical species. To this end, we introduce and develop a stochastic formulation of nonequilibrium thermodynamics of chemical reaction systems based on a master equation defined on the space of microscopic chemical states and on appropriate definitions of entropy and entropy production. The system is in contact with a heat reservoir and is placed out of equilibrium by the contact with particle reservoirs. In our approach, the fluxes of various types, such as the heat and particle fluxes, play a fundamental role in characterizing the nonequilibrium chemical state. We show that the rate of entropy production in the stationary nonequilibrium state is a bilinear form in the affinities and the fluxes of reaction, which are expressed in terms of rate constants and transition rates, respectively. We also show how the description in terms of microscopic states can be reduced to a description in terms of the numbers of particles of each species, from which follows the chemical master equation. As an example, we calculate the rate of entropy production of the first and second Schlögl reaction models.

Thermodynamic Analysis of Chemically Reacting Mixtures-Comparison of First and Second Order Models

Recently, a method based on non-equilibrium continuum thermodynamics which derives thermodynamically consistent reaction rate models together with thermodynamic constraints on their parameters was analyzed using a triangular reaction scheme. The scheme was kinetically of the first order. Here, the analysis is further developed for several first and second order schemes to gain a deeper insight into the thermodynamic consistency of rate equations and relationships between chemical thermodynamic and kinetics. It is shown that the thermodynamic constraints on the so-called proper rate coefficient are usually simple sign restrictions consistent with the supposed reaction directions. Constraints on the so-called coupling rate coefficients are more complex and weaker. This means more freedom in kinetic coupling between reaction steps in a scheme, i.e., in the kinetic effects of other reactions on the rate of some reaction in a reacting system. When compared with traditional mass-action rate equations, the method allows a reduction in the number of traditional rate constants to be evaluated from data, i.e., a reduction in the dimensionality of the parameter estimation problem. This is due to identifying relationships between mass-action rate constants (relationships which also include thermodynamic equilibrium constants) which have so far been unknown.

Mesoscopic kinetic basis of macroscopic chemical thermodynamics: A mathematical theory

Gibbs' macroscopic chemical thermodynamics is one of the most important theories in chemistry. Generalizing it to mesoscaled nonequilibrium systems is essential to biophysics. The nonequilibrium stochastic thermodynamics of chemical reaction kinetics suggested a free energy balance equation dF^{(meso)}/dt=E_{in}-e_{p} in which the free energy input rate E_{in} and dissipation rate e_{p} are both non-negative, and E_{in}≤e_{p}. We prove that in the macroscopic limit by merely allowing the molecular numbers to be infinite, the generalized mesoscopic free energy F^{(meso)} converges to φ^{ss}, the large deviation rate function for the stationary distributions. This generalized macroscopic free energy φ^{ss} now satisfies a balance equation dφ^{ss}(x)/dt=cmf(x)-σ(x), in which x represents chemical concentration. The chemical motive force cmf(x) and entropy production rate σ(x) are both non-negative, and cmf(x)≤σ(x). The balance equation is valid generally in isothermal driven systems and is different from mechanical energy conservation and the first law; it is actually an unknown form of the second law. Consequences of the emergent thermodynamic quantities and equalities are further discussed. The emergent "law" is independent of underlying kinetic details. Our theory provides an example showing how a macroscopic law emerges from a level below.