1. Lecture 5: Mean-field theory and the Gross-Pitaevskii equation
In the previous lecture we saw two applications of the variational principle to the class of product states and fermionic gaussian states, respectively. In both cases we obtained an effective equation for the ground-state properties involving only a single effective particle (a single spin in the first case and a single majorana fermion in the second case).
In this lecture we continue our study of mean-field theory in the bosonic setting, in particular to the description of Bose-Einstein condensates. Here we find a similar result: we’ll obtain a nonlinear effective equation for the condensate in terms of a single effective particle degree of freedom. Before we do this we need to review some of the formalism for the description of quantum fields and coherent states.
The pdf version of the notes can be found here.
2. Quantum fields
Recall that in lecture 3 we showed that the second-quantised many particle Schrödinger equation for the abstract state vector admitted the form
where the operator is given by
where and are the annihilation and creation operators, respectively, for a complete set of single-particle wave functions . (We consider the bosonic and fermionic cases in parallel, so that and either satisfy the CCR or CAR, respectively.)
It is convenient to introduce the following field operators
These operators satisfy the simple (anti)commutation relations, following from the completeness relation:
where the subscript on the brackets denotes commutator/anticommutator depending on whether the system is bosonic/fermionic, respectively.
In these terms the hamiltonian operator may be rewritten as
One example that will concern us in the sequel is that of a dilute gas of bosons dominated by low energy two-particle -wave collisions moving in an external potential. In this case we have that is given by , and , where , so that
3. Coherent states
We now recall the definition of a coherent state. Suppose that and are the annihilation and creation operators, respectively, for a single quantum harmonic oscillator satisfying , etc. Let be the vacuum state which satisfies
Let and define the displacement operator
This is a unitary operator. One can use the Baker-Campbell-Hausdorff expansion to show that
A coherent state is then defined to be
These states have the property that
The set of coherent states supply an overcomplete basis for hilbert space; we have the completeness relation
Any coherent state is an example of a gaussian state. However the set of gaussian states is larger than that of coherent states: there exist squeezed states which are also gaussian and which aren’t coherent.
The generalisation to many harmonic oscillators proceeds as follows. Let and , , be the annihilation and creation operators, respectively, for a quantum harmonic oscillators satisfying , . We think of the harmonic oscillators as being arranged in a line.
Let and define
The set of states forms a convenient variational class. To see what a variational calculation looks like in this setting we study the following
4. Example 1: the Bose-Hubbard model
The Bose-Hubbard model is an effective model for bosons hopping in a lattice with repulsive onsite interactions. It can be taken as a model for, e.g., trapped cold atoms moving in the presence of a periodic potential. The hamiltonian is given by
The expectation value of in a coherent state is given by
Taking the limit and assuming translation invariance gives the following expression for the energy density
The variation is then over the single complex number :
The minimum occurs at and
This solution reflects the fact that if the chemical potential is large enough then it is energetically favourable for the system to support a finite density of particles. Unfortunately the translation-invariant version of our class which we used here is a little too trivial to give good information about the Bose-Hubbard model. The non-translation invariant case is more interesting and gives a better approximation. However, we can no longer find the minimum in compact form.
5. Coherent states for quantum fields
Let and be the bosonic field annihilation and creation operators. We define the quantum field displacement operator via
where is a complex function. We use the field displacement operator to define a field coherent state via
Note that a field coherent state obeys
etc. Exercise: prove this using a Taylor series for the exponential and the commutation relations for the field operators. We now take as our variational class the set
where the set of functions in the definition of are taken to be twice-differentiable.
6. Example 2: Bose-Einstein condensates
We apply the variational principle to given in (6) modelling a dilute gas of bosons using the class ; we obtain the expectation value
To apply the variational principle we need to extremise by minimising over all . This is easy to carry out using variational calculus: the Euler-Lagrange equations give
which reduce to
This nonlinear Schrödinger equation is known as the (time-independent) Gross-Pitaevskii equation (GPE). The single-particle solution models a Bose-Einstein condensate in an external potential. The nonlinear term models the averaged, or mean-field, interaction with the other particles.
Exercise: make the Thomas-Fermi approximation in the large- limit by neglecting the kinetic energy derivative term in the GPE. What solution do you obtain for ? Compare this with the lattice calculation of the previous section.