the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# Review article: Large fluctuations in non-equilibrium physics

### Giovanni Jona-Lasinio

Non-equilibrium is dominant in geophysical and climate phenomena. However the study of non-equilibrium is much more difficult than equilibrium, and the relevance of probabilistic simplified models has been emphasized. Large deviation rates have been used recently in climate science. In this paper, after recalling progress during the last decades in understanding the role of large deviations in a class of non-equilibrium systems, we point out differences between equilibrium and non-equilibrium. For example, in non-equilibrium (a) large deviation rates may be extensive but not simply additive. (b) In non-equilibrium there are generically long-range space correlations, so large deviation rates are non-local. (c) Singularities in large deviation rates denote the existence of phase transitions often not possible in equilibrium. To exemplify, we shall refer to lattice gas models like the symmetric simple exclusion process and other models which are playing an important role in the understanding of non-equilibrium physics. The reasons why all this may be of interest in climate physics will be briefly indicated.

This paper is an enlarged version of a seminar in the series
*Perspectives in Climate Sciences 2021* and provides an
introduction to the main ideas of the so-called macroscopic
fluctuation theory (MFT). This theory developed mainly for
non-equilibrium diffusive systems is very well supported by
mathematical models and provides a methodological approach that can be
followed for other non-equilibrium systems.

The theory is based on variational principles determining for each model the large deviation rates (LDRs) of thermodynamic variables like a density or a time average like the time-averaged current flowing through the system in stationary states. Models in climate science are considerably more complex than those for which MFT has been proved to hold. However, MFT may provide a guide for more complex problems, and it shows that there are considerable differences in fluctuations from an equilibrium or a non-equilibrium state.

In the last decades, the theory of large deviations has become a main tool in statistical mechanics especially in the study of non-equilibrium. MFT has been formulated and verified in probabilistic models of lattice gases which macroscopically lead to hydrodynamic diffusion equations, which in turn represent laws of large numbers. Therefore, large deviations from hydrodynamic behavior have been studied. Here, we shall follow the formulation given in Bertini et al. (2001, 2002) and the reviews in Bertini et al. (2007, 2015). The reader is advised to refer also to Derrida (2007, 2011). Although the theory has been developed especially for diffusions, several conclusions seem to hold more generally, e.g., for reaction–diffusion systems (Jona-Lasinio et al., 1993). Reaction–diffusion systems include important energy balance models in geophysics (Sellers, 1969; Ghil, 1976).

Non-equilibrium is dominant in geophysical and climate phenomena, and we refer to the following papers for the relevance of large deviation estimates in this domain (Galfi et al., 2019; Galfi and Lucarini, 2021; Galfi et al., 2021); see also Lucarini et al. (2022). In Galfi et al. (2019), a simplified model of the general circulation of the atmosphere is adopted and assumed to be a chaotic dynamical system so that probabilistic concepts and methods can be applied. Stochastic climate models were introduced by Hasselmann (1976). For a recent assessment of the Hasselmann program, we refer to Lucarini and Chekroun (2023).

Non-equilibrium includes an enormous variety of phenomena, so we cannot hope to formulate a unique theory having a generality comparable to classical thermodynamics. We have to restrict to subclasses of problems. One difficulty is to define suitable thermodynamic functionals in far from equilibrium situations. Large fluctuations have offered a way out, as large deviation rates provide genuine thermodynamic functionals in non-equilibrium stationary states.

Studying irreversible processes is much more difficult than understanding equilibrium phenomena. In equilibrium statistical mechanics, we do not have to solve any equation of motion, and the Gibbs distribution provides the basis for the calculation of macroscopic quantities and their fluctuations. In non-equilibrium, we cannot bypass the dynamics even in the study of stationary states, which we may consider as the simplest beyond equilibrium. Examples are the heat flow in an iron rod whose endpoints are thermostated at different temperatures or the stationary flow of electrical current in a given potential difference.

For such states, the fluctuations exhibit novel and rich features with respect to the equilibrium situation. As experimentally observed (Dorfman et al., 1994), the space correlations of the thermodynamic variables generically extend to macroscopic distances, which means, for instance, that the fluctuations of the density in different points of the system are not independent. This is reflected in the structure of large deviation rates, which in non-equilibrium stationary states are in general non-local in space.

In this paper, we consider a class of systems that behave macroscopically as diffusions. The systems considered locally deviate only slightly from equilibrium so that small gradients and linear response to external fields are reasonable approximations. Microscopically, this implies that the system reaches a local equilibrium in a time which is short compared to the times typical of macroscopic evolution. So what characterizes situations in which this description applies is a separation of scales both in space and time. Far from equilibrium states are those which exhibit differences over the size of the whole system. In other words, far from equilibrium is obtained from building up small local differences. Local equilibrium is the first assumption we make. For the relevance of local equilibrium in climate models, see Lorenz (1967), Peixoto and Oort (1992).

This restriction allows us to understand some typical phenomena induced by non-equilibrium in somewhat ideal cases which however give a hint of what may happen in more realistic cases.

The content of the paper is as follows: in the next section, we recall the use by Einstein of the Boltzmann relation between entropy and probability to estimate the probability of a thermodynamic fluctuation in equilibrium. In Sect. 3, we describe the essentials of the macroscopic fluctuation theory in the version reviewed in Bertini et al. (2015). Section 4 is devoted to the illustration of the differences between equilibrium and non-equilibrium through explicit calculations for the symmetric simple exclusion process (SSEP). Final comments are in Sect. 5. Our exposition is not very detailed, and if the reader wants to understand the underpinning calculations, we suggest consulting the reviews by Bertini et al. (2015), Derrida (2007, 2011), Mallick (2015) and references therein.

