Landau damping in the multiscale Vlasov theory

Vlasov kinetic theory is extended by adopting an extra one particle distribution function as an additional state variable characterizing the micro-turbulence internal structure. The extended Vlasov equation keeps the reversibility, the Hamiltonian structure, and the entropy conservation of the original Vlasov equation. In the setting of the extended Vlasov theory we then argue that the Fokker-Planck type damping in the velocity dependence of the extra distribution function induces the Landau damping. The same type of extension is made also in the setting of fluid mechanics.


1.
Introduction. The physical system under investigation is a gas composed of particles interacting via long range forces. We call this gas a Vlasov gas rather than plasma in order to emphasize that the interaction forces that we consider in this paper are not associated with any type of field (e.g. the electromagnetic field) that has its own dynamics (e.g. Maxwell's electromagnetic field theory). Because of the presence of long range interactions, we expect to see in the Vlasov gas a collective behavior and thus large pair and even higher order correlations. In the position space, the correlations manifest themselves in the formation of fine scale spatial structures, in the velocity space, in the micro-turbulence. It is thus surprising that the Vlasov kinetic theory, with only one particle distribution function serving as the state variable and with the mean-field consideration of long range interactions, has proved to provide a good theoretical description of the behavior observed in experiments. In this paper we put into focus one particular experimental observation that became known as Landau damping. It is the time irreversible evolution of the one particle distribution function toward spatial homogeneity. Before being experimentally observed [21], this behavior was predicted [17] from an analysis of solutions of the Vlasov kinetic equation. Landau damping is particularly surprising since the Vlasov kinetic equation is time reversible. How can the reversible time evolution equation have time irreversible solutions? The recent mathematical analysis in [24], extending Landau's investigation of the linearized Vlasov equation to the complete nonlinear setting, points out to the analogy, both in the physics and the mathematics involved, with the general conundrum of statistical mechanics, namely with the question of how does the reversible Hamiltonian time evolution of an isolated system of 10 23 particles imply the observed time irreversible approach to thermodynamic equilibrium states. Inspired by this analogy, we investigate in this paper the Landau damping by using some of the methods developed originally in statistical mechanics.
Following Boltzmann, the tendency of solutions of the microscopic time evolution equations to increase their irregularity is mathematically expressed in the emergence of dissipative and irreversible time evolution of appropriately chosen mesoscopic state variables. This strategy can be seen as a way to regularize the complex microscopic time evolution. The loss of regularity and the time irreversible weak convergence seen in the microscopic time evolution becomes after regularization an explicitly visible regular property of the time evolution of appropriately chosen mesoscopic state variables in which many microscopic details are ignored. The weak convergence is transformed in the regularization into the strong convergence. In the example of the Boltzmann gas (gas particles interact only via very short range repulsive forces), the particle trajectories become very complex (their irregularity increases) due to strong interactions experienced by the gas particles during collisions. Boltzmann's regularization of the Boltzmann gas dynamics consists in choosing only one particle distribution function as the state variable and taking into account the influence of collisions in the Boltzmann collision term that is dissipative and time irreversible. Solutions to the Boltzmann equation are regular from the mathematical point of view and show approach to the Maxwell distribution in velocities and to the spatial homogeneity, both observed in experiments.
One possible way to regularize the Vlasov equation is to add to its right hand side a self-diffusion term that drives solutions to the spatial homogeneity. For such regularized Vlasov equation the Landau damping becomes a regular and explicitly visible property of its solution. In order to see the physics that is behind the added self-diffusion term, we step out of the confinement of the Vlasov one particle kinetic theory and formulate (in Section 4) the time evolution of the Vlasov gas on a different, more detailed, level. Two questions arise: first, how do we formulate such extended Vlasov theory, and second, how does the Landau damping present itself in solutions of the governing equations of the extended theory.
The paper is organized as follows. In Section 2 we present a general nondissipative one particle kinetic equation and recall three properties of their solutions (reversibility, Hamiltonian structure, and entropy conservation). A special case, corresponding to the Boltzmann gas, is reviewed in Section 3. We introduce the Boltzmann regularization (based on considering collisions as the source of dissipation). In Section 4, we construct extended Vlasov kinetic equation (whose solutions still possess the three properties identified in Section 2 for one particle kinetic equations), regularize them (based on considering the decay of micro-turbulence as the source of dissipation) and show then how this dissipation induces the Landau damping in the one particle distribution function. In Section 5 we make the same type of investigation of the Vlasov gas but with the level of fluid mechanics replacing the level of kinetic theory. 2. Nondissipative one particle kinetic theory. Let f (r, v) be one particle distribution function, r and v denote the position coordinate and momentum of one particle. Its time evolution is governed by where E(f ) is a real valued function of f having the physical meaning of energy, E f denotes derivative of E with respect to f , L (f ) is the Poisson bivector expressing mathematically the kinematics of f . It is given by associating it with the Poisson bracket where A and B are real valued as sufficiently regular functions of f . We recall that {A, B} (f ) given in (2) holds (which can either be verified directly or by associating it with the Lie group of canonical transformations of (r, v) -see [22]). As for the boundary conditions, we require that all the integrals over the boundary arising in by parts integration disappear. Written explicitly, the kinetic equation (1) takes the form From the mathematical point of view, (1) is a general formulation of the Hamiltonian dynamics, i.e. a covector constructed as a gradient of a potential (having the physical meaning of energy) is transformed into a vector (serving as a vector field generating the time evolution) by a Poisson bivector. From the physical point of view, (1) expresses the mechanics. It may not be the classical canonical mechanics (in which the role of f is played by (r, v) and L (rv) = 0 1 −1 0 ) but still it is the time evolution possessing the physically essential features of mechanics. The noncanonical formulation (1) of mechanics has been introduced by Arnold [2] in the context of fluid mechanics (putting the Euler equation, that has already been proven to represent the Hamiltonian dynamics by Clebsch in [4], into the form (1)). The Poisson bracket (2) has been introduced in [22]. Now we turn to important qualitative properties of solutions to (1). It follows directly from the general formulation (1) that solutions to the kinetic equation (1) possess the following three properties: Conservation of the energy The skewsymmetry of L (f ) implies immediately the conservation of the energy The energy conservation is indeed a cornerstone of the physics involved in mechanics.

Conservation of entropy
Another conservation appears as a consequence of the degeneracy of L (f ) . The degeneracy is defined as follows. If there exists a potential C (f ) (f ) (in addition to the potential It can be directly verified that potentials of the form C (f ) (f ) = dr dvc(f (r, v)), where c is a sufficiently regular function c : R → R, are all Casimirs of the Poisson bracket (2). Equation (5) represents thus an infinite number of conservation laws that supplement the energy conservation law (4). We note that in the completely microscopic classical mechanics of particles these additional conservation laws are absent since the Poisson bivector arising there is nondegenerate. The essential difference between the microscopic particle mechanics and mesoscopic mechanics (as e.g. the kinetic theory represented by (1)) is that the former keeps all the details while the latter ignores some of them. The ignored details have not however disappeared from the time evolution of the mesoscopic state variables ( i.e. the state variables that we have decided to keep; in kinetic theory we keep the one particle distribution function). How shall we express their influence? It is thermodynamics that answers this question. The influence of the ignored details on the time evolution of the mesoscopic state variables is expressed in the time evolution of a new potential, called entropy (we shall use the symbol S to denote it). The entropy S tends to its maximum allowed by constraints. In the case of the nondissipative mesoscopic dynamics that we discuss in this section, the entropy S is a potential that remains unchanged. In the nondissipative mesoscopic dynamics the ignored details manifest themselves in the emergence of a companion conservation law. Since the Casimirs C are conserved (see (5)), we choose the entropy to be one of the Casimirs C(f ). To continue and select one among the infinite number of Casimirs requires an investigation of the dissipative dynamics (in Sections 3 and 4).

Time reversibility
So far, we have made no restriction on the choice of the energy E (f ) (f ). Now, we make one. We assume that the energy is invariant with respect to the transformation v → −v. If this is the case then we see immediately that (1) is time reversible in the sense that the transformation t → −t is fully compensated by the transformation v → −v (i.e. the trajectory f (r, v, t) corresponding to the initial condition f 0 (r, v) is the same as the trajectory f (r, −v, −t) corresponding to the initial condition f 0 (r, −v)).
3. Boltzmann kinetic theory. In this section we recall several well known results about the Boltzmann kinetic equation. Our objective is to prepare the setting for discussing the Vlasov kinetic equation in Section 4 and also to highlight similarities and differences between the Boltzmann and the Vlasov kinetic equations.
3.1. Nondissipative Boltzmann equation. The physical system under investigation in this section is the Boltzmann gas (we use the superscript (f B) to denote the quantities that arise in this investigation). The energy E (f B) (f ) is given by where is the kinetic energy, m denotes the mass of one gas particle, and The potential V (f B) is a very short-range hard-core type two-particle repulsive potential, f 2 (f ; r 1 , v 1 , r 2 , v 2 ) is the two particle distribution function expressed in terms of the one particle distribution function f . The Boltzmann gas particles are free, they are not subjected to any force, except when two particles become very close. In such case, the two particles are subjected to a large repulsive force. The kinetic equation governing the time evolution of the Boltzmann gas is thus (1) with the energy (6), i.e.
What can we say about solutions to (9)? First, we note that the four properties proven in Section 2 hold provided f 2 (f ; r 1 , v 1 , r 2 , v 2 ) is sufficiently regular and invariant with respect to (v 1 , v 2 ) → (−v 1 , −v 2 ) function and C is a Casimir of the bracket (2). The potential N (B) (f ) = dr dvf (r, v) is a Casimir that has the physical interpretation of the number of moles. Regarding more detailed properties of solutions to (9), we recall that experimental observations of the Boltzmann gas show that, in the absence of external forces, the Boltzmann gas approaches (as t → ∞) a distribution function f eq (r, v) that is independent of r and Maxwellian (i.e. proportional to e − v 2 2mk B T , where k B is the Boltzmann constant and T is the equilibrium temperature) in the dependence on v. Is this experimental observation seen in solutions to (9)?
To begin with answering this question, we first recall (see [24]) that in the absence of the potential energy (8) in the energy (6), solutions to (9) converge weakly in large time to a spatially homogeneous distribution that is equal to the space average of the initial distribution function f 0 (r, v). We expect that this property will still hold for the full kinetic equation (9) but, to the best of our knowledge, such property of solutions to (9), as well as the Maxwellian dependence on v in the asymptotic solutions to (9), has not been rigorously proven. There are many results coming from direct simulations of systems composed of (up to 10 5 ) particles interacting via potential of the type (8) but such results, while addressing the time evolution of the Boltzmann gas, are not results about solutions to the kinetic equation (9).
Nevertheless, we may assume that solutions to (9) show time irreversible tendency towards spatial homogeneity and Maxwell distribution in v. This tendency is not seen in the properties (10).
In the further investigation of solutions to (9) we shall follow Boltzmann and make the time irreversible tendency of solutions to (9) manifestly visible by modifying Eq.(9) into a new equation whose solutions have the properties (10) except that the last line in (10) is absent and in the second line the Casimir C (f ) (f ) is a one specific Casimir called the Boltzmann entropy and the equality is changed into an inequality. It is this inequality and the time irreversible nature of the equation that brings into the visibility the time irreversible tendency of solutions to (9). The modified (regularized) kinetic equation (9) is then called a Boltzmann equation. The transformation of (9) into the Boltzmann equation is in fact a first step in the investigation of solutions to (9). The Boltzmann equation is simpler than (9), it is easier to solve it, and its solutions are expected (the expectation is based on the physical insight leading from (9 to the Boltzmann equation) to approximate well solutions to (9). This is the Boltzmann strategy (recalled below in the context of the Boltzmann gas) that we shall also follow in the context of the Vlasov gas in Section 4.

3.2.
Regularization. The physical insight on which the Boltzmann regularization of (9) is based is that the very rapid and very large changes in momenta of gas particles during collisions are responsible for the spatial homogenization and also for the Maxwell distribution in velocities seen in the one particle distribution function. Details of the trajectories during collisions are ignored in the Boltzmann analysis.
The first step in the Boltzmann modification is a complete omission of the hardcore repulsive potential energy (8). Such omission is justified if the kinematics is modified in such a way that the particles are prevented to be close enough to make the potential energy (8) different from zero. This is easy to do if we regard the particles as point particles so that, even without changing kinematics, the potential energy (8) completely disappears. The total energy of the Boltzmann gas is thus only the kinetic energy of the gas particles. But even with this idealized particles, we still have to admit that collisions occur and that they influence the behavior of the gas. Having no mechanical basis to express their influence (recall that the second term on the right hand side of (1) is now absent) we turn to chemical kinetics. Collisions are seen in the Boltzmann analysis as chemical reactions. The momentum v serves as the variable parameterizing the species. The chemical reaction representing the binary collision is a reaction in which two species, one with the parameter v and the other with the parameter v 1 , interact. The outcome of the reaction are two new species, one with the parameters v and v 1 . The mechanical nature of this "chemical reaction" is expressed in two constraints, having the physical meaning of the energy the momentum conservations in binary collisions. With this viewpoint of collisions, the second term on the right hand side of (9) becomes replaced with the so called Boltzmann collision term ∂f ∂t diss .
With such replacement, the kinetic equation (1) turns into the Boltzmann kinetic For the mathematical expression of ∂f ∂t diss we turn now to the mass action law in chemical kinetics (see [11]): where Ξ (B) (f, f * ). called a dissipation potential, is given by the thermodynamic force X (B) (called affinity in the context of chemical kinetics) by , playing the role of the rate coefficients, equals zero except when (11) holds and then it is positive, and S (B) (f ) is a Casimir (see (5)). The mechanical origin of binary collisions is retained in (13) only in the constraints (11) expressing mechanical conservation laws, and, indirectly, in the choice of S (B) (f ) that we shall discuss below. It can directly be verified that the collision term is invariant with respect to the transformation v → −v and thus it is the term that breaks the time reversibility of the equation (9).

3.3.
Properties of solutions. Now we turn to solutions to the Boltzmann equation (12). We begin by noting that These three properties are directly seen in (14). From (16) we conclude that solutions to (12) approach, as t → ∞, solutions to X (B) = 0. We shall denote the distribution functions that arise as solutions to X (B) = 0 by the symbol f leq (r, v), the submanifold that they form by the symbol M coll where M (kt) denotes the state space of kinetic theory.
So far, the Casimir S (B) (f ) is unspecified. Now we specify it. There are in fact three routes that we can follow to specify it. All three routes lead to the same Boltzmann entropy (that is one particular On the first route we ask the question of what is X (B) . Since -see (15) -X (B) is directly related to S (B) (f ), we ask in fact the question of what is S (B) (f )) for which solutions to X (B) = 0 have the experimentally observed Maxwellian dependence on v. It is easy to see that the answer is: On the second route we look for the entropy for which the fundamental thermodynamic relation implied by (12) is the fundamental thermodynamic relation representing in equilibrium thermodynamics the ideal gas. This is because the physical system whose time evolution is represented by the Boltzmann kinetic equation (12) is the ideal gas. We briefly recall (see more in [12]) the way how to derive the fundamental thermodynamic relation implied by a kinetic equation whose solutions have the properties (16). First, we construct the thermodynamic potential where T is the equilibrium thermodynamic temperature and µ the equilibrium thermodynamic chemical potential. Next, we solve the equation Φ The solution, denoted f eq (r, v, T, µ) and called an equilibrium distribution, is approached as t → ∞. Finally, the fundamental thermodynamic relation P = P (T, µ) implied by (12) (P is the equilibrium thermodynamic pressure) is where V is the volume of the region in which the gas under investigation is confined. It is easy to verify that in order that (19) is the fundamental thermodynamic relation of the ideal gas, S (B) (f ) has to be the Boltzmann entropy. The third route is the route taken by Ludwig Boltzmann in his original introduction of the Boltzmann equation. The collision term ∂f ∂t diss emerges first from an analysis of the mechanics of binary collisions. Boltzmann then investigates solutions of his equation by looking for a potential whose time evolution obeys the inequality in the third equation in (16). Such search amounts in fact to casting the collision term, derived initially from mechanics of binary collisions, into the mass-action-law form (13).
From the above discussion we see that the properties (16) suffice to show agreement with the experimental observation of the approach to the Maxwellian distribution in velocities but they do not, by themselves, suffice to prove the observed approach to the spatially homogeneous distribution. To prove it, we can ether follow Refs. [8], [5] and investigate further solutions of (12) or follow the spirit of the Boltzmann passage from (9) to (12) and make another passage from (12) to a modified Boltzmann equation in which the approach to the spatially homogeneous distribution becomes manifestly displayed.
Grad [8], Desvillettes and Villani [5] investigate the influence of the first term on the right hand side of (12) (i.e. the free flow term that leaves the entropy S (B) unchanged) on solutions of the Boltzmann equation. They have shown that solutions to Boltzmann equation, (12), approach rapidly a neighborhood of the locally Maxwellian manifold M coll (see (17)). Due to the influence of the free flow term, however, the solutions touch it only at the end of the time evolution, where f does not change in time and is spatially homogeneous, i.e. the distribution f reaches the total Maxwellian distribution, which is the minimum of the thermodynamic potential Φ (B) (see (18)).
The second way to investigate the tendency to spatial homogenization in solutions to the Boltzmann equation (12) is motivated by another experimental observation made on the Boltzmann gas. According to this observation the time evolution of the Boltzmann gas becomes, as t → ∞, well described by fluid mechanics that is a mesoscopic theory in which more microscopic details are ignored than they are ignored in the kinetic theory. The proof consists of showing first that solutions to (12) become, as t → ∞, well approximated by solutions to the governing equations of fluid mechanics and then, in the second part, showing that solutions to the governing equations of fluid mechanics approach, again as t → ∞, spatially homogeneous solutions. Regarding the first part of the proof, there are several methods (the most well known among them is the Chapman-Enskog method) that can be used to reduce kinetic theory to fluid mechanics. All the methods display well the underlying physics but their mathematical rigor is inferior to the mathematical rigor with which the results in [5] are derived. Another way to transform the Boltzmann equation (12) into an equation manifestly displaying approach to spatially homogeneous solution can be an adaptation to the Boltzmann equation of the modification introduced in the next section for the Vlasov equation. We shall return to this possibility later in this paper.
In order to find still more information about solutions to the Boltzmann equation, we can restrict ourselves to particular stages in the time evolution and correspondingly simplify the Boltzmann equation. For example, in the final stages of the approach to the thermodynamic equilibrium states f eq (r, v) the thermodynamic force X (B) (see (15)) is small and we can thus approximate the dissipation potential Ξ (B) (see (14)) by Ξ (B) (f, f * ) = dr dv dv W (B) (X (B) ) 2 . If moreover the changes of the momenta in the binary collisions are small, we can approximate X (B) by X (B) = ∂f * ∂v . With these approximations, ∂f ∂t diss takes the Fokker-Planck form where Λ (B) is an operator that is degenerate (in order that the energy is preserved during the time evolution) and positive definite if acting outside of its null-space, and f * = Φ f . An additional physical insight into the Boltzmann regularization can be gained by regarding the kinetic equations (9) and (12) as Liouville equations corresponding to equations governing the time evolution of one quasi-particle and then interpreting the regularization as a modification of the quasi-particle dynamics. We now proceed to identify the quasi-particle dynamics.
We begin by noting that entropy conservation in the nondissipative time evolution allows us to write the kinetic equation (9) in the form and the Boltzmann equation (12) Note that this dissipation is very similar to the Landau collision integral [18], §41.
In the case of the general Boltzmann collision term (13) we still can write the Boltzmann kinetic equation in the form (22) (see [9]) if we use the identity where ϕ is a sufficiently regular mapping ϕ : R N → R N , y(η) = y − ηy ; y (η) = y + (1 − η)y used in particular in [16] Next, we regard the above two kinetic equations as Liouville equations corresponding toṙ in the case of (21) andṙ in the case of the Boltzmann equation (22). Both the quasi-particle-dynamics equations (24) and (25) involve the distribution function f and have to be therefore considered always together with their corresponding Liouville equations (21), in the case of (24), and (22), in the case of (25). The regularization of the kinetic equation (9) into the Boltzmann equation (12) we can now see as replacing the force T ∂f * ∂r in the quasi-particle dynamics (24) with the friction force Λ (B) ∂f * ∂v . Note that this friction conserves the energy (see e.g. [18], §41), it is a friction in the sense that it causes an increase of the entropy. 4. Vlasov kinetic theory. Now we turn to the main subject of this paper. We consider the Vlasov gas and the Vlasov kinetic equation governing its time evolution. The Vlasov equation has still the form (1) with the Poisson bivector given by (2) and the energy where V (f V ) pot (|r 2 − r 1 |) is the two-particle long range attractive potential. Written explicitly, the Vlasov kinetic equation has the form We see that the equations governing the time evolution of the Boltzmann gas (Eq.(9)) and the equation governing the time evolution of the Vlasov gas (Eq. (27)) are the same except that the potential energy in the case of the Boltzmann gas is short range and repulsive and in the case of the Vlasov gas long range and attractive. As for the solutions to (27), we note, as we did in the case of (9), that (27) is a particular realization of the general kinetic equation (1) and thus the four properties of solutions derived in Section 2 where C (V ) is a Casimir of the bracket (2), and N (V ) (f ) = dr dvf (r, v) is another Casimir, that has the physical interpretation of the number of moles, hold also for (27). Regarding more specific properties, we shall concentrate on the experimentally observed [21] tendency to spatial homogenization known as Landau damping. It has been rigorously proven in [24] that solutions to (27), both with and without the second term on its right hand side, converge weakly, as t → ∞, to spatially homogeneous distribution. Inspired by the Boltzmann regularization of (9), our objective in this paper is to modify (27) into a new kinetic equation in which the tendency to the spatial homogenization is manifestly displayed.

4.1.
Extension. As in the Boltzmann regularization, our first task is to attempt to recognize the physics that is behind the property of solution of the governing equation that is in the focus of our interest. In the case of the Vlasov equation it is the Landau damping (i.e. the time irreversible tendency towards the spatial homogenization). There does not seem to be a general agreement about the physics behind the Landau damping (see e.g. Section 1.2 in [24]) but there does seem to be an agreement that, due to the phase mixing, the particles experience a complex motion that can be regarded as micro-turbulence (i.e. the turbulence seen on the level of kinetic theory with the one particle distribution function serving as the state variable) [28,30].
In the rest of this paper we extend the Vlasov kinetic theory. By an extension we mean first of all an enlargement of the state space. How do we choose the extra state variables? In principle, we choose them as the state variables through which we can best and in the most succinct way express the new features that are missed in the original non extended theory. In practice, the extra state variables are often found in the vector field of the original dynamical theory (as it is done in the Grad type extensions of the fluid mechanics in [25], [15], [29] by adopting hydrodynamic fluxes as the extra state variables) or in the investigation of the more detailed (more microscopic) nature of the macroscopic system under investigation (as it is done in fluid mechanics of complex fluids for example in [16], [3]). The former type of extension is, from the physical point of view, an extension by introducing an extra inertia, the latter type of the extension is, from the physical point of view, an involvement of the internal structure whose time evolution cannot be decoupled from the time evolution of the interest in the original dynamical theory.
Having chosen the extra state variables, we face the problem of establishing equations governing their time evolution. In the Grad type extensions this is typically done by suggesting closures of the infinite Grad hierarchy. In the internal-structure type extensions the time evolution equations are typically found, first, by investigating the time evolution of the internal structure (e.g. the time evolution of a dumbbell suspended in a fluid -see e.g. [3]), and second, by investigation the coupling between the time evolution of the original and the time evolution of the extra state variables (expressed, for example in the context of fluid mechanics of complex fluids in the dependence of the stress tensor on the extra state variables). A complementary and a more systematic way of finding the time evolution equations of an extended theory is to require that both the original and the extended dynamical theories share the mathematical structure guaranteeing agreement of their solutions with certain basic experimental observations.
In our extension of the Vlasov kinetic theory we shall follow the internal-structure type extension and we shall require that the time evolution in both the original and the extended theory is Hamiltonian. The extended kinetic theory that we shall formulate below involves an extra state variable that is chosen to be the one particle distribution function characterizing the internal structure induced by the micro-turbulence. In the nutshell, we replace the distribution function f (r, v) with a pair of distribution functions We interpret φ(r, v) as the average (or regularized) distribution function and ψ(r, v) as a distribution function characterizing the internal structure. The energy (26) acquires in the extended theory a new term expressing the extra kinetic energy involved in the small scale micro-turbulent motion. The extended Vlasov kinetic equation (i.e. the equation governing the time evolution of (29) is Hamiltonian and its solutions satisfy the properties (28) (appropriately adapted to the extended theory). The regularization (discussed in Section 4.2) of the extended Vlasov kinetic equation is then made by adding to it a Fokker-Planck like dissipation of the microturbulent structure. In the course of the time evolution, this dissipation induces then the Landau damping of φ(r, v).
We now proceed to introduce the equation governing the time evolution of (29), i.e. the equation extending the Vlasov equation (27). The physical requirement on which we base our construction is the preservation of the Hamiltonian structure and of the properties (28) of solutions of the original Vlasov equation (27). We therefore begin by extending the Poisson bracket (2) to the Poisson bracket expressing kinematics of (29). In order to see well both the physics and the mathematics that is behind the extension, we shall make the extension by following two routes. On the first route (Section 4.1.1), the extension is motivated by the geometrical content of the bracket (2), on the second route (Section 4.1.2), by the Reynolds approach to the turbulence in fluid mechanics. 4.1.1. Geometrical approach. From the physical point of view, we assume that the kinematics is still determined by the kinematics of φ(r, v) expressed in the Poisson bracket (2) (in which f is replaced by φ) and that the newly adopted distribution function ψ(r, v) is passively advected. From the mathematical point of view, we shall make therefore the extension with the help of the concept of semidirect product (see e.g. [23], [6]).
Poisson bracket (2) can be regarded as the Lie-Poisson bracket on the dual of the Lie algebra of Hamiltonian vector fields on the cotangent bundle T * Q, Q ⊂ R 3 of one-particle classical mechanics, see [23]. The vector fields where h is the one particle Hamiltonian, generate the Lie group G of canonical transformations on T * Q. Hence, a natural extension of the dynamics generated by (2) is a geometric extension (semidirect product) of the underlying Lie group. The vector field (30), which depends on the choice of the Hamiltonian h, generates motion of a particle within the cotangent bundle and that motion can be regarded as a Lie group. The vector field form an algebra called Lie algebra of G and denoted by G. Bracket (2) can be seen as the Lie-Poisson bracket governing evolution in the dual of G, where fields of differential forms live. The distribution function φ can be then introduced as the phase-space divergence of the differential forms, and bracket (2) then generates evolution equations for φ.
The Lie group G can be extended by adding an another group passively advected by G. That is the semidirect product, and according to Eq. 3.1 in [22] the Poisson bracket constructed from the semidirect product G V , which gives evolution on the Lie algebra dual of the semidirect product, is The new density ψ(r, v) is the density expressing the added (passively advected) group. Let us refer to this bracket as to the extended kinetic theory Poisson bracket. The bracket [•, •] stands for the Poisson bracket of classical mechanics, which governs evolution on T Q , Bracket (31) can be then rewritten more explicitly as The evolution equations of φ and ψ can be obtained by casting {A, T Φ} (ext) into the form which leads to where φ * = Φ φ and ψ * = Φ ψ . These are the governing equations of the extended Vlasov theory. They are Hamiltonian as the original Vlasov equations and their solutions have the properties where C (ext) is a Casimir of the Poisson bracket (33) and N (ext) is the Casimir having the physical interpretation of the number of moles, that are the same as (28). Regarding the energy E (ext) (φ, ψ), the physical interpretation of (29) leads us to The new contribution E (turb) (φ, ψ) to the total energy E (ext) (φ, ψ) is the energy associated with the micro-turbulence. It is therefore mainly an extra kinetic energy. The total kinetic energy is thus a sum of the overall kinetic energy dr dvφ(r, v) v 2 2m , that depends only on the overall distribution function φ, and the extra kinetic energy E (turb) (φ, ψ) of the fine scale micro-turbulent motion. For the qualitative analysis of solutions to (35), presented in Section 4.2, we do not need a specific form of E (turb) (φ, ψ). We assume that the contribution of the micro-turbulence to the potential energy of the Vlasov gas is negligible. The potential energy given by the second term in (37) is thus assumed to represent the complete potential energy also in the extended Vlasov theory.

4.1.2.
Reynolds approach adapted to kinetic theory. Our objective is to adapt the classical Reynolds approach to the turbulence (developed originally by Reynolds [27] inside fluid mechanics) to the micro-turbulence of the Vlasov gas (formulated now inside kinetic theory). We anticipate that in the Reynolds-type reformulation of the Vlasov equation we will be able to see the Landau damping similarly as the extra Reynolds stress (arising due to the turbulence) is seen in the Reynolds analysis.
We begin by recalling the essential steps that are made in the original Reynolds analysis. The starting point is the Navier-Stokes equation governing the time evolution of the velocity field that plays the role of the state variable. In the first step, this equation is extended to a more microscopic level of description (i.e. to a level on which more details are taken into account). This is done by promoting the velocity field to the status of the random variable and extending the Navier-Stokes equation to the corresponding to it Langevin equation. Alternatively, the extension of the velocity field can be done by representing the random velocity field as the average field together with an infinite hierarchy of moments that are physically interpreted as fields characterizing the internal structure. The Langevin equation governing the time evolution of the random velocity field then provides the equations governing the time evolution of this extended set of state variables. The subsequent passage from the infinite to a finite number of moments requires a closure of the infinite hierarchy of the time evolution equations.
We take the above description of the Reynolds analysis as a motivation rather than a recipe. The starting point in our Reynolds-type investigation of the microturbulence is the Vlasov equation (27). Our first step is not the promotion of the one particle distribution function to the status of random variable but rather its replacement by a pair of distribution functions (29). Both distribution functions φ(r, v) and ψ(r, v) are regarded as being independent. The equations governing their time evolution are then constructed by requiring that the Hamiltonian structure of the original Vlasov equation (27) is preserved in the extended time evolution equations.
If we put the Reynolds analysis into the setting of the Hamiltonian dynamics and recall that the state variable in the Hamiltonian formulation of the Euler equation is the momentum field (see e.g. [22], [14]) and that the velocity field is its conjugate, we see that the essential step in the Reynolds analysis, namely the separation into an average and fluctuations, is not made for the state variable but for its conjugate. This means that in the setting of kinetic theory we shall replace the distribution function f (r, v) with a pair of distribution functions (29) in such a way that In order to arrive at the Poisson bracket expressing the kinematics of (29), we proceed as follows. First, we introduce two uncoupled and identical Vlasov equations, one for the distribution function f and the other for the distribution function g. The Poisson bracket expressing the kinematics of (f, g) is simply where {A, B} (f ) is the bracket (2) and {A, B} (g) is the same bracket but with f replaced by g. Now, we look for a one-to-one transformation (f, g) (φ, ψ) for which f * = φ * + ψ * (which is the relation (38)) and g * = φ * . We see easily that such transformation is If we now apply this transformation to (39) we obtain where {A, B} (φ) is the Poisson bracket (2) with f replaced by φ, {A, B} (ψ) is the Poisson bracket (2) with f replaced by ψ, and Since the transformation (40) is one-to-one and (39) is a Poisson bracket then also the bracket (33) is a Poisson bracket (in particular, we are certain that the Jacobi identity holds for (33)). If we now compare the bracket (41)  involves two extra terms: ∂ ∂r T ψ ∂ψ * ∂v + ∂ ∂v T ψ ∂ψ * ∂r .

4.2.
Regularization. In order to get more detailed properties of solutions to the extended Vlasov equation (35) (more detailed than (36)), we can either follow Desvillettes and Villani [24] and adopt their rigorous mathematical analysis of the Vlasov equation (27) to the extended Vlasov equation (35) or we can follow Reynolds and, first, regularize (35), and then investigate solutions of the regularized equation. The regularization of (44) consists of supplying (35) with an explicit dissipation term expressing the physics of the micro-turbulence. We shall follow the latter route. We assume that the thermodynamic force X (V ) (compare with (15)) that drives the decay of the micro-turbulence is effective velocity ∂ψ * ∂v . We limit ourselves to the quadratic dissipation potential where Λ (V ) is a symmetric, degenerate (in order the total energy remains conserved, i.e. Λ (V ) ∂E (ext) ψ ∂v = 0) and positive definite if acting outside its null-space. The regularized extended Vlasov kinetic equation becomes The dissipative term in the last line is −Ξ Solutions to (44) have the following properties: To continue the investigation of solutions to (44), we regard it as a pair of Liouville equations corresponding to two quasi-particle time evolution equations (compare with (25)). The equation governing the time evolution of the "average" quasi-particle areṙ and the equations governing the time evolution of the "turbulence" quasi-particle areṙ Now we use our physical insight into the micro-turbulence and estimate the relative importance of the terms appearing on the right hand side of (46) and (47).
From the assumptions that we made about the energy E (ext) (see (37)), we can assume that the extra force T ∂ψ * ∂r brought about by the micro-turbulence is small relative to the force T ∂φ * ∂r brought about by the interaction potential V φV pot . On the contrary, the extra velocity of the "average" quasi-particle that is proportional to ∂ψ * ∂v and that arises due to the micro-turbulence, is not, in general, small relative to the velocity T ∂φ * ∂v . Moreover, we can assume that the "turbulent" quasi-particle has a small mass and thus we we can limit ourselves to the motion with zero acceleration. This means (if we omit the force T ∂φ * ∂r ) that If we now insert this equation to the first equation in (46), we obtaiṅ where is the diffusion coefficient of the Landau damping. With (49) and with omitting the extra force T ∂φ * ∂r in the second equation in (46), the equation governing the time evolution of the "average" distribution function φ(r, v) becomes This is the Vlasov equation (27) equipped with the diffusion term (the last term on the right hand side of (51)) that manifestly displays the Landau damping.
5. Spatial homogenization in the extended fluid mechanics. As we noted in Introduction, it is already surprising that the one particle kinetic theory has been found to be microscopic enough to provide an appropriate setting for describing dynamics of the Vlasov gas. We therefore do not expect that such description can be made in the even less microscopic setting of the classical fluid mechanics. Nevertheless, we can still ask the question as to whether there exist extensions of the classical fluid mechanics which predict spatial homogenization (in the sense that we have seen above in the setting of kinetic theory) and could thus become candidates for the setting suitable for the Vlasov gas. In this section we define what we mean by Landau damping in fluid mechanics (in the last paragraph before Section 5.1) and introduce two extensions predicting it. We begin by recalling the Hamiltonian formulation of the classical fluid mechanics. The fields playing the role of the state variables are (ρ(r), u(r), s(r)) (52) denoting the mass, momentum, and entropy density respectively. Their kinematics is expressed mathematically in the Poisson bracket (see e.g. [22]) We use hereafter the summation convention; i = 1, 2, 3; α = 1, 2, 3. The time evolution equationsȦ = {A, E} (hyd) ; ∀A, where E(ρ, u, s) is the energy, are the familiar equations of the Euler fluid mechanics with where σ = uE u (56) and p = −e + ρE ρ + sE s + < u, E u > (57) The last equation in (54) is a consequence of the three equations above it (it is a companion local conservation law) but we have added it for the completeness.
Next, we investigate solutions to (54). First, we note that the properties (10) (also (28), or (36)) hold also for solutions to (54). The thermodynamic potential Φ (N SF ) is the same as (18)  The energy E (N SF ) is assumed to be invariant with respect to the transformation (ρ, u, s) → (ρ, −u, s).
In order to investigate solutions to (54) in more details, we can again take two routes. Either we follow example of [24] and take the route of a direct and rigorous investigation or we take an indirect route inspired by the methods of statistical mechanics. As in the previous sections, we choose the latter. Assuming that the forces driving fluids to equilibrium, denoted X (N SF ) , are proportional to the gradient of velocity (the Navier-Stokes force) and the gradient of temperature (the Fourier force), introducing dissipation potential Ξ (N SF ) that is a quadratic function of X (N SF ) , and supplying the Euler hydrodynamic equations (54) with the corresponding to it dissipative terms, we arrive at the governing equations of the Navier-Stokes-Fourier fluid mechanics. It can be shown (e.g. by following the Chapman-Enskog method) that solutions to the Navier-Stokes-Fourier equations approximate well asymptotic solutions to the Boltzmann equation (12). We can thus regard the Navier-Stokes-Fourier dissipation in the setting of fluid mechanics as an analog of the Boltzmann collision dissipation in the setting of kinetic theory and anticipate that the approach to the spatially homogeneous equilibrium in the Boltzmann kinetic theory proceeds in the same way also in the Navier-Stokes-Fourier fluid mechanics.
More specifically, we introduce two manifolds and M inc = {f ∈ M (hyd) |ρ = const., div v = 0} (59) where M (hyd) denotes the state space of the classical hydrodynamics. The manifold M N SF is composed of hydrodynamic fields for which the conjugates of u(r) and e(r) are spatially homogeneous but the field ρ(r) can still change with r. The manifold M N SF plays in the Navier-Stokes-Fourier hydrodynamics an analogical role as the manifold M coll (see (17)) plays in the Boltzmann kinetic theory. The manifold M inc is the manifold on which the fluid has the spatially homogeneous mass density and div v = 0. The thermodynamic equilibrium state is the state for which the thermodynamic potential Φ (N SF ) reaches its minimum, it is an element of both M N SF and M inc . We conjecture that the approach to the thermodynamic equilibrium in the Navier-Stokes-Fourier hydrodynamics proceeds in the same way as in the Boltzmann kinetic theory (see the third paragraph after (18)). We thus conjecture that solutions to the Navier-Stokes-Fourier hydrodynamic equations (54) approach rapidly the vicinity of M N SF and M inc but never touch them (due to the influence of the Euler term -the right hand side of (54) -that is analogical to the free flow term in (12)) except at the end of the time evolution where all hydrodynamic fields are independent of t and r.
Results reported in [5], [20], [19], [7], [1] and in [13] are related to the conjecture. The relation to the properties of solution to the Boltzmann equation found in [5] is based on the relation between solutions to (54) and asymptotic solutions to (12) that can be established by following, for example, the well known Chapman-Enskog type investigation of solutions to (12). The authors of Refs. [20], [19] have proven in the case of isentropic flow and the authors of Ref. [7] in the case of non-isentropic flows that in the low Mach number limit solutions to the Navier-Stokes-Fourier equations approach solutions to the incompressible Navier-Stokes-Fourier equations (i.e. the Navier-Stokes-Fourier equations constrained to the manifold M inc ). Thermodynamic interpretation of these result has been suggested in [1]. The authors of Ref. [13] have shown that, in the high Mach number limit, solutions to the Navier-Stokes-Fourier equations approach solutions to a diffusion equation in the mass density ρ(r) which then leads to ρ = const.
Now we turn to the Vlasov gas. We do not expect the setting of the classical fluid mechanics recalled above to be suitable for describing its dynamics. Appropriate extensions of (52)-(57) can however become suitable. What shall we consider to be the requirement for the suitability? We suggest it to be the appearance of the Landau damping in solutions to the governing equations. In the setting of kinetic theory, as we have seen in the previous sections, the Landau damping is the tendency to the spatial homogenization in the absence of the Boltzmann collision term. We therefore define the Landau damping in fluid mechanics as a tendency to the spatial homogenization in the absence of the Navier-Stokes-Fourier dissipation. Our strategy in constructing the extensions will be the same as in kinetic theory. We select extra state variables characterizing the micro-turbulence internal structure, construct Hamiltonian time evolution equations governing their time evolution, introduce an appropriate dissipation of the extra state variables, and finally show that the dissipation induces the self-diffusion.

5.1.
Weakly nonlocal extension. The most straightforward way (we can say a "minimalist" way) to take into account long range interactions in the setting of fluid mechanics is to use energy that is a nonlocal function of (52). With such energy, the role of the extra state variables is played by higher order gradients of the classical state variables (52). Nonlocal and Hamiltonian extensions of the classical fluid mechanics have already been considered in [10] (in the context of fluids that are in the vicinity of gas-liquid phase transition) and in [31] (in the context of fluids involving self-diffusion). We shall therefore limit ourselves only to pointing out the similarity with Section 4.
If we use in the equationȦ = {A, E} (hyd) ; ∀A, (see the text that follows (53)) the energy E that depends on spatial derivatives of ρ(r) and u(r), we still obtain the time evolution equations (54) but with modified fluxes. In particular, the mass