We review phase space techniques based on the Wigner representation that provide an approximate description of dilute ultra-cold Bose gases. In this approach the quantum field evolution can be represented using equations of motion of a similar form to the Gross-Pitaevskii equation but with stochastic modifications that include quantum effects in a controlled degree of approximation. These techniques provide a practical quantitative description of both equilibrium and dynamical properties of Bose gas systems. We develop versions of the formalism appropriate at zero temperature, where quantum fluctuations can be important, and at finite temperature where thermal fluctuations dominate. The numerical techniques necessary for implementing the formalism are discussed in detail, together with methods for extracting observables of interest. Numerous applications to a wide range of phenomena are presented.

We consider the ground state properties of a trapped dipolar condensate under the influence of quantum fluctuations. We show that this system can undergo a phase transition from a low density condensate state to a high density droplet state, which is stabilized by quantum fluctuations. The energetically favored state depends on the geometry of the confining potential, the number of atoms and the two-body interactions. We develop a simple variational ansatz and validate it against full numerical solutions. We produce a phase diagram for the system and present results relevant to current experiments with dysprosium and erbium condensates.

We extend the projected Gross-Pitaevskii equation formalism of Davis et al. ͓Phys. Rev. Lett. 87, 160402 ͑2001͔͒ to the experimentally relevant case of thermal Bose gases in harmonic potentials and outline a robust and accurate numerical scheme that can efficiently simulate this system. We apply this method to investigate the equilibrium properties of the harmonically trapped three-dimensional projected Gross-Pitaevskii equation at finite temperature and consider the dependence of condensate fraction, position, and momentum distributions and density fluctuations on temperature. We apply the scheme to simulate an evaporative cooling process in which the preferential removal of high-energy particles leads to the growth of a Bose-Einstein condensate. We show that a condensate fraction can be inferred during the dynamics even in this nonequilibrium situation.