The first explicit large deviation theory in equilibrium states is
presumably Einstein's theory of opalescence (Einstein, 1910); see also
Landau and Lifshitz (1980). He uses the Boltzmann–Planck formula connecting entropy
and probability to estimate probabilities of fluctuations of
thermodynamic variables, e.g., of densities. Therefore, entropy appears
as a large deviation rate. The small parameter representing the
intensity of the noise in this case is the Boltzmann constant
(Bertini et al., 2015), $\mathit{\kappa}=R/N$, *R* is the gas constant, and *N* is Avogadro's
number, so that its inverse is proportional to the typical number of
the microscopic degrees of freedom of the system. Einstein's paper,
which is entirely correct far from the critical point, unfortunately
contains an improper application to the critical point where long-range correlations are present. This paper was criticized later by
Ornstein and Zernike (1914) but provided an important source of
inspiration for non-equilibrium.

Starting from what he calls the Boltzmann principle,

it is interesting to quote from Einstein (1910) the following:

Wis commonly equated with the number of different possible ways (complexions) in which the state considered – which is incompletely defined in the sense of a molecular theory, by parameters of a system – can conceivably be realized. To be able to calculateW, one needs a complete theory of the system under consideration. If considered from a phenomenological point of view Eq. (1) appears devoid of content.However, Boltzmann's principle does acquire some content independent of any elementary theory if one assumes and generalizes from molecular kinetics the proposition that the irreversibility of physical processes is only apparent.

It follows from Eq. (1) that

$$\begin{array}{}\text{(2)}& W=\text{const}\cdot {e}^{{\scriptscriptstyle \frac{N}{R}}S}.\end{array}$$The order of magnitude of the constant is determined by taking into account that for the state of maximum entropy (entropy

S_{0}),Wis of the order of magnitude 1, so that we then have, with order of magnitude accuracy,$$\begin{array}{}\text{(3)}& W={e}^{{\scriptscriptstyle \frac{N}{R}}(S-{S}_{\mathrm{0}})}.\end{array}$$From this, we can conclude that the probability d

Wthat the quantitiesλ_{1}….λ_{k}lie betweenλ_{1}andλ_{1}+dλ_{1}......λ_{k}, andλ_{k}+dλ_{k}is given, in order of magnitude, by the equation$$\begin{array}{}\text{(4)}& \mathrm{d}W={e}^{{\scriptscriptstyle \frac{N}{R}}(S-{S}_{\mathrm{0}})}\mathrm{d}{\mathit{\lambda}}_{\mathrm{1}}\mathrm{\dots}\mathrm{d}{\mathit{\lambda}}_{k},\end{array}$$in the case when the system is determined only incompletely by

λ_{1}….λ_{k}.

In 1931, Onsager (1931a, b), in the same vein as Einstein (1910) that was quoted, made use of the Boltzmann formula in the study of fluctuations in non-equilibrium phenomena under the condition of small global deviations from equilibrium. The theory was developed further by Onsager and Machlup in Onsager and Machlup (1953) where fluctuations of time trajectories of thermodynamic variables were considered under the same hypotheses. In the next section, we discuss how it is possible to give a meaning to a formula like Eq. (3) in the more general situation of stationary states non-necessarily close to equilibrium. Typically we think of systems in contact with thermostats at different temperatures and/or reservoirs characterized by different chemical potentials and under the action of external fields. The MFT represents a step forward with respect to Onsager and the Onsager–Machlup theory.

To study the fluctuations in states far from equilibrium, let us
analyze the meaning of the difference *S*−*S*_{0} in Eq. (3). In an
isolated system, energy is conserved so that, if the volume remains
constant, $S-{S}_{\mathrm{0}}=-\frac{F-{F}_{\mathrm{0}}}{T}$ where *F* is the Helmholtz
free energy. The quantity *F*_{0}−*F*, it represents the minimal work to
bring the equilibrium state to the state corresponding to *F* at
constant temperature and volume (Landau and Lifshitz, 1980) where *F*_{0} is the
equilibrium free energy different for different systems.

The concept of minimal work is meaningful also in non-equilibrium and can be taken as a generalization of the free energy. This is essentially the starting point of the macroscopic fluctuation theory. In the following section, we shall present the basic ideas of the MFT, stating explicitly the main assumptions, following mainly Bertini et al. (2015).

The MFT, as we mentioned, was inspired by non-equilibrium microscopic probabilistic models, or the so-called lattice gases, in particular the symmetric simple exclusion process (SSEP) for which the macroscopic dynamics are diffusive and can be proved rigorously. Also the probabilities of large deviations can be obtained and the rates computed. For a general introduction to lattice gases, we refer to Spohn (1991).

## 3.1 Macroscopic fluctuations in stationary states

The macroscopic dynamics of diffusive systems are described by
hydrodynamic equations often provided by conservation laws and
constitutive equations; that is equations expressing the current in
terms of the thermodynamic variables. More precisely on the basis of
a local equilibrium assumption, at the macroscopic level, the system is
completely described by a local density *ρ*(*t*,*x*), which may have
several components and the local current *j*(*t*,*x*). The evolution
equations are of the form

where we omit the space variable $x\in \mathrm{\Lambda}\subset {R}^{d}$.

For diffusive systems, the constitutive equations take the following form:

where the *diffusion coefficient* *D*(*ρ*) and the
*mobility* *χ*(*ρ*) are *d*×*d* symmetric and positive
definite matrices, and *E* is an external field acting on the bulk.

