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 theorem along with the exact solution by characteristics of GHD, and gives a recursive procedure able to generate $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 analyze correlations in the partitioning protocol, extract large-time asymptotics, and, in free models, derive all Euler-scale $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.

Generalized hydrodynamics (GHD) is a large-scale theory for the dynamics of many-body integrable systems. It consists of an infinite set of conservation laws for quasi-particles traveling with effective ("dressed") velocities that depend on the local state. We show that these equations can be recast into a geometric dynamical problem. They are conservation equations with state-independent quasi-particle velocities, in a space equipped with a family of metrics, parametrized by the quasi-particles' type and speed, that depend on the local state. In the classical hard rod or soliton gas picture, these metrics measure the free length of space as perceived by quasi-particles, in the quantum picture, they weigh space with the density of states available to them. Using this geometric construction, we find a general solution to the initial value problem of GHD, in terms of a set of integral equations where time appears explicitly. These integral equations are solvable by iteration and provide an extremely efficient solution algorithm for GHD.

Based on the method of hydrodynamic projections we derive a concise formula for the Drude weight of the repulsive Lieb-Liniger $\delta$-Bose gas. Our formula contains only quantities which are obtainable from the thermodynamic Bethe ansatz. The Drude weight is an infinite-dimensional matrix, or bilinear functional: it is bilinear in the currents, and each current may refer to a general linear combination of the conserved charges of the model. As a by-product we obtain the dynamical two-point correlation functions involving charge and current densities at small wavelengths and long times, and in addition the scaled covariance matrix of charge transfer. We expect that our formulas extend to other integrable quantum models.

Using an extension of the concept of twist field in QFT to space-time (external) symmetries, we study conical twist fields in two-dimensional integrable QFT. These create conical singularities of arbitrary excess angle. We show that, upon appropriate identification between the excess angle and the number of sheets, they have the same conformal dimension as branch-point twist fields commonly used to represent partition functions on Riemann surfaces, and that both fields have closely related form factors. However, we show that conical twist fields are truly different from branch-point twist fields. They generate different operator product expansions (short distance expansions) and form factor expansions (large distance expansions). In fact, we verify in free field theories, by re-summing form factors, that the conical twist fields operator product expansions are correctly reproduced. We propose that conical twist fields are the correct fields in order to understand null polygonal Wilson loops/gluon scattering amplitudes of planar maximally supersymmetric Yang-Mills theory.

We define a simple model of conformal field theory in random space-time environments, which we refer to as stochastic conformal field theory. This model accounts for the effects of dilute random impurities in strongly interacting critical many-body systems. On one hand, surprisingly, although impurities are separated by macroscopic distances, we find that the infinite-time steady state is factorized on microscopic lengths, a signature of the emergence of localization. The stationary state also displays vanishing energy current and strong uncorrelated spatial fluctuations of local observables. On the other hand, at finite times, the transient shows a crossover from ballistic to diffusive energy propagation. In this regime and a Markovian limit, concentrating on current-generating initial states with a temperature imbalance, we show that the energy current and density satisfy simple dissipative hydrodynamic equations. We describe the space-time scales at which non-equilibrium currents exist. We show that a light-cone effect subsists in the presence of impurities although a momentum burst propagates transiently on a diffusive scale only.

The theory of generalized hydrodynamics (GHD) was recently developed as a new tool for the study of inhomogeneous time evolution in many-body interacting systems with infinitely many conserved charges. In this letter, we show that it supersedes the widely used conventional hydrodynamics (CHD) of one-dimensional Bose gases. We illustrate this by studying "nonlinear sound waves" emanating from initial density accumulations in the Lieb-Liniger model. We show that, at zero temperature and in the absence of shocks, GHD reduces to CHD, thus for the first time justifying its use from purely hydrodynamic principles. We show that sharp profiles, which appear in finite times in CHD, immediately dissolve into a higher hierarchy of reductions of GHD, with no sustained shock. CHD thereon fails to capture the correct hydrodynamics. We establish the correct hydrodynamic equations, which are finite-dimensional reductions of GHD characterized by multiple, disjoint Fermi seas. We further verify that at nonzero temperature, CHD fails at all nonzero times. Finally, we numerically confirm the emergence of hydrodynamics at zero temperature by comparing its predictions with a full quantum simulation performed using the NRG-TSA-ABACUS algorithm. The analysis is performed in the full interaction range, and is not restricted to either weak- or strong-repulsion regimes.

