Abstract
We propose a general formalism, within large-deviation theory, giving access to the exact statistics of fluctuations of ballistically transported conserved quantities in homogeneous, stationary states. The formalism is expected to apply to any system with an Euler hydrodynamic description, classical or quantum, integrable or not, in or out of equilibrium. We express the exact scaled cumulant generating function (or full counting statistics) for any (quasi-)local conserved quantity in terms of the flux Jacobian. We show that the “extended fluctuation relations” of Bernard and Doyon follow from the linearity of the hydrodynamic equations, forming a marker of “freeness” much like the absence of hydrodynamic diffusion does. We show how an extension of the formalism gives exact exponential behaviours of spatio-temporal two-point functions of twist fields, with applications to order-parameter dynamical correlations in arbitrary homogeneous, stationary state. We explain in what situations the large-deviation principle at the basis of the results fail, and discuss how this connects with nonlinear fluctuating hydrodynamics. Applying the formalism to conformal hydrodynamics, we evaluate the exact cumulants of energy transport in quantum critical systems of arbitrary dimension at low but nonzero temperatures, observing a phase transition for Lorentz boosts at the sound velocity.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Far-from-equilibrium physics has seen a large amount of theoretical and experimental developments in recent years [1,2,3,4,5]. A distinctive feature of such non-equilibrium states is entropy production and the breaking of time-reversal invariance. This is associated with the existence nonzero currents describing transport of various quantities, such as particles, charge or energy. These effects are ubiquitous in nature and of fundamental importance. Despite this, there exists no fully satisfying non-equilibrium parallel to equilibrium thermodynamics, that is, organising principles for the behaviour and statistics of non-equilibrium currents. A promising avenue is the study of fluctuations, and a widely popular approach is the use of large-deviation theory (LDT) [6, 7], which studies the rare but significant fluctuations around almost-sure values of macroscopic quantities. In particular, fluctuations in transport encode many universal properties of non-equilibrium physics (see [8,9,10,11,12]). LDT offers a conceptual link with thermodynamics and gives a general understanding of a wide class of non-equilibrium phenomena.
In the context of non-equilibrium transport, it is natural to focus on the LDT for the total transfer of quantities, at long times, between two or more macroscopic regions. As an illustration, consider the total energy that has passed from the left to the right halves of an infinitely large system after a long time t; one can enquire about the distribution of this random variable as \(t\rightarrow \infty \). Of particular interest are the cumulants, scaled by 1 / t, of the transferred quantity in the large-t limit. These are finite if the expected large-deviation principle holds, and are encoded within the scaled cumulant generating function (SCGF), or full counting statistics, a non-equilibrium counterpart to the equilibrium free energy. The LDT for transport has been studied in many systems, with the SCGF often calculated exactly. For quantum transport of free fermions, the SCGF of U(1) charges is given by the celebrated Levitov–Lesovik formula [13, 14], which has applications in mesoscopic physics. Free-particle advanced techniques have been used [15,16,17,18,19,20] (see also [21]), and exact results exist in certain integrable impurity models [22] and in general \(1+1\)-dimensional conformal field theory (CFT) [23, 24] (see the review [25]). In these studies, non-equilibrium currents are generated by the partitioning protocol [26,27,28] (see also [25] and references therein). In this protocol, the left and right halves of the system are initially in different ensembles, e.g., different temperatures. The system is then left to evolve for a long time, and the final state of the system is observed in finite regions around the origin. For systems that admit ballistic transport, one generically obtains a non-equilibrium steady state where steady transport occurs. There are also many exact results for the SCGF in classical stochastic many-body systems such as “exclusion processes” (see, e.g., [10,11,12]). Many techniques have been used, and a successful framework is macroscopic fluctuation theory (MFT) [29,30,31,32,33], based on a hydrodynamic description and taking as input diffusion coefficients. Exact results for fluctuations in open or stochastic quantum systems have also been obtained (see, e.g., [34,35,36,37,38,39], and, recently, in certain cellular automata [40]).
In this paper, we propose a general theory for studying the LDT of ballistic transport in homogeneous, stationary, maximal-entropy states (MES) of deterministic many-body systems. MES are temporally and spatially homogeneous states (thermodynamic ensembles) reached after relaxation processes have occurred [5]. In such states, entropy is maximised with respect to all available local conservation laws. They include thermal Gibbs states, Galilean or relativistic boosts thereof, and, in integrable systems, generalised Gibbs ensembles [41,42,43,44] (experimentally observed [45]). MES may admit ballistic currents if there are conserved charges which are odd under time reversal, such as the momentum; in such cases, these are non-equilibrium steady states. In particular, they include the non-equilibrium steady states emerging in the partitioning protocol [25, 46, 47].
In the set-up we consider, time evolution is deterministic, and initial states are fluctuating. In a hydrodynamic description [48], the presence of ballistic transport leads to nontrivial Euler-scale hydrodynamic equations, and to Euler-scale linear fluctuating hydrodynamics. We use aspects of these (see [49, App A], [50, 51]), in order to propose a framework for the LDT of ballistic transport, expressing exact SCGFs in terms of the flux Jacobian (or linearised Euler matrix). Our proposal somewhat parallels MFT, being based on a hydrodynamic description, but differs from MFT in that it only necessitates the Euler scale—diffusion and other higher-order contributions would give subleading corrections. It can be applied to a large variety of systems with Euler hydrodynamics, quantum or classical, integrable or not, interacting or not.
The theory is based on biasing the measure by a total time-integrated current, thus accessing rare fluctuations and explicitly generating the scaled cumulants. Such a bias is a widely used technique in stochastic systems (sometimes referred to as exponential tilting, or s-ensemble), and how it gives rise to a new stochastic dynamics is referred to as the (classical or quantum) generalised Doob transform (see [52, 53] and references therein). Here, instead, we find the exact modification of the initial state that reproduces this bias.
The proposal provides an organising principle for all results for exact SCGF in homogeneous free-particle models and \(1+1\)-dimensional CFT; it is a nonlinear generalisation of principles found for these systems [21, 23,24,25]. We clarify the origin of the “extended fluctuation relations” of Bernard and Doyon [21], showing that they arise when the flux Jacobian, in the coordinates of the conserved densities, is state independent (the hydrodynamic equations are linear). We propose that this is a property of a many-body system that characterises it as being “free” by opposition to “interacting”, much like the absence of diffusion is [54].
The theory generalises to the total integrated currents along arbitrary rays (space-time points \(\{(x,t) : x/t=\xi \}\)). We show how this gives the exact exponential asymptotics of spatio-temporal two-point correlation functions of twist fields, and thus of order-parameter dynamical correlation lengths in homogeneous, stationary states.
In certain situations, we find divergent scaled cumulants, signalling that the fluctuations scale in a different fashion. Seeing the SCGF as a function of the state’s parameters (or of the bias), this may be interpreted as a “dynamical phase transition”, of the type seen in other contexts (see, e.g., [55,56,57,58,59,60,61,62,63]). We explain why this occurs from hydrodynamic principles, and we propose that this may connect with the breaking of Gaussianity found in nonlinear fluctuating hydrodynamics [49, 50, 64].
Of immediate importance is the application to integrable systems, a detailed examination of which is completed in a separate work [65]. There, we exploit generalised hydrodynamics [46, 47, 66,67,68] (experimentally verified [69]) and confirm the proposal in the hard rod gas by comparing with Monte Carlo simulations. In the present paper, we apply the theory to non-integrable quantum critical systems of dimension higher than 1, using conformal hydrodynamics. Non-equilibrium steady states for energy transport were first studied in such systems in [70, 71] and [72,73,74]. Here we obtain exact fluctuation results. We write explicitly the first few scaled cumulants for energy transport as functions of the rest-frame temperature and the relativistic boost, and differential equations for the SCGF, which we solve numerically. We observe a dynamical phase transition in thermal states boosted to the sound velocity.
The paper is organised as follows. In Sect. 2, we explain the context and review the main aspects of large-deviation theory for transport. In Sect. 3, we present our main results, explain the main idea of the derivation, discuss the extended fluctuation relations, give the generalisation to arbitrary rays and the application to twist-field correlation functions, and discuss dynamical phase transitions. In Sect. 4 we present the application to conformal hydrodynamics, and in Sect. 5, we present conclusions and open questions. Finally, in “Appendix A” we provide the main derivation of the general results; in “Appendix B”, we review basic aspects of Euler hydrodynamics, including the solution to the Riemann problem in free (linear) hydrodynamics and the normal modes of conformal hydrodynamics; and in “Appendix C”, we discuss the multi-parameter SCGF and present related general arguments.
2 Maximal-Entropy States and Large-Deviation Theory
In this section, we first describe the context in which the main results are proposed to apply: all systems that possess an Euler hydrodynamic description. We specify what we believe would be the general properties of many-body systems that are expected to be necessary for our results to hold. We then recall the main aspects of large-deviation theory that we need, including the scale cumulant generating function (SCGF).
2.1 Systems and States of Interest
We consider an infinite-length one-dimensional many-body system with a dynamics that is homogeneous in both space and time, and with local interactions. This can be a classical or quantum lattice model, field theory or gas, and the dynamics may or may not be generated by a Hamiltonian. We believe the results apply to both deterministic and stochastic dynamics, although we concentrate on the former,Footnote 1 where randomness lies in the initial state. The assumption of one-dimensionality may be partially lifted by applying the results below to effectively one-dimensional transport in higher-dimensional systems; we discuss this in the context of higher-dimensional conformal hydrodynamics in Sect. 4.
The model is assumed to admit a certain number of homogeneous conserved charges \(Q_i = \int _{{\mathbb {R}}}{\mathrm{d}}x\,\mathcal {q}_i(x,t)\). These are dynamical observables satisfying \({\mathrm{d}}Q_i/{\mathrm{d}}t=0\). They are assumed to have associated conservation laws
indexed by i. Here \(\mathcal {q}_i(x,t)\) and \(\mathcal {j}_i(x,t)\) are the charge density and current, respectively, at space-time point (x, t), which are local (supported on a finite region containing x) or quasi-local (an appropriate extension [43, 75]) observables at that point (we likewise say that \(Q_i\) are local or quasi-local charges). The conservation laws follow purely from the dynamics of the model. If the system is Hamiltonian, then the conserved charges include the Hamiltonian, and we assume them to be in involution (they commute with each other as well as with the Hamiltonian).
States are statistical in nature; fluctuations in deterministic dynamical systems may come from various sources, such as fluctuating initial conditions or quantum fluctuations. Following the general philosophy of the \(C^*\)-algebra description of quantum and classical statistical mechanics, states are fully described by the set of all expectations, here denoted using the bracket notation \(\langle \cdots \rangle \), of local observables; this space may be completed under various norms (see, e.g., [43, 76,77,78]).
We consider the manifold of homogeneous and stationary (i.e., invariant under space and time translations) maximal-entropy states (MES). Formally, MES are characterised by as many Lagrange parameters \(\beta ^i\) (indexed by i) as there are conserved quantities \(Q_i\), and have probability measure or density matrix proportional to
We will denote expectations in a MES by \(\langle \cdots \rangle _{{\underline{\beta }}}\), where \({\underline{\beta }}\) is the vector of all \(\beta ^i\)’s, coordinates for the MES manifold (when it is clearer, we will also use the notation \(\beta ^\bullet \) to represent the set of \(\beta ^i\)’s for all i).
The form \({\mathrm{e}}^{-\sum _i \beta ^i Q_i}\) for the probability measure or density matrix is formal in most situations. In the simple cases where the series \(\sum _i \beta ^i Q_i\) truncates and the charges \(Q_i\) are local, there is a variety of ways to make it rigorous in the context of \(C^*\) algebras: as an infinite-volume limit; via the Kubo–Martin–Schwinger (KMS) relation (in the quantum case), or the Dobrushin–Lanford–Ruelle (DLR) equations (in the classical case); by a precise notion of entropy maximisation; or by considering appropriate tangents to a manifold of states (see [76,77,78] for discussions).
When the series \(\sum _i \beta ^i Q_i\) does not truncate, or the charges involved are quasi-local, a formulation, also based on tangents to a manifold of states, is developed in [43]. It accounts in particular for “generalised thermalisation” in integrable systems. It is based on considering conserved charges \(Q_i\) as vectors in the tangent space to the MES manifold. In this formulation, the state is essentially of the form
for some conserved charge \(Q(\underline{\beta })\), which may be local or quasi-local. The Lagrange parameters \(\beta _i\), associated with various charges \(Q_i\), are defined by the deformation equation
where the symmetric inner product on the space of (quasi-)local fields is
and
refers to the connected, symmetrised correlation function.Footnote 2 According to [43], state (2.3) gives rise to a Hilbert space (which in many cases is expected to be countable-dimensional) induced by (2.5). In order to construct this Hilbert space, one completes the space of local observables with respect to the natural topology from (2.5). This requires the adjunction of elements which are not local in the usual sense, but which still satisfy a weak property of locality. The subspace that is invariant under time evolution is the space of conserved charges. This subspace might be finite- or infinite-dimensional, and the conserved charges might not all be local, depending on what is allowed by the dynamics. It is this subspace that is identified with the tangent space at that point in the MES manifold. In this subspace, one can choose a complete set of conserved charges \(Q_i\), which span a dense subspace. We believe that the results below hold if the set \(Q_i\) is complete in this sense.
In one dimension, away from the ground state, there is no spontaneous symmetry breaking, and thus, MES describe “pure phases”. As such, they are extremal states [76,77,78], and therefore cluster at large distances. Clustering of two-point functions can be shown to be exponential on general grounds in Gibbs states of local quantum Hamiltonians [79]. In states of form (2.3) with Q quasi-local, clustering is still expected to hold. However, because of quasi-locality, it may be weaker than exponential, with “long-range” algebraically decaying correlations; see, for instance, the analysis provided in [17]. In any case, below we assume clustering to be strong enough so that integrals of connected correlation functions converge.
The manifold of MES forms the basis of the emergent hydrodynamics in slowly varying, long-wavelength states [48], in that it gives rise to the “equations of state”: the relation between average currents and average densities.
In generic Galilean invariant quantum and classical gases, the set of conserved quantities contain the particle number, the energy and the momentum, and the MES are simply the Gibbs ensembles and Galilean boosts thereof. Similar statements hold for generic relativistic gases. In these cases, the Euler hydrodynamics is the standard one, Galilean or relativistic. In integrable systems, there are infinitely many conserved quantities, including quasi-local ones [75, 80, 81], and the MES manifold is infinite-dimensional. In these cases, the MES are referred to as generalised Gibbs ensembles [41,42,43,44], and the Euler hydrodynamics, referred to as generalised hydrodynamics, was developed in [46, 47, 66] (and was recently verified experimentally [69]). In a large family of integrable models, the Lagrange parameters \(\beta ^i\)’s are more appropriately represented by a function on a “spectral space”, a space of available stable “quasi-particles”, via the thermodynamic Bethe ansatz [82,83,84] (see, e.g., [75, 85]); thus in these cases, the MES manifold is a manifold of functions, whose full description is, however, not yet known in most models.
The non-equilibrium steady states emerging from the partitioning protocol [26,27,28] also are expected to be MES. In this protocol, two semi-infinite, separate halves (seen as two baths) of the system are initially in different states, often taken to be different thermal states at different temperatures (possibly with different boosts). Suppose the set of Lagrange parameters are \(\beta ^i_{\mathrm{l}}\) (left) and \(\beta ^i_{\mathrm{r}}\) (right). The two halves are then connected to each other and let to evolve for a long time. In any finite region around the connection point, a steady state develops at infinite times, with Lagrange parameters \(\beta ^i\) depending on the initial left and right sets, \(\beta ^i({\underline{\beta }}_{\mathrm{l}},\,{\underline{\beta }}_{\mathrm{r}})\). In particular, if ballistic transport is supported, non-equilibrium currents may emerge. In fluid dynamics, this is known as the Riemann problem [86]. The explicit steady states may be constructed in a number of models (see the review [25] and the papers [46, 47]).
Remark 2.1
MES are essentially the “invariant equilibrium states” as first developed by Gibbs. Gibbs’ theory for gases usually does not include the momentum and its associated intensive “potential” as thermodynamic quantities, but these may be simply re-introduced by Galilean (or relativistic) boosts. Away from the rest frame, the state is no longer time-reversal invariant, hence the appellation “equilibrium state” is not appropriate, and “maximal-entropy states” seems better suited to describe the full set. In more general systems, such as in integrable systems, there are usually infinitely many conserved charges that break time-reversal invariance, which cannot be accounted for simply by Galilean or relativistic boosts. A further difference from conventional Gibbs states is the possibility, as mentioned above, that quasi-local charges involved in a MES give rise to correlations with slower decay than exponential.
Remark 2.2
Physically, the MES manifold is the set of states that are homogeneous, stationary and clustering and that describe averages of observables on finite regions of infinite space emerging after relaxation. That is, these are all states that occur after evolving for a long time from generic initial states in infinite volume, in accordance with local relaxation in isolated, thermodynamically large systems [5]. The set of conservation laws restricts the state’s evolution, and MES are thus states where entropy is maximised with respect to all available local (and quasi-local) conservation laws. An interesting question, which to our knowledge has not been settled, is as to if the set of homogeneous, stationary and clustering states includes other states than MES and their convex combinations. In quantum systems, an argument for the negative is via the (generalised) eigenstate thermalisation hypothesis (ETH): the space of homogeneous, stationary states is spanned by the eigenstates (the diagonal ensemble); by the ETH, this is the convex space of linear combinations of density matrices (2.2) over, generically, many sets of Lagrange parameters; only the extremal points are clustering, thus clustering imposes a single set of Lagrange parameters \(\{\beta _i\}\).
2.2 Large Deviations in Transport
In LDT, one concentrates on fluctuating quantities \(J^{(t)}\), which are extensive with respect to some parameter t, and whose densities \(J^{(t)}/t\) take almost-sure values \(\bar{j}\) in the limit \(t\rightarrow \infty \). For the purpose of studying non-equilibrium transport, it is natural to focus on the LDT of the total transfer of a conserved quantity after time t, between two regions of a system. In higher-dimensional systems, we assume transport to occur in a single direction of space, so that the system is effectively one-dimensional. Transfer occurs say from the left, \(x<0\), to the right, \(x>0\). The conserved quantity Q is one of the \(Q_i\)’s, with index \(i=i_*\),
The total transfer of Q after time t, for the purpose of the LDT, is the total current passing by the origin,
The quantity \(J^{(t)}\) is a random variable, whose probability distribution \({\mathbf {P}}\) is determined by the state of interest, which we take to be a MES for some \({\underline{\beta }}\). According to the large-deviation principle, such extensive quantities have probability distributions that are exponentially peaked at the almost-sure value (here the scaling in time has exponent 1, which is the one relevant here),Footnote 3
The almost-sure value \(\bar{j}\) is simply the average current in the state of interest. The function I(j) controlling this exponential is referred to as the large-deviation rate function. It describes the probabilities of rare but significant events where the quantity \(J^{(t)}\) deviates by large amounts from \(t\bar{j}\). The rate function is the Legendre–Frenchel transform of the scaled cumulant generating function (SCGF) \(F(\lambda )\) for \(J^{(t)}\),
Here \(c_n\) are the cumulants, scaled by time. Up to a conventional minus sign, the function \(F(\lambda )\) can be interpreted as a non-equilibrium equivalent of the equilibrium specific free energy, with non-equilibrium partition function \(Z = \langle {\mathrm{e}}^{\lambda J^{(t)}}\rangle _{{\underline{\beta }}} \asymp {\mathrm{e}}^{tF(\lambda )}\). Crucially, the large-deviation principle (2.9) implies that all cumulants of the random variable \(J^{(t)}\) scale like t at large t. An important question in LDT is the evaluation of the exact large-deviation function I(j), or the exact SCGF \(F(\lambda )\).
Remark 2.3
In the LDT for quantum systems, formulating the problem of the fluctuations of the total transfer of a charge Q by a direct interpretation of the above formulae is not physically natural: the operator \(\int _0^t {\mathrm{d}}s\,\mathcal {j}(0,s)\) is not a natural quantum observable on which von Neumann measurements can be made, as it involves a time integral. Instead, one needs a formulation that takes into account properly the quantum nature of the system (see, e.g., [9]). This can be via a two-time von Neumann measurement procedure where the charge difference between the left and right halves of the system is measured at time 0, the system is let to evolve, and the charge difference is again measured at time t. It can also be via some indirect measurement scheme, for instance where the current passing by the origin is coupled to an external device on which von Neumann measurements are made (see, e.g., [13]). There are indications that suggest that different measurement schemes lead, in the large-deviation limit and for ballistic transport, to the same result, and that this result is in agreement with (2.8), (2.9), and (2.10) (see, e.g., the review [25]). In this paper, we assume this to be the case.
Remark 2.4
The most interesting applications of the LDT for transport are in non-equilibrium steady states. A crucial condition for the theory that we develop below to be applicable is that the MES be established over the full period of time integration in the variable \(J^{(t)}\) in (2.8). Specifically, given some non-equilibrium protocol, if there exist strong correlations between baths, and system at long enough times, these affect the fluctuations of \(J^{(t)}\). Fluctuations of this type depart from those of a MES and invalidate the results of our theory. Importantly, we expect the theory to hold in the partitioning protocol for deterministic systems (see the discussion and numerical results in [65]).
3 Main Results: Ballistic Large-Deviation Theory
In this section, we give the exact expressions for SCGF for ballistic transport in terms of objects from linear fluctuating hydrodynamics. We present an overview of the derivation, with the details given in “Appendix A”. We then provide various consequences, and we extend the results to fluctuations along arbitrary rays, with applications to correlation functions of twist fields.
3.1 Scaled Cumulant Generating Function for Ballistic Transport
It is expected that there is a set of averages of all local or quasi-local densities, which we will denote by \({\mathsf {q}}_i = \langle \mathcal {q}_i(0,0)\rangle _{{\underline{\beta }}}\), which provide a good system of coordinates for the MES manifold—that is, the map \({\underline{\beta }} \mapsto \underline{{\mathsf {q}}}\) is bijective (from an appropriate space of \({\underline{\beta }}\)). In (a large family of) integrable systems, this set is, again, more appropriately represented by a function on a spectral space (the “quasi-particle density” of the thermodynamic Bethe ansatz [82,83,84]).
Consider, then, the averages of the currents, \({\mathsf {j}}_i = \langle \mathcal {j}_i(0,0)\rangle _{\underline{\beta }}\), as functions of the state coordinates \(\underline{{\mathsf {q}}}\). These model-dependent functions—the fluxes—are the equations of state of the model. Construct the flux Jacobian (or linearised Euler matrix)
This is a (model-specific) matrix that is a function of the MES, hence a function of the coordinates \(\underline{\beta }\), or \(\underline{{\mathsf {q}}}\). The flux Jacobian is at the basis of the Euler hydrodynamic theory (see “Appendix B”) and is a fundamental part of what is often referred to as linear fluctuating hydrodynamics.
Consider the MES characterised by some coordinates \({\underline{\beta }}\). We define a flow \(\lambda \mapsto {\underline{\beta }}(\lambda )\) on the manifold of MES, starting on this state \({\underline{\beta }}(0)={\underline{\beta }}\) and for \(\lambda \) lying in some interval of \({{\mathbb {R}}}\), by the differential equation
where \({\mathsf {A}}(\lambda )\) is the flux Jacobian in the state with Lagrange parameters \({\underline{\beta }}(\lambda )\). On the right-hand side, the sign of the flux Jacobian \({\mathsf {A}}\) is the matrix obtained by diagonalising \({\mathsf {A}}\) and taking the sign of its eigenvalues,
where \(v^{\mathrm{eff}}=\mathrm{diag}(v^{\mathrm{eff}}_1,v^{\mathrm{eff}}_2,\ldots )\) is the diagonal matrix of eigenvalues [the \(v^{\mathrm{eff}}_i\)’s are the effective velocities attached to the normal modes of the hydrodynamics (see (B.6))]. Recall that \(i_*\) is the index corresponding to the conserved charge of interest (Eq. (2.7)).
By the chain rule, Eq. (3.2) is equivalent to the following flow in the conserved density coordinates,
The static correlation matrix is defined by
and \({\mathsf {C}}(\lambda )\) in (3.4) is the flow-dependent static correlation matrix, evaluated in the state with Lagrange parameters \({\underline{\beta }}(\lambda )\).
Our main result, shown in “Appendix A”, is as follows. We identify the flow parameter \(\lambda \) with the conjugate parameter in (2.10), and we have an expression for the SCGF \(F(\lambda )\) directly in terms of the current along the flow:
where \({\mathsf {j}}(\lambda ) = {\mathsf {j}}_{i_*}(\lambda ) = \langle \mathcal {j}_{i_*}(0,0)\rangle _{\underline{\beta }(\lambda )}\). Thus, the knowledge of the Euler hydrodynamics (giving the flux Jacobian \({\mathsf {A}}\), and the currents \({\mathsf {j}}_i\), as functions of the state) is sufficient in order to obtain the SCGF. As we also show in “Appendix A.4”, this agrees with, and largely generalises, the result for the second cumulant \(c_2\) which follows from the current–current sum rule written in [50, 51].
If \(F(\lambda )\) is strictly convex and everywhere differentiable, then the Legendre–Frenchel transform reduces to the Legendre transform, and it is a simple matter to obtain the large-deviation function as:
Remark 3.1
One can, formally, generalise to SCGFs with multiple parameters \(\lambda ^j\) associated with all currents \(\mathcal {j}_j\), with
In (3.4) and (3.2), we make the replacements \(\lambda {\mapsto }\lambda ^j\) and \(i_*{\mapsto } j\). This, however, requires multiple differentiability with respect to the parameters \(\lambda _j\), which are nontrivial relations on the matrix \({\mathsf {A}}\) as a function of these parameters. In this paper, we do not investigate this aspect; however, see “Appendix C”.
Remark 3.2
In certain cases, where the generating function of the currents separates into a sum of functions of the normal modes, we obtain a more explicit expression for \(F(\lambda )\), developed in “Appendix C”. This agrees with the general expression found in integrable systems (see [65]).
3.2 Derivation: Biasing the Measure
The derivation of the main results (3.6) with (3.2) is provided in “Appendix A”. The main argument is to bias the measure in a particular way and show that the bias generates a flow on the MES manifold. This latter fact can be shown either using the theory of pseudolocal charges [43, 75], or from a strong version of the hydrodynamic projection principles [48,49,50,51]. In order to specify the explicit flow in terms of the flux Jacobian, one needs certain basic results from linear fluctuating hydrodynamics. Here we give the main lines of the proof. We note that the basic techniques leading to form (3.6) for the SCGF were introduced in [23, 24] in the context of \(1+1\)-dimensional conformal field theory and more generally in [21]; in fact, result (3.6) with (3.2) may be seen as a nonlinear generalisation of the results found there, as is made clearer in the next subsection.
Besides assumptions which are expected to be valid quite generally in local many-body systems, the leading assumption of physical relevance is that of sufficiently strong clustering of local observables at long times. That is, multi-point connected correlation functions of local fields, in particular of local currents, vanish at large time separations, in a way that makes them integrable.
Contrary to the case of clustering of correlation functions at large space separations in MES, we do not know of any strong mathematical result on clustering of correlation functions in time. However, some insights can be gained from fluctuating hydrodynamics (see [49, 50, 64]). There, it is phenomenologically argued that correlation functions will cluster exponentially along any ray, in space-time, which is not co-propagating with one of the fluid’s normal modes. Thus, strong clustering of correlation functions in time should occur if there are no normal modes with vanishing velocity. The breaking of this condition leads to a failure of the large-deviation principle (2.9). We discuss the physics and potential consequences of this in Sect. 3.5.
Let us modify the measure for the state \(\langle \cdots \rangle _{{\underline{\beta }}}\) by a time integral of the current \(\mathcal {j}(0,t)\). That is, let us construct a family of states \(\langle \cdots \rangle ^{(\lambda )}\), parametrised by \(\lambda \in {{\mathbb {R}}}\), with \(\langle \cdots \rangle ^{(0)}=\langle \cdots \rangle _{{\underline{\beta }}}\), obtained by modifying the state \(\langle \cdots \rangle _{{\underline{\beta }}}\) by the insertion of the time-integrated local current \(\mathcal {j}(0,t)\) of the charge Q at the origin,
(the symmetrisation guarantees that averages of Hermitian observables are real numbers in the quantum case). This is well defined as a formal expansion in \(\lambda \) if connected correlation functions vanish fast enough at large time separation, and it is in fact expected to be well defined for real values of \(\lambda \) in an interval containing the origin. In particular, we have
where the connected correlation function \(\langle \cdot ,\cdot \rangle ^{(\lambda ), \mathrm c}\) is defined as in (2.6) but for the state \(\langle \cdots \rangle ^{(\lambda )}\).
As a loose interpretation, the insertion of the exponential of the time integral of the current can be seen as “biasing” the dynamics, changing the weights of trajectories in order to make rare events “typical” and access their probabilities. The biasing by a time-integrated current is natural and has been used widely in the study of large deviations in stochastic dynamics or open quantum systems. In this context, one attempts to relate it to a change of the stochastic dynamics or of the Lindbladian, something referred to as the generalised (classical or quantum) Doob transformation [52, 53]. By contrast, here we relate it to a change of the distribution ruling the initial state, and crucially, the new distribution is still a MES—we obtain a function \({\underline{\beta }}(\lambda )\) with \({\underline{\beta }}(0)={\underline{\beta }}\). We determine this change solely from the Euler hydrodynamics of the system.
Before determining \({\underline{\beta }}(\lambda )\), we explain how the \(\lambda \)-dependent state is useful in order to fix \(F(\lambda )\) (see, e.g., [21]). Taking definition (2.10), we find, after a convenient shift in the time integration region allowed by stationarity of the state and after symmetrisation,
Under an appropriate assumption of sufficiently fast clustering in time, when t is large we expect the main contribution to the s integral to come from the central region away from the boundaries \(s=\pm t/2\), where the state is stationary. That is, we may write
where we used (3.9). The first equality is because the main contribution is from the central region and the second because this contribution is time independent. Integrating on \(\lambda \) with the condition \(F(0)=0\), we obtain (3.6).
The state resulting from the \(\lambda \)-bias is manifestly stationary. Also, by the fact that \(\mathcal {j}(x,t)\) is part of a conservation law, and by clustering at large time differences (uniformly on finite intervals of position differences), we have \(\int _{-\infty }^\infty {\mathrm{d}}t\,\mathcal {j}(0,t) = \int _{-\infty }^\infty {\mathrm{d}}t\,\mathcal {j}(x,t) + \lim _{t\rightarrow \infty } \int _0^x {\mathrm{d}}y\,( \mathcal {q}(y,t)- \mathcal {q}(y,-t))= \int _{-\infty }^\infty {\mathrm{d}}t\,\mathcal {j}(x,t)\) for any x (where the last equality holds inside connected correlation functions), whence the resulting state is homogeneous. It is also possible to argue that the state is clustering. Is this state a MES (see Remark 2.2)? We show in “Appendix A” that indeed it is. We show in two ways—either from the theory of pseudolocal charges or from hydrodynamic projection principles—that infinitesimal \(\lambda \) modifications, Eq. (3.10), generate tangents to the MES manifold. Since at \(\lambda =0\) the state lies on the MES manifold, then it stays on it. As a consequence, \(\lambda \mapsto \langle \cdots \rangle ^{(\lambda )}\) forms a path lying within this manifold. Therefore, there exists \({\underline{\beta }}(\lambda )\) with \({\underline{\beta }}(0)={\underline{\beta }}\) and
This is the crucial observation of the method.
Using (3.10), we then have
for any local or quasi-local observable \({{\mathcal {O}}}\). We may obtain an equation for the coordinates \({\mathsf {q}}_i(\lambda )\) by specifying \({{\mathcal {O}}}\) in (3.14) to be the available conserved densities of the model, \(\mathcal {q}_i(0,0)\). The left-hand side is therefore the time derivative of the state coordinates. The right-hand side is a time-integrated two-point function of conserved densities and currents, and this is a function of the state and hence can be seen as a function of the state coordinates \(\underline{{\mathsf {q}}}(\lambda )\). Equation (3.14) therefore fully specifies the path by giving its tangent at the point \(\underline{{\mathsf {q}}}(\lambda )\) in terms of a function of \(\underline{{\mathsf {q}}}(\lambda )\). A result from linear fluctuating hydrodynamics [48,49,50,51] is that, in an appropriate Euler scaling limit [48, 51, 87] (see “Appendix A”), where in particular both x and t are large in fixed ratio, we have
where \(\delta (x-\mathsf At) = {\mathsf {M}}\delta (x-v^{\mathrm{eff}}t){\mathsf {M}}^{-1}\) [(see (3.3)]. Integrating over time, the result is independent of x, and we obtain
This is the crucial technical step in the derivation. Equation (3.16), and a generalisation of it necessary in order to show that the \(\lambda \)-bias keeps the state within the MES manifold, is shown more rigorously in “Appendix A” from hydrodynamic projection. Combining (3.16) and (3.14), we indeed find (3.4), which implies (3.2).
3.3 Constant Flux Jacobian and Extended Fluctuation Relations
The form of the flux Jacobian \({\mathsf {A}}_i^{~ j}\) depends on the fluid coordinate system chosen. As the name suggests, the flux Jacobian transforms as a Jacobian: a covariant (contravariant) vector in its first (second) index. Seen as a matrix, this is a similarity transformation, which is generically coordinate dependent. Therefore, the only coordinate-independent information within the flux Jacobian is its spectrum, the elements of the diagonal matrix \(v^{\mathrm{eff}}\) in (B.6) [and in (3.3)]. However, there is more information about the physical system within the flux Jacobian. Indeed, the physical system provides a favoured, special set of coordinate systems: the densities of conserved charges. These are specified by the model up to \({{\mathbb {R}}}\)-linear transformations, but \({{\mathbb {R}}}\)-linear transformations form a subset of the set of all coordinate transformations. Hence, one can define a “natural” flux Jacobian as the flux Jacobian in a system of coordinates given by the conserved densities. This is unique up to \({{\mathbb {R}}}\)-linear transformations, which are coordinate-independent similarity transformations.
In some cases, for instance in non-interacting models and \(1+1\)-dimensional conformal field theory, the natural flux Jacobian is independent of the state (a property which is indeed invariant under \({{\mathbb {R}}}\)-linear transformations). Equivalently, the Euler hydrodynamic equations are linear (see (B.3)). In many ways, this can be considered as a hydrodynamic system without interactions.
In such cases, it is a simple matter to solve for flow (3.2):
That is, the flow corresponds to a shift of the Lagrange parameters proportional to \(\lambda \). In particular, we have from (3.6)
Let us consider non-equilibrium steady states emerging from the partitioning protocol, or the Riemann problem (see the discussion at the end of Sect. 2.1). Consider the Lagrange parameters \(\beta ^i\) characterising the steady state, as functions of the set of Lagrange parameters on the left and right halves of the initial state, \(\beta ^i({\underline{\beta }}_{\mathrm{l}},\,{\underline{\beta }}_{\mathrm{r}})\).
In the cases of a natural flux Jacobian that is independent of the state, it is a simple matter to solve the Riemann problem and to evaluate the steady state in the region around the connection point. We show in “Appendix B.2” that this solution leads to the relation
This, combined with (3.18), is the fully general statement of the extended fluctuation relations, first introduced by Bernard and Doyon [21]. That is, according to the extended fluctuation relations, the biasing of the measure necessary to generate transport cumulants can be performed by linear shifts of the Lagrange parameters in the initial baths of the partitioning protocol. Such linear shifts generate fluctuation statistics of the initial state, hence the extended fluctuation relations indicate that, in free models, the statistics of transport fluctuations is directly obtained from that of the initial-state fluctuations. This appears to be physically sensible, as without interactions, initial-state fluctuations are not affected during transport.
The statement of the extended fluctuation relations [21] was obtained by extracting principles found in [23, 24] in the context of energy and charge transport in \(1+1\)-dimensional CFT, and was argued to hold also in free-particle models, later confirmed by various explicit calculations [16, 17, 19]. This shows that the present formalism fully agrees with these results and that, effectively, it generalises the method to nonlinear Euler hydrodynamics. As a consistency check, it is also a simple matter to see that the solution presented in [65], for interacting integrable models (and based on the present formalism), indeed reproduces the extended fluctuation relations when specialised to models without interactions.
3.4 Fluctuations Along Rays and Dynamical Correlation Functions of Twist Fields
The proposal of Sect. 3.1 can be generalised to the statistics of the component of currents perpendicular to other space-time paths instead of the time-directed paths with constant space coordinates. Of particular interest is the application to dynamical correlation functions of twist fields, including order and disorder fields, in thermal states and other MES.
Consider, instead of \(J^{(t)}\) defined in (2.8), the quantity
determined by the path \(\vec \ell = \{s\mapsto \vec {\ell }(s) = (\hat{x}(s),\hat{t}(s))\in {{\mathbb {R}}}^2:s\in [0,1]\}\), where \({\mathrm{d}}\vec {\ell }(s) = ({\mathrm{d}}\hat{x}(s),{\mathrm{d}}\hat{t}(s))={\mathrm{d}}s\,({\mathrm{d}}\hat{x}/{\mathrm{d}}s,{\mathrm{d}}\hat{t}/ {\mathrm{d}}s)\) is the infinitesimal tangent to the path, \(\vec {\mathcal {j}} = (\mathcal {j},\mathcal {q})\) is the conserved current vector, and
Suppose (without loss of generality) that the path has end points \(\vec \ell (0) = (0,0)\) and \(\vec \ell (1) = (x,t)\), and denote by \(\ell = \sqrt{x^2+t^2}\) the Euclidean distance between the end points. By current conservation, quantity (3.20) is independent of the path chosen that connects (0, 0) to (x, t). We may therefore choose it to be the segment of ray \(x/t=\tan \theta \) (with \(\theta \in [0,2\pi )\)) determined by
We are interested again in the large-\(\ell \), scaled statistics, and so we must evaluate the expectation value
and the generator
Extending the arguments presented in “Appendix A”, or the simpler derivation of Sect. 3.2, it is a simple matter to derive, for path (3.22), the flow [recall (2.7)]
which generalises (3.2) to the case \(\theta \ne 0\). With this flow, the result takes a form that generalises (3.6),
Consider the limit \(\theta \rightarrow \pi /2\), where the path is “horizontal”, lying on a the time slice \(t=0\). In this case, flow (3.25) does not depend on the flux Jacobian anymore and is immediately solvable, simply effecting a shift of the Lagrange parameter \(\beta ^{i_*}\) proportional to \(\lambda \), that is,
From (3.26), we have in this case \(\partial F(\lambda ;\pi /2)/\partial \lambda = -{\mathsf {q}}(\lambda ;\pi /2)\), which, along with \(F(0;\pi /2)=0\) and (3.27), allows us to identify \(F(\lambda ;\pi /2)\) with a free energy difference. That is, we obtain
where \(\Delta f(\lambda )\) is the specific (dimensionless) free energy difference
with the specific free energy being \(f({\underline{\beta }}) = -\log Z({\underline{\beta }})\), where \(Z({\underline{\beta }})\) is the partition function for Lagrange parameters \({\underline{\beta }}\).
The above results have perhaps their most interesting application to the evaluation of correlation functions of twist fields. Let us introduce the “height fields” \(\varphi (x,t)\) defined as \(\mathcal {q}(x,t) = \partial _x \varphi (x,t)\) and \(\mathcal {j}(x,t) = -\partial _t \varphi (x,t)\), which automatically solves the continuity relation (2.1). Differences of height fields \(\varphi (x_2,t) - \varphi (x_1,t)\) count (for \(x_2>x_1\), say) the quantity of charge present in \([x_1,x_2]\) at time t. Exponential of height fields \({{\mathcal {O}}}(x,t) \propto e^{\lambda \varphi (x,t)}\) are a certain type of fields that have been studied in a variety of cases in the literature, and are referred to as twist fields.Footnote 4 If the charge Q is associated with an internal symmetry, then they are local in the general sense used in many-body quantum physics (they commute with the energy density at equal times). In particular, with U(1) symmetry, the observable \(e^{\lambda \varphi (x,t)}\) can be used to represent order parameters in many-body models: in the free Dirac fermion and in the Thirring model (or sine-Gordon model), it naturally occurs by bosonisation and has applications to the transverse field Ising model and XXZ chains (see, e.g., [88, 89]). Certain classes of such twist fields can also be used to study entanglement entropy in free-particle models [90]. Twist fields associated with space-time symmetries have also been studied recently [91]; they do not possess the conventional locality property of many-body physics anymore, although there is still path independence (a field \(\varphi (x,t)\) can be defined independently of the path \(\vec \ell \) chosen).
Interestingly, SCGF (2.10) gives rise to the leading exponential behaviour of the dynamical two-point functions of twist fields:
(recall footnote 3). This is expected to hold in arbitrary maximal-entropy states of arbitrary many-body systems, integrable or not, including thermal states and, in integrable systems, GGEs. The evaluation of exponential behaviours of dynamical two-point correlation functions of order parameters—which give a “dynamical correlation lengths”—is a notoriously difficult problem, for which there are only partial solutions (see, e.g., [89]). Formula (3.30) provides the first exact result in interacting models.
For equal time, non-dynamical correlations (\(\theta =\pi /2\)), (3.30) with (3.28) givesFootnote 5
In this case, some exact results already exist that confirm (3.31). A specialisation to the \({{\mathbb {Z}}}_2\) twist fields of free Majorana fermions of (3.31) was proposed for arbitrary GGEs in [92], and shown to agree with results derived from special quantum quenches in the Ising model [93]. In the context of entanglement entropy, taking into account the twist field interpretation of the Rényi entanglement entropy as expressed in [90], the exact results of [94,95,96,97,98] can be interpreted as giving the leading exponential decay of permutation twist fields correlation functions in GGEs of interacting integrable systems, which also agree withFootnote 6 (3.31).
Remark 3.3
From the path independence of variable (3.20), we may choose different paths than ray (3.22) and try to evaluate the related flow. For instance, we may choose the piecewise straight path which goes first in the time direction \((\hat{x}(s),\hat{t}(s)) = (0,2st),\;s\in [0,1/2]\), and then in the space direction \((\hat{x}(s),\hat{t}(s)) = ((2s-1)x,t),\;s\in [1/2,1]\). Evaluating the flow equation from this is, however, more complicated. It might be tempting to think that the result will separate into two contributions, one from each straight piece, \(F(\lambda ;\theta )\) being proportional to a sum of \(F(\lambda ;0) = F(\lambda )\) (the SCGF calculated in Sect. 3.1) and \(F(\lambda ;\pi /2) = -\Delta f(\lambda )\) (the difference of free energy densities). However, this is generically incorrect: for instance, because of ballistic transport, the space integral \(\int _0^x\,{\mathrm{d}}\hat{x}\,\langle \mathcal {q}(\hat{x},t),\mathcal {q}_i(0,0)\rangle _{{\underline{\beta }}}^{\mathrm{c}}\) does not necessarily vanish in the limit \(x,t\rightarrow \infty \). Likewise, the evaluation of \(F(\lambda ;\theta )\) for any choice of path that is not straight (at the scale set by \(\ell \rightarrow \infty \)) may receive contributions from correlations between separated portions of the path if they are connected by ballistic transport of normal modes. See Sect. 3.5.
Remark 3.4
It is a simple matter to generalise (3.28) and (3.31) to states which are inhomogeneous at Euler scales, and described by fluid cells with position-dependent Lagrange parameters \({\underline{\beta }}(\hat{x})\). By using the idea of local entropy maximisation, one simply expects each “fluid cell” at \(\hat{x}\) to produce a contribution proportional to its free energy difference \(\Delta f_{\hat{x}}(\lambda )\), and thus
An ansatz of this form appeared, for the order-parameter correlation functions in the Ising model, in [99], and was verified against direct numerical calculations. This is also connected to formulae for entanglement entropies in inhomogeneous states [100]. Again, this is expected to be valid for arbitrary twist fields and in interacting models as well. However, for dynamical correlation functions in inhomogeneous, non-stationary states, because of correlations produced by ballistically propagating modes, we do not expect the simple generalisation of (3.26) and (3.30) to similar integrals over space-time paths to be correct. The theory developed in [87] for charge-density dynamical correlations in inhomogeneous, non-stationary states might be useful for solving this problem.
3.5 Divergence of Scaled Cumulants and Non-Gaussianity
We now explain in what situations the limit defining the scaled cumulants in SCGF (2.10) may be divergent [in which case the large-deviation principle expressed in (2.9) fails], and what the meaning of this may be.
We first note that the scaled cumulants of (2.10) are time-integrated, connected, multi-point correlation functions of local current observables, as expressed in (A.1). They exist if correlation functions of local currents cluster fast enough at large time separations. More generally, the cumulants on arbitrary rays, (3.23) and (3.24) defined in Sect. 3.4, exist if correlation functions cluster fast enough at large separations along rays in space-time. This is the main assumption behind the results we have presented: that of strong enough clustering. Clustering in space can be shown on quite general grounds in Gibbs states of local Hamiltonians in one dimension [79] (in higher dimensions, it may fail as thermal phase transitions are possible). However, clustering along nontrivial rays is more subtle. In what situations may it be broken?
Consider, for instance, the current–density correlation. In the Euler scaling limit (see “Appendix A”), it can be expressed (formally) as (3.15). It is clear that this is zero at \(x/t=\xi \) if there is no effective velocity \(v^{\mathrm{eff}}_i\) that takes the value \(\xi \). This means that, far along this ray, the current–density connected correlation function is expected to vanish exponentially fast. Likewise, current–current correlation functions,
vanish exponentially fast if \(v^{\mathrm{eff}}_i\ne \xi \;\forall \;i\). In fact, from the strong hydrodynamic projection principle (A.8), and the slightly stronger version [87, Eq 3.35], this holds more generally for local observables. However, if there is an effective velocity at the value \(\xi \), then the Euler-scale expression diverges. This is generically associated with algebraic instead of exponential clustering of correlation functions.
There is therefore a link between the strength of the correlation on the ray \(\xi \) and the presence or not of an effective velocity with value \(\xi \). Recall that the effective velocity is the velocity of ballistic propagation of the normal modes of the fluid (such as pressure waves in air—sound waves). It is physically natural that ballistically propagating normal modes create strong correlations along their paths. These are sometimes referred to as “sound peaks”, or “heat peaks”. Such strong correlation also occurs naturally in rarefaction waves [86]: there, the state at ray \(\xi \) is such that \(v^{\mathrm{eff}}_j=\xi \) for some j. As a consequence, the large-deviation principle (2.9) does not hold if there exists such a normal mode propagating along the ray—a “co-propagating mode”. In this situation, cumulants may be divergent.
In fact, it is possible to argue from our explicit results that co-propagating modes have a strong effect. The calculation of cumulants requires us to take \(\lambda \) derivatives. Consider the case \(\xi =0\) for simplicity. Clearly, from (3.6), the first cumulant (the average current) is expected to be finite. For higher cumulants, the flow equation (3.2) is ambiguous at \(\lambda =0\) if there is a co-propagating mode (an eigenvalue of \({\mathsf {A}}\) that vanishes), because of the ambiguity of the sign function at 0. Despite this, the second cumulant can be evaluated explicitly and its final expression (A.34) is unambiguous (due to continuity of the absolute value). However, this phenomenon does not hold at higher orders: the sign function and its derivatives must be evaluated at 0, giving ambiguities and even apparent divergencies, because naively the sign function has a divergent derivative at 0. Therefore, the present theory does not provide access to cumulants when there are co-propagating modes, but does indicate that divergencies may occur.
In fact, the special points—either in parameter space or in \(\lambda \) space—where co-propagating modes exist may be interpreted as a points where a “dynamical phase transition” occurs. The concept of dynamical phase transitions is widely studied in the non-equilibrium large-deviation theory of stochastic models (see, for instance, [55,56,57,58,59,60,61,62,63]). Our theory gives explicit results for such phase transitions. Indeed, it predicts a sudden change in the fluctuation spectrum for a dynamical quantity as a parameter of the state is modified, or as the bias \(\lambda \) is modified, through such points. On both sides of the phase transition point, the cumulants will take different forms, as \({{\,\mathrm{sgn}\,}}(A)\) changes discontinuously. See, for instance, our results in higher-dimensional conformal field theory (4.21), (4.22).
One subtle point is worth mentioning. It seems as though in order for this phase transition to occur, the co-propagating mode should be “isolated”. We define a non-isolated mode as a mode whose effective velocity is part of a continuum of effective velocities, in a state where the associated modes are smoothly populated and smoothly coupled to the charge whose transport we study. From the results [65], it appears as though if a mode is non-isolated, then no discontinuity appears due to the ensuing smoothing, and no divergence emerges. Indeed, this is what we observe for transport of generic local conserved quantities in integrable systems [65], where there is a continuum of quasi-particle velocities. There are situations where isolated effective velocities may be present in integrable systems, for instance for spin transport in the XXZ spin chain (see [101]), and it might be possible to study transport of charges that couple to a single quasi-particle velocity.
What does the divergence of scaled cumulants mean? We propose that the leading Gaussian form of the fluctuation spectrum may be broken in states with a co-propagating mode. Recall that the large-deviation principle—which says that all cumulants scale with t in (2.10)—is an extension of the law of large numbers: subtracting the average \(t\bar{j}\) and the fluctuations of \((J^{(t)}-t\bar{j})/\sqrt{t}\) are Gaussian at large t, with nonzero second cumulant and vanishing higher cumulants [as the higher cumulants of this variable receive a scaling \(t^{-n/2}\) instead of the \(t^{-1}\) used in (2.10)]. The divergence of the scaled cumulants of large-deviation theory \(c_n,\;n\ge 3\) suggests that the cumulants of \((J^{(t)}-t\bar{j})/\sqrt{t}\) might no longer be vanishing, thus breaking Gaussianity (of course, the cumulants of this quantity, including the second, might simply not exist, another way of breaking Gaussianity).
Crucially, this may have a connection with nonlinear fluctuating hydrodynamics [49, 50, 64]. Nonlinear fluctuating hydrodynamics can be used to describe the broadening of correlation peaks occurring along the ballistic rays of normal modes, and correlations in rarefactions waves. It is observed that the first-order (linear) expansion of noisy hydrodynamic equations, leading to Gaussian fluctuations, vanishes if there is a co-propagating mode, and the next order needs to be taken, leading to fluctuations in the KPZ class. Satisfactorily, the present theory also has peculiarities when co-propagating modes are present that point towards the breaking of Gaussianity. In particular, with [65] and the above discussion, this suggests that, as already predicted in [49], the KPZ fluctuations of nonlinear fluctuating hydrodynamics do not generically occur in integrable systems, except perhaps for very specific variables and states, such as the spin in spin transport problems of the Heisenberg chain (as [101] indicates that there is an isolated mode),Footnote 7 or perhaps observables able to isolate quasi-particle velocities. The present theory does not yet confirm the KPZ class of fluctuations predicted by nonlinear fluctuating hydrodynamics. We hope to develop these ideas in a future work.
4 Application to Conformal Hydrodynamics in Arbitrary Dimensions
The proposal of Sect. 3.1 applies to all systems in the class described in Sect. 2.1. Although the focus was on one dimension, the formalism applies as well to effectively one-dimensional set-ups in higher-dimensional systems. The goal of this section is to provide a non-integrable example of the formalism and to show how dimensional reduction is performed.
There are many non-integrable systems which admit ballistic transport. One family of examples is relativistic or Galilean quantum and classical field theory, in arbitrary dimension. Higher-dimensional relativistic conformal field theory is particularly interesting, as it makes predictions for quantum systems tuned to quantum critical points at small but nonzero temperatures. In dimensions higher than one, very few results are available. Further, the equations of state—giving conformal hydrodynamics—are almost completely determined by the symmetries, making the present formalism immediately applicable.
Here we study the important example of energy transport in conformal hydrodynamics of arbitrary dimension, obtaining explicitly flow (3.2), expressions for the cumulants \(c_2\), \(c_3\) and \(c_4\), and numerically evaluating the SCGF \(F(\lambda )\). All these are, to our knowledge, new results.Footnote 8 We further note that conformal hydrodynamics was used in order to obtain exact non-equilibrium steady states of quantum critical systems in the partitioning protocol [70, 71, 73, 74]. The properties of \(F(\lambda )\) in such states will be analysed in a future publication.
4.1 Reduction to a One-Dimensional Hydrodynamic Problem
Consider conformal field theory in \(d> 1\) dimensions of space, and recall that the energy–momentum tensor \(\mathcal {T}^{\mu \nu }\) satisfies \(\mathcal {T}^{\mu \nu } = \mathcal {T}^{\nu \mu }\) (Lorentz invariance), \(\mathcal {T}^\mu _{\ \mu }=0\) (scale invariance) and \(\partial _\mu \mathcal {T}^{\mu \nu } = 0\) (conservation of energy and momentum). Assume the system not to be integrable—this is the generic situation and excludes free field theory. Then, the full MES manifold is that of boosted thermal states. For simplicity, consider transport in the direction \(x^1\) (with \(x^0\) the time coordinate), and the associated momentum operator \(P = \int {\mathrm{d}}^{d} x\,\mathcal {T}^{01}(x)\). In this case, we may restrict to the space of thermal states boosted in that direction [70], with density matrices
and corresponding state denoted by \(\langle \cdots \rangle ^{(d)}_{\beta _1,\beta _2}\). Here \(H = \int {\mathrm{d}}^{d} x\,\mathcal {T}^{00}(x)\) is the Hamiltonian, the rest-frame temperature is \(T_\mathrm{rest}=\beta _\mathrm{rest}^{-1}\), and the Lorentz boost is of rapidity \(\theta \). By relativistic and conformal invariance, expectation values of energy–momentum components take the form
where \(u^\mu = (\cosh \theta , \sinh \theta , 0,\ldots ,0)^\mu \) and a is a model-dependent positive constant.
In order to make the connection with the formalism developed, we need to render the system effectively one-dimensional. This can be done by integration over the transverse directions. Specifically, we assume the transverse space \(S_\perp \), with coordinates \(x^\perp = (x^2,\ldots ,x^d)\), to be of \(d-1\)-dimensional hyperarea \(V_\perp \) and to be periodic in all its coordinates, with equal periods. Let us denote by \(\langle \cdots \rangle ^{(d,V_\perp )}_{\beta _1,\beta _2}\) the resulting state with density matrix of form (4.1); in particular, the limit of infinite transverse hyperarea reproduces the infinite-volume results,
We then define one-dimensional densities as
(the equality \(\mathcal {j}_1^{(V_\perp )}(x,t) = \mathcal {q}_2^{(V_\perp )}(x,t)\) comes from \(\mathcal {T}^{01}(x,x^\perp ,t) = \mathcal {T}^{10}(x,x^\perp ,t)\), due to Lorentz invariance). It is a simple matter to see that
and correlation functions cluster at large longitudinal distances. Therefore this is an effectively one-dimensional system, with two conserved charges, and assuming that the transverse direction does not give rise to additional thermodynamic degrees of freedom (this in particular assumes no turbulent instabilities (see the discussion in [70])), one-dimensional Euler hydrodynamics apply. We may then ask about the SCGF for the energy current \(\mathcal {j}_1^{(V_\perp )}\) as defined in (2.10), and the general discussion and results of Sect. 3 hold. In fact, it is convenient to divide the SCGF by \(V_\perp \), and so we consider
By clustering in \(d+1\)-dimensional space-time, all cumulants generated by \(F(\lambda )\) have a finite limit as \(V_\perp \rightarrow \infty \), and we look for
Clearly, this is obtained by taking the large-\(V_\perp \) limit of the solution presented in Sect. 3.1 for the effectively one-dimensional system, that is, of expression (3.6) divided by \(V_\perp \), with (3.2). It is then sufficient to know the following large-\(V_\perp \) limits of conserved densities and currents:
and to solve for flow (3.2) with the flux Jacobian given by
Using (4.3) as well as homogeneity in the transverse direction, the limits in (4.8) are given exactly by the expression on the right-hand side of (4.2) for \(\mu ,\nu \in \{0,1\}\). Thus, we have fully reduced the problem of evaluating (4.7) to a one-dimensional hydrodynamic problem.
4.2 Exact SCGF and Cumulants
Using (4.2), we have more explicitly
The flux Jacobian takes the form
As the flux Jacobian is not state independent if \(d>1\), the extended fluctuation relations do not hold in higher-dimensional CFT.Footnote 9 The matrix \({{\,\mathrm{sgn}\,}}({\mathsf {A}})\) is obtained by diagonalising \({\mathsf {A}}\) and taking the sign of the eigenvalues. The eigenvalues of \({\mathsf {A}}\) are:
with the speed of sound of conformal hydrodynamics given by \(v_{\mathrm{s}} = 1/{\sqrt{d}}\). Clearly, \({{\,\mathrm{sgn}\,}}(v^{\mathrm{eff}}_\pm ) = {{\,\mathrm{sgn}\,}}(\sinh 2\theta \pm \alpha ) = {{\,\mathrm{sgn}\,}}(\theta \pm \theta _{\mathrm{s}})\) with the sound rapidity \(\theta _{\mathrm{s}}\) defined by \(\tanh \theta _{\mathrm{s}} = v_{\mathrm{s}}\). Define \( \xi _1 = {{\,\mathrm{sgn}\,}}{v^{\mathrm{eff}}_+} + {{\,\mathrm{sgn}\,}}{v^{\mathrm{eff}}_-}\) and \( \xi _2 = {{\,\mathrm{sgn}\,}}{v^{\mathrm{eff}}_+} - {{\,\mathrm{sgn}\,}}{v^{\mathrm{eff}}_-}\). They take the following values:
That is, \(\xi _1\) is nonzero for supersonic rapidities and zero otherwise, and \(\xi _2\) is nonzero for infrasonic rapidities and zero otherwise. We then have
Consider the SCGF for energy transport. We fix \(i_{*} = 1\) in (3.2) and let the Lagrange multipliers become \(\lambda \) dependent. Then using \(\beta _1 = \beta _\mathrm{rest}\cosh \theta \) and \(\beta _2 =- \beta _\mathrm{rest}\sinh \theta \), we obtain
We may evaluate the SCGF according to (3.6) by integrating current (4.11)
Cumulants are simply obtained by taking derivatives with respect to \(\lambda \) and setting \(\lambda =0\).
The derivatives may be readily evaluated:
Reading off the required identities from (4.17), setting \(\lambda = 0\) and using the explicit form of \(\xi _{1/2}\) given by (4.15), we obtain the second cumulant:
where \(\theta = \theta (0)\) and \(\beta _\mathrm{rest}= \beta _\mathrm{rest}(0)\). This process easily generates cumulants of \(n^{\text {th}}\) order. As an example, using the Mathematica symbolic calculation software, we calculate \(c_3\),
and \(c_4\):
One can determine the SCGF itself by numerically solving (4.17) to find \(\beta _\mathrm{rest}(\lambda )\) and \(\theta (\lambda )\), then inserting these expressions into (4.18) and numerically integrating (3.6) to obtain \(F(\lambda )\). We show the results for \(d=2\) and \(d=3\) in Fig. 1. We verify that the resulting functions are convex, as they should by the general theory [86]. From these plots, many insights can be drawn out. In particular, it is possible to verify that the powerful fluctuation relations of Gallavotti–Cohen type [21, 103,104,105,106,107,108,109] hold. We leave an in-depth analysis to a future work.
Finally, we note that the phenomenon discussed in Sect. 3.5 can be explicitly seen here. Consider a thermal state boosted the sound velocity, \(\theta = \pm \theta _{\mathrm{s}}\). In this case, as a consequence of the discontinuities in \(\xi _{1,2}\) in (4.15), the derivatives of \(\beta _{\mathrm{rest}}(\lambda )\) and \(\theta (\lambda )\) with respect to \(\lambda \) have discontinuities at \(\lambda =0\). This implies that the third derivative of \(F(\lambda )\) does not exist at \(\lambda =0\). That is, the scaled cumulants \(c_n\) for \(n\ge 3\) do not exist, and the large-deviation principle is broken; these higher cumulants are expected to be divergent. Intuitively, when an object moves in a medium exactly at the speed of sound, there is a build up of linear waves generated. At the macroscopic scale, this appears to increase correlations of transported energy to such an extent so as to modify the scaling of higher-order cumulants with time. There thus appear to be a dynamical phase transition. On both sides of the phase transition point, the cumulants take different form, as is clear from (4.20) and (4.21).
5 Conclusion
In this paper, we show how to calculate the scaled cumulant generating function (or full counting statistics) for transport of any conserved quantity in stationary, homogeneous, clustering states of many-body systems, in or out of equilibrium. The technique is based on large-deviation theory, and the result is expressed in terms of quantities readily available from the Euler hydrodynamics description of the system. This can be seen as a nonlinear generalisation of the construction in [23, 24] for \(1+1\)-dimensional conformal field theory to interacting integrable [65] and non-integrable models. We show that the extended fluctuation relations proposed in [21] hold whenever the Euler hydrodynamics is linear. We extend the theory to arbitrary rays and make the connection with spatio-temporal correlation functions of twist fields, which have applications to order-parameter correlations. We also explain in what situations the theory may break and “dynamical phase transitions” may occur, making connection with nonlinear fluctuating hydrodynamics. Finally, we give the example of conformal hydrodynamics in arbitrary dimensions, obtaining the first exact results for energy transport cumulants in spatial dimensionality higher than 1. We observe, in this example, a breaking in thermal states boosted at the sound velocity.
Future works would include an in-depth study of fluctuations in higher-dimensional conformal hydrodynamics, especially in the non-equilibrium steady states constructed in [70, 71, 73, 74], as well as the analysis in other non-integrable models where ballistic transport exist, such as anharmonic chains or one-dimensional hard rods with alternating masses [110]. An understanding of how the Gallavotti–Cohen-type fluctuation relations emerge from our theory in the general setting of the Riemann problem of Euler hydrodynamics is also lacking. It would be interesting to extend the ideas developed here to include diffusion, and to make potential connections with macroscopic fluctuation theory. A more in-depth study of the exact formulae for dynamical correlation functions for order parameters and other twist field in stationary, homogeneous states would also be needed. A full connection with nonlinear fluctuating hydrodynamics, and fluctuations in the KPZ class, would be very interesting. It would be nice to see if the formalism can be appropriately extended in order to include the exact logarithmic large-deviation results found recently [20] (based on [111]).
Finally, as we already remarked, the extended fluctuation relations form a marker of “freeness”—they hold in free-particle models and \(1+1\)-dimensional conformal field theory. In these cases, Euler hydrodynamics is linear—equivalently, the natural flux Jacobian is state independent—and fluctuations in transport are directly related to initial-state fluctuations. This is in contrast to nonlinear Euler hydrodynamic systems, where the nonlinear evolution affects the structure of transport fluctuations as per the theory developed here. Hydrodynamic diffusion has also been argued to be a signal for the lack of interactions [54], which has been confirmed for integrable models [67, 68], and we observe that all known models where extended fluctuation relations hold also have vanishing hydrodynamic diffusion, and vice versa. Is there a relation between state independence of the natural flux Jacobian, and the vanishing of the diffusion matrix (although these two objects operate at different hydrodynamic scales)? Interestingly, this potential relation is further brought to light by a recent result [112], which can be interpreted as connecting, in the context of integrable systems, state differentiation of the flux Jacobian to the diagonal elements of the diffusion matrix. Why these two different hydrodynamic scales may be connected in this way remains to be explained.
Notes
In the stochastic case, the dynamics should be Markovian, and the various concepts used here have natural stochastic correspondents; for instance, conserved quantities should be interpreted as martingales.
Following standard notation, \(A(t)\asymp B(t)\) means \(\lim _{t\rightarrow \infty } (\log A(t))/(\log B(t))=1\).
These are in fact twist fields associated with continuous symmetries, where the conservation law arises from the associated Noether current. More generally, twist fields may be associated with discrete symmetries as well, in which case, however, there is no Noether current, hence no obvious flux Jacobian.
Results (3.31) and (3.32) hold as well for twist fields that are associated with discrete symmetries, as \(\Delta f(\lambda )\) is well defined in these cases also. Hence, the results are applicable to the study of the entanglement entropy in interacting models. They give an exact formula in integrable models using the thermodynamic Bethe ansatz [83]: in a state described by the source term \(w(\theta )\), for the n-copy permutation twist-field two-point function \(\langle {{\mathcal {T}}}_n(0,0)\bar{{\mathcal {T}}}_n(x,0)\rangle \), we have \(\Delta f_n = \int {\mathrm{d}}p(\theta )\,\log \left( \frac{1+{\mathrm{e}}^{-\epsilon _n(\theta )}}{(1+{\mathrm{e}}^{-\epsilon _1(\theta )})^n}\right) \), with pseudoenergy \(\epsilon _n(\theta )\) having source term \(nw(\theta )\) implementing the n-times larger imaginary time direction induced by the twist property, and \(p(\theta )\) being the momentum function.
In the isotropic Heisenberg chain, this was numerically observed in [102], after the first preprint of this paper appeared.
This invalidates the conjecture made in [70] for the SCGF, which was based on the extended fluctuation relations.
A more precise derivation can be obtained by rewriting the conservation laws in integral form, and assuming that, in the appropriate Euler scaling limit, (B.1) holds uniformly enough.
In this notation, we have in particular \({\mathsf {q}}_i^{~j} = \delta _i^j\) and \({\mathsf {j}}_i^{~j} = {\mathsf {A}}_i^{~j}\).
This solution is smooth except for contact discontinuities. It does not contain shocks, hence does not absorb (or generate) entropy, and is therefore expected to be the physically relevant solution.
References
Zwanzig, R.: Nonequilibrium Statistical Mechanics. Oxford University Press, New York (2001)
Kinoshita, T., Wenger, T., Weiss, D.S.: A quantum Newton’s cradle. Nature 440, 900 (2006)
Jezouin, S., Parmentier, F.D., Anthore, A., Gennser, U., Cavanna, A., Jin, Y., Pierre, F.: Quantum limit of heat flow across a single electronic channel. Science 342, 601 (2013)
Brantut, J.-P., Grenier, C., Meineke, J., Stadler, D., Krinner, S., Kollath, C., Esslinger, T., Georges, A.: A thermoelectric heat engine with ultracold atoms. Science 342, 713 (2013)
Eisert, J., Friesdorf, M., Gogolin, C.: Quantum many-body systems out of equilibrium. Nat. Phys. 11, 124 (2015)
Touchette, H.: The large deviation approach to statistical mechanics. Phys. Rep. 478, 1 (2009)
Touchette, H., Harris, R.J.: Large deviation approach to nonequilibrium systems. In: Klages, R., Just, W., Jarzynski, C. (eds.) Nonequilibrium Statistical Physics of Small Systems: Fluctuation Relations and Beyond. Wiley, Weinheim (2013)
Marconi, U.M.B., Puglisi, A., Rondoni, L., Vulpiani, A.: Fluctuation-dissipation: response theory in statistical physics. Phys. Rep. 461, 111 (2008)
Esposito, M., Harbola, U., Mukamel, S.: Nonequilibrium fluctuations, fluctuation theorems, and counting statistics in quantum systems. Rev. Mod. Phys. 81, 1665 (2009)
Derrida, B.: Non-equilibrium steady states: fluctuations and large deviations of the density and of the current. J. Stat. Mech. Theory Exp. 2007, P07023 (2007)
Derrida, B.: Microscopic versus macroscopic approaches to non-equilibrium systems. J. Stat. Mech. Theory Exp. 2011, P01030 (2011)
Lazarescu, A.: The physicist’s companion to current fluctuations: one-dimensional bulk-driven lattice gases. J. Phys. A 48, 503001 (2015)
Levitov, S., Lesovik, G.: Charge distribution in quantum shot noise. JETP Lett. 58, 230 (1993)
Avron, J.E., Bachmann, S., Graf, G.M., Klich, I.: Fredholm determinants and the statistics of charge transport. Commun. Math. Phys 280, 807 (2008)
Bernard, D., Doyon, B.: Full counting statistics in the resonant-level model. J. Math. Phys. 53, 122302 (2012)
Saito, K., Dhar, A.: Fluctuation theorem in quantum heat conduction. Phys. Rev. Lett. 99, 180601 (2007)
Doyon, B., Lucas, A., Schalm, K., Bhaseen, M.J.: Non-equilibrium steady states in the Klein–Gordon theory. J. Phys. A 48, 095002 (2015)
Gawȩdzki, K., Tauber, C.: Nonequilibrium transport through quantum-wire junctions and boundary defects for free massless bosonic fields. Nucl. Phys. B 896, 138 (2015)
Yoshimura, T.: Full counting statistics in the free Dirac theory. J. Phys. A 51, 475002 (2018)
Moriya, H., Nagao, R., Sasamoto, T.: Exact large deviation function of spin current for the one dimensional XX spin chain with domain wall initial condition. J. Stat. Mech. 2019, 063105 (2019)
Bernard, D., Doyon, B.: Time-reversal symmetry and fluctuation relations in non-equilibrium quantum steady states. J. Phys. A 46, 372001 (2013)
Komnik, A., Saleur, H.: Quantum fluctuation theorem in an interacting setup: point contacts in fractional quantum Hall edge state devices. Phys. Rev. Lett. 107, 100601 (2011)
Bernard, D., Doyon, B.: Energy flow in non-equilibrium conformal field theory. J. Phys. A 45, 362001 (2012)
Bernard, D., Doyon, B.: Non-equilibrium steady states in conformal field theory. Ann. Henri Poincaré 16, 113 (2015)
Bernard, D., Doyon, B.: Conformal field theory out of equilibrium: a review. J. Stat. Mech. Theor. Exp. 2016, 064005 (2016)
Spohn, H., Lebowitz, J.L.: Stationary non-equilibrium states of infinite harmonic systems. Commun. Math. Phys. 54, 97 (1977)
Ruelle, D.: Natural nonequilibrium states in quantum statistical mechanics. J. Stat. Phys. 98, 57 (2000)
Jakšić, V., Pillet, C.-A.: Mathematical theory of non-equilibrium quantum statistical mechanics. J. Stat. Phys. 108, 787 (2002)
Derrida, B., Lebowitz, J.L., Speer, E.R.: Free energy functional for nonequilibrium systems: an exactly solvable case. Phys. Rev. Lett. 87, 150601 (2001)
Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., Landim, C.: Fluctuations in stationary nonequilibrium states of irreversible processes. Phys. Rev. Lett. 87, 040601 (2001)
Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., Landim, C.: Macroscopic fluctuation theory for stationary non-equilibrium states. J. Stat. Phys. 107, 635 (2002)
Bodineau, T., Derrida, B.: Current fluctuations in nonequilibrium diffusive systems: an additivity principle. Phys. Rev. Lett. 92, 180601 (2004)
Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., Landim, C.: Macroscopic fluctuation theory. Rev. Mod. Phys. 87, 593 (2015)
Buča, B., Prosen, T.: Exactly solvable counting statistics in open weakly coupled interacting spin systems. Phys. Rev. Lett. 112, 06720 (2014)
Žnidarič, M.: Large-deviation statistics of a diffusive quantum spin chain and the additivity principle. Phys. Rev. E 89, 042140 (2014)
Žnidarič, M.: Exact large-deviation statistics for a nonequilibrium quantum spin chain. Phys. Rev. Lett. 112, 040602 (2014)
Carollo, F., Garrahan, J.P., Lesanovsky, I., Pérez-Espigares, C.: Fluctuating hydrodynamics, current fluctuations, and hyperuniformity in boundary-driven open quantum chains. Phys. Rev. E 96, 052118 (2017)
Bauer, M., Bernard, D., Jin, T.: Stochastic dissipative quantum spin chains (I): quantum fluctuating discrete hydrodynamics. SciPost Phys. 3, 033 (2017)
Bauer, M., Bernard, D., Jin, T.: Equilibrium fluctuations in maximally noisy extended quantum systems. SciPost Phys. 6, 045 (2019)
Buča, B., Garrahan, J. P., Prosen, T., Vanicat, M.: Exact large deviation statistics and trajectory phase transition of a deterministic boundary driven cellular automaton (2019). Preprint arXiv:1901.00845
Rigol, M., Dunjko, V., Yurovsky, V., Olshanii, M.: Relaxation in a completely integrable many-body quantum system: an ab initio study of the dynamics of the highly excited states of 1D lattice hard-core bosons. Phys. Rev. Lett. 98, 050405 (2007)
Ilievski, E., De Nardis, J., Wouters, B., Caux, J.-S., Essler, F.H.L., Prosen, T.: Complete generalized Gibbs ensembles in an interacting theory. Phys. Rev. Lett. 115, 157201 (2015)
Doyon, B.: Thermalization and pseudolocality in extended quantum systems. Commun. Math. Phys. 351, 155 (2017)
Essler, F.H.L., Fagotti, M.: Quench dynamics and relaxation in isolated integrable quantum spin chains. J. Stat. Mech. Theory Exp. 2016, 064002 (2016)
Langen, T., Erne, S., Geiger, R., Rauer, B., Schweigler, T., Kuhnert, M., Rohringer, W., Mazets, I.E., Gasenzer, T., Schmiedmayer, J.: Experimental observation of a generalized Gibbs ensemble. Science 348, 207 (2015)
Castro-Alvaredo, O.A., Doyon, B., Yoshimura, T.: Emergent hydrodynamics in integrable quantum systems out of equilibrium. Phys. Rev. X 6, 041065 (2016)
Bertini, B., Collura, M., De Nardis, J., Fagotti, M.: Transport in out-of-equilibrium \(XXZ\) chains: exact profiles of charges and currents. Phys. Rev. Lett. 117, 207201 (2016)
Spohn, H.: Large Scale Dynamics of Interacting Particles. Springer, Heidelberg (1991)
Spohn, H.: Nonlinear fluctuating hydrodynamics for anharmonic chains. J. Stat. Phys. 155, 1191 (2014)
Mendl, C.B., Spohn, H.: Current fluctuations for anharmonic chains in thermal equilibrium. J. Stat. Mech. Theor. Exp. 2015, P03007 (2015)
Doyon, B., Spohn, H.: Drude weight for the Lieb-Liniger Bose gas. SciPost Phys. 3, 039 (2017)
Chetrite, R., Touchette, H.: Nonequilibrium Markov processes conditioned on large deviations. Ann. Henri Poincaré 16, 2005 (2015)
Carollo, F., Garrahan, J.P., Lesanovsky, I., Pérez-Espigares, C.: Making rare events typical in Markovian open quantum systems. Phys. Rev. A 98, 010103 (2018)
Spohn, H.: Interacting and noninteracting integrable systems. J. Math. Phys. 59, 091402 (2018)
Bodineau, T., Derrida, B.: Distribution of current in nonequilibrium diffusive systems and phase transitions. Phys. Rev. E 72, 066110 (2005)
Garrahan, J.P., Jack, R.L., Lecomte, V., Pitard, E., vanDuijvendijk, K., van Wijland, F.: Dynamical first-order phase transition in kinetically constrained models of glasses. Phys. Rev. Lett. 98, 195702 (2007)
Garrahan, J.P., Jack, R.L., Lecomte, V., Pitard, E., vanDuijvendijk, K., van Wijland, F.: First-order dynamical phase transition in models of glasses: an approach based on ensembles of histories. J. Phys. A Math. Theor. 42, 075007 (2009)
Espigares, C.P., Garrido, P.L., Hurtado, P.I.: Dynamical phase transition for current statistics in a simple driven diffusive system. Phys. Rev. E 87, 032115 (2013)
Jack, R.L., Sollich, P.: Large deviations of the dynamical activity in the East model: analysing structure in biased trajectories. J. Phys. A: Math. Theor. 47, 015003 (2014)
Hurtado, P.I., Espigares, C.P., del Pozo, J.J., Garrido, P.L.: Thermodynamics of currents in nonequilibrium diffusive systems: theory and simulation. J. Stat. Phys. 154, 214 (2014)
Tsobgni Nyawo, P., Touchette, H.: A minimal model of dynamical phase transition. Eur. Phys. Lett. 116, 50009 (2016)
Tizón-Escamilla, N., Pérez-Espigares, C., Garrido, P.L., Hurtado, P.I.: Order and symmetry-breaking in the fluctuations of driven systems. Phys. Rev. Lett. 119, 090602 (2017)
Lazarescu, A.: Generic dynamical phase transition in one-dimensional bulk-driven lattice gases with exclusion. J. Phys. A Math. Theor. 50, 254004 (2017)
Chen, Z., de Gier, J., Hiki, I., Sasamoto, T.: Exact confirmation of 1D nonlinear fluctuating hydrodynamics for a two-species exclusion process. Phys. Rev. Lett. 120, 240601 (2018)
Myers, J., Bhaseen, M. J., Harris, R. J., Doyon, B.: Transport fluctuations in integrable models out of equilibrium (2018). Preprint arXiv:1812.02082
Doyon, B., Yoshimura, T.: A note on generalized hydrodynamics: inhomogeneous fields and other concepts. SciPost Phys. 2, 014 (2017)
De Nardis, J., Bernard, D., Doyon, B.: Hydrodynamic diffusion in integrable systems. Phys. Rev. Lett. 121, 160603 (2018)
De Nardis, J., Bernard, D., Doyon, B.: Diffusion in generalized hydrodynamics and quasiparticle scattering. SciPost Phys. 6, 049 (2019)
Schemmer, M., Bouchoule, I., Doyon, B., Dubail, J.: Generalized hydrodynamics on an atom chip. Phys. Rev. Lett. 122, 090601 (2019)
Bhaseen, M.J., Doyon, B., Lucas, A., Schalm, K.: Far from equilibrium energy flow in quantum critical systems. Nat. Phys. 11, 509 (2015)
Chang, H.-C., Karch, A., Yarom, A.: An ansatz for one dimensional steady state configurations. J. Stat. Mech. Theory Exp. 2014, P06018 (2014)
Pourhasan, R.: Non-equilibrium steady state in the hydro regime. J. High Energy Phys. 02, 005 (2016)
Lucas, A., Schalm, K., Doyon, B., Bhaseen, M.J.: Shock waves, rarefaction waves and non-equilibrium steady states in quantum critical systems. Phys. Rev. D 94, 025004 (2016)
Spillane, M., Herzog, C.P.: Relativistic hydrodynamics and non-equilibrium steady states. J. Stat. Mech. Theory Exp. 2016, 103208 (2016)
Ilievski, E., Medenjak, M., Prosen, T., Zadnik, L.: Quasilocal charges in integrable lattice systems. J. Stat. Mech. 2016, 064008 (2016)
Israel, R.B.: Convexity in the Theory of Lattice Gases. Princeton University Press, Princeton (1979)
Bratteli, O., Robinson, D.W.: Operator Algebras and Quantum Statistical Mechanics 1. Springer, Berlin (1987)
Bratteli, O., Robinson, D.W.: Operator Algebras and Quantum Statistical Mechanics 2. Springer, Berlin (1987)
Araki, H.: Gibbs states of a one dimensional quantum lattice. Commun. Math. Phys. 14, 120 (1969)
Prosen, T.: Quantum invariants of motion in a generic many-body system. J. Phys. A. 31, L645 (1998)
Prosen, T.: Ergodic properties of a generic nonintegrable quantum many-body system in the thermodynamic limit. Phys. Rev. E 60, 3949 (1999)
Yang, C.N., Yang, C.P.: Thermodynamics of a one-dimensional system of bosons with repulsive delta-function interaction. J. Math. Phys. 10, 1115 (1969)
Zamolodchikov, A.: Thermodynamic Bethe ansatz in relativistic models: scaling 3-state potts and Lee-Yang models. Nucl. Phys. B 342, 695 (1990)
Takahashi, M.: Thermodynamics of One-Dimensional Solvable Models. Cambridge University Press, Cambridge (1999)
Ilievski, E., Quinn, E., Caux, J.-S.: From interacting particles to equilibrium statistical ensembles. Phys. Rev. B 95, 115128 (2017)
Bressan, A.: Hyperbolic conservation laws: an illustrated tutorial. In: Modelling and Optimisation of Flows on Networks. Lecture Notes in Mathematics, vol. 2062, p. 157. Springer, Berlin (2013)
Doyon, B.: Exact large-scale correlations in integrable systems out of equilibrium. SciPost Phys. 5, 054 (2018)
Tsvelik, A.: Quantum Field Theory in Condensed Matter Physics. Cambridge University Press, Cambridge (2003)
Sachdev, S.: Quantum Phase Transitions. Cambridge University Press, Cambridge (2011)
Cardy, J.L., Castro Alvaredo, O.A., Doyon, B.: Form factors of branch-point twist fields in quantum integrable models and entanglement entropy. J. Stat. Phys. 130, 129 (2008)
Castro-Alvaredo, O.A., Doyon, B., Fioravanti, D.: Conical twist fields and null polygonal Wilson loops. Nucl. Phys. B 931, 146 (2018)
Chen, Y., Doyon, B.: Form factors in equilibrium and non-equilibrium mixed states of the ising model. J. Stat. Mech. 2014, P09021 (2014)
Calabrese, P., Essler, F.H.L., Fagotti, M.: Quantum quench in the transverse field Ising chain II: stationary state properties. J. Stat. Mech. 2012, P07022 (2012)
Alba, V., Calabrese, P.: Entanglement and thermodynamics after a quantum quench in integrable systems. PNAS 114, 7947 (2017)
Alba, V., Calabrese, P.: Entanglement dynamics after quantum quenches in generic integrable systems. SciPost Phys. 4, 017 (2018)
Alba, V., Calabrese, P.: Quench action and Rényi entropies in integrable systems. Phys. Rev. B 96, 115421 (2017)
Alba, V., Calabrese, P.: Rényi entropies after releasing the Néel state in the XXZ spin-chain. J. Stat. Mech. 2017, 113105 (2017)
Mestyán, M., Alba, V., Calabrese, P.: Rényi entropies of generic thermodynamic macrostates in integrable systems. J. Stat. Mech. 2018, 083104 (2018)
Bastianello, A., Calabrese, P.: Spreading of entanglement and correlations after a quench with intertwined quasiparticles. SciPost Phys. 5, 033 (2018)
Bertini, B., Fagotti, M., Piroli, L., Calabrese, P.: Entanglement evolution and generalised hydrodynamics: noninteracting systems. J. Phys. A 51, 39LT01 (2018)
Piroli, L., De Nardis, J., Collura, M., Bertini, B., Fagotti, M.: Transport in out-of-equilibrium XXZ chains: non-ballistic behavior and correlation functions. Phys. Rev. B 96, 115124 (2017)
Ljubotina, M., Znidaric, M., Prosen, T.: Kardar–Parisi–Zhang physics in the quantum Heisenberg magnet (2019). arXiv:1903.01329
Gallavotti, G., Cohen, E.G.D.: Dynamical ensembles in nonequilibrium statistical mechanics. Phys. Rev. Lett. 74, 2694 (1995)
Jarzynski, C.: Nonequilibrium equality for free energy differences. Phys. Rev. Lett. 78, 2690 (1997)
Jarzynski, C.: Equilibrium free-energy differences from non-equilibrium measurements: a master-equation approach. Phys. Rev. E 56, 5018 (1997)
Crooks, G.: Nonequilibrium measurements of free energy differences for microscopically reversible Markovian systems. J. Stat. Phys. 90, 1481 (1998)
Kurchan, J.: Fluctuation theorem for stochastic dynamics. J. Phys. A 31, 3719 (1998)
Lebowitz, J.L., Spohn, H.: A Gallavotti–Cohen-type symmetry in the large deviation functional for stochastic dynamics. J. Stat. Phys. 95, 333 (1999)
Jarzynski, C., Wójcik, D.K.: Classical and quantum fluctuation theorems for heat exchange. Phys. Rev. Lett. 92, 230602 (2004)
Brunet, E., Derrida, B., Gerschenfeld, A.: Fluctuations of the heat flux of a one-dimensional hard particle gas. Eur. Phys. Lett. 90, 20004 (2010)
Eisler, V., Rácz, Z.: Full counting statistics in a propagating quantum front and random matrix spectra. Phys. Rev. Lett. 110, 060602 (2013)
Gopalakrishnan, S., Huse, D.A., Khemani, V., Vasseur, R.: Hydrodynamics of operator spreading and quasiparticle diffusion in interacting integrable systems. Phys. Rev. B 98, 220303 (2018)
Ruelle, D.: Statistical Mechanics. Benjamin, New-York - Amsterdam (1969)
Bastianello, A., Doyon, B., Watts, G., Yoshimura, T.: Generalized hydrodynamics of classical integrable field theory: the sinh-Gordon model. SciPost Phys. 4, 45 (2018)
Lieb, E.H., Robinson, D.W.: The finite group velocity of quantum spin systems. Commun. Math. Phys. 28, 251 (1972)
Bravyi, S., Hastings, M.B., Verstraete, F.: Lieb-Robinson bounds and the generation of correlations and the topological quantum order. Phys. Rev. Lett. 97, 050401 (2006)
Acknowledgements
We thank M. Joe Bhaseen, Denis Bernard, Juan P. Garrahan, Krzysztof Gawedzki, Rosemary J. Harris, Tomohiro Sasamoto, Herbert Spohn, Pierpaolo Vivo, Takato Yoshimura and Claudio Zeni, and concerning Sect. 3.4, Vincenzo Alba, Pasquale Calabree and Martón Kormos, for discussions, comments and suggestions. BD acknowledges funding from the Royal Society under a Leverhulme Trust Senior Research Fellowship, “Emergent hydrodynamics in integrable systems: non-equilibrium theory”, ref. SRF\R1\180103, and from the EPSRC through a standard grant, “Entanglement Measures, Twist Fields, and Partition Functions in Quantum Field Theory”, ref. EP/P006132/1. JM acknowledges funding from the EPSRC Centre for Doctoral Training in Cross-Disciplinary Approaches to Non-Equilibrium Systems (CANES) under Grant EP/L015854/1. BD is grateful to École Normale Supérieure de Paris for an invited professorship (2018) and acknowledges hospitality and funding from the Erwin Schrödinger Institute in Vienna (Austria) during the programme “Quantum Paths” (2018), and the Galileo Galilei Institute in Florence (Italy) during the programme “Entanglement in Quantum Systems” (2018). BD’s research was supported in part by Perimeter Institute for Theoretical Physics. Research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development and by the Province of Ontario through the Ministry of Research, Innovation and Science. Both authors thank the Centre for Non-Equilibrium Science (CNES) and the Thomas Young Centre (TYC).
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Christian Maes.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendices
Derivation of the Main Result (3.6) with (3.2)
The derivation is based mainly on the assumptions of strong enough clustering of correlation functions of local observables, both in space and in time, along with a standard result from hydrodynamics, Eq. (A.5), which can be seen as a weak version of the hydrodynamic projection principle. Clustering is to be strong enough. For instance, exponential clustering at large spacial separations can be shown rigorously in extremal KMS states associated with local Hamiltonians (see, e.g., [77,78,79, 113]). Strong enough clustering in time is more difficult to prove, but expected to hold generically in many systems and states—as explained in Sect. 3.5, it is broken in some situations, leading to a breaking of the large-deviation principle (2.9).
For completeness, we present various ways of proving the results, which involve different assumptions and principles. In one way of proving the emergence of a flow on the MES manifold, we make use of a theorem from [43] which shows that the pseudolocal charges, a concept originally introduced in [80, 81] (see [75]), form a Hilbert space that describe the tangent to the MES manifold as per (2.4). We require not only strong clustering in time, but also strong clustering along all rays \(\xi =x/t\) in a neighbourhood of the \(\xi =0\), expected to hold if no effective velocity lies within such a neighbourhood. In [43], the particular context of quantum statistical mechanics, specifically the quasi-local \(C^*\) algebras, is taken, and a specific definition of MES (which are in [43] referred to as GGEs) is used. This definition is in agreement with the one used here, in particular with (2.4), if we assume the tangent spaces, which may be different at different points along the flow, to have the same countable basis all along the flow, the \(Q_i\)’s. In generic, non-integrable system, we would in fact expect all tangent spaces to be finite-dimensional, but this is difficult to prove.
In another way of proving the emergence of a flow on the MES manifold, we instead make use of a stronger version of the hydrodynamic projection principle (Eq. (A.8)). Hydrodynamic projection principles are nontrivial, but have been used successfully in the context of statistical fluid dynamics [48,49,50,51], see also [87] where an even stronger version of (A.8) is explicitly used to derive correlation results verified numerically in [114].
In addition, we make various technical assumptions, such as appropriate boundedness and differentiability assumptions for correlation functions. The precise specification of the assumptions of clustering, boundedness and differentiability is possible but would require a full mathematical framework—such as that of \(C^*\) algebras; we hope to come back to such matters in future works.
1.1 Statement of the Problem
Consider expression (2.10) for the scaled cumulants \(c_n\). Let us assume that correlation functions of the current \(\mathcal {j}(0,t)\) cluster strongly enough at large time differences. Using stationarity of the state, standard arguments show that the scaled cumulants exist and can be written in the form
where the limits are taken, in order, on \(t_{n-1},t_{n-2},\ldots ,t_{1}\), as written. Here the many-point connected, symmetrised correlation function on the right-hand side is a natural generalisation of (2.6): we define \(\langle {{\mathcal {O}}}_1,\ldots ,{{\mathcal {O}}}_n\rangle ^{(\mathrm c)}_{{\underline{\beta }}}\) as the connected part of the expectation value of the normalised iterated anti-commutators
(where \(\{{{\mathcal {O}}},{{\mathcal {O}}}'\} = {{\mathcal {O}}}{{\mathcal {O}}}' + {{\mathcal {O}}}'{{\mathcal {O}}}\) is the anti-commutator).
Let us assume more generally that correlation functions involving the currents \(\mathcal {j}(0,t)\) and local observables cluster strongly enough at large time differences. Then standard arguments show that (i) the state \(\langle \cdots \rangle ^{(\lambda )}\) defined by the series expansion in \(\lambda \) of (3.9) has a nonzero radius of convergence for any local observable \({{\mathcal {O}}}(x,t)\) (the radius may depend on the observable), (ii) for any local observable \({{\mathcal {O}}}(x,t)\), (3.10) holds within the convergence region, (iii) the series expansion in (2.10) for \(F(\lambda )\) has a nonzero radius of convergence, and (iv) result (3.12) holds within the convergence region.
There are then two ways—which are equivalent under an appropriate assumption of clustering at large time differences—to define precisely the function \(F(\lambda )\), in order to prove (3.6) with (3.2):
- I.
We may define \(F(\lambda )\) as the solution to (3.12) with \(F(0)=0\), where the state \(\langle \cdots \rangle ^{(\lambda )}\) is defined as the solution to (3.10) with \(\langle \cdots \rangle ^{(0)} = \langle \cdots \rangle _{{\underline{\beta }}}\).
- II.
We may take the explicit form of cumulants (A.1), along with the second equation in (2.10), as our working expression.
Below, we give simple proofs (3.6) with (3.2) under both definitions. The arguments presented in Sect. 3.2 follow the proof that is natural under definition I; this proof is conceptual, but neglects technical difficulties about the manifold structure of MES. The proof under definition II is more explicit.
1.2 Hydrodynamic Projections
The main assumption for this section is the hydrodynamic projection principle. This principle has various implementations. We first use a very weak version of it, which expresses Fourier transforms of dynamical two-point functions of conserved densities in terms of the static correlation matrix \({\mathsf {C}}\) and the flux Jacobian \({\mathsf {A}}\), in the long-wavelength, longtime limit. This is in fact a standard result in linear fluctuating hydrodynamics.
Define
as well as the “Euler scaling limit” of long wavelengths and long times,
The weak version of the hydrodynamics projection principle that we assume is that the time dependence of \(S_{ij}(k,t)\) takes the form
We need one additional assumption. As we will present two alternative arguments, this additional assumption takes two different forms, depending on the argument used.
In one form, we require that there exist a \(k\in (0,1)\) such that clustering of correlation functions of local observables in time is strong enough along all rays \(\xi =x/t\) in the interval \([-k,k]\),
We believe exponential clustering with uniform exponent would be sufficient, but it is not necessary.
In the other form, a stronger version of the hydrodynamic projection principle is required to hold: Fourier transforms of two-point functions involving a conserved density and an arbitrary local observable are expressed as Fourier transforms of conserved density two-point functions. Define, for local observables \({{\mathcal {O}}}(x,t)\),
and the Euler scaling limit \(S_{i{{\mathcal {O}}}}^{\mathrm{Eul}}(kt)\) as in (A.4). Then, the stronger version is
where the inner product is defined in (2.5) and \({\mathsf {C}}\) is assumed to be invertible. Equation (A.8) represents the idea that correlations, at large scales, are produced by propagation of ballistically transported quantities, and thus, two-point functions of the local observables are fully determined by evaluating their overlap with the conserved charges and propagating conserved densities. The static correlation matrix provides the metric in the space of conserved densities to use for the completeness relation. Below, we provide a proof of the weaker version (A.5) as a consequence of assumption (A.8), for completeness. We note that assumption (A.8), with (A.5), immediately implies
With these assumptions, the main result of this subsection is to show that
This implies in particular (3.16), but is more general.
Proof of (A.10)
Note that (A.9), specialised to \({{\mathcal {O}}}(x,t) = \mathcal {q}_j(x,t)\), is expression (A.5), our main assumption (weak hydrodynamic projection principle). Below, we show (A.10) from (A.9). Under the strong hydrodynamic projection principle (A.8), this completes the proof. Under the weaker one (A.5), we need one additional step: we need to show that there exists a conserved density \(\tilde{\mathcal {q}}_i(x,t) = \sum _j {\mathsf {V}}_{i}^{~j}\mathcal {q}_i(x,t)\) such that
for all local observables \({{\mathcal {O}}}(x,t)\). Once this is shown, specialising it to \({{\mathcal {O}}}(x,t)=\mathcal {q}_j(x,t)\) and using (A.10) for \({{\mathcal {O}}}(x,t) = \mathcal {q}_j(x,t)\), we determine that \({\mathsf {V}}_{i}^{~j} = {{\,\mathrm{sgn}\,}}({\mathsf {A}})_{i}^{~j}\), which completes the proof.
We first prove (A.11). The main idea is to show that \(\int _{-\infty }^\infty {\mathrm{d}}t\,\mathcal {j}_i(0,t)\) is a conserved pseudolocal charge. This is useful, as [43, Defs 5.4, 6.2] says that the tangent space to a MES (as per (2.4)) is the conserved subspace of the Hilbert space completion of the inner product induced by (2.6), and [43, Thm 5.7] shows that this Hilbert space completion is in bijection with the space of pseudolocal charges. The precise definition of pseudolocal charges is given in [43, Def 5.1] (we use the “two-sided pseudolocal charges” (see [43, Sect 5.2])). Making contact with this definition, we must consider \(t\mapsto J^{(t)} = \int _{-t/2}^{t/2} {\mathrm{d}}s\,\mathcal {j}_i(0,t)\) (with an unimportant shift with respect to (2.8)) to form an infinite sequence of operators with increasing t, supported on increasing intervals whose length is proportional to t. The definition [43, Def 5.1] asks for the support to be strict, while, under time evolution, it is known by the Lieb–Robinson theorem [115, 116] that the support is only exponentially accurate. But by exponential accuracy, it is possible, using the techniques of [116], to approximate the time-evolved fields by observables with strict supports, and modify the sequence in order for the strict support to grow linearly and converge to the same object. Here for simplicity we assume this has been done and that \(J^{(t)}\) has strict support growing linearly with t. We must then check three requirements for the large-t limit of \(J^{(t)}\) to form, in the sense of [43, Def 5.1], a pseudolocal charge: (I) that \(\langle J^{(t)},J^{(t)}\rangle ^{\mathrm{c}}_{{\underline{\beta }}}\) has a growth that is bounded linearly in t; (II) that \(\lim _{t\rightarrow \infty }\langle J^{(t)},{{\mathcal {O}}}(x,0)\rangle ^{\mathrm{c}}\) exists for all local observables \({{\mathcal {O}}}(x,0)\); and (III) that the result of the latter is independent of x. For the latter point, in fact, a strong enough independence must hold not just in the limit: that there exists \(k\in (0,1)\) such that the supremum of the difference \(|\langle J^{(t)},{{\mathcal {O}}}(x,0)\rangle ^{\mathrm{c}}-\langle J^{(t)},{{\mathcal {O}}}(y,0)\rangle ^{\mathrm{c}}|\) within the region \(x,y\in [-kt,kt]\) tends to zero as \(t\rightarrow \infty \). As in the large-t limit it is clear that \(J^{(t)}\) is conserved, with these three requirements, (A.11) follows from [43, Thm 5.7].
The first two requirements are immediate from strong enough clustering in time. The last one follows from the conservation laws (2.1), and strong enough clustering along all rays \(\xi =x/t\) with \(|\xi |\in [0,k]\), as expressed in (A.6). Indeed, the conservation laws allow us bound the supremum of the difference by
If clustering is strong enough along all rays \(\xi =x/t\in [-k,k]\) (and using the dominated convergence theorem), the large-t limit vanishes. This completes the proof of (A.11).
We now prove (A.10) from (A.9). Let us take some \(t>0\), and consider a function f(x) whose derivative is \(f'(x) \propto e^{-\mu x^2}\) for some \(\mu >0\), with the normalisation condition
Then using (A.9) as well as the dominated convergence theorem, we have
The equation also holds with \(t\mapsto -t\). By the conservation laws,
By clustering in space, we can use integration by parts, and we find
Again by the conservation laws,
Inserting this in the right-hand side of (A.16) and using (A.14), we obtain
Using the first equation in (A.13) to evaluate the first term on the right-hand side, and simplifying the integral in the second term, we find
where, using the second equation in (A.13),
The second term on the right-hand side of (A.19) can be evaluated by (A.14). Simplifying by using (A.20) and recalling that \(t>0\), we obtain (A.10). \(\square \)
Proof of (A.5) from (A.8)
First note that by symmetry (in particular, recall definition (2.6)), assumption (A.8) implies
where \(S_{{{\mathcal {O}}}i}^{\mathrm{Eul}}(kt)\) is the Euler scaling limit of \(S_{{{\mathcal {O}}}i}(k,t) = \int _{-\infty }^\infty {\mathrm{d}}x\,{\mathrm{e}}^{-{\mathrm{i}}kx}\langle {{\mathcal {O}}}(x,t)\mathcal {q}_i(0,0)\rangle ^{\mathrm{c}}_{{\underline{\beta }}}\). Using the conservation laws (2.1) and integration by parts, we obtain
From (A.21), we then find in the Euler scaling limit, assuming that we can exchange the derivative and the limit,
whose solution, with the initial condition \(S_{ij}^{\mathrm{Eul}}(0) = {\mathsf {C}}_{ij}\), is (A.5). \(\square \)
1.3 \(\lambda \)-Flow
We define the “Lie derivative” \({{\mathcal {L}}} \langle \cdots \rangle _{{\underline{\beta }}}\), at the point \(\langle \cdots \rangle _{{\underline{\beta }}}\) in the MES manifold characterised by the Lagrange parameters \({\underline{\beta }}\), by
Leibniz’s rule fixes its form on products of expectations, and in particular on connected, symmetrised correlation functions, a combinatoric analysis from the definition given around (A.2) gives
As a consequence of (A.10),
and, by Leibniz’s rule,
Let us show (3.6) with (3.2) using definition I of “Appendix A.1”. First (A.26) shows that the Lie derivative on the MES manifold lies within the tangent space of the manifold. Therefore, by (3.10), the flow \(\lambda \mapsto \langle \cdots \rangle ^{(\lambda )}\) is that along the direction set by the Lie derivative \({{\mathcal {L}}}\langle \cdots \rangle _{{\underline{\beta }}}\), starting at some point \({\underline{\beta }}(0)\), which lies entirely on the MES manifold. Hence, it can be characterised by \(\lambda \mapsto {\underline{\beta }}(\lambda )\), which solves (3.2):
Finally, the differential equation (3.12) with \(F(0)=0\) shows (3.6). \(\square \)
On the other hand, let us show (3.6) with (3.2) using definition II of “Appendix A.1”, from the explicit expressions of the cumulants \(c_n\) given in (A.1). It is sufficient to show that
where \({\underline{\beta }}(\lambda )\) solves (3.2). We show this by induction. The induction is on the statement that
In particular, this implies (A.29). If (A.30) holds for \(n=m\), then we have
Assuming that we can exchange the limits and integrals with the \(\lambda \)-derivative, this gives
where on the third line we used (3.2) (along with Leibniz’s rule), on the fourth (A.27), and on the fifth (A.25). This shows (A.30) for \(n=m+1\). Since (A.30) holds by definition for \(n=1\), this completes the proof. \(\square \)
1.4 Second Cumulant
The second cumulant was shown in [50, 51] to satisfy the sum rule
Using the Euler-scale expression (A.5), this gives
where \(|A|={{\,\mathrm{sgn}\,}}(A)A\). On the other hand, from our result (3.6) with (3.2), we have
where in the last step we used (B.5). This agrees with (A.34).
Euler Hydrodynamics
1.1 Standard Results
Consider a non-stationary, inhomogeneous state \(\langle \cdots \rangle \) of the system. In Euler hydrodynamics, one assumes the every local average at space-time point x, t can be approximated by a local MES, which depends on x, t but not on the observable whose average is taken,
(since the MES is homogeneous and stationary, one can put the observable at 0, 0 on the right-hand side). Physically, there are fluid cells, which are large compared to microscopic scales but small compared to variation scales of the states, in which the state has, to a good approximation, maximised entropy, and is very nearly homogeneous and stationary. Writing the conservation laws (2.1) in average form within the state \(\langle \cdots \rangle \) and using the approximation (B.1), one obtainsFootnote 10
where \({\mathsf {q}}_i(x,t) = \langle \mathcal {q}_i(0,0)\rangle _{{\underline{\beta }}(x,t)}\) and \({\mathsf {j}}_i(x,t) = \langle \mathcal {q}_i(0,0)\rangle _{{\underline{\beta }}(x,t)}\). This can be rewritten using the flux Jacobian (3.1),
Since \({\mathsf {q}}_i(x,t)\) form a system of coordinates for the MES at x, t, these are equations of motion of the space-time-dependent MES—they are the Euler hydrodynamic equations corresponding to the dynamical system of interest.
It is known [48,49,50,51] (see also the explicit proof in [46]) that there exists a generating function g for the currents,
The function g is a generating function for the average currents in MES, much like the specific free energy is for the average densities. Equation (B.4) is a consequence of the fact that \(-\partial {\mathsf {j}}_i/\partial \beta ^j = \int {\mathrm{d}}x\,\langle \mathcal {q}_j(x,0),\mathcal {j}_i(0,0)\rangle _{{\underline{\beta }}}^{\mathrm{c}}\) is symmetric under \(i\leftrightarrow j\), which can be shown by using the conservation laws and integration by parts. Note that changing variables to \(\underline{{\mathsf {q}}}\), this symmetry also immediately implies the important relation
involving the static correlation matrix \({\mathsf {C}} = \frac{\partial \underline{{\mathsf {q}}}}{\partial {\underline{\beta }}}\) (see (3.5)). In particular, this implies that with respect to the inner product defined by the \({\mathsf {C}}\) matrix, \(( {\mathsf {a}},{\mathsf {b}}) = \sum _{ij} {\mathsf {a}}^{i}{\mathsf {b}}^j{\mathsf {C}}_{ij}\), which is non-degenerate if \({\mathsf {C}}\) is invertible (\({\mathsf {C}}\) is always nonnegative), the matrix \({\mathsf {A}}\) is Hermitian.Footnote 11 Hence, if \({\mathsf {C}}\) is invertible, then \({\mathsf {A}}\) is diagonalisable by a similarity transformation.
The normal coordinates of Euler hydrodynamics are a different system of coordinates, \(\underline{{\mathsf {q}}}\mapsto {\underline{n}}\), which diagonalise the flux Jacobian. That is, there is a diagonal matrix \(v^{\mathrm{eff}}\), with diagonal elements \(v^{\mathrm{eff}}_i\), such that
The quantities \(v_i^{\mathrm{eff}}\) are the “effective velocities” of the normal modes in the fluid, nonlinear versions of the sound velocity. Changing coordinates, one then obtains
It is also possible to rewrite all equations in terms of the coordinates given by the Lagrange parameters \({\underline{\beta }}\). Using (B.3) and (B.5), we find
Further, we also have (if \({\mathsf {C}}\) is invertible)
and using (B.5), we get
1.2 The Riemann Problem in Linear Hydrodynamics
Consider a linear Euler hydrodynamics, with \({\mathsf {A}}_i^{~ j}\) independent of the state. We can use (B.8) with flux Jacobian independent of x, t:
Let us assume that \({\mathsf {A}}_j^{~i}\) is diagonalisable by a similarity transformation (\({\mathsf {C}}\) is invertible). Consider its right eigenvectors:
Then we find normal modes by the linear coordinate transformation
and we have
The Riemann problem is that of solving the Euler hydrodynamic equations with initial conditions giving by two distinct homogeneous states, one on the left, one on the right. We set
Since both the initial condition in the Riemann problem and the Euler equation are invariant under simultaneous scaling \((x,t)\mapsto (\mu x\,\mu t)\), we may assume the solution to the Riemann problem to have this symmetry as well: all functions of space time are functions of \(\xi = x/ t\) only. We therefore obtain
The solution isFootnote 12
where \(n_{k;{\mathrm{l}}}\) and \(n_{k;{\mathrm{r}}}\) are the normal modes in the states specified by the Lagrange parameters \({\underline{\beta }}_{\mathrm{l}}\) and \({\underline{\beta }}_{\mathrm{r}}\), respectively. Let us consider the ray \(\xi =0\), which is the relevant one for the transport statistics problem, and denote \(n_k(0;{\underline{\beta }}_{\mathrm{l}},{\underline{\beta }}_{\mathrm{r}}) = n_k({\underline{\beta }}_{\mathrm{l}},{\underline{\beta }}_{\mathrm{r}})\). Let us consider flow (3.2) for generating the transport statistics in the state specified by \(n_k({\underline{\beta }}_{\mathrm{l}},{\underline{\beta }}_{\mathrm{r}})\), and denote by \(n_k({\underline{\beta }}_{\mathrm{l}},{\underline{\beta }}_{\mathrm{r}};\lambda )\) the normal modes along the flow. Let us finally consider solution (3.17) to the flow problem in the free case. Using (B.13) and (B.12), this is
Using (B.17), we obtain
which shows the general expression of the extended fluctuation relations (3.19).
1.3 Some Aspects of Conformal Hydrodynamics
It is a simple matter to solve the diagonalisation problem for the flux Jacobian (4.13) of conformal hydrodynamics in arbitrary dimension. We find that the following combinations of the rest-frame temperature \(T_\mathrm{rest}\) and the boost \(\theta \) form normal modes:
with effective velocities \(v^{\mathrm{eff}}_\pm \), respectively, as given in (4.14).
It is also a simple matter to find the current generating functions (B.4). One can check that the currents \({\mathsf {j}}_1\) and \({\mathsf {j}}_2\) as given in (4.11) and (4.12) are generated as per (B.4) by the function
Multi-parameter SCGF and Normal Mode Decompositions
In certain cases, even beyond free models, we can obtain an explicit expression for \(F(\lambda )\), where the integral in (3.6) is performed in terms of the normal modes of the Euler hydrodynamics of the model. This result holds in hydrodynamic theories where the generating function for the currents, Eq. (B.4), has a property of separation into normal modes, Eq. (C.1). We do not know yet the full range of theories with this property, but it includes generalised hydrodynamics [46].
Consider the generating function for currents (B.4). Arguments (see below) suggest that this function may in some situations separate into a sum of functions of the normal coordinates:
This decomposition holds in generalised hydrodynamics [46], but it is clear, from (B.20) and (B.21), that it holds in conformal hydrodynamics if and only if \(d=1\) (the “trivial”, linear case). In the cases where it holds, the SCGF is given by
where \(I_\lambda = [0,\lambda )\) if \(\lambda >0\) and \((\lambda ,0]\) if \(\lambda <0\), and the sets \(\lambda _\star ^{(i,\pm )}\) are the turning points of the sign of the effective velocity,
This parallels what is found in generalised hydrodynamics [65].
In order to show (C.2), we calculate, using (C.1) and (B.4),
where in passing from the second to the third line, we used (B.10). Therefore, assuming we are away from the turning points of \(v^{\mathrm{eff}}_i(\lambda )\), we have
in agreement with (3.6). At the turning points of \(v^{\mathrm{eff}}_i\), there are additional delta-function terms. One can verify that the second term in (C.2) exactly cancels these terms.
In fact, using (C.11), we can also write the multi-parameter SCGF \(F({\underline{\lambda }})\) (see Remark 3.1) as a function of the normal coordinates \({\underline{n}} = {\underline{n}}(\lambda )\) in a similar fashion, as
where the sets \( n_\star ^{(i,{\underline{a}})}\) are the turning points of the sign of the effective velocity,
The constant \(F_0\) is such that at the original state \({\underline{n}}(\lambda =0)\), we recover \(F(0)=0\).
We now provide an argument for decomposition (C.1). Let us consider g in (B.4) as a function of the normal coordinates \({\underline{n}}\). We argue, under certain assumptions (which may be hard to verify), that \(\partial g/\partial n_i\) is independent of \(n_j\) for \(j\ne i\). This would imply decomposition (C.1).
For this argument, we consider the multi-parameter SCGF. Let us assume that there is a differentiable multi-parameter flow
and a differentiable SCGF \(F({\underline{\lambda }})\) for the transport of all charges \(Q_{j}\), each associated with \(\lambda _j\), as per Remark 3.1. We combine (3.6) with (C.8) in this general situation, in order to obtain, assuming the matrix \(\partial {\underline{\beta }}/\partial {\underline{\lambda }}\) to be invertible,
where \(\frac{\partial F}{\partial {\underline{\beta }}}\) and \(\frac{\partial g}{\partial {{\underline{\beta }}}}\) are to be seen as line vectors. Changing variables, we have
and using (B.10) this gives
We do not expect differentiability of F at the points where \({{\,\mathrm{sgn}\,}}(v^{\mathrm{eff}})\) changes. However, away from these points, it is natural to assume that F is differentiable. Consider therefore taking in (C.11) another derivative with respect to the normal modes. Away from the points where the effective velocity changes sign, \({{\,\mathrm{sgn}\,}}(v^{\mathrm{eff}}_j)\) has zero derivative, and we obtain
Since the left-hand side is symmetric (by differentiability of F), we find
With \(i\ne j\), in states where the signs of \(v^{\mathrm{eff}}_i\) and \(v^{\mathrm{eff}}_j\) are different, this implies
Suppose that for each \(i>j\), there exists a neighbourhood of states such that \({{\,\mathrm{sgn}\,}}(v^{\mathrm{eff}}_i)\ne {{\,\mathrm{sgn}\,}}(v^{\mathrm{eff}}_j)\). Then (C.14) will hold in all these neighbourhoods, for the corresponding (i, j). Suppose also that for all j and all i, the function \(\partial g/\partial n_j\) is analytic in \(n_i\) (in appropriate neighbourhoods of \(n_i\) such that \({\underline{n}}\) lies in the manifold of MES). Then, by analytic continuation, one would have (C.14) for all \({\underline{n}}\), and therefore (C.1).
One can verify that the multi-parameter flow (C.8) for the flux Jacobian of higher-dimensional CFT (4.13) is not consistent: it does not lead to differentiable Lagrange parameters as functions of the many parameters \(\lambda _i\), at least in the region \(|\theta |<\theta _{\mathrm{s}}\) (and we note that in the region \(|\theta |>\theta _{\mathrm{s}}\), the effective velocities have the same sign). This is, technically, where the above argument fails in this case. It would be interesting to further study this situation.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Doyon, B., Myers, J. Fluctuations in Ballistic Transport from Euler Hydrodynamics. Ann. Henri Poincaré 21, 255–302 (2020). https://doi.org/10.1007/s00023-019-00860-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00023-019-00860-w