Detailed balance, local detailed balance, and global potential for stochastic chemical reaction networks

Detailed balance of a chemical reaction network can be defined in several different ways. Here we investigate the relationship among four types of detailed balance conditions: deterministic, stochastic, local, and zero-order local detailed balance. We show that the four types of detailed balance are equivalent when different reactions lead to different species changes and are not equivalent when some different reactions lead to the same species change. Under the condition of local detailed balance, we further show that the system has a global potential defined over the whole space, which plays a central role in the large deviation theory and the Freidlin–Wentzell-type metastability theory of chemical reaction networks. Finally, we provide a new sufficient condition for stochastic detailed balance, which is applied to construct a class of high-dimensional chemical reaction networks that both satisfies stochastic detailed balance and displays multistability.

Applications of large deviation theory in geophysical fluid dynamics and climate science

The climate is a complex, chaotic system with many degrees of freedom. Attaining a deeper level of understanding of climate dynamics is an urgent scientific challenge, given the evolving climate crisis. In statistical physics, many-particle systems are studied using Large Deviation Theory (LDT). A great potential exists for applying LDT to problems in geophysical fluid dynamics and climate science. In particular, LDT allows for understanding the properties of persistent deviations of climatic fields from long-term averages and for associating them to low-frequency, large-scale patterns. Additionally, LDT can be used in conjunction with rare event algorithms to explore rarely visited regions of the phase space. These applications are of key importance to improve our understanding of high-impact weather and climate events. Furthermore, LDT provides tools for evaluating the probability of noise-induced transitions between metastable climate states. This is, in turn, essential for understanding the global stability properties of the system. The goal of this review is manifold. First, we provide an introduction to LDT. We then present the existing literature. Finally, we propose possible lines of future investigations. We hope that this paper will prepare the ground for studies applying LDT to solve problems encountered in climate science and geophysical fluid dynamics.

Rare Event Analysis for Minimum Hellinger Distance Estimators via Large Deviation Theory

Hellinger distance has been widely used to derive objective functions that are alternatives to maximum likelihood methods. While the asymptotic distributions of these estimators have been well investigated, the probabilities of rare events induced by them are largely unknown. In this article, we analyze these rare event probabilities using large deviation theory under a potential model misspecification, in both one and higher dimensions. We show that these probabilities decay exponentially, characterizing their decay via a “rate function” which is expressed as a convex conjugate of a limiting cumulant generating function. In the analysis of the lower bound, in particular, certain geometric considerations arise that facilitate an explicit representation, also in the case when the limiting generating function is nondifferentiable. Our analysis involves the modulus of continuity properties of the affinity, which may be of independent interest.

Generalized Euler-Lotka equation for correlated cell divisions

Cell division times in microbial populations display significant fluctuations. These fluctuations impact the population growth rate in a non-trivial way. If fluctuations are uncorrelated among different cells, the population growth rate is predicted by the Euler-Lotka equation, which is a classic result in mathematical biology. However, cell division times can present significant correlations, due to physical properties of cells that are passed from mothers to daughters. In this paper, we derive an equation remarkably similar to the Euler-Lotka equation which is valid in the presence of correlations. Our exact result is based on large deviation theory and does not require particularly strong assumptions on the underlying dynamics. We apply our theory to a phenomenological model of bacterial cell division. We find that the discrepancy between the growth rate predicted by the Euler-Lotka equation and our generalized version is relatively small, but large enough to be measurable in experiments.