1. Introduction
When the average collision time in ultracold dilute gases made of bosonic atoms is much larger than the relevant dynamics characteristic time scale, it is possible to have a model based on the Landau–Vlasov equation [
1]. The Landau–Vlasov equation is obtained from the Boltzmann–Vlasov equation [
1,
2,
3,
4] neglecting the collision operator. The dynamics of ultracold bosonic systems, e.g., in the crossover from collisionless to collisional regimes needs the Boltzmann–Vlasov equation [
5]. Hydrodynamic equations [
6,
7] are useful tools in the collisional case, for instance, for Bose–Einstein condensates [
6] or the superfluid Fermi gas in the BCS–BEC crossover [
8].
Under a very strong transverse confinement, a bosonic gas is in a quasi one-dimensional (1D) configuration. Experimental achievement of quasi-1D systems is realized in ultracold atoms trapped in optical potentials with harmonic transverse confinement energies much larger than the temperature or chemical potential [
9]. The collisionless regime is enhanced in the quasi-1D configuration. Indeed, in 1D binary elastic collisions, particles exchange their energies completely, thus there is no sensible effects from these collisions between identical particles. Consequently, no thermalization is possible, as verified in ultracold bosonic atoms trapped in 1D optical lattices [
10,
11]. For these dilute 1D bosonic systems, the Landau–Vlasov equation is applicable, provided the gas does not contain a quasi-condensate and that it is not in the Tonks–Girardeau regime, with fermionic properties [
10,
12]. We are following the terminology of ultracold atoms community [
12] (and references therein) when referring to the Landau–Vlasov equation. Namely, it is collisionless so that it has no “Landau collision operator”, as would be implied in the context of plasma physics.
Recently [
12], the linear stability of solutions of the 1D Landau–Vlasov equation was investigated by means of well-known methods from plasma physics, namely, the Landau or Laplace transform approach. In this method, the time-evolution of perturbations around the equilibrium distribution function is treated as an initial-value problem. The linear Landau damping rate (or growth rate, for unstable equilibria) is therefore determined upon the adequate analysis in the complex plane (Landau contour). The similarity between the Landau–Vlasov equation and the Vlasov–Poisson system describing collisionless electrostatic plasmas provides a stimulating scenario for the application of plasma techniques in a seemingly uncorrelated area such as in the study of ultracold bosonic gases.
In this context, the present work is dedicated to the discussion of time-independent solutions and stationary wave solutions for the 1D Landau–Vlasov equation. In plasmas, stationary solutions for the Vlasov–Poisson system can be derived starting from Jeans’s theorem according to which the particle distribution function satisfying Vlasov’s equation should be a function of the constants of motion. In the time-independent case, the particle energy is such a constant of motion or invariant, as treated in the original work [
13] by Bernstein, Greene, and Kruskal (BGK). By construction, these so-called BGK modes are exact nonlinear plasma oscillations which do not present damping or growth. The BGK approach where the energy is the central dynamical variable can be adapted for the derivation of phase-space hole structures [
14,
15,
16,
17] and, to a more limited extent, to quantum plasmas [
18].
In spite of the more popular view in terms of the surfing electron interpretation [
19], an alternative, more rigorous interpretation of Landau damping is in terms of the phase mixing superposition of Case–Van Kampen modes [
20]. Introduced by Van Kampen [
21] and demonstrated by Case [
22] to form a complete orthogonal set for the linearized Vlasov–Poisson system, the stationary wave or Case–Van Kampen modes have been discussed in a variety of contexts. For instance, in plasmas with an ionic background slowly varying in time [
23]; in multidimensional non-uniform plasmas [
24]; for nonlinear waves [
25]; extended Fermi systems [
26]; and electromagnetic [
27], collisional [
28], and quantum [
29] plasmas. As discussed in
Section 3, the Case–Van Kampen modes are also available for an ultracold boson gas described by the Landau–Vlasov equation.
This work is organized as follows. In
Section 2, we revisit the 1D Landau–Vlasov equation, which was derived and discussed in detail in [
12].
Section 3 considers BGK modes and
Section 4 the Case–Van Kampen modes for the 1D Landau–Vlasov equation.
Section 5 is reserved to the conclusions and final remarks.
2. The One-Dimensional Landau–Vlasov Equation
The one-dimensional (1D) Landau–Vlasov equation is given [
12] by
where
is the 1D probability distribution function,
is the external confinement potential,
is the axial local density and
is the renormalized 1D interaction strength, in which
ℏ is the reduced Planck constant,
m is the atomic mass,
is the s-wave scattering length of the interaction between atoms, and
is the characteristic transverse width occupied by the dilute bosonic gas. The integrals are taken from minus to plus infinity except if explicitly stated. The normalization
is adopted, where
N is the total number of bosons. In some cases we will use the harmonic external potential
although this choice is not decisive for the following treatment.
The validity conditions of the 1D Landau–Vlasov equation are
where
is the energy associated with the transverse confinement,
and
Equations (
6)–(
8) are resp. Equations (10), (18) and (19) in [
12], assuring a quasi-1D configuration where the ultracold dilute bosonic gas is neither in a quasi-condensate or Tonks–Girardeau regime.
3. Bernstein–Greene–Kruskal Modes
In the stationary case where
everywhere, the general solution to Equation (
1) is
where
f is an arbitrary function of the energy function
with the total potential
This holds for arbitrary external potential, as long as it is time-independent. The same reasoning applies to the BGK solution for the stationary Vlasov–Poisson system, with some differences. The energy function in the plasma problem contains the electrostatic potential, while in the bosons problem
H depends on the particle distribution function itself, through the interaction potential
where
is a functional of
f, viz. Equation (
2). Moreover, there is nothing similar to Poisson’s equation to be self-consistently solved, but only the normalization condition (
4). In this context, therefore, it is not an exaggeration to consider the stationary Vlasov–Landau equation to be much simpler than the stationary Vlasov–Poisson system. Nevertheless, concrete applications require a detailed analysis, as shown in the next examples.
3.1. Maxwell–Boltzmann Distribution
The functional form of
is entirely free, which is in accordance with the collisionless assumption so that no particular equilibrium (e.g., the Bose–Einstein distribution) is preferred. Suppose there is a Maxwell–Boltzmann distribution
where
has the role of inverse temperature in energy units,
A is a normalization constant to be determined and
is a numerical factor included for convenience.
From Equations (
2) and (
10)–(
12), the 1D particle number density is
The axial local density appears in both sides of Equation (
13). Nevertheless, in this example the determining equation can be easily disentangled according to
where the Lambert
W function or product log function is defined [
30] as the solution of
, in the domain
. By construction, the solution is analytically exact. The exact total potential (
11) is also entirely available. For simplicity, a repulsive interaction (
) is assumed, so that
in Equation (
14) is automatically a real, positive definite quantity.
The last step concerns the determination of the normalization constant
A. For instance, for the harmonic potential in Equation (
5), it is convenient to introduce the rescaled variables
in terms of the characteristic length
. The normalization condition (
4) yields
which cannot be analytically solved for
. Nevertheless, given the rescaled coupling constant
one can readily numerically obtain
, as shown in
Figure 1. It should be remarked that
has very small values in today’s experiments [
12,
31], which allows to approximate
for a generic argument
in Equation (
16), yielding
in this approximation. For instance, with
(Li atom) as in [
12] together with typical values
[
32,
33] one has
. The validity conditions (
7) and (
8) are safely meet for these parameters. Moreover, from Equation (
6) one would need an energy of transverse confinement
with
.
It is interesting to rewrite the validity condition (
7) using the approximations
and
so that
. From Equation (
13), one has the estimate
. Together with
, one has that Equation (
7) becomes
which has an evident thermodynamic meaning. Under the same approximation, the Lambert function in Equation (
14) can be safely replaced by
so that the number density assumes the Maxwellian form
To summarize and without any approximation within the Landau–Vlasov model, for the Maxwell–Boltzmann distribution (
12) one has the exact number density (
14), subject to
which determines the normalization constant
A given an arbitrary external potential
.
3.2. Water Bag Distribution
In a non-equilibrium situation we are free to have any function of the total energy as a suitable particle distribution function. A second example is provided by a completely degenerate Fermi–Dirac-like distribution
where
is the step function,
A is a normalization constant, and
is a energy parameter which would be the Fermi energy in a Fermi gas. Moreover, we assume
; otherwise, some quantities become complex valued in the following. However, in the context of a bosonic gas,
is just a measure of the energy spread, precisely as in the water bag model for plasmas [
34]. By construction, Equation (
19) shows an exact stationary solution of the Landau–Vlasov equation.
Integration in momentum space implicitly gives the 1D number density
Note for real
one has
and hence automatically
, supposing
(repulsive interaction). The always non-negative solution of Equation (
20) is
The exact total potential (
11) is also immediately available, for arbitrary external potential.
The last step is the determination of the normalization constant
A, once a specific external potential is chosen. For the harmonic potential (
5), it is convenient to introduce the rescaling
in terms of the characteristic length
. The dimensionless 1D number density becomes
and the normalization condition
yields
where
yields
in dimensionless variables. The integral in Equation (
24) can be analytically done, but in terms of transcendental functions which are not useful to show here. For specific values of
, one can numerically determine
and hence the necessary normalization, as depicted in
Figure 2. In the limit of very small
one has
. The resulting dimensionless 1D number density is shown in
Figure 3.
Similarly to the Maxwellian case, it is possible to rewrite the validity condition (
7) using the approximations
and
. From Equation (
20) one has the estimate
. Together with
yields
with an evident thermodynamic meaning.
To summarize and without any approximation within the Landau–Vlasov model, for the water bag, completely degenerate Fermi–Dirac-like distribution (
19) one has the exact number density (
21), subject to
which determines the normalization constant
A given an arbitrary external potential
. For the sake of introducing the next section, notice the correspondence between the Van Kampen mode decomposition and the water bag distribution which is shown to be granted in the limit of an infinite number of bags [
35].
4. Case–Van Kampen Modes
The Case–Van Kampen modes are the normal modes in a plasma [
21,
22]. For the Landau–Vlasov equation, Case–Van Kampen modes can be derived starting from the assumption
. Experimentally, a 1D configuration of ultracold atoms free of external confinement can be produced by means of a 1D ring geometry with large radius [
11] or considering two high barriers at the edges of a straight axial arrangement [
36]. In both cases, spatial periodic conditions can be assumed. The 1D Landau–Vlasov equation reduces to
To proceed, we set
where
is the equilibrium distribution subject to
where
is the equilibrium 1D number density and
is a first-order perturbation. Due to the periodic boundary conditions, it is meaningful to Fourier transform according to
where
k is a multiple of a fundamental wavenumber. Linearizing the Landau–Vlasov Equation (
26) the result is
Complementary to the Landau approach where the linearized Landau–Vlasov equation is taken as an initial value problem analyzed by Laplace transform methods, the Case–Van Kampen approach assumes
where
is a real, arbitrary constant and
is the phase speed. Inserting from Equation (
31) into Equation (
30) gives the eigenvalue problem
Following [
21,
22], it is convenient to assume the normalization
So that the integral Equation (
32) simplifies to
In a distributional sense, the solution for Equation (
34) is
where
℘ denotes the Cauchy principal value symbol,
is a function to be determined and where
is the Dirac delta.
As can be readily verified, the naive solution (without principal value and with
) can not be made compatible with the normalization (
33). Indeed, to comply with the normalization one needs
where the integral is taken in the principal value sense. The final result is
As apparent, these Case–Van Kampen modes are not damped (they are stationary waves) and have a singular character. Following Case, by means of the introduction of adequate adjoint solutions, it is possible to demonstrate that the Equation (
37) provides a complete set, in the sense that all solutions to the linearized Landau–Vlasov equation can be expressible as a linear combination of these modes. More precisely, for simplicity we have discussed only the class 1a among the four classes of eigenfunctions in Case’s terminology, as detailed in the original article [
22] and textbooks [
37]. To demonstrate the completeness of the full set of Case–Van Kampen eigenfunctions makes necessary to introduce an auxiliary (or adjoint) equation with a different set of eigenfunctions, orthogonal to those in the original set, except when the eigenvalues coincide. The complete analysis is not trivial but entirely similar to the Vlasov–Poisson case, shown in [
22,
37], for instance.
Note that in spite of the singular character of the Case–Van Kampen modes, they can be used to compute well behaved physical quantities. For instance, we can consider the 1D number density perturbation
where
is an arbitrary weight function. Allowing a superposition law taking into account a frequency spread is valid in the context of the linearized Landau–Vlasov equation. Restricting to the
k-th Fourier component, applying Equations (
31) and (
37), we have the well behaved function
after using the property
. In this case, the first and last terms of Equation (
37) cancel upon integration (the same occurs for Vlasov–Poisson plasmas [
37]).
As a simple illustration, the Gaussian weight function
produces from integration of Equation (
39) the density perturbation
This is an example of the well known fact that although the isolated Case–Van Kampen eigenmodes are stationary waves, they can produce damped macroscopic objects, taking into account phase mixing. Consistently, the monochromatic limit is not damped.