Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0705.1226 · Physical Review A

Semiclassical field method for the equilibrium Bose gas and application to thermal vortices in two dimensions

Luca Giorgetti, Iacopo Carusotto, Yvan Castin

Published Jul 16, 2007• 32 citations
arXivHTMLPDF

Abstract

We develop a semiclassical field method for the study of the weakly interacting Bose gas at finite temperature which, contrary to the usual classical field model, does not suffer from an ultraviolet cutoff dependence. We apply the method to the study of thermal vortices in spatially homogeneous, two-dimensional systems. We present numerical results for the vortex density and the vortex pair distribution function. Insight in the physics of the system is obtained by comparing the numerical results with the predictions of simple analytical models. In particular, we calculate the activation energy required to form a vortex pair at low temperature.

DOI
10.1103/physreva.76.013613

Related papers

Vortex pairing in two-dimensional Bose gases

Christopher J. Foster, P. B. Blakie, Matthew J. Davis

Recent experiments on ultracold Bose gases in two dimensions have provided evidence for the existence of the Berezinskii-Kosterlitz-Thouless (BKT) phase via analysis of the interference between two independent systems. In this work we study the two-dimensional quantum degenerate Bose gas at finite temperature using the projected Gross-Pitaevskii equation classical field method. Although this describes the highly occupied modes of the gas below a momentum cutoff, we have developed a method to incorporate the higher momentum states in our model. We concentrate on finite-sized homogeneous systems in order to simplify the analysis of the vortex pairing. We determine the dependence of the condensate fraction on temperature and compare this to the calculated superfluid fraction. By measuring the first order correlation function we determine the boundary of the Bose-Einstein condensate and BKT phases, and find it is consistent with the superfluid fraction decreasing to zero. We reveal the characteristic unbinding of vortex pairs above the BKT transition via a coarse-graining procedure. Finally, we model the procedure used in experiments to infer system correlations [Hadzibabic et al., Nature 441, 1118 (2006)], and quantify its level of agreement with directly calculated in situ correlation functions.

arXiv:0912.1675•Physical Review A•Feb 24, 2010•64 citations

PHASE TRANSITIONS IN ULTRA-COLD TWO-DIMENSIONAL BOSE GASES

David A. W. Hutchinson, P. B. Blakie

We briefly review the theory of Bose-Einstein condensation in the two-dimensional trapped Bose gas and, in particular the relationship to the theory of the homogeneous two-dimensional gas and the Berezinskii-Kosterlitz-Thouless phase. We obtain a phase diagram for the trapped two-dimensional gas, finding a critical temperature above which the free energy of a state with a pair of vortices of opposite circulation is lower than that for a vortex-free Bose-Einstein condensed ground state. We identify three distinct phases which are, in order of increasing temperature, a phase coherent Bose-Einstein condensate, a vortex pair plasma with fluctuating condensate phase and a thermal Bose gas. The thermal activation of vortex-antivortex pair formation is confirmed using finite-temperature classical field simulations.

arXiv:cond-mat/0602614•International Journal of Modern Physics B•Dec 15, 2006•7 citations

Observation of Thermally Activated Vortex Pairs in a Quasi-2D Bose Gas

Jae-yoon Choi, Sang Won Seo, Yong-il Shin

We measure the in-plane distribution of thermally activated vortices in a trapped quasi-2D Bose gas, where we enhance the visibility of density-depleted vortex cores by radially compressing the sample before releasing the trap. The pairing of vortices is revealed by the two-vortex spatial correlation function obtained from the vortex distribution. The vortex density decreases gradually as temperature is lowered, and below a certain temperature, a vortex-free region emerges in the center of the sample. This shows the crossover from a Berezinskii-Kosterlitz-Thouless phase containing vortex-pair excitations to a vortex-free Bose-Einstein condensate in a finite-size 2D system.

arXiv:1211.5649•Physical Review Letters•Apr 25, 2013•51 citations

Simulations of thermal Bose fields in the classical limit

Matthew J. Davis, S. A. Morgan, Keith R. Burnett

We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.

arXiv:cond-mat/0201571•Physical Review A•Nov 25, 2002•103 citations

An exact stochastic field method for the interacting Bose gas at thermal equilibrium

Iacopo Carusotto, Yvan Castin

We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.

arXiv:cond-mat/0108042•Journal of Physics B Atomic Molecular and Optical Physics•Nov 23, 2001•37 citations

Thermal and superfluid properties of the two-dimensional Bose gas

Rémi Desbuquois, Jean Dalibard

