Arxiv:2007.04861V2 [Cond-Mat.Quant-Gas] 23 Jul 2020
Total Page:16
File Type:pdf, Size:1020Kb
Thermalisation of a trapped one-dimensional Bose gas via diffusion Alvise Bastianello,1 Andrea De Luca,2 Benjamin Doyon,3 and Jacopo De Nardis4 1Institute for Theoretical Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands 2Laboratoire de Physique Th´eorique et Mod´elisation(CNRS UMR 8089), Universit´ede Cergy-Pontoise, F-95302 Cergy-Pontoise, France 3Department of Mathematics, King's College London, Strand WC2R 2LS, London, U.K. 4Department of Physics and Astronomy, University of Ghent, Krijgslaan 281, 9000 Gent, Belgium For a decade the fate of a one-dimensional gas of interacting bosons in an external trapping potential remained mysterious. We here show that whenever the underlying integrability of the gas is broken by the presence of the external potential, the inevitable diffusive rearrangements between the quasiparticles, quantified by the diffusion constants of the gas, eventually lead the system to thermalise at late times. We show that the full thermalising dynamics can be described by the generalised hydrodynamics with diffusion and force terms, and we compare these predictions with numerical simulations. Finally, we provide an explanation for the slow thermalisation rates observed in numerical and experimental settings: the hydrodynamics of integrable models is characterised by a continuity of modes, which can have arbitrarily small diffusion coefficients. As a consequence, the approach to thermalisation can display pre-thermal plateau and relaxation dynamics with long polynomial finite-time corrections. In the past two decades, low-temperature gases of of the external force that breaks integrability and the dif- bosonic or fermionic atoms emerged as the best exper- fusion spreading that redistributes the quasiparticle mo- imental platform where to study many-body phenomena menta and increases the local thermodynamic entropy. and where to probe their non-equilibrium dynamics [1{3]. Previous works have studied the hydrodynamics of an A great deal of attention has been dedicated to the prob- integrable gas in the presence of the external force but lem of thermalisation in isolated systems, such as a single without accounting for diffusive terms [20{24]. Here, we gas of bosonic atoms interacting via contact repulsion. show that diffusive corrections play an essential role in The experiment of 2007 dubbed \quantum Newton cra- the late-time dynamics. dle" [4] represented a turning point. The apparent lack Lieb-Liniger gas in an external potential. Our of thermalisation at late times for a non-equilibrium one- model of interest is the Lieb-Liniger model for a gas of dimensional cold atomic gas in an external trap started bosons interacting with contact repulsion. The Hamilto- numerous research directions on the role of the under- nian is given by lying integrability in the isolated dynamics. It was then 2 later clarified that integrable models do not thermalise to ~ y y y H^LL = dx @x (x)@x (x)+g (x) (x) (x) (x); standard Gibbs ensembles, but to Generalises Gibbs ones ˆ 2m [1,5,6], where entropy is maximised given the constraints (1) y of all local and quasi-local conserved quantities. As the with corresponding bosonic operators [ (x); (y)] = integrability is weakly broken by the external trap, one δ(x − y). In the following we set the mass of the bosons 2 expects therefore to crossover from a fast relaxation to a such that ~ =2m = 1 and g = c. The model is known to generalised Gibbs state, to a slow pre-thermal decay to describe state-of-the-art cold atom experiments [4, 25{ the final thermal ensemble. In the past years, the study of 32] and also to be integrable [33]. As a consequence, all pre-thermalisation dynamics has received many contribu- physical initial states [34{39] eventually relax to a GGE − P Q^ β tions, and for homogeneous systems, a quite comprehen- density matrix% ^ ∼ e j j j , where the set of Lagrange sive understanding has now been achieved [7{17]. How- multipliers βj is fixed by the initial condition. ever, a full description of the (pre) thermalisation mech- Its large scale dynamics can be accurately described anism in inhomogenous systems is still lacking. While via the so-called Generalised Hydrodynamics (GHD), in- a recent experimental work showed that multiple one- troduced in the past years [40, 41] and the resulting arXiv:2007.04861v2 [cond-mat.quant-gas] 23 Jul 2020 dimensional gases coupled to each other or where tran- predictions have even received an experimental confir- sitions to higher dimension are allowed, do indeed ther- mation [31]. In contrast with standard hydrodynam- malise at late times [18, 19], it remains to be understood ics, GHD is built on the infinite number of local and ^ what is the fate of a single isolated one-dimensional gas. quasi-local conserved quantities Qj = dx q^j(x) charac- terising integrable systems. Equivalently,´ each stationary In this letter, we show that in the presence of an ex- state of the system can be described by a set of stable ternal force and of diffusive spreading, a single gas de- quasiparticles with density ρθ, in one-to-one correspon- scribed by an integrable Hamiltonian and by an external dence with the expectation values of conserved quantities trapping potential does thermalise, due to the interplay dθρθ hj,θ = hq^ji, with hj,θ, the single-particle eigenval- ´ 2 interacting systems, describes both a diffusive spreading of the quasiparticle motion and a re-distribution of quasi- particle velocities. Both quantities are non-linear func- tional of the densities ρθ(x; t), fixed by the equation of state and by the explicit form of the diffusion kernel. The effect of an external trap is described by the fol- lowing Hamiltonian H^ = H^ + dx V^ (x) (2) trap LL ˆ where the potential V^ (x) = V (x)^q0(x) and the parti- y cle densityq ^0(x) = (x) (x). The contribution of the trapping potential to the dynamics of the quasiparticle densities can be treated perturbatively in the degree of smoothness of the external potential V (x). It was first shown in [20] that at the Euler scale this amounts to adding the simple force term (@xV (x))@θρθ. We here FIG. 1. Left: Quasiparticle density ρθ(x; t) at t = L claim that the addition of this single term also gives the for a Lieb-Liniger gas (2) inside the trap V (x) = 1=2(1 − correct hydrodynamic equation up to diffusive scales: cos 2πx=L), in a circle with length L = 1 (top) and L = 25:6 eff α (down) (see text for other parameter's values). Right: Log- @tρθ = −@x(vθ ρθ) + @x(Dθ @xρα) + (@xV (x))@θρθ: (3) plot of the thermodynamic entropy density S(t)=L (top) con- verging to the thermal entropy for the given quench and of That is, this equation holds up to order O(L−3) cor- 4 L y 2 the integrated correlator hhj jii = 0 dxh( ) i=L. The rections, with L the minimal length of spatial varia- inset is a log-plot to show the approach´ to the thermal value tions of the system and of the external potential V (x). at late times. From a scaling analysis, at this diffusive order possible corrections to Eq. (3) could come from three different 2 sources: i) second derivative in the potential @xV (x); ii) second power of its first derivative (@ V (x))2; and ues of each quasiparticle with rapidity θ. x iii) mixed derivatives, @xV (x)@xρα. We now argue that Given the derivative of the scattering shift between two none of these terms contributes. First, we recall that quasiparticles with different rapidities T , any single- θ,α the definitions of densities of chargesq ^j(x) are ambigu- particle function hj,θ is dressed via the linear transfor- ous by the addition of total derivative terms, and as dr −1 0 maton hj = (1 − T n) · h, with nθ = 2πρθ=(kθ) the shown in [49], we can make use of this freedom to en- fermionic filling function and (1−T n)θ,α = δθ,α −Tθ,αnα force PT (parity and time) symmetry. With PT sym- the dressing matrix (here and in the following we de- metry, we can argue that terms proportional to @2V 0 x note the derivative respect to rapidities as @θhθ = hθ). are excluded in (3). Specifically, consider the Heisen- The effective quasiparticle momentum kθ and energy berg time-evolution equation for any conserved density " are obtained by integrating the dressed derivatives: θ q^j(0), say at the origin. Taylor expanding the potential, 0 −1 bare 0 k = (1 − T n) · (k ) (and similarly for the en- V (x) = V (0) + xV 0(0) + O(x2), we argue that the O(x2) ergy). For the LL model, T = 2c=((θ − α)2 + c2), and 2 θ,α does not contribute to Eq. (3), as dx x [^q0(x); q^j(0)] bare 0 bare 0 (kθ ) = 1, ("θ ) = 2θ. Because of the interactions, has a vanishing expectation on any´ stationary state by the effective energy and momentum of quasiparticles de- PT symmetry. Second, no mixed derivatives appear as pend non-trivially on the local thermodynamic state via the charge densityq ^0(x) is ultra-local. That is, in terms the dressing of excitations. of the field (x) it does not contain any derivative, and Without external forces, the large scale evolution of as a consequence the number of particles does not gen- any local density of conserved quantities can be writ- erate any flow in space [ dxq^0(x); qi] = 0. This implies ten in terms of the evolution of the local density of that diffusion is not modified´ by the presence of the ex- quasiparticles ρθ(x; t).