These equations have to be supplemented by appropriate boundary
conditions on ∂Λ, the boundary of Λ. If
*λ*(*t*,*x*), $x\in \partial \mathrm{\Lambda}$ is the chemical potential of
the external reservoirs, the boundary conditions read

where *f*(*ρ*) is the equilibrium free energy
density. Non-equilibrium depends on the boundary conditions and the
external field. It may happen that the two equilibrate each other, in
which case we have a non-homogeneous equilibrium state. Much of what
we shall say applies to equations which are not in divergence form.

We assume that the microscopic evolution is given by a Markov process
*X*_{t} which represents the configuration of the system at time
*t*. The stationary non-equilibrium state (SNS) is described by a
stationary, i.e., invariant with respect to time shifts, probability
distribution *P*_{st} over the trajectories of *X*_{t}.

The macroscopic equations are supposed to derive from an underlying
microscopic dynamics through an appropriate scaling limit where the
microscopic time is divided by a factor ≈*N*^{2} and space is
divided by ≈*N* for *N* tending to ∞ where *N* is the
number of degrees of freedom of the system.

The hydrodynamic equations represent laws of large numbers with
respect to the probability measure *P*_{st} conditioned on an initial
state *X*_{0}. The initial conditions are determined by *X*_{0}. Of course
many microscopic configurations give rise to the same value of
*ρ*_{0}(*x*). In general *ρ*_{t}(*x*) is an appropriate limit as the
number of degrees of freedom diverges.

Classically we should start from molecules interacting with realistic forces and evolving with Newtonian dynamics. This is beyond the reach of present-day mathematical theory, and much simpler models have to be adopted in the reasonable hope that some essential features are adequately captured.

We further assume that the stationary measure *P*_{st} admits a
principle of large deviations describing the fluctuations from the
typical hydrodynamic behavior. This means that the probability that
the macroscopic variable *ρ*_{t} deviates from the solutions of the
hydrodynamic equations and is close to some trajectory
${\widehat{\mathit{\rho}}}_{t}$ is exponentially small and of the following form:

where $\mathcal{J}\left(\widehat{\mathit{\rho}}\right)$ is a positive functional which
vanishes if ${\widehat{\mathit{\rho}}}_{t}$ is a solution of
Eq. (5), and $V\left({\widehat{\mathit{\rho}}}_{{t}_{\mathrm{1}}}\right)$ is the cost to produce the
initial value ${\widehat{\mathit{\rho}}}_{{t}_{\mathrm{1}}}$. The parameter *ϵ* is a
scaling factor of the order of the ratio between the microscopic
length scale (typical intermolecular distance) and the macroscopic
one. The factor *ϵ*^{−d} is of the order of the number of
particles in a macroscopic volume. The role of Avogadro's number in
Eq. (3) is played here by *ϵ*^{d}. We normalize
$V\left(\widehat{\mathit{\rho}}\right)$ so that $V\left(\widehat{\mathit{\rho}}\right)=\mathrm{0}$ when $\widehat{\mathit{\rho}}=\overline{\mathit{\rho}}$ is
the stationary solution. Therefore $\mathcal{J}\left(\widehat{\mathit{\rho}}\right)$
represents the extra cost necessary in order that the system follows
the trajectory ${\widehat{\mathit{\rho}}}_{t}$. Finally
$\mathit{\rho}\left({X}_{t}\right)\approx {\widehat{\mathit{\rho}}}_{t}$ means closeness in some metric. *V*
is the generalization to the infinite dimensional case of the
Freidlin–Wentzell *quasi-potential* (Freidlin and Wentzell, 2012) and is the large
deviation rate of the stationary probability measure. The exponent in
Eq. (8) due to the factor *ϵ*^{−d} is extensive in
space. When we consider large deviations of a time-averaged quantity
like the current (Bertini et al., 2005b), we get an extra time factor.

Let us denote by *θ* the time inversion operator defined by
$\mathit{\theta}{X}_{t}={X}_{-t}={X}_{t}^{*}$. The probability measure ${P}_{\mathrm{st}}^{*}$
describing the evolution of the time-reversed process ${X}_{t}^{*}$ is
given by the composition of *P*_{st} and *θ*^{−1} that is

We assume that the time-reversed process also admits a hydrodynamic description. This hypothesis is physically very reasonable: by acting on a system from the outside, we can invert the currents flowing through the system and for example have heat passing from lower temperature to higher temperature.

At the level of large deviations, Eq. (9) implies

where ${\widehat{\mathit{\rho}}}_{{t}_{\mathrm{1}}}$, ${\widehat{\mathit{\rho}}}_{{t}_{\mathrm{2}}}$ represents the initial and final
points of the trajectory and $V\left({\widehat{\mathit{\rho}}}_{{t}_{i}}\right)$ represents the costs associated with the creation of the fluctuations ${\widehat{\mathit{\rho}}}_{{t}_{i}}$
starting from the stationary non-equilibrium state (SNS). The
functional 𝒥^{*} vanishes on the solutions of the
hydrodynamics associated to the time-reversed process.

The physical situation we are considering is the following. The
system is in the stationary state $\overline{\mathit{\rho}}$ at $t=-\mathrm{\infty}$, but at
*t*=0 we find it in the state ${\widehat{\mathit{\rho}}}_{\mathrm{0}}$. We want to determine
the most probable trajectory followed in the spontaneous creation of
this fluctuation. According to Eq. (8), this trajectory is the one
that minimizes 𝒥 among all trajectories connecting
$\overline{\mathit{\rho}}$ to ${\widehat{\mathit{\rho}}}_{\mathrm{0}}$ in the time interval $[-\mathrm{\infty},\mathrm{0}]$. From
Eq. (10) we have

