Fluctuations in Ballistic Transport from Euler Hydrodynamics

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.


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 → ∞. 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 largedeviation 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 = ξ}). 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.

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).

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, 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 = R dx q i (x, t). These are dynamical observables satisfying dQ i /dt = 0. They are assumed to have associated conservation laws indexed by i. Here q i (x, t) and 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 · · · , 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 β 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 · · · β , where β is the vector of all β i 's, coordinates for the MES manifold (when it is clearer, we will also use the notation β • to represent the set of β i 's for all i).
The form e − i β i Qi for the probability measure or density matrix is formal in most situations. In the simple cases where the series i β 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 i β 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(β), which may be local or quasi-local. The Lagrange parameters β 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. 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 quasilocality, 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 β i 's are more appropriately represented by a function on a "spectral space", a space of available stable "quasiparticles", 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 β i l (left) and β i 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 β i depending on the initial left and right sets, β i (β l , β 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 {β i }.
valuesj in the limit t → ∞. 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 * , (2.7) 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 P is determined by the state of interest, which we take to be a MES for some β. According to the large-deviation principle, such extensive quantities have probability distributions that are exponentially peaked at the almostsure value (here the scaling in time has exponent 1, which is the one relevant here), 3 P(J (t) = tj) e −tI(j) , I(j) = 0, I(j) > 0 (j =j). (2.9) The almost-sure valuej is simply the average current in the state of interest. The function I(j) controlling this exponential is referred to as the largedeviation rate function. It describes the probabilities of rare but significant events where the quantity J (t) deviates by large amounts from tj. The rate function is the Legendre-Frenchel transform of the scaled cumulant generating function (SCGF) F (λ) for J (t) , Here c n are the cumulants, scaled by time. Up to a conventional minus sign, the function F (λ) can be interpreted as a non-equilibrium equivalent of the equilibrium specific free energy, with non-equilibrium partition function Z = e λJ (t) β e tF (λ) . 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 (λ).

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 t 0 ds 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]).

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.

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 q i = q i (0, 0) β , which provide a good system of coordinates for the MES manifold-that is, the map β → q is bijective (from an appropriate space of β). 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, j i = j i (0, 0) β , as functions of the state coordinates 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 β, or 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 β. We define a flow λ → β(λ) on the manifold of MES, starting on this state β(0) = β and for λ lying in some interval of R, by the differential equation where A(λ) is the flux Jacobian in the state with Lagrange parameters β(λ). On the right-hand side, the sign of the flux Jacobian A is the matrix obtained by diagonalising A and taking the sign of its eigenvalues, 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 C(λ) in (3.4) is the flow-dependent static correlation matrix, evaluated in the state with Lagrange parameters β(λ).
Our main result, shown in "Appendix A", is as follows. We identify the flow parameter λ with the conjugate parameter in (2.10), and we have an expression for the SCGF F (λ) directly in terms of the current along the flow: . Thus, the knowledge of the Euler hydrodynamics (giving the flux Jacobian A, and the currents 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 (λ) 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 λ j associated with all currents j j , with In (3.4) and (3.2), we make the replacements λ →λ j and i * →j. This, however, requires multiple differentiability with respect to the parameters λ j , which are nontrivial relations on the matrix 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 (λ), developed in "Appendix C". This agrees with the general expression found in integrable systems (see [65]).

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, multipoint 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 · · · β by a time integral of the current j(0, t). That is, let us construct a family of states · · · (λ) , parametrised by λ ∈ R, with · · · (0) = · · · β , obtained by modifying the state · · · β by the insertion of the time-integrated local current 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 λ 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 λ in an interval containing the origin. In particular, we have where the connected correlation function ·, · (λ),c is defined as in (2.6) but for the state · · · (λ) . 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 β(λ) with β(0) = β. We determine this change solely from the Euler hydrodynamics of the system.
Before determining β(λ), we explain how the λ-dependent state is useful in order to fix F (λ) (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 = ±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 λ with the condition F (0) = 0, we obtain (3.6). The state resulting from the λ-bias is manifestly stationary. Also, by the fact that 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 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 λ modifications, Eq. (3.10), generate tangents to the MES manifold. Since at λ = 0 the state lies on the MES manifold, then it stays on it. As a consequence, λ → · · · (λ) forms a path lying within this manifold. Therefore, there exists β(λ) with β(0) = β and This is the crucial observation of the method. Using (3.10), we then have for any local or quasi-local observable O. We may obtain an equation for the coordinates q i (λ) by specifying O in (3.14) to be the available conserved densities of the model, 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 q(λ). Equation (3.14) therefore fully specifies the path by giving its tangent at the point q(λ) in terms of a function of q(λ). 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 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 λ-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).

Constant Flux Jacobian and Extended Fluctuation Relations
The form of the flux Jacobian A j i 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 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 R-linear transformations, but 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 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 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 λ. 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 β i characterising the steady state, as functions of the set of Lagrange parameters on the left and right halves of the initial state, β i (β l , β 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.

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 Suppose (without loss of generality) that the path has end points (0) = (0, 0) and (1) = (x, t), and denote by = √ 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 θ (with θ ∈ [0, 2π)) determined by (s) = (sx, st) = s (sin θ, cos θ). (3.22) We are interested again in the large-, scaled statistics, and so we must evaluate the expectation value (3.23) 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 which generalises (3.2) to the case θ = 0. With this flow, the result takes a form that generalises (3.6), Consider the limit θ → π/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 β i * proportional to λ, that is, From (3.26), we have in this case ∂F (λ; π/2)/∂λ = −q(λ; π/2), which, along with F (0; π/2) = 0 and (3.27), allows us to identify F (λ; π/2) with a free energy difference. That is, we obtain where Δf (λ) is the specific (dimensionless) free energy difference with the specific free energy being f (β) = − log Z(β), where Z(β) is the partition function for Lagrange parameters β.
The above results have perhaps their most interesting application to the evaluation of correlation functions of twist fields. Let us introduce the "height fields" ϕ(x, t) defined as q(x, t) = ∂ x ϕ(x, t) and j(x, t) = −∂ t ϕ(x, t), which automatically solves the continuity relation (2.1). Differences of height fields ,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. 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 λϕ(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 ϕ(x, t) can be defined independently of the path chosen).
Interestingly, SCGF (2.10) gives rise to the leading exponential behaviour of the dynamical two-point functions of twist fields: . 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 (θ = π/2), (3.30) with (3.28) gives 5 (3.31) In this case, some exact results already exist that confirm (3.31). A specialisation to the 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 with 6 (3.31). is, however, more complicated. It might be tempting to think that the result will separate into two contributions, one from each straight piece, F (λ; θ) being proportional to a sum of F (λ; 0) = F (λ) (the SCGF calculated in Sect. 3.1) and F (λ; π/2) = −Δf (λ) (the difference of free energy densities). However, this is generically incorrect: for instance, because of ballistic transport, the space integral x 0 dx q(x, t), q i (0, 0) c β does not necessarily vanish in the limit x, t → ∞. Likewise, the evaluation of F (λ; θ) for any choice of path that is not straight (at the scale set by → ∞) 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. 5 BD acknowledges discussions with V. Alba at the Perimeter Institute, September 2017, that led to both (3.31) and (3.32). 6 Results (3.31) and (3.32) hold as well for twist fields that are associated with discrete symmetries, as Δf (λ) 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(θ), for the n-copy permutation twist-field two-point function Tn(0, 0)Tn(x, 0) , we have Δfn = dp(θ) log 1+e − n(θ) (1+e − 1 (θ) ) n , with pseudoenergy n(θ) having source term nw(θ) implementing the n-times larger imaginary time direction induced by the twist property, and p(θ) being the momentum function.

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 β(x). By using the idea of local entropy maximisation, one simply expects each "fluid cell" atx to produce a contribution proportional to its free energy difference Δfx(λ), 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.

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 = ξ if there is no effective velocity v eff i that takes the value ξ. 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 eff i = ξ ∀ 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 ξ, 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 ξ and the presence or not of an effective velocity with value ξ. 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 ξ is such that v eff j = ξ 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 λ derivatives. Consider the case ξ = 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 λ = 0 if there is a co-propagating mode (an eigenvalue of 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 λ spacewhere 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 λ is modified, through such points. On both sides of the phase transition point, the cumulants will take different forms, as 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 quasiparticle 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 tj and the fluctuations of (J (t) −tj)/ √ 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 ≥ 3 suggests that the cumulants of (J (t) − tj)/ √ 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), 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.

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 (λ). All these are, to our knowledge, new results. 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 (λ) in such states will be analysed in a future publication.

Reduction to a One-Dimensional Hydrodynamic Problem
Consider conformal field theory in d > 1 dimensions of space, and recall that the energy-momentum tensor T μν satisfies T μν = T νμ (Lorentz invariance), T μ μ = 0 (scale invariance) and ∂ μ T μν = 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 = d d x T 01 (x). In this case, we may restrict to the space of thermal states boosted in that direction [70], with density matrices e −β 1 H−β 2 P = e −βrest(cosh θ H−sinh θ P ) (4.1) and corresponding state denoted by · · · (d) β1,β2 . Here H = d d x T 00 (x) is the Hamiltonian, the rest-frame temperature is T rest = β −1 rest , and the Lorentz boost is of rapidity θ. By relativistic and conformal invariance, expectation values of energy-momentum components take the form where u μ = (cosh θ, sinh θ, 0, . . . , 0) μ 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 ⊥ , with coordinates x ⊥ = (x 2 , . . . , x d ), to be of d − 1-dimensional hyperarea V ⊥ and to be periodic in all its coordinates, with equal periods. Let us denote by · · · (d,V ⊥ ) β1,β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 , 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 j 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 ⊥ , and so we consider By clustering in d + 1-dimensional space-time, all cumulants generated by F (λ) have a finite limit as V ⊥ → ∞, and we look for Clearly, this is obtained by taking the large-V ⊥ limit of the solution presented in Sect. 3.1 for the effectively one-dimensional system, that is, of expression (3.6) divided by V ⊥ , with (3.2). It is then sufficient to know the following large-V ⊥ limits of conserved densities and currents: (4.8) 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 Ann. Henri Poincaré μ, ν ∈ {0, 1}. Thus, we have fully reduced the problem of evaluating (4.7) to a one-dimensional hydrodynamic problem.

Exact SCGF and Cumulants
Using (4.2), we have more explicitly The flux Jacobian takes the form (4.13) As the flux Jacobian is not state independent if d > 1, the extended fluctuation relations do not hold in higher-dimensional CFT. 9 The matrix sgn(A) is obtained by diagonalising A and taking the sign of the eigenvalues. The eigenvalues of A are: with the speed of sound of conformal hydrodynamics given by v s = 1/ √ d. Clearly, sgn(v eff ± ) = sgn(sinh 2θ ± α) = sgn(θ ± θ s ) with the sound rapidity θ s defined by tanh θ s = v s . Define ξ 1 = sgn v eff + + sgn v eff − and ξ 2 = sgn v eff + − sgn v eff − . They take the following values: That is, ξ 1 is nonzero for supersonic rapidities and zero otherwise, and ξ 2 is nonzero for infrasonic rapidities and zero otherwise. We then have sgn A = 1 2α Consider the SCGF for energy transport. We fix i * = 1 in (3.2) and let the Lagrange multipliers become λ dependent. Then using β 1 = β rest cosh θ and β 2 = −β rest sinh θ, 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 λ and setting λ = 0. The derivatives may be readily evaluated: Reading off the required identities from (4.17), setting λ = 0 and using the explicit form of ξ 1/2 given by (4.15), we obtain the second cumulant: where θ = θ(0) and β rest = β rest (0). This process easily generates cumulants of n 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 β rest (λ) and θ(λ), then inserting these expressions into (4.18) and numerically integrating (3.6) to obtain F (λ). 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, θ = ±θ s . In this case, as a consequence of the discontinuities in ξ 1,2 in (4.15), the derivatives of β rest (λ) and θ(λ) with respect to λ have discontinuities at λ = 0. This implies that the third derivative of F (λ) does not exist at λ = 0. That is, the

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 higherdimensional 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. Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

A. 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 ξ = x/t in a neighbourhood of the ξ = 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.

A.1. Statement of the Problem
Consider expression (2.10) for the scaled cumulants c n . Let us assume that correlation functions of the current 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 Let us assume more generally that correlation functions involving the currents j(0, t) and local observables cluster strongly enough at large time differences. Then standard arguments show that (i) the state · · · (λ) defined by the series expansion in λ of (3.9) has a nonzero radius of convergence for any local observable O(x, t) (the radius may depend on the observable), (ii) for any local observable O(x, t), (3.10) holds within the convergence region, (iii) the series expansion in (2.10) for F (λ) 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 (λ), in order to prove (3.6) with (3.2): I. We may define F (λ) as the solution to (3.12) with F (0) = 0, where the state · · · (λ) is defined as the solution to (3.10) with · · · (0) = · · · β . 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.

A.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 C and the flux Jacobian 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 ∈ (0, 1) such that clustering of correlation functions of local observables in time is strong enough along all rays ξ = 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 O(x, t), and the Euler scaling limit S Eul iO (kt) as in (A.4). Then, the stronger version is where the inner product is defined in (2.5) and C is assumed to be invertible. Equation ( This implies in particular (3.16), but is more general.
Proof of (A. 10). Note that (A.9), specialised to O(x, t) = 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 densitỹ for all local observables O(x, t). Once this is shown, specialising it to O(x, t) = q j (x, t) and using (A.10) for O(x, t) = q j (x, t), we determine that V j i = sgn(A) j i , which completes the proof. We first prove (A.11). The main idea is to show that ∞ −∞ dt 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 → J (t) = t/2 −t/2 ds 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 J (t) , J (t) c β has a growth that is bounded linearly in t; (II) that lim t→∞ J (t) , O(x, 0) c exists for all local observables 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 ∈ (0, 1) such that the supremum of the difference | J (t) , O(x, 0) c − J (t) , O(y, 0) c | within the region x, y ∈ [−kt, kt] tends to zero as t → ∞. 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 ξ = x/t with |ξ| ∈ [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 ξ = x/t ∈ [−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) ∝ e −μx 2 for some μ > 0, with the normalisation condition Then using (A.9) as well as the dominated convergence theorem, we have Vol. 21 (2020) Fluctuations in Ballistic Transport 287 The equation also holds with t → −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).
Proof of (A.5) from (A. 8). First note that by symmetry (in particular, recall definition (2.6)), assumption (A.8) implies 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 Eul ij (0) = C ij , is (A.5).

A.3. λ-Flow
We define the "Lie derivative" L · · · β , at the point · · · β in the MES manifold characterised by the Lagrange parameters β, 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 λ → · · · (λ) is that along the direction set by the Lie derivative L · · · β , starting at some point β(0), which lies entirely on the MES manifold. Hence, it can be characterised by λ → β(λ), which solves (3.2): (A.28) Finally, the differential equation (3.12) with F (0) = 0 shows (3.6).
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 β(λ) solves (3.2). We show this by induction. The induction is on the statement that (A. 30) 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 λ-derivative, this gives

A.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| = 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).

B.1. Standard Results
Consider a non-stationary, inhomogeneous state · · · 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 · · · and using the approximation (B.1), one obtains 10 t) . This can be rewritten using the flux Jacobian (3.1), Since 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 −∂j i /∂β j = dx q j (x, 0), j i (0, 0) c β is symmetric under i ↔ j, which can be shown by using the conservation laws and integration by parts. Note that changing variables to q, this symmetry also immediately implies the important relation involving the static correlation matrix C = ∂q ∂β (see (3.5)). In particular, this implies that with respect to the inner product defined by the C matrix, (a, b) = ij a i b j C ij , which is non-degenerate if C is invertible (C is always nonnegative), the matrix A is Hermitian. 11 Hence, if C is invertible, then A is diagonalisable by a similarity transformation.
The normal coordinates of Euler hydrodynamics are a different system of coordinates, q → n, which diagonalise the flux Jacobian. That is, there is a diagonal matrix v eff , with diagonal elements v eff i , such that The quantities v eff i 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 β. Using (B.3) and (B.5), we find (B.9) 11 In this notation, we have in particular q j i = δ j i and j j i = A j i . and using (B.5), we get v eff = ∂n ∂β A T ∂β ∂n . (B.10)

B.2. The Riemann Problem in Linear Hydrodynamics
Consider a linear Euler hydrodynamics, with A j i independent of the state. We can use (B.8) with flux Jacobian independent of x, t: Let us assume that A i j is diagonalisable by a similarity transformation (C is invertible). Consider its right eigenvectors: Then we find normal modes by the linear coordinate transformation and we have 14) 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 The solution is 12 where n k;l and n k;r are the normal modes in the states specified by the Lagrange parameters β l and β r , respectively. Let us consider the ray ξ = 0, which is the relevant one for the transport statistics problem, and denote n k (0; β l , β r ) = n k (β l , β r ). Let us consider flow (3.2) for generating the transport statistics in the state specified by n k (β l , β r ), and denote by n k (β l , β r ; λ) 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 n k (β l , β r ; λ) = n k (β l , β r ) − λ sgn(v eff k )w i * k . (B.18) Using (B.17), we obtain which shows the general expression of the extended fluctuation relations (3.19).

