1. Introduction
The Boltzmann transport equation (BTE) [
1,
2,
3,
4] furnishes a semi-classical description of the effects of particle transport, including the influence of external forces on the effective single-particle distribution function
of a small cell in phase phase, centered at position
r and momentum
p. For a system with identical particles of mass
m, the Boltzmann equation reads:
Here,
is the number of particles in a cell of phase volume
, centered at position
r and momentum
p [
3]. In addition,
is the force field acting on the particles in the fluid. The term on the right-hand side is added to describe the effect of collisions between particles. It is a statistical term and requires knowledge of the statistics that the particles obey, like the Maxwell–Boltzmann, Fermi–Dirac or Bose–Einstein distributions. In his famous “Stoßzahlansatz” (or hypothesis of molecular chaos), Boltzmann obtained an explicit form for it. In modern notation, for example for an interacting Fermi gas, where a particle from a state with momentum
p is scattered to a state with momentum
, whereas a second particle is scattered from a momentum
q to a momentum
, the collision term reads:
where
is the normalized transition probability from the two-particle state with momenta
to the state labeled by
. Clearly, solving this widely-studied equation is a very difficult task. It might be hoped that symmetries could be helpful. The equation without the collision term is known as the Vlasov equation [
5]. The relationship with Landau damping and a physicists’ derivation can be found in [
6,
7]. In this work, we shall explore a class of symmetries of the (collisionless) BTE.
Throughout, we shall restrict to the
space dimension. (By analogy with other constructions of local scale symmetries (see [
8,
9,
10,
11] and especially [
12] and references therein), we expect a straightforward extension of the results reported here to
. Since we shall construct here a finite-dimensional Lie algebra of dynamical conformal symmetries of the
collisionless BTE, one should indeed expect that an extension to
exists. That symmetry algebra should contain three generators
, along with a vector of generators
and also spatial rotations.) We start from a non-standard representation, isomorphic to the infinite-dimensional Lie algebra of conformal transformations in
dimensions. (For the sake of clarity, we shall adopt the following convention of terminology: the infinite-dimensional Lie algebra
will be called a (centerless) “conformal Virasoro algebra”. Its maximal finite-dimensional sub-algebra
will be called a “conformal algebra”) This Lie algebra is spanned by the generators
and can be defined from the commutators [
9,
12]:
where μ is a parameter. An explicit realization in terms of time-space transformation is [
9,
12]:
such that
can be interpreted as a velocity (“speed of light/sound”) and where
are constants. (The contraction
of Equation (
3) produces the non-semi-simple “altern-Virasoro algebra”
(but without central charges). Its maximal finite-dimensional sub-algebra is the conformal Galilean algebra
cga(1) [
9,
13]; see also [
8,
11]. The
cga(d) is non-isomorphic to either the standard Galilei algebra or else the Schrödinger algebra.) Writing
and
, where the generators
satisfy
, it can be seen that, provided
, the above Lie algebra Equation (
2) is isomorphic to a pair of Virasoro algebras
with a vanishing central charge. However, this isomorphism does not imply that physical systems described by two different representations of the conformal Virasoro algebra, or the conformal algebra, with commutators Equation (
2), were trivially related. For example, it is well known that if one uses the generators of the standard representation of conformal invariance or else the non-standard representation Equation (
4) in order to find co-variant two-point functions, the resulting scaling forms are different [
9].
Now, consider the maximal finite-dimensional sub-algebra
, which for
, in turn, is isomorphic to the direct sum
. The explicit realization follows from from Equation (
3):
Here, the generators
describe time- and space-translations,
is a (conformal) Galilei transformation (since the commutator
does not vanish and does not give a central element of the Lie algebra Equation (
2), its structure is fundamentally different from algebras containing the usual Galilei algebra as a sub-algebra),
gives the dynamical scaling
of
(with
), such that the so-called “dynamical exponent”
, since both time and space are re-scaled in the same way, and, finally,
give “special” conformal transformations. In the context of statistical mechanics of conformally-invariant phase transitions, one characterizes co-variant quasi-primary scaling operators through the invariant parameters
, where
x is the scaling dimension.
Finally, the finite-dimensional representation Equation (
4) acts as a dynamical symmetry on the equation of motion:
in the sense that a solution
ϕ of
is mapped onto another solution of the same equation. Indeed, it is easily checked that:
and
It follows that for fields
ϕ with scaling dimensions
, the algebra Equation (
4) really leaves the solution space of Equation (
5) invariant.
In order to return to the Boltzmann equation, we consider Equation (
5) in the form:
where
is interpreted as a single-particle distribution function and where we consider
v as an additional variable. Equation (
7) is a simple Boltzmann equation, without an external force, without a collision term and in one space dimension. From Equation (
6), with
v fixed (and normalized to
), its solution space is conformally invariant. (With respect to Equation (
5),
was replaced. This change must also be made in the generators Equation (
4) and commutators Equation (
2)). In
Section 2, we shall generalize the above representation of the conformal algebra to the situation with
v as a further variable. In
Section 3, we shall further extend this to the case when an external force
, possibly depending on time, spatial position and velocity, is included. The aim of these calculations is to determine which situations of potential physical interest with a non-trivial conformal symmetry might be identified. This explorative study aims at identifying lines for further study, which might lead later to a more comprehensive understanding of the possible symmetries of Boltzmann equations. Taking into account the collision term is left for future work. We shall concentrate on the
space dimension throughout. Conclusions and final comments are given in
Section 4.
2. Collisionless Boltzmann Equation without External Forces
In our construction of conformal dynamical symmetries of the 1D collisionless BTE, we shall often meet Lie algebras of a certain structure. These will be isomorphic to the two-dimensional conformal algebra.
Proposition 1. The Lie algebra defined by the commutators:where are constants, is isomorphic to the pair of centerless Viraso algebras .
Proof. For either
or
, this is either evident or else has already been seen in
Section 1. In the other case, consider the change of basis
and
, where
are two families of commuting generators of
and α and β are constants, such that
. It then follows
and
. ☐
This implies in particular the isomorphism of the maximal finite-dimensional sub-algebras, or “conformal algebras” in the terminology chosen here. By definition, this “conformal algebra” obeys the commutators Equation (
8), but with
.
Our construction of dynamical symmetries of the Equation (
7) follows the lines of the construction of local scale-invariance in time-dependent critical phenomena [
9]. The physically-motivated requirements are: First of all it, is clear that the equation is invariant under time-translations:
Some kind of dynamical scaling must be present, as well. Its most general form is:
Whenever the dynamical exponent
, we shall find an explicit dependence on
v. In general, we look for a family of generators
, for which we make the ansatz:
We shall find
from the following three conditions (throughout, we use the notations
):
must be a symmetry for the Equation (
7); hence,
. This gives:
The generator
is assumed to be in the Cartan sub-algebra; hence,
. It follows:
The action of
is as a lowering operator; hence,
. It follows:
These conditions, combined with the following initial conditions:
must be sufficient for the determination of all admissible forms of
.
In the special case
, we have
and find the most general form of
as a symmetry of Equation (
7) as follows (the requirement that
close into the Lie algebra
fixes
):
and
with a certain set of undetermined constants.
For conformal invariance, a family of generators
must also be found. Its construction is straightforward if the explicit form of
is known. Really,
must act as a raising operator, in both hierarchies, such that [
9]:
which implies that
. However, the usual realization of
as space translations does not work, since if we set all undetermined constants in Equation (
19) to zero, one would have
. It is better to work with the form:
as we shall do from now on.
We first consider the special case, when all of the constants in the expression Equation (
19) for
vanish:
Case A: A12 = A110 = A100 = B110 = B100 = D0 = 0.
Proposition 2. The six generators:span a representation of the conformal algebra Equation (2), which acts as dynamical symmetry algebra of the Equation (7), for arbitrary dynamical exponent z.
Proof. It is readily checked that the generator Equation (
26) satisfies the commutation relations (
2), with
. On the other hand, for any
, one has:
which establishes the asserted dynamical symmetry. ☐
Next, we treat the general case, when all of the constants are non-zero:
Case B: A12 ≠ 0, A110 ≠ 0, A100 ≠ 0, B110 ≠ 0, B100 ≠ 0, D0 ≠ 0.
Then, the generators are modified as follows:
Now, computing:
we conclude that the cases
and
must be treated separately.
Case B1: A12 = 0. It follows that the constants in Equation (
19) are given by:
Proposition 3. Let and be arbitrary constants. Then, the six generators:span a representation of the conformal algebra (the above result of Case A is recovered upon setting ). These generators give more symmetries of Equation (7). Proof. From the above, the commutator Equation (
2) is readily verified, with
. For the dynamical symmetries, one checks the commutators:
which proves the assertion. ☐
In contrast to the previous Case A, the representation acting only on , but keeping v as a constant parameter, can no longer be obtained by simply setting . Rather, one must set first, and only then, the limit is well-defined.
Case B2: A12 ≠ 0, A110 ≠ 0, B110 ≠ 0, A100 ≠ 0, B100 ≠ 0, D0 ≠ 0.
It turns out that for , the algebra also can be closed, but only if and (then, all other constants also vanish).
Proposition 4. Let z be an arbitrary constant. Then, the generators , where:close into a Lie algebra, with the following non-zero commutation relations:with . The algebra is isomorphic to the usual conformal algebra Equation (2) and further extends the dynamical symmetries of Equation (7). Proof. The commutation relation is directly verified. The isomorphism with the conformal algebra follows from Proposition 1. The requirement to have a symmetry algebra of Equation (
7) implies a relation between the constants
(called
in Proposition 1) and μ, namely
. In this case at hand, we have
,
. It is then verified that
and:
which proves that these are dynamical symmetries of (
7). ☐
We now ask whether the finite-dimensional representation Equations (
26), (30) and (31), with
, acting on functions
, and having a dynamical exponent
, can be extended to representations of an infinite-dimensional conformal Virasoro algebra. The answer turns out to be negative:
Proposition 5. The representation Equations (26), (30) and (31) of the finite-dimensional conformal algebra with commutator Equation (8) cannot be extended to representations of an infinite-dimensional conformal Virasoro algebra with commutator Equation (8) when .
Similar no-go results have been found before for variants of representations of the Schrödinger and conformal Galilean algebras [
14]. On the other hand, for
, extensions to a representation of a conformal Virasoro algebra with
exist [
15].
Proof. Since for the finite-dimensional representations Equations (
26), (30) and (31), we have:
we suppose that this must be valid for all admissible
. Now, using the condition Equation (
17) for
, a conformal Virasoro algebra should contain a new generator
. Starting from the most general form,
we find that the coefficients are obtained from:
where
are unknown functions of their arguments, but do no longer depend on the time
t. We want to satisfy
. However, when calculating:
we see that closure is not possible for
. Indeed, although the dependence on
of the functions
can be chosen to satisfy the above closure condition, the
t-dependence cannot be absorbed into these functions. Hence, our new representation Equations (
26), (30) and (31) of the conformal algebra Equation (
8) are necessarily finite-dimensional. ☐
3. Symmetry Algebra of Collisionless Boltzmann Equation with an Extra Force Term
We write the collisionless Boltzmann equation in the form:
We want to determine the admissible forms of an external force
, such that Equation (33) is invariant under a representation of the conformal algebra Equation (
8). The unknown representation must include the “force” term and, in particular, for
, it should coincide with the representations of conformal algebra obtained in the previous section.
The idea of the construction is similar to the one used in
Section 2. First, we impose invariance under basic symmetries:
From invariance under time translation
, it follows:
From invariance under dynamical scaling
, we obtain that:
if
satisfies the equation
, with solution:
where
is an arbitrary function of the scaling variable
.
It turns out that for the following calculations, it is more convenient to make a change of independent variables
. In the new variables, the generator of dynamical scaling just reads:
Next, in order to be specific, we make the following ansatz for the analogue of space translations. (Indeed, we might also require to find
from the conditions to be (i) a symmetry of Boltzmann equation and (ii) to form a closed Lie algebra with the other basic symmetries
. Such requirements lead to a system of differential equations, and the ansatz Equation (38) is a particular solution of this system, which has the special property that the Boltzmann operator can be linearly expressed
by the generators. We believe this to be a natural auxiliary hypothesis):
In the same coordinate system, the collisionless Boltzmann equation becomes:
Here, some comments are in order. In the structure of Boltzmann Equation (39), as well as in the form Equation (38) of the modified space translations,
enters an unknown function
. Therefore, the form of
cannot be found only from its commutator with the other generators
, but the constraints form the entire conformal algebra must be used, as well as the requirement that
and
are dynamical symmetries of Equation (
39):
In fact, commuting the unknown generators
with
and
, we can fix the
t- and
r-dependence of the yet undetermined functions that occur in them:
In particular, looking for a representation of the analog of the extended Galilei algebra
, we find that the unknown functions
must satisfy the system:
Because of Equation (46), one must distinguish two cases:
, when can be arbitrary
, when is a constant.
In the second case, taking Equations (44) and (45) together, we obtain an equation for
. It is:
and has in general two independent solution:
. It follows that, for a given arbitrary value of
, we have in general two distinct realizations of the analogue of Galilei transformation; and consequently, also two realizations of the analogue of the Galilei algebra. By construction, these are Lie algebras of symmetries of the collisionless Boltzmann Equation (39) (with
):
Next, we include the generators of special conformal transformation
and
to the extended Galilei algebras Equation (48) just constructed. We must also satisfy the other commutators of the conformal algebra Equation (8). Furthermore, the generators of the representation we are going to construct are dynamical symmetries of the collisionless Boltzmann equation (we use the commutators
and
in order to establish a relation between the constants
and
). We find:
for the eigenvalue and:
The system of Equations (43)–(45) and (50)–(71) must give a solution for the unknown functions . Of course, it is possible that several of the above equations are equivalent. Because of this fact, although the above system might look to be over-determined, we have not yet been able to produce an explicit solution without making an auxiliary assumption. A classification of all solutions of the above system is left as an open problem. We shall now describe some examples of solutions of this large system.
Example 1: Let Φ(u) = 0. This case seems to be quite simple, provided it is compatible with our system. From Equation (43), we obtain:
Using this value of
from Equations (52)–(54) and (56)–(58), we directly obtain:
When we substitute the above results in relation Equation (42), we also find:
One can now verify that the above results for the functions
and
satisfy all equations of the above system. Now, we can finally write the algebra generators:
We return to the original variables via the change
, done through the substitutions
and
. Finally, we have the following representation of a conformal symmetry algebra of the collisionless Boltzmann Equation (33):
Proposition 6. The generator Equation (77) close into the following Lie algebra:for ∈
and for an arbitrary dynamical exponent z. They give a representation of the finite-dimensional conformal algebra, which acts as a dynamical symmetry algebra of the Boltzmann equation in the form: Proof. The commutation relation Equation (78) is directly checked. From the commutators
and:
it is seen that they generate dynamical symmetries. ☐
Example 2: Let k = 0. In this case,
left arbitrary, which leads to
from Equation (43) and:
Then, from Equation (42), we obtain:
However, when substituting in Equations (64)–(67), taking also into account that , we find that . Then, it is easy to check that the condition Equations (68)–(71) are fulfilled. This allows us to formulate:
Proposition 7. Let . Consider the generators:where and satisfy: For any triplet , which gives a solution of the system Equations (83) and (84), the generator Equation (82) close into the following Lie algebra:for ∈
and for an arbitrary constant z. Equation (82) is a representation of the finite-dimensional conformal algebra and acts as the dynamical symmetry algebra of the Vlasov–Boltzmann equation, with a quite general “force” term: Proof. The commutators are satisfied for
and
if condition Equations (83) and (84) are fulfilled. Under the same conditions, the symmetries are proven by the relations:
☐
In particular, if we implement the physical requirement that the “force” term should depend only on the positions
r, that is
, we can compute explicitly the representation of the algebra Equation (82). To do this, one must find a solution of the system:
The solution of the second equation is relatively simple, even for an arbitrary
z:
The solution of the Equation (86) for an arbitrary
z can be expressed in terms of hypergeometric functions, but we shall not give its explicit form here. However, for
, the system Equations (86) and (87) have an elementary solution:
Substituting this into the generator Equation (82) for
gives a finite-dimensional representation of the dynamical conformal symmetry of a collisionless Boltzmann equation of the form:
4. Conclusions
In this work, we have described the results of the first exploration of dynamical symmetries of collisionless Vlasov–Boltzmann transport equations. Our main finding is that these equations admit conformal dynamical symmetries, although it does not seem to be possible to extend this to infinite-dimensional conformal Virasoro symmetries, not even in the case of space dimensions. These conformal symmetries are new representations of the conformal algebra and are inequivalent to the standard representation, which is habitually used in conformal field-theory descriptions of equilibrium critical phenomena. Our first class of new symmetries was found by admitting the momentum p (or equivalently, the velocity ) as an additional independent variable, leading to the representations Equations (26), (30) and (31). The second class of symmetries also allowed for external driving forces , and it has been one of the questions of which types of forces should be compatible with conformal invariance. As an example, we have seen that time-independent forces , with an arbitrary scaling function φ, are admissible and lead to the general representation Equation (82). However, the solutions of the associated system of equations for the coefficients have not yet been classified and the complete content of these representations remains to be worked out in the future.
Some intuition can be gleaned from some examples. We have written down the explicit representations for the force
, with an
arbitrary Equation (77) and for
Equations (82), (86), (87) with an arbitrary
. In the later case, which could be related to physical situations, we have given the explicit representation of the conformal algebra for
, when
, Equations (82) and (89). Having identified these symmetries, the next step would be to use these to find either exact solutions [
16] or else to use the algebra representations for fixing the form of co-variant
n-point correlation functions, in analogy to time-dependent critical phenomena; see, e.g., [
12].
The results derived here can be used as a starting point to derive forms of the transition rates
w in the collision terms, which would be compatible with the dynamical symmetries of the collision-free equations. This kind of approach would be analogous to the one used for finding dynamical symmetries of non-linear Schrödinger equations; see, e.g., [
17,
18]. We hope to return to this elsewhere.