
Using the theory of generalized hydrodynamics (GHD), we derive exact
Eulerscale dynamical twopoint correlation functions of conserved densities
and currents in inhomogeneous, nonstationary states of manybody integrable
systems with weak spacetime variations. This extends previous works to
inhomogeneous and nonstationary situations. Using GHD projection operators, we
further derive formulae for Eulerscale twopoint functions of arbitrary local
fields, purely from the data of their homogeneous onepoint functions. These
are new also in homogeneous generalized Gibbs ensembles. The technique is based
on combining a fluctuationdissipation 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 sinhGordon model and the LiebLiniger 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 LeclairMussardotype infinite
formfactor series in integrable quantum field theory, and exact Eulerscale
twopoint functions of exponential fields in the sinhGordon model and of
powers of the density field in the LiebLiniger model. We also analyze
correlations in the partitioning protocol, extract largetime asymptotics, and,
in free models, derive all Eulerscale $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 quasiparticle
formulation of GHD remains valid for radiative modes, even though these do not
display particlelike 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 finiteaverage
observables; this set is larger for GGEs with higher conserved charges. We
concentrate on the sinhGordon model, which only has radiative modes, and study
transport in the domainwall initial problem as well as Eulerscale
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 largescale theory for the dynamics of
manybody integrable systems. It consists of an infinite set of conservation
laws for quasiparticles 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
stateindependent quasiparticle velocities, in a space equipped with a family
of metrics, parametrized by the quasiparticles' 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 quasiparticles, 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 LiebLiniger $\delta$Bose gas. Our
formula contains only quantities which are obtainable from the thermodynamic
Bethe ansatz. The Drude weight is an infinitedimensional 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
byproduct we obtain the dynamical twopoint 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 spacetime
(external) symmetries, we study conical twist fields in twodimensional
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 branchpoint 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 branchpoint 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 resumming 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
YangMills theory.

We define a simple model of conformal field theory in random spacetime
environments, which we refer to as stochastic conformal field theory. This
model accounts for the effects of dilute random impurities in strongly
interacting critical manybody systems. On one hand, surprisingly, although
impurities are separated by macroscopic distances, we find that the
infinitetime 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 currentgenerating initial states with a
temperature imbalance, we show that the energy current and density satisfy
simple dissipative hydrodynamic equations. We describe the spacetime scales at
which nonequilibrium currents exist. We show that a lightcone 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 manybody interacting
systems with infinitely many conserved charges. In this letter, we show that it
supersedes the widely used conventional hydrodynamics (CHD) of onedimensional
Bose gases. We illustrate this by studying "nonlinear sound waves" emanating
from initial density accumulations in the LiebLiniger 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 finitedimensional 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 NRGTSAABACUS
algorithm. The analysis is performed in the full interaction range, and is not
restricted to either weak or strongrepulsion regimes.

The evaluation of vacuum expectation values (VEVs) in massive integrable
quantum field theory (QFT) is a nontrivial renormalizationgroup "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 branchpoint 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 formfactor
series which truncate at twoparticle level in free QFT. We verify the method
by studying U(1) twist fields in free models, which are simply related to the
branchpoint 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 formfactor observations by Bianchini
and CastroAlvaredo [Nucl. Phys. B913 (2016)], the standard powerlaw
shortdistance 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 nearcritical harmonic chains, including log log
corrections.

Inspired by recent results on the nonequilibrium dynamics of manybody
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
spacetime scale, the fluid of hard rods satisfies Eulertype 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 NavierStokes (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 nonunitary but
${\cal PT}$invariant quantum field theory in two spacetime dimensions. In
addition to unbroken $\mathcal{PT}$symmetry and a positive energy spectrum, we
assume standard properties of quantum field theory including a local
energymomentum 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 nonnegative 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}}=c24\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 quantumclassical equivalence: we demonstrate the equivalence
between the equations of GHD, and the Eulerscale hydrodynamic equations of a
new family of classical gases which generalize the gas of hard rods. In this
family, the "quasiparticles", upon colliding, jump forward or backward by a
distance that depends on their velocities, generalizing the jump forward by the
rods' length of the fixedvelocity tracer upon elastic collision of two hard
rods. Such velocitydependent 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 LiebLiniger model realized in coldatom experiments.

