text:
# Large amplitude spin waves in ultra-cold gases

## 1 Introduction

Recent experiments  have renewed the interest in spin waves in dilute quantum gases, in conditions where they had not been observed before. A first novelty is that these waves were neither purely nuclear nor electronic, but involved two hyperfine atomic levels where the nuclear and electronic spin are coupled. This difference is actually minor: the exchange effets from which the waves originate are unaffected by the precise nature of the levels - in other words one can, without any loss of generality, assimilate any pair of atomic states to up and down states of a fictitious spin. The second, more important, difference is that the waves were observed with very large amplitudes, leading to an almost complete apparent segregation of the atoms in each internal state, and involving situations where the gas is very far from local equilibrium; we study such situations in the present article.

In degenerate liquids, the existence of spin waves has been known for many years . In dilute gases, which are non-degenerate, it took about twenty more years  to realize that they should also sustain similar waves, for both bosonic and fermionic systems. This came as a surprise to some since, before, these waves had been mostly associated with the Landau formalism for Fermi liquids, and with the notion of “molecular field” resulting from exchange and interactions with nearest neighbors. In a liquid, a test particle remains trapped in a sort of local cavity and interacts constantly with several neighbors, so that the notion of mean field emerges naturally from the averaging over the effect of several neighbors. By contrast, in a gas, the motion of a particle is often described as free flights along straight lines, interrupted by short collisional processes, so that the physics involved is clearly very different. It nevertheless turns out that, when averaged over all possible collisions at each point of space, the final result of the exchange interactions in a gas and a liquid are very similar (see, for example, Ref.  for a general discussion). Soon after these theoretical predictions, experimental observations demonstrated the existence of spin waves in spin polarized Hydrogen gas , 3He gas , and dilute 3He-4He solutions .

The physical difference between gases and liquids has its counterpart in the different theoretical approaches used to derive kinetic equations. Historically, the first derivation was that of Boltzmann, with his Stosszahlansatz (or molecular chaos Ansatz): collisions are considered as point processes (with no duration and no spatial extent) taking place between particles which are completely uncorrelated before collision. As a consequence, one can add at each point of space the effect of all possible binary collisions between uncorrelated particles, as if one was adding the effect of many “beam to beam collision experiments” in atomic physics. This provides the famous Boltzmann equation, with a relatively complicated collision term on the right hand side, while the drift term on the left hand side corresponds merely to a completely free flight of particles between collisions. In the approach introduced later by Landau for the study of degenerate Fermi liquids , quasi-particles are never free, since their motion is constantly guided by a mean field created by the neighbor particles. Therefore interactions modify the drift term while, on the other hand, the collision term is generally treated phenomenologically, by a simple relaxation approximation. Both point of view have their advantages. Boltzmann’s point of view does not require the use of a pseudo-potential and therefore allows a more microscopic treatment of collisions, for instance a more precise inclusion of lateral collisions including their full angular end energy dependence. On the other hand, it does not encompass quantum degenerate systems, does not introduce the powerful notion of quasi-particles, and generally speaking remains limited to dilute interacting gases.

The connection between the two points of view is provided by collisions in the forward direction. In two body collision theory, it is well known that wave interference in the forward direction is responsible for particle absorption (optical theorem). In addition, more physics may be involved in this forward interference effect: for spinless particles, retardation effects are introduced (see for instance ); for particles with spin, in addition, identical spin rotation effects (ISRE) due to particle indistinguishability  also take place. In many-body transport theory, both these effects have their equivalents. The cumulated retardation effects during collisions give rise to an average force, or equivalently to a scalar mean field (see for instance  and references therein); in the kinetic equation, this corresponds to terms which are quadratic in the distribution function, as the collision term, but contain gradients (spatial gradients of the density as well as momentum gradients of the particle distribution); in other words, one gets zero sound Landau type mean field terms which are naturally grouped with the drift term on the l.h.s. of the kinetic equation. As for the identical spin rotation effects, they average to a spin molecular field which is the equivalent of the molecular field considered by Silin  and Leggett ; the ISRE is already present to zero order gradient expansion, as opposed to the scalar retardation terms , which explains why the dispersion relation of spin waves is different from the usual sound wave dispersion. A detailed study of forward scattering in a binary collision can therefore lead to a microscopic understanding of the origin of the Landau mean field (lateral scattering introduces collisional damping). From a practical point of view, it remains true that the mean field treatment of the drift term, necessarily associated with the use of a pseudo-potential, is often more compact and elegant than a detailed study of the collisional properties in the forward direction; it is nevertheless interesting to check precisely to what extent they are equivalent in various situations.

Our purpose in this article is twofold. First, for particles with spins, the equivalence in question has been verified in the literature only to zero order in the gradient expansion, in other words only for the local terms appearing in the kinetic equation. Here we wish to derive the kinetic equation to first (non-local) order in the gradient expansion. Moreover, we will consider the situation where the scattering may depend on the relative spin orientations (even if this dependence was relatively weak in the experiments of the JILA group ). We will therefore examine in detail the properties of collisional interference in the forward direction in this more general case; the details of the calculations are given in two appendices (gradient expansion of the interaction term within the $`S`$-matrix collision Ansatz). Second, we also wish to study situation where the gas is far from local equilibrium. In fact, most of the work on spin waves in the literature deals with hydrodynamic situations where, at every point of space, the gas is close to local equilibrium, so that a simple hydrodynamic expression of the spin current can be used; of course, this reduces the number of variables, but at the price of an approximation which is not necessarily justified. In fact, in the experiments in question, the gas sometimes evolves very far from an hydrodynamic regime, which is not so surprising since the spin waves have a large amplitude and lead to an almost complete segregation of the two spin species. For a study of spin waves in the collisionless regime (mean-free-path large compared to characteristic lengths), but still close to equilibrium (linear regime), see Ref. .

## 2 Kinetic equation

We study an ensemble of identical atoms with two internal levels, obeying either Bose or Fermi statistics. As mentioned in the introduction we can, without any loss of generality, assimilate these two levels to two spin levels, whatever their real physical origin is. For instance, hyperfine atomic structure states, even in a situation of intermediate magnetic decoupling, are possible; the only important thing is that the quantum states describing the internal variables should be orthogonal. We now derive a kinetic equation for a dilute gas of such atoms.

### 2.1 Collisional methods; Yvon-Snider and Lhuillier-Laloë equations

In classical statistical mechanics, Boltzmann’s intuitive method for deriving a kinetic equation is based on the study of individual collisions; one adds the effect of all possible binary collisions taking place at each point of space on the single particle distribution in phase-space $`f_1(𝐫,𝐩)`$. Another point of view starts from the infinite hierarchy of BBGKY equations and closes its first equation, which relates $`f_1`$ to the two particle distribution function $`f_2`$, by studying the evolution of $`f_2`$ during a binary collision and expressing it as some functional of $`f_1`$ (see for instance Ref. ).

The same lines can be followed in quantum mechanics, but $`f_1(𝐫,𝐩)`$ has to be replaced by the reduced single-particle density operator $`\widehat{\rho }_1`$. Since we assume here that the particles have two internal states, which we treat as spin states, the single-particle density operator $`\widehat{\rho }_1`$ acts in the product space of orbital and spin variables of one particle. The operatorial kinetic equation has the form:

$$\frac{d\widehat{\rho }_1}{dt}+\frac{1}{i\mathrm{}}[\widehat{\rho }_1,\widehat{H}_1]_{}\frac{d}{dt}|_{coll}\widehat{\rho }_1$$
(1)
where:

$$\widehat{H}_1=\frac{𝐩^2}{2m}\widehat{1}+\widehat{V}^{ext}$$
(2)
is the single-particle Hamiltonian, $`𝐩`$ the momentum of the particle, $`m`$ its mass, $`\widehat{1}`$ the unity operator in spin space and $`\widehat{V}^{ext}`$ the operator describing the external forces acting on the particle, which may be spin dependent. The r.h.s. of eq. (1) takes into account the effects of binary interactions between the particles; it contains the result of the collisional approximation mentioned above, and will be discussed in more detail below. The equation is only valid for a dilute gas ($`n^{1/3}a`$, where $`n`$ is the number density of the gas and $`a`$ the scattering length) and on time/length scales much greater than the duration of a collision/scattering length.

The next step is to introduce the Wigner transform  $`\widehat{\rho }_W(𝐫,𝐩)`$ of $`\widehat{\rho }_1`$ with respect to orbital variables:

$$\widehat{\rho }_W(𝐫,𝐩)(2\pi \mathrm{})^3d^3𝐫^{}e^{i𝐩.𝐫^{}/\mathrm{}}𝐫\frac{𝐫^{}}{2}|\widehat{\rho }_1|𝐫+\frac{𝐫^{}}{2}$$
(3)
It is a classical function of position and momentum, but still an operator - or a $`2\times 2`$ matrix - in spin space. In the following, we will treat the orbital degrees of freedom semi-classically, but not the spin degrees of freedom. For the sake of simplicity, from now on we will drop the index $`W`$, assuming that the dependence on $`𝐫`$ and $`𝐩`$ is sufficient to signal a Wigner distribution. One then performs a Wigner transform of equation (1), which involves taking the Wigner transform of products of operators. This is possible by using Groenewold’s formula , which provides the result as a infinite series of gradient expansion: the first term is merely the product of the Wigner transforms of the two operators, followed by product of gradients with respect to $`𝐫`$ and $`𝐩`$, followed by higher order gradients, etc. The Wigner transform of the l.h.s. of equation (1) then gives:

$`_t\widehat{\rho }(𝐫,𝐩)`$ $`+{\displaystyle \frac{𝐩}{m}}_𝐫\widehat{\rho }(𝐫,𝐩)+{\displaystyle \frac{1}{i\mathrm{}}}[\widehat{\rho }(𝐫,𝐩),\widehat{V}^{ext}(𝐫)]_{}`$
$`{\displaystyle \frac{1}{2}}[_𝐩\widehat{\rho }(𝐫,𝐩),_𝐫\widehat{V}^{ext}(𝐫)]_++\mathrm{}`$ (4)
The first two terms are exact, whereas the next two are only the first orders terms in a gradient expansion involving the external potential. The commutator is obtained to zero order of the gradient expansion while the anticommutator occurs only to first order, followed by higher order terms symbolized by the dots. The anticommutator corresponds to the classical force term (describing the effect of a trap, for example), which we wish to retain in our calculations; as a consequence, we have to include up to first order gradient expansion in the Wigner transform of the r.h.s. of equation (1) .

The non-trivial part in the derivation of a kinetic equation begins when one gives a precise expression to the formal r.h.s. of equation (1). This is the subject of many studies in the literature and even the subject of books; here, we limit ourselves to a simplified discussion of two different approaches which are appropriate for the study of spin waves. We consider only collisions taking place in elastic channels, where the number of atoms in each internal spin state is conserved, and begin with some simple considerations on elementary collision theory. For each channel, at very low temperatures, collisions occur in the extreme quantum regime where the typical $`s`$-wave scattering length $`a`$ is small compared to the de Broglie thermal wavelength $`\lambda _T`$. In this case, collisions are well described by the isotropic $`s`$-wave, and the $`T`$ matrix can be approximated by:

$$T_k(\widehat{k},\widehat{k}^{})=\frac{\mathrm{}^2}{2\pi ^2m}f_k(\theta )=\frac{\mathrm{}^2}{2\pi ^2m}(aika^2+..)$$
(5)
(terms of order $`a^2`$ should be retained in order to satisfy explicitly the unitarity of the $`S`$-matrix in the optical theorem); in this equation, $`m`$ is the mass of the particles, $`f_k(\theta )`$ the scattering amplitude which, for $`s`$-wave scattering, is independent of the polar angle $`\theta `$, and $`k`$ the collision wave vector.

