This post is based on my joint work with Linda E. Reichl.
Let us recall from the previous chapters 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:
(Uehling-Ulenbeck term),
- condensate atoms – excited atoms:
(condensate growth term).
The above model has been remarkably successful in describing several Bose-Einstein Condensate (BEC) dynamical problems, in both the hydrodynamic and collisionless regimes, with extremely broad temperature ranges, and even the description of the condensate growth. It also provides a clear picture of what physical processes are taking place in a partially condensed Bose gas and gives information about the superfluid properties of the system. However, as pointed out by Linda Reichl et. al. (Phys. Rev. A, 2013 & 2014), the description of the thermal component of the BEC is simplified in this system due to the use an approximate particlelike Hartree-Fock excitation spectrum; that makes predictions for transport properties at low temperature inaccurate, for instance, sound speeds at . As a consequence, according to the breakthrough idea of Linda Reichl, a collision operator is missing in the model!

Figure 1: The three collision operators
In order to understand the intuition of the derivation of . We start with the collision operator for classical particles. Suppose that we need to study the density of a spatially homogeneous gas, whose distribution function is
, which indicates the probability of finding a particle with momentum
and time
. Suppose the two particles with momenta
and
collide and they change their momenta into
and
. That means we lose two particles with momenta
and
, and gain two new particles with momenta
and
. The probability of losing two particles with momenta
and
is
and of gaining two particles with momenta
and
is
. As a result, the collision operator has to contain
. Since the collision conserves mass and momentum, we also have
and
, in which
is the energy of the particle with momentum
. As a consequence, the collision operator should be proportional to
. The classical Boltzmann equation then reads
where the collision kernel is omitted for the sake of simplicity.
We will try to apply the above argument for a gas of fermions. For fermions, we cannot consider the collisions of each quantum particle as in the case of classical gas. We then modify the argument as follows: Suppose that two momentum states with momenta and
collide and after the collision some particles will move to the new momentum states
and
. As a result, we will have a term proportional to
in the collision operator. However, fermions follow the Pauli Exclusion Principle, which says that if the momentum states
and
are already full, then particles cannot move from
and
to
and
. We therefore modify the collision as
where the factor indicates that momentum state
is occupied. The quantum Boltzmann equation for fermions can be written as follows
The question now is how to adapt the above arguments for fermions to the case of the operator for bosons. Suppose again that two momentum states with momenta
and
collide but the momentum state
is hidden inside the condensate: we have a collision between
and the condensate. After the collision some particles will move to the new momentum states
and
. Similar as above, in the collision operator, we will have a term proportional to
. If we apply the same argument as above, the collision operator is proportional to
However, different from fermions, bosons do not follow the Pauli Exclusion Principle. Therefore, instead of multiplying with , we use the factor
and modify the collision as
The conservation of momentum and energy in this case are and
. The collision operator should contain
Denoting this case as , we observe that there are still other two cases
and
. Due to the symmetry of
and
, we can combine them into
The form of then follows. The interaction is of the type
. The operator
describes the process that two momentum states with momenta
and
collide and after the collision some bosons will move to the new momentum states
and
. The same argument as above also gives the form of
. The
operator describes the
interaction. Let us look carefully into the structure of
: We call this a four-excitation process due to the fact that there are four particles in each collision. In a four-excitation process, besides the
interaction, there should be another
one. The missing collision operator is, therefore, of the type
, which describes the process of three momentum states
,
and
collide and particles are moving into a new momentum state
.
The new kinetic equation should be (see Figure 3)

Figure 2: The Bose-Einstein Condensate (BEC) and the excited atoms:
The new model has been tested with several experiments. For example, the speed and lifetime of the first and second sound modes in mono-atomic BECs computed by the new model, with the new collision operator , are in agreement with many experimental results. So,
should be the missing piece of the picture.
In this joint work with Linda Reichl at Department of Physics and Center for Complex Quantum System, UT Austin, we revisit a method developed by Peletmiskii and Yatsenko to derive kinetic equations for quantum gases.