Inhomogeneous Bose Gas

File:Superfluid to Mott insulator transition-bec3-first-expt

The physics of a BEC happens not just in momentum space, but also in position space, and it is useful to analyze it accordingly. With a trapping potential applied, the Hamiltonian is

$H=\int d^{3}r\psi ^{\dagger }(r)[{\frac {-\hbar ^{2}}{2m}}+V_{trap}]\psi (r)+{\frac {1}{2}}\int d^{3}r\int d^{3}r'\psi ^{\dagger }(r)\psi ^{\dagger }(r^{\prime })U(r-r^{\prime })\psi (r^{\prime })\psi (r)$

This must be approximated, in the spirit of Bogolubov's momentum space approximation, to obtain a useful solution. We thus replace

${\hat {\psi }}(r,t)=\psi (r,t)+{\tilde {\psi }}(r,t)\,,$

where $\psi (r,t)$ is an expectation, and ${\tilde {\psi }}(r,t)$ captures the quantum (+ thermal) fluctuations. We further assume that the interaction potential is a delta function, and the resulting equation is a nonlinear Schrodinger equation (also known as a Gross-Pitaevskii equation):

$i\hbar {\frac {\partial \psi }{\partial t}}=\left[{-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V_{trap}+U_{0}N|\psi (r,t)|^{2}}\right]\psi (r,t)$

The $|\psi (r,t)|^{2}$ term adds an energy proportional to the density due to the interactions. That density-dependence implies that increasing the density comes with an energy cost and decreasing density lowers the energy. Since the density is multiplied by the wave function, there is a non-linear energy dependence on the density. Therefore, lowering the density in one region and raising it in another costs energy and so this term makes the condensate try to have a uniform density distribution (the total number of particles is fixed so the integrated density cannot change).

The density-dependent term has an interesting interplay with the other two terms. Consider the effect on a condensate in a harmonic trap for instance. If it were not for the interactions, the BEC density distribution would simply be the Gaussian ground state of the harmonic oscillator. However, the central region of that distribution is very dense, so the density-dependent term will try to flatten out the distribution and push atoms out towards the wings. It can flatten out the distribution near the center very well, but the trap potential makes it impossible to push atoms out too far, so the atom distribution stays roughly flat over a finite region set be the interaction energy and trap potential. At the edge of the flat region of the density distribution is a crossover to the wings of the distribution. The density cost makes the cloud want to make this crossover region as small as possible. However, making it very small would take a lot of curvature in the wave function, causing the kinetic energy term to increase. Therefore the distribution at the edges reaches a compromise between interaction energy and kinetic energy curvature. The length scale of the crossover region set by this compromise (assuming repulsive interactions) is known as the healing length $\xi$,

$\xi =(8\pi an)^{-1/2}\,,$

arising from

${\frac {\hbar ^{2}}{2m\xi ^{2}}}=U_{0}h\,.$

File:Superfluid to Mott insulator transition-bec3-healing

If the interactions are really strong, the kinetic energy term can be neglected, because the interactions will keep the density constant in its spatial distribution. Such an approximation is the Thomas-Fermi approximation, giving an equation for the wavefunction,

$[V_{ext}+U_{0}|\psi (r,t)|^{2}-\mu ]\psi (r)=0\,,$

giving the solution

$|\psi (r,t)|^{2}={\frac {1}{U_{0}}}(\mu -V(r)\,.$

The wavefunction is essentially just the potential filled up to the chemical potential level, inverted. For a quadaratic potential, $V(r)=m\omega ^{2}r^{2}/2$, the chemical potential is

$\mu ={\frac {\hbar \omega }{2}}(15N_{a}/a_{osc})^{2/5}\,,$

where $x=15N_{a}/a_{osc}$ is a common term worth identifying, and $a_{osc}={\sqrt {\hbar ^{2}/m\omega }}$ is a characteristic length scale of the oscillator, its zero point motion. Defining $R=a_{osc}x^{1/5}$, we may find $\mu =m\omega ^{2}R^{2}/2$. This explains the profile of the condensate data obtained in experiments:

File:Superfluid to Mott insulator transition-bec3-bec-peak

\noindent Note that the size of the ground stat BEC is much larger than the zero-point motion of the harmonic oscillator. This is due to the pressure of the repulsive interactions. The Gross-Pitaevskii interaction gives not only the ground state wavefunction, but also the dynamics of the system. For example, it predicts soliton formation: stable wavefunctions with a size scale determined by a balance of the kinetic energy and the internal interactions. This requires, however, an attractive potential. Such soliton formation can nevertheless be seen in BEC's, with tight traps (see recent Paris experiments).

Length and energy scales in BEC

\begin{itemize}

• Size of atom: $a=3$ nm
• Separation between atoms $n^{1/3}=200$ nm
• Matter wavelength $\lambda _{dB}=1$ $\mu$m
• Size of confinement $a_{osc}=30$ $\mu$m

Note that

$a\ll n^{1/3}\leq \lambda _{dB}<2\pi \xi <a_{osc}\,.$

For a gas, $a\ll n^{1/3}$. For a BEC, $\lambda _{dB}\geq n^{1/3}$ in addition. The corresponding energy scales are also useful to identify. Let $\hbar \omega =a_{osc}$. Then:

$k_{B}T_{s-wave}\gg k_{B}T_{c}\geq k_{B}T>U_{int}>\hbar \omega$

The interaction energy scale $U_{int}\sim (h^{2}/m)na$, corresponding to the healing length.

Back to: Quantum gases