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.
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.
In this work we investigate the unique properties of ultracold Bose gases in one and two dimensions. In two dimensions, we present simulations of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition using the projected Gross-Pitaevskii equation formalism. We confirm the presence of the BKT phase by observing algebraic decay in the first-order correlation function, unbinding of vortex pairs, and behaviour of the superfluid fraction. In one dimension we use a transverse Gaussian ansatz and the Lagrangian approach to derive an effective equation of motion for the one dimensional Bose gas at zero temperature.
This thesis concerns ultracold two-dimensional atomic gases and the associated novel superfluid phases in an optical lattice subjected to staggered rotations. In the first part, we study the Berezinskii-Kosterlitz-Thouless (BKT) transition in a two-dimensional Bose gas. Motivated by recent experiments on the cold gas in the BKT regime, we employ a modified Popov theory and a renormalization group theory to describe the system. The two-dimensional renormalization group flow equations display unique features, which we interpret as precursors of the BKT physics. We further calculate higher correlation functions close to the BKT transition, which can serve as observables for detecting effects beyond mean-field. Finally, we obtain the density profile under the conditions of the experiment performed in the group of Dalibard in Paris, France. In the second part of the thesis, we study a recent experimental proposal of a time-dependent optical lattice with staggered rotations. We show that in the tight-binding regime, the system can be effectively described by a (Bose-)Hubbard model subjected to an artificial staggered magnetic field. The staggered magnetic field (flux) in the lattice leads to conical points in the band structure which possesses distinct minima for different strength, and can give rise to Dirac-like dispersions. For bosons, we show that besides the homogenous superfluid phase, a staggered-vortex superfluid phase is stabilized for large flux values. We show that the staggered flux renormalizes the phase boundary and find a multi-critical point separating the two superfluid phases and the Mott insulating phase. Finally, the full phase diagram of the generalized Bose-Hubbard model at zero temperature is determined. We then generalize the model for a mixture of ultracold bosons and spin-1/2 fermions. By using the Bogoliubov theory for the bosonic sector, we show that the bosonic superfluid mediates attractive interactions between fermions up to the nearest-neighbor site. The resulting model is an extended Hubbard model in the presence of a staggered magnetic field, where the on-site and the nearest-neighbor attractive interactions can be tuned independently. We then employ the Bardeen-Cooper-Schrieffer theory for superconductivity, by considering an on-site pairing order parameter competing with a nearest-neighbor resonating-valence-bond order parameter. Due to conical points in the energy band, the system possesses quantum critical property and the resulting phase diagram displays features which are strikingly similar to that of high-critical-temperature superconductors or heavy-fermion superconductors.
Current experiments on atomic gases in highly anisotropic traps present the opportunity to study in detail the low temperature phases of two-dimensional inhomogeneous systems. Although, in an ideal gas, the trapping potential favors Bose-Einstein condensation at finite temperature, interactions tend to destabilize the condensate, leading to a superfluid Kosterlitz-Thouless-Berezinskii phase with a finite superfluid mass density but no long-range order, as in homogeneous fluids. The transition in homogeneous systems is conveniently described in terms of dissociation of topological defects (vortex-antivortex pairs). However, trapped two-dimensional gases are more directly approached by generalizing the microscopic theory of the homogeneous gas. In this paper, we first derive, via a diagrammatic expansion, the scaling structure near the phase transition in a homogeneous system, and then study the effects of a trapping potential in the local density approximation. We find that a weakly interacting trapped gas undergoes a Kosterlitz-Thouless-Berezinskii transition from the normal state at a temperature slightly below the Bose-Einstein transition temperature of the ideal gas. The characteristic finite superfluid mass density of a homogeneous system just below the transition becomes strongly suppressed in a trapped gas.
This paper develops and implements the stochastic projected Gross-Pitaevskii equation for spin-1 Bose gases, addressing key considerations for numerical simulations. As an application of the theory we explore equilibrium phases in a two-dimensional spin-1 gas, where quasi-long-range order emerges via a Berezinskii-Kosterlitz-Thouless transition. Our analysis includes the definition of superfluid densities for both mass and spin degrees of freedom in a manner suitable for implementation within a stochastic projected Gross-Pitaevskii equation simulation. We present a finite-temperature phase diagram for the ferromagnetic spin-1 Bose gas and identify three distinct superfluid phases: two exhibiting conventional Berezinskii-Kosterlitz-Thouless-like behavior and an additional phase that simultaneously supports independent mass and spin superflows. As temperature increases, the stability region of this additional phase shrinks. We provide a characterization of the phase transitions through consideration of the spin-component densities and the unbinding of multiple types of vortices. This work provides a foundation for further studies of nonequilibrium and finite-temperature phenomena in spinor Bose gases.
Abstract Characterizing the superconducting and superfluid transitions in two-dimensional (2D) many-body systems is of broad interest and remains a fundamental issue. In this study, we establish the condensate fraction scaling as a highly efficient tool to achieve that and accordingly propose efficient schemes to accurately determine the associated Berezinskii−Kosterlitz−Thouless (BKT) transitions. Using the 2D attractive Fermi-Hubbard model as a testbed and applying numerically exact auxiliary-field quantum Monte Carlo simulations, we access unprecedented system sizes (up to 64 × 64 = 4096 lattice sites) and perform a comprehensive analysis for the temperature dependence and finite-size scaling of condensate fraction across the BKT transition. We demonstrate that this quantity exhibits algebraic scaling below the transition and exponential scaling above it, with significantly reduced finite-size effects comparing to the extensively studied on-site pairing correlator. This greatly improves the determination of BKT transition using moderate system sizes. We also extract finite-size BKT transition temperature from condensate fraction, and confirm its logarithmic correction on system size. Based on the accurately determined transition, we reveal that the specific heat displays an anomaly, showing a peak at a temperature slightly above BKT transition. Our findings should be generally applicable to 2D fermionic and bosonic systems hosting superconductivity or superfluidity.
Spinor Bose gases may exhibit both phase and ferromagnetic order. Consequently, they may support superflow of both mass and spin currents. In this thesis we utilize a stochastic Gross-Pitaevskii model to investigate the superfluid properties of a two-dimensional spin-1 Bose gas with ferromagnetic interactions. Owing to the low dimensionality, the superfluidity of this system can not be attributed to the presence of a Bose-Einstein condensate, with thermal fluctuations precluding long-range order. Rather, mass and spin superfluidity emerge via respective Berezinskii-Kosterlitz-Thouless (BKT) transitions. Such transitions are characterized by the emergence of quasi-long-range order, driven by the binding of topological defect pairs. First, we consider an easy-plane ferromagnetic gas. We find the system superfluidity may be understood by direct application of existing BKT theory. In particular, mass and spin superfluidity are respectively associated with quasi-long-range order in the global phase and local magnetization. The former arises from the binding of mass vortex pairs, while the latter arises from the binding of polar-core spin vortex pairs. Following this, we introduce axial magnetization, and consider a gas in the broken-axisymmetric phase. We find axial magnetization acts to couple mass and spin superflows via a mechanism analogous to the superfluid drag exhibited by two-component superfluids. This coupling renders invalid the direct applicability of existing BKT theory: here mass and spin superfluidity may each be associated with both phase and magnetic order. We remedy this with a generalization making concrete the relation between superfluidity and quasi-long-range order in this system. In addition, through variations of the quadratic Zeeman energy and axial magnetization we identify three distinct superfluid phases, which serve as finite-temperature generalizations of the polar, broken-axisymmetric, and easy-axis ferromagnetic ground states.