We devise an iterative scheme for numerically calculating dynamical two-point correlation functions in integrable many-body systems, in the Eulerian scaling limit. Expressions for these were originally derived in Ref. [1] by combining the fluctuation-dissipation principle with generalized hydrodynamics. Crucially, the scheme is able to address non-stationary, inhomogeneous situations, when motion occurs at the Euler-scale of hydrodynamics. In such situations, in interacting systems, the simple correlations due to fluid modes propagating with the flow receive subtle corrections, which we test. Using our scheme, we study the spreading of correlations in several integrable models from inhomogeneous initial states. For the classical hard-rod model we compare our results with Monte-Carlo simulations and observe excellent agreement at long time-scales, thus providing the first demonstration of validity for the expressions derived in Ref. [1]. We also observe the onset of the Euler-scale limit for the dynamical correlations.
We devise an iterative scheme for numerically calculating dynamical two-point correlation functions in integrable many-body systems, in the Eulerian scaling limit. Expressions for these were originally derived in Ref. [1] by combining the fluctuation-dissipation principle with generalized hydrodynamics. Crucially, the scheme is able to address non-stationary, inhomogeneous situations, when motion occurs at the Euler-scale of hydrodynamics. In such situations, in interacting systems, the simple correlations due to fluid modes propagating with the flow receive subtle corrections, which we test. Using our scheme, we study the spreading of correlations in several integrable models from inhomogeneous initial states. For the classical hard-rod model we compare our results with Monte-Carlo simulations and observe excellent agreement at long time-scales, thus providing the first demonstration of validity for the expressions derived in Ref. [1]. We also observe the onset of the Euler-scale limit for the dynamical correlations.
Using the theory of generalized hydrodynamics (GHD), we derive exact Euler-scale dynamical two-point correlation functions of conserved densities and currents in inhomogeneous, non-stationary states of many-body integrable systems with weak space-time variations. This extends previous works to inhomogeneous and non-stationary situations. Using GHD projection operators, we further derive formulae for Euler-scale two-point functions of arbitrary local fields, purely from the data of their homogeneous one-point functions. These are new also in homogeneous generalized Gibbs ensembles. The technique is based on combining a fluctuation-dissipation principle along with the exact solution by characteristics of GHD, and gives a recursive procedure able to generate n n -point correlation functions. Owing to the universality of GHD, the results are expected to apply to quantum and classical integrable field theory such as the sinh-Gordon model and the Lieb-Liniger model, spin chains such as the XXZ and Hubbard models, and solvable classical gases such as the hard rod gas and soliton gases. In particular, we find Leclair-Mussardo-type infinite form-factor series in integrable quantum field theory, and exact Euler-scale two-point functions of exponential fields in the sinh-Gordon model and of powers of the density field in the Lieb-Liniger model. We also analyse correlations in the partitioning protocol, extract large-time asymptotics, and, in free models, derive all Euler-scale n n -point functions.
At large scales of space and time, the nonequilibrium dynamics of local observables in extensive many-body systems is well described by hydrodynamics. At the Euler scale, one assumes that each mesoscopic region independently reaches a state of maximal entropy under the constraints given by the available conservation laws. Away from phase transitions, maximal entropy states show exponential correlation decay, and independence of fluid cells might be assumed to subsist during the course of time evolution. We show that this picture is incorrect: under ballistic scaling, regions separated by macroscopic distances "develop long-range correlations as time passes." These correlations take a universal form that only depends on the Euler hydrodynamics of the model. They are rooted in the large-scale motion of interacting fluid modes and are the dominant long-range correlations developing in time from long-wavelength, entropy-maximized states. They require "the presence of interaction" and "at least two different fluid modes" and are of a fundamentally different nature from well-known long-range correlations occurring from diffusive spreading or from quasiparticle excitations produced in far-from-equilibrium quenches. We provide a universal theoretical framework to exactly evaluate them, an adaptation of the macroscopic fluctuation theory to the Euler scale. We verify our exact predictions in the hard-rod gas, by comparing with numerical simulations and finding excellent agreement.
We derive an exact formula for the scaled cumulant generating function of the time-integrated current associated to an arbitrary ballistically transported conserved charge. Our results rely on the Euler-scale description of interacting, many-body, integrable models out of equilibrium given by the generalized hydrodynamics, and on the large deviation theory. Crucially, our findings extend previous studies by accounting for inhomogeneous and dynamical initial states in interacting systems. We present exact expressions for the first three cumulants of the time-integrated current. Considering the non-interacting limit of our general expression for the scaled cumulant generating function, we further show that for the partitioning protocol initial state our result coincides with previous results of the literature. Given the universality of the generalized hydrodynamics, the expression obtained for the scaled cumulant generating function is applicable to any interacting integrable model obeying the hydrodynamic equations, both classical and quantum.
We consider a one-dimensional gas of hard rods, one of the simplest examples of an interacting integrable model. It is well known that the hydrodynamics of such integrable models can be understood by viewing the system as a gas of quasiparticles. Here, we explore the dynamics of individual quasiparticles for a variety of initial conditions of the background gas. The mean, variance, and two-time correlations are computed exactly and lead to a picture of quasiparticles as drifting Brownian particles. For the case of a homogeneous background, we show that the motion of two tagged quasiparticles is strongly correlated, and they move like a rigid rod at late times. Apart from a microscopic derivation based on the mapping to point particles, we provide an alternate derivation which emphasizes that quasiparticle fluctuations are related to initial phase-space fluctuations, which are carried over in time by Euler scale dynamics. For the homogeneous state, we use the Brownian motion picture to develop a Dean-Kawasaki-type fluctuating hydrodynamic theory, formally having the same structure as that derived recently by Ferrari and Olla. We discuss differences with existing proposals on the hydrodynamics of hard rods and some puzzles.
We introduce a new universal framework describing fluctuations and correlations in quantum and classical many-body systems, at the Euler hydrodynamic scale of space and time. The framework adapts the ideas of the conventional macroscopic fluctuation theory (MFT) to systems that support ballistic transport. The resulting “ballistic MFT” (BMFT) is solely based on the Euler hydrodynamics data of the many-body system. Within this framework, mesoscopic observables are classical random variables depending only on the fluctuating conserved densities, and Euler-scale fluctuations are obtained by deterministically transporting thermodynamic fluctuations via the Euler hydrodynamics. Using the BMFT, we show that long-range correlations in space generically develop over time from long-wavelength inhomogeneous initial states in interacting models. This result, which we verify by numerical calculations, challenges the long-held paradigm that at the Euler scale, fluid cells may be considered uncorrelated. We also show that the Gallavotti-Cohen fluctuation theorem for non-equilibrium ballistic transport follows purely from time-reversal invariance of the Euler hydrodynamics. We check the validity of the BMFT by applying it to integrable systems, and in particular the hard-rod gas, with extensive simulations that confirm our analytical results.
One–dimensional integrable and integrable systems display, on macroscopic scales, a universal form of transport known as (GHD). In its standard Euler-scale formulation, GHD mirrors the equations of a two-dimensional compressible fluid but ignores fluctuations and becomes numerically unwieldy as soon as integrability-breaking perturbations are introduced. We show that GHD can be efficiently simulated as a gas of —a natural generalization of hard-rod particles—whose trajectories are efficiently mapped onto those of point particles. This representation (i) provides a transparent route to incorporate integrability-breaking terms, and (ii) automatically embeds the exact fluctuating-hydrodynamics extension of GHD. The resulting framework enables fast, large-scale simulations of quasi-integrable systems even in the presence of complicated integrability-breaking perturbations. It also manifests the pivotal role of two-point correlations in systems confined by external potentials: we demonstrate that situations where local one-point observables appear thermalized can nevertheless sustain long-lived, far-from-equilibrium correlations for arbitrarily long times, signaling that, differently from what was previously stated, true thermalization is reached at diffusive time-scales.