This post is based on my paper with Avy Soffer.

Let us recall from Chapter 1 that the BEC occupies the lowest quantum state, at that point macroscopic quantum phenomena become apparent. Above the BEC are excited atoms, which occupy higher quantum states. The are two types of interactions:

- excited atoms – excited atoms,
- condensate atoms – excited atoms.

*“The behaviour of condensates at finite temperatures is a frontier of manybody physics, the experimental exploration of which has so far concentrated mainly on the initial formation of condensates.” (James R. Anglin & Wolfgang Ketterle, Nature, 2002)*

Up to now, under some physical constraints, we have assumed that the interactions between excited atoms are much smaller than the interactions between condensate atoms and excited atoms, which are described by the operator.

In this chapter, we will consider the full system in which both types of interactions are taken into account. Let us denote the operator that describes the interactions between excited atoms by . We then have

- excited atoms – excited atoms: (Uehling-Ulenbeck term),
- condensate atoms – excited atoms: (condensate growth term).

As usual, the dynamics of the condensate itself follows the nonlinear Schrodinger equation (NLS).

Finally, the full system should contain the NLS and a kinetic equation of and . This can be seen in the following Figure:

And the full system for the excited atoms density distribution and the condensate wave function reads (we will come back to the derivation of this equation in my joint work with Linda Reichl in Chapter 6):

where are some constants. is the condensate density, satisfies

and is the Bogoliubov dispersion law

for , are some positive constants, is the mass of the particles, is the interaction coupling constant.

Let us try to understand each term from the equation. From Chapter 1, we know that the condensate growth term describes the interaction between excited atoms and the condensate. Notice that depends on and vanishes when . The Uehling-Ulenbeck (Boltzmann-Norheim) term describes the interaction between excited atoms. Following the same argument for the quantum Boltzmann equation for fermions, we suppose that two microscopic boxes with momenta and collide. After the collision, the particles will move to the new boxes and . For the fermion case, the collision contains the term

where the factor indicates the fullness of the box . From Chapter 1, we also know that bosons do not follow the Pauli Exclusion Principle: the bosons in are very friendly and they are trying to attract their friends to their home. As a result, we use the factor and modify the collision as

Taking into account the conservation of momentum and energy, we obtain

which gives the form of

The operator describes the interaction or the -wave interaction between excited atoms.

Above the BEC critical temperature, the density of the condensate is , then . The homogeneous version of Equation (1) is reduced to the Uehling-Ulenbeck equation

which has a blow-up positive radial solution in the norm if the mass of the initial data is too concentrated around the origin due to a result by Miguel Escobedo and Juan JL Velazquez (Invent. Math. 2015). Let us mention that when the temperature is above the BEC critical temperature, the energy is of the form . The collision of two microscopic boxes of particles with momenta and changes the momenta into and ; and the conservation laws read:

Since , , , belong to the sphere centered at with radius , the collision operator in (2) can be expressed as a integration on a sphere, following the strategy for the classical Boltzmann operator. In (1), is approximated by the Bogoliubov dispersion law (2) , which means that the collision operators are integrals on much more complicated manifolds. New estimates on these energy manifolds are required. Moreover, (3) conserves the mass of the solution, while the full equation (1) does not.

In this joint work with Avy Soffer (Rutgers), we first decouple the QB and NLS equations, then show a global existence and uniqueness result for strong, positive, radial solutions to the spatially homogeneous system.

In the next chapter, we will come back to the condensate growth term to study the convergence to equilibrium.