Atomic physics (2005)
.pdf
236 Magnetic trapping, evaporative cooling and Bose–Einstein condensation
Hence, the number density of atoms n (r) = N |ψ (r)|2 in the harmonic potential has the form of an inverted parabola:
n (r) = n0 |
1 |
|
x2 |
|
|
y2 |
|
z2 |
(10.38) |
||
|
|
|
|
|
|
|
|
, |
|||
− Rx2 |
− Ry2 |
|
|||||||||
|
|
− Rz2 |
|
||||||||
where n0, the number density at the centre of the condensate, is |
|||||||||||
n0 = |
Nµ |
. |
|
|
|
(10.39) |
|||||
g |
|
|
|
||||||||
The condensate has an ellipsoidal shape and the density goes to zero at points on the axes given by x = ±Rx, y = ±Ry and z = ±Rz , defined by
1 |
M ω2R2 = µ , |
(10.40) |
2 x x
and similarly for Ry and Rz . In the Thomas–Fermi regime, the atoms fill up the trap to the level of the chemical potential as illustrated in Fig. 10.11, just like water in a trough. The chemical potential µ is determined by the normalisation condition
1 = |ψ|2 |
dx dy dz = |
µ 8π |
(10.41) |
|||||||||
|
|
|
|
RxRy Rz . |
||||||||
g |
15 |
|||||||||||
A useful form for µ is |
|
|
|
|
|
|
|
|
|
|
||
µ = |
|
× |
1 |
|
15Na |
2/5 . |
(10.42) |
|||||
ω |
||||||||||||
2 |
|
|||||||||||
|
|
aho |
|
|
|
|||||||
The mean oscillation frequency is defined as ω = (ωxωy ωz )1/3 and aho is calculated using this frequency in eqn 10.34. Typical values for an Io e trap are given in the following table together with the important properties of a BEC of sodium, calculated from the key formulae in eqns 10.42, 10.40 and 10.39 (in that order).
Scattering length (for Na) |
|
|
|
|
|
a |
2.9 nm |
||
Radial oscillation frequency (ωx = ωy ) |
ωx/2π |
250 Hz |
|||||||
Axial oscillation frequency |
|
|
ωz /2π |
16 Hz |
|||||
Average oscillation frequency |
|
|
|
|
|
|
|
100 Hz |
|
|
|
|
ω/2π |
||||||
Zero-point energy (in temperature units) |
1 |
|
|
|
2.4 nK |
||||
ω/kB |
|||||||||
Harmonic oscillator length (for |
|
|
) |
2 |
aho |
2.1 µm |
|||
ω |
|
||||||||
Number of atoms in condensate |
|
|
|
N0 |
106 |
||||
Chemical potential |
|
|
|
µ |
130 nK |
||||
Radial size of condensate |
|
|
Rx = Ry |
15 µm |
|||||
Axial size of condensate |
|
|
|
Rz |
95 µm |
||||
Peak density of condensate |
6 |
|
|
|
n0 |
2 × 1014 cm−3 |
|||
Critical temperature (for 4 × 10 |
|
atoms) |
|
TC |
760 nK |
||||
Critical density (for 4 × 106 atoms) |
|
nC |
4 × 1013 cm−3 |
||||||
10.6 A Bose–Einstein condensate 237
(a)
(b)
Fig. 10.11 In the Thomas–Fermi regime the condensate has the same shape as the confining potential. (a) A harmonic potential. (b) The density of the atoms in a harmonic trap has an inverted-parabolic shape (along all three axes).
The critical temperature and density at the onset of Bose condensation are not properties of the condensate, but are calculated for a cloud of 4 × 106 atoms (from eqn 10.19); this would lead to a condensate of roughly N0 106 atoms after evaporative cooling to T /TC 0.5 (where most atoms are in the condensate, see eqn F.16). Both µ and TC have a weak dependence on N (eqns 10.19 and 10.42) and a similar (but not the same) dependence on ω, so these quantities have a similar relative magnitude for many cases. Note, however, that µ depends on
(a) |
(b) |
(c) |
Fig. 10.12 This sequence of images shows a Bose–Einstein condensate being born out of a cloud of evaporatively-cooled atoms in a magnetic trap. Each image was taken after a time-of-flight expansion. The cloud changes its size and shape as it undergoes a phase transition: (a) a thermal cloud just above the critical temperature TC has a spherical shape (isotropic expansion); (b) a cloud of atoms at 0.9 TC has a Bose-condensed fraction in the centre surrounded by a halo of thermal atoms; and (c) well below the critical temperature (< 0.5 TC) most of the atoms are in the condensate (lowest energy state of the trap). These images come from a system that does not have the same aspect ratio as an Io e trap, but they illustrate the anisotropic expansion of the condensate wavefunction. See Fig. 10.13 for dimensions and other details. Data provided by Nathan Smith, Physics department, University of Oxford.
240 Magnetic trapping, evaporative cooling and Bose–Einstein condensation
travel in the gas has great significance for superfluidity. For motion slower than this speed the condensate flows smoothly around obstacles without exciting any particles out of the ground state of the quantum gases. This type of flow does not dissipate any energy and so it is frictionless and the gas is superfluid.
10.7.2Healing length
The Thomas–Fermi approximation neglects the kinetic energy term in the Schr¨odinger equation. This leads to a physically unrealistic sharp edge at the surface of the condensate (see Fig. 10.11)—such a discontinuity in the gradient would make 2ψ infinite. Therefore we have to take kinetic energy into account at the boundary. To determine the shortest distance ξ over which the wavefunction can change we equate the kinetic term (that contains 2ψ 2/(2M ξ2)) to the energy scale of the system given by the chemical potential. Atoms with energy higher than µ leave the condensate. Using n0 = Nµ/g (from eqn 10.39) and eqn 10.30 for g, we find that
2 |
gn0 |
= |
4π 2an0 |
|
(10.45) |
|
|
µ = |
|
|
. |
||
2M ξ2 |
N |
M |
||||
√
Hence ξ = 1/ 8πan0, e.g. ξ = 0.3 µm for a sodium condensate with n0 = 2 × 1014 cm−3. Typically, ξ Rx and smoothing of the wavefunction only occurs in a thin boundary layer, and these surface e ects give only small corrections to results calculated using the Thomas–Fermi approximation. This so-called healing length also determines the size of the vortices that form in a superfluid when the confining potential rotates (or a fast moving object passes through it). In these little ‘whirlpools’ the wavefunction goes to zero at the centre, and ξ determines the distance over which the density rises back up to the value in the bulk of the condensate, i.e. this healing length is the distance over which the superfluid recovers from a sharp change.
10.7.3The coherence of a Bose–Einstein condensate
Figure 10.14 shows the result of a remarkable experiment carried out by the group led by Wolfgang Ketterle at MIT. They created two separate condensates of sodium at the same time. After the trapping potential was turned o the repulsion between the atoms caused the two clouds to expand and overlap with each other (as in the time-of-flight technique used to observe Bose–Einstein condensation, see Fig. 10.12). The two condensates interfere to give the fringes shown in the figure; there are no atoms at certain positions where the matter waves from the two sources interfere destructively—these atoms do not disappear, but they are redistributed to positions in the fringe pattern where the matter waves add constructively. Such interference is well known in optics; however, there is a very interesting di erence between this experiment
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Exercises for Chapter 10 |
|
|
245 |
||||||
|
(b) In thermodynamics the chemical potential is |
approximations leads to the same zeroth- |
|||||||||||||||||||||||||||||||
|
the energy required to remove a particle from |
order approximation for the chemical poten- |
|||||||||||||||||||||||||||||||
|
the system µ = ∂E/∂ |
N |
, where E is the total |
tial given above and |
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
energy of the system. Show that E = 7 Nµ. |
i |
d |
(δψ (t)) |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
(10.8) |
Expansion of a non-interacting condensate |
dt |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
Although experiments are not carried out with |
|
|
2k2 |
|
|
|
|
2 |
|
2 |
|
|
|
|||||||||||||||||||
|
a non-interacting gas it is instructive to consider |
= |
2M |
δψ (t) + g |ψ0| 2δψ (t) + gψ0 δψ |
|
|
(t) . |
||||||||||||||||||||||||||
|
what happens when a = 0. In this case the con- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(10.47) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
densate has the same size as the ground state |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
of a quantum harmonic oscillator (for any N0) |
(b) Show that equating terms with the same time |
|||||||||||||||||||||||||||||||
|
dependence leads to two coupled equations |
||||||||||||||||||||||||||||||||
|
and the initial momentum along each direction |
||||||||||||||||||||||||||||||||
|
can be estimated from the uncertainty principle. |
for u and v that, in matrix form, are |
|
|
|
|
|||||||||||||||||||||||||||
|
For atoms of the same mass as sodium (but with |
k + 2g ψ0 2 |
|
µ |
|
gψ02 |
|
|
|
|
u |
||||||||||||||||||||||
|
a = 0) released from a trap with a radial os- |
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
g(ψ| |
0 )2| |
− |
|
k + 2g |ψ0|2 − µ v |
||||||||||||||||||||||||
|
cillation frequency of 250 Hz and a frequency of |
|
|
|
|
||||||||||||||||||||||||||||
|
16 Hz for axial motion, estimate roughly the time |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u |
||||||||||||||
|
of flight at which the cloud is spherical. |
|
|
|
|
|
|
|
|
|
|
|
|
|
= ω −v , |
||||||||||||||||||
(10.9) |
Excitations of a Bose condensate |
|
|
where k = 2k2/2M. |
|
|
|
|
|
|
|||||||||||||||||||||||
|
The vibrational modes of a condensate can be |
(c) Hence show that u and v are solutions of the |
|||||||||||||||||||||||||||||||
|
viewed as compression waves that form a stand- |
||||||||||||||||||||||||||||||||
|
matrix equation |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
ing wave within the condensate; |
|
hence these |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
k + µ − ω |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
modes have frequencies of the order of the speed |
|
|
|
gψ02 |
u |
= 0 . |
||||||||||||||||||||||||||
|
of sound divided by the size of the condensate |
|
g(ψ0 )2 |
|
k + µ + ω v |
|
|
|
|
||||||||||||||||||||||||
|
vs/R. Show that this collective motion of the |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
condensate occurs at a comparable frequency to |
From the determinant of this matrix, find |
|||||||||||||||||||||||||||||||
|
the oscillation of individual atoms in the mag- |
the relation between the angular frequency of |
|||||||||||||||||||||||||||||||
|
netic trap. |
|
|
|
|
|
|
|
|
|
the small oscillations ω and the magnitude of |
||||||||||||||||||||||
(10.10) |
Derivation of the speed of sound |
|
|
their wavevector k (the dispersion relation). |
|||||||||||||||||||||||||||||
|
|
Show that for low energies this gives the same |
|||||||||||||||||||||||||||||||
|
The time-dependent Schr¨odinger equation for the |
||||||||||||||||||||||||||||||||
|
expression for the speed of sound ω/k found |
||||||||||||||||||||||||||||||||
|
wavefunction of an atom in a Bose–Einstein con- |
||||||||||||||||||||||||||||||||
|
in Section 10.7.1.50 |
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
densate in a uniform potential is |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
dψ |
|
|
2 |
2 |
ψ + g |ψ| |
2 |
|
(10.11) Attractive interactions |
|
|
|
|
|
|
|
||||||||||||||
|
|
i |
|
= |
− |
|
|
|
|
ψ , |
In certain hyperfine states, the scattering length |
||||||||||||||||||||||
|
|
dt |
2M |
|
|
||||||||||||||||||||||||||||
|
where, for simplicity, the potential has been taken |
a of alkali metal atoms changes with the applied |
|||||||||||||||||||||||||||||||
|
magnetic field and this feature has been used to |
||||||||||||||||||||||||||||||||
|
as zero (V = 0). The wavefunction ψ = ψ0e−iµt/ |
||||||||||||||||||||||||||||||||
|
satisfies this equation with a chemical potential |
perform experiments in which the atoms have at- |
|||||||||||||||||||||||||||||||
|
tractive interactions a < 0. |
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
µ = g |ψ0|2 . |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
Show that eqn 10.33 can be written in the form |
||||||||||||||||||||||||
|
The trial wavefunction with small fluctuations |
|
|
|
|
4 |
|
E |
= x−2 |
+ x2 + Gx−3 . |
|
|
|
|
|
||||||||||||||||||
|
can be written as |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
3 ω |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
ψ = |
ψ |
+ uei(kx−ωt) + v e−i(kx−ωt) e−iµt/ |
By plotting graphs for various values of the pa- |
|||||||||||||||||||||||||||||
|
5 |
0 |
−iµt/ |
+ δψ (t) , |
|
|
|
6 |
rameter G between 0 and −1, estimate the lowest |
||||||||||||||||||||||||
|
= ψ0e |
|
|
|
|
|
|
|
value of G for which there exists a minimum in |
||||||||||||||||||||||||
|
where the amplitudes |u| and |v| are small com- |
the energy as a function of x. For an atomic |
|||||||||||||||||||||||||||||||
|
species with a scattering length of a = |
− |
5 nm |
||||||||||||||||||||||||||||||
|
pared to ψ . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
| |
0| |
|
|
|
|
|
|
|
|
|
|
in a trap where aho = 2 µm, estimate the max- |
|||||||||||||||||||
|
(a) Show that substituting this function into the |
imum number of atoms that a Bose condensate |
|||||||||||||||||||||||||||||||
|
Schr¨odinger |
|
equation |
and making suitable |
can contain without collapsing. |
|
|
|
|
|
|
||||||||||||||||||||||
50After problem devised by Professor Keith Burnett, Physics graduate class, University of Oxford.