The right-hand side is minimal if ${\mathcal{J}}_{[\mathrm{0},\mathrm{\infty}]}^{*}\left(\mathit{\theta}\widehat{\mathit{\rho}}\right)=\mathrm{0}$; that is if $\mathit{\theta}\widehat{\mathit{\rho}}$ is a solution of the time-reversed hydrodynamics. The existence of such a relaxation solution is due to the fact that the stationary solution $\overline{\mathit{\rho}}$ is attractive also for the time-reversed hydrodynamics. We have therefore the following generalization of Onsager–Machlup to non-equilibrium stationary states:

In a SNS the spontaneous emergence of a macroscopic fluctuation takes place most likely following a trajectory which is the time reversal of the relaxation path according to the time-reversed hydrodynamics.

The above statement follows from assuming the existence of a time-reversed dynamics and from our general hypotheses. In equilibrium, a fluctuation emerges following the time reversal of the relaxation trajectory. As illustrated in Gabrielli et al. (1996, 1999), this property may hold even if the microscopic dynamics do not satisfy detailed balance. Therefore Onsager symmetry and the above non-equilibrium generalization can be true under rather general conditions: this is possible because going from the microscopic to the macroscopic level there is a loss of information.

From Eqs. (10) or (11) we have that the free energy is related to 𝒥 by

where the infimum is taken over all trajectories connecting $\overline{\mathit{\rho}}$
to *ρ*.

## 3.2 Density fluctuations

The functional 𝒥 for diffusive systems has the Freidlin–Wentzell form (Freidlin and Wentzell, 2012) generalized to the infinite dimensional situation.

where the kernel *K*(*ρ*) is the elliptic operator defined on
functions $\mathit{\pi}:\mathrm{\Lambda}\to \mathbf{R}$ vanishing at the boundary
∂Λ by

This operator is the generalization of the Onsager matrix. Interpreting $\mathcal{J}(\mathit{\rho},{\partial}_{t}\mathit{\rho})$ as a Lagrangian, there corresponds by Legendre duality the Hamiltonian

where *H* is the momentum conjugate to *ρ*; that is
$H=\frac{\mathit{\delta}\mathcal{L}}{\mathit{\delta}\left({\partial}_{t}\mathit{\rho}\right)}$. The scalar
product 〈,〉 means integration over *x*.

The associated Hamilton equations are

These equations with appropriate boundary conditions are the
variational equations to calculate the optimal trajectory creating
*ρ*. They characterize the MFT. They are difficult to solve for
generic dependence on *ρ* of the transport coefficients
*D*(*ρ*),*χ*(*ρ*). The variational problem has been solved for
constant *D* and quadratic *χ*; see for example Imparato et al. (2009). For
the special model *zero-range* where *V* turns out to be local,
see e.g., Bertini et al. (2001, 2002).

The quasi-potential or non-equilibrium free energy satisfies the associated Hamilton–Jacobi equation:

As we shall see, the expression of *V*(*ρ*) for SSEP can be obtained
by solving Eq. (17).

## 3.3 Current fluctuations

For lattice gases, the following expression has been derived (Bertini et al., 2005b, 2006) for the joint fluctuations of density and current:

where

Here, *j* is the actual value of the current, which is connected to
*ρ* by the continuity equation ${\partial}_{t}\mathit{\rho}+\mathrm{\nabla}\cdot j=\mathrm{0}$,
while *J*(*ρ*) is the hydrodynamic current for the given value of
*ρ*. For a simple interpretation of the exponent, think of an
electric circuit. In this case, *χ*^{−1} is the resistance, and the
double integral is the energy dissipated by *j*(*t*)−*J*(*ρ*(*t*)).

By minimizing first with respect to the current *j*, it is possible to
show that

where $\overline{\mathit{\rho}}$ is the stationary solution.

In a stationary state, it is natural to consider the fluctuations of the local time-averaged current,

For fluctuations of $\stackrel{\mathrm{\u0303}}{j}$, the following large deviation principle has been derived:

where

The set 𝒜_{τ,j} is the set of all currents *j* such
that ${\mathit{\tau}}^{-\mathrm{1}}{\int}_{\mathrm{0}}^{\mathit{\tau}}\phantom{\rule{-0.125em}{0ex}}\mathrm{d}t\phantom{\rule{0.33em}{0ex}}j(t,x)=\stackrel{\mathrm{\u0303}}{j}\left(x\right)$.

This is a more general form of a large fluctuation principle proposed
by Bodineau and Derrida (2004) called an *additivity principle*. Suppose we split a one-dimensional system in two
segments of different length *L*_{1} and *L*_{2}. In this case, we must
specify the boundary condition in the intermediate point; that is a
value *ρ* of the density. The additivity principle takes the
following form:

where *P*_{L} is the probability corresponding to the length *L*, and
*ρ*_{a} and *ρ*_{b} are the boundary values of the density. This principle
is correct for various models and equivalent to Eq. (23). However
in this approach there is no time dependence, and it is assumed that in
the variational calculation is enough to consider minimizers that are
independent of time, while in Eq. (23) we admit time
dependence. This is a non-trivial difference because as it has been
clarified a phase transition may be involved. In such a case, the
method of Bodineau and Derrida (2004) would underestimate the probability of a current
fluctuation. Such a transition has been proved to exist in the model
of Kipnis–Marchioro–Presutti (Kipnis et al., 1982) in equilibrium
(Bertini et al., 2006; Bodineau and Derrida, 2005) and found numerically in Hurtado and Garrido (2011).

## 3.4 Phase transitions

