In this chapter we present several dynamical simulations that make use of the numerical methods discussed in the previous chapter. Some of these are model simulations that are not directly linked to experiment but are designed to investigate some aspect of the dynamical behaviour. Others are performed with the express purpose of explaining specific experimental data. Both kinds of simulation serve to illustrate the range of nonequilibrium phenomena that can be studied in ultracold Bose gases using the ZNG equations.

A dynamical simulation is typically initiated in one of two ways. Either an appropriate nonequilibrium initial state is imposed on the system, or the system, initially in equilibrium, is dynamically excited by the application of an external perturbation. The latter parallels the procedure used experimentally to study small-amplitude collective excitations and usually amounts to some parametric modulation of the trapping potential. However, this approach may not always be feasible if the excitation phase requires a prohibitively long simulation time. In this case, the best one can do is to specify some initial nonequilibrium state that represents the experimental situation as closely as possible. This is not an issue in model simulations, where we are at liberty to specify the initial state in whatever way serves our purpose.

In Sections 12.1–12.3 we present three examples of model simulations. All essentially check some aspect of the numerical procedures. By studying the equilibration of an initial nonequilibrium state in Section 12.1, we confirm that the total number of atoms is conserved to a very good approximation during the course of the evolution. This is a nontrivial result since, as explained at the end of subsection 11.3.2, the numbers of condensate and thermal atoms change in quite different ways.

In contrast with Chapter 2, in this chapter we include the dynamics of the thermal cloud. As noted in Chapter 1, we treat the noncondensate atoms using the simplest microscopic model approximation that captures the important physics. In particular, we consider only temperatures high enough (T ≥ 0.4TBEC) that the noncondensate atoms can be described by a particle-like Hartree–Fock (HF) spectrum. To extend the analysis to very low temperatures is in principle straightforward (see Chapter 7). However, the details are more complicated since the excitations of the thermal cloud take on a collective aspect (i.e. a Bogoliubov-type quasiparticle spectrum must be used). In trapped Bose gases, the HF single-particle spectrum gives a good approximation down to much lower temperatures than in the case of uniform Bose gases, as first emphasized by Giorgini et al. (1997).

In Section 3.1, we derive a generalized form of the Gross–Pitaevskii equation for the Bose order parameter Φ(r, t) that is valid at finite temperatures. It involves terms that are coupled to the noncondensate component (the thermal cloud) and thus its solution in general requires one to know the equations of motion for the dynamics of the noncondensate atoms. In Section 3.2 we restrict ourselves to finite temperatures high enough that the noncondensate atoms can be described by a quantum kinetic equation for the single-particle distribution function f(p, r, t). A detailed microscopic derivation of this kind of kinetic equation is given in Chapters 6 and 7 using the Kadanoff–Baym Green's function formalism.

A characteristic feature of a Bose-condensed gas is that the kinetic equation governing f(p, r, t) involves a collision integral C12[f, Φ] describing collisions between condensate and noncondensate atoms.

The time-dependent behaviour of Bose–Einstein condensed clouds, such as collective modes and the expansion of a cloud when released from a trap, is an important source of information about the physical nature of the condensate. In addition, the spectrum of elementary excitations of the condensate is an essential ingredient in calculations of thermodynamic properties. In this chapter we treat the dynamics of a condensate at zero temperature starting from a time-dependent generalization of the Gross–Pitaevskii equation used in Chapter 6 to describe static properties. From this equation one may derive equations very similar to those of classical hydrodynamics, which we shall use to calculate properties of collective modes.

We begin in Sec. 7.1 by describing the time-dependent Gross–Pitaevskii equation and deriving the hydrodynamic equations, which we then use to determine the excitation spectrum of a homogeneous Bose gas (Sec. 7.2). Subsequently, we consider modes in trapped clouds (Sec. 7.3) within the hydrodynamic approach, and also describe the method of collective coordinates and the related variational method. In Sec. 7.4 we consider surface modes of oscillation, which resemble gravity waves on a liquid surface. The variational approach is used in Sec. 7.5 to treat the free expansion of a condensate upon release from a trap. Finally, in Sec. 7.6 we discuss solitons, which are exact one-dimensional non-linear solutions of the time-dependent Gross–Pitaevskii equation.

