The stochastic Gross-Pitaevskii equation: II

We provide a derivation of a more accurate version of the stochastic Gross-Pitaevskii equation, as introduced by Gardiner et al (2002 J. Phys. B: At. Mol. Opt. Phys. 35 1555). This derivation does not rely on the concept of local energy and momentum conservation and is based on a quasiclassical Wigner function representation of a 'high temperature' master equation for a Bose gas, which includes only modes below an energy cut-off ER that are sufficiently highly occupied (the condensate band). The modes above this cutoff (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provides noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross-Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation and by the feasibility of its numerical implementation.

Experimental tests of quantum nonlinear dynamics in atom optics

Cold atoms in optical potentials provide an ideal test bed to explore quantum nonlinear dynamics. Atoms are prepared in a magneto-optic trap or as a dilute Bose-Einstein condensate and subjected to a far detuned optical standing wave that is modulated. They exhibit a wide range of dynamics, some of which can be explained by classical theory while other aspects show the underlying quantum nature of the system. The atoms have a mixed phase space containing regions of regular motion which appear as distinct peaks in the atomic momentum distribution embedded in a sea of chaos. The action of the atoms is of the order of Planck's constant, making quantum effects significant. This tutorial presents a detailed description of experiments measuring the evolution of atoms in time-dependent optical potentials. Experimental methods are developed providing means for the observation and selective loading of regions of regular motion. The dependence of the atomic dynamics on the system parameters is explored and distinct changes in the atomic momentum distribution are observed which are explained by the applicable quantum and classical theory. The observation of a bifurcation sequence is reported and explained using classical perturbation theory. Experimental methods for the accurate control of the momentum of an ensemble of atoms are developed. They use phase space resonances and chaotic transients providing novel ensemble atomic beamsplitters. The divergence between quantum and classical nonlinear dynamics is manifest in the experimental observation of dynamical tunnelling. It involves no potential barrier. However a constant of motion other than energy still forbids classically this quantum allowed motion. Atoms coherently tunnel back and forth between their initial state of oscillatory motion and the state 180 out of phase with the initial state.

Quantum dynamics with stochastic gauge simulations

The general idea of a stochastic gauge representation is introduced and compared with more traditional phase-space expansions, like the Wigner expansion. Stochastic gauges can be used to obtain an infinite class of positive-definite stochastic time-evolution equations, equivalent to master equations, for many systems including quantum time evolution. The method is illustrated with a variety of simple examples ranging from astrophysical molecular hydrogen production, through to the topical problem of Bose-Einstein condensation in an optical trap and the resulting quantum dynamics.

Properties of the stochastic Gross-Pitaevskii equation: finite temperature Ehrenfest relations and the optimal plane wave representation

We present Ehrenfest relations for the high temperature stochastic Gross-Pitaevskii equation description of a trapped Bose gas, including the effect of growth noise and the energy cutoff. A condition for neglecting the cutoff terms in the Ehrenfest relations is found which is more stringent than the usual validity condition of the truncated Wigner or classical field method-that all modes are highly occupied. The condition requires a small overlap of the nonlinear interaction term with the lowest energy single particle state of the noncondensate band, and gives a means to constrain dynamical artefacts arising from the energy cutoff in numerical simulations. We apply the formalism to two simple test problems: (i) simulation of the Kohn mode oscillation for a trapped Bose gas at zero temperature, and (ii) computing the equilibrium properties of a finite temperature Bose gas within the classical field method. The examples indicate ways to control the effects of the cutoff, and that there is an optimal choice of plane wave basis for a given cutoff energy. This basis gives the best reproduction of the single particle spectrum, the condensate fraction and the position and momentum densities.

Atom counting in ultracold gases using photoionization and ion detection

We analyze photoionization and ion detection as a means of accurately counting ultracold atoms. We show that it is possible to count clouds containing many thousands of atoms with accuracies better than N-1/2 with current technology. This allows the direct probing of sub-Poissonian number statistics of atomic samples. The scheme can also be used for efficient single-atom detection with high spatiotemporal resolution. All aspects of a realistic detection scheme are considered, and we discuss experimental situations in which such a scheme could be implemented.

Observing a coherent superposition of an atom and a molecule

We demonstrate that it is possible, in principle, to perform a Ramsey-type interference experiment to exhibit a coherent superposition of a single atom and a diatomic molecule. This gedanken experiment, based on the techniques of Aharonov and Susskind [Phys. Rev. 155, 1428 (1967)], explicitly violates the commonly accepted superselection rule that forbids coherent superpositions of eigenstates of differing atom number. A Bose-Einstein condensate plays the role of a reference frame that allows for coherent operations analogous to Ramsey pulses. We also investigate an analogous gedanken experiment to exhibit a coherent superposition of a single boson and a fermion, violating the commonly accepted superselection rule forbidding coherent superpositions of states of differing particle statistics. In this case, the reference frame is realized by a multimode state of many fermions. This latter case reproduces all of the relevant features of Ramsey interferometry, including Ramsey fringes over many repetitions of the experiment. However, the apparent inability of this proposed experiment to produce well-defined relative phases between two distinct systems each described by a coherent superposition of a boson and a fermion demonstrates that there are additional, outstanding requirements to fully lift the univalence superselection rule.

Projected Gross-Pitaevskii equation for harmonically confined Bose gases at finite temperature

We extend the projected Gross-Pitaevskii equation formalism of Davis [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.

Critical Temperature of a Trapped Bose Gas: Comparison of Theory and Experiment

We apply the projected Gross-Pitaevskii equation (PGPE) formalism to the experimental problem of the shift in critical temperature T-c of a harmonically confined Bose gas as reported in Gerbier , Phys. Rev. Lett. 92, 030405 (2004). The PGPE method includes critical fluctuations and we find the results differ from various mean-field theories, and are in best agreement with experimental data. To unequivocally observe beyond mean-field effects, however, the experimental precision must either improve by an order of magnitude, or consider more strongly interacting systems. This is the first application of a classical field method to make quantitative comparison with experiment.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

Self-consistent theory of atomic Fermi gases with a Feshbach resonance at the superfluid transition

A self-consistent theory is derived to describe the BCS-Bose-Einstein-condensate crossover for a strongly interacting Fermi gas with a Feshbach resonance. In the theory the fluctuation of the dressed molecules, consisting of both preformed Cooper pairs and bare Feshbach molecules, has been included within a self-consistent T-matrix approximation, beyond the Nozieres and Schmitt-Rink strategy considered by Ohashi and Griffin. The resulting self-consistent equations are solved numerically to investigate the normal-state properties of the crossover at various resonance widths. It is found that the superfluid transition temperature T-c increases monotonically at all widths as the effective interaction between atoms becomes more attractive. Furthermore, a residue factor Z(m) of the molecule's Green function and a complex effective mass have been determined to characterize the fraction and lifetime of Feshbach molecules at T-c. Our many-body calculations of Z(m) agree qualitatively well with recent measurments of the gas of Li-6 atoms near the broad resonance at 834 G. The crossover from narrow to broad resonances has also been studied.

Coherent molecular bound states of bosons and fermions near a Feshback resonance

We analyze molecular bound states of atomic quantum gases near a Feshbach resonance. A simple, renormalizable field theoretic model is shown to have exact solutions in the two-body sector, whose binding energy agrees well with observed experimental results in both Bosonic and Fermionic cases. These solutions, which interpolate between BEC and BCS theories, also provide a more general variational ansatz for resonant superfluidity and related problems.

Multiquantum eigenstates of a linear chain of coupled qubits

We present a technique to identify exact analytic expressions for the multiquantum eigenstates of a linear chain of coupled qubits. A choice of Hilbert subspaces is described that allows an exact solution of the stationary Schrodinger equation without imposing periodic boundary conditions and without neglecting end effects, fully including the dipole-dipole nearest-neighbor interaction between the atoms. The treatment is valid for an arbitrary coherent excitation in the atomic system, any number of atoms, any size of the chain relative to the resonant wavelength and arbitrary initial conditions of the atomic system. The procedure we develop is general enough to be adopted for the study of excitation in an arbitrary array of atoms including spin chains and one-dimensional Bose-Einstein condensates.

Ultra-cold fermions in optical lattices

We summarize recent theoretical results for the signatures of strongly correlated ultra-cold fermions in optical lattices. In particular, we focus on collective mode calculations, where a sharp decrease in collective mode frequency is predicted at the onset of the Mott metal-insulator transition; and correlation functions at finite temperature, where we employ a new exact technique that applies the stochastic gauge technique with a Gaussian operator basis.

Fringe spacing and phase of interfering matter waves

We experimentally investigate the outcoupling of atoms from Bose-Einstein condensates using two radio-frequency (rf) fields in the presence of gravity. We show that the fringe separation in the resulting interference pattern derives entirely from the energy difference between the two rf fields and not the gravitational potential difference between the two resonances. We subsequently demonstrate how the phase and polarization of the rf radiation directly control the phase of the matter wave interference and provide a semiclassical interpretation of the results.