The formation of dispersive shock waves in one-dimensional Bose gas represents a limitation of Generalized Hydrodynamics (GHD) due to the coarse-grained nature of the theory. Nevertheless, GHD accurately captures long-wavelength behavior, thus indicating an implicit knowledge of the underlying microscopic physics. Such representations are already known through the Whitham modulation theory, where dispersionless equations describe the evolution of the slowly varying shock wave parameters. Here we study the correspondence between Whitham's approach to the Gross-Pitaevskii equation and GHD in the semiclassical limit and beyond. Our findings enable the recovery of the shock wave solution directly from GHD simulations, which we demonstrate for both zero and finite temperature. Additionally, we study how free expansion protocols affect the shock wave density and their implications for experimental detection. The combined picture of Whitham and GHD lends itself to additional physical interpretation regarding the formation of shock waves. Published by the American Physical Society 2024
The formation of dispersive shock waves in one-dimensional Bose gas represents a limitation of Generalized Hydrodynamics (GHD) due to the coarse-grained nature of the theory. Nevertheless, GHD accurately captures long-wavelength behavior, thus indicating an implicit knowledge of the underlying microscopic physics. Such representations are already known through the Whitham modulation theory, where dispersionless equations describe the evolution of the slowly varying shock wave parameters. Here we study the correspondence between Whitham's approach to the Gross-Pitaevskii equation and GHD in the semiclassical limit and beyond. Our findings enable the recovery of the shock wave solution directly from GHD simulations, which we demonstrate for both zero and finite temperature. Additionally, we study how free expansion protocols affect the shock wave density and their implications for experimental detection. The combined picture of Whitham and GHD lends itself to additional physical interpretation regarding the formation of shock waves. Published by the American Physical Society 2024
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. It offers a full and efficient framework to capture the long-wavelengths nonequilibrium dynamics in experimentally realizable one-dimensional systems, like Bose gases with delta repulsion. In this letter, we study the nonlinear sound waves emanating from initial density accumulations in the zero-temperature Bose gas. We numerically confirm the emergence of GHD by comparing its predictions with a full quantum simulation performed using the NRG-TSA-ABACUS algorithm. We prove that at zero temperature, GHD is equivalent to conventional Galilean hydrodynamics with nontrivial equations of state. We further predict, at long wavelengths and large densities, the appearance of sharp profiles in finite times, that immediately dissolve (no sustained shocks). The analysis is performed in the full interaction range, and is not restricted to either weak- or strong-repulsion regimes.
Abstract The problem of wave breaking during its propagation in the Bose–Einstein condensate to a stationary medium is considered for the case when the initial profile at the breaking instant can be approximated by a power function of the form (–x)1/n. The evolution of the wave is described by the Gross–Pitaevskii equation so that a dispersive shock wave is formed as a result of breaking; this wave can be represented using the Gurevich–Pitaevskii approach as a modulated periodic solution to the Gross–Pitaevskii equation, and the evolution of the modulation parameters is described by the Whitham equations obtained by averaging the conservation laws over fast oscillations in the wave. The solution to the Whitham modulation equations is obtained in closed form for n = 2, 3, and the velocities of the dispersion shock wave edges for asymptotically long evolution times are determined for arbitrary integers n > 1. The problem considered here can be applied for describing the generation of dispersion shock waves observed in experiments with the Bose–Einstein condensate.
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 results of benchmarks of the recent theory of generalized hydrodynamics (GHD) in the context of the dynamics of the 1D Bose gas described by the Lieb-Liniger model are reported. Through comparisons against several alternative, well-established computational approaches, the range of applicability of GHD in some of the most challenging out-of-equilibrium scenarios, such as propagation of quantum shock waves and collisional relaxation, is verified and validated and, additionally, some of the methods limitations are pointed out.
We study the formation and the subsequent dynamics of shock waves in repulsive one-dimensional Bose gases during the free expansion of a density hump. By building coherent Fermi states for interacting Bethe fermions, we define a quantum fluctuating initial state expressed in terms of universal quantities, namely the density and the Luttinger parameter. In the integrable case, this fluctuating state is then evolved by generalized hydrodynamics (GHD) and, differently from nonfluctuating initial states, it develops density ripples on top of the hydrodynamic mean value. Our analysis gives a general theory of quantum ripples and wave breaking in integrable and quasi-integrable one-dimensional liquids and clarifies the role of the interaction strength. In particular, for strongly/intermediately interacting bosons, we find quantum ripples originating from low-energy modes at the Fermi surface interfering when transported by GHD. In the low-coupling limit, near the quasicondensate regime, we find instead that density ripples have a semiclassical nature, and their description requires information on the curvature of the Fermi surface.
Generalized hydrodynamics (GHD) is a recent theoretical approach that is becoming a go-to tool for characterizing out-of-equilibrium phenomena in integrable and near-integrable quantum many-body systems. Here, we benchmark its performance against an array of alternative theoretical methods, for an interacting one-dimensional Bose gas described by the Lieb-Liniger model. In particular, we study the evolution of both a localized density bump and dip, along with a quantum Newton's cradle setup, for various interaction strengths and initial equilibrium temperatures. We find that GHD generally performs very well at sufficiently high temperatures or strong interactions. For low temperatures and weak interactions, we highlight situations where GHD, while not capturing interference phenomena on short lengthscales, can describe a coarse-grained behaviour based on convolution averaging that mimics finite imaging resolution in ultracold atom experiments. In a quantum Newton's cradle setup based on a double-well to single-well trap quench, we find that GHD with diffusive corrections demonstrates excellent agreement with the predictions of a classical field approach.