The physical properties of a homogeneous system at the thermodynamic equilibrium are strongly constrained by its dimensionality. In this respect, the two dimensional Bose gas constitutes a particularly interesting system: while it is impossible to observe a long range order at non-zero temperatures, there exists nevertheless a phase transition to a superfluid state at low temperatures. Furthermore, owing to the reduced dimensionality, the equation of state of the weakly interacting two-dimensional Bose gas is scale invariant. In this thesis, we present an experimental study of the two-dimensional Bose gas. We measure its equation of state with two different methods, and find a good agreement with analytic and numerical predictions. These results confirm as well the scale invariance of the system. Furthermore, one of the methods allows for a determination of the equation of state with a single adjustable parameter. We then characterize the superfluid response of the system., by showing evidence for a dissipationless response of the system to a moving perturbation. Finally, we analyze the fluctuations dynamics of the 2D Bose gas, which confirms both the suppression of density fluctuations in the superfluid phase, and the dominating contribution of phonons to the phase fluctuations.

•Jan 1, 2013•0 citations
We develop a semiclassical field method for the study of the weakly interacting Bose gas at finite temperature which, contrary to the usual classical field model, does not suffer from an ultraviolet cutoff dependence. We apply the method to the study of thermal vortices in spatially homogeneous, two-dimensional systems. We present numerical results for the vortex density and the vortex pair distribution function. Insight in the physics of the system is obtained by comparing the numerical results with the predictions of simple analytical models. In particular, we calculate the activation energy required to form a vortex pair at low temperature.
Recent experiments on ultracold Bose gases in two dimensions have provided evidence for the existence of the Berezinskii-Kosterlitz-Thouless (BKT) phase via analysis of the interference between two independent systems. In this work we study the two-dimensional quantum degenerate Bose gas at finite temperature using the projected Gross-Pitaevskii equation classical field method. Although this describes the highly occupied modes of the gas below a momentum cutoff, we have developed a method to incorporate the higher momentum states in our model. We concentrate on finite-sized homogeneous systems in order to simplify the analysis of the vortex pairing. We determine the dependence of the condensate fraction on temperature and compare this to the calculated superfluid fraction. By measuring the first order correlation function we determine the boundary of the Bose-Einstein condensate and BKT phases, and find it is consistent with the superfluid fraction decreasing to zero. We reveal the characteristic unbinding of vortex pairs above the BKT transition via a coarse-graining procedure. Finally, we model the procedure used in experiments to infer system correlations [Hadzibabic et al., Nature 441, 1118 (2006)], and quantify its level of agreement with directly calculated in situ correlation functions.
We briefly review the theory of Bose-Einstein condensation in the two-dimensional trapped Bose gas and, in particular the relationship to the theory of the homogeneous two-dimensional gas and the Berezinskii-Kosterlitz-Thouless phase. We obtain a phase diagram for the trapped two-dimensional gas, finding a critical temperature above which the free energy of a state with a pair of vortices of opposite circulation is lower than that for a vortex-free Bose-Einstein condensed ground state. We identify three distinct phases which are, in order of increasing temperature, a phase coherent Bose-Einstein condensate, a vortex pair plasma with fluctuating condensate phase and a thermal Bose gas. The thermal activation of vortex-antivortex pair formation is confirmed using finite-temperature classical field simulations.
We measure the in-plane distribution of thermally activated vortices in a trapped quasi-2D Bose gas, where we enhance the visibility of density-depleted vortex cores by radially compressing the sample before releasing the trap. The pairing of vortices is revealed by the two-vortex spatial correlation function obtained from the vortex distribution. The vortex density decreases gradually as temperature is lowered, and below a certain temperature, a vortex-free region emerges in the center of the sample. This shows the crossover from a Berezinskii-Kosterlitz-Thouless phase containing vortex-pair excitations to a vortex-free Bose-Einstein condensate in a finite-size 2D system.
We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.
We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.
The physical properties of a homogeneous system at the thermodynamic equilibrium are strongly constrained by its dimensionality. In this respect, the two dimensional Bose gas constitutes a particularly interesting system: while it is impossible to observe a long range order at non-zero temperatures, there exists nevertheless a phase transition to a superfluid state at low temperatures. Furthermore, owing to the reduced dimensionality, the equation of state of the weakly interacting two-dimensional Bose gas is scale invariant. In this thesis, we present an experimental study of the two-dimensional Bose gas. We measure its equation of state with two different methods, and find a good agreement with analytic and numerical predictions. These results confirm as well the scale invariance of the system. Furthermore, one of the methods allows for a determination of the equation of state with a single adjustable parameter. We then characterize the superfluid response of the system., by showing evidence for a dissipationless response of the system to a moving perturbation. Finally, we analyze the fluctuations dynamics of the 2D Bose gas, which confirms both the suppression of density fluctuations in the superfluid phase, and the dominating contribution of phonons to the phase fluctuations.