Contributor info: Prof. Benjamin Doyon

Giuseppe del Vecchio del Vecchio, Benjamin Doyon, Paola Ruggiero

SciPost Phys. Core 7, 005 (2024) · published 14 February 2024 | Toggle abstract · pdf

The large-scale behaviour of entanglement entropy in finite-density states, in and out of equilibrium, can be understood using the physical picture of particle pairs. However, the full theoretical origin of this picture is not fully established yet. In this work, we clarify this picture by investigating entanglement entropy using its connection with the large-deviation theory for thermodynamic and hydrodynamic fluctuations. We apply the universal framework of Ballistic Fluctuation Theory (BFT), based the Euler hydrodynamics of the model, to correlation functions of branch-point twist fields, the starting point for computing Rényi entanglement entropies within the replica approach. Focusing on free fermionic systems in order to illustrate the ideas, we show that both the equilibrium behavior and the dynamics of Rényi entanglement entropies can be fully derived from the BFT. In particular, we emphasise that long-range correlations develop after quantum quenches, and accounting for these explain the structure of the entanglement growth. We further show that this growth is related to fluctuations of charge transport, generalising to quantum quenches the relation between charge fluctuations and entanglement observed earlier. The general ideas we introduce suggest that the large-scale behaviour of entanglement has its origin within hydrodynamic fluctuations.

Benjamin Doyon, Gabriele Perfetto, Tomohiro Sasamoto, Takato Yoshimura

SciPost Phys. 15, 136 (2023) · published 4 October 2023 | Toggle abstract · pdf

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.

Benjamin Doyon, Joseph Durnin, Takato Yoshimura

SciPost Phys. 13, 072 (2022) · published 28 September 2022 | Toggle abstract · pdf

We introduce an extension of the generalised $T\bar{T}$-deformation described by Smirnov-Zamolodchikov, to include the complete set of extensive charges. We show that this gives deformations of S-matrices beyond CDD factors, generating arbitrary functional dependence on momenta. We further derive from basic principles of statistical mechanics the flow equations for the free energy and all free energy fluxes. From this follows, without invoking the microscopic Bethe ansatz or other methods from integrability, that the thermodynamics of the deformed models are described by the integral equations of the thermodynamic Bethe-Ansatz, and that the exact average currents take the form expected from generalised hydrodynamics, both in the classical and quantum realms.

Olalla A. Castro-Alvaredo, Cecilia De Fazio, Benjamin Doyon, Aleksandra A. Ziółkowska

SciPost Phys. 12, 115 (2022) · published 31 March 2022 | Toggle abstract · pdf

In the context of quantum field theory (QFT), unstable particles are associated with complex-valued poles of two-body scattering matrices in the unphysical sheet of rapidity space. The Breit-Wigner formula relates this pole to the mass and life-time of the particle, observed in scattering events. In this paper, we uncover new, dynamical signatures of unstable excitations and show that they have a strong effect on the non-equilibrium properties of QFT. Focusing on a 1+1D integrable model, and using the theory of Generalized Hydrodynamics, we study the formation and decay of unstable particles by analysing the release of hot matter into a low-temperature environment. We observe the formation of tails and the decay of the emitted nonlinear waves, in sharp contrast to the situation without unstable excitations. We also uncover a new phenomenon by which a wave of a stable population of unstable particles may persist without decay for long times. We expect these signatures of the presence of unstable particles to have a large degree of universality. Our study shows that the out-of-equilibrium dynamics of many-body systems can be strongly affected not only by the spectrum, but also by excitations with finite life-times.

Gabriele Perfetto, Benjamin Doyon

SciPost Phys. 10, 116 (2021) · published 27 May 2021 | Toggle abstract · pdf

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.

Frederik S. Møller, Gabriele Perfetto, Benjamin Doyon, Jörg Schmiedmayer

SciPost Phys. Core 3, 016 (2020) · published 22 December 2020 | Toggle abstract · pdf

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.

Isabelle Bouchoule, Benjamin Doyon, Jerome Dubail

SciPost Phys. 9, 044 (2020) · published 2 October 2020 | Toggle abstract · pdf

We theoretically investigate the effects of atom losses in the one-dimensional (1D) Bose gas with repulsive contact interactions, a famous quantum integrable system also known as the Lieb-Liniger gas. The generic case of K-body losses (K = 1,2,3,...) is considered. We assume that the loss rate is much smaller than the rate of intrinsic relaxation of the system, so that at any time the state of the system is captured by its rapidity distribution (or, equivalently, by a Generalized Gibbs Ensemble). We give the equation governing the time evolution of the rapidity distribution and we propose a general numerical procedure to solve it. In the asymptotic regimes of vanishing repulsion -- where the gas behaves like an ideal Bose gas -- and hard-core repulsion -- where the gas is mapped to a non-interacting Fermi gas -- we derive analytic formulas. In the latter case, our analytic result shows that losses affect the rapidity distribution in a non-trivial way, the time derivative of the rapidity distribution being both non-linear and non-local in rapidity space.

