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 nonequilibrium states is entropy production and the breaking of time-reversal invariance. This is associated with the existence non-zero 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 nonequilibrium phenomena.
In the context of nonequilibrium 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 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+1dimensional conformal field theory (CFT) [23,24], see the review [25]. These systems admit ballistic transport, and nonequilibrium currents are generated by the partitioning protocol [26,27,28] (see also [25] and references therein), where an imbalance exists in the initial condition. 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 nonequilibrium steady states. In particular, they include the nonequilibrium steady states emerging in the partitioning protocol.
In the setup we consider, time evolution is deterministic, and initial states are fluctuating. In a hydrodynamic description [46], 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 [61, App A], [47] and [48], 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 [49,50] 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 [51].
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. [52,53,54,55,56,57,58,59,60]. We explain why this occurs from hydrodynamic principles, and how this may connect with the breaking of Gaussianity found in nonlinear fluctuating hydrodynamics [61,47,62].
Of immediate importance is the application to integrable systems, a detailed examination of which is completed in a separate work [63]. There, we exploit generalised hydrodynamics [64,65,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. Nonequilibrium 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 section 2, we explain the context and review the main aspects of large deviation theory for transport. In section 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 section 4 we present the application to conformal hydrodynamics, and in section 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 both in 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 both to 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 section 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 [75,43]) 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. [76,77,43].
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 e − i β i Q i . 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 Q i for the probability measure or density matrix is formal in most situations. When 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] for discussions. A formulation related to the latter, which is developed in [43], accounts for "generalised thermalisation" in integrable systems, where the series i β i Q i may be infinite and the charges quasi-local. It is based on considering conserved charges Q i as vectors in the tangent space to the MES manifold. A charge Q i deforms the state according to the equation where the symmetric inner product on the space of local fields is and refers to the connected, symmetrised correlation function 2 . According to [43], the tangent space is the Hilbert space, induced by (2.3), of pseudolocal conserved charges [78,79,75], and we believe that the results below hold if the set Q i is complete in the sense that it spans a dense subset of this Hilbert space. 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], 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 [80]. 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 [46], 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, 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 [64,65,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 "quasi-particles", via the thermodynamic Bethe ansatz [81,82,83], see e.g. [75,84]; thus in these cases the MES manifold is a manifold of functions, whose full description is however not yet known in most models.
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.
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 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.

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 in the limit t → ∞. For the purpose of studying nonequilibrium 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 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 almost-sure value (here the scaling in time has exponent 1, which is the one relevant here) 3 , P(J (t) = tj) e −tI(j) , I() = 0, I(j) > 0 (j =). (2.7) The almost-sure value is simply the average current in the state of interest. The function I(j) controlling this exponential is referred to as the large-deviation function. It describes the probabilities of rare but significant events where the quantity J (t) deviates by large amounts from t. Its Legendre-Frenchel transform, F (λ), is the scaled cumulant generating function (SCGF) 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 nonequilibrium equivalent of the equilibrium specific free energy, with nonequilibrium partition function Z = e λJ (t) β e tF (λ) . Crucially, the largedeviation principle (2.7) 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.6), (2.7), and (2.8), see e.g. the review [25]. In this paper, we assume this to be the case.

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 [81,82,83]).
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, is the diagonal matrix of eigenvalues (the v 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.5).
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.8), and we have an expression for the SCGF F (λ) directly in terms of the current along the flow: where j(λ) = j i * (λ) = j i * (0, 0) β(λ) . 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 [47,48]. 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 to 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 opbtain a more explicit expression for F (λ), developed in Appendix C. This agrees with the general expression found in integrable systems, see [63].

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 [75,43], or from a strong version of the hydrodynamic projection principles [46,61,47,48]. 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 the 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, the 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 manybody 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. Although this vanishing is expected generically, contrary to that of connected correlation functions at large space separations, it is nevertheless not guaranteed in MES. Its breaking leads to a failure of the large deviation principle (2.7). We discuss the physics and potential consequences of this in subsection 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.4) 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 timeintegrated 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 [49,50]. 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 the definition (2.8), 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, 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, we have for any x, 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, hence can be seen as a function of the state coordinates q(λ). Eq. (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 [46,61,47,48] is that, in an appropriate Euler scaling limit [46,48,85] (see appendix A), where in particular both x and t are large in fixed ratio, we have 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 λ-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 the 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 nonequilibrium steady states emerging from the partitioning protocol [26,27,28]. 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, and if ballistic transport is supported, nonequilibrium currents may emerge. Let us denote the Lagrange parameters β i characterising this steady state by β i (β l , β r ). In fluid dynamics, this is known as the Riemann problem [86].
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+1dimensional 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 [63], 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 subsection 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.6), the quantity , dt(s)) = ds (dx/ds, dt/ds) is the infinitesimal tangent to the path, j = (j, q) is the conserved current vector, and 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, the 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 Extending the arguments presented in appendix A, or the simpler derivation of section 3.2, it is a simple matter to derive, for the path (3.22), the flow (recall (2.5)) 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, the 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 ϕ(x 2 , t) − ϕ(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 O(x, t) ∝ e λϕ(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 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 bosonization and has applications to the transverse field Ising model and XXZ chains, see e.g. [87,88]. Certain classes of such twist fields can also be used to study entanglement entropy in free-particle models [89]. Twist fields associated to spacetime symmetries have also been studied recently [90]; 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, the SCGF (2.8) 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. [88]). Formula 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 [91], and shown to agree with results derived from special quantum quenches in the Ising model [92]. In the context of entanglement entropy, taking into account the twist field interpretation of the Rényi entanglement entropy as expressed in [89], the exact results of [93,94,95,96,97] 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).  1)x, t), s ∈ [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 (λ; θ) being proportional to a sum of F (λ; 0) = F (λ) (the SCGF calculated in subsection 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 subsection 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 β(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 [98], and was verified against direct numerical calculations. This is also connected to formulae for entanglement entropies in inhomogeneous states [99]. 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 [85] 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 the SCGF (2.8) may be divergent (in which case the large deviation principle expressed in (2.7) fails), and what the meaning of this may be. We first note that the scaled cumulants of (2.8) are time-integrated, connected, multipoint 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 subsection 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 [80] (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 [85,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 occur 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.7) does not hold if there exists such a normal mode propagating along the ray -a "co-propagating mode". In fact it is possible to argue from our explicit results that this is the case. 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. From (3.2), the second cumulant also is expected to be finite, although its exact value is ambiguous (in the present theory) if there is a co-propagating mode, because of the ambiguity of the sign function at 0. However, for the third cumulant, we need to take another derivative, which does not exist if there is such a co-propagating mode that couples to the current of interest (i.e. such that the corresponding element of the A is nonzero), because of the discontinuity of the sign function. In this case, all scaled cumulants c n for n ≥ 3 do not exist. Because a discontinuity points to a diverging derivative, this indicates that these higher cumulants are actually divergent.
The divergence of scaled cumulants may be interpreted as a "dynamical phase transition", a concept widely studied in the nonequilibrium large-deviation theory of stochastic models, see for instance [52,53,54,55,56,57,58,59,60]. Indeed, we see a change in the fluctuation spectrum for a dynamical quantity as a parameter of the state is modified, or as the bias λ is modified. On both sides of the phase transition point, the cumulants will take different forms, as sgn(A) changes discontinuously.
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. We claim that if a mode is non-isolated, then no discontinuity appears due to the ensuing smoothing, and no divergence emerges. This is what we observe for transport of generic local conserved quantities in integrable systems [63], 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 [100], 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.8) -is an extension of the law of large numbers: subtracting the average t, the fluctuations of (J (t) − t)/ √ 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.8)). The divergence of the scaled cumulants of largedeviation theory c n , n ≥ 3 suggest that the cumulants of (J (t) − t)/ √ t might no longer be vanishing, thus breaking Gaussianity.
Crucially, this may have a connection with nonlinear fluctuating hydrodynamics [61,47,62]. 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. Nonlinear fluctuating hydrodynamics describes more than the particular current fluctuations studied here. Nevertheless, the present theory is in agreement with it, and provides an additional confirmation of the breaking of Gaussianity. In particular, with [63] and the above discussion, this shows that, as already predicted in [61], 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 XXZ chain (as [100] indicates that there is an isolated mode), 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 subsection 3.1 applies to all systems in the class described in subsection 2.1. Although the focus was on one dimension, the formalism applies as well to effectively one-dimensional setups 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 the 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 7 . We further note that conformal hydrodynamics was used in order to obtain exact nonequilibrium 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 energymomentum 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) 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 the form (4.1); in particular, the limit of infinite transverse hyperarea reproduces the infinite-volume results, We then define one-dimensional densities as (x, t) comes from T 01 (x, x ⊥ , t) = T 10 (x, x ⊥ , 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]), onedimensional Euler hydrodynamics apply. We may then ask about the SCGF for the energy current j as defined in (2.8), and the general discussion and results of section 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 subsection 3.1 for the effectively one-dimensional system; that is, of the 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: and to solve for the 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 µ, ν ∈ {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 As the flux Jacobian is not state-independent if d > 1, the extended fluctuation relations do not hold in higher-dimensional CFT 8 . 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 the 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: , |θ| > θ s cosh(θ) ((3d + 1) cosh(2θ) − (d + 1)), |θ| < θ s (4.20) 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 [101,102,103,104,105,106,107,21] hold. We leave an in-depth analysis to a future work. Finally, we note that the phenomenon discussed in subsection 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 scaled cumulants c n for n ≥ 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).

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 [63] 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 spacio-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 nonequilibrium 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 [108]. 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 [109]).
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 initialstate 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 [51], 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 [110], 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

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, equation (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. [80,111,77]. Strong enough clustering in time is more difficult to prove, but expected to hold generically in many systems and statesas explained in subsection 3.5, it is broken in some situations, leading to a breaking of the large deviation principle (2.7).
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 [78,79] (see [75]), form a Hilbert space that describe the tangent to the MES manifold as per (2.2). 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 MESs (which are in [43] referred to as GGEs) is used. This definition is in agreement with the one used here, in particular with (2.2), 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, equation (A.8). Hydrodynamic projection principles are nontrivial, but have been used successfully in the context of statistical fluid dynamics [46,61,47,48], see also [85] where an even stronger version of (A.8) is explicitly used to derive correlation results verified numerically in [112].
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 the expression (2.8) 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 where the limits are taken, in order, on t n−1 , t n−2 , . . . , t 1 , as written. Here the manypoint connected, symmetrised correlation function on the right-hand side is a natural generalisation of (2. 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.8) for F (λ) has a nonzero radius of convergence, and (iv) the 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 the cumulants (A.1), along with the second equation in (2.8), as our working expression.
Below we give simple proofs (3.6) with (3.2) under both definitions. The arguments presented in subsection 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, long-time limit. This is in fact a standard result in linear fluctuating hydrodynamics.
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. 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.2)) is the conserved subspace of the Hilbert space completion of the inner product induced by (2.4), 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.6)) 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 [113,114] that the support is only exponentially accurate. But by exponential accuracy, it is possible, using the techniques of [114], 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] tend 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 (A.14) 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 (A.18) 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 19) 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.4)), assumption (A.8) implies where S Eul Oi (kt) is the Euler scaling limit of S Oi (k, t) = ∞ −∞ dx e −ikx O(x, t)q i (0, 0) c β . 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), 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 subsection 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.4 Second cumulant
The second cumulant was shown in [47,48] 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 Euler hydrodynamics 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, that 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 9 where q i (x, t) = q i (0, 0) β(x,t) and j i (x, t) = q i (0, 0) β(x,t) . This can be re-written 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 [46,61,47,48] (see also the explicit proof in [64]) 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)). 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 re-write all equations in terms of the coordinates given by the Lagrange parameters β. Using (B.3) and (B.5), we find

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. Consider its right-eigenvectors: Then we find normal modes by the linear coordinate transformation and we have ∂ t n k (x, t) + v eff k ∂ x n k (x, t) = 0. (B.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 10 n k (ξ; β l , β r ) = n k;l (ξ < v eff k ) n k;r (ξ > v eff k ) (B.17) 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 the 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 the 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.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 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) below. We do not know yet the full range of theories with this property, but it includes generalised hydrodynamics [64]. Consider the generating function for the 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 [64], 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 λ = [0, λ) if λ > 0 and (λ, 0] if λ < 0, and the sets λ (i,±) are the turning points of the sign of the effective velocity, This parallels what is found in generalised hydrodynamics [63]. 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) below, 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 the 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 the 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 to λ 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 differenitable. 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 higherdimensional 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.