We report diffusion Monte Carlo (DMC) and path integral Monte Carlo (PIMC) calculations of the properties of a one-dimensional (1D) Bose quantum fluid. The equation of state, the superfluid fraction ${\ensuremath{\rho}}_{S}/{\ensuremath{\rho}}_{0}$, the one-body density matrix
n(x), the pair distribution function
g(x), and the static structure factor
S(q) are evaluated. The aim is to test Luttinger liquid (LL) predictions for 1D fluids over a wide range of fluid density and LL parameter
K. The 1D Bose fluid examined is a single chain of
4He atoms confined to a line in the center of a narrow nanopore. The atoms cannot exchange positions in the nanopore, the criterion for 1D. The fluid density is varied from the spinodal density where the 1D liquid is unstable to droplet formation to the density of bulk liquid
4He. In this range,
K varies from
K>2 at low density, where a robust superfluid is predicted, to
K<0.5, where fragile 1D superflow and solidlike peaks in
S(q) are predicted. For uniform pore walls, the ${\ensuremath{\rho}}_{S}/{\ensuremath{\rho}}_{0}$ scales as predicted by LL theory. The
n(x) and
g(x) show long range oscillations and decay with
x as predicted by LL theory. The amplitude of the oscillations is large at high density (small
K) and small at low density (large
K). The
K values obtained from different properties agree well verifying the internal structure of LL theory. In the presence of disorder, the ${\ensuremath{\rho}}_{S}/{\ensuremath{\rho}}_{0}$ does not scale as predicted by LL theory. A single
vJ parameter in the LL theory that recovers LL scaling was not found. The one body density matrix (OBDM) in disorder is well predicted by LL theory. The ``dynamical'' superfluid fraction, ${\ensuremath{\rho}}_{S}^{D}/{\ensuremath{\rho}}_{0}$, is determined. The physics of the deviation from LL theory in disorder and the ``dynamical'' ${\ensuremath{\rho}}_{S}^{D}/{\ensuremath{\rho}}_{0}$ are discussed.
We investigate the superfluid-insulator transition of one-dimensional interacting bosons in both deep and shallow periodic potentials. We compare a theoretical analysis based on quantum Monte Carlo simulations in continuum space and Luttinger liquid approach with experiments on ultracold atoms with tunable interactions and optical lattice depth. Experiments and theory are in excellent agreement. Our study provides a quantitative determination of the critical parameters for the Mott transition and defines the regimes of validity of widely used approximate models, namely, the Bose-Hubbard and sine-Gordon models.
We evaluate the Bose-Einstein condensate density and the superfluid fraction of bosons in a periodic external potential using path-integral Monte Carlo methods. The periodic lattice consists of a cubic cell containing a potential well that is replicated along one dimension using periodic boundary conditions. The aim is to describe bosons in a one-dimensional optical lattice or helium confined in a periodic porous medium. The one-body density matrix is evaluated and diagonalized numerically to obtain the single boson natural orbitals and the occupation of these orbitals. The condensate fraction is obtained as the fraction of bosons in the orbital that has the highest occupation. The superfluid density is obtained from the winding number. From the condensate orbital and superfluid fraction, we investigate (1) the impact of the periodic external potential on the spatial distribution of the condensate, and (2) the correlation of localizing the condensate into separated parts and the loss of superflow along the lattice. For high-density systems, as the well depth increases, the condensate becomes depleted in the wells and confined to the plateaus between successive wells, as in pores between necks in a porous medium. For low-density systems, as the well depth increases the BEC is localized at the center of the wells (tight binding) and depleted between the wells. In both cases, the localization of the condensate suppresses superflow leading to a superfluid-insulator crossover. The impact of the external potential on the temperature dependence of the superfluidity is also investigated. The external potential suppresses the superfluid fraction at all temperatures, with a superfluid fraction significantly less than one at low temperature. The addition of an external potential does not, however, significantly reduce the transition temperature.
We present results for the long-distance asymptotics of correlation functions of mesoscopic one-dimensional systems with periodic and open (Dirichlet) boundary conditions, as well as at finite temperature in the thermodynamic limit. The results are obtained using Haldane's harmonic-fluid approach (also known as 'bosonization'), and are valid for both bosons and fermions, in weakly and strongly interacting regimes. The harmonic-fluid approach and the method of computing the correlation functions using conformal transformations are explained in great detail. As an application relevant to one-dimensional systems of cold atomic gases, we consider the model of bosons interacting with a zero-range potential. The Luttinger-liquid parameters are obtained from the exact solution by solving the Bethe-ansatz equations in finite-size systems. The range of applicability of the approach is discussed, and the prefactor of the one-body density matrix of bosons is fixed by finding an appropriate parametrization of the weak-coupling result. The formula thus obtained is shown to be accurate, when compared with recent diffusion Monte Carlo calculations, within less than 10%. The experimental implications of these results for Bragg scattering experiments at low and high momenta are also discussed.
Ultracold bosonic and fermionic quantum gases confined to quasi-one-dimensional (1D) geometry are promising candidates for probing fundamental concepts of Luttinger liquid (LL) physics. They can also be exploited for devising applications in quantum information processing and precision measurements. Here, we focus on 1D dipolar Bose gases, where evidence of super-strong coupling behavior has been demonstrated by analyzing the low-energy static and dynamical structures of the fluid at zero temperature by a combined reptation quantum Monte Carlo (RQMC) and bosonization approach. Fingerprints of LL behavior emerge in the whole crossover from the already strongly interacting Tonks–Girardeau at low density to a dipolar density wave regime at high density. We have also shown that a LL framework can be effectively set up and utilized to describe this strongly correlated crossover physics in the case of confined 1D geometries after using the results for the homogeneous system in LL hydrodynamic equations within a local density approximation. This leads to the prediction of observable quantities such as the frequencies of the collective modes of the trapped dipolar gas under the more realistic conditions that could be found in ongoing experiments. The present paper provides a description of the theoretical framework in which the above results have been worked out, making available all the detailed derivations of the hydrodynamic Luttinger equations for the inhomogeneous trapped gas and of the correlation functions for the homogeneous system.
Abstract We present a comprehensive study shedding light on how thermal fluctuations affect correlations in a Bose gas with contact repulsive interactions in one spatial dimension. The pair correlation function, the static structure factor, and the one-body density matrix are calculated as a function of the interaction strength and temperature with the exact ab-initio Path Integral Monte Carlo method. We explore all possible gas regimes from weak to strong interactions and from low to high temperatures. We provide a detailed comparison with a number of theories, such as perturbative (Bogoliubov and decoherent classical), effective (Luttinger liquid) and exact (ground-state and thermal Bethe Ansatz) ones. Our Monte Carlo results exhibit an excellent agreement with the tractable limits and provide a fundamental benchmark for future observations which can be achieved in atomic gases, cavity quantum-electrodynamic and superconducting-circuit platforms.