In this chapter we consider selected topics in the theory of trapped gases at non-zero temperature when the effects of interactions are taken into account. The task is to extend the considerations of Chapters 8 and 10 to allow for the trapping potential. In Sec. 11.1 we begin by discussing energy scales, and then calculate the transition temperature and thermodynamic properties. We show that at temperatures of the order of Tc the effect of interactions on thermodynamic properties of clouds in a harmonic trap is determined by the dimensionless parameter N1/6a/ā. Here ā, which is defined in Eq. (6.24), is the geometric mean of the oscillator lengths for the three principal axes of the trap. Generally this quantity is small, and therefore under many circumstances the effects of interactions are small. At low temperatures, thermodynamic properties may be evaluated in terms of the spectrum of elementary excitations of the cloud, which we considered in Secs. 7.2, 7.3, and 8.2. At higher temperatures it is necessary to take into account thermal depletion of the condensate, and useful approximations for thermodynamic functions may be obtained using the Hartree–Fock theory as a starting point.

The remainder of the chapter is devoted to non-equilibrium phenomena. As we have seen in Secs. 10.3–10.5, two ingredients in the description of collective modes and other non-equilibrium properties of uniform gases are the two-component nature of condensed Bose systems, and collisions between excitations.

The topic of Bose–Einstein condensation in a uniform, non-interacting gas of bosons is treated in most textbooks on statistical mechanics. In the present chapter we discuss the properties of a non-interacting Bose gas in a trap. We shall calculate equilibrium properties of systems in a semi-classical approximation, in which the energy spectrum is treated as a continuum. For this approach to be valid the temperature must be large compared with Δ∊/k, where Δ∊ denotes the separation between neighbouring energy levels. As is well known, at temperatures below the Bose–Einstein condensation temperature, the lowest-energy state is not properly accounted for if one simply replaces sums by integrals, and it must be included explicitly.

The statistical distribution function is discussed in Sec. 2.1, as is the single-particle density of states, which is a key ingredient in the calculations of thermodynamic properties. Calculations of the transition temperature and the fraction of particles in the condensate are described in Sec. 2.2. In Sec. 2.3 the semi-classical distribution function is introduced. From this we obtain the density profile and the velocity distribution of particles, and use these to determine the shape of an anisotropic cloud after free expansion. Thermodynamic properties of Bose gases are calculated as functions of the temperature in Sec. 2.4. Finally, Sec. 2.5 treats corrections to the transition temperature due to a finite particle number.

The laser cooling mechanisms described in Chapter 4 operate irrespective of the statistics of the atom, and they can therefore be used to cool Fermi species. The statistics of a neutral atom is determined by the number of neutrons in the nucleus, which must be odd for a fermionic atom. Since alkali atoms have odd atomic number Z, their fermionic isotopes have even mass number A. Such isotopes are relatively less abundant than those with odd A since they have both an unpaired neutron and an unpaired proton, which increases their energy by virtue of the odd–even effect. In early experiments, 40K and 6Li atoms were cooled to about one-quarter of the Fermi temperature. More recently, fermionic alkali atoms have been cooled to temperatures well below one-tenth of their Fermi temperature, and in addition a degenerate gas of the fermionic species 173Yb (with I = 5/2) has been prepared.

In the classical limit, at low densities and/or high temperatures, clouds of fermions and bosons behave alike. The factor governing the importance of quantum degeneracy is the phase-space density ϖ introduced in Eq. (2.24), and in the classical limit ϖ ≪ 1. When ϖ becomes comparable with unity, gases become degenerate: bosons condense in the lowest single-particle state, while fermions tend towards a state with a filled Fermi sea. As one would expect on dimensional grounds, the degeneracy temperature for fermions – the Fermi temperature TF – is given by the same expression as the Bose–Einstein transition temperature for bosons, apart from a numerical factor of order unity.

The electric field intensity of a standing-wave laser field is periodic in space. Due to the ac Stark effect, this gives rise to a spatially periodic potential acting on an atom, as explained in Chapter 4 (see, e.g., Eq. (4.31)). This is the physical principle behind the generation of optical lattices. By superimposing a number of different laser beams it is possible to generate potentials which are periodic in one, two or three dimensions. The suggestion that standing light waves may be used to confine the motion of atoms dates back to 1968 and is due to Letokhov. The first experimental realization of an optical lattice was achieved in 1987 for a classical gas of cesium atoms.

