Numerical method for evolving the Projected Gross-Pitaevskii equation

In this paper we describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for a Bose gas in a harmonic oscillator potential. The central difficulty in solving this equation is the requirement that the classical field is restricted to a small set of prescribed modes that constitute the low energy classical region of the system. We present a scheme, using a Hermite-polynomial based spectral representation, that precisely implements this mode restriction and allows an efficient and accurate solution of the PGPE. We show equilibrium and non-equilibrium results from the application of the PGPE to an anisotropic trapped three-dimensional Bose gas.Comment: 12 pages, 5 figures. To appear in Phys. Rev. E. Convergence results added, a few minor changes made and typos fixe

Quantum turbulence and correlations in Bose-Einstein condensate collisions

We investigate numerically simulated collisions between experimentally realistic Bose-Einstein condensate wavepackets, within a regime where highly populated scattering haloes are formed. The theoretical basis for this work is the truncated Wigner method, for which we present a detailed derivation, paying particular attention to its validity regime for colliding condensates. This paper is an extension of our previous Letter [A. A. Norrie, R. J. Ballagh, and C. W. Gardiner, Phys. Rev. Lett. 94, 040401 (2005)] and we investigate both single-trajectory solutions, which reveal the presence of quantum turbulence in the scattering halo, and ensembles of trajectories, which we use to calculate quantum-mechanical correlation functions of the field

Classical Region of a Trapped Bose Gas

The classical region of a Bose gas consists of all single-particle modes that have a high average occupation and are well-described by a classical field. Highly-occupied modes only occur in massive Bose gases at ultra-cold temperatures, in contrast to the photon case where there are highly-occupied modes at all temperatures. For the Bose gas the number of these modes is dependent on the temperature, the total number of particles and their interaction strength. In this paper we characterize the classical region of a harmonically trapped Bose gas over a wide parameter regime. We use a Hartree-Fock approach to account for the effects of interactions, which we observe to significantly change the classical region as compared to the idealized case. We compare our results to full classical field calculations and show that the Hartree-Fock approach provides a qualitatively accurate description of classical region for the interacting gas.Comment: 6 pages, 5 figures; updated to include new results with interaction

Microcanonical temperature for a classical field: application to Bose-Einstein condensation

We show that the projected Gross-Pitaevskii equation (PGPE) can be mapped exactly onto Hamilton's equations of motion for classical position and momentum variables. Making use of this mapping, we adapt techniques developed in statistical mechanics to calculate the temperature and chemical potential of a classical Bose field in the microcanonical ensemble. We apply the method to simulations of the PGPE, which can be used to represent the highly occupied modes of Bose condensed gases at finite temperature. The method is rigorous, valid beyond the realms of perturbation theory, and agrees with an earlier method of temperature measurement for the same system. Using this method we show that the critical temperature for condensation in a homogeneous Bose gas on a lattice with a UV cutoff increases with the interaction strength. We discuss how to determine the temperature shift for the Bose gas in the continuum limit using this type of calculation, and obtain a result in agreement with more sophisticated Monte Carlo simulations. We also consider the behaviour of the specific heat.Comment: v1: 9 pages, 5 figures, revtex 4. v2: additional text in response to referee's comments, now 11 pages, to appear in Phys. Rev.

Simulations of thermal Bose fields in the classical limit

We demonstrate that the time-dependent projected Gross-Pitaevskii equation derived earlier [Davis, et al., J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. We find that this equation will evolve randomised initial wave functions to equilibrium, and compare our numerical data to the predictions of a gapless, second-order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)]. We find that we can determine the temperature of the equilibrium state when this theory is valid. Outside the range of perturbation theory we describe how to measure the temperature of our simulations. We also determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. As the Gross-Pitaevskii equation is non-perturbative, we expect that it can describe the correct thermal behaviour of a Bose gas as long as all relevant modes are highly occupied.Comment: 15 pages, 12 figures, revtex4, follow up to Phys. Rev. Lett. 87 160402 (2001). v2: Modified after referee comments. Extra data added to two figures, section on temperature determination expande