**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*

and

These operators satisfy the simple (anti)commutation relations, following from the completeness relation:

,

and

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

where

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

and

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.