Benjamin Doyon

SciPost Phys. Lect. Notes 18 (2020) · published 28 August 2020 | Toggle abstract · pdf

These are lecture notes for a series of lectures given at the Les Houches Summer School on Integrability in Atomic and Condensed Matter Physics, 30 July to 24 August 2018. The same series of lectures has also been given at the Tokyo Institute of Technology, October 2019. I overview in a pedagogical fashion the main aspects of the theory of generalised hydrodynamics, a hydrodynamic theory for quantum and classical many-body integrable systems. Only very basic knowledge of hydrodynamics and integrable systems is assumed.

Jason Myers, M. J. Bhaseen, Rosemary J. Harris, Benjamin Doyon

SciPost Phys. 8, 007 (2020) · published 20 January 2020 | Toggle abstract · pdf

We propose exact results for the full counting statistics, or the scaled cumulant generating function, pertaining to the transfer of arbitrary conserved quantities across an interface in homogeneous integrable models out of equilibrium. We do this by combining insights from generalised hydrodynamics with a theory of large deviations in ballistic transport. The results are applicable to a wide variety of physical systems, including the Lieb-Liniger gas and the Heisenberg chain. We confirm the predictions in non-equilibrium steady states obtained by the partitioning protocol, by comparing with Monte Carlo simulations of this protocol in the classical hard rod gas. We verify numerically that the exact results obey the correct non-equilibrium fluctuation relations with the appropriate initial conditions.

Jean-Sébastien Caux, Benjamin Doyon, Jérôme Dubail, Robert Konik, Takato Yoshimura

SciPost Phys. 6, 070 (2019) · published 20 June 2019 | Toggle abstract · pdf

Describing and understanding the motion of quantum gases out of equilibrium is one of the most important modern challenges for theorists. In the groundbreaking Quantum Newton Cradle experiment [Kinoshita, Wenger and Weiss, Nature 440, 900, 2006], quasi-one-dimensional cold atom gases were observed with unprecedented accuracy, providing impetus for many developments on the effects of low dimensionality in out-of-equilibrium physics. But it is only recently that the theory of generalized hydrodynamics has provided the adequate tools for a numerically efficient description. Using it, we give a complete numerical study of the time evolution of an ultracold atomic gas in this setup, in an interacting parameter regime close to that of the original experiment. We evaluate the full evolving phase-space distribution of particles. We simulate oscillations due to the harmonic trap, the collision of clouds without thermalization, and observe a small elongation of the actual oscillation period and cloud deformations due to many-body dephasing. We also analyze the effects of weak anharmonicity. In the experiment, measurements are made after release from the one-dimensional trap. We evaluate the gas density curves after such a release, characterizing the actual time necessary for reaching the asymptotic state where the integrable quasi-particle momentum distribution function emerges.

Jacopo De Nardis, Denis Bernard, Benjamin Doyon

SciPost Phys. 6, 049 (2019) · published 25 April 2019 | Toggle abstract · pdf

We extend beyond the Euler scales the hydrodynamic theory for quantum and classical integrable models developed in recent years, accounting for diffusive dynamics and local entropy production. We review how the diffusive scale can be reached via a gradient expansion of the expectation values of the conserved fields and how the coefficients of the expansion can be computed via integrated steady-state two-point correlation functions, emphasising that PT-symmetry can fully fix the inherent ambiguity in the definition of conserved fields at the diffusive scale. We develop a form factor expansion to compute such correlation functions and we show that, while the dynamics at the Euler scale is completely determined by the density of single quasiparticle excitations on top of the local steady state, diffusion is due to scattering processes among quasiparticles, which are only present in truly interacting systems. We then show that only two-quasiparticle scattering processes contribute to the diffusive dynamics. Finally we employ the theory to compute the exact spin diffusion constant of a gapped XXZ spin-1/2 chain at finite temperature and half-filling, where we show that spin transport is purely diffusive.

Benjamin Doyon

SciPost Phys. 5, 054 (2018) · published 28 November 2018 | Toggle abstract · pdf

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$-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.

Alvise Bastianello, Benjamin Doyon, Gerard Watts, Takato Yoshimura

SciPost Phys. 4, 045 (2018) · published 30 June 2018 | Toggle abstract · pdf

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.

Benjamin Doyon, Herbert Spohn

SciPost Phys. 3, 039 (2017) · published 15 December 2017 | Toggle abstract · pdf

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.

Benjamin Doyon, Takato Yoshimura

SciPost Phys. 2, 014 (2017) · published 22 April 2017 | Toggle abstract · pdf

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.