In general, the appearance of singularities in the large deviation
rates denotes the presence of a non-equilibrium phase
transition. Actually the variational principle of Bodineau and Derrida
may provide several time-independent solutions which in fact have been
found in models discussed in Baek et al. (2018) representing different
phases. There is another type of phase transition whose appearance is
signaled by the non-differentiability of the quasi-potential
*V*(*ρ*). This type of transition has been found in the weakly
asymmetric exclusion process (WASEP), that is, in presence of an
external field, for sufficiently high values of the field
(Bertini et al., 2010). The existence of non-equilibrium phase transitions,
often impossible in equilibrium, is a generic feature which has to be
taken in account when analyzing a phenomenon.

## 3.5 Long-range space correlations

We are concerned with *macroscopic correlations*, which are a
generic feature of non-equilibrium non-linear models. Microscopic
space correlations which decay as a summable power law disappear at
the macroscopic level.

We introduce the *pressure* functional as the Legendre transform
of the quasi-potential *V*:

By Legendre duality, we have the change of variable formulae
$h=\frac{\mathit{\delta}V}{\mathit{\delta}\mathit{\rho}}$, $\mathit{\rho}=\frac{\mathit{\delta}G}{\mathit{\delta}h}$, so
that the Hamilton–Jacobi equation (Eq. 17) can then be
rewritten in terms of *G* as

where *h* vanishes at the boundary of Λ. As for equilibrium
systems, *G* is the generating functional of the correlation
functions; see e.g., Amit (1978); Di Castro and Raimondi (2015).

We define

By expanding Eq. (25) around the stationary state one obtains, after non-trivial manipulations and combinatorics, recursive equations for the ${C}_{n}({x}_{\mathrm{1}},\mathrm{\dots},{x}_{n})$; see Bertini et al. (2009). We discuss the pair correlation function by splitting it into the local equilibrium part and a possibly non-local term as follows:

where

and $\overline{\mathit{\rho}}$ is the stationary solution. We then obtain from the
general equations the following equation for *B*:

where ℒ^{†} is the formal adjoint of the elliptic
operator $\mathcal{L}={L}_{x}+{L}_{y}$ given, using the usual convention
that repeated indices are summed, by the following:

and

where
$\overline{J}=J\left(\overline{\mathit{\rho}}\right)=-D\left(\overline{\mathit{\rho}}\right(x\left)\right)\mathrm{\nabla}\overline{\mathit{\rho}}\left(x\right)+\mathit{\chi}\left(\overline{\mathit{\rho}}\right(x\left)\right)E\left(x\right)$ is the macroscopic current in the stationary
state. In order that Eq. (29) may have a non-trivial solution,
*α* must be different from zero, which is generically the case
when *χ* and *D* have a non-linear dependence on *ρ*. For
non-trivial *α*, long-range correlations appear. In particular, for
the SSEP where $\mathit{\chi}\left(\mathit{\rho}\right)=\mathit{\rho}(\mathrm{1}-\mathit{\rho})$,

Here, Δ^{−1} is the Green function of the Dirichlet Laplacian.

In Basile and Jona-Lasinio (2004), it was shown that long-range correlations may appear also in equilibrium in reaction–diffusion dynamics when microscopic time-reversal invariance is strongly violated.

The simple exclusion process SSEP is the most studied system in far from equilibrium situations and has a role similar to the Ising model in the study of phase transitions. The SSEP in one-dimensional lattice is a process in which particles perform symmetric random walks subject to hard core exclusion. In non-equilibrium, the boundary conditions at the two ends of the lattice are different; an external field may act on the system, or both, so that a current is flowing through the system.

Large deviations functions have also been calculated via the MFT for other models, like the zero-range process (Bertini et al., 2002) and the Kipnis–Marchioro–Presutti model (Bertini et al., 2005a).

## 4.1 Non-locality of the quasi-potential *V* in non-equilibrium stationary states

We consider the variational problem defining *V*(*ρ*) for the
one-dimensional simple exclusion process characterized by *D*=1 and
$\mathit{\chi}=\mathit{\rho}(\mathrm{1}-\mathit{\rho})$. By performing the change of variable and
inserting it in the Hamilton–Jacobi equation (Eq. 17),

for some functional *ϕ*(*x*;*ρ*) to be determined satisfying the
boundary conditions $\mathit{\varphi}(\pm \mathrm{1})=\mathrm{log}\mathit{\rho}(\pm \mathrm{1})/[\mathrm{1}-\mathit{\rho}(\pm \mathrm{1}\left)\right]$.
We obtain a solution of the variational problem if we solve the
following ordinary differential equation which relates the functional
$\mathit{\varphi}\left(x\right)=\mathit{\varphi}(x;\mathit{\rho})$ to *ρ*:

This equation admits a unique monotone solution, which is the relevant one for the quasi-potential. A computation shows that

where *F*(*ρ*) is the equilibrium free energy

This expression was first obtained by Derrida, Lebowitz and Speer solving the microscopic model (Derrida et al., 2002). They also proved that if one splits the system in two parts, the rates obey an additivity rule more complicated than a simple sum and similar to the additivity principle of Bodineau and Derrida (2004). The above macroscopic calculation via the Hamilton–Jacobi equation was done in Bertini et al. (2002).

If an external field is present, the large deviation rate has been calculated by Enaud and Derrida (2004).

## 4.2 Non-stationary states

The time evolution depends on the initial condition. Basic work on non-stationary problems for SSEP is due to
Derrida and Gerschenfeld (2009a, b) with interesting developments in recent work by
Mallick–Moriya–Sasamoto (Mallick et al., 2022). They studied the evolution of a
step initial condition like in Fig. 1. By considering the time integral
of the local current through the origin ${Q}_{t}={\int}_{\mathrm{0}}^{t}\phantom{\rule{0.125em}{0ex}}j(\mathrm{0},{t}^{\prime})\phantom{\rule{0.125em}{0ex}}\mathrm{d}{t}^{\prime}$
where *t* is the microscopic time, we expect in a diffusive
one-dimensional regime a law of large numbers for the quantity
${Q}_{t}/\sqrt{t}$ for large *t* and a large deviation principle,