B.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 rest and the boost θ form normal modes: with effective velocities v eff ± , 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 j 1 and j 2 as given in (4.11)

C. Multi-parameter SCGF and Normal Mode Decompositions
In certain cases, even beyond free models, we can obtain an explicit expression for F (λ), 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 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 eff i (λ), we have in agreement with (3.6). At the turning points of v 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 (λ) (see Remark 3.1) as a function of the normal coordinates n = n(λ) in a similar fashion, as The constant F 0 is such that at the original state n(λ = 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 n. We argue, under certain assumptions (which may be hard to verify), that ∂g/∂n i is independent of n j for j = 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 (λ) for the transport of all charges Q j , each associated with λ j , as per Remark 3.1. We combine (3.6) with (C.8) in this general situation, in order to obtain, assuming the matrix ∂β/∂λ to be invertible, We do not expect differentiability of F at the points where sgn(v 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, sgn(v eff j ) has zero derivative, and we obtain ∂ 2 F ∂n i ∂n j = ∂ 2 g ∂n i ∂n j sgn(v eff j ). (C.12) Since the left-hand side is symmetric (by differentiability of F ), we find Suppose that for each i > j, there exists a neighbourhood of states such that sgn(v eff i ) = sgn(v 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 ∂g/∂n j is analytic in n i (in appropriate neighbourhoods of n i such that n lies in the manifold of MES). Then, by analytic continuation, one would have (C.14) for all 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 λ i , at least in the region |θ| < θ s (and we note that in the region |θ| > θ 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.