We predict the existence of a dip below unity in the second-order coherence function of a partially condensed ideal Bose gas in harmonic confinement, signaling the anticorrelation of density fluctuations in the sample. The dip in the second-order coh
erence function is revealed in a canonical-ensemble calculation, corresponding to a system with fixed total number of particles. In a grand-canonical ensemble description, this dip is obscured by the occupation-number fluctuation catastrophe of the ideal Bose gas. The anticorrelation is most pronounced in highly anisotropic trap geometries containing small particle numbers. We explain the fundamental physical mechanism which underlies this phenomenon, and its relevance to experiments on interacting Bose gases.
The classical-field formalism has been widely applied in the calculation of normal correlation functions, and the characterization of condensation, in finite-temperature Bose gases. Here we discuss the extension of this method to the calculation of m
ore general correlations, including the so-called anomalous correlations of the field, without recourse to symmetry-breaking assumptions. Our method is based on the introduction of U(1)-symmetric classical-field variables analogous to the modified quantum ladder operators of number-conserving approaches to the degenerate Bose gas, and allows us to rigorously quantify the anomalous and non-Gaussian character of the field fluctuations. We compare our results for anomalous correlation functions with the predictions of mean-field theories, and demonstrate that the nonlinear classical-field dynamics incorporate a full description of many-body processes which modify the effective mean-field potentials experienced by condensate and noncondensate atoms. We discuss the role of these processes in shaping the condensate mode, and thereby demonstrate the consistency of the Penrose-Onsager definition of the condensate orbital in the classical-field equilibrium. We consider the contribution of various noncondensate-field correlations to the overall suppression of density fluctuations and interactions in the field, and demonstrate the distinct roles of phase and density fluctuations in the transition of the field to the normal phase.
The coherence properties of degenerate Bose gases have usually been expressed in terms of spatial correlation functions, neglecting the rich information encoded in their temporal behavior. In this paper we show, using a Hamiltonian classical-field fo
rmalism, that temporal correlations can be used to characterize familiar properties of a finite-temperature degenerate Bose gas. The temporal coherence of a Bose-Einstein condensate is limited only by the slow diffusion of its phase, and thus the presence of a condensate is indicated by a sharp feature in the temporal power spectrum of the field. We show that the condensate mode can be obtained by averaging the field for a short time in an appropriate phase-rotating frame, and that for a wide range of temperatures, the condensate obtained in this approach agrees well with that defined by the Penrose-Onsager criterion based on one-body (spatial) correlations. For time periods long compared to the phase diffusion time, the field will average to zero, as we would expect from the overall U(1) symmetry of the Hamiltonian. We identify the emergence of the first moment on short time scales with the concept of U(1) symmetry breaking that is central to traditional mean-field theories of Bose condensation. We demonstrate that the short-time averaging procedure constitutes a general analog of the anomalous averaging operation of symmetry-broken theories by calculating the anomalous thermal density of the field, which we find to have form and temperature dependence consistent with the results of mean-field theories.
We perform finite-temperature dynamical simulations of the arrest of a rotating Bose-Einstein condensate by a fixed trap anisotropy, using a Hamiltonian classical-field method. We consider a quasi-two-dimensional condensate containing a single vortex
in equilibrium with a rotating thermal cloud. Introducing an elliptical deformation of the trapping potential leads to the loss of angular momentum from the system. We identify the condensate and the complementary thermal component of the nonequilibrium field, and compare the evolution of their angular momenta and angular velocities. By varying the trap anisotropy we alter the relative efficiencies of the vortex-cloud and cloud-trap coupling. For strong trap anisotropies the angular momentum of the thermal cloud may be entirely depleted before the vortex begins to decay. For weak trap anisotropies, the thermal cloud exhibits a long-lived steady state in which it rotates at an intermediate angular velocity.
We consider a finite-temperature Bose-Einstein condensate in a quasi-two-dimensional trap containing a single precessing vortex. We find that such a configuration arises naturally as an ergodic equilibrium of the projected Gross-Pitaevskii equation,
when constrained to a finite conserved angular momentum. In an isotropic trapping potential the condensation of the classical field into an off-axis vortex state breaks the rotational symmetry of the system. We present a methodology to identify the condensate and the Goldstone mode associated with the broken rotational symmetry in the classical-field model. We also examine the variation in vortex trajectories and thermodynamic parameters of the field as the energy of the microcanonical field simulation is varied.
We present a quantum mechanical treatment of the mechanical stirring of Bose-Einstein condensates using classical field techniques. In our approach the condensate and excited modes are described using a Hamiltonian classical field method in which the
atom number and (rotating frame) energy are strictly conserved. We simulate a T = 0 quasi-2D condensate perturbed by a rotating anisotropic trapping potential. Vacuum fluctuations in the initial state provide an irreducible mechanism for breaking the initial symmetries of the condensate and seeding the subsequent dynamical instability. Highly turbulent motion develops and we quantify the emergence of a rotating thermal component that provides the dissipation necessary for the nucleation and motional-damping of vortices in the condensate. Vortex lattice formation is not observed, rather the vortices assemble into a spatially disordered vortex liquid state. We discuss methods we have developed to identify the condensate in the presence of an irregular distribution of vortices, determine the thermodynamic parameters of the thermal component, and extract damping rates from the classical field trajectories.