Define the cumulant generating function
$\langle {e}^{\mathit{\lambda}Q}\rangle \simeq {e}^{\surd t\mathit{\mu}\left(\mathit{\lambda}\right)}$ where *λ* is a
real parameter and 〈,〉 now stands for stochastic average. The
functions Φ(*q*) and *μ*(*λ*) are Legendre transform of each
other. This problem was studied by Derrida and Gerschenfeld (2009a, b).

The variational equations are

where *χ*(*ρ*) is the mobility. Notice that the equation for *H*
is backward in time, the boundary conditions being
$H(t,x)=\mathit{\lambda}\mathit{\theta}\left(x\right)$ and
$H(\mathrm{0},x)=\mathit{\lambda}\mathit{\theta}\left(x\right)+{f}^{\prime}\left(\mathit{\rho}\right(\mathrm{0},x\left)\right)-{f}^{\prime}\left(\overline{\mathit{\rho}}\right))$,
*f*(*ρ*) is the equilibrium free energy density.

In Derrida and Gerschenfeld (2009b), the authors calculated the moment generating function with quenched and fluctuating initial conditions, microscopically with the Bethe ansatz and macroscopically with the MFT. The variational equations could be solved only in special cases. Recently Mallick–Moriya–Sasamoto (Mallick et al., 2022) have discovered that the following non-local transformations, which generalize the Cole–Hopf transformation, allow to map the variational equations for the SSEP to the Ablowitz–Kaup–Newell–Segur (AKNS) equations:

The AKNS equations,

were solved for the SSEP, that is $\mathit{\chi}=\mathit{\rho}(\mathrm{1}-\mathit{\rho})$, with the
inverse scattering method. These transformations are valid for
quadratic *χ*. The boundary conditions for the step initial
condition are

From their solution, they obtained the moment generating function

where

For integrability results in the Kipnis–Marchioro–Presutti model of heat transfer, see the recent work by Bettelheim–Smith–Meerson (Bettelheim et al., 2022) and references therein.

The MFT shows that once macroscopic evolution equations like hydrodynamics are available and a separation of scales holds, a self-consistent macroscopic description of non-equilibrium phenomena can be obtained through a study of rare fluctuations. The origin of the probabilistic behavior may be due to the influence of a smaller scale on a larger one or to chaotic properties of the underlying dynamical system.

The discovery that a purely macroscopic theory could reproduce in the case of SSEP the large deviation function for non-equilibrium stationary states obtained from a microscopic calculation in Derrida et al. (2002) has been an important support and a source of inspiration for the MFT. The work of Derrida and Gerschenfeld (2009a, b) extended the theory to time-dependent evolutions, and its correctness received an important mathematical support by the recent work of Mallick et al. (2022).

We have illustrated the theory in an idealized case: we have
considered (i) purely diffusive systems and (ii) simplified stochastic
models. However as the Ising model allowed us to understand a lot
about phase transitions and the critical point, we believe that the
SSEP and other solvable models are providing a guide to what may
happen out of equilibrium. Furthermore the MFT applies to all variants
of SSEP or of other diffusive models that macroscopically have the
same transport coefficients *D*,*χ*.

The general approach of MFT has been extended to systems with more than one conservation law (Bernardin, 2008) and to some reaction diffusion processes like the Glauber–Kawasaki dynamics (Jona-Lasinio et al., 1993; Bodineau and Lagouge, 2010). It is reasonable to expect that it will be applicable to more general cases, which may lead to a different structure of the large deviation rates.

The experience so far indicates that the phenomenon of long-range space correlations is not limited to purely diffusive systems and is rather general in non-equilibrium.

In climate science, the models are comparatively more complicated. However the application of large deviation theory to rather complex models as in Galfi et al. (2019), Galfi and Lucarini (2021), Galfi et al. (2021) and Lucarini et al. (2022) is encouraging.

Formally, the equations of the MFT can be derived also from assuming an extension to non-equilibrium of the so-called fluctuating hydrodynamics (FHs). The idea of fluctuating hydrodynamics goes back to Landau (Landau and Lifshitz, 1980) who considered only linear hydrodynamic equations near equilibrium, while the most interesting phenomena are generated by non-linearities far from equilibrium. It consists in adding to the macroscopic equations a noise term. In the case of hydrodynamic equations in divergence form, we add to the current a fluctuating term,

where, conditionally on *ρ*, *ξ* is a Gaussian random term with
variance as follows:

The hydrodynamic equation takes the following form:

However we have to emphasize that such equations, in the case of non-linear hydrodynamics, are very singular and need renormalization, and so far there is no mathematical theory applicable even in one dimension.

From the previous equations, we obtain

the same as formula (18).

It would be interesting to provide a rigorous foundation to fluctuating hydrodynamics; however this requires as a first step to give a clear mathematical meaning to the stochastic partial differential equations on which it is founded.

No data sets were used in this article.

The author has declared that there are no competing interests.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This article is part of the special issue “Interdisciplinary perspectives on climate sciences – highlighting past and current scientific achievements”. It is not associated with a conference.

I wish to acknowledge my long-standing collaboration on the topics of this paper with Lorenzo Bertini, Alberto De Sole, Davide Gabrielli and Claudio Landim. I thank Vera Melinda Galfi and Valerio Lucarini for a critical reading of the paper and for very useful comments.

