Review article: Large ﬂuctuations in non-equilibrium physics

. Non-equilibrium is dominant in geophysical and climate phenomena. However the study of non-equilibrium is much more difﬁcult than equilibrium, and the relevance of probabilistic simpliﬁed 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 brieﬂy indicated.


Introduction
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 ver-ified 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. (2001Bertini et al. ( , 2002) ) and the reviews in Bertini et al. (2007Bertini et al. ( , 2015)).The reader is advised to refer also to Derrida (2007Derrida ( , 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 Published by Copernicus Publications on behalf of the European Geosciences Union & the American Geophysical Union.
G. Jona-Lasinio: Large fluctuations in non-equilibrium physics 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 (2007Derrida ( , 2011)), Mallick (2015) and references therein.

Macroscopic fluctuations in equilibrium
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), κ = 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 longrange 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: W is 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 calculate W , 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 The order of magnitude of the constant is determined by taking into account that for the state of maximum entropy (entropy S 0 ), W is of the order of magnitude 1, so that we then have, with order of magnitude accuracy, From this, we can conclude that the probability dW that 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 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 0 = − F −F 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 nonequilibrium 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).

Macroscopic fluctuation theory (MFT)
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).

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 ∈ ⊂ 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 ∈ ∂ is the chemical potential of the external reservoirs, the boundary conditions read where f (ρ) is the equilibrium free energy density.Nonequilibrium 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 ρt is exponentially small and of the following form: where J ( ρ) is a positive functional which vanishes if ρt is a solution of Eq. ( 5), and V ( ρt 1 ) is the cost to produce the initial value ρt 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 ( ρ) so that V ( ρ) = 0 when ρ = ρ is the stationary solution.Therefore J ( ρ) represents the extra cost necessary in order that the system follows the trajectory ρt .Finally ρ(X t ) ≈ ρ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 θ X t = X −t = X * t .The probability measure P * 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 ρt 1 , ρt 2 represents the initial and final points of the trajectory and V ( ρt i ) represents the costs associated with the creation of the fluctuations ρt i starting from the stationary non-equilibrium state (SNS).The functional J * vanishes on the solutions of the hydrodynamics associated to the timereversed process.
The physical situation we are considering is the following.The system is in the stationary state ρ at t = −∞, but at t = 0 we find it in the state ρ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 J among all trajectories connecting ρ to ρ0 in the time interval [−∞, 0].From Eq. ( 10) we have The right-hand side is minimal if J * [0,∞] (θ ρ) = 0; that is if θ ρ is a solution of the time-reversed hydrodynamics.The existence of such a relaxation solution is due to the fact that the stationary solution ρ 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. (1996Gabrielli et al. ( , 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 J by where the infimum is taken over all trajectories connecting ρ to ρ.

Density fluctuations
The functional J 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 π : → R vanishing at the boundary ∂ by This operator is the generalization of the Onsager matrix.Interpreting J (ρ, ∂ t ρ) as a Lagrangian, there corresponds by Legendre duality the Hamiltonian where H is the momentum conjugate to ρ; that is H = δL δ(∂ t ρ) .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. (2001Bertini et al. ( , 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).

Current fluctuations
For lattice gases, the following expression has been derived (Bertini et al., 2005b(Bertini et al., , 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 ∂ t ρ + ∇ • j = 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 ρ is the stationary solution.
In a stationary state, it is natural to consider the fluctuations of the local time-averaged current, For fluctuations of j , the following large deviation principle has been derived: where The set A τ,j is the set of all currents j such that τ −1 τ 0 dt j (t, x) = j (x).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 nontrivial 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).

Phase transitions
In general, the appearance of singularities in the large deviation rates denotes the presence of a non-equilibrium phase https://doi.org/10.5194/npg-30-253-2023 Nonlin.Processes Geophys., 30, 253-262, 2023 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.

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 = δV δρ , ρ = δG δ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 1 , . .., 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 ρ is the stationary solution.We then obtain from the general equations the following equation for B: where L † is the formal adjoint of the elliptic operator L = L x + L y given, using the usual convention that repeated indices are summed, by the following: and where J = J ( ρ) = −D( ρ(x))∇ ρ(x) + χ( ρ(x))E(x) 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 α, longrange correlations appear.In particular, for the SSEP where 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 timereversal invariance is strongly violated.

The simple exclusion process (SSEP)
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).

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 χ = ρ(1 − ρ).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 φ(±1) = log ρ(±1)/[1 − ρ(±1)].We obtain a solution of the variational problem if we solve the following ordinary differential equation which relates the functional φ(x) = φ(x; ρ) 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).

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 = t 0 j (0, t ) dt where t is the microscopic time, we expect in a diffusive one-dimensional regime a law of large numbers for the quantity Q t / √ t for large t and a large deviation principle, 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) = λθ (x) and H (0, x) = λθ (x)+f (ρ(0, x))−f ( ρ)), 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 χ = ρ(1 − ρ), with the inverse scattering method.These transformations are valid for quadratic χ.The boundary conditions for the step initial condition are u(0, x) = ωδ(x) v(t, x) = δ(x). (43) From their solution, they obtained the moment generating function 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.

Figure 1 .
Figure 1.Step initial condition with a density ρ a at the left of the origin and ρ b at the right of the origin.