A first obvious remark is that the intensity of the spherical scattered wave is proportional to the square of $`a`$, while interference effects in the forward direction between the incoming plane wave and the spherical scattered wave are proportional to $`a`$ itself. Generally speaking, the relative magnitude of forward and lateral scattering effects is of order $`\lambda _T/a`$, which is large at low energies. This is why, here, we concentrate on forward scattering effects and their relation to mean field corrections to the drift terms. Another remark is that, because waves interfere only if they are in phase, phases and $`i`$ factors are important here. It turns out that, in quantum mechanics as well as in optics , the summation over many scatterers in the forward direction introduces an additional $`i`$ factor (which is actually the origin of the $`i`$ factor in the optical theorem). Hence, in every scattering channel, equation (5) shows that no first order $`a`$ interference effect occurs in the forward direction, the first contribution arising from the second term in $`ika^2`$; to first order in $`a`$, only the phase of the forward wave is changed, not its intensity. But, if several scattering channels are open, the situation is different, because these different phase shifts in the forward scattered waves can introduce spin rotation. Actually, this is precisely the origin of the first order identical spin rotation effect in the forward direction .

#### 2.1.1 Spin-independent interactions

We now briefly summarize two approaches that have been developed in the context of particles with spin: the Yvon-Snider approximation, which emphasizes more the relation to the BBGKY hierarchy, and closes the infinite set of equations with an approximation suitable for a binary collision; the Lhuillier-Laloë (or $`S`$ matrix) approach, which is closer to the initial approach of Boltzmann and makes use of the $`S`$ collision matrix. We begin with the simplest case, when the interactions are independent of spin.

##### Yvon-Snider

An interesting method is that originally proposed by Yvon , then independently by Snider ; it provides a result which, in the literature, is often called the Waldmann-Snider equation; the main idea is to express the collision term as the trace over a collision partner (particle $`2`$) of a commutator containing the binary interaction potential $`V_{12}`$ and a product of single particle operators modified by the unitary transformation associated with the Møller collision operator $`\mathrm{\Omega }\mathrm{\Omega }^{(+)}`$:

$$\frac{d}{dt}|_{coll}\widehat{\rho }_1(1)=(i\mathrm{})^1Tr_2\left\{[V_{12},\mathrm{\Omega }\widehat{\rho }_1(1)\widehat{\rho }_1(2)\mathrm{\Omega }^{}]_{}\right\}$$
(6)
This equation is valid for dilute gases of distinguishable particles; for non-degenerate gases ($`n\lambda _T^31`$, where $`\lambda _T`$ is the thermal wavelength of the particles) of bosons or fermions, it generalizes into:

$$\frac{d}{dt}|_{coll}\widehat{\rho }_1(1)=(i\mathrm{})^1Tr_2\left\{\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}[V_{12},\mathrm{\Omega }\widehat{\rho }_1(1)\widehat{\rho }_1(2)\mathrm{\Omega }^{}]_{}\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\right\}$$
(7)
where $`\widehat{P}_{ex.}`$ is the exchange operator between particles $`1`$ and $`2`$, and where:

$$ϵ=1\text{for bosons, }ϵ=1\text{for fermions}$$
(8)
Wigner transforms can then be applied to these expressions to complete the derivation of the closed kinetic equation for $`\widehat{\rho }(𝐫,𝐩)`$. For instance, expanding (6) in a Groenewold gradient expansion of the products provides, to zero order in gradients, the straightforward Boltzmann expression for the collision term:

$`{\displaystyle \frac{d}{dt}}|_{coll}\widehat{\rho }(𝐫_1,𝐩_1)`$ $`={\displaystyle }d^3q{\displaystyle \frac{q}{m}}[\sigma _T(k)\widehat{\rho }(𝐫_1,𝐩_1)f(𝐫_1,𝐩_2)`$
$`{\displaystyle }d^2\widehat{k}^{}\sigma _k(\theta )\widehat{\rho }(𝐫_1,𝐩_1^{})f(𝐫_1,𝐩_2^{})]`$ (9)
where $`f(𝐫,𝐩)`$ is the spin trace of $`\widehat{\rho }(𝐫,𝐩)`$:

$$f(𝐫,𝐩)=\text{Tr}_S\left\{\widehat{\rho }(𝐫,𝐩)\right\}$$
(10)
and where $`\widehat{k}`$ is the unit vector in the direction of $`𝐤`$ and $`\theta `$ the angle between $`𝐤`$ and $`𝐤^{^{}}`$. Here we use the usual notation in Boltzmann theory:

$$\begin{array}{cc}𝐪=2\mathrm{}𝐤& 𝐩_2=𝐩_1𝐪\\ 𝐩_1^{}=𝐩(𝐪/2)+\mathrm{}k\widehat{k}^{}& 𝐩_2^{}=𝐩(𝐪/2)\mathrm{}k\widehat{k}^{}\end{array}$$
(11)
as well as the definition of the differential and total cross sections:

$$\begin{array}{cc}\sigma _k(\theta )=(4\pi ^4m^2/\mathrm{}^4)\left|T(𝐤,𝐤^{})\right|^2;& \sigma _T(k)=d^2\widehat{k}^{}\sigma _k(\theta )\end{array}$$
(12)
Higher order terms will introduce scalar mean field terms, which we do not write here explicitly for concision; they can be found in the appendix of ref. .

Now, if we introduce spin and statistics, we have to use (7) instead of (6). To zero order in gradients, this introduces the additional exchange terms to the right hand side of (9):

$$\begin{array}{c}\hfill \frac{ϵ}{2}d^3q\frac{q}{m}\{i\tau _{fwd.}^{ex.}(k)[\widehat{\rho }(𝐫_1,𝐩_1),\widehat{\rho }(𝐫_1,𝐩_2)]_{}\\ \hfill d^2\widehat{k}^{}i\tau _k^{ex.}(\theta )[\widehat{\rho }(𝐫_1,𝐩_1^{}),\widehat{\rho }(𝐫_1,𝐩_2^{})]_{}+\\ \hfill +d^2\widehat{k}^{}\sigma _k^{ex.}(\theta )([\widehat{\rho }(𝐫_1,𝐩_1),\widehat{\rho }(𝐫_1,𝐩_2)]_+[\widehat{\rho }(𝐫_1,𝐩_1^{}),\widehat{\rho }(𝐫_1,𝐩_2^{})]_+)\}\end{array}$$
(13)
where the “generalized cross sections” are defined by:

$$\tau _{fwd.}^{ex.}(k)=(8\pi ^3m/\mathrm{}^2k)\mathrm{Re}\left\{T(𝐤,𝐤)\right\}$$
(14)
and:

$$\sigma _k^{ex.}(\theta )i\tau _k^{ex.}(\theta )=(4\pi ^4m^2/\mathrm{}^4)T(𝐤,𝐤^{^{}})T^{}(𝐤,𝐤^{^{}})$$
(15)
In (13), the two terms in the first two lines correspond respectively to the ISRE in the forward and lateral direction; the two terms in the third line correspond to exchange effects changing the values of the total and lateral cross sections, as discussed in more detail in . For simplicity, we do not write the first order gradient terms, but more details on the limit of validity of the Yvon-Snider equation and of the various terms which it contains can be found in ,  or  (in particular the Appendix of this last reference). Inserting the low energy limit (5) is trivial; we write the explicit results in the next section.

##### Lhuillier-Laloë

Another point of view, more directly in the spirit of the initial Ansatz of Boltzmann, was employed by Lhuillier and Laloë (LL) ; it is based on the use of the $`S`$ collision matrix. When collisions are treated only as “closed processes” (one ignores particles “in the middle of a collision”), it is indeed possible to use the $`S`$ matrix to relate exactly the single particle operators after and before collision. One gets the following expression for the single particle density operator after collision $`\widehat{\rho }_1^{}`$ as a function of the same operator $`\widehat{\rho }_1`$ before collision:

$$\widehat{\rho }_1^{}(1)=Tr_2\left\{\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\widehat{S}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\widehat{S}^{}\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\right\}$$
(16)
(the trace acts over orbital and spin variables of the collision partner $`2`$) where the $`S`$ matrix is related to the collision $`T`$ matrix by:

$$1:𝐤_f;2:𝐤_f|\widehat{S}|1:𝐤_i;2:𝐤_i=\delta (𝐤_f𝐤_i)i\frac{\pi m}{\mathrm{}^2k_i}\delta (k_fk_i)T_{k_i}$$
(17)
($`\mathrm{}k`$ is the relative momentum of the two particles). One then approximates the rate of change of $`\widehat{\rho }_1`$ by the variation $`(\widehat{\rho }_1^{}\widehat{\rho }_1)/\mathrm{\Delta }t`$ during a short time interval 111We will see that $`\mathrm{\Delta }t`$ is a time which emerges naturally from the calculation through the occurrence of square of delta functions of the energy (see Appendix I). $`\mathrm{\Delta }t`$, which is however much longer than the duration of a collision, to obtain:

$$\frac{d}{dt}|_{coll}\widehat{\rho }_1(1)\frac{1}{\mathrm{\Delta }t}Tr_2\left\{\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\left[\widehat{S}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\widehat{S}^{}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\right]\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\right\}$$
(18)
One should note that the $`S`$ matrix Ansatz of equation (18) is only valid for a non-degenerate gas where:

$$Tr_2\left\{\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\right\}=\widehat{\rho }_1(1)+ϵ\widehat{\rho }_1(1)^2\widehat{\rho }_1(1)$$
(19)
As above, the kinetic equation is obtained from equation (18) by performing a Wigner transform followed by a gradient expansion. In Appendix I we first calculate the zero order gradient expansion, which is given by equation (76) ; it is easy to see that it coincides exactly with the sum of (9) and (13), which shows that the Yvon-Snider and $`S`$ matrix approximations are exactly equivalent to this order. The first order gradient terms are also given in the same Appendix; here, we only give their low energy limit by using equation (5), which provides the following kinetic equation:

$`_t\widehat{\rho }(𝐫_1,𝐩_1)+{\displaystyle \frac{𝐩_1}{m}}_𝐫\widehat{\rho }(𝐫_1,𝐩_1)+{\displaystyle \frac{1}{i\mathrm{}}}[\widehat{\rho }(𝐫_1,𝐩_1),\widehat{V}^{ext}(𝐫_1)+ϵg\widehat{n}(𝐫_1)]_{}`$
$`{\displaystyle \frac{1}{2}}[_𝐩\widehat{\rho }(𝐫_1,𝐩_1),_𝐫(\widehat{V}^{ext}(𝐫_1)+gn(𝐫_1)\widehat{1}+ϵg\widehat{n}(𝐫_1))]_+=I_{coll}[\widehat{\rho }]`$ (20)
where $`g4\pi \mathrm{}^2a/m`$ and $`\widehat{n}(𝐫_1)`$ is the local density operator integrated over velocities:

$$\widehat{n}(𝐫)d^3p\widehat{\rho }(𝐫,𝐩)$$
(21)
In this equation, as usual, forward scattering terms linear in the coupling constant $`ga`$ have been included in the l.h.s. The local term of the collision integral, which regroups terms proportional to $`a^2`$, is given by:

$`I_{coll}[\widehat{\rho }]`$ $`={\displaystyle }d^3q{\displaystyle \frac{q}{m}}{\displaystyle }d^2\widehat{k}^{}a^2(\widehat{\rho }(𝐫_1,𝐩_1)f(𝐫_1,𝐩_2)\widehat{\rho }(𝐫_1,𝐩_1^{})f(𝐫_1,𝐩_2^{})`$
$`+{\displaystyle \frac{ϵ}{2}}[\widehat{\rho }(𝐫_1,𝐩_1),\widehat{\rho }(𝐫_1,𝐩_2)]_+{\displaystyle \frac{ϵ}{2}}[\widehat{\rho }(𝐫_1,𝐩_1^{}),\widehat{\rho }(𝐫_1,𝐩_2^{})]_+)`$ (22)
where $`𝐩_2`$, $`𝐩_1^{}`$ and $`𝐩_2^{}`$ are defined in (11). In addition to the anticommutator in the l.h.s. of (20), obtained to first order in the gradient expansion, we recover exactly the low-energy limit of the LL equation  with, in the l.h.s., the ISRE in the forward direction (the ISRE in the lateral directions disappears in the low energy limit); the second line of (22) provides the statistical terms introduced by particle indistinguishability into the total and differential collision cross section.