This paper was edited by Vera Melinda Galfi and reviewed by Valerio Lucarini and one anonymous referee.

Amit, D. J.: Field Theory the Renormalization Group, and Critical Phenomena, McGraw-Hill, 1978, ISBN 9780070015753, 1978. a

Baek, Y., Kafri, Y., and Lecomte, V.: Dynamical phase transitions in the current distribution of driven diffusive channels, J. Phys. A, 51, 105001, https://doi.org/10.1088/1751-8121/aaa8f9, 2018. a

Basile, G. and Jona-Lasinio, G.: Equilibrium states with macroscopic correlations, Int. J. Mod. Phys. B, 18, 479, https://doi.org/10.1142/s0217979204024094, 2004. a

Bernardin, C.: Stationary non-equilibrium propeties for a heat conduction model, Phys. Rev. E, 78, 021134, https://doi.org/10.1103/physreve.78.021134, 2008. a

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Fluctuations in stationary nonequilibrium states of irreversible processes, Phys. Rev. Lett., 87, 040601, https://doi.org/10.1103/physrevlett.87.040601, 2001. a, b

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Macroscopic fluctuation theory for stationary non-equilibrium states, J. Stat. Phys., 107, 635–675, 2002. a, b, c, d

Bertini, L., Gabrielli, D., and Lebowitz, J. L.: Large Deviations for a Stochastic Model of Heat Flow, J. Stat. Phys., 121, 843–885, https://doi.org/10.1007/s10955-005-5527-2, 2005a. a

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Current fluctuations in stochastic lattice gases, Phys. Rev. Lett., 94, 030601, https://doi.org/10.1103/physrevlett.94.030601, 2005b. a, b

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Non-equilibrium current fluctuations in stochastic lattice gases, J. Stat. Phys., 123, 237–276, 2006. a, b

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Stochastic interacting particle systems out of equilibrium, J. Stat. Mech. Theory Exp., 2007, P07014, https://doi.org/10.1088/1742-5468/2007/07/p07014, 2007. a

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Towards a non-equilibrium thermodynamics: a selfcontained macroscopic description of driven diffusive systems, J. Stat. Phys., 135, 857–872, 2009. a

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Lagrangian phase transitions in non-equilibrium thermodynamic systems, J. Stat. Mech. Theory Exp., 2010, L11001, https://doi.org/10.1088/1742-5468/2010/11/l11001, 2010. a

Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Macroscopic Fluctuation Theory, Rev. Mod. Phys., 87, 593–636, 2015. a, b, c, d, e

Bettelheim, E., Smith, N. R., and Meerson, B.: Full statistics of nonstationary heat transfer in the Kipnis–Marchioro–Presutti model, ArXiv, arXiv:2204.06278, 2022. a

Bodineau, T. and Derrida, B.: Current fluctuations in non-equilibrium diffusive systems: an additivity principle, Phys. Rev. Lett., 92, 180601, https://doi.org/10.1103/physrevlett.92.180601, 2004. a, b, c

Bodineau, T. and Derrida, B.: Distribution of current in non-equilibrium diffusive systems and phase transitions, Phys. Rev. E, 72, 066110, https://doi.org/10.1103/physreve.72.066110, 2005. a

Bodineau, T. and Lagouge, M.: Current large deviations in a driven dissipative model, J. Stat. Phys., 139, 201, https://doi.org/10.1007/s10955-010-9934-7, 2010. a

Derrida, B.: Non-equilibrium steady states: fluctuations and large deviations of the density and of the current, J. Stat. Mech. Theory Exp., 2007, P07023, https://doi.org/10.1088/1742-5468/2007/07/p07023, 2007. a, b

Derrida, B.: Microscopic versus macroscopic approaches to non-equilibrium systems, J. Stat. Mech. Theory Exp., 2011, P01030, https://doi.org/10.1088/1742-5468/2011/01/p01030, 2011. a, b

Derrida, B. and Gerschenfeld, A.: Current fluctuations in the one-dimensional symmetric simple exclusion process with step initial condition, J. Stat. Phys., 136, 1, https://doi.org/10.1007/s10955-009-9772-7, 2009a. a, b, c

Derrida, B. and Gerschenfeld, A.: Current fluctuations in the one-dimensional diffusive systems with step initial density profile, J. Stat. Phys., 137, 978, https://doi.org/10.1007/s10955-009-9830-1, 2009b. a, b, c, d

Derrida, B., Lebowitz, J. L., and Speer, E. R.: Large deviation of the density profile in the steady state of the open symmetric exclusion process, J. Stat. Phys., 107, 599–634, 2002. a, b

Di Castro, C. and Raimondi, R.: Statistical Mechanics and Applications in Condensed Matter, Cambridge University Press, ISBN 9781139600286, https://doi.org/10.1017/CBO9781139600286, 2015. a

Dorfman, J. R., Kirkpatrick, T. R., and Sengers, J. V.: Generic Long Range Correlations in Molecular Fluids, Ann. Rev. Phys. Chem., 45, 213–239, 1994. a

Einstein, A.: Theorie der Opaleszenz von homogenen Flüssigkeiten und Flüssigkeitsgemischen in der Nähe des kritischen Zustandes, Annalen der Physik, 33, 1275–1298, 1910 (English translation in The collected papers of Albert Einstein, 3, 231–249, Princeton University Press, 1993). a, b, c

Enaud, C. and Derrida, B.: Large Deviation Functional of the Weakly Asymmetric Exclusion Process, J. Stat. Phys., 114, 537–562, 2004. a

