Abstract The Toda lattice is an integrable system and its natural space-time stationary states are the generalized Gibbs ensembles (GGE). Of particular physical interest are then the space-time correlations of the conserved fields. To leading order they scale ballistically. We report on the exact solution of the respective generalized hydrodynamic equations linearized around a GGE as background state. Thereby we obtain a concise formula for the family of scaling functions.
Abstract The Toda lattice is an integrable system and its natural space-time stationary states are the generalized Gibbs ensembles (GGE). Of particular physical interest are then the space-time correlations of the conserved fields. To leading order they scale ballistically. We report on the exact solution of the respective generalized hydrodynamic equations linearized around a GGE as background state. Thereby we obtain a concise formula for the family of scaling functions.
Abstract In recent years, a lot of effort has been put in describing the hydrodynamic behavior of integrable systems. In this paper, we describe such picture for the Volterra lattice. Specifically, we are able to explicitly compute the susceptibility matrix and the current-field correlation matrix in terms of the density of states of the Volterra lattice endowed with a Generalized Gibbs ensemble. Furthermore, we apply the theory of linear Generalized Hydrodynamics (GHDs) to describe the Euler scale behavior of the correlation functions. We anticipate that the solution to the GHDs equations develops shocks at ξ 0 = x t ; so this linear approximation does not fully describe the behavior of correlation functions. Intrigued by this fact, we performed several numerical investigations which show that, exactly when the solution to the hydrodynamic equations develops shock, the correlation functions show an highly oscillatory behavior. In view of this empirical observation, we believe that at this point ξ 0 the diffusive contribution are not sub-leading corrections to the ballistic transport, but they are of the same order.
Abstract In many-body systems, the dynamics is governed, at large scales of space and time, by the hydrodynamic principle of projection onto the conserved densities admitted by the model. This is formalised as local relaxation of fluctuations in the Ballistic Macroscopic Fluctuation Theory, and is a nonlinear version of the Boltzmann-Gibbs principle. We use it to derive a projection formula, expressing n -point connected correlation functions (cumulants) of local observables at different space-time points, in terms of those of conserved densities. This applies in every d≥1 d ≥ 1 spatial dimensions and under the ballistic scaling of space and time, both in and out of equilibrium. It generalises the well-known linear-response principle for 2-point functions. For higher-point functions, one needs to account for nonlinear fluctuations of conserved densities and, correspondingly, higher derivatives of local averages. Using Malyshev’s formula for the cumulant expansion, and keeping the leading order, the result is a nonlinear projection, expressed as a sum of products of correlation functions of conserved densities with equilibrium multivariances as coefficients. The sum is combinatorially organised via certain covers of the set of space-time points, which we call “minimal connected covers”. We use this in order to get general, explicit formulas for two- and three-point functions in stationary states, expressed in terms of thermodynamic and Euler-scale data.
Abstract The generalized Gibbs ensemble (GGE) was introduced ten years ago to describe observables in isolated integrable quantum systems after equilibration. Since then, the GGE has been demonstrated to be a powerful tool to predict the outcome of the relaxation dynamics of few-body observables in a variety of integrable models, a process we call generalized thermalization. This review discusses several fundamental aspects of the GGE and generalized thermalization in integrable systems. In particular, we focus on questions such as: which observables equilibrate to the GGE predictions and who should play the role of the bath; what conserved quantities can be used to construct the GGE; what are the differences between generalized thermalization in noninteracting systems and in interacting systems mappable to noninteracting ones; why is it that the GGE works when traditional ensembles of statistical mechanics fail. Despite a lot of interest in these questions in recent years, no definite answers have been given. We review results for the XX model and for the transverse field Ising model. For the latter model, we also report original results and show that the GGE describes spin–spin correlations over the entire system. This makes apparent that there is no need to trace out a part of the system in real space for equilibration to occur and for the GGE to apply. In the past, a spectral decomposition of the weights of various statistical ensembles revealed that generalized eigenstate thermalization occurs in the XX model (hard-core bosons). Namely, eigenstates of the Hamiltonian with similar distributions of conserved quantities have similar expectation values of few-spin observables. Here we show that generalized eigenstate thermalization also occurs in the transverse field Ising model.
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.
Using generalized hydrodynamics (GHD), we develop the Euler hydrodynamics of classical integrable field theory. Classical field GHD is based on a known formalism for Gibbs ensembles of classical fields, that resembles the thermodynamic Bethe ansatz of quantum models, which we extend to generalized Gibbs ensembles (GGEs). In general, GHD must take into account both solitonic and radiative modes of classical fields. We observe that the quasi-particle formulation of GHD remains valid for radiative modes, even though these do not display particle-like properties in their precise dynamics. We point out that because of a UV catastrophe similar to that of black body radiation, radiative modes suffer from divergences that restrict the set of finite-average observables; this set is larger for GGEs with higher conserved charges. We concentrate on the sinh-Gordon model, which only has radiative modes, and study transport in the domain-wall initial problem as well as Euler-scale correlations in GGEs. We confirm a variety of exact GHD predictions, including those coming from hydrodynamic projection theory, by comparing with Metropolis numerical evaluations.