#### 2.1.2 Spin-dependent interactions

We now study the gas with spin-dependent interactions and, for a moment, consider the particles as distinguishable. In the limit of slow elastic collisions, we need only consider four $`s`$-wave scattering lengths: $`a_{11}`$, $`a_{22}`$, $`a_{12}^{(d)}`$ and $`a_{12}^{(t)}`$. The two first describe collisions processes between particles that are in the same internal state; when they are in othogonal internal states, two different processes occur for distinguishable particles: a direct collision process without energy transfer described by $`a_{12}^{(d)}`$, and a transfer collision process described by $`a_{12}^{(t)}`$ \- similar processes occur in the theory of spin exchange collisions . The $`T`$ matrix is now a $`4\times 4`$ matrix in spin space; to lowest order in the scattering lengths, its matrix elements with distinguishable particles are given by:

$$1:\alpha ;2:\beta |\widehat{T}_k|1:\gamma ;2:\delta =\frac{\mathrm{}^2}{2\pi ^2m}\left(\begin{array}{cccc}\stackrel{~}{a}_{11}(k)& 0& 0& 0\\ 0& \stackrel{~}{a}_{22}(k)& 0& 0\\ 0& 0& \stackrel{~}{a}_{12}^{(d)}(k)& \stackrel{~}{a}_{12}^{(t)}(k)\\ 0& 0& \stackrel{~}{a}_{12}^{(t)}(k)& \stackrel{~}{a}_{12}^{(d)}(k)\end{array}\right)$$
(23)
where $`\alpha `$, $`\beta `$, $`\gamma `$ and $`\delta `$ are the internal state indexes taking values $`\{1,2\}`$ and where:

$$\stackrel{~}{a}_{\alpha \beta }(k)=a_{\alpha \beta }\left[1ika_{\alpha \beta }\right]$$
(24)
Equation (23) replaces (5) when the interactions depend on spin.

##### Kinetic equation

The derivation of the kinetic equation then proceeds along the same lines as in the previous section: we use (16) again to symmetrize the density operator and take quantum statistics into account, but now insert (23) instead of (17) in it, calculate the partial trace, and obtain an interaction term for the single particle density operator; we then Wigner transform the result and use a gradient expansion limited to first order. This calculation is similar to that of Appendix I for spin-independent interactions, but of course more complicated; it provides more general expressions for the mean field terms and the collision integrals. For conciseness, we limit ourselves to the mean field terms, which will actually be sufficient for our calculations below. In this particular case, the terms can actually be obtained more conveniently by an operatorial method, different from that of Appendix I, as shown in Appendix II. We only reproduce the final kinetic equation:

$`_t\widehat{\rho }(𝐫,𝐩)`$ $`+{\displaystyle \frac{𝐩}{m}}_𝐫\widehat{\rho }(𝐫,𝐩)+{\displaystyle \frac{1}{i\mathrm{}}}[\widehat{\rho }(𝐫,𝐩),\widehat{U}(𝐫)]_{}`$
$`{\displaystyle \frac{1}{2}}[_𝐩\widehat{\rho }(𝐫,𝐩),_𝐫\widehat{U}(𝐫)]_+=I_{coll}[\widehat{\rho }]`$ (25)
where the effective single-particle potential $`\widehat{U}`$ is:

$$\widehat{U}=\widehat{V}^{ext}+\left(\begin{array}{cc}(1+ϵ)g_{22}n_2+g_{12}n_1& ϵg_{12}n_{21}\\ ϵg_{12}n_{12}& (1+ϵ)g_{11}n_1+g_{12}n_2\end{array}\right)$$
(26)
with $`g_{\alpha \beta }4\pi \mathrm{}^2a_{\alpha \beta }/m`$, where $`a_{12}a_{12}^{(d)}+ϵa_{12}^{(t)}`$; it turns out that only this combination of the direct and transfer coupling constants enters the calculation (for more details, see Appendix II).

This kinetic equation resembles the Landau-Silin equation  for a normal Fermi liquid in a magnetic field. The second term in the l.h.s. of (25) is the usual free drift term; the anticommutator is a force term including both the effect of an external potential and of the mean field; the commutator is a spin precession term containing effective magnetic fields from different origins (differential Zeeman effect, spin mean field, etc.).

##### Evolution of the density and spin density in phase space

The structure of the distribution function, which is still a matrix in the space of internal variables, is more explicit if we decompose $`\widehat{\rho }`$ and $`\widehat{U}`$ in the basis formed by Pauli matrices $`\widehat{𝝈}`$ and the unit $`2\times 2`$ matrix:

$`\widehat{\rho }(𝐫,𝐩,t)`$ $`={\displaystyle \frac{1}{2}}\left(f(𝐫,𝐩,t)\widehat{1}+𝐌(𝐫,𝐩,t)\widehat{𝝈}\right)`$
$`\widehat{U}(𝐫,t)`$ $`=U_0(𝐫,t)\widehat{1}+𝐔(𝐫,t)\widehat{𝝈}`$ (27)
When written in terms of the phase-space density $`f`$ and spin density $`𝐌`$, equation (25) becomes:

$`_tf+{\displaystyle \frac{𝐩}{m}}_𝐫f_𝐫U_0_𝐩f_r𝐔_𝐩𝐌`$ $`=I_{coll}^{(f)}[f,𝐌]`$
$`_t𝐌+{\displaystyle \frac{𝐩}{m}}_𝐫𝐌_𝐫U_0_𝐩𝐌_𝐫𝐔_𝐩f{\displaystyle \frac{2𝐔}{\mathrm{}}}\times 𝐌`$ $`=I_{coll}^{(M)}[f,𝐌]`$ (28)
In these equations, $`U_0`$ plays the role of the overall trapping potential and $`𝐔`$ is the local effective magnetic field around which the spins rotate. The effective trapping potential $`U_0`$ is given by:

$`U_0(𝐫,t)`$ $`={\displaystyle \frac{V_2^{ext}+V_1^{ext}}{2}}+\left[{\displaystyle \frac{1+ϵ}{2}}(g_{22}+g_{11})+g_{12}\right]{\displaystyle \frac{n}{2}}`$
$`+{\displaystyle \frac{1+ϵ}{2}}(g_{22}g_{11}){\displaystyle \frac{m_{}}{2}}`$ (29)
and the effective magnetic field is:

$$\begin{array}{cc}𝐔(𝐫,t)=& \frac{\mathrm{}\mathrm{\Omega }(x,t)}{2}𝐞_{}+ϵ\frac{g_{12}𝐦(x,t)}{2}\hfill \\ \mathrm{}\mathrm{\Omega }(𝐫,t)=& V_2^{ext}V_1^{ext}+\frac{1+ϵ}{2}\left[(g_{22}g_{11})n+(g_{22}+g_{11}2g_{12})m_{}\right]\hfill \end{array}$$
(30)
where the local density and spin polarization density are defined as:

$`n=n_2+n_1`$ $`;m_{}=n_2n_1`$
$`m_{,1}=n_{21}+n_{12}`$ $`;m_{,2}=i(n_{21}n_{12})`$ (31)
The basis in spin space is denoted by $`\{e_{,1}`$; $`e_{,2}`$; $`e_{}\}`$. The two first vectors define the transverse plane, the last defines the longitudinal direction.

The collision integral in (25), i.e. the equivalent of eq. (22) when the interactions are spin-dependent, is not written explicitly; to simplify the calculations, and because we will mostly study experimental situations where lateral collisions do not play a dominant role, we limit ourselves to a simple relaxation time approximation:

$$I_{coll}^{(f)}\frac{f(𝐫,𝐩,t)f^{eq}(𝐫,𝐩,t)}{\tau }\text{}I_{coll}^{(M)}\frac{𝐌(𝐫,𝐩,t)𝐌^{eq}(𝐫,𝐩,t)}{\tau }$$
(32)
where $`f^{eq}`$ (resp. $`𝐌^{eq}`$) is the local equilibrium (spin) density in phase-space, and $`\tau `$ is a relaxation time. Actually, one could allow for different relaxation times for $`f`$, $`𝐌_{}`$ and $`𝐌_{}`$. Nevertheless, in the case of a non-degenerate gas with spin-independent interactions, it is known (Ref.  and references therein) that these relaxation times are equal; below we consider a non-degenerate gas of 87Rb, for which the three scattering lengths $`a_{11}`$, $`a_{22}`$ and $`a_{12}`$ are very close , so that in practice we will we only need a single relaxation time $`\tau `$ given approximately by:

$$\tau \frac{1}{4\pi a_{12}^2n}\sqrt{\frac{m}{k_BT}}$$
(33)
Finally, the kinetic equations we will use below to discuss spin waves are:

$`_tf+{\displaystyle \frac{𝐩}{m}}_𝐫f_𝐫U_0_𝐩f_𝐫𝐔_𝐩𝐌`$ $`{\displaystyle \frac{ff^{eq}}{\tau }}`$
$`_t𝐌+{\displaystyle \frac{𝐩}{m}}_𝐫𝐌_𝐫U_0_𝐩𝐌_𝐫𝐔_𝐩f{\displaystyle \frac{2𝐔}{\mathrm{}}}\times 𝐌`$ $`{\displaystyle \frac{𝐌𝐌^{eq}}{\tau }}`$ (34)
The main differences with the LL kinetic equation derived in Ref.  are that the full mean field is included (not just the spin mean field), and that the equation is obtained for general spin-dependent interactions. Nevertheless, our result is also less general, since the collisions are assumed to occur only in the $`s`$-wave channel and that the collision integral is treated at the relaxation time approximation level (if necessary, these two restrictions could be lifted without any special difficulty).

### 2.2 Mean field

Another way to derive the modified drift term is to directly use a mean field approximation within field theory. The following second quantized Hamiltonian density is used to describe a trapped ultra-cold gas of bosons or fermions:

$``$ $`=_1+_{int}`$
$`={\displaystyle \underset{\alpha =1,2}{}}\left({\displaystyle \frac{\mathrm{}^2}{2m}}\psi _\alpha ^{}(𝐫)\psi _\alpha (𝐫)+V_\alpha ^{ext}(𝐫)\psi _\alpha ^{}(𝐫)\psi _\alpha (𝐫)\right)`$
$`+{\displaystyle \frac{1}{2}}{\displaystyle \underset{\alpha ,\beta }{}}g_{\alpha \beta }\psi _\alpha ^{}(𝐫)\psi _\beta ^{}(𝐫)\psi _\beta (𝐫)\psi _\alpha (𝐫)`$ (35)
The annihilation and creation field operators of a particle in the internal state $`\alpha `$ are $`\psi _\alpha (𝐫)`$ and $`\psi _\alpha ^{}(𝐫)`$. They have bosonic or fermionic equal-time commutation or anticommutation relations. The interactions only depend on the three coupling constants $`g_{11}`$, $`g_{12}`$ and $`g_{22}`$; the real potential $`V`$-matrix elements have been replaced by the appropriate $`T`$-matrix elements at low-energy, according to an usual procedure in the study of cold gases (Fermi pseudo-potential method , or $`VT`$ renormalization procedure in Ref. ).

We wish to obtain the time evolution of the single-particle density matrix $`\rho _{\alpha \beta }(𝐫,𝐫^{},t)=\psi _\alpha ^{}(𝐫^{},t)\psi _\beta (𝐫,t)`$, where $`\mathrm{}`$ denotes the expectation value with respect to the time-dependent $`N`$-body density operator of the system. To this end we write down the Heisenberg equation of motion for $`\psi _\alpha ^{}(𝐫^{},t)\psi _\beta (𝐫,t)`$ using the Hamiltonian (35). The non-interacting part of the Hamiltonian density $`_1`$ gives the single-particle Liouville-von Neumann equation:

$`i\mathrm{}_t\rho _{\alpha \beta }(𝐫,𝐫^{},t)`$ $`={\displaystyle \frac{\mathrm{}^2}{2m}}\left(_𝐫^2_𝐫^{}^2\right)\rho _{\alpha \beta }(𝐫,𝐫^{},t)`$
$`+\left(V_\alpha ^{ext}(𝐫)V_\beta ^{ext}(𝐫^{})\right)\rho _{\alpha \beta }(𝐫,𝐫^{},t)`$ (36)
The interacting part $`_{int}`$ couples the single-particle density matrix with the expectation value of the product of four field operators (two-particle density matrix). If we are only interested in time scales large compared to the duration of a collision, we can decompose the averages of such products into averages of products of two operators 222This result can be obtained for instance from the Wick theorem. It is valid because the gas is not Bose condensed; otherwise, the mean field could not be obtained in that way, due to the absence of the exchange term.. In the context of the equation of motion for the single-particle density matrix, this approximation is often called the random phase approximation (RPA) . The interacting part of the equation of motion then becomes:

$$\underset{\delta =1,2}{}V_{\alpha \delta }^{mf}(𝐫,t)\rho _{\delta \beta }(𝐫,𝐫^{},t)\underset{\delta =1,2}{}\rho _{\alpha \delta }(𝐫,𝐫^{},t)V_{\delta \beta }^{mf}(𝐫^{},t)$$
(37)
where the mean field potential $`2\times 2`$ matrix is given by:

$$V_{\alpha \beta }^{mf}(𝐫)=\delta _{\alpha ,\beta }\underset{\gamma =1,2}{}g_{\alpha \gamma }\psi _\gamma ^{}(𝐫)\psi _\gamma (𝐫)+ϵg_{\alpha \beta }\psi _\alpha ^{}(𝐫)\psi _\beta (𝐫)$$
(38)
which contains a direct term and an exchange mean field term, proportional to $`ϵ`$. Within this approximation, each particle evolves under the influence of an effective single-particle potential $`\widehat{U}(𝐫,t)`$, which is the sum of the external potential $`\widehat{V}^{ext}(𝐫)`$ and the time-dependent mean field $`\widehat{V}^{mf}(𝐫,t)`$.

In order to compare this result with that of the preceding section, we note that $`n_\alpha (𝐫,t)=\psi _\alpha ^{}(𝐫,t)\psi _\alpha (𝐫,t)`$ is the average density of atoms in internal state $`\alpha =1,2`$ and the off-diagonal elements $`n_{12}(𝐫,t)=n_{21}^{}(𝐫,t)=\psi _1^{}(𝐫,t)\psi _2(𝐫,t)`$ are the coherences between the internal states, so that the effective potential $`\widehat{U}`$ reproduces exactly the one obtained in the preceding section, equation (26). We therefore arrive at the following operatorial kinetic equation:

$$\frac{d\widehat{\rho }_1}{dt}+\frac{1}{i\mathrm{}}[\widehat{\rho }_1,\widehat{H}_1+\widehat{V}^{mf}]_{}=0$$
(39)
which is the quantum equivalent of the Vlasov equation (see Ref.  for example); in this equation, $`\widehat{H}_1`$ is the single particle hamiltonian defined in (2) and $`\widehat{V}^{mf}`$ is defined by (38). Obviously, if we Wigner transform this operatorial equation, we will obtain a kinetic equation that is equivalent to (25), but without the collision integral. The effect of lateral collisions can then be treated phenomenologically by adding by hand a term corresponding to a relaxation time approximation, as usual ; then equations (34) are recovered and the mean field treatment of interactions becomes exactly equivalent to these simplified equations.

### 2.3 Discussion

We have used three different methods to derive a kinetic equation for a gas of atoms with two internal levels: the $`S`$ matrix (or LL) Ansatz, the Yvon-Snider Ansatz (YS) and the mean field approximation; we now summarize a comparison between the results.

To zero order in the gradients, LL and YS give exactly the same equation, containing a Boltzmann-type collision integral and a spin mean field. As for the mean field approximation, it provides only the spin mean field, obtained in the limit of low-energy collisions because of the use of the pseudopotential; this limit is in full agreement with the results of the two other methods. We note in passing that this spin mean field occurs for a non-condensed Bose gas, but not in a Bose-Einstein condensate (BEC) at very low temperature, because the exchange term is absent when a single quantum state is involved. In this case, there is however a different mechanism for creating spin waves, the so-called “quantum torque” of purely kinetic origin ; it arises naturally because, in quantum mechanics, any gradient of the phase corresponds to a kinetic energy for each component which, in turn, affects the evolution of the relative phase and therefore of the transverse spin orientation.

To first order in the gradients and for forward scattering, LL, YS and the mean field approximation all give the same scalar mean field term (zero sound), again in the limit of low-energy collisions for the third method. We have checked that the complete expression of the LL and YS first order gradient terms coincide exactly for forward scattering, but not for lateral scattering; whether or not the two methods are completely equivalent to this order is left as an open question; for more details see .

In a non-degenerate dilute gas, the spin mean field generally dominates over the scalar mean field, because it appears to lower order in the gradient expansion; it is therefore not surprising that zero sound type collective modes should not propagate as easily as spin waves in dilute non-degenerate gases.

To what extent collisional methods and mean field approximations are equivalent is not a trivial question. Generally speaking, mean field theory is assumed to be valid for systems with either long range interactions (Vlasov’s equation for a plasma ) or in strongly interacting degenerate systems described in terms of quasi-particles (Landau’s equation for a Fermi liquid ); for a dilute gas, neither of these conditions is met. Our conclusion is nevertheless that the equivalence is perfect to first order in the scattering length $`a`$ but does not hold to higher orders in $`a`$. We finally remark that, when describing a non-degenerate dilute atomic gas in terms of quasi-particles (see the work of Bashkin , for example), a quasi-particle is just a particle whose kinetic energy $`p^2/2m`$ is shifted by a local mean field term, so that the kinetic drift term of the quasi-particle transport equation is left unchanged.

## 3 Spin oscillations and internal conversion

Our purpose now is to use the kinetic equation obtained in the preceding section to discuss spin waves in ultra-cold trapped atomic gases. We first briefly review a recent experiment realized at JILA  and the theoretical work stimulated by it . We then discuss the non-hydrodynamic character of the observed spin waves, which allows a comparison between the different theoretical treatments. Finally, we describe the decay of the spin waves and discuss the relevance of decoherence.

### 3.1 The JILA experiment

In a recent experiment performed at JILA , Lewandowsky *et al* have studied bosonic 87Rb atoms with two hyperfine states of interest (denoted by $`1`$ and $`2`$), confined in an axially symmetric magnetic trap elongated along the $`Ox`$ direction. Their experiment is described in ; the temperature $`T`$ of the gas in this experiment is about twice the critical temperature for Bose-Einstein condensation, so that the gas is not strongly degenerate and can be treated reasonably well as a Boltzmann gas; however, since the de Broglie thermal wavelength is much larger than the average scattering length, collisions occur in the full quantum regime. Initially, the gas is at equilibrium with only state $`1`$ populated. A $`\pi /2`$ radio frequency pulse is then applied, which suddenly puts all the atoms into the same coherent superposition of states $`1`$ and $`2`$. The subsequent evolution of the system along the axial direction is then monitored by measuring optically the local densities $`n_1`$ and $`n_2`$ of atoms in each internal state. Experimentally, one observes that the system “segregates”: after about $`100`$ ms, atoms in state $`1`$ are mostly found away from the center of the trap, while atoms in state $`2`$ move towards the center, the total density $`n=n_1+n_2`$ remaining practically unchanged. After about $`200`$ ms, the local densities of each species return to equilibrium.

In Ref. , this phenomenon was explained in terms of an internal conversion resulting from the identical spin rotation effect (ISRE), which is at the origin of the spin mean field; similar considerations were almost simultaneously provided by two other groups . Qualitatively, this spin oscillation can be understood as follows:

(i) the field gradient creates an inhomogeneous spin precession, so that the gas develops a gradient of transverse spin polarization

(ii) the thermal motion of the atoms then creates correlations between transverse spin and velocity; therefore, a particle moving with a given velocity at point $`x`$ gets a spin polarization which is not parallel to the average local polarization

(iii) the ISRE then makes its spin polarization leave the transverse plane and get a longitudinal component, in a direction which depends on the sign of the velocity

(iv) this is equivalent to a velocity dependent internal conversion which, at some later time, results in an efficient spatial separation of atoms in the two internal states.

To study this phenomenon quantitatively, we derive from (34) an effective one-dimensional kinetic equation in terms of the local density in phase space $`f(x,p,t)`$ and spin density $`𝐌(x,p,t)`$, by assuming that radial equilibrium is quickly established (the radial trap frequency $`\omega _{rad}/2\pi `$ is much larger than the axial trap frequency $`\omega _{ax}/2\pi `$). Integrating (34) over radial coordinates $`y,z`$ and momenta $`p_y,p_z`$ provides the following equations:

$`_tf+{\displaystyle \frac{p}{m}}_xf_xU_0_pf_x𝐔_p𝐌`$ $`(ff^{eq})/\tau `$ (40)
$`_t𝐌+{\displaystyle \frac{p}{m}}_x𝐌_xU_0_p𝐌_x𝐔_pf`$ $`{\displaystyle \frac{2𝐔}{\mathrm{}}}\times 𝐌`$
$`(𝐌𝐌^{eq})/\tau `$ (41)
(every quantity depending initially on three-dimensional coordinates $`𝐫`$ and momenta $`𝐩`$ has been integrated over radial coordinates and momenta, so that it now depends only on the coordinate $`x`$ and the momentum $`pp_x`$). In the process of radial averaging, the coupling constants $`g_{\alpha \beta }`$ are renormalized by a factor $`1/2`$, as discussed by Levitov . The three-dimensional effective trapping potential (29) therefore becomes:

$`U_0(x,t)`$ $`={\displaystyle \frac{V_2^{ext}+V_1^{ext}}{2}}+\left[{\displaystyle \frac{1+ϵ}{2}}(g_{22}+g_{11})+g_{12}\right]{\displaystyle \frac{n}{4}}`$
$`+{\displaystyle \frac{1+ϵ}{2}}(g_{22}g_{11}){\displaystyle \frac{m_{}}{4}}`$ (42)
and the three-dimensional effective magnetic field (30) is changed into:

$$𝐔(x,t)=\frac{\mathrm{}\mathrm{\Omega }(x,t)}{2}𝐞_{}+ϵ\frac{g_{12}𝐦(x,t)}{4}$$
(43)
with:

$$\mathrm{}\mathrm{\Omega }(x,t)=V_2^{ext}V_1^{ext}+\frac{1+ϵ}{2}\left[(g_{22}g_{11})\frac{n}{2}+(g_{22}+g_{11}2g_{12})\frac{m_{}}{2}\right]$$

where $`𝐞_{}`$ is the unit vector in the longitudinal direction in spin space.

A few simplifying assumptions are appropriate for the experimental conditions of Ref. . The confining energy $`(V_2^{ext}+V_1^{ext})/2`$ is of order $`k_BT13`$ kHz$`\times h`$ which is much larger than the mean field interaction energy $`gn(0)140`$ Hz$`\times h`$. This allows us to keep only the confining energy of the harmonic trap in the effective trapping potential $`U_0`$:

$$U_0(x)\frac{V_1^{ext}(x)+V_2^{ext}(x)}{2}=\frac{1}{2}m\omega _{ax}^2x^2$$
(44)
The differential trapping energy $`V_2^{ext}V_1^{ext}10`$ Hz$`\times h`$ is even weaker than the mean field interaction energy. This ensures that the force terms $`_x𝐔_p(f\text{ or }𝐌)`$ appearing in equations (40,41) are negligible. On the contrary, when the local effective magnetic field $`𝐔`$ does not appear under a spatial gradient, as in the term $`2𝐔\times 𝐌/\mathrm{}`$, it can not be neglected. This effective magnetic field $`𝐔`$ is made of two terms (see equation (43)): one is an effective *external* magnetic field $`\mathrm{}𝛀/2`$; the other is an *exchange* magnetic field or spin mean field $`g_{12}𝐦/4`$ which results from the ISRE. The effective external magnetic field is the sum of the contributions of a differential Zeeman and of a differential mean field:

$$\mathrm{}𝛀(x)=\mathrm{}\mathrm{\Omega }(x)𝐞_{}[V_2^{ext}(x)V_1^{ext}(x)+(g_{22}g_{11})n(x)/2]𝐞_{}$$
(45)
where, following ref. , we have assumed that $`2g_{12}g_{11}+g_{22}`$ for simplicity. The average value over the sample of the effective external magnetic field can be removed by going to a uniformly rotating frame (Larmor frame). For the numerical simulations, we need to know the $`x`$ dependence of this effective magnetic field. As the profile $`n(x)`$ will be shown to be a Gaussian (as for a non-interacting non-degenerate gas) which does not vary in time, we model the experimentally measured effective external magnetic field by:

$$𝛀(x)=\delta \mathrm{\Omega }\mathrm{exp}(m\omega _{ax}^2x^2/2k_BT)𝐞_{}$$
(46)
where the parameter $`\delta \mathrm{\Omega }`$ is the variation of $`\mathrm{\Omega }(x)`$ between the center and the edge of the cloud.

Taking into account all approximations mentioned above, we finally arrive at equations:

$`_tf+{\displaystyle \frac{p}{m}}_xfm\omega _{ax}^2x_pf`$ $`(ff^{eq})/\tau `$ (47)
$`_t𝐌+{\displaystyle \frac{p}{m}}_x𝐌m\omega _{ax}^2x_p𝐌`$ $`(𝛀+{\displaystyle \frac{g_{12}𝐦}{2\mathrm{}}})\times 𝐌`$
$`(𝐌𝐌^{eq})/\tau `$ (48)
The initial equilibrium Maxwell-Boltzmann distribution $`f(x,p)`$ solves the kinetic equation (47), so that the dynamics after the $`\pi /2`$ pulse can be expressed in terms of $`𝐌`$ only. In addition, the density $`n(x)`$ is obtained by integrating the equilibrium Maxwell-Boltzmann distribution over momentum $`p`$, which shows that it is time independent; from now on, we will only consider equation (48). The spin distribution immediately after the $`\pi /2`$ pulse is assumed to be:

$$𝐌(x,p,t=0)=f(x,p)𝐞_{,1}$$
(49)
Equation (48) can be solved numerically with values of the parameters taken from Ref. : the trap frequencies are $`\omega _{ax}/2\pi =7`$ Hz and $`\omega _{rad}/2\pi =230`$ Hz; the average time between collisions is $`\tau 10`$ ms; the temperature is $`T0.6`$ $`\mu `$K; the density at the center of the trap is $`n(0)=1.8\times 10^{13}`$ cm-3; $`\delta \mathrm{\Omega }/2\pi `$ is typically $`12`$ Hz, and the scattering lengths are $`a_{11}=100.9a_0`$, $`a_{22}=95.6a_0`$ and $`a_{12}=98.2a_0`$ where $`a_0`$ is the Bohr radius. The solution at the center of trap is plotted in Figure 1; it shows very good agreement with the experimental observations, without any adjustable parameter.

### 3.2 Large, non-hydrodynamic, spin waves

Ordinary spin waves in non-degenerate gases usually occur either in the hydrodynamic regime or/and in a regime of small amplitudes . One can then assume a small departure of the spin distribution (from either local or global equilibrium); this is similar to the first order Chapman-Enskog gradient expansion, which leads to the Navier-Stokes equations (see for instance ). The result is the well known Leggett equations , transposed from degenerate Fermi liquids to dilute non-degenerate gases :

$$\begin{array}{cc}_t𝐦+_x𝐣=𝛀\times 𝐦& \\ _t𝐣(𝛀+\frac{g_{12}𝐦}{2\mathrm{}})\times 𝐣+\frac{k_BT}{m}_x𝐦+\omega _{ax}^2x𝐦& \frac{𝐣}{\tau }\end{array}$$
(50)
where the spin current along $`Ox`$ is defined by:

$$𝐣(x,t)𝑑p\frac{p}{m}𝐌(x,p,t)$$
(51)
The momentum dependence has disappeared from these equations, which contain only position variables. These equations can describe spin waves, not only in the hydrodynamic regime, but also in the collisionless regime as, for example, in liquid 3He (Silin spin waves ) or in H$``$ gas . For a discussion of the hydrodynamic-like description of spin waves in the collisionless regime of non-degenerate dilute gases, see Refs. .

The spin waves observed at JILA were of large amplitude and occurred in the intermediate regime between hydrodynamic and collisionless, where $`\omega _{ax}\tau 1`$, so that there is no a priori reason to believe in the validity of the Leggett equations; we therefore need to resort to a numerical solution of the kinetic equation in terms of both position and momentum variables. Figure 2 shows the results of this calculation, and the spin density distribution (longitudinal and transverse) in $`p`$-space for different times. For comparison, the local equilibrium spin distribution is also plotted (dashed line); we see that the spin distribution in $`p`$-space can indeed get very distorted, which illustrates a strongly non-hydrodynamic situation.

For very short times, $`t40`$ ms, the transverse and longitudinal spin distributions are still close to local equilibrium ($`M_{||}^{eq}=0`$), so that the Leggett equations are valid. For $`t40`$ ms, this is no longer true, so that one would lose a significant part of the physics by not retaining the full momentum dependence of the distribution. In other words, studying the two first moments of the distribution ($`𝐦`$ and $`𝐣`$) is not equivalent to studying the full distribution $`𝐌(p)`$. What actually takes place is a phenomenon analogous to a density “shock wave” rather than a “sound wave”. Eventually, for $`t160`$ ms (end of the “segregation”), since equilibrium is almost reached, the Leggett equations become valid again.

The initial part of the internal conversion (or state separation) can be treated analytically in several approximations, allowing for comparison between different approaches such as those of Refs. . In Ref. , we solved the kinetic equation analytically, for times $`t\tau `$, by expanding in Taylor series in time and obtained:

$$\frac{m_{||}(x,t)}{n(x)}=\frac{g_{12}n(x)}{2\mathrm{}}[\frac{k_BT}{m}\mathrm{\Omega }^{\prime \prime }(x)2\omega _{ax}^2x\mathrm{\Omega }^{}(x)]\frac{t^4}{4!}$$
(52)
where $`\mathrm{\Omega }^{}`$ and $`\mathrm{\Omega }^{\prime \prime }`$ are the first and second spatial derivative of $`\mathrm{\Omega }(x)`$.

This formula predicts an initial quartic in $`t`$ behavior, which correctly reproduces the numerical results; each power of $`t`$ corresponds to one of the four physical processes (i) to (iv) described in § 3.1.

Another approach to calculate the short time behavior is that used by Williams *et al* . They use the Leggett equations which are valid for $`t\tau `$ because the spin distribution is at equilibrium at $`t=0`$. Solving the Leggett equations in the small time limit, where $`𝐣/\tau `$ is negligible, they recover (52).

A third approach is to consider, following Oktel and Levitov ), the Leggett equations in the hydrodynamic regime. The approximation consists in assuming that the spin current remains close to its stationary value, which allows to neglect the term $`_t𝐣`$ in equation (50). Using the resulting equations, we obtain a different behavior for small times:

$$\frac{m_{||}(x,t)}{n(x)}\frac{g_{12}n(x)}{2\mathrm{}}[\frac{k_BT}{m}\mathrm{\Omega }^{\prime \prime }(x)2\omega _{ax}^2x\mathrm{\Omega }^{}(x)]\frac{(\tau t)^2}{4!}$$
(53)
(the effect of the hydrodynamic approximation is to replace $`t^2`$ by $`t\tau `$). The $`t^4`$ and $`t^2`$ results are plotted in Fig. 3; we see that times smaller than a time between collisions $`\tau `$ can not be dealt with in this approach. Nevertheless, after each atom has made one collision on average, the spin current is close to its stationary value and the hydrodynamic approximation is reasonably valid. Therefore, when $`t>\tau 10`$ ms, we get:

$$\frac{m_{||}(x,t)}{n(x)}\frac{g_{12}n(x)}{2\mathrm{}}[\frac{k_BT}{m}\mathrm{\Omega }^{\prime \prime }(x)2\omega _{ax}^2x\mathrm{\Omega }^{}(x)]\frac{\tau ^2(tt_0)^2}{4!}$$
(54)
(the parameter $`t_0`$ includes the accumulated effect of the retardation of the spin current $`𝐣`$ for small time). The $`(tt_0)^2`$ behavior is plotted in Fig. 4 (the reason why, after $`40`$ ms, the numerical solution of the kinetic equation is not well approximated by the $`(tt_0)^2`$ law is that the spin distribution gets very distorted, so that the Leggett equations are no longer valid).

### 3.3 Decay of the phenomenon

Up to this point, we have studied mostly the initial part of the phenomenon, when the amplitude of the spin wave grows and the two species tend to separate from each other; we now discuss its final part, when the spin oscillations decay and, eventually, the system returns to equilibrium. A first remark is that it is relatively easy to obtain an order of magnitude of the time at which the transition between the two regimes occurs, i.e. the time at which the maximum of the separation takes places. For short times, we have seen that this separation is an indirect result of a differential precession of the transverse component of the spin at different points along the $`Ox`$ axis; the thermal motion of the atoms then creates correlations between the transverse orientation of the spins and the sign of the $`x`$ component of the velocities. But, for long times, the differential rotation of the spins will be so large that each sign of this velocity will become associated with a widely open fan of transverse orientation, with almost zero average; clearly, the selective internal conversion effect will then also average to almost zero. Let us consider the simple situation of a homogeneous gas in a box of size $`L`$ with a linear gradient of external effective magnetic field $`\mathrm{\Omega }^{}`$. We focus, for instance, on the situation at the center of the box. At time $`t`$, atoms initially at $`x=v\times t`$ will cross this point with a transverse orientation rotated by an angle of the order of:

$$\mathrm{\Omega }^{}\times x\times t=\mathrm{\Omega }^{}vt^2$$
(55)
Therefore, at times greater than:

$$t_m\frac{1}{\sqrt{\mathrm{\Omega }^{}v}}$$
(56)
the transverse directions of the spins average out so that the apparent segregation effect no longer takes place. In the case of the trapped gas, $`\mathrm{\Omega }^{}\delta \mathrm{\Omega }/L`$ and $`v\omega _{ax}L`$ where $`L`$ is the size of the cloud, and we obtain:

$$t_m\frac{2\pi }{\sqrt{\delta \mathrm{\Omega }\omega _{ax}}}$$
(57)
Typical values give $`t_m100`$ ms in accordance with the results of the simulation. After this maximum, the system tends to return to equilibrium during a time period which we call the decay of the phenomenon.

Fig.1 shows that the decay of the longitudinal and transverse spin polarizations are significantly different. The longitudinal spin polarization $`m_{}(x,t)`$ returns to equilibrium in two steps: it first quickly decreases from its maximum value to almost zero in a time scale of the order of $`50`$ ms; it then goes to zero much more slowly, on a time scale of the order of $`200`$ ms. In order to understand this, we note in Fig.1 that, at the maximum of separation, the wave packets (or the clouds) associated with each internal state overlap only little; the first step is then easily interpreted as the free motion of the clouds under the restoring force of the trap. It takes a quarter of a period $`T_{ax}/4=\pi /2\omega _{ax}=36`$ ms to go from the maximum of oscillation to the center of the trap. Actually, as soon as the two clouds start to overlap, they interact through a repulsive mean field ( $`g>0`$), so that it takes a little more time than a quarter of a period to reach the situation where the two clouds overlap significantly 333That the two clouds interact mainly trough the mean field and not trough (lateral) collisions is confirmed by the numerical solution of the kinetic equation. Indeed, we checked that removing the collision integral does not change the decay time scale. The role of the collision integral is merely to damp revivals of the longitudinal spin waves, which are observed when the collision integral is discarded.. At this point, the subsequent evolution of the system can be understood as the mutual diffusion of two gases 444Actually, an instability of the transverse spin polarization is possible when there is a strong longitudinal spin polarization gradient (which is the case when the two clouds overlap again) and a spin mean field. This is known as Castaing’s instability . In Ref. , we checked that it does not play a role in the experiment done at JILA .; this happens $`150`$ ms after the initial pulse and can again be described by the Leggett equations (we checked on the $`p`$ distribution that the system is close to local equilibrium). The relaxation time scale is then given by the spin diffusion time $`t_{diff}L^2/D`$ (where $`L2\sqrt{k_BT/m\omega _{ax}^2}`$ is the size of the cloud and $`D=k_BT\tau /m`$ is the spin diffusion constant). Using the same parameters as in the simulation, we obtain $`t_{diff}200`$ ms, in good agreement with the numerical result.