Freidlin, M. and Wentzell, A.: Random Perturbations of Dynamical Systems, 3rd edn., Springer, https://doi.org/10.1007/978-3-642-25847-3_8, 2012. a, b

Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Onsager Reciprocity Relations without Microscopic Reversibility, Phys. Rev. Lett., 77, 1202, https://doi.org/10.1103/physrevlett.77.1202, 1996. a

Gabrielli, D., Jona-Lasinio, G., and Landim, C.: Onsager Symmetry from Microscopic TP Invariance, J. Stat. Phys., 96, 639, https://doi.org/10.1023/a:1004550307453, 1999. a

Galfi, V. M. and Lucarini, V.: Fingerprinting heatwaves and cold spells and assessing their response to climate change using large deviation theory, Phys. Rev. Lett., 127, 058701, https://doi.org/10.1103/physrevlett.127.058701, 2021. a, b

Galfi, V. M., Lucarini, V., and Wouters, J.: A Large Deviation Theory-based Analysis of Heat Waves and Cold Spells in a Simplified Model of the General Circulation of the Atmosphere, J. Stat. Mech., 2019, 033404, https://doi.org/10.1088/1742-5468/ab02e8, arXiv:1807.08261, 2019. a, b, c

Galfi, V. M., Lucarini, V., Ragone, F., and Wouters, J.: Applications of large deviation theory in geophysical fluid dynamics and climate science, Riv. Nuovo Cim., 44, 291363, https://doi.org/10.1007/s40766-021-00020z, 2021. a, b

Ghil, M.: Climate stability for a Sellers type model, J. Atmos. Sci., 33, 3–20, 1976. a

Hasselmann, K.: Stochastic climate models Part I. Theory, Tellus, 28, 473–485, 1976. a

Hurtado, P. I. and Garrido, P. L.: Spontaneous symmetry breaking at the fluctuating level, Phys. Rev. Lett., 107, 180601, https://doi.org/10.1103/physrevlett.107.180601, 2011. a

Imparato, A., Lecomte, V., van and Wijland, F.: Equilibrium-like fluctuations in some boundary-driven open diffusive systems, Phys. Rev. E, 80, 011131, https://doi.org/10.1103/physreve.80.011131, 2009. a

Jona-Lasinio, G., Landim, C., and Vares, M. E.: Large deviations for a reaction-diffusion model, Probab. Theory Related Fields, 97, 339–361, https://doi.org/10.1007/bf01195070, 1993. a, b

Kipnis, C., Marchioro, C., and Presutti, E.: Heat flow in an axactly solvable model, J. Stat. Phys., 27, 65–74, https://doi.org/10.1007/bf01011740, 1982. a

Landau, L. D. and Lifshitz, E. M.: Statistical Physics, 3rd edn., chap. XII, ISBN 9780750633727, 1980. a, b, c

Lorenz, E. N.: The Nature and Theory of the General Circulation of the Atmosphere, WMO Publication, 218, World Meteorological Organization, Geneva, 218, 115, 1967. a

Lucarini, V. and Chekroun, M.: Hasselmann's Program and Beyond: New Theoretical Tools for Understanding the Climate Crisis, ArXiv, arXiv:2303.12009, 2023. a

Lucarini, V., Serdukova, L., and Margazoglou, G.: Lévy noise versus Gaussian-noise-induced transitions in the Ghil–Sellers energy balance model, Nonlin. Processes Geophys., 29, 183–205, https://doi.org/10.5194/npg-29-183-2022, 2022. a, b

Mallick, K.: The Exclusion Process: A paradigm for non-equilibrium behaviour, Physica A, 418, 17–48, 2015. a

Mallick, K., Moriya, H., and Sasamoto, T.: Exact solution of the macroscopic fluctuation theory for the symmetric exclusion process, ArXiv, arXiv:2202.05213, 2022. a, b, c

Onsager, L.: Reciprocal Relations in Irreversible Processes. I, Phys. Rev., 37, 405–426, https://doi.org/10.1103/physrev.37.405, 1931a. a

Onsager, L.: Reciprocal Relations in Irreversible Processes. II, Phys. Rev., 38, 2265–2279, https://doi.org/10.1103/physrev.38.2265, 1931b. a

Onsager, L. and Machlup, S.: Fluctuations and irreversible processes, Phys. Rev., 91, 1505–1512, 1953. a

Ornstein, L. S. and Zernike, F.: Accidental deviations of density and opalescence at the critical point of a single substance, Proc. Acad. Sci. (Amsterdam), 17, 793–806, 1914. a

Peixoto, J. P. and Oort, A. H.: Physics of Climate, AIP Press, New York, ISBN 9780883187128, 1992. a

Sellers, W. D.: A global climatic model based on the energy balance of the earth atmosphere, J. Appl. Meteorol., 8, 392–400, 1969. a

Spohn, H.: Large Scale Dyhamics of Interacting Particles, Springer, Berlin, ISBN 9783642843716, https://doi.org/10.1007/978-3-642-84371-6, 1991. a

- Abstract
- Introduction
- Macroscopic fluctuations in equilibrium
- Macroscopic fluctuation theory (MFT)
- The simple exclusion process (SSEP)
- Final remarks
- Appendix A: Fluctuating hydrodynamics
- Data availability
- Competing interests
- Disclaimer
- Special issue statement
- Acknowledgements
- Review statement
- References

- Abstract
- Introduction
- Macroscopic fluctuations in equilibrium
- Macroscopic fluctuation theory (MFT)
- The simple exclusion process (SSEP)
- Final remarks
- Appendix A: Fluctuating hydrodynamics
- Data availability
- Competing interests
- Disclaimer
- Special issue statement
- Acknowledgements
- Review statement
- References