The study of atoms in such potentials has many different facets. At the simplest level, it is possible to study the energy band structure of atoms moving in these potentials and to explore experimentally a number of effects that are difficult to observe for electrons in the periodic lattice of a solid. Interactions between atoms introduce qualitatively new effects. Within mean-field theory, which applies when the number of atoms in the vicinity of a single minimum of the potential is sufficiently large, one finds that interatomic interactions give rise to novel features in the band structure. These include multivaluedness of the energy for a given band and states possessing a periodicity different from that of the optical lattice.

The experimental discovery of Bose–Einstein condensation in trapped atomic clouds opened up the exploration of quantum phenomena in a qualitatively new regime. Our aim in the present work is to provide an introduction to this rapidly developing field.

The study of Bose–Einstein condensation in dilute gases draws on many different subfields of physics. Atomic physics provides the basic methods for creating and manipulating these systems, and the physical data required to characterize them. Because interactions between atoms play a key role in the behaviour of ultracold atomic clouds, concepts and methods from condensed matter physics are used extensively. Investigations of spatial and temporal correlations of particles provide links to quantum optics, where related studies have been made for photons. Trapped atomic clouds have some similarities to atomic nuclei, and insights from nuclear physics have been helpful in understanding their properties.

In presenting this diverse range of topics we have attempted to explain physical phenomena in terms of basic principles. In order to make the presentation self-contained, while keeping the length of the book within reasonable bounds, we have been forced to select some subjects and omit others. For similar reasons and because there now exist review articles with extensive bibliographies, the lists of references following each chapter are far from exhaustive.

This book originated in a set of lecture notes written for a graduate-level one-semester course on Bose–Einstein condensation at the University of Copenhagen.

In preceding chapters we have explored properties of Bose–Einstein condensates with a single macroscopically occupied quantum state, and spin degrees of freedom of the atoms were assumed to play no role. In the present chapter we extend the theory to systems in which two or more quantum states are macroscopically occupied.

The simplest example of such a multi-component system is a mixture of two different species of bosons, for example two isotopes of the same element, or two different atoms. The theory of such systems can be developed along the same lines as that for one-component systems developed in earlier chapters, and we do this in Sec. 12.1.

Since alkali atoms have spin, it is also possible to make mixtures of the same isotope, but in different internal spin states. This was first done experimentally by the JILA group, who made a mixture of 87Rb atoms in hyperfine states F = 2, mF = 2 and F = 1, mF = −1. Mixtures of hyperfine states of the same isotope differ from mixtures of distinct isotopes because atoms can undergo transitions between hyperfine states, while transitions that convert one isotope into another may be neglected. Transitions between different hyperfine states can influence equilibrium properties markedly if the interaction energy per particle is comparable with or larger than the energy difference between hyperfine levels. In magnetic traps it is difficult to achieve such conditions, since the trapping potential depends on the particular hyperfine state.

The ability to vary the force constants of a trapping potential makes it possible to create very elongated or highly flattened clouds of atoms. This opens up the study of Bose–Einstein condensation in lower dimensions, since motion in one or more directions may then effectively be frozen out at sufficiently low temperatures. In a homogeneous system Bose–Einstein condensation cannot take place at non-zero temperature in one or two dimensions, but in traps the situation is different because the trapping potential changes the energy dependence of the density of states. This introduces a wealth of new phenomena associated with lower dimensions which have been explored both theoretically and experimentally. A general review may be found in the lecture notes Ref.

For a system in thermal equilibrium, the condition for motion in a particular direction to be frozen out is that the energy difference between the ground state and the lowest excited state for the motion must be much greater than the thermal energy kT. This energy difference is ħωi if interactions are unimportant for the motion in the i direction. If the interaction energy nU0 is large compared with ħωi and the trap is harmonic, the lowest excited state is a sound mode with wavelength comparable to the spatial extent of the cloud in the i direction.

Bose–Einstein condensates of particles behave in many ways like coherent radiation fields, and the realization of Bose–Einstein condensation in dilute gases has opened up the experimental study of interactions between coherent matter waves. In addition, the existence of these dilute trapped quantum gases has prompted a re-examination of a number of theoretical issues.

