Publications

Order by: publication date number of citations

• A note on generalized hydrodynamics: inhomogeneous fields and other concepts

  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.

  156 citations

• Diffusion in generalized hydrodynamics and quasiparticle scattering

  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.

  153 citations

• Lecture notes on Generalised Hydrodynamics

  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.

  152 citations

• Drude Weight for the Lieb-Liniger Bose Gas

  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.

  138 citations

• Walking, Weak first-order transitions, and Complex CFTs II. Two-dimensional Potts model at $Q>4$

  Victor Gorbenko, Slava Rychkov, Bernardo Zan

  SciPost Phys. 5, 050 (2018) · published 22 November 2018 | Toggle abstract · pdf

  We study complex CFTs describing fixed points of the two-dimensional $Q$-state Potts model with $Q>4$. Their existence is closely related to the weak first-order phase transition and walking RG behavior present in the real Potts model at $Q>4$. The Potts model, apart from its own significance, serves as an ideal playground for testing this very general relation. Cluster formulation provides nonperturbative definition for a continuous range of parameter $Q$, while Coulomb gas description and connection to minimal models provide some conformal data of the complex CFTs. We use one and two-loop conformal perturbation theory around complex CFTs to compute various properties of the real walking RG flow. These properties, such as drifting scaling dimensions, appear to be common features of the QFTs with walking RG flows, and can serve as a smoking gun for detecting walking in Monte Carlo simulations. The complex CFTs discussed in this work are perfectly well defined, and can in principle be seen in Monte Carlo simulations with complexified coupling constants. In particular, we predict a pair of $S_5$-symmetric complex CFTs with central charges $c\approx 1.138 \pm 0.021 i$ describing the fixed points of a 5-state dilute Potts model with complexified temperature and vacancy fugacity.

  91 citations

• Exact large-scale correlations in integrable systems out of equilibrium

  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.

  89 citations

• Universality of the SAT-UNSAT (jamming) threshold in non-convex continuous constraint satisfaction problems

  Silvio Franz, Giorgio Parisi, Maksim Sevelev, Pierfrancesco Urbani, Francesco Zamponi

  SciPost Phys. 2, 019 (2017) · published 2 June 2017 | Toggle abstract · pdf

  Random constraint satisfaction problems (CSP) have been studied extensively using statistical physics techniques. They provide a benchmark to study average case scenarios instead of the worst case one. The interplay between statistical physics of disordered systems and computer science has brought new light into the realm of computational complexity theory, by introducing the notion of clustering of solutions, related to replica symmetry breaking. However, the class of problems in which clustering has been studied often involve discrete degrees of freedom: standard random CSPs are random K-SAT (aka disordered Ising models) or random coloring problems (aka disordered Potts models). In this work we consider instead problems that involve continuous degrees of freedom. The simplest prototype of these problems is the perceptron. Here we discuss in detail the full phase diagram of the model. In the regions of parameter space where the problem is non-convex, leading to multiple disconnected clusters of solutions, the solution is critical at the SAT/UNSAT threshold and lies in the same universality class of the jamming transition of soft spheres. We show how the critical behavior at the satisfiability threshold emerges, and we compute the critical exponents associated to the approach to the transition from both the SAT and UNSAT phase. We conjecture that there is a large universality class of non-convex continuous CSPs whose SAT-UNSAT threshold is described by the same scaling solution.

  75 citations

• RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles

  Nicholas P. Bailey, Trond S. Ingebrigtsen, Jesper Schmidt Hansen, Arno A. Veldhorst, Lasse Bøhling, Claire A. Lemarchand, Andreas E. Olsen, Andreas K. Bacher, Lorenzo Costigliola, Ulf R. Pedersen, Heine Larsen, Jeppe C. Dyre, Thomas B. Schrøder

  SciPost Phys. 3, 038 (2017) · published 14 December 2017 | Toggle abstract · pdf

  RUMD is a general purpose, high-performance molecular dynamics (MD) simulation package running on graphical processing units (GPU's). RUMD addresses the challenge of utilizing the many-core nature of modern GPU hardware when simulating small to medium system sizes (roughly from a few thousand up to hundred thousand particles). It has a performance that is comparable to other GPU-MD codes at large system sizes and substantially better at smaller sizes.RUMD is open-source and consists of a library written in C++ and the CUDA extension to C, an easy-to-use Python interface, and a set of tools for set-up and post-simulation data analysis. The paper describes RUMD's main features, optimizations and performance benchmarks.

  69 citations

• Twisted and untwisted negativity spectrum of free fermions

  Hassan Shapourian, Paola Ruggiero, Shinsei Ryu, Pasquale Calabrese

  SciPost Phys. 7, 037 (2019) · published 25 September 2019 | Toggle abstract · pdf

  A basic diagnostic of entanglement in mixed quantum states is known as the positive partial transpose (PT) criterion. Such criterion is based on the observation that the spectrum of the partially transposed density matrix of an entangled state contains negative eigenvalues, in turn, used to define an entanglement measure called the logarithmic negativity. Despite the great success of logarithmic negativity in characterizing bosonic many-body systems, generalizing the operation of PT to fermionic systems remained a technical challenge until recently when a more natural definition of PT for fermions that accounts for the Fermi statistics has been put forward. In this paper, we study the many-body spectrum of the reduced density matrix of two adjacent intervals for one-dimensional free fermions after applying the fermionic PT. We show that in general there is a freedom in the definition of such operation which leads to two different definitions of PT: the resulting density matrix is Hermitian in one case, while it becomes pseudo-Hermitian in the other case. Using the path-integral formalism, we analytically compute the leading order term of the moments in both cases and derive the distribution of the corresponding eigenvalues over the complex plane. We further verify our analytical findings by checking them against numerical lattice calculations.

  60 citations

• Nonlinear modes disentangle glassy and Goldstone modes in structural glasses

  Luka Gartner, Edan Lerner

  SciPost Phys. 1, 016 (2016) · published 31 December 2016 | Toggle abstract · pdf

  One outstanding problem in the physics of glassy solids is understanding the statistics and properties of the low-energy excitations that stem from the disorder that characterizes these systems' microstructure. In this work we introduce a family of algebraic equations whose solutions represent collective displacement directions (modes) in the multi-dimensional configuration space of a structural glass. We explain why solutions of the algebraic equations, coined nonlinear glassy modes, are quasi-localized low-energy excitations. We present an iterative method to solve the algebraic equations, and use it to study the energetic and structural properties of a selected subset of their solutions constructed by starting from a normal mode analysis of the potential energy of a model glass. Our key result is that the structure and energies associated with harmonic glassy vibrational modes and their nonlinear counterparts converge in the limit of very low frequencies. As nonlinear modes never suffer hybridizations, our result implies that the presented theoretical framework constitutes a robust alternative definition of `soft glassy modes' in the thermodynamic limit, in which Goldstone modes overwhelm and destroy the identity of low-frequency harmonic glassy modes.

  57 citations