Quantitative numerical analyses of interacting dilute Bose-Einstein condensates are most often based on semiclassical approximations. Since the complex-valued field-theoretic action of the Bose gas does not offer itself to the direct application of standard Monte Carlo techniques, simulations beyond their scope by now almost exclusively rely on quantum-mechanical techniques. Here we explore an alternative approach based on a Langevin-type sampling in an extended state space, known as the complex Langevin (CL) algorithm. While the use of the CL technique has a long-standing history in high-energy physics, in particular in the simulation of QCD at finite baryon density, applications to ultracold atoms are still in their infancy. Here we examine the applicability of the CL approach for a one- and two-component, three-dimensional nonrelativistic Bose gas in thermal equilibrium, below and above the Bose-Einstein phase transition. By comparison with analytic descriptions at the Gaussian level, including Bogoliubov and Hartree-Fock theory, we find that the method allows computing reliably and efficiently observables in the regime of experimentally accessible parameters. Close to the transition, quantum corrections lead to a shift of the critical temperature which we reproduce within the numerical range known in the literature. With this work, we aim to provide a first test of CL as a potential out-of-the-box tool for the simulation of experimentally realistic situations, including trapping geometries and multicomponent (or multiple species) models.
Quantitative numerical analyses of interacting dilute Bose-Einstein condensates are most often based on semiclassical approximations. Since the complex-valued field-theoretic action of the Bose gas does not offer itself to the direct application of standard Monte Carlo techniques, simulations beyond their scope by now almost exclusively rely on quantum-mechanical techniques. Here we explore an alternative approach based on a Langevin-type sampling in an extended state space, known as the complex Langevin (CL) algorithm. While the use of the CL technique has a long-standing history in high-energy physics, in particular in the simulation of QCD at finite baryon density, applications to ultracold atoms are still in their infancy. Here we examine the applicability of the CL approach for a one- and two-component, three-dimensional nonrelativistic Bose gas in thermal equilibrium, below and above the Bose-Einstein phase transition. By comparison with analytic descriptions at the Gaussian level, including Bogoliubov and Hartree-Fock theory, we find that the method allows computing reliably and efficiently observables in the regime of experimentally accessible parameters. Close to the transition, quantum corrections lead to a shift of the critical temperature which we reproduce within the numerical range known in the literature. With this work, we aim to provide a first test of CL as a potential out-of-the-box tool for the simulation of experimentally realistic situations, including trapping geometries and multicomponent (or multiple species) models.
Numerical simulations of the full quantum properties of interacting many-body systems by means of field-theoretic Monte Carlo techniques are often limited due to a sign problem. Here we simulate properties of a dilute two-dimensional Bose gas in the vicinity of the Berezinskii-Kosterlitz-Thouless (BKT) transition by means of the complex Langevin (CL) algorithm, thereby extending our previous CL study of the three-dimensional Bose gas to the lower-dimensional case. The purpose of the paper is twofold. On the one hand, it adds to benchmarking of the CL method and thus contributes to further exploring the range of applicability of the method. With the respective results, the universality of the equation of state is recovered, as well as the long-wave-length power-law dependence of the single-particle momentum spectrum below the BKT transition. Analysis of the rotational part of the current density corroborates vortex unbinding in crossing the transition. Beyond these measures of consistency we compute quantum corrections to the critical density and chemical potential in the weakly coupled regime. Our results show a shift of these quantities to lower values as compared to those obtained from classical field theory. It points in the opposite direction as compared to the shift of the critical density found by means of the path-integral Monte Carlo method at larger values of the coupling. Our simulations widen the perspective for precision comparisons with experiment.
In this thesis, we perform ab-initio lattice simulations of the nonrelativistic and relativistic Bose-Einstein condensates under strong external fields such as the electric field, the magnetic field, and the rotation.All the lattice actions discussed in this thesis are complex, so that the conventional Monte Carlo method su↵ers from the notorious sign problem.The path integral measure is complex due to the complex action, so that the important sampling algorithm cannot be applied.To overcome this di culty, we adopt the complex Langevin method, which has been developed in the context of finite density quantum chromodynamics to attack the fermion sign problem.In particular, we perform ab-initio simulation of the quantum vortex nucleation by the magnetic field or the rotation.In dilute and low temperature system, quantum and thermal fluctuations are negligible.The system can be remarkably well described by the mean field approximation.However, when quantum or thermal fluctuation becomes large, it is highly nontrivial how the quantized vortices behave.Around the critical values of temperature, chemical potential, magnetic field, or angular velocity, the fluctuation grows and the mean-field description inevitably breaks down.To understand the quantum vortex nucleation in such situations, the ab-initio simulation is necessary as discussed in this thesis.The direct evidence of the quantum vortex nucleation is the quantization of circulation.We show that the circulation is quantized in the superconducting (superfluid) phase far from the critical chemical potential, but it is not just above the critical chemical potential.We observe that the quantized circulation is blurred by the quantum fluctuation as the chemical potential is getting close to its critical value.In our simulation, the quantum fluctuation of vortex number is observed, while the averaged circulation is clearly quantized.At first glance, these two facts seem to be incompatible.However, we show that the fluctuation of vortex number behaves as Gaussian and, as a result of cancellation, the average circulation becomes exact integer.iii
This article reviews recent advances in the theory of the three-dimensional dilute homogeneous Bose gas at zero and finite temperature. Effective-field-theory methods are used to formulate a systematic perturbative framework that can be used to calculate the properties of a system at T=0. The perturbative expansion of these properties is essentially an expansion in the gas parameter na3, where a is the s-wave scattering length and n is the number density. In particular, the leading quantum corrections to the ground-state energy density, the condensate depletion, and long-wavelength collective excitations are rederived in an efficient and economical manner. Nonuniversal effects are also discussed. These effects are higher-order corrections that depend on properties of the interatomic potential other than the scattering length, such as the effective range. The article critically examines various approaches to the dilute Bose gas in equilibrium at finite temperature. These include the Bogoliubov approximation, the Popov approximation, the Hartree-Fock-Bogoliubov approximation, the \ensuremath{\Phi}-derivable approach, optimized perturbation theory, and renormalization-group techniques. The article ends with a look at recent calculations of the critical temperature of the dilute Bose gas, which include 1/N techniques, lattice simulations, self-consistent calculations, and variational perturbation theory.
The quantum many-body problem is at the heart of some of the most formidable open problems in modern physics, such as high-Tc superconductivity or the behaviour of neutron stars. Ultracold atomic systems can now be used to simulate model Hamiltonians of condensed matter or nuclear physics, in very well-controlled environnement. In this thesis, we have developed a general method to probe the thermodynamics of homogeneous quantum systems using trapped atomic gases. These measurements are directly compared to the predictions of theories of the quantum many-body problem. We have applied this technique to the spin-1/2 Fermi gas and the Bose with short-range interactions. Using fermionic 6Li, we explored a part of the wide parameter space by changing the interaction strength, the spin-population imbalance or the temperature of the gas. This system exhibits remarkably rich physics, such as normal/superfluid phase transitions (that can be of thermal or quantum character) or a Fermi liquid-type behaviour of the normal phase. We have also used this method to probe a Bose gas of 7Li atoms close to a Feshbach resonance. We have measured the Equation of State of the Bose gas as a function of interactions at very low temperature. For the first time, we measured quantitatively the Lee-Huang-Yang beyond mean-field correction to the ground-state energy of the system, first predicted in 1957. We compared the experimental results to Quantum Monte-Carlo calculations. We have extended this study using out-of-equilibrium measurements of the Bose gas in the strongly interacting regime, which gives a first hint on the properties of the hypothetical unitary Bose gas.
Using a finite-temperature path integral Monte Carlo simulation (PIMC) method and finite-size scaling, we have investigated the interaction-induced shift of the phase-transition temperature for Bose-Einstein condensation of homogeneous weakly interacting Bose gases in three dimensions, which is given by a proposed analytical expression ${T}_{c}={T}_{c}^{0}{1+{c}_{1}a{n}^{1∕3}+[{c}_{2}^{\ensuremath{'}}\phantom{\rule{0.2em}{0ex}}\mathrm{ln}(a{n}^{1∕3})+{c}_{2}^{\ensuremath{''}}]{a}^{2}{n}^{2∕3}+\mathrm{O}({a}^{3}n)}$, where Tc0 is the critical temperature for an ideal gas, a is the s-wave scattering length, and n is the number density. We have used smaller number densities and more time slices than in the previous PIMC simulations [Gruter et al., Phys. Rev. Lett. 79, 3549 (1997)] in order to understand the difference in the value of the coefficient c1 between their results and the (apparently) other reliable results in the literature. Our results show that ${({T}_{c}\ensuremath{-}{T}_{c}^{0})∕{T}_{c}^{0}}∕(a{n}^{1∕3})$ depends strongly on the interaction strength an1∕3 while the previous PIMC results are considerably flatter and smaller than our results. We obtain ${c}_{1}=1.32\ifmmode\pm\else\textpm\fi{}0.14$, in agreement with results from recent Monte Carlo methods of three-dimensional O(2) scalar ${\ensuremath{\phi}}^{4}$ field theory and variational perturbation theory.
The mean-field properties of trapped, two-dimensional interacting Bose gases of finite size are surveyed numerically. The coupled nonlinear equations of the Gross-Pitaevskii (GP) equation and the Hartree-Fock equation for the condensate and the thermal component, respectively, are calculated by the hybrid method, in which the GP equation and the low-lying excitations less than cut-off energy are calculated quantum mechanically using the basis representation, while the high-lying states are obtained within the local density approximation. Consistent solutions are obtained for the whole range of temperatures including the critical temperature by determining the thermodynamic chemical potential µ from the condition that the total atomic number must be the desired value, instead of identifying µ with the eigenvalue of the GP equation. Comparisons with the results of the semiclassical two-fluid model are presented.