In Sec. 13.1 we discuss Josephson tunnelling of a condensate between two wells and the role of fluctuations in particle number and phase. The number and phase variables play a key role in the description of the classic interference experiment, in which two clouds of atoms are allowed to expand and overlap (Sec. 13.2). Rather surprisingly, an interference pattern is produced even though initially the two clouds are completely isolated. Density fluctuations in a Bose gas are studied in Sec. 13.3, where we relate atomic clock shifts to the two-particle correlation function. The ability to manipulate gaseous Bose–Einstein condensates by lasers has made possible the study of coherent matter wave optics and in Sec. 13.4 we describe applications of these techniques to observe solitons, Bragg scattering, and non-linear mixing of matter waves. How to characterize Bose–Einstein condensation in terms of the density matrix is the subject of Sec. 13.5, where we also consider fragmented condensates.

An important feature of cold atomic vapours is that particle separations, which are typically of order 102 nm, are usually an order of magnitude larger than the length scales associated with the atom–atom interaction. Consequently, the two-body interaction between atoms dominates, and three-and higher-body interactions are unimportant. Moreover, since the atoms have low velocities, many properties of these systems may be calculated in terms of a single parameter, the scattering length.

An alkali atom in its electronic ground state has several different hyperfine states, as we have seen in Secs. 3.1 and 3.2. Interatomic interactions give rise to transitions between these states and, as we described in Sec. 4.6, such processes are a major mechanism for loss of trapped atoms. In a scattering process, the internal states of the particles in the initial or final states are described by a set of quantum numbers, such as those for the spin, the atomic species, and their state of excitation. We shall refer to a possible choice of these quantum numbers as a channel. At the temperatures of interest for Bose–Einstein condensation, atoms are in their electronic ground states, and the only relevant internal states are therefore the hyperfine states. Because of the existence of several hyperfine states for a single atom, the scattering of cold alkali atoms is a multi-channel problem.

The advent of the laser opened the way to the development of powerful methods for producing and manipulating cold atomic gases. To set the stage we describe a typical experiment, which is shown schematically in Fig. 4.1. A beam of sodium atoms emerges from an oven at a temperature of about 600 K, corresponding to a speed of about 800 m s−1, and is then passed through a so-called Zeeman slower, in which the velocity of the atoms is reduced to about 30 m s−1, corresponding to a temperature of about 1 K. In the Zeeman slower, a laser beam propagates in the direction opposite that of the atomic beam, and the radiation force produced by absorption of photons retards the atoms. Due to the Doppler effect, the frequency of the atomic transition in the laboratory frame is not generally constant, since the atomic velocity varies. However, by applying an inhomogeneous magnetic field designed so that the Doppler and Zeeman effects cancel, the frequency of the transition in the rest frame of the atom may be held fixed. On emerging from the Zeeman slower the atoms are slow enough to be captured by a magneto-optical trap (MOT), where they are further cooled by interactions with laser light to temperatures of order 100 μK. Another way of compensating for the changing Doppler shift is to increase the laser frequency in time, which is referred to as ‘chirping’.

A number of atomic properties play a key role in experiments on cold atomic gases, and we shall discuss them briefly in the present chapter with particular emphasis on alkali atoms. Basic atomic structure is the subject of Sec. 3.1. Two effects exploited to trap and cool atoms are the influence of a magnetic field on atomic energy levels, and the response of an atom to radiation. In Sec. 3.2 we describe the combined influence of the hyperfine interaction and the Zeeman effect on the energy levels of an atom, and in Sec. 3.3 we review the calculation of the atomic polarizability. In Sec. 3.4 we summarize and compare some energy scales.

Atomic structure

The total spin of a Bose particle must be an integer, and therefore a boson made up of fermions must contain an even number of them. Neutral atoms contain equal numbers of electrons and protons, and therefore the statistics that an atom obeys is determined solely by the number of neutrons N: if N is even, the atom is a boson, and if it is odd, a fermion. Since the alkalis have odd atomic number Z, boson alkali atoms have odd mass numbers A. Likewise for atoms with even Z, bosonic isotopes have even A. In Table 3.1 we list N, Z, and the nuclear spin quantum number I for some alkali atoms and hydrogen.