We now discuss the transverse spin polarization $`𝐦_{}(x,t)`$. As seen on Fig.1, it oscillates many times before going to zero. After $`150`$ ms, when the Leggett equations are valid again, we know  that spin waves occur in the hydrodynamic regime. The spectrum derived from these equations predicts damped transverse spin waves . The damping time is of the order of the spin diffusion time $`t_{diff}200`$ ms, which gives an order of magnitude for the time it takes the transverse spin polarization to go from its value at the maximum of separation to zero. This is in agreement with the relaxation time obtained by numerically solving the kinetic equation (Fig.1). We finally note that the frequency of the transverse spin waves is essentially given by the external effective magnetic field precession frequency $`\delta \mathrm{\Omega }/2\pi =12`$ Hz.

### 3.4 “Ghost” wave packet and decoherence

We now come back to the time at which the separation is maximum. In order to simplify the discussion, we consider two clouds of atoms, one corresponding to state $`1`$ and located around $`x=d/2`$, the other corresponding to state $`2`$ located around $`x=d/2`$, where $`d`$ is the distance between the center of the two clouds; we note $`\delta x`$ the width of each of these wave packets, assuming that $`\delta x<d`$. In this case, elementary quantum mechanics predicts that the local transverse spin polarization vanishes everywhere (this is because the local spin density corresponds to an operator which is local in ordinary space). Nevertheless, the superposition of the two wave packets remains coherent so that, rigorously speaking, the system is not equivalent to a classical mixture of two gases, in opposition to what we have assumed in the previous section.

One may wonder how this coherence translates in terms of the Wigner distribution. The answer is well known: the presence of the coherence is contained in a so-called “ghost” wave packet, which exists around the middle point between the two real wave packets (or clouds). For spinless particles, it “carries interferences with it” , and reconstructs them as soon as the real packets overlap again. For particles with spin, when the two separated wave packets have opposite spin orientations, the ghost wave packet centered at $`x=0`$ appears only in the non-diagonal spin matrix elements (transverse spin components of the Wigner distribution). In both cases, the wave packet rapidly oscillates as a function of momentum $`p`$, taking positive and negative values, with an oscillation period $`2\pi \mathrm{}/d`$; the further apart the two clouds, the faster the oscillation. On the other hand, for a statistical mixture of the internal states, the ghost wave packet does not exist. Such oscillating wave packets are well-known in the context of macroscopic superpositions of states .

The next question is to what extent the ghost wave packet is preserved by the time evolution of the kinetic equation. Since the Wigner formalism in itself implies no approximation (it is strictly equivalent to the use of operators), the question arises only because our kinetic equation does not provide the exact quantum evolution: it was actually obtained from the operatorial (quantum) equation for the single-particle density operator by truncating the gradient expansion of the Wigner transform. In other words, we made a semi-classical approximation. But, for the ghost wave packet, the $`p`$ oscillation period is $`2\pi \mathrm{}/d`$ while the width of the wave packet is $`\delta x`$, so that the parameter of the gradient expansion is $`d/\delta x`$, which may be larger than $`1`$: therefore, our kinetic equation does not necessarily remain correct for this case; the semi-classical gradient expansion may smooth the oscillations of the ghost wave packet, and artificially introduce decoherence. Another related source of artificial decoherence555Concerning the second source of artificial decoherence, we think that it is not a problem in our simulations. We checked that, when removing the collision integral so that the kinetic equation becomes reversible, the transverse spin polarization disappears when the two wave packets do not overlap but reappears when they mix again. might come from the discretization of phase-space (necessary for the numerical solution of the kinetic equation), which may introduce a lattice spacing in momentum space larger than $`\delta p`$.

Up to this point, no real physical decoherence mechanism was included in our discussion: we have assumed that the system is perfectly isolated. But this is not the case in practice. In order to obtain an estimate of the physical decoherence time, we proceed by analogy with a well-understood situation in quantum optics: a macroscopic superposition of two coherent states of the electro-magnetic field in a cavity. In that case, it is known that the coherence is lost as soon as one photon on average escapes the cavity and goes in the environment . We therefore assume here that for our trapped gas, a source of decoherence comes from the fact that atoms can leave the trap, because of 3-body collisions for example. The time it takes one atom to leave the trap is $`1/\alpha _{3b}n(0)^2N`$ where $`\alpha _{3b}4\times 10^{29}`$ cm6/s is the 3-body recombination rate constant , $`N10^6`$ is the total number of atoms and $`n(0)2\times 10^{13}`$ cm-3 is the density at the center of the trap. This gives $`0.1`$ ms as a rough upper bound to the physical decoherence time, which is still much shorter than the characteristic time scale of the spin state segregation. Our conclusion is therefore that, when the wave packets recombine, coherence is not likely to play a role in these experiments; this justifies our discussion in the preceding section.

## 4 Conclusion

We have derived a kinetic equation for a Boltzmann gas with two internal levels using different methods. A first conclusion which emerges from this work is that, for dilute gases at low temperature, several significantly different theoretical approaches lead essentially to the same result. Here we have discussed the Yvon-Snider method, with its peculiar way to close the BBGKY hierarchy, the $`S`$ matrix method which is very close to the spirit of the original Boltzmann equation, as well as the popular mean field method. In the latter case, the equivalence actually holds only in the low temperature limit, when forward scattering dominates over lateral scattering by a factor of the order of $`\lambda _T/a`$. Whether or not it is possible to find experimental situations where the results of these theories are significantly different, and can be tested experimentally, is still an open question (for instance, whether or not ISRE in lateral collisions may become dominant).

A second conclusion is that a correct treatment of the correlations between internal variables and velocities at each point of space may be important; this is illustrated by the figures contained in figure 2. In other words, hydrodynamic equations are not always appropriate. What is remarkable is that the effects of local field inhomogeneity and velocities combine non-linearly to create this almost complete separation of the two internal states; the effect is not a segregation of states but rather an internal conversion that depends on the direction of the velocity. It is probably even more remarkable that this effect was discovered experimentally, without any theoretical prediction as a guide to the appropriate experimental conditions, and appeared basically as a $`100\%`$ effect from the beginning. Clearly, ISRE is not a small quantum correction to the effect of collisions, but may dominate the entire dynamics of the quantum gas!

It would be interesting to extend the present work in a few directions, for instance exploring the influence of a temperature gradient on the spin current in hydrodynamical situations, as predicted in . For bosons, it would also be important to understand better how spin waves evolve progressively from the non-condensed regime, which we have studied in this article, to the condensed low temperature regime; the theory recently developed in  seems to provide an appropriate tool for this purpose.

Acknowledgments: Le LKB est UMR 8552 du CNRS, de l’ENS et de l’Université P. et M. Curie, Paris. Part of this work was made during a summer visit at the ECT\* center of theoretical physics near Trento (Italy), to which the authors are very grateful. They also express their gratitude to J. Williams for many stimulating discussions at ECT\* and, later, by e-mail.

APPENDICES: INTERACTION TERM OF THE S-MATRIX ANSATZ

## Appendix I: spin-independent interactions

In this appendix, we give the detailed calculation of the Wigner transform of the collision term of the $`S`$ matrix Ansatz, equation (18). We will closely follow the notations of a similar calculation done in the Appendix of Ref. , where the Wigner transform of the Yvon-Snider Ansatz, equation (6) without internal levels and without exchange, is obtained to first order in the gradient expansion. We define the on-shell $`T`$ matrix $`𝒯`$ by:

$$\widehat{S}=\widehat{1}i2\pi \widehat{𝒯}$$
(58)
and rewrite the core of the Ansatz in equation (18) as:

$`\widehat{S}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\widehat{S}^{}\widehat{\rho }_1(1)\widehat{\rho }_1(2)`$ $`=i2\pi \widehat{𝒯}\widehat{\rho }_1(1)\widehat{\rho }_1(2)+\text{c.c.}`$
$`+(2\pi )^2\widehat{𝒯}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\widehat{𝒯}^{}`$ (59)
where c.c. is the complex conjugate of the preceding term. The first two terms of the preceding equation are linear in the $`T`$ matrix, the last is quadratic. The Wigner transform $`\widehat{F}_T`$ of the linear in $`T`$ matrix operator between curly brackets in (18) is:

$$\widehat{F}_T(𝐑,𝐫,𝐏,𝐩)=i(2\pi )^5\mathrm{}^6d^3Kd^3\kappa e^{i𝐊𝐑}e^{i\kappa 𝐫}$$

$$\times 𝐊_+,𝐤_+|\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\widehat{𝒯}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}|𝐊_{},𝐤_{}+\text{c.c.}$$
(60)
where:

$$𝐊_\pm =\frac{𝐏}{\mathrm{}}\pm \frac{𝐊}{2};𝐤_\pm =\frac{𝐩}{\mathrm{}}\pm \frac{𝐤}{2}$$
(61)
With two closure relations for relative wavevectors $`𝐤_1`$ and $`𝐤_2`$, the r.h.s. of (60) becomes equal to:

$$\frac{1}{2i(2\pi )^5\mathrm{}^6}d^3Kd^3\kappa d^3k_1d^3k_2e^{i𝐊𝐑}e^{i𝜿𝐫}$$

$$\times \left(T(𝐤_+,𝐤_1)+ϵT(𝐤_+,𝐤_1)\widehat{P}_S\right)\delta (E_{k_+}E_{k_1})$$

$$\times 𝐊_+,𝐤_1|\widehat{\rho }_1(1)\widehat{\rho }_1(2)|𝐊_{},𝐤_2\left(\delta (𝐤_2𝐤_{})+ϵ\delta (𝐤_2+𝐤_{})\widehat{P}_S\right)+\text{c.c.}$$
(62)
where $`\widehat{P}_S`$ is the exchange operator in spin space. Now using the inverse Wigner transform formula, we can introduce into (62) the Wigner transform $`\widehat{\rho }(𝐫,𝐩)`$ of $`\widehat{\rho }_1`$ and replace the matrix element of the product of $`\widehat{\rho }`$’s by:

$$\begin{array}{c}\mathrm{}^6d^3R^{}d^3r^{}e^{i𝐊𝐑^{}}e^{i(𝐤_2𝐤_1)𝐫^{}}\hfill \\ \times \widehat{\rho }(𝐑^{}+\frac{𝐫^{}}{2},\frac{𝐏+\mathrm{}𝐤_1+\mathrm{}𝐤_2}{2})\widehat{\rho }(𝐑^{}\frac{𝐫^{}}{2},\frac{𝐏\mathrm{}𝐤_1\mathrm{}𝐤_2}{2})\hfill \end{array}$$
(63)
When inserting the result into (62), a delta function appears:

$$d^3Ke^{i𝐊(𝐑𝐑^{})}=(2\pi )^3\delta (𝐑𝐑^{})$$
(64)
Tracing over the orbital and the spin space of particle $`2`$, we finally obtain the following expression for the Wigner transform $`\widehat{I}_T(𝐫_1,𝐩_1)`$ for the terms linear in $`T`$ in the Ansatz (18):