The evaluation of vacuum expectation values (VEVs) in massive integrable quantum field theory (QFT) is a nontrivial renormalization-group "connection problem" -- relating large and short distance asymptotics -- and is in general unsolved. This is particularly relevant in the context of entanglement entropy, where VEVs of branch-point twist fields give universal saturation predictions. We propose a new method to compute VEVs of twist fields associated to continuous symmetries in QFT. The method is based on a differential equation in the continuous symmetry parameter, and gives VEVs as infinite form-factor series which truncate at two-particle level in free QFT. We verify the method by studying U(1) twist fields in free models, which are simply related to the branch-point twist fields. We provide the first exact formulae for the VEVs of such fields in the massive uncompactified free boson model, checking against an independent calculation based on angular quantization. We show that logarithmic terms, overlooked in the original work of Callan and Wilczek [Phys. Lett. B333 (1994)], appear both in the massless and in the massive situations. This implies that, in agreement with numerical form-factor observations by Bianchini and Castro-Alvaredo [Nucl. Phys. B913 (2016)], the standard power-law short-distance behavior is corrected by a logarithmic factor. We discuss how this gives universal formulae for the saturation of entanglement entropy of a single interval in near-critical harmonic chains, including log log corrections.

Inspired by recent results on the non-equilibrium dynamics of many-body quantum systems, we study the classical hard rod problem in one dimension with initial domain wall condition. Hard rods are an integrable system, in the sense that for each velocity the density of particles is locally conserved. It was proven by Boldrighini, Dobrushin and Sukhov (1983) that on the hydrodynamic space-time scale, the fluid of hard rods satisfies Euler-type equations which comprise all conservation laws. We provide the general solution to these equations on the line, with an initial condition where the left and right halves are, asymptotically, in different states. The solution is interpreted as being composed of a continuum of contact discontinuities, one for each velocity. This is a classical counterpart of the transport problem solved recently in quantum integrable systems. We provide two independent derivations, one paralleling that in the quantum case, the other using explicitly the microscopic dynamics. Finally, we discuss the Navier-Stokes (viscous) corrections, and study its effect on the broadening of the contact discontinuity and on entropy production.

We show irreversibility of the renormalization group flow in non-unitary but ${\cal PT}$-invariant quantum field theory in two space-time dimensions. In addition to unbroken $\mathcal{PT}$-symmetry and a positive energy spectrum, we assume standard properties of quantum field theory including a local energy-momentum tensor and relativistic invariance. This generalizes Zamolodchikov's $c$-theorem to ${\cal PT}$-symmetric hamiltonians. Our proof follows closely Zamolodchikov's arguments. We show that a function $c_{\mathrm{eff}}(s)$ of the renormalization group parameter $s$ exists which is non-negative and monotonically decreasing along renormalization group flows. Its value at a critical point is the "effective central charge" entering the specific free energy. At least in rational models, this equals $c_{\mathrm{eff}}=c-24\Delta$, where $c$ is the central charge and $\Delta$ is the lowest primary field dimension in the conformal field theory which describes the critical point.

Dynamical equations in generalized hydrodynamics (GHD), a hydrodynamic theory for integrable quantum systems at the Euler scale, take a rather simple form, even though an infinite number of conserved charges are taken into account. We show a remarkable quantum-classical equivalence: we demonstrate the equivalence between the equations of GHD, and the Euler-scale hydrodynamic equations of a new family of classical gases which generalize the gas of hard rods. In this family, the "quasi-particles", upon colliding, jump forward or backward by a distance that depends on their velocities, generalizing the jump forward by the rods' length of the fixed-velocity tracer upon elastic collision of two hard rods. Such velocity-dependent position shifts are characteristics of classical soliton scattering. The emerging hydrodynamics of a quantum integrable model is therefore that of the classical gas of its solitons. This provides a "molecular dynamics" for GHD which is numerically efficient and flexible. This is directly applicable, for instance, to the study of inhomogeneous dynamics in integrable quantum chains and in the Lieb-Liniger model realized in cold-atom experiments.