Generalized hydrodynamics (GHD) was proposed recently as a formulation of
hydrodynamics for integrable systems, taking into account infinitelymany
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 nonequilibrium
dynamics of onedimensional 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 Eulerlike 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 socalled 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
translationinvariant states on hypercubic quantum lattices of any
dimensionality (including quantum chains) and finiterange 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
manybody quantum systems, it is paramount to develop powerful methods that
encode the emergent physics. Up to now, the strong dichotomy observed between
integrable and nonintegrable evolutions made an overarching theory difficult
to build, especially for transport phenomena where spacetime profiles are
drastically different. We present a novel framework for studying transport in
integrable systems: hydrodynamics with infinitelymany conservation laws. This
bridges the conceptual gap between integrable and nonintegrable quantum
dynamics, and gives powerful tools for accurate studies of spacetime profiles.
We apply it to the description of energy transport between heat baths, and
provide a full description of the currentcarrying nonequilibrium steady state
and the transition regions in a family of models including the LiebLiniger
model of interacting Bose gases, realized in experiments.

We study the nonequilibrium 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 onesoliton
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
DarbouxB\"acklund transformations, GelfandLevitanMarchenko 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 reexamine the emergence of a universal nonequilibrium 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 twoshock and rarefactionshock 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 finitedomain
marginals of critical distributions of statistical models in infinite volume.
The algorithm corrects the problem of the longrange 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 infiniteplane prediction. We accurately
reproduce planar two, three and fourpoint functions of spin and energy
operators. We also define a lattice stressenergy 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
nonequilibrium 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 nonequilibrium steady states. Using this we review the theory for energy
transport in homogeneous onedimensional critical systems, including the
complete description of its large deviations and the resulting (extended)
fluctuation relations. We generalize some of these ideas to onedimensional
critical charge transport and to the presence of defects, as well as beyond
onedimensional criticality. We describe nonequilibrium transport in
freeparticle models, where connections are made with generalized Gibbs
ensembles, and in higherdimensional and nonintegrable quantum field theories,
where the use of the powerful hydrodynamic ideas for nonequilibrium 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,
largedeviation theory or hydrodynamics.

We consider the logarithmic negativity, a measure of bipartite entanglement,
in a general unitary 1+1dimensional massive quantum field theory, not
necessarily integrable. We compute the negativity between a finite region of
length $r$ and an adjacent semiinfinite region, and that between two
semiinfinite 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 subleading
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 nonequilibrium conformal field theory.
We study nonequilibrium steady states in the context of onedimensional
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
onedimensional 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 nonunitary
minimal conformal field theory: the LeeYang model. We are particularly
interested in the features of the bipartite entanglement entropy for this
model and on building blocks thereof, namely twist field form factors.
Nonunitarity selects out a new type of twist field as the operator whose
twopoint function (appropriately normalized) yields the entanglement entropy.
We compute this twopoint 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 nexttoleading 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
modeldependent for nonunitary models.

We study the dynamics of the entanglement in one dimensional critical quantum
systems after a local quench in which two independently thermalized
semiinfinite 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 lightcone 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 onedimensional critical model whose ground state breaks
conformal invariance (such as in those described by nonunitary conformal field
theories), behaves as $S_n \sim \frac{c_{\mathrm{eff}}(n+1)}{6n} \log \ell$,
where $c_{\mathrm{eff}}=c24\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 nonunitary spin chain and an analytical
computation using the corner transfer matrix method for a nonunitary 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 nonunitary
models.

Let an infinite, homogeneous, manybody quantum system be unitarily evolved
for a long time from a state where two halves are independently thermalized.
One says that a nonequilibrium 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
LiebRobinson bound, we derive, under a certain regularity condition, a lower
bound for the nonequilibrium steadystate 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
nonintegrable models, and "generalized sound velocities", which encode
generalized Gibbs thermalization in integrable models. These are bounded by the
LiebRobinson 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 onedimensional criticality.

We construct nonequilibrium steady states in the KleinGordon theory in
arbitrary space dimension $d$ following a local quench. We consider the
approach where two independently thermalized semiinfinite systems, with
temperatures $T_{\rm L}$ and $T_{\rm R}$, are connected along a
$d1$dimensional hypersurface. A currentcarrying steady state, described by
thermally distributed modes with temperatures $T_{\rm L}$ and $T_{\rm R}$ for
left and rightmoving modes, respectively, emerges at late times. The
nonequilibrium density matrix is the exponential of a nonlocal conserved
charge. We obtain exact results for the average energy current and the complete
distribution of energy current fluctuations. The latter shows that the
longtime 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 stressenergy tensor, approach the
steady state with a powerlaw 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 KleinGordon
analysis we also include exact results for free Dirac fermions.