$`\widehat{I}_T(𝐫_1,𝐩_1)`$ $`={\displaystyle \frac{1}{2i(2\pi )^2\mathrm{\Delta }t}}{\displaystyle d^3qd^3rd^3\kappa d^3k_1^{}d^3k_2^{}d^3r^{}}`$
$`\times e^{i𝜿𝐫}e^{i(𝐤_2^{}𝐤_1^{})𝐫^{}}\delta (E_{k_+}E_{k_1^{}})`$
$`\times [\delta (𝐤_2^{}𝐤_{})T(𝐤_+,𝐤_1^{})\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })f(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })`$
$`+ϵ\delta (𝐤_2^{}𝐤_{})T(𝐤_+,𝐤_1^{})\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })`$
$`+ϵ\delta (𝐤_2^{}+𝐤_{})T(𝐤_+,𝐤_1^{})\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })`$
$`+\delta (𝐤_2^{}+𝐤_{})T(𝐤_+,𝐤_1^{})f(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })]`$
$`+\text{c.c.}`$ (65)
with notation (11) and with:

$$\begin{array}{cc}𝐫=𝐫_1𝐫_2\text{ ;}& 𝐤_\pm =\frac{𝐩}{\mathrm{}}\pm \frac{𝜿}{2}\\ 𝐩_1^{\prime \prime }=𝐩_1\frac{𝐪}{2}+\mathrm{}\frac{𝐤_1^{}+𝐤_2^{}}{2}\text{ ;}& 𝐩_2^{\prime \prime }=𝐩_1\frac{𝐪}{2}\mathrm{}\frac{𝐤_1^{}+𝐤_2^{}}{2}\end{array}$$
(66)
The distribution $`f`$ is the spin trace of $`\widehat{\rho }`$, see equation (10).

The same kind of calculation can be done for the term quadratic in the $`T`$ matrix in (18). It gives the Wigner transform $`\widehat{I}_{T^2}(𝐫_1,𝐩_1)`$ as:

$`\widehat{I}_{T^2}(𝐫_1,𝐩_1)`$ $`={\displaystyle \frac{1}{4\pi \mathrm{\Delta }t}}{\displaystyle d^3qd^3rd^3\kappa d^3k_1^{}d^3k_2^{}d^3r^{}}`$
$`\times e^{i𝜿𝐫}e^{i(𝐤_2^{}𝐤_1^{})𝐫^{}}\delta (E_{k_+}E_{k_1^{}})\delta (E_k_{}E_{k_2^{}})`$
$`\times [T(𝐤_+,𝐤_1^{})T(𝐤_{},𝐤_2^{})^{}\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })f(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })`$
$`+ϵT(𝐤_+,𝐤_1^{})T(𝐤_{},𝐤_2^{})^{}\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_2^{\prime \prime })\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_1^{\prime \prime })`$
$`+ϵT(𝐤_+,𝐤_1^{})T(𝐤_{},𝐤_2^{})^{}\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })`$
$`+T(𝐤_+,𝐤_1^{})T(𝐤_{},𝐤_2^{})^{}f(𝐫_1{\displaystyle \frac{𝐫𝐫^{}}{2}},𝐩_1^{\prime \prime })\widehat{\rho }(𝐫_1{\displaystyle \frac{𝐫+𝐫^{}}{2}},𝐩_2^{\prime \prime })]`$ (67)
As in Ref. , we now assume that $`\widehat{\rho }(𝐫,𝐩)`$ varies slowly in space over microscopic distances and expand the product of $`\widehat{\rho }`$ and $`f`$ in (65) and (67) according to

$`\widehat{\rho }(𝐫_1,𝐩_1^{\prime \prime })f(𝐫_1,𝐩_2^{\prime \prime }){\displaystyle \frac{𝐫}{2}}_{𝐫_1}\left[\widehat{\rho }(𝐫_1,𝐩_1^{\prime \prime })f(𝐫_1,𝐩_2^{\prime \prime })\right]`$
$`+{\displaystyle \frac{𝐫^{}}{2}}\left[\widehat{\rho }(𝐫_1,𝐩_1^{\prime \prime })_{𝐫_1}f(𝐫_1,𝐩_2^{\prime \prime })\left(_{𝐫_1}\widehat{\rho }(𝐫_1,𝐩_1^{\prime \prime })\right)f(𝐫_1,𝐩_2^{\prime \prime })+\mathrm{}\right]`$ (68)
The first term in (68) corresponds to the local term, the two which follow to first order in gradients non-local terms. A very similar equation exists for the product of $`\widehat{\rho }`$’s in (65) and (67).

##### (1) Local term

We now calculate the zeroth order in gradients. The following three integrals occur:

$`{\displaystyle d^3re^{i𝜿𝐫}}=(2\pi )^3\delta (𝜿);{\displaystyle d^3r^{}e^{i(𝐤_2^{}𝐤_1^{})𝐫^{}}}=(2\pi )^3\delta (𝐤_2^{}𝐤_1^{})`$
$`{\displaystyle d^3r^{}e^{i(𝐤_2^{}+𝐤_1^{})𝐫^{}}}=(2\pi )^3\delta (𝐤_2^{}+𝐤_1^{})`$ (69)
The linear in $`T`$ matrix terms become:

$`\widehat{I}_T(𝐫_1,𝐩_1)`$ $`={\displaystyle \frac{(2\pi )^4}{2i\mathrm{\Delta }t}}{\displaystyle d^3qd^3k_1^{}\delta (E_kE_{k_1^{}})}`$
$`\times [\delta (𝐤_1^{}𝐤)T(𝐤,𝐤_1^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤_1^{})f(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤_1^{})`$
$`+ϵ\delta (𝐤_1^{}𝐤)T(𝐤,𝐤_1^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤_1^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤_1^{})`$
$`+ϵ\delta (𝐤_1^{}+𝐤)T(𝐤,𝐤_1^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤_1^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤_1^{})`$
$`+\delta (𝐤_1^{}+𝐤)T(𝐤,𝐤_1^{})f(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤_1^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤_1^{})]`$
$`+\text{c.c.}`$ (70)
In the preceding equation, we note the appearance of the square of delta functions of the energy $`\delta (E_k)\delta (𝐤)\left[\delta (E_k)\right]^2`$. They are handled by the following well-known simplification from scattering theory:

$$\left[\delta (E_k)\right]^2=\frac{\mathrm{\Delta }t}{2\pi \mathrm{}}\delta (E_k)$$
(71)
where $`\mathrm{\Delta }t`$ is a time larger than the duration of a collision (see Ref. , for example), which simplifies with the one introduced in the $`S`$ matrix Ansatz (18). The same manipulations can be done with the quadratic in $`T`$ matrix term, they lead to:

$`\widehat{I}_{T^2}(𝐫_1,𝐩_1)`$ $`={\displaystyle \frac{(2\pi )^4}{2\mathrm{}}}{\displaystyle d^3qd^3k^{}\delta (E_kE_k^{})}`$
$`\times [T(𝐤,𝐤^{})T(𝐤,𝐤^{})^{}\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤^{})f(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤^{})`$
$`+ϵT(𝐤,𝐤^{})T(𝐤,𝐤^{})^{}\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤^{})`$
$`+ϵT(𝐤,𝐤^{})T(𝐤,𝐤^{})^{}\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤^{})`$
$`+T(𝐤,𝐤^{})T(𝐤,𝐤^{})^{}f(𝐫_1,𝐩_1\mathrm{}𝐤+\mathrm{}𝐤^{})\widehat{\rho }(𝐫_1,𝐩_1\mathrm{}𝐤\mathrm{}𝐤^{})]`$ (72)
We will now use different properties of the $`T`$ matrix. In addition to the differential $`\sigma _k(\widehat{k},\widehat{k}^{})`$ and total $`\sigma _T(k)`$ cross sections defined in (12), we introduce, following Ref. , the following cross sections:

$`T(𝐤,𝐤)`$ $`={\displaystyle \frac{\mathrm{}^2k}{i8\pi ^3m}}\left(\sigma _{fwd.}^{ex.}(k)i\tau _{fwd.}^{ex.}(k)\right)`$
$`T(𝐤^{},𝐤)T(𝐤^{},𝐤)^{}`$ $`={\displaystyle \frac{\mathrm{}^4}{4\pi ^4m^2}}\left(\sigma _k^{ex.}(\widehat{k},\widehat{k}^{})i\tau _k^{ex.}(\widehat{k},\widehat{k}^{})\right)`$ (73)
Unitarity of the $`S`$ matrix implies the optical theorem:

$$\frac{T(𝐤,𝐤)}{2i}+c.c.=\mathrm{Im}T(𝐤,𝐤)=\frac{\mathrm{}^2k}{(2\pi )^3m}\sigma _T(k)$$
(74)
and the rotational invariance of the interaction Hamiltonian can be used to show that:

$$T(𝐤,𝐤^{})=T(𝐤,𝐤^{})$$
(75)
Performing the integral over the length of wavevector $`𝐤^{}`$ in (70) and (72) and using the previous properties of the $`T`$ matrix, we obtain the Wigner transform to zero order of the $`S`$ matrix Ansatz $`\widehat{I}_W=\widehat{I}_T+\widehat{I}_{T^2}`$:

$`\widehat{I}_W(𝐫_1,𝐩_1)={\displaystyle }d^3q{\displaystyle \frac{q}{m}}\{[\sigma _T(k)\widehat{\rho }(𝐫_1,𝐩_1)f(𝐫_1,𝐩_2)`$
$`+{\displaystyle \frac{ϵ}{2}}(\sigma _{fwd.}^{ex.}(k)[\widehat{\rho }(𝐫_1,𝐩_1),\widehat{\rho }(𝐫_1,𝐩_2)]_++i\tau _{fwd.}^{ex.}(k)[\widehat{\rho }(𝐫_1,𝐩_1),\widehat{\rho }(𝐫_1,𝐩_2)]_{})]`$
$`{\displaystyle }d^2\widehat{k}^{}[\sigma _k(\widehat{k},\widehat{k}^{})\widehat{\rho }(𝐫_1,𝐩_1^{})f(𝐫_1,𝐩_2^{})+{\displaystyle \frac{ϵ}{2}}(\sigma _k^{ex.}(\widehat{k},\widehat{k}^{})[\widehat{\rho }(𝐫_1,𝐩_1^{}),\widehat{\rho }(𝐫_1,𝐩_2^{})]_+`$
$`+i\tau _k^{ex.}(\widehat{k},\widehat{k}^{})[\widehat{\rho }(𝐫_1,𝐩_1^{}),\widehat{\rho }(𝐫_1,𝐩_2^{})]_{})]\}`$ (76)
with the notations (11). This local term is the LL “collision integral” . It contains terms linear in $`T`$ matrix, usually written on the l.h.s. of the kinetic equation, as well as terms quadratic in $`T`$ matrix.

The low-energy limit of the cross sections is obtained from the low-energy expression of the $`T`$ matrix (5):

$$\begin{array}{cc}\sigma _T(k)4\pi a^2& \sigma _k(\widehat{k},\widehat{k}^{})a^2\\ \sigma _{fwd.}^{ex.}(k)4\pi a^2& \sigma _k^{ex.}(\widehat{k},\widehat{k}^{})a^2\\ \tau _{fwd.}^{ex.}(k)4\pi a/k& \tau _k^{ex.}(\widehat{k},\widehat{k}^{})0\end{array}$$

when $`k0`$, see Ref. . When introduced into (76), the collision integral of equation (22) is obtained, as well as the spin mean field contained in the commutator in equation (20). Other mean field terms (contained in the anticommutator) appear to first order of the gradient expansion.

##### (2) First order terms

We retain only first order terms that are in addition linear in $`T`$ matrix.

###### (2-a) $`𝐫`$ gradients

The first-order terms introduce the gradient of a delta function:

$$d^3r𝐫\text{ }e^{i𝜿𝐫}=i(2\pi )^3_𝜿\delta (𝜿)$$

which implies taking the derivative with respect to $`𝜿`$ of the function under the integral. We obtain:

$`{\displaystyle \frac{(2\pi )^3}{2}}`$ $`{\displaystyle }d^3q\{\mathrm{Re}T(𝐤,𝐤)\text{ }_𝐩_𝐫\left[\widehat{\rho }(𝐫_1,𝐩_1)f(𝐫_1,𝐩_2)\right]`$
$`_𝐩\left[\mathrm{Re}T(𝐤,𝐤)\right]\text{ }_𝐫\left[\widehat{\rho }(𝐫_1,𝐩_1)f(𝐫_1,𝐩_2)\right]`$
$`+{\displaystyle \frac{ϵ}{2}}T(𝐤,𝐤)\text{ }_𝐩_𝐫\left[\widehat{\rho }(𝐫_1,𝐩_1)\widehat{\rho }(𝐫_1,𝐩_2)\right]`$
$`{\displaystyle \frac{ϵ}{2}}_𝐩\left[T(𝐤,𝐤)\right]_𝐫\left[\widehat{\rho }(𝐫_1,𝐩_1)\widehat{\rho }(𝐫_1,𝐩_2)\right]`$
$`+{\displaystyle \frac{ϵ}{2}}T(𝐤,𝐤)^{}\text{ }_𝐩_𝐫\left[\widehat{\rho }(𝐫_1,𝐩_2)\widehat{\rho }(𝐫_1,𝐩_1)\right]`$
$`{\displaystyle \frac{ϵ}{2}}_𝐩\left[T(𝐤,𝐤)^{}\right]_𝐫\left[\widehat{\rho }(𝐫_1,𝐩_2)\widehat{\rho }(𝐫_1,𝐩_1)\right]\}`$ (77)
###### (2-b) $`𝐫^{}`$ gradients

The following integrals occur:

$`{\displaystyle d^3r^{}e^{i(𝐤_2^{}𝐤_1^{})𝐫^{}}𝐫^{}}=i(2\pi )^3_{𝐤_1^{}}\delta (𝐤_1^{}𝐤_2^{})`$
$`{\displaystyle d^3r^{}e^{i(𝐤_2^{}+𝐤_1^{})𝐫^{}}𝐫^{}}=i(2\pi )^3_{𝐤_1^{}}\delta (𝐤_2^{}+𝐤_1^{})`$ (78)
so that we obtain:

$`{\displaystyle \frac{(2\pi )^3}{2}}`$ $`{\displaystyle }d^3q_𝐩\{\mathrm{Re}T(𝐤,𝐤)\text{ }[\widehat{\rho }(𝐫_1,𝐩_1)_𝐫f(𝐫_1,𝐩_2)f(𝐫_1,𝐩_2)_𝐫\widehat{\rho }(𝐫_1,𝐩_1)]`$
$`+{\displaystyle \frac{ϵ}{2}}T(𝐤,𝐤)\text{ }\left[\widehat{\rho }(𝐫_1,𝐩_1)_𝐫\widehat{\rho }(𝐫_1,𝐩_2)\widehat{\rho }(𝐫_1,𝐩_2)_𝐫\widehat{\rho }(𝐫_1,𝐩_1)\right]`$
$`+{\displaystyle \frac{ϵ}{2}}T(𝐤,𝐤)^{}\text{ }[_𝐫\widehat{\rho }(𝐫_1,𝐩_2)\text{ }\widehat{\rho }(𝐫_1,𝐩_1)_𝐫\widehat{\rho }(𝐫_1,𝐩_1)\text{ }\widehat{\rho }(𝐫_1,𝐩_2)]\}`$ (79)
The expression of the $`T`$ matrix to lowest order in the scattering length (see equation (5))

$$(2\pi )^3T(𝐤,𝐤)g$$

is independent of $`𝐤`$. Using the preceding expression in (77) and (79), we obtain the following first order term:

$$\frac{1}{2}[_𝐩\widehat{\rho }(𝐫_1,𝐩_1),_𝐫(gn(𝐫_1)\widehat{1}+ϵg\widehat{n}(𝐫_1))]_+$$

This is the mean field term appearing in the anticommutator in the kinetic equation (20).

## Appendix II: spin-dependent interactions (forward scattering)

In this appendix, we show how to obtain the full mean field in the case of spin-dependent interactions, using the $`S`$-matrix Ansatz. We will limit ourselves to collisions at low-energy, using the $`T`$ matrix only at lowest order in the scattering lengths (equation (23)):

$$\alpha ;\beta |\widehat{T}_k|\gamma ;\delta =\delta _{\alpha ,\gamma }\delta _{\beta ,\delta }\frac{g_{\alpha \beta }^{(d)}}{(2\pi )^3}+(1\delta _{\alpha ,\beta })\delta _{\alpha ,\delta }\delta _{\beta ,\gamma }\frac{g_{\alpha \beta }^{(t)}}{(2\pi )^3}$$
(80)
where $`g_{\alpha \alpha }^{(d)}g_{\alpha \alpha }`$. In this limit the $`T`$ matrix elements are real.

As we are interested in mean field terms, we only keep terms linear in the $`T`$ matrix in the r.h.s. of the $`S`$ matrix Ansatz, equation (18):

$$\frac{2\pi }{i\mathrm{\Delta }t}Tr_2\left\{\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\left[\widehat{𝒯}\widehat{\rho }_1(1)\widehat{\rho }_1(2)\widehat{\rho }_1(1)\widehat{\rho }_1(2)\widehat{𝒯}\right]\frac{\widehat{1}+ϵ\widehat{P}_{ex.}}{\sqrt{2}}\right\}$$
(81)
where we introduced the $`𝒯`$ matrix (the on-shell $`T`$ matrix), whose definition is given in Appendix I, equation (58). Using the properties of the trace on particle $`2`$ and the fact that the exchange operator commutes with the $`𝒯`$ matrix, we can rewrite the preceding equation in the form of a commutator:

$$\frac{1}{i\mathrm{}}[\frac{2\pi \mathrm{}}{\mathrm{\Delta }t}Tr_2\left\{(\widehat{1}+ϵ\widehat{P}_{ex.})\widehat{𝒯}\widehat{\rho }_1(2)\right\},\widehat{\rho }_1(1)]_{}$$
(82)
It seems therefore natural that the first operator in the commutator should play the role of an effective single-particle Hamiltonian.

We now consider a general matrix element of the preceding commutator. To simplify the notation, we only write the first term of the commutator:

$$1:𝐩_1,\alpha |\frac{2\pi \mathrm{}}{\mathrm{\Delta }t}Tr_2\left\{(\widehat{1}+ϵ\widehat{P}_{ex.})\widehat{𝒯}\widehat{\rho }_1(2)\right\}\widehat{\rho }_1(1)|1:𝐩_1^{},\alpha ^{}$$
(83)
We calculate this matrix element by introducing three closure relations of the form:

$$\widehat{1}=d^3p\underset{\beta }{}|𝐩,\beta 𝐩,\beta |$$
(84)
and by using the expression (80) of the $`T`$ matrix. We obtain:

$`{\displaystyle \frac{2\pi \mathrm{}}{(2\pi )^3\mathrm{\Delta }t}}`$ $`\times {\displaystyle }d^3p_2{\displaystyle }d^3p_3\delta (E_{(𝐩_1𝐩_2)/2}E_{(𝐩_3𝐩_4)/2})`$
$`\times {\displaystyle \underset{\beta }{}}\{g_{\alpha \beta }^{(d)}𝐩_3,\alpha |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}𝐩_4,\beta |\widehat{\rho }_1|𝐩_2,\beta `$
$`+(1\delta _{\alpha ,\beta })g_{\alpha \beta }^{(t)}𝐩_3,\beta |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}𝐩_4,\alpha |\widehat{\rho }_1|𝐩_2,\beta `$
$`+ϵg_{\alpha \beta }^{(d)}𝐩_3,\beta |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}𝐩_4,\alpha |\widehat{\rho }_1|𝐩_2,\beta `$
$`+ϵ(1\delta _{\alpha ,\beta })g_{\alpha \beta }^{(t)}𝐩_3,\alpha |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}𝐩_4,\beta |\widehat{\rho }_1|𝐩_2,\beta \}`$ (85)
where $`𝐩_4=𝐩_1+𝐩_2𝐩_3`$. By defining the three following coupling constants:

$$g_{\alpha \alpha }=g_{\alpha \alpha }^{(d)};g_{\alpha \beta }=g_{\alpha \beta }^{(d)}+ϵg_{\alpha ,\beta }^{(t)}\text{if }\alpha \beta $$
(86)
we can rewrite equation (85) as:

$`{\displaystyle \frac{2\pi \mathrm{}}{\mathrm{\Delta }t}}`$ $`\times {\displaystyle }d^3p_2{\displaystyle }d^3p_3\delta (E_{(𝐩_1𝐩_2)/2}E_{(𝐩_3𝐩_4)/2})`$
$`\times {\displaystyle \underset{\beta }{}}{\displaystyle \frac{g_{\alpha \beta }}{(2\pi )^3}}\{𝐩_4,\beta |\widehat{\rho }_1|𝐩_2,\beta 𝐩_3,\alpha |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}`$
$`+ϵ𝐩_4,\alpha |\widehat{\rho }_1|𝐩_2,\beta 𝐩_3,\beta |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}\}`$ (87)
We therefore see that, for both statistics (bosonic or fermionic), only three coupling constants ($`g_{11}`$, $`g_{22}`$ and $`g_{12}`$) are involved.

Now, taking the Wigner transform of the operator (82) and expanding in gradients (see Appendix I) introduces another delta function of energy conservation, which is present at each order in gradients. This allows to use the “$`\delta (E)^2`$” simplification (equation (71) in Appendix I) to make the $`\mathrm{\Delta }t`$ disappear. The net result is that:

$$\frac{2\pi \mathrm{}}{\mathrm{\Delta }t}\delta (E_{(𝐩_1𝐩_2)/2}E_{(𝐩_3𝐩_4)/2})$$
(88)
gets replaced by $`1`$. Formally the matrix element of equation (83) becomes:

$`{\displaystyle d^3p_2d^3p_3\underset{\beta }{}\frac{g_{\alpha \beta }}{(2\pi )^3}}`$ $`\times \{𝐩_4,\beta |\widehat{\rho }_1|𝐩_2,\beta 𝐩_3,\alpha |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}`$
$`+ϵ𝐩_4,\alpha |\widehat{\rho }_1|𝐩_2,\beta 𝐩_3,\beta |\widehat{\rho }_1|𝐩_1^{},\alpha ^{}\}`$ (89)
This result can be rewritten:

$$1:𝐩_1,\alpha |\widehat{V}^{mf}(1)\widehat{\rho }_1(1)|1:𝐩_1^{},\alpha ^{}$$
(90)
where:

$$V_{\alpha \beta }^{mf}(𝐫)=\delta _{\alpha \beta }\underset{\gamma =1,2}{}g_{\alpha \gamma }𝐫,\gamma |\widehat{\rho }_1|𝐫,\gamma +ϵg_{\alpha \beta }𝐫,\alpha |\widehat{\rho }_1|𝐫,\beta $$
(91)
is the mean field potential. Therefore, equation (82) can be written formally:

$$\frac{1}{i\mathrm{}}[\widehat{V}^{mf}(1),\widehat{\rho }_1(1)]_{}$$
(92)
This equation is valid provided it is used to compute the Wigner transform and do a subsequent gradient expansion.

In the body of the article, we define the effective single-particle Hamiltonian $`\widehat{U}`$ which is the sum of the external potential $`\widehat{V}^{ext}`$ and of the mean field potential $`\widehat{V}^{mf}`$. The operatorial kinetic equation is now:

$$\frac{d\widehat{\rho }_1}{dt}+\left(i\mathrm{}\right)^1[\widehat{\rho }_1,\widehat{H}_1+\widehat{V}^{mf}]_{}=0$$
(93)
The kinetic equation is then derived following the procedure of Appendix I (Wigner transform and gradient expansion up to first order); the result is:

$`_t\widehat{\rho }(𝐫,𝐩)`$ $`+{\displaystyle \frac{𝐩}{m}}_𝐫\widehat{\rho }(𝐫,𝐩)+{\displaystyle \frac{1}{i\mathrm{}}}[\widehat{\rho }(𝐫,𝐩),\widehat{U}(𝐫)]_{}`$
$`{\displaystyle \frac{1}{2}}[_𝐩\widehat{\rho }(𝐫,𝐩),_𝐫\widehat{U}(𝐫)]_+=0`$ (94)
which is equal to the l.h.s. of (25).