Generalized hydrodynamics (GHD) was proposed recently as a formulation of hydrodynamics for integrable systems, taking into account infinitely-many conservation laws. In this note we further develop the theory in various directions. By extending GHD to all commuting flows of the integrable model, we provide a full description of how to take into account weakly varying force fields, temperature fields and other inhomogeneous external fields within GHD. We expect this can be used, for instance, to characterize the non-equilibrium dynamics of one-dimensional Bose gases in trap potentials. We further show how the equations of state at the core of GHD follow from the continuity relation for entropy, and we show how to recover Euler-like equations and discuss possible viscosity terms.

Recently, it was understood that modified concepts of locality played an important role in the study of extended quantum systems out of equilibrium, in particular in so-called generalized Gibbs ensembles. In this paper, we rigorously study pseudolocal charges and their involvement in time evolutions and in the thermalization process of arbitrary states with strong enough clustering properties. We show that the densities of pseudolocal charges form a Hilbert space, with inner product determined by thermodynamic susceptibilities. Using this, we define the family of pseudolocal states, which are determined by pseudolocal charges. This family includes thermal Gibbs states at high enough temperatures, as well as (a precise definition of) generalized Gibbs ensembles. We prove that the family of pseudolocal states is preserved by finite time evolution, and that, under certain conditions, the stationary state emerging at infinite time is a generalized Gibbs ensemble with respect to the evolution dynamics. If the evolution dynamics does not admit any conserved pseudolocal charges other than the evolution Hamiltonian, we show that any stationary pseudolocal state with respect to this dynamics is a thermal Gibbs state, and that Gibbs thermalization occurs. The framework is that of translation-invariant states on hypercubic quantum lattices of any dimensionality (including quantum chains) and finite-range Hamiltonians, and does not involve integrability.

Understanding the general principles underlying strongly interacting quantum states out of equilibrium is one of the most important tasks of current theoretical physics. With experiments accessing the intricate dynamics of many-body quantum systems, it is paramount to develop powerful methods that encode the emergent physics. Up to now, the strong dichotomy observed between integrable and non-integrable evolutions made an overarching theory difficult to build, especially for transport phenomena where space-time profiles are drastically different. We present a novel framework for studying transport in integrable systems: hydrodynamics with infinitely-many conservation laws. This bridges the conceptual gap between integrable and non-integrable quantum dynamics, and gives powerful tools for accurate studies of space-time profiles. We apply it to the description of energy transport between heat baths, and provide a full description of the current-carrying non-equilibrium steady state and the transition regions in a family of models including the Lieb-Liniger model of interacting Bose gases, realized in experiments.

We study the non-equilibrium dynamics obtained by an abrupt change (a {\em quench}) in the parameters of an integrable classical field theory, the nonlinear Schr\"odinger equation. We first consider explicit one-soliton examples, which we fully describe by solving the direct part of the inverse scattering problem. We then develop some aspects of the general theory using elements of the inverse scattering method. For this purpose, we introduce the {\em quench map} which acts on the space of scattering data and represents the change of parameter with fixed field configuration (initial condition). We describe some of its analytic properties by implementing a higher level version of the inverse scattering method, and we discuss the applications of Darboux-B\"acklund transformations, Gelfand-Levitan-Marchenko equations and the Rosales series solution to a related, dual quench problem. Finally, we comment on the interplay between quantum and classical tools around the theme of quenches and on the usefulness of the quantization of our classical approach to the quantum quench problem.

We re-examine the emergence of a universal non-equilibrium steady state following a local quench between quantum critical heat baths in spatial dimensions greater than one. We show that energy transport proceeds by the formation of an instantaneous shock wave and a broadening rarefaction wave on either side of the interface, and not by two shock waves as previously proposed. For small temperature differences the universal steady state energy currents of the two-shock and rarefaction-shock solutions coincide. Over a broad range of parameters, the difference in the energy flow across the interface between these two solutions is at the level of two percent. The properties of the energy flow remain fully universal and independent of the microscopic theory. We briefly discuss the width of the shock wave in a viscous fluid, the effects of momentum relaxation, and the generalization to charged fluids.

In this paper we propose a Monte Carlo method for generating finite-domain marginals of critical distributions of statistical models in infinite volume. The algorithm corrects the problem of the long-range effects of boundaries associated to generating critical distributions on finite lattices. It uses the advantage of scale invariance combined with ideas of the renormalization group in order to construct a type of "holographic" boundary condition that encodes the presence of an infinite volume beyond it. We check the quality of the distribution obtained in the case of the planar Ising model by comparing various observables with their infinite-plane prediction. We accurately reproduce planar two-, three- and four-point functions of spin and energy operators. We also define a lattice stress-energy tensor, and numerically obtain the associated conformal Ward identities and the Ising central charge.

We provide a pedagogical review of the main ideas and results in non-equilibrium conformal field theory and connected subjects. These concern the understanding of quantum transport and its statistics at and near critical points. Starting with phenomenological considerations, we explain the general framework, illustrated by the example of the Heisenberg quantum chain. We then introduce the main concepts underlying conformal field theory (CFT), the emergence of critical ballistic transport, and the CFT scattering construction of non-equilibrium steady states. Using this we review the theory for energy transport in homogeneous one-dimensional critical systems, including the complete description of its large deviations and the resulting (extended) fluctuation relations. We generalize some of these ideas to one-dimensional critical charge transport and to the presence of defects, as well as beyond one-dimensional criticality. We describe non-equilibrium transport in free-particle models, where connections are made with generalized Gibbs ensembles, and in higher-dimensional and non-integrable quantum field theories, where the use of the powerful hydrodynamic ideas for non-equilibrium steady states is explained. We finish with a list of open questions. The review does not assume any advanced prior knowledge of conformal field theory, large-deviation theory or hydrodynamics.

We consider the logarithmic negativity, a measure of bipartite entanglement, in a general unitary 1+1-dimensional massive quantum field theory, not necessarily integrable. We compute the negativity between a finite region of length $r$ and an adjacent semi-infinite region, and that between two semi-infinite regions separated by a distance $r$. We show that the former saturates to a finite value, and that the latter tends to zero, as $r\rightarrow\infty$. We show that in both cases, the leading corrections are exponential decays in $r$ (described by modified Bessel functions) that are solely controlled by the mass spectrum of the model, independently of its scattering matrix. This implies that, like the entanglement entropy, the logarithmic negativity displays a very high level of universality, allowing one to extract information about the mass spectrum. Further, a study of sub-leading terms shows that, unlike the entanglement entropy, a large-$r$ analysis of the negativity allows for the detection of bound states.

We develop a hydrodynamic approach to non-equilibrium conformal field theory. We study non-equilibrium steady states in the context of one-dimensional conformal field theory perturbed by the $T\bar T$ irrelevant operator. By direct quantum computation, we show, to first order in the coupling, that a relativistic hydrodynamic emerges, which is a simple modification of one-dimensional conformal fluids. We show that it describes the steady state and its approach, and we provide the main characteristics of the steady state, which lies between two shock waves. The velocities of these shocks are modified by the perturbation and equal the sound velocities of the asymptotic baths. Pushing further this approach, we are led to conjecture that the approach to the steady state is generically controlled by the power law $t^{-1/2}$, and that the widths of the shocks increase with time according to $t^{1/3}$.

In this paper we study the simplest massive 1+1 dimensional integrable quantum field theory which can be described as a perturbation of a non-unitary minimal conformal field theory: the Lee-Yang model. We are particularly interested in the features of the bi-partite entanglement entropy for this model and on building blocks thereof, namely twist field form factors. Non-unitarity selects out a new type of twist field as the operator whose two-point function (appropriately normalized) yields the entanglement entropy. We compute this two-point function both from a form factor expansion and by means of perturbed conformal field theory. We find good agreement with CFT predictions put forward in a recent work involving the present authors. In particular, our results are consistent with a scaling of the entanglement entropy given by $\frac{c_{\text{eff}}}{3}\log \ell$ where $c_{\text{eff}}$ is the effective central charge of the theory (a positive number related to the central charge) and $\ell$ is the size of the region. Furthermore the form factor expansion of twist fields allows us to explore the large region limit of the entanglement entropy and find the next-to-leading order correction to saturation. We find that this correction is very different from its counterpart in unitary models. Whereas in the latter case, it had a form depending only on few parameters of the model (the particle spectrum), it appears to be much more model-dependent for non-unitary models.

We study the dynamics of the entanglement in one dimensional critical quantum systems after a local quench in which two independently thermalized semi-infinite halves are joined to form a homogeneous infinite system and left to evolve unitarily. We show that under certain conditions a nonequilibrium steady state (NESS) is reached instantaneously as soon as the entanglement interval is within the light-cone emanating from the contact point. In this steady state, the exact expressions for the entanglement entropy and the logarithmic negativity are in agreement with the steady state density matrix being a boosted thermal state, as expected. We derive various general identities: relating the negativity after the quench with unequal left and right initial temperatures, with that with equal left and right temperatures; and relating these with the negativity in equilibrium thermal states. In certain regimes the resulting expressions can be analytically evaluated. Immediately after the interval interesects the lightcone, we find logarithmic growth. For a very long interval, we find that the negativity approaches a plateau after sufficiently long times, different from its NESS value. This provides a theoretical framework explaining recently obtained numerical results.

In this letter we show that the R\'enyi entanglement entropy of a region of large size $\ell$ in a one-dimensional critical model whose ground state breaks conformal invariance (such as in those described by non-unitary conformal field theories), behaves as $S_n \sim \frac{c_{\mathrm{eff}}(n+1)}{6n} \log \ell$, where $c_{\mathrm{eff}}=c-24\Delta>0$ is the effective central charge, $c$ (which may be negative) is the central charge of the conformal field theory and $\Delta\neq 0$ is the lowest holomorphic conformal dimension in the theory. We also obtain results for models with boundaries, and with a large but finite correlation length, and we show that if the lowest conformal eigenspace is logarithmic ($L_0 = \Delta I + N$ with $N$ nilpotent), then there is an additional term proportional to $\log(\log \ell)$. These results generalize the well known expressions for unitary models. We provide a general proof, and report on numerical evidence for a non-unitary spin chain and an analytical computation using the corner transfer matrix method for a non-unitary lattice model. We use a new algebraic technique for studying the branching that arises within the replica approach, and find a new expression for the entanglement entropy in terms of correlation functions of twist fields for non-unitary models.

Let an infinite, homogeneous, many-body quantum system be unitarily evolved for a long time from a state where two halves are independently thermalized. One says that a non-equilibrium steady state emerges if there are nonzero steady currents in the central region. In particular, their presence is a signature of ballistic transport. We analyze the consequences of the current observable being a conserved density; near equilibrium this is known to give rise to linear wave propagation and a nonzero Drude peak. Using the Lieb-Robinson bound, we derive, under a certain regularity condition, a lower bound for the non-equilibrium steady-state current determined by equilibrium averages. This shows and quantifies the presence of ballistic transport far from equilibrium. The inequality suggests the definition of "nonlinear sound velocities", which specialize to the sound velocity near equilibrium in non-integrable models, and "generalized sound velocities", which encode generalized Gibbs thermalization in integrable models. These are bounded by the Lieb-Robinson velocity. The inequality also gives rise to a bound on the energy current noise in the case of pure energy transport. We show that the inequality is satisfied in many models where exact results are available, and that it is saturated at one-dimensional criticality.

We construct non-equilibrium steady states in the Klein-Gordon theory in arbitrary space dimension $d$ following a local quench. We consider the approach where two independently thermalized semi-infinite systems, with temperatures $T_{\rm L}$ and $T_{\rm R}$, are connected along a $d-1$-dimensional hypersurface. A current-carrying steady state, described by thermally distributed modes with temperatures $T_{\rm L}$ and $T_{\rm R}$ for left and right-moving modes, respectively, emerges at late times. The non-equilibrium density matrix is the exponential of a non-local conserved charge. We obtain exact results for the average energy current and the complete distribution of energy current fluctuations. The latter shows that the long-time energy transfer can be described by a continuum of independent Poisson processes, for which we provide the exact weights. We further describe the full time evolution of local observables following the quench. Averages of generic local observables, including the stress-energy tensor, approach the steady state with a power-law in time, where the exponent depends on the initial conditions at the connection hypersurface. We describe boundary conditions and special operators for which the steady state is reached instantaneously on the connection hypersurface. A semiclassical analysis of freely propagating modes yields the average energy current at large distances and late times. We conclude by comparing and contrasting our findings with results for interacting theories and provide an estimate for the timescale governing the crossover to hydrodynamics. As a modification of our Klein-Gordon analysis we also include exact results for free Dirac fermions.