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

Trade-off between coherence and dissipation for excitable phase oscillators

Thermodynamic uncertainty relation (TUR) bounds coherence in stochastic oscillatory systems. In this paper, we show that both dynamical and thermodynamic bounds play important roles for the excitable oscillators, e.g., neurons. Firstly, we investigate the trade-off between coherence and dissipation both in the sub- and superthreshold regions for a single excitable unit, where both the TUR and the saddle-node on an invariant circle (SNIC) bounds constrain the fluctuation of interspike intervals. Secondly, we show that the widely studied phenomenon called coherence resonance, where there exists a noise strength to make the oscillatory responses of the system most coherent, is also bounded by the TUR in the one-dimensional excitable phase model. Finally, we study the coherence-dissipation relation in ensembles of strongly coupled excitable oscillators.

Thermodynamic Bounds on Generalized Transport: From Single-Molecule to Bulk Observables

We prove that the transport of any differentiable scalar observable in d-dimensional nonequilibrium systems is bounded from above by the total entropy production scaled by the amount the observation "stretches" microscopic coordinates. The result-a time-integrated generalized speed limit-reflects the thermodynamic cost of transport of observables, and places underdamped and overdamped stochastic dynamics on equal footing with deterministic motion. Our work allows for stochastic thermodynamics to make contact with bulk experiments, and fills an important gap in thermodynamic inference, since microscopic dynamics is, at least for short times, underdamped. Requiring only averages but not sample-to-sample fluctuations, the proven transport bound is practical and applicable not only to single-molecule but also bulk experiments where only averages are observed, which we demonstrate by examples. Our results may facilitate thermodynamic inference on molecular machines without an obvious directionality from bulk observations of transients probed, e.g., in time-resolved x-ray scattering.

Stochastic differential equation for a system coupled to a thermostatic bath via an arbitrary interaction Hamiltonian

The conventional Langevin equation offers a mathematically convenient framework for investigating open stochastic systems interacting with their environment or a bath. However, it is not suitable for a wide variety of systems whose dynamics rely on the nature of the environmental interaction, as the equation does not incorporate any specific information regarding that interaction. Here, we present a stochastic differential equation (SDE) for an open system coupled to a thermostatic bath via an arbitrary interaction Hamiltonian. This SDE encodes the interaction information to a fictitious potential (mean force) and a position-dependent damping coefficient. Surprisingly, we find that the conventional Langevin equation can be recovered in the presence of arbitrary strong interactions given two conditions: translational invariance of the potential and mutual independence of baths. Our results provide a comprehensive framework for studying open stochastic systems with an arbitrary interaction Hamiltonian and yield deeper insight into why various experiments fit the conventional Langevin description regardless of the strength or type of interaction.

Role of interactions in nonequilibrium transformations

For arbitrary nonequilibrium transformations in complex systems, we show that the distance between the current state and a target state can be decomposed into two terms: one corresponding to an independent estimate of the distance, and another corresponding to interactions, quantified using the relative mutual information between the variables. This decomposition is a special case of a more general decomposition involving successive orders of correlation or interactions among the degrees of freedom of the system. To illustrate its practical significance, we study the thermal relaxation of two interacting, optically trapped colloidal particles, where increasing pairwise interaction strength is shown to prolong the longevity of the time-dependent nonequilibrium state. Additionally, we study a system with both pairwise and triplet interactions, where our approach identifies their distinct contributions to the transformation. In more general setups where it is possible to control the strength of different orders of interactions, our findings provide a way to disentangle their effects and identify interactions that facilitate the transformation.

Information geometric bound on general chemical reaction networks

We investigate the convergence of chemical reaction networks (CRNs), aiming to establish an upper bound on their reaction rates. The nonlinear characteristics and discrete composition of CRNs pose significant challenges in this endeavor. To circumvent these complexities, we adopt an information geometric perspective, utilizing the natural gradient to formulate a nonlinear system. This system effectively determines an upper bound for the dynamics of CRNs. We corroborate our methodology through numerical simulations, which reveal that our constructed system converges more rapidly than CRNs within a particular class of reactions. This class is defined by the count of chemicals, the highest stoichiometric coefficients in the reactions, and the total number of reactions involved. Further, we juxtapose our approach with traditional methods, illustrating that the latter falls short in providing an upper bound for CRN reaction rates. Although our investigation centers on CRNs, the widespread presence of hypergraphs across various disciplines, ranging from natural sciences to engineering, indicates potential wider applications of our method, including in the realm of information science.

Persistence time bound for subdiffusion based on multidimensional thermodynamic uncertainty relation: Application to an analytically solvable model

The thermodynamic uncertainty relation (TUR) is an inequality showing the tradeoff relationship between the relative fluctuation of current observables and thermodynamic costs. It is one of the most important results of stochastic thermodynamics. There are various applications for TUR, one of which is the recent finding of thermodynamic constraints on the time window in which anomalous diffusion of Brownian particles can occur, including subdiffusion and superdiffusion, which are slower and faster than normal diffusion, respectively. These constraints are quite nontrivial because they are not generally derived from the asymptotic normal-diffusive behavior of the anomalous diffusion itself. In this study, we applied multidimensional TUR to the subdiffusion of Brownian particles obeying multivariate Langevin dynamics with a translationally invariant Hamiltonian in equilibrium. Multidimensional TUR is an improved TUR that includes information on another observable in addition to the one currently being considered. The use of an additional observable yields tighter bounds on the current fluctuation than those obtained using TUR. As an example, we demonstrated our theory using the one-dimensional Rouse model, which is known as a simple and analytically tractable model of polymer chains. We demonstrated that we improved the bounds for the persistence time of subdiffusion of the Rouse model, which became tighter as a more correlated observable with the current was used.

Improving the Cramér-Rao bound with the detailed fluctuation theorem

In some nonequilibrium systems, the distribution of entropy production p(Σ) satisfies the detailed fluctuation theorem (DFT) p(Σ)/p(-Σ)=exp(Σ). When the distribution p(Σ) shows a time dependence, the celebrated Cramér-Rao (CR) bound asserts that the mean entropy production rate is upper bounded in terms of the variance of Σ and the Fisher information with respect to time. In this paper we employ the DFT to derive an upper bound for the mean entropy production rate that improves the CR bound. We show that this new bound serves as an accurate approximation for the entropy production rate in the heat exchange problem mediated by a weakly coupled bosonic mode. The bound is saturated for the same setup when mediated by a weakly coupled qubit.

Thermodynamic Bounds on Correlation Times

We derive a variational expression for the correlation time of physical observables in steady-state diffusive systems. As a consequence of this variational expression, we obtain lower bounds on the correlation time, which provide speed limits on the self-averaging of observables. In equilibrium, the bound takes the form of a trade-off relation between the long- and short-time fluctuations of an observable. Out of equilibrium, the trade-off can be violated, leading to an acceleration of self-averaging. We relate this violation to the steady-state entropy production rate, as well as the geometric structure of the irreversible currents, giving rise to two complementary speed limits. One of these can be formulated as a lower estimate on the entropy production from the measurement of time-symmetric observables. Using an illustrating example, we show the intricate behavior of the correlation time out of equilibrium for different classes of observables and how this can be used to partially infer dissipation even if no time-reversal symmetry breaking can be observed in the trajectories of the observable.

The fundamental theorem of natural selection in optimization and games

The Fisher's fundamental theorem of natural selection (FTNS) is a matter of longstanding debate in the community of mathematical biologists. Many researchers proposed different clarifications and mathematical reconstructions of the Fisher's original statement. The present study is motivated by our opinion that the controversy can be resolved by examining the Fisher's statement within the framework of two mathematical theories that are inspired by the Darwinian formalism: evolutionary game theory (EGT) and evolutionary optimization (EO). We present four rigorous formulations (some of them previously reported) of FTNS in four different setups that come from EGT and EO. Our study demonstrates that FTNS in its original form is correct only in certain setups. In order to be recognized as a universal law, the Fisher's statement should be: (a) clarified and completed and (b) relaxed by replacing the words "is equal to" with "does not exceed". Moreover, the real meaning of FTNS can be best understood from the information-geometric point of view. Such an approach shows that FTNS imposes an upper geometric bound on information flows in evolutionary systems. In this light, FTNS appears to be a statement about the intrinsic time scale of an evolutionary system. This leads to a novel insight: FTNS is an analogue of the time-energy uncertainty relation in physics. This further emphasizes a close relation with results on speed limits in stochastic thermodynamics.

Trade-offs between number fluctuations and response in nonequilibrium chemical reaction networks

We study the response of chemical reaction networks driven far from equilibrium to logarithmic perturbations of reaction rates. The response of the mean number of a chemical species is observed to be quantitively limited by number fluctuations and the maximum thermodynamic driving force. We prove these trade-offs for linear chemical reaction networks and a class of nonlinear chemical reaction networks with a single chemical species. Numerical results for several model systems support the conclusion that these trade-offs continue to hold for a broad class of chemical reaction networks, though their precise form appears to sensitively depend on the deficiency of the network.

Fast Functionalization with High Performance in the Autonomous Information Engine

Mandal and Jarzynski have proposed a fully autonomous information heat engine, consisting of a demon, a mass, and a memory register interacting with a thermal reservoir. This device converts thermal energy into mechanical work by writing information to a memory register or, conversely, erasing information by consuming mechanical work. Here, we derive a speed limit inequality between the relaxation time of state transformation and the distance between the initial and final distributions, where the combination of the dynamical activity and entropy production plays an important role. Such inequality provides a hint that a speed-performance trade-off relation exists between the relaxation time to a functional state and the average production. To obtain fast functionalization while maintaining the performance, we show that the relaxation dynamics of the information heat engine can be accelerated significantly by devising an optimal initial state of the demon. Our design principle is inspired by the so-called Mpemba effect, where water freezes faster when initially heated.

Topological Speed Limit

Any physical system evolves at a finite speed that is constrained not only by the energetic cost but also by the topological structure of the underlying dynamics. In this Letter, by considering such structural information, we derive a unified topological speed limit for the evolution of physical states using an optimal transport approach. We prove that the minimum time required for changing states is lower bounded by the discrete Wasserstein distance, which encodes the topological information of the system, and the time-averaged velocity. The bound obtained is tight and applicable to a wide range of dynamics, from deterministic to stochastic, and classical to quantum systems. In addition, the bound provides insight into the design principles of the optimal process that attains the maximum speed. We demonstrate the application of our results to chemical reaction networks and interacting many-body quantum systems.

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.

Speed Limit for a Highly Irreversible Process and Tight Finite-Time Landauer's Bound

Landauer's bound is the minimum thermodynamic cost for erasing one bit of information. As this bound is achievable only for quasistatic processes, finite-time operation incurs additional energetic costs. We find a tight finite-time Landauer's bound by establishing a general form of the classical speed limit. This tight bound well captures the divergent behavior associated with the additional cost of a highly irreversible process, which scales differently from a nearly irreversible process. We also find an optimal dynamics which saturates the equality of the bound. We demonstrate the validity of this bound via discrete one-bit and coarse-grained bit systems. Our Letter implies that more heat dissipation than expected occurs during high-speed irreversible computation.

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.