## A geometric approximation to the Euler equations : the Vlasov-Monge-Ampere system

### November 27, 2006

Abstract
This paper studies the Vlasov-Monge-Ampere system ( $VMA$  ), a fully non-linear version of the Vlasov-Poisson system ( $VP$  ) where the (real) Monge-Ampere equation $det\frac{{\partial }^{2}\Psi }{\partial {x}_{i}\partial {x}_{j}}=\rho$  substitutes for the usual Poisson equation. This system can be derived as a geometric approximation of the Euler equations of incompressible fluid mechanics in the spirit of Arnold and Ebin. Global existence of weak solutions and local existence of smooth solutions are obtained. Links between the $VMA$  system, the $VP$  system and the Euler equations are established through rigorous asymptotic analysis.

1 Introduction

The classical Vlasov-Poisson ( $VP$  ) system describes the evolution of an electronic cloud in a neutralizing uniform background through the following equations
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+\xi \cdot {\nabla }_{x}f+{\nabla }_{x}\phi \cdot {\nabla }_{\xi }f=0\end{array}$ (1)
 $\begin{array}{ccc}& & {\epsilon }^{2}\Delta \phi =\rho -1,\end{array}$ (2)
where $f\left(t,x,\xi \right)\ge 0$  denotes the electronic density at time $t\ge 0$  , point $x\in {\mathbb{R}}^{d}$  , velocity $\xi \in {\mathbb{R}}^{d}$  (usually $d=3$  ), $\rho \left(t,x\right)\ge 0$  denotes the ’macroscopic’ density
 $\begin{array}{ccc}\rho \left(t,x\right)={\int }_{{\mathbb{R}}^{d}}f\left(t,x,\xi \right)d\xi ,& & \end{array}$ (3)
and $\phi \left(t,x\right)$  denotes the electric potential at time $t$  and point $x$  generated, through the Poisson equation ( 2 ), where $\epsilon$  is a coupling constant, by the difference between the electronic density $\rho \left(t,x\right)$  and the neutralizing background density, which is supposed to be uniform and normalized to unity. Standard notations $\nabla =\left({\partial }_{1},...,{\partial }_{d}\right)$  and $\Delta ={\partial }_{1}^{2}+...+{\partial }_{d}^{2}$  have been used and $\cdot$  stands for the inner product in ${\mathbb{R}}^{d}$  . The mathematical theory of the $VP$  system is now well understood. In particular, existence of global smooth solutions in three space dimensions has been proved in [? (see also [?, [?). In the present paper, a fully nonlinear version of the $VP$  system is addressed :
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+\xi \cdot {\nabla }_{x}f+{\nabla }_{x}\phi \cdot {\nabla }_{\xi }f=0\end{array}$ (4)
 $\begin{array}{ccc}& & det\left(\mathbb{I}+{\epsilon }^{2}{D}^{2}\phi \right)=\rho ,\end{array}$ (5)
where the (real) Monge-Ampere equation ( 5 ) substitutes for the Poisson equation ( 2 ). Here, ${D}^{2}\phi \left(t,x\right)$  stands for the $d×d$  symmetric matrix made of all second order $x-$  partial derivatives of $\phi$  , $\mathbb{I}$  stands for the $d×d$  identity matrix and $det$  for the determinant of a square matrix. The occurrence of the Monge-Ampere equation in mathematical modeling is not very common.
Notice, however, that a very similar system can be found in meteorology with Hoskins’ semi-geostrophic equations (cf. [?, [? and the included references).
In a simplified two dimensional setting, the semi-geostrophic equations read
 $\begin{array}{ccc}& & \frac{\partial \rho }{\partial t}+\left\{\phi ,\rho \right\}=0\end{array}$ (6)
 $\begin{array}{ccc}& & det\left(\mathbb{I}+{\epsilon }^{2}{D}^{2}\phi \right)=\rho ,\end{array}$ (7)
where $\left\{\cdot ,\cdot \right\}$  denotes the usual Poisson bracket.
Formally, as the coupling constant $\epsilon$  is small, the $VP$  and $VMA$  equations asymptotically approach each other up to order $O\left({\epsilon }^{4}\right)$  . Indeed, linearizing the determinant about the identity matrix leads to
 $\begin{array}{ccc}det\left(\mathbb{I}+{\epsilon }^{2}{D}^{2}\phi \right)=1+{\epsilon }^{2}\Delta \phi +O\left({\epsilon }^{4}\right).& & \end{array}$ (8)
The formal limit, as $\epsilon =0$  , reads
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+\xi \cdot {\nabla }_{x}f+{\nabla }_{x}\phi \cdot {\nabla }_{\xi }f=0\end{array}$ (9)
 $\begin{array}{ccc}& & \rho =1,\end{array}$ (10)
where constraint ( 10 ) substitutes for both the Poisson and the Monge-Ampere equations. The limit system ( 9 , 10 ), that we call constrained Vlasov system, can be seen as a ’kinetic’ extension of the Euler equations of classical incompressible fluid mechanics,
 $\begin{array}{ccc}& & {\partial }_{t}v+\left(v\cdot \nabla \right)v=-\nabla p\end{array}$ (11)
 $\begin{array}{ccc}& & \nabla \cdot v=0,\end{array}$ (12)
where $v\left(t,x\right)\in {\mathbb{R}}^{d}$  and $p\left(t,x\right)\in \mathbb{R}$  respectively are the velocity and the pressure of the fluid at time $t$  and position $x$  . Indeed, any smooth solution $\left(v,p\right)$  provides a ’monokinetic’ solution to the constrained Vlasov system ( 9 , 10 ), defined by $f\left(t,x,\xi \right)=\delta \left(\xi -v\left(t,x\right)\right),\phi =-p.$  Here a monokinetic solution means a delta-valued solution in the $\xi$  variable.
In addition, the constrained Vlasov system ( 9 , 10 ) turns out to be a natural extension (or $\Gamma$  limit) of the Euler equations from both geometrical and variational reasons, as explained in section  2 In a similar way, there is a monokinetic version of the $VP$  system, the so-called (pressureless) Euler-Poisson ( $EP$  ) system, which reads
 $\begin{array}{ccc}& & {\partial }_{t}v+\left(v\cdot \nabla \right)v=\nabla \phi \end{array}$ (13)
 $\begin{array}{ccc}& & {\partial }_{t}\rho +\nabla \cdot \left(\rho v\right)=0\end{array}$ (14)
 $\begin{array}{ccc}& & {\epsilon }^{2}\Delta \phi =\rho -1.\end{array}$ (15)
A rigorous asymptotic analysis of the $VMA$  system as $\epsilon \to 0$  will be provided (sections  5.1 and  5.2 ), in the case when the initial electronic density
 $\begin{array}{ccc}f\left(t=0,x,\xi \right)={f}^{0}\left(x,\xi \right)& & \end{array}$ (16)
is asymptotically monokinetic, namely approaching $\delta \left(\xi -{v}_{0}\left(x\right)\right),$  for some smooth divergence free velocity field ${v}_{0}$  , as $\epsilon$  tends to zero. Before this asymptotic analysis, we want to explain the geometric origin of the $VMA$  system. It has been known, since Arnold’s celebrated work (cf. [?), that the Euler equations (formally) describe geodesics curves along a suitable group of volume preserving maps, lengths being measured in the ${L}^{2}$  sense.
We will show (section  2 ) that the $VMA$  system just describes approximate geodesics obtained through a very natural penalty method, where $\epsilon$  stands for the penalty parameter. For this geometric interpretation to be valid, the Monge-Ampere equation ( 5 ) must be understood in the following weak sense:
for each fixed $t$  , $\phi \left(t,\cdot \right)$  is the unique (up to an additive constant) function such that $\Psi \left(x\right)={x}^{2}/2+{\epsilon }^{2}\phi \left(t,x\right)$  is convex in $x$  and
 $\begin{array}{ccc}\forall g\in {C}^{0}\left({\mathbb{R}}^{d}\right),{\int }_{{\mathbb{R}}^{d}}g\left(\nabla \Psi \left(x\right)\right)\rho \left(t,x\right)dx={\int }_{\Omega }g\left(y\right)dy,& & \end{array}$ (17)
where $\Omega$  is a fixed bounded open convex set where the neutralizing background of the electrons is assumed to be located. (This definition is made precise in section  2.3 .) Notice that, by construction, $\nabla \Psi$  must be valued in the closure of $\Omega$  and, therefore, the potential $\phi$  enjoys the following property $|x+{\epsilon }^{2}{\nabla }_{x}\phi \left(t,x\right)|\le {sup}_{y\in \Omega }|y|<+\infty .$  There is no similar bound for the electrostatic potential of the classical $VP$  system. Thus, in some sense, the $VMA$  system can be seen as a nonlinearly saturated version of the $VP$  system.
Beyond the geometric derivation of the $VMA$  system, our main analytic results are as follows:
• The $VMA$  system admits global energy preserving weak solutions.
• The $VMA$  system admits local strong solutions in periodic domains.
• For well prepared, nearly monokinetic initial data, the solutions of the $VMA$  system converge when $\epsilon$  goes to 0 to those of the Euler equations.
• In this asymptotic, the $EP$  system is a higher order approximation of the $VMA$  system.
The paper is organized as follows: in section 2, we first recall the geometric nature of the Euler equations, then we explain why the constrained Vlasov system ( 9 , 10 ) is a natural extension of the Euler equations from a variational point of view, finally we introduce the concept of approximate geodesics for volume preserving maps, and derive the $VMA$  system. Section 3 is devoted to the proof of existence of global energy preserving weak solutions. In section 4, we prove existence of local strong solutions, in the case of a periodic domain.
Finally, in section 5, we study the asymptotic behavior of the $VMA$  system as $\epsilon$  goes to 0.

2 The geometric origin of the Vlasov-Monge-Ampere system

2.1 The Euler equations

The motion of an incompressible fluid in a domain $\Omega \subset {\mathbb{R}}^{d}$  is classically described by the Euler equations $\left(E\right)$  :
 $\begin{array}{ccc}& & {\partial }_{t}v+\left(v\cdot \nabla \right)v=-\nabla p\end{array}$ (18)
 $\begin{array}{ccc}& & \nabla \cdot v=0,\end{array}$ (19)
with $t\in \mathbb{R}$  , $x\in \Omega$  , where $v=v\left(t,x\right)$  stands for the velocity field and $p=p\left(t,x\right)$  for the scalar pressure field. These equations have a nice geometrical interpretation going back to Arnold (see [?). Introducing $G\left(\Omega \right)$  the group of all volume preserving diffeomorphisms of $\Omega$  with jacobian determinant equal to 1, and measuring lengths in the ${L}^{2}$  sense, we may define (at least formally) geodesic curves along $G\left(\Omega \right)$  . It turns out that the Euler equations just describe these curves. For the same reasons, the Euler equations can be seen as the optimality equations for the corresponding minimization problem:
given two maps chosen in $G\left(\Omega \right)$  , find an ${L}^{2}-$  shortest path between them along $G\left(\Omega \right)$  . It was shown by Shnirelman [? (see also [? and [?) that, in the case when $\Omega$  is the unit cube in ${\mathbb{R}}^{3}$  , there are many maps for which there are no such shortest paths. Beyond this negative result, [? established that minimizing paths are more appropriately described by doubly stochastic measures. These measures (also called polymorphisms) generalize volume preserving maps in the following way: a doubly stochastic measure $\mu \left(dx,dy\right)$  is a (Borel) probability measure on $\Omega ×\Omega$  with two projections on each copy of $\Omega$  both equal to the (normalized) Lebesgue measure. It is known -see [?, for instancethat any such $\mu$  can be weakly approximated by a sequence ${\mu }_{n}\left(dx,dy\right)=\delta \left(x-{g}_{n}\left(y\right)\right)dy$  where each ${g}_{n}$  is a volume preserving map of $\Omega$  .
In [? it was shown that, in the case considered by Shnirelman for which there is no classical shortest path, minimizing paths along $G\left(\Omega \right)$  converge to paths of doubly stochastic measures $t\to \mu \left(t;dx,dy\right)$  governed by the following extension of the Euler equations
 $\begin{array}{ccc}& & {\partial }_{t}\mu +{\nabla }_{x}\cdot \left(\mu v\right)=0,\end{array}$ (20)
 $\begin{array}{ccc}& & {\partial }_{t}\left(v\mu \right)+{\nabla }_{x}\cdot \left(\mu v\otimes v\right)+\mu {\nabla }_{x}p=0,\end{array}$ (21)
where $v=v\left(t;x,y\right)$  and $p=p\left(t,x\right)$  can be respectively seen as the velocity field and the pressure field attached to $\mu$  . (Notice that the velocity field $v$  generally depends on the extra variable $y$  and is not a classical but rather a multivalued velocity field.) These equations are just a reformulation of the constrained Vlasov system ( 9 , 10 ). Indeed, it can be checked, under appropriate regularity assumptions, that the kinetic measure $f$  defined by
 $\begin{array}{ccc}f\left(t;dx,d\xi \right)={\int }_{y\in \Omega }\delta \left(\xi -v\left(t;x,y\right)\right)\mu \left(t;dx,dy\right)& & \end{array}$ (22)
solves ( 9 , 10 ) when $\left(\mu ,v,p\right)$  solves ( 20 , 21 ). Thus we conclude that the constrained Vlasov system ( 9 , 10 ) is a natural variational extension of the Euler equations.

2.2 Approximate geodesics

A general strategy to define approximate geodesics along a manifold $M$  (in our case $M=G\left(\Omega \right)$  ) embedded in a Hilbert space $H$  (here $H={L}^{2}\left(\Omega ,{\mathbb{R}}^{d}\right)$  ) is to introduce a penalty parameter $\epsilon >0$  and the following $unconstrained$  dynamical system in $H$
 $\begin{array}{ccc}{\partial }_{tt}X+\frac{1}{2{\epsilon }^{2}}{\nabla }_{X}\left({d}^{2}\left(X,M\right)\right)\right)=0.& & \end{array}$ (23)
In this equation, the unknown $t\to X\left(t\right)$  is a curve in $H$  , $d\left(X,M\right)$  is the distance (in $H$  ) of $X$  to the manifold $M$  , i.e. in our case as $M=G\left(\Omega \right)$  ,
 $\begin{array}{ccc}d\left(X,G\left(\Omega \right)\right)={inf}_{g\in G\left(\Omega \right)}\parallel X-g{\parallel }_{H},& & \end{array}$ (24)
and, finally, ${\nabla }_{X}$  denotes the gradient operator in $H$  . This penalty approach has been used for the Euler equations by the first author in [?. It is similar-but not identicalto Ebin’s slightly compressible flow theory [?, and is a natural extension of the theory of constrained finite dimensional mechanical systems [?. The penalized system is formally hamiltonian in variables $\left(X,{\partial }_{t}X\right)$  with Hamiltonian (or energy) given by:
 $\begin{array}{ccc}E=\frac{1}{2}\parallel {\partial }_{t}X{\parallel }_{H}^{2}+\frac{1}{2{\epsilon }^{2}}{d}^{2}\left(X,G\left(\Omega \right)\right).& & \end{array}$
(Multiplying equation ( 23 ) by ${\partial }_{t}X$  , we formally get that the energy is conserved.) Therefore it is plausible that the map $X\left(t\right)$  will remain close to $G\left(\Omega \right)$  if properly initialized at $t=0$  . A formal computation shows that, given a point $X$  for which there is a unique closest point ${\pi }_{X}$  to $X$  in the $H$  closure of $G\left(\Omega \right)$  , we have:
 $\begin{array}{ccc}{\nabla }_{X}\left(d\left(X,G\right)\right)=\frac{1}{d\left(X,G\right)}\left(X-{\pi }_{X}\right).& & \end{array}$ (25)
Thus the equation ( 23 ) formally becomes:
 $\begin{array}{ccc}{\partial }_{tt}X+\frac{1}{{\epsilon }^{2}}\left(X-{\pi }_{X}\right)=0.& & \end{array}$ (26)
To understand why solutions to such a system may approach geodesics along $G\left(\Omega \right)$  as $\epsilon$  goes to 0, just recall that, in the simple framework of a surface $S$  embedded in the 3 dimensional Euclidean space, a geodesic $t\to s\left(t\right)$  along $S$  is characterized by the fact that for every $t$  , the plane defined by $\left\{\stackrel{˙}{s}\left(t\right),\stackrel{¨}{s}\left(t\right)\right\}$  is orthogonal to $S$  . In our case, ${\partial }_{tt}X\left(t\right)$  is nearly orthogonal to $G\left(\Omega \right)$  thanks to ( 26 ), meanwhile $X\left(t\right)$  remains close to $G\left(\Omega \right)$  .
The approximate geodesic equation was introduced in [? in order to allow a spatial approximation of $G\left(\Omega \right)$  by the group of permutations of $N$  points ${A}_{j}$  chosen to form a discrete grid on $\Omega$  . On such a discrete group, the concept of geodesics becomes unclear meanwhile approximate geodesics still make sense.
They can be interpreted as trajectories of a cloud of $N$  particles ${X}_{i}$  moving in the Euclidean space ${\mathbb{R}}^{dN}$  , which substitutes for $H$  . These particles solve the following coupled system of harmonic oscillators ${\epsilon }^{2}\frac{{d}^{2}{X}_{i}}{d{t}^{2}}+{X}_{i}-{A}_{{\sigma }_{i}}=0,$  where $\sigma$  is a time dependent permutation minimizing, at each fixed time $t$  , $\Sigma {|{X}_{i}-{A}_{\sigma \left(i\right)}|}^{2}$  among all other permutations of the first $N$  integers. The convergence of this discrete model to the incompressible Euler equations for well prepared initial data was proved in [?. In order to study the continuous version ( 26 ), a specific study of the projection problem ( 24 ) is needed.

2.3 The polar decomposition Theorem

Let us first recall a general measure theoretic definition:
Definition 2.1 Let $A$  and $B$  be two topological spaces, let $\rho$  be a Borel finite measure of $A$  and $X$  a Borel map $A\to B$  , we call the push-forward of $\rho$  by $X$  and note $X#d\rho$  the Borel measure $\eta$  on $B$  defined by
 $\begin{array}{ccc}\forall f\in {C}^{0}\left(B\right),{\int }_{B}f\left(y\right)d\eta \left(y\right)={\int }_{A}f\left(X\left(x\right)\right)d\rho \left(x\right).& & \end{array}$
Let us now consider the case of a bounded open subset $\Omega$  of the Euclidean space ${\mathbb{R}}^{d}$  equipped with the Lebesgue measure that we denote $dx$  . We say that a Borel map $s:\overline{\Omega }\to \overline{\Omega }$  is volume (or Lebesgue measure) preserving if $s#dx=dx$  , i.e. if for all $g\in {C}^{0}\left(\overline{\Omega }\right)$  one has ${\int }_{\Omega }g\left(x\right)dx={\int }_{\Omega }g\left(s\left(x\right)\right)dx,$  or equivalently, for any Borel subset $B$  of $\overline{\Omega }$  one has $|{s}^{-1}\left(B\right)|=|B|$  . The set of all measure preserving maps of $\Omega$  is a closed subset of the Hilbert space $H={L}^{2}\left(\Omega ,{\mathbb{R}}^{d}\right)$  and will be denoted by $S\left(\Omega \right)$  . Notice that $S\left(\Omega \right)$  is only a semi-group for the composition rule and contains the group of volume preserving diffeomorphisms $G\left(\Omega \right)$  . It is known [? that, at least in the case when $\Omega$  is convex and $d\ge 2$  , $S\left(\Omega \right)$  is exactly the closure of $G\left(\Omega \right)$  in ${L}^{2}\left(\Omega ,{\mathbb{R}}^{d}\right)$  , which implies $d\left(.,G\left(\Omega \right)\right)=d\left(.,S\left(\Omega \right)\right).$  The polar decomposition Theorem for maps [? (extended to Riemannian manifolds in [?) will be crucial for our analysis of the $VMA$  system:
Theorem 2.2 Let $\Omega$  be a bounded convex open subset of ${\mathbb{R}}^{d}$  , let $X\in {L}^{2}\left(\Omega ;{\mathbb{R}}^{d}\right)$  and ${\rho }_{X}=X#dx,$  where $dx$  is the Lebesgue measure on $\Omega .$  Assume ${\rho }_{X}$  to be a Lebesgue integrable function, or, equivalently, $X$  to satisfy the non-degeneracy condition:
 $\begin{array}{ccc}\forall E\subset {\mathbb{R}}^{d}\text{Borel},|E|=0⇒|{X}^{-1}\left(E\right)|=0.& & \end{array}$ (27)
Then there exists a unique pair $\left(\nabla {\Phi }_{X},{\pi }_{X}\right)$  where ${\Phi }_{X}$  is a convex function and ${\pi }_{X}\in S\left(\Omega \right)$  , such that
 $\begin{array}{ccc}X=\nabla {\Phi }_{X}\circ {\pi }_{X}.& & \end{array}$ (28)
In this ’polar decomposition’, ${\pi }_{X}$  is also characterized as the unique closest point to $X$  on $S\left(\Omega \right)$  in the ${L}^{2}$  sense and ${\Phi }_{X}$  is characterized to be (up to an additive constant) the unique convex function on $\Omega$  satisfying
 $\begin{array}{ccc}{\int }_{{\mathbb{R}}^{d}}g\left(x\right)d{\rho }_{X}={\int }_{\Omega }g\left(X\left(y\right)\right)dy={\int }_{\Omega }g\left(\nabla {\Phi }_{X}\left(y\right)\right)dy,& & \end{array}$ (29)
for any $g\in {C}^{0}\left({\mathbb{R}}^{d}\right)$  such that $|g\left(x\right)|\le C\left(1+|x{|}^{2}\right)$  .
In addition, the Legendre-Fenchel transform ${\Psi }_{X}$  of ${\Phi }_{X}$  defined by
 $\begin{array}{ccc}{\Psi }_{X}\left(x\right)={sup}_{y\in \Omega }\left\{x\cdot y-{\Phi }_{X}\left(y\right)\right\}& & \end{array}$ (30)
is Lipschitz continuous on ${\mathbb{R}}^{d}$  , with Lipschitz constant bounded by ${sup}_{x\in \Omega }|x|$  and has the following properties :
$\nabla {\Psi }_{X}\left(x\right)\in \Omega$  holds true for ${\rho }_{X}$  a.e. $x$  ,
 $\begin{array}{ccc}{\int }_{{\mathbb{R}}^{d}}g\left(\nabla {\Psi }_{X}\right){\rho }_{X}\left(x\right)dx={\int }_{\Omega }g\left(\nabla {\Psi }_{X}\left(X\left(x\right)\right)\right)dx={\int }_{\Omega }g\left(x\right)dx& & \end{array}$ (31)
for any $g\in {C}^{0}\left(\overline{\Omega }\right)$  , and
 $\begin{array}{ccc}& & \nabla {\Phi }_{X}\left(\nabla {\Psi }_{X}\left(x\right)\right)=x{\rho }_{X}\left(x\right)dxa.e,\end{array}$ (32)
 $\begin{array}{ccc}& & \nabla {\Psi }_{X}\left(\nabla {\Phi }_{X}\left(y\right)\right)=ydya.e,\end{array}$ (33)
 $\begin{array}{ccc}& & {\pi }_{X}\left(y\right)=\nabla {\Psi }_{X}\left(X\left(y\right)\right)dya.e.\end{array}$ (34)
We make here several remarks on Theorem  2.2 :
Link with the Monge-Ampere equation We can interpret ( 29 ) as a weak version of the Monge-Ampere equation:
 $\begin{array}{ccc}{\rho }_{X}\left(\nabla \Phi \right)det{D}^{2}\Phi =1& & \end{array}$
and ( 31 ) can be seen as a weak version of another Monge-Ampere equation:
 $\begin{array}{ccc}& & det{D}^{2}\Psi ={\rho }_{X}\end{array}$
 $\begin{array}{ccc}& & \nabla \Psi \text{maps supp}\left({\rho }_{X}\right)\text{in}\Omega .\end{array}$
The pair $\left({\Phi }_{X},{\Psi }_{X}\right)$  depends in fact only of $\Omega$  and the measure ${\rho }_{X}=X#dx,$  and if condition ( 27 ) fails, then existence and uniqueness of the projection ${\pi }_{X}$  may fail, but existence and uniqueness of $\nabla {\Phi }_{X}$  remain true.
Theorem  2.2 and the subsequent remarks allow us to introduce the following notation that will be used throughout the paper:
Definition 2.3 Let $\Omega$  be a fixed bounded convex open set of ${\mathbb{R}}^{d}$  , let $\rho$  be a positive measure on ${\mathbb{R}}^{d}$  of total mass $|\Omega |$  , absolutely continuous w.r.t the Lebesgue measure and such that $\int \left(1+|x{|}^{2}\right)d\rho \left(x\right)<+\infty$  . We call $\Phi \left[\Omega ,\rho \right]$  , or, in short, $\Phi \left[\rho \right]$  , the unique up to a constant convex function on $\Omega$  satisfying
 $\begin{array}{ccc}\forall g\in {C}^{0}\left({\mathbb{R}}^{d}\right)\cap {L}^{1}\left(d\rho \right),{\int }_{{\mathbb{R}}^{d}}g\left(x\right)d\rho \left(x\right)={\int }_{\Omega }g\left(\nabla \Phi \left[\Omega ,\rho \right]\left(y\right)\right)dy.& & \end{array}$ (35)
We call $\Psi \left[\Omega ,\rho \right]$  its Legendre-Fenchel transform satisfying
 $\begin{array}{ccc}\forall g\in {C}^{0}\left({\mathbb{R}}^{d}\right)\cap {L}^{1}\left(\Omega ,dy\right),{\int }_{{\mathbb{R}}^{d}}g\left(\nabla \Psi \left[\Omega ,\rho \right]\left(x\right)\right)d\rho \left(x\right)={\int }_{\Omega }g\left(y\right)dy.& & \end{array}$ (36)
If no confusion is possible we may write $\Phi$  (resp. $\Psi$  ) instead of $\Phi \left[\Omega ,\rho \right]$  (resp.
$\Psi \left[\Omega ,\rho \right]\right).$  We will use some additional results from [?. The first one establishes the continuity of the polar decomposition:
Theorem 2.4 Let $\rho$  be a Lebesgue integrable positive measure on ${\mathbb{R}}^{d}$  , with total mass $\Omega$  , such that $\int \left(1+|x{|}^{2}\right)d\rho <+\infty$  . Let ${\rho }_{n}$  be a sequence of Lebesgue integrable positive measures on ${\mathbb{R}}^{d}$  , with total mass $\Omega$  , such that $\forall n$  , $\int \left(1+|x{|}^{2}\right)d{\rho }_{n}<+\infty$  . Let ${\Phi }_{n}=\Phi \left[\Omega ,{\rho }_{n}\right]$  and ${\Psi }_{n}=\Psi \left[\Omega ,{\rho }_{n}\right]$  be as in Definition  2.3 . If for any $f\in {C}^{0}\left({\mathbb{R}}^{d}\right)$  such that $|f\left(x\right)|\le C\left(1+|x{|}^{2}\right),$  $\int fd{\rho }_{n}$  converges to $\int fd\rho$  , then
• ${\Phi }_{n}$  converges to $\Phi \left[\Omega ,\rho \right]$  uniformly on each compact set of $\Omega$  and strongly in ${W}^{1,1}\left(\Omega \right),$
• ${\Psi }_{n}$  converges to $\Psi \left[\Omega ,\rho \right]$  uniformly on each compact set of ${\mathbb{R}}^{d}$  and strongly in ${W}^{1,1}\left(K\right)$  for every $K$  compact in ${\mathbb{R}}^{d}$  .
The second one provides a ’dual’ definition of the distance between a map $X$  and the semi-group $S\left(\Omega \right)$  :
Theorem 2.5 Let $X\in {L}^{2}\left(\Omega ;{\mathbb{R}}^{d}\right)$  and $\rho =X#dx,$  where $dx$  is the Lebesgue measure on $\Omega .$  Assume $\rho$  to be a Lebesgue integrable function. Then $\frac{1}{2}{d}^{2}\left(X,S\left(\Omega \right)\right)=\int \left(|x{|}^{2}/2-\Psi \left[\Omega ,\rho \right]\left(x\right)\right)\rho \left(x\right)dx+{\int }_{\Omega }\left(|y{|}^{2}/2-\Phi \left[\Omega ,\rho \right]\left(y\right)\right)dy$  $={sup}_{u,v}\int \left(|x{|}^{2}/2-u\left(x\right)\right)\rho \left(x\right)dx+{\int }_{\Omega }\left(|y{|}^{2}/2-v\left(y\right)\right)dy,$  where the supremum if performed over all pairs $\left(u,v\right)$  of continuous functions on ${\mathbb{R}}^{d}$  such that $u\left(x\right)+v\left(y\right)\ge x\cdot y$  pointwise.

2.4 The Vlasov-Monge-Ampere system

Let us now derive the $VMA$  system as the kinetic formulation of the approximate geodesic equation ( 26 ). First, from the polar decomposition Theorem  2.2 , equation ( 26 ) reads
 $\begin{array}{ccc}{\partial }_{tt}X\left(t,x\right)=\nabla \phi \left(t,X\left(t,x\right)\right),& & \end{array}$ (37)
where
 $\begin{array}{ccc}\nabla \phi \left(t,x\right)=\frac{\nabla \Psi \left[\Omega ,\rho \left(t,\cdot \right)\right]\left(x\right)-x}{{\epsilon }^{2}}& & \end{array}$ (38)
and $\Psi \left[\Omega ,\rho \right]$  is as in Definition ( 2.3 ). This means that $\nabla \phi$  satisfies ( 5 ) in a weak form with the additional condition that the range of $x\to x+{\epsilon }^{2}\nabla \phi \left(t,x\right)$  is contained in $\overline{\Omega }$  .
Next, let ${f}^{0}\ge 0$  be a given initial density function, that we assume to be in ${L}^{\infty }\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)$  , compactly supported and satisfying the compatibility condition
 $\begin{array}{ccc}\int {f}^{0}\left(x,\xi \right)dxd\xi =|\Omega |.& & \end{array}$ (39)
For each $t\ge 0$  , let us define $\left(x,\xi \right)\to f\left(t,x,\xi \right)$  to be ${f}^{0}$  pushed forward by the following ODE
 $\begin{array}{ccc}& & {\partial }_{t}X\left(t,x,\xi \right)=\Xi \left(t,x,\xi \right)\end{array}$ (40)
 $\begin{array}{ccc}& & {\partial }_{t}\Xi \left(t,x,\xi \right)=\left(\nabla \phi \right)\left(X\left(t,x,\xi \right)\right)\end{array}$ (41)
 $\begin{array}{ccc}& & \left(X,\Xi \right)\left(t=0,x,\xi \right)=\left(x,\xi \right).\end{array}$ (42)
Then $f$  satisfies the following kinetic (or Liouville) equation
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+{\nabla }_{\xi }\cdot \left(\nabla \phi f\right)=0\end{array}$ (43)
 $\begin{array}{ccc}& & f\left(0,\cdot ,\cdot \right)={f}^{0},\end{array}$ (44)
which must be understood in the following weak sense
 $\begin{array}{ccc}& & \forall g\in {C}_{c}^{\infty }\left(\left[0,+\infty \right)×{\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right),\end{array}$
 $\begin{array}{ccc}& & {\int }_{0}^{\infty }dt{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\left(\frac{\partial g}{\partial t}+\xi \cdot {\nabla }_{x}g+\nabla \phi \cdot {\nabla }_{\xi }g\right)fdxd\xi \end{array}$
 $\begin{array}{ccc}& & =-{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}{f}_{0}\left(x,\xi \right)g\left(t=0,x,\xi \right)dxd\xi .\end{array}$ (45)
This linear Liouville equation is nonlinearly coupled to equation ( 38 ), where $\rho$  is linked to $f$  by equation ( ??  ). Finally, we have defined, through ( 38 , 43 , 44 ), the weak formulation of the $VMA$  initial value problem.
The energy of the system is defined by
 $\begin{array}{ccc}E\left(t\right)=& & \frac{1}{2}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}f\left(t,x,\xi \right)|\xi {|}^{2}dxd\xi \end{array}$
 $\begin{array}{ccc}+& & \frac{1}{2{\epsilon }^{2}}{\int }_{{\mathbb{R}}^{d}}\rho \left(t,x\right){|\nabla \Psi \left[\Omega ,\rho \right]\left(t,x\right)-x|}^{2}dx.\end{array}$ (46)

3 Existence of global renormalized weak solutions

The main result of this section is as follows:
Theorem 3.1 Let $\left(x,\xi \right)\to {f}^{0}\left(x,\xi \right)\ge 0$  be in ${L}^{\infty }\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)$  , with compact support in both $x$  and $\xi$  , satisfying condition ( 39 ).
Then the $VMA$  system ( 38 , 43 , 44 ) admits a global weak solution $\left(f,\rho ,\Psi \right)$  , with $f\in {L}^{\infty }\left({\mathbb{R}}^{+}×{\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)$  and $\left(\rho ,\nabla \psi \right)\in {L}^{\infty }\left(\left[0,T\right]×{\mathbb{R}}^{d}\right)$  for all $T>0$  . In addition, each such weak solution enjoys the following properties:
• $f$  is a continuous function of $t$  , valued in ${L}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)$  , for every $1\le p<\infty$
• the density $\rho$  is a continuous function of $t$  , valued in ${L}^{p}\left({\mathbb{R}}^{d}\right)$  , for every $1\le p<\infty$  ,
• the support of $f\left(t,\cdot ,\cdot \right)$  in $\left(x,\xi \right)$  is compact, with a diameter growing no more than linearly in $t$  .
• the total energy defined by ( 46 ) is conserved,
• the ’renormalization’ property (in the sense of [?) $\frac{\partial g\left(f\right)}{\partial t}+{\nabla }_{x}\cdot \left(\xi g\left(f\right)\right)+{\nabla }_{\xi }\cdot \left(\nabla \phi g\left(f\right)\right)=0$  holds true for all $g\in {C}^{1}\left(\mathbb{R}\right)$  ,
• the trajectories of ( 41 , 42 ) are uniquely defined for almost every initial condition $\left(x,\xi \right)$  ,
• $t\to f\left(t,\cdot ,\cdot \right)$  is just ${f}^{0}$  pushed forward along the trajectories of ( 41 , 42 ).
Proof of Theorem  3.1 :
We build a sequence of approximate solutions $\left({f}_{h},{\Psi }_{h}{\right)}_{h>0}$  by time discretization and let the time step $h$  go to zero. To handle the limiting process, the non-linear terms will be treated with the help of Theorem  2.4 . More precisely if one can extract a subsequence such that, for every $t$  , ${f}_{h}\left(t,\cdot ,\cdot \right)$  converges weakly, then we can deduce from Theorem  2.4 that the corresponding sequence $\nabla {\Psi }_{h}\left(t,\cdot \right)$  will converge strongly, and this will allow us pass to the limit in the nonlinear term.

3.1 Construction of a sequence of approximate solutions

We consider $\eta \in {C}_{c}^{\infty }\left({\mathbb{R}}^{d}\right)$  such that $\eta \ge 0$  , ${\int }_{{\mathbb{R}}^{d}}\eta =1$  and ${\eta }_{h}=\frac{1}{{h}^{d}}\eta \left(\frac{\cdot }{h}\right)$  .
We then seek approximate solutions as solutions of the approximate problem
 $\begin{array}{ccc}& & \frac{\partial {f}_{h}}{\partial t}+\xi \cdot {\nabla }_{x}{f}_{h}+\frac{\nabla {\Psi }_{h}\left(x\right)-x}{{\epsilon }^{2}}\cdot {\nabla }_{\xi }{f}_{h}=0\end{array}$ (47)
 $\begin{array}{ccc}& & {f}_{h}\left(0,x,\xi \right)={f}_{h}^{0}\left(x,\xi \right)={f}_{0}{*}_{x,\xi }{\eta }_{h}\otimes {\eta }_{h}\end{array}$ (48)
 $\begin{array}{ccc}& & {\Psi }_{h}\left(t\right)={\eta }_{h}*\Psi \left[\Omega ,\rho \left(t=nh\right)\right]\text{for}t\in \left[nh,\left(n+1\right)h\left[.\end{array}$ (49)
$\nabla {\Psi }_{h}$  being a smooth function of space this regularized equation admits a unique solution that one builds by the method of characteristics. Since the flow is divergence-free in the phase space, the solution ${f}_{h}$  satisfies
 $\begin{array}{ccc}\forall p\in \left[1,+\infty \right],\parallel {f}_{h}\left(t\right){\parallel }_{{L}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)}=\parallel {f}_{h}\left(0\right){\parallel }_{{L}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)}.& & \end{array}$ (50)
By construction (through Theorem  2.2 ), $\nabla {\Psi }_{h}$  is valued in the convex bounded set $\overline{\Omega }$  . Suppose that ${f}^{0}\left(x,\xi \right)$  vanishes outside of the set $\left\{{x}^{2}+{\epsilon }^{2}{\xi }^{2}\le {C}^{2}\right\}$  for some constant $C>0$  fixed and denote $R={sup}_{y\in \Omega }|y|$  . Then we have
Lemma 3.2 $\forall t\ge 0,{f}_{h}\left(t,\cdot ,\cdot \right)$  is supported in $\left\{\sqrt{{x}^{2}+{\epsilon }^{2}{\xi }^{2}}\le C+Rt/\epsilon \right\}$  .
Proof : We just write ${\epsilon }^{2}{\partial }_{tt}X+X=\nabla {\Psi }_{h}\left(X\right)$  in complex notation $-i\epsilon {\partial }_{t}Z+Z=F$  , where $Z=X+i\epsilon {\partial }_{t}X$  and $F=\nabla {\Psi }_{h}\left(X\right)$  , which is bounded by $R$  . This leads to $Z\left(t\right)=Z\left(0\right)exp\left(-it/\epsilon \right)+i{\epsilon }^{-1}{\int }_{0}^{t}exp\left(-i\left(t-s\right)/\epsilon \right)F\left(s\right)ds$  ant the desired bound easily follows. Notice here a sharp contrast with the classical $VP$  system, for which the $\xi -$  support of the solutions cannot be controlled so easily (except in the one dimensional case). $\square$  Convergence of the sequence of approximate solutions Using ( 50 ) and Lemma  3.2 there exists, for any $1  , up to the extraction of a subsequence, $f\in {L}^{p}\left(\left[0,T\right]×{\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)$  such that ${f}_{h}$  converges weakly to $f$  as $h\to 0$  .
It remains to show that the product ${f}_{h}\nabla {\Psi }_{h}$  converges to the good limit.
For this we need strong convergence of $\nabla {\Psi }_{h}.$  We already know that $\nabla {\Psi }_{h}\in {L}^{\infty }\left(\left[0,T\right]×{\mathbb{R}}^{d}\right).$  We claim that for all $t>0$  , $\nabla {\Psi }_{h}\left(t,\cdot \right)$  converges strongly to $\nabla \Psi \left(t,\cdot \right)$  in ${L}_{loc}^{q}\left({\mathbb{R}}^{d}\right)$  , $\forall q\in \left[1,+\infty \left[$  . Indeed, such a strong convergence of $\nabla {\Psi }_{h}$  follows from Theorem  2.4 provided that we have for all $t>0$  ,
 $\begin{array}{ccc}{\int }_{{\mathbb{R}}^{d}}g\left(x\right){\rho }_{h}\left(t,x\right)dx\to {\int }_{{\mathbb{R}}^{d}}g\left(x\right)\rho \left(t,x\right)dx,& & \end{array}$ (51)
for any $g\in {C}^{0}\left({\mathbb{R}}^{d}\right)$  such that $\int \left(1+|x{|}^{2}\right)g\left(x\right)dx<+\infty$  . Note first that from Lemma  3.2 , we can restrict ourselves here to test functions $g$  that are compactly supported. Then we show that the sequence ${\rho }_{h}$  is relatively compact in $C\left(\left[0,T\right],{L}^{p}\left({\mathbb{R}}^{d}\right)-w\right).$  This is done by the following lemma:
Lemma 3.3 For all $T>0$  , for all $p$  with $1\le p<\infty$  the sequence ${f}_{h}$  (resp. ${\rho }_{h}$  ) satisfies
• ${f}_{h}$  (resp. ${\rho }_{h}$  ) is a bounded sequence in ${L}^{\infty }\left(\left[0,T\right];{L}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)\right)$  (resp.
in ${L}^{\infty }\left(\left[0,T\right];{L}^{p}\left({\mathbb{R}}^{d}\right)\right)$  ,
• ${\partial }_{t}{f}_{h}$  (resp. ${\partial }_{t}{\rho }_{h}$  ) is a bounded sequence in ${L}^{\infty }\left(\left[0,T\right];{W}^{-1,p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)\right)\right)$  , (resp. in ${L}^{\infty }\left(\left[0,T\right];{W}^{-1,p}\left({\mathbb{R}}^{d}\right)\right)$  ,
and one can extract from ${f}_{h}$  (resp. from ${\rho }_{h}$  ) a subsequence converging in $C\left(\left[0,T\right],{L}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)-w\right)$  (resp. in $C\left(\left[0,T\right],{L}^{p}\left({\mathbb{R}}^{d}\right)-w\right)$  ).
Proof: the first point uses equation ( 50 ) and Lemma  3.2 . The second point uses equation ( 43 ) and the identity:
${\partial }_{t}{\rho }_{h}=-{\nabla }_{x}\cdot {\int }_{{\mathbb{R}}^{d}}\xi {f}_{h}d\xi ,$  with the fact that the ${f}_{h}$  are uniformly compactly supported in $x$  and $\xi$  (Lemma  3.2 ); the last point is a classical result of functional analysis (see [?for example). $\square$  This lemma and Lemma  3.2 yield ( 51 ). Then using Theorem  2.4 , with $\rho$  the limit of a subsequence of ${\rho }_{h}$  , we have convergence of the sequence $\nabla {\Psi }_{h}$  to $\nabla \Psi \left[\Omega ,\rho \right]$  in $C\left(\left[0,T\right],{L}^{p}\left({\mathbb{R}}^{d}\right)\right)$  . We have extracted a subsequence ${f}_{h}$  such that
• ${f}_{h}$  converges in $C\left(\left[0,T\right],{L}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)-w\right)$  for every $1\le p<\infty$  .
• ${\rho }_{h}$  converges in $C\left(\left[0,T\right],{L}^{p}\left({\mathbb{R}}^{d}\right)-w\right)$  for every $1\le p<\infty$  .
• $\nabla {\Psi }_{h}\left(t,\cdot \right)$  converges in ${L}^{p}\left({\mathbb{R}}^{d}\right)$  for every $t$  and for every $1\le p<\infty$  .
Thus the limit $\left(f,\nabla \Psi \right)$  satisfies equations ( 43 - 44 ) and the first part of Theorem  3.1 is proved.

3.2 Conservation of energy

We now give a rigorous proof of the conservation of energy following an argument going back to F. Otto (in an unpublished work on the semi-geostrophic equations). We recall the definition of the energy as
 $\begin{array}{ccc}E\left(t\right)=\frac{1}{2}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}f\left(t,x,\xi \right)|\xi {|}^{2}dxd\xi +\frac{1}{2{\epsilon }^{2}}{\int }_{{\mathbb{R}}^{d}}\rho \left(t,x\right)|\nabla \Psi \left(t,x\right)-x{|}^{2}dx.& & \end{array}$
We call the first term the kinetic energy ${E}_{c}$  and the second term, multiplied by ${\epsilon }^{2}$  , the (normalized) potential energy ${E}_{p}.$  We have
Proposition 3.4 Let $f$  be any solution of ( 43 ) such that on every interval $\left[0,T\right]$  , $f\left(t,\cdot ,\cdot \right)$  is uniformly compactly supported in $|x|,|\xi |\le R\left(T\right)$  for some function $R\left(T\right)$  . Then the energy of the solution $f$  is conserved.
Proof: From Theorem  2.5 , we know that ${E}_{p}\left(t\right)=\int \left(|x{|}^{2}/2-\Psi \left(t,x\right)\right)\rho \left(t,x\right)dx+{\int }_{\Omega }\left(|y{|}^{2}/2-\Phi \left(t,y\right)\right)dy$  $={sup}_{u,v}\int \left(|x{|}^{2}/2-u\left(x\right)\right)\rho \left(t,x\right)dx+{\int }_{\Omega }\left(|y{|}^{2}/2-v\left(y\right)\right)dy,$  where the supremum if performed over all pairs $\left(u,v\right)$  of continuous functions on ${\mathbb{R}}^{d}$  such that $u\left(x\right)+v\left(y\right)\ge x\cdot y$  pointwise. Thus for each $t,{t}_{0}\in \mathbb{R}+$  , we have
 $\begin{array}{ccc}{E}_{p}\left(t\right)& & \ge \int \left(|x{|}^{2}/2-\Psi \left({t}_{0},x\right)\right)\rho \left(t,x\right)dx+{\int }_{\Omega }\left(|y{|}^{2}/2-\Phi \left({t}_{0},y\right)\right)dy,\end{array}$
and this implies
 $\begin{array}{ccc}{E}_{p}\left(t\right)-{E}_{p}\left({t}_{0}\right)& & \ge {\int }_{{\mathbb{R}}^{d}}\left(|x{|}^{2}/2-\Psi \left({t}_{0},x\right)\right)\left(\rho \left(t,x\right)-\rho \left({t}_{0},x\right)\right)dx\end{array}$
 $\begin{array}{ccc}& & ={\int }_{{t}_{0}}^{t}{\int }_{{\mathbb{R}}^{d}}{\partial }_{t}\rho \left(s,x\right)\left(|x{|}^{2}/2-\Psi \left({t}_{0},x\right)\right)dxds\end{array}$
 $\begin{array}{ccc}& & ={\int }_{{t}_{0}}^{t}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\xi f\left(s,x,\xi \right)\left(x-\nabla \Psi \left({t}_{0},x\right)\right)dxd\xi ds.\end{array}$
Notice that the product in the second line is licit since ${\partial }_{t}\rho$  is in ${W}^{-1,p}$  for any $1\le p<\infty$  , $f\left(t,\cdot ,\cdot \right)$  and therefore $\rho \left(t,\cdot \right)$  are compactly supported in space uniformly on $\left[0,T\right]$  , and $\Psi -|x{|}^{2}/2$  is in ${W}_{loc}^{1,\infty }$  . Exchanging ${t}_{0}$  and $t$  we would have found
 $\begin{array}{ccc}{E}_{p}\left({t}_{0}\right)-{E}_{p}\left(t\right)\ge {\int }_{t}^{{t}_{0}}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\xi f\left(s,x,\xi \right)\left(x-\nabla \Psi \left(t,x\right)\right)dxd\xi ds,& & \end{array}$
moreover we have for the kinetic energy
 $\begin{array}{ccc}{\epsilon }^{2}\left({E}_{c}\left(t\right)-{E}_{c}\left({t}_{0}\right)\right)={\int }_{{t}_{0}}^{t}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\xi f\left(t,x,\xi \right)\cdot \left(\nabla \Psi \left(s,x\right)-x\right)dxd\xi ds.& & \end{array}$
Dividing by $t-{t}_{0},t>{t}_{0}$  we find
 $\begin{array}{ccc}& & {\epsilon }^{2}\frac{E\left(t\right)-E\left({t}_{0}\right)}{t-{t}_{0}}\end{array}$
 $\begin{array}{ccc}& & \ge \frac{1}{t-{t}_{0}}{\int }_{{t}_{0}}^{t}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\xi f\left(s,x,\xi \right)\cdot \left(\nabla \Psi \left(s,x\right)-\nabla \Psi \left({t}_{0},x\right)\right)dxd\xi ds\end{array}$
and
 $\begin{array}{ccc}& & {\epsilon }^{2}\frac{E\left(t\right)-E\left({t}_{0}\right)}{t-{t}_{0}}\end{array}$
 $\begin{array}{ccc}& & \le \frac{1}{t-{t}_{0}}{\int }_{{t}_{0}}^{t}{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\xi f\left(s,x,\xi \right)\cdot \left(\nabla \Psi \left(t,x\right)-\nabla \Psi \left(s,x\right)\right)dxd\xi ds.\end{array}$
We know from  3.1 that $\nabla \Psi \left(t,.\right)$  converges strongly in ${L}_{loc}^{p}\left({\mathbb{R}}^{d}\right),1\le p<\infty$  to $\nabla \Psi \left({t}_{0},.\right)$  as $t$  goes to ${t}_{0}$  , and so the right hand sides of the above inequalities converges to 0 and we conclude that ${lim}_{t>{t}_{0}}\frac{E\left(t\right)-E\left({t}_{0}\right)}{t-{t}_{0}}=0.$  We could take $t<{t}_{0}$  and find the same result. Finally we conclude that
 $\begin{array}{ccc}\frac{dE}{dt}\equiv 0.& & \end{array}$
$\square$

3.3 Renormalized solutions and existence of characteristics

The study of renormalized solutions for transport equations has been introduced in [? for vector fields in ${W}^{1,1}$  with bounded divergence. These results have been extended by Bouchut [? to the case of Vlasov-type equations with acceleration field in $BV$  (A recent result of L. Ambrosio, [?, has extended the existence of renormalized solutions to transport equations with vector fields in $BV$  and with bounded divergence). The fact that solutions of ( 43 ,  44 ) are renormalized solutions is an immediate consequence of the following theorem:
Theorem 3.5 $\left($  F. Bouchut $\right)$  Let $f\in {L}^{\infty }\left(\right]0,T\left[,{L}_{loc}^{\infty }\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)\right)$  satisfy
 $\begin{array}{ccc}\frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+{\nabla }_{\xi }\cdot \left(E\left(t,x\right)f\right)=0,& & \end{array}$
with $E\left(t,x\right)\in {L}^{1}\left(\right]0,T\left[;{L}_{loc}^{1}\left({\mathbb{R}}^{d}\right)\right)\cap {L}^{1}\left(\right]0,T\left[;B{V}_{loc}\left({\mathbb{R}}^{d}\right)\right)$  , then, for any $g\in {C}^{1}\left(\mathbb{R}\right)$  , $\frac{\partial g\left(f\right)}{\partial t}+{\nabla }_{x}\cdot \left(\xi g\left(f\right)\right)+{\nabla }_{\xi }\cdot \left(E\left(t,x\right)g\left(f\right)\right)=0,$  and for every $1\le p<\infty$  , $f$  belongs to $C\left(\right]0,T\left[,{L}_{loc}^{p}\left({\mathbb{R}}^{d}×{\mathbb{R}}^{d}\right)\right)$  .
In our case the BV bound on the acceleration $\nabla \Psi$  is a direct consequence of the fact that $\Psi$  is a globally Lipschitz convex function. This result implies the strong time continuity results for $f$  and $\rho$  in Theorem  3.1 . Finally, as in [?, it can be deduced from the renormalization property that 1) for almost every initial condition $\left(x,\xi \right)$  , there is a unique trajectory solving ( 41 , 42 ), 2) $t\to f\left(t\right)$  is just ${f}^{0}$  pushed forward along these trajectories.
A complete proof is given in appendix.
Remark: From the renormalization property it follows that, once the potential $\Psi \left(t,x\right)$  is known, there exists a unique solution to ( 43 ) in ${L}_{t,x,\xi }^{\infty }$  .
Of course, this does not imply at all the uniqueness of weak solutions to the Vlasov-Monge-Ampere system! This paragraph ends the proof of Theorem  3.1 .

4 Strong solutions

In this section we show existence of strong solutions over a finite time intervall. To do so, we need regularity estimates for solutions of Monge-Ampere equation. We will get rid of the difficulties that may arise at the free boundary of the set $\left\{\rho >0\right\}$  by considering the periodic case. Note that for the Vlasov-Poisson system, existence of global smooth solutions has been proved (see [?); in the present case, due to the non-linearity of the Monge-Ampere equation, we were only able to obtain a result for finite time.

4.1 The periodic Vlasov-Monge-Ampere system

Polar factorization of maps in a periodic domain The polar decomposition Theorem has been generalized in [? to general Riemannian manifolds, while the particular case of the flat torus ${\mathbb{T}}^{d}={\mathbb{R}}^{d}/{\mathbb{Z}}^{d}$  had been addressed in [?.
Definition 4.1 We say that a mapping $Y:{\mathbb{R}}^{d}\to {\mathbb{R}}^{d}$  is ${\mathbb{Z}}^{d}$  additive if the mapping $x\to Y\left(x\right)-x$  is ${\mathbb{Z}}^{d}$  periodic. The set of all measurable ${\mathbb{Z}}^{d}$  additive mappings is denoted by $\mathcal{P}$  . For each $x\in {\mathbb{R}}^{d}$  we call $\stackrel{^}{x}$  the class of $x$  in ${\mathbb{R}}^{d}/{\mathbb{Z}}^{d}$  , and for any $X\in \mathcal{P}$  , $\stackrel{^}{X}$  the mapping of ${\mathbb{T}}^{d}$  into itself defined by $\forall x\in {\mathbb{R}}^{d},\stackrel{^}{X}\left(\stackrel{^}{x}\right)=\stackrel{^}{X}\left(x\right).$
We may say if no confusion is possible additive instead of ${\mathbb{Z}}^{d}$  additive. Then the following theorem can be deduced from the results of [? and [?:
Theorem 4.2 Let $X:{\mathbb{R}}^{d}\to {\mathbb{R}}^{d}$  be additive and assume that ${\rho }_{X}=X#dx$  has a density in ${L}^{1}\left(\left[0,1{\right]}^{d}\right)$  . Then there exists a unique pair $\left(\nabla {\Phi }_{X},{\pi }_{X}\right)$  such that $X=\nabla {\Phi }_{X}\circ {\pi }_{X}$  where ${\Phi }_{X}$  is a convex function and ${\Phi }_{X}\left(x\right)-|x{|}^{2}/2$  is ${\mathbb{Z}}^{d}$  periodic, ${\pi }_{X}:{\mathbb{R}}^{d}\to {\mathbb{R}}^{d}$  is additive and ${\stackrel{^}{\pi }}_{X}$  is Lebesgue measure preserving in ${\mathbb{T}}^{d}$  . Moreover we have $\parallel X-{\pi }_{X}{\parallel }_{{L}^{2}\left(\left[0,1{\right]}^{d}\right)}=\parallel \stackrel{^}{X}-{\stackrel{^}{\pi }}_{X}{\parallel }_{{L}^{2}\left({\mathbb{T}}^{d}\right)}$  and, ${\Psi }_{X}$  denoting the Legendre transform of ${\Phi }_{X}$  , we have ${\pi }_{X}=\nabla {\Psi }_{X}\circ X.$
Remark : The pair $\left({\Phi }_{X},{\Psi }_{X}\right)$  is uniquely defined by the density ${\rho }_{X}=X#dx$  .
Notice that the periodicity of ${\Phi }_{X}\left(x\right)-|x{|}^{2}/2$  implies that $\nabla {\Phi }_{X}$  and $\nabla {\Psi }_{X}$  are ${\mathbb{Z}}^{d}$  additive, and that ${\Psi }_{X}-|x{|}^{2}/2$  is also ${\mathbb{Z}}^{d}$  periodic. As in the previous case, we introduce the following notation:
Definition 4.3 Let $\rho$  be a probability measure on ${\mathbb{T}}^{d}$  , with density in ${L}^{1}\left({\mathbb{T}}^{d}\right)$  . We denote $\Phi \left[\rho \right]$  (resp. $\Psi \left[\rho \right]$  ) the unique up to a constant convex function such that
 $\begin{array}{ccc}& & \Phi \left[\rho \right]-|x{|}^{2}/2\text{is}{\mathbb{Z}}^{d}\text{periodic},\end{array}$ (52)
 $\begin{array}{ccc}& & \forall f\in {C}^{0}\left({\mathbb{T}}^{d}\right),{\int }_{{\mathbb{T}}^{d}}f\left(\stackrel{^}{\nabla }\Phi \left[\rho \right]\left(x\right)\right)dx=\int fd\rho \end{array}$ (53)
(resp. its Legendre fenchel transform).
$\Psi \left[\rho \right]$  will thus be a generalized solution of the Monge-Ampere equation $det{D}^{2}\Psi =\rho .$  Next the results of Caffarelli ([?, [?, [?) on the regularity of solutions of the Monge-Ampere equation yield the following theorem:
Theorem 4.4 Let $\rho >0$  be a ${C}^{\alpha }\left({\mathbb{T}}^{d}\right)$  probability density on ${\mathbb{T}}^{d}$  , for some $\alpha \in \right]0,1\left[$  .
Then $\Psi =\Psi \left[\rho \right]$  (see Definition  4.3 ) is a classical solution of
 $\begin{array}{ccc}det{D}^{2}\Psi =\rho & & \end{array}$
and satisfies:
 $\begin{array}{ccc}& & \parallel \nabla \Psi \left(x\right)-x{\parallel }_{{L}^{\infty }}\le C\left(d\right)=\sqrt{d}/2\end{array}$
 $\begin{array}{ccc}& & \parallel {D}^{2}\Psi {\parallel }_{{C}^{\alpha }}\le K\left(m,M,\parallel \rho {\parallel }_{{C}^{\alpha }}\right)\end{array}$
where $m=inf\rho$  and $M=sup\rho$  .
This theorem is an adaptation of the regularity results stated above, whose complete proof is given in appendix.
The periodic Vlasov-Monge-Ampere system We now seek $f:\left(t,x,\xi \right)\in \left({\mathbb{T}}^{d}×{\mathbb{R}}^{d}×\left[0,T\right]\right)\to f\left(t,x,\xi \right)\in {\mathbb{R}}^{+}$  , for some $T>0$  , solution of the initial value problem for the periodic Vlasov-Monge-Ampere $\left(VM{A}_{p}\right)$  system
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+\frac{1}{{\epsilon }^{2}}{\nabla }_{\xi }\cdot \left(\left(\nabla \Psi \left[\rho \right]\left(x\right)-x\right)f\right)=0\end{array}$ (54)
 $\begin{array}{ccc}& & f\left(0,\cdot ,\cdot \right)={f}^{0},\end{array}$ (55)
for a given ${f}^{0}$  satisfying the compatibility condition
 $\begin{array}{ccc}\int {f}^{0}\left(x,\xi \right)dxd\xi =1.& & \end{array}$ (56)
The macroscopic density $\rho$  is still related to $f$  by equation ( 3 ), and $\Psi \left[\rho \right]$  is as in Definition  4.3 .

4.2 Existence of local strong solutions

We mention first that the proof of existence of global weak solutions adapts with minor changes to the periodic case, and that the obtain for the periodic $\left(VM{A}_{p}\right)$  system the same result as Theorem  3.1 .
Our result in this section is the following:
Theorem 4.5 Let ${f}_{0}\in {W}^{1,\infty }\left({\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right),$  be such that:
 $\begin{array}{ccc}& & \exists {C}_{0}>0:{f}_{0}\equiv 0\text{for}|\xi |\ge {C}_{0},\end{array}$ (57)
 $\begin{array}{ccc}& & \exists m>0:{\rho }_{0}\left(x\right)={\int }_{{\mathbb{R}}^{d}}{f}_{0}\left(x,\xi \right)d\xi \ge m\forall x\in {\mathbb{T}}^{d},\end{array}$ (58)
then there exists $T>0$  and a solution $f$  to the $VM{A}_{p}$  system ( 54 , 55 ), in the space ${W}^{1,\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)$  .
Proof of Theorem  4.5 : First we deduce from Theorem  4.4 :
Corollary 4.6 Let $\rho ,\Psi =\Psi \left[\rho \right]$  be as in Theorem  4.4 . Then, we have $\parallel {D}^{2}\Psi {\parallel }_{{L}^{\infty }\left({\mathbb{T}}^{d}\right)}\le C\left(m,M,\parallel \nabla \rho {\parallel }_{{L}^{\infty }\left({\mathbb{T}}^{d}\right)}\right),$  and we can define $K\left(m,M,l\right)=sup\left\{\parallel {D}^{2}\Psi \left[\rho \right]{\parallel }_{{L}^{\infty }\left({\mathbb{T}}^{d}\right)};\parallel \nabla \rho {\parallel }_{{L}^{\infty }\left({\mathbb{T}}^{d}\right)}\le l,m\le \rho \le M\right\}<\infty .$
We see that in order to use Theorem  4.4 we need $\rho$  to be bounded away from 0. In the following lemma, we show that under suitable assumptions on the initial data, it is possible to enforce locally in time the condition $0
Lemma 4.7 Let $f\in {L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)$  satisfy
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+{\nabla }_{\xi }\cdot \left(E\left(t,x\right)f\right)=0\end{array}$
 $\begin{array}{ccc}& & f\left(0,.,.\right)={f}^{0},\end{array}$
with $E\in {L}^{1}\left(\left[0,T\right],BV\left({\mathbb{T}}^{d}\right)\right)$  and
 $\begin{array}{ccc}\parallel E{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)}\le F,& & \end{array}$
let the initial condition ${f}_{0}$  be such that $a\left(x,\xi \right)\le f\left(0,x,\xi \right)\le b\left(x,\xi \right),$  with ${\rho }_{a}\left(x\right)=\int a\left(x,\xi \right)d\xi \ge m>0$  and ${\rho }_{b}\left(x\right)=\int b\left(x,\xi \right)d\xi \le M<\infty$  and $a,b$  satisfying
 $\begin{array}{ccc}|{\nabla }_{x,\xi }\left(a,b\right)|\le \frac{c}{1+|\xi {|}^{d+2}}.& & \end{array}$
Then there exists a constant $R>0$  depending on $m,M,c,F$  , such that
 $\begin{array}{ccc}\left({\rho }_{a}\left(x\right)-Rt\right)\le \rho \left(t,x\right)\le \left({\rho }_{b}\left(x\right)+Rt\right).& & \end{array}$
The proof of the lemma is given in appendix.

4.2.1 Construction of approximate solutions

Let us consider $\left(t,x\right)\to E\left(t,x\right)$  a smooth vector-field on ${\mathbb{T}}^{d}$  , and write
 $\begin{array}{ccc}{T}_{E}\left(f\right)=\frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+{\nabla }_{\xi }\cdot \left(Ef\right).& & \end{array}$
If $f$  satisfies ${T}_{E}\left(f\right)=0$  , we have
 $\begin{array}{ccc}& & {T}_{E}{\nabla }_{x}f=-\left({\nabla }_{x}E\right)\cdot {\nabla }_{\xi }f\end{array}$
 $\begin{array}{ccc}& & {T}_{E}{\nabla }_{\xi }f=-{\nabla }_{x}f\end{array}$
 $\begin{array}{ccc}& & {T}_{E}{\partial }_{t}f=-{\partial }_{t}E\cdot {\nabla }_{\xi }f,\end{array}$
and therefore
 $\begin{array}{ccc}\frac{d}{dt}\parallel {\nabla }_{x,\xi }f{\parallel }_{{L}^{\infty }}\le \parallel {\nabla }_{x,\xi }f{\parallel }_{{L}^{\infty }}\left(1+\parallel {\nabla }_{x}E{\parallel }_{{L}^{\infty }}\right)& & \end{array}$ (59)
which implies
 $\begin{array}{ccc}\parallel {\nabla }_{x,\xi }f\left(t\right){\parallel }_{{L}^{\infty }}\le \parallel {\nabla }_{x,\xi }f\left(t=0\right){\parallel }_{{L}^{\infty }}exp\left({\int }_{0}^{t}\left(1+\parallel {\nabla }_{x}E\left(s\right){\parallel }_{{L}^{\infty }}\right)ds\right).& & \end{array}$
Now let ${f}_{0}$  be given as in Theorem  4.5 , satisfying ( 57 , 58 ). Thanks to Lemma  4.7 it is possible to find ${t}_{1},m,M$  such that for any $f$  satisfying
 $\begin{array}{ccc}& & {T}_{E}\left(f\right)=0\end{array}$
 $\begin{array}{ccc}& & f\left(t=0\right)={f}_{0},\end{array}$
for any field $E\in {L}^{1}\left(\left[0,{t}_{1}\right],BV\left({\mathbb{T}}^{d}\right)\right)$  satisfying $\parallel E{\parallel }_{{L}^{\infty }\left(\left[0,{t}_{1}\right]×{\mathbb{T}}^{d}\right)}\le \sqrt{d}/\left(2{\epsilon }^{2}\right)$  , we have
 $\begin{array}{ccc}& & m\le \rho \left(t,\cdot \right)\le M,\forall t\in \left[0,{t}_{1}\right]\end{array}$ (60)
 $\begin{array}{ccc}& & |\xi {|}_{max}\le {C}_{1}=10{C}_{0},\end{array}$ (61)
with $f$  supported in $\left\{|\xi |\le |\xi {|}_{max}\right\}$  and with ${C}_{0}$  as in Theorem  4.5 , so that we have for $0\le t\le {t}_{1}$  :
 $\begin{array}{ccc}\parallel \nabla \rho {\parallel }_{{L}^{\infty }}\le {\omega }_{d}{C}_{1}^{d}\parallel {\nabla }_{x}f{\parallel }_{{L}^{\infty }},& & \end{array}$ (62)
${\omega }_{d}$  being the volume of the unit ball of ${\mathbb{R}}^{d}$  . Then we construct a family of approximate solutions $\left({f}_{h},{\Psi }_{h}\right)$  to ( 54 ), in the same spirit as we did for weak solutions, by solving
 $\begin{array}{ccc}& & \frac{\partial {f}_{h}}{\partial t}+\xi \cdot {\nabla }_{x}{f}_{h}+\frac{\nabla {\Psi }_{h}\left(x\right)-x}{{\epsilon }^{2}}\cdot {\nabla }_{\xi }{f}_{h}=0\end{array}$
 $\begin{array}{ccc}& & {f}_{h}\left(t=0\right)={f}_{0}\end{array}$
 $\begin{array}{ccc}& & {\Psi }_{h}\left(t\right)=\Psi \left[\rho \left(t=nh\right)\right]\text{for}t\in \left[nh,\left(n+1\right)h\left[.\end{array}$
Note that we have neither mollified the term $\nabla \Psi$  nor the initial condition and that $\parallel \nabla {\Psi }_{h}-x{\parallel }_{{L}^{\infty }}\le C\left(d\right)=\sqrt{d}/2$  . Now choose $l=10\parallel {\nabla }_{x,\xi }{f}_{0}{\parallel }_{{L}^{\infty }}{\omega }_{d}{C}_{1}^{d}$  , if for some $t=nh\le {t}_{1}-h$  we have
 $\begin{array}{ccc}\parallel {\nabla }_{x,\xi }{f}^{h}\left(t=nh\right){\parallel }_{{L}^{\infty }}\le \frac{l}{{\omega }_{d}{C}_{1}^{d}}& & \end{array}$
this implies, thanks to ( 62 ), that
 $\begin{array}{ccc}\parallel {\nabla }_{x}{\rho }^{h}\left(t=nh\right){\parallel }_{{L}^{\infty }}\le l,& & \end{array}$
and conditions ( 60 , 61 ) are satisfied because $t\le {t}_{1}$  . Then if we denote $K=K\left(m,M,l\right)$  as in Corollary  4.6 , we have for $nh\le t
 $\begin{array}{ccc}\frac{d}{dt}\parallel {\nabla }_{x,\xi }{f}^{h}{\parallel }_{{L}^{\infty }}\le \left(K+1\right)\parallel {\nabla }_{x,\xi }{f}^{h}{\parallel }_{{L}^{\infty }},& & \end{array}$
and then
 $\begin{array}{ccc}\parallel {\nabla }_{x,\xi }{f}^{h}\left(t=nh+h\right){\parallel }_{{L}^{\infty }}\le \parallel {\nabla }_{x,\xi }{f}^{h}\left(t=nh+h\right){\parallel }_{{L}^{\infty }}exp\left(K+1\right)h.& & \end{array}$
So if we define $T$  as
 $\begin{array}{ccc}T=min\left\{{t}_{1},{t}_{2}\right\},& & \end{array}$
with $exp\left(\left(K+1\right){t}_{2}\right)=10$  , we have for $0\le t\le T$  ,
 $\begin{array}{ccc}& & \parallel {\nabla }_{x,\xi }{f}^{h}{\parallel }_{{L}^{\infty }}\le 10\parallel {\nabla }_{x,\xi }{f}_{0}{\parallel }_{{L}^{\infty }}\end{array}$
 $\begin{array}{ccc}& & \parallel \nabla {\rho }^{h}{\parallel }_{{L}^{\infty }}\le l\end{array}$
 $\begin{array}{ccc}& & m\le \rho \le M\end{array}$
 $\begin{array}{ccc}& & \parallel {D}^{2}{\Psi }^{h}{\parallel }_{{L}^{\infty }}\le K.\end{array}$
Thus we can extract a subsequence converging to a strong solution of ( 54 , 55 ). Then we argue as in section 2 to show that all terms converge to the correct limit. This ends the proof of Theorem  4.5 .
$\square$

5 Asymptotic analysis

5.1 Convergence to the Euler equation

In this section we justify that the Vlasov-Monge-Ampere system describes approximate geodesics on volume preserving transformations: indeed we will show that weak solutions of this system converge to a solution of the incompressible Euler equations $\left(E\right)$  as the parameter $\epsilon$  goes to 0, at least for well prepared initial data. We restrict ourselves to the space periodic case, the macroscopic density $\rho$  is still defined by ( 3 ) and the convex potentials $\Phi \left[\rho \right],\Psi \left[\rho \right]$  are still as in Definition  4.3 .
For sake of simplicity, we slightly modify our notations and introduce the following rescaled potentials
 $\begin{array}{ccc}& & \stackrel{~}{\phi }\left[\rho \right]=\frac{|x{|}^{2}/2-\Psi \left[\rho \right]}{\epsilon },\end{array}$
 $\begin{array}{ccc}& & \phi \left[\rho \right]=\frac{\Phi \left[\rho \right]-|x{|}^{2}/2}{\epsilon },\end{array}$
so that
 $\begin{array}{ccc}\nabla \phi \left[\rho \right]=\nabla \stackrel{~}{\phi }\left[\rho \right]\circ \nabla \Phi \left[\rho \right],& & \end{array}$
and the $VM{A}_{p}$  system takes the following form:
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+\xi \cdot {\nabla }_{x}f-\frac{\nabla \stackrel{~}{\phi }\left[\rho \right]}{\epsilon }\cdot {\nabla }_{\xi }f=0\end{array}$ (63)
 $\begin{array}{ccc}& & f\left(0,\cdot ,\cdot \right)={f}_{0}.\end{array}$ (64)
The energy is given by
 $\begin{array}{ccc}E\left(t\right)& & =\frac{1}{2}\int f\left(t,x,\xi \right)|\xi {|}^{2}dxd\xi +\frac{1}{2}\int |\nabla \phi {|}^{2}dx\end{array}$ (65)
 $\begin{array}{ccc}& & =\frac{1}{2}\int f\left(t,x,\xi \right)|\xi {|}^{2}dxd\xi +\frac{1}{2}\int \rho |\nabla \stackrel{~}{\phi }{|}^{2}dx.\end{array}$
It has been shown in section  3.2 that the energy is conserved. The Euler equations for incompressible fluids $\left(E\right)$  reads:
 $\begin{array}{ccc}& & {\partial }_{t}v+v\cdot \nabla v=-\nabla p\end{array}$ (66)
 $\begin{array}{ccc}& & \nabla \cdot v=0.\end{array}$ (67)
We shall here consider a smooth solution of $E$  and a weak solution of $VM{A}_{p}$  , with ‘well prepared initial data’, meaning that the initial data of both systems are close a time 0. Then we will show that as time evolves, both solutions stay close to each other.
Theorem 5.1 Let $f$  be a weak solution of ( 63 ,  64 ) with finite energy, let $\left(t,x\right)\to v\left(t,x\right)$  be a smooth ${C}^{2}\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)$  solution of ( 66 , 67 ) for $t\in \left[0,T\right],$  and $p\left(t,x\right)$  the corresponding pressure, let
 $\begin{array}{ccc}{H}_{\epsilon }\left(t\right)=\frac{1}{2}\int f\left(t,x,\xi \right)|\xi -v\left(t,x\right){|}^{2}dxd\xi +\frac{1}{2}\int |\nabla \phi \left(t,x\right){|}^{2}dx,& & \end{array}$
then
 $\begin{array}{ccc}{H}_{\epsilon }\left(t\right)\le Cexp\left(Ct\right)\left({H}_{\epsilon }\left(0\right)+{\epsilon }^{2}\right),\forall t\in \left[0,T\right].& & \end{array}$
$C$  depends only on $T,{sup}_{0\le s\le T}\left\{\parallel v\left(s,.\right),p\left(s,.\right),{\partial }_{t}p\left(s,.\right),\nabla p\left(s,.\right){\parallel }_{{W}^{1,\infty }\left({\mathbb{T}}^{d}\right)}\right\}$  .
Remark 1: This estimate is enough to compare the weak solutions $f$  to the $VM{A}_{p}$  system (for well prepared initial data) and the smooth solutions $v$  of the Euler equations. For instance, $\int f\left(t=0,x,\xi \right)d\xi \equiv 1$  implies $\phi \left(t=0,x\right)\equiv 0$  and therefore, $\int |\xi -v\left(t=0,x\right){|}^{2}f\left(t=0,x,\xi \right)dxd\xi \le {C}_{0}{\epsilon }^{2}$  implies ${sup}_{t\in \left[0,T\right]}\int |\xi -v\left(t,x\right){|}^{2}f\left(t,x,\xi \right)dxd\xi \le {C}_{T}{\epsilon }^{2},$  where ${C}_{T}$  depends only on ${C}_{0}$  , $T$  and $v$  .
Remark 2: We see that we consider nearly monokinetic initial data for the $VM{A}_{p}$  system.
Proof of Theorem  5.1 We shall show that
 $\begin{array}{ccc}\frac{d}{dt}{H}_{\epsilon }=& & -\int f\left(t,x,\xi \right)\left(\xi -v\right)\nabla v\left(\xi -v\right)\end{array}$
 $\begin{array}{ccc}& & +\int f\left(t,x,\xi \right)\frac{1}{\epsilon }v\cdot \nabla \stackrel{~}{\phi }\end{array}$
 $\begin{array}{ccc}& & -\int f\left(t,x,\xi \right)\left(v-\xi \right)\cdot \nabla p,\end{array}$ (68)
where we will use the notation
 $\begin{array}{ccc}u\nabla vw={\sum }_{i,j=1}^{d}{u}^{i}{\partial }_{i}{v}^{j}{w}^{j}.& & \end{array}$
The proof of this identity is postponed to the end of the section.
Now we look at all terms of the right hand side. All the constants that we denote by $C$  are controlled as in Theorem  5.1 . We set
 $\begin{array}{ccc}& & {T}_{1}=-\int f\left(t,x,\xi \right)\left(\xi -v\right)\nabla v\left(\xi -v\right),\end{array}$
 $\begin{array}{ccc}& & {T}_{2}=\int f\left(t,x,\xi \right)\frac{1}{\epsilon }v\cdot \nabla \stackrel{~}{\phi },\end{array}$
 $\begin{array}{ccc}& & {T}_{3}=-\int f\left(t,x,\xi \right)\left(v-\xi \right)\cdot \nabla p.\end{array}$
First we have ${T}_{1}\le C{H}_{\epsilon }.$  For ${T}_{2}$  we have
 $\begin{array}{ccc}{T}_{2}=\frac{1}{\epsilon }\int \rho v\cdot \nabla \stackrel{~}{\phi }=& & \frac{1}{\epsilon }\int v\left(\nabla \Phi \left[\rho \right]\right)\cdot \nabla \stackrel{~}{\phi }\left(\nabla \Phi \left[\rho \right]\right)\end{array}$
 $\begin{array}{ccc}=& & \frac{1}{\epsilon }\int v\left(x+\epsilon \nabla \phi \right)\cdot \nabla \phi \end{array}$
 $\begin{array}{ccc}=& & \frac{1}{\epsilon }\int v\cdot \nabla \phi +\left(v\left(x+\epsilon \nabla \phi \right)-v\left(x\right)\right)\cdot \nabla \phi \end{array}$
 $\begin{array}{ccc}\le & & 0+C\int {|\nabla \phi |}^{2}\le C{H}_{\epsilon },\end{array}$
we have used that $v$  is divergence-free thus its integral against any gradient is zero. Next we have the following lemma:
Lemma 5.2 Let $G:{\mathbb{T}}^{d}\to \mathbb{R}$  be Lipschitz continuous such that ${\int }_{{\mathbb{T}}^{d}}G=0$  , then for all $R>0$  , one has $|\int \rho G|\le \frac{1}{2}\parallel \nabla G{\parallel }_{{L}^{\infty }}\left(\frac{1}{R}{\epsilon }^{2}+R{H}_{\epsilon }\right).$
Proof: We just write a Taylor expansion of $G$  :
 $\begin{array}{ccc}& & |\int \left(\rho -1\right)G|=|\int \left(G\left(x+\epsilon \nabla \phi \right)-G\left(x\right)|\end{array}$
 $\begin{array}{ccc}\le & & \epsilon \parallel \nabla G{\parallel }_{{L}^{\infty }}\parallel \nabla \phi {\parallel }_{{L}^{1}}\le \epsilon \parallel \nabla G{\parallel }_{{L}^{\infty }}{H}_{\epsilon }^{1/2}\le \frac{1}{2}\parallel \nabla G{\parallel }_{{L}^{\infty }}\left(\frac{1}{R}{\epsilon }^{2}+R{H}_{\epsilon }\right).\end{array}$
$\square$  Again, since $v$  is divergence-free, $\int v\cdot \nabla p=0,$  thus from Lemma  5.2 we have
 $\begin{array}{ccc}-\int \rho v\cdot \nabla p\le C\left({\epsilon }^{2}+{H}_{\epsilon }\right).& & \end{array}$
We remind that ${\partial }_{t}\rho \left(t,x\right)=-{\nabla }_{x}\cdot \int f\left(t,x,\xi \right)\xi d\xi .$  Since it costs no generality to suppose that for all $t\in \left[0,T\right]$  , $\int p\left(t,x\right)dx\equiv 0$  , we obtain that
 $\begin{array}{ccc}\int f\left(t,x,\xi \right)\xi \cdot \nabla p& =& \int \frac{\partial \rho }{\partial t}p\end{array}$
 $\begin{array}{ccc}& =& \frac{d}{dt}\int \rho p-\int \rho \frac{\partial p}{\partial t}\end{array}$
 $\begin{array}{ccc}& \le & C\left({\epsilon }^{2}+{H}_{\epsilon }\right)-\frac{dQ}{dt}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
again using Lemma  5.2 , where $Q\left(t\right)=-\int \rho p.$  Thus
 $\begin{array}{ccc}{T}_{3}\le C\left({H}_{\epsilon }+{\epsilon }^{2}\right)-\frac{dQ}{dt}& & \end{array}$
and we have the following inequality:
 $\begin{array}{ccc}\frac{d}{dt}\left({H}_{\epsilon }+Q\right)\le C{H}_{\epsilon }+O\left({\epsilon }^{2}\right).& & \end{array}$ (69)
Moreover, using Lemma  5.2 ,
 $\begin{array}{ccc}|Q\left(t\right)|\le C{\epsilon }^{2}+{H}_{\epsilon }\left(t\right)/2,& & \end{array}$ (70)
thus
 $\begin{array}{ccc}{H}_{\epsilon }+Q\ge {H}_{\epsilon }/2-C{\epsilon }^{2},& & \end{array}$ (71)
and we can transform ( 69 ) in
 $\begin{array}{ccc}\frac{d}{dt}\left({H}_{\epsilon }+Q\right)\le C\left({H}_{\epsilon }+Q\right)+C{\epsilon }^{2}.& & \end{array}$ (72)
Gronwall’s lemma then yields ${H}_{\epsilon }\left(t\right)+Q\left(t\right)\le \left({H}_{\epsilon }\left(0\right)+Q\left(0\right)+Ct{\epsilon }^{2}\right)exp\left(Ct\right).$  Using again ( 70 ) we obtain
 $\begin{array}{ccc}{H}_{\epsilon }\left(t\right)\le C\left({H}_{\epsilon }\left(0\right)+{\epsilon }^{2}\right)exp\left(Ct\right),& & \end{array}$ (73)
which achieves the proof of Theorem  5.1 .
$\square$  Proof of identity ( 68 ):
We first notice that, for all $g\in {C}^{1}\left(\mathbb{R}×{\mathbb{T}}^{d}\right)$  , we have:
$\frac{d}{dt}\int \rho \left(t,x\right)g\left(t,x\right)dx=\int \int f\left(t,x,\xi \right)\left({\partial }_{t}g\left(t,x\right)+\xi \cdot \nabla g\left(t,x\right)\right)d\xi dx.$  We also use the conservation of energy defined by ( 65 ). Then we get
 $\begin{array}{ccc}\frac{d}{dt}{H}_{\epsilon }=& & \frac{d}{dt}\frac{1}{2}\int f\left(t,x,\xi \right)\left(|v{|}^{2}-2\xi \cdot v\right)dxd\xi \end{array}$
 $\begin{array}{ccc}=& & \int f\left(t,x,\xi \right)\left({\partial }_{t}v\cdot v-{\partial }_{t}v\cdot \xi \right)-\frac{1}{2}\int {\nabla }_{x}\cdot \left(f\left(t,x,\xi \right)\xi \right)\left(|v{|}^{2}-2\xi \cdot v\right)\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\int {\nabla }_{\xi }\cdot \left(\frac{1}{\epsilon }\nabla \stackrel{~}{\phi }f\left(t,x,\xi \right)\right)\left(|v{|}^{2}-2\xi \cdot v\right).\end{array}$
Integrating by part, we get
 $\begin{array}{ccc}\frac{d}{dt}{H}_{\epsilon }=& & \int f\left(t,x,\xi \right)\left({\partial }_{t}v\cdot v-{\partial }_{t}v\cdot \xi \right)+\int f\left(t,x,\xi \right)\xi \nabla v\left(v-\xi \right)\end{array}$
 $\begin{array}{ccc}& & +\int f\left(t,x,\xi \right)\frac{1}{\epsilon }\nabla \stackrel{~}{\phi }\cdot v.\end{array}$
The first two terms can be rewritten as
 $\begin{array}{ccc}& & \int f\left(t,x,\xi \right)\left({\partial }_{t}v\cdot v-{\partial }_{t}v\cdot \xi \right)+\int f\left(t,x,\xi \right)\xi \nabla v\left(v-\xi \right)\end{array}$
 $\begin{array}{ccc}=& & -\int f\left(t,x,\xi \right)\left(v-\xi \right)\nabla v\left(v-\xi \right)+\int f\left(t,x,\xi \right){\partial }_{t}v\cdot \left(v-\xi \right)\end{array}$
 $\begin{array}{ccc}& & +\int f\left(t,x,\xi \right)v\nabla v\left(v-\xi \right)\end{array}$
 $\begin{array}{ccc}=& & -\int f\left(t,x,\xi \right)\left(v-\xi \right)\nabla v\left(v-\xi \right)+\int f\left(t,x,\xi \right)\left(v-\xi \right)\cdot \left({\partial }_{t}v+v\cdot \nabla v\right),\end{array}$
and finally using equation ( 66 ) we conclude.
$\square$

5.2 Comparison with the Euler-Poisson system

Here we show that, as mentioned in the introduction, the Euler-Poisson ( $EP$  ) system is a more accurate approximation to the Vlasov Monge-Ampere system than the Euler equations, as $\epsilon$  goes to zero.
The $EP$  system Let us recall that the (pressureless) Euler-Poisson system describes the motion of a continuum of electrons on a neutralizing background of ions through electrostatic interaction. Let $\overline{v}$  and $\overline{\rho }$  be the velocity and density of electrons. Let $\overline{\phi }$  be the (rescaled) electric potential. Under proper scaling, these functions of $x\in {\mathbb{R}}^{d}$  and $t>0$  satisfy the Euler-Poisson system:
 $\begin{array}{ccc}& & {\partial }_{t}\overline{v}+\overline{v}\cdot \nabla \overline{v}=-\frac{1}{\epsilon }\nabla \overline{\phi }\end{array}$ (74)
 $\begin{array}{ccc}& & {\partial }_{t}\overline{\rho }+\nabla \cdot \left(\overline{\rho }\overline{v}\right)=0\end{array}$ (75)
 $\begin{array}{ccc}& & 1-\epsilon \Delta \overline{\phi }=\overline{\rho }.\end{array}$ (76)
The so-called ’quasi-neutral’ limit $\epsilon \to 0$  of similar systems has been studied for example in [? and [?, and convergence results have been established using pseudo-differentials energy estimates. For well-prepared initial data, solutions of $EP$  are expected to behave as solutions of Euler incompressible equations.
This fact is proved by the second author in his PhD thesis ([?, Chap 2), see also [?. We give here the complete result that we will use herafter. We will denote by ${\overline{v}}^{\epsilon }$  (resp. ${f}^{\epsilon }\right)$  the solutions of the $EP$  (resp. $VM{A}_{p}$  ) system with parameter $\epsilon$  .
Theorem 5.3 Let $v$  be a solution of ( 66 , 67 ) on $\left[0,T\right]×{\mathbb{T}}^{d}$  , with initial data ${v}_{0}$  , and satisfying $v\in {L}^{\infty }\left(\left[0,T\right],{H}^{s}\left({\mathbb{T}}^{d}\right)\right)$  for some $s\ge {s}_{0}\left(d\right)$  .
There for some ${s}^{\prime }>0$  , ${s}^{\prime }  , if $\left({\overline{v}}_{0}^{\epsilon },{\overline{\rho }}_{0}^{\epsilon }\right)$  is such that the sequences
 $\begin{array}{ccc}\frac{{\overline{v}}_{0}^{\epsilon }-{v}_{0}}{\epsilon },\frac{{\overline{\rho }}_{0}^{\epsilon }-1}{{\epsilon }^{2}}& & \end{array}$
are bounded in ${H}^{{s}^{\prime }}\left({\mathbb{T}}^{d}\right)$  , then there exists ${T}_{\epsilon }>0$  with ${liminf}_{\epsilon \to 0}{T}_{\epsilon }\ge T$  and a sequence $\left({\overline{v}}^{\epsilon },{\overline{\rho }}^{\epsilon }\right)$  of solutions to the $EP$  system on $\left[0,{T}_{\epsilon }\left[$  with initial data $\left({\overline{v}}_{0}^{\epsilon },{\overline{\rho }}_{0}^{\epsilon }\right)$  , belonging to ${L}^{\infty }\left(\left[0,{T}_{\epsilon }\right],{H}^{{s}^{\prime }}\left({\mathbb{T}}^{d}\right)\right)$  . Moreover, for $\epsilon$  small enough, the sequences
 $\begin{array}{ccc}\frac{{\overline{v}}^{\epsilon }-v}{\epsilon },\frac{{\overline{\rho }}^{\epsilon }-1}{{\epsilon }^{2}}& & \end{array}$
are bounded in ${L}^{\infty }\left(\left[0,T\right],{H}^{{s}^{\prime }}\left({\mathbb{T}}^{d}\right)\right)$  . Finally, ${s}^{\prime }$  goes to $+\infty$  as $s$  goes to $+\infty$  .
Assumptions Here we consider $v$  a solution to $E$  ( 66 ,  67 ) with initial data ${v}_{0}$  , a sequence ${f}^{\epsilon }$  of solutions of $VM{A}_{p}$  ( 63 , 64 ) with initial data ${f}_{0}^{\epsilon }$  , and a sequence $\left({\overline{v}}^{\epsilon },{\overline{\rho }}^{\epsilon }\right)$  solutions of $EP$  ( 74 ,  75 ,  76 ) with initial data $\left({\overline{v}}_{0}^{\epsilon },{\overline{\rho }}_{0}^{\epsilon }\right)$  . We still define ${H}_{\epsilon }$  as in Theorem  5.1 :
 $\begin{array}{ccc}{H}_{\epsilon }\left(t\right)=\frac{1}{2}\int {f}^{\epsilon }\left(t,x,\xi \right)|\xi -v\left(t,x\right){|}^{2}dxd\xi +\frac{1}{2}\int |\nabla {\phi }^{\epsilon }{|}^{2}dx.& & \end{array}$
We introduce the following assumptions:
• H0 $v$  solution of $E$  satisfies, for some ${C}_{0}>0$  , $\parallel v{\parallel }_{{L}^{\infty }\left(\left[0,T\right]{H}^{s}\left({\mathbb{T}}^{d}\right)\right)}\le {C}_{0}$  , and $s$  is large enough so that ${s}^{\prime }$  in Theorem  5.3 satisfies $s\ge {s}^{\prime }>\left[\frac{d}{2}\right]+2$  .
• H1 The sequence $\left({\overline{v}}_{0}^{\epsilon },{\overline{\rho }}_{0}^{\epsilon }\right)$  of initial data of $EP$  is such that, for some ${C}_{1}>0$  ,  $\begin{array}{ccc}{sup}_{\epsilon >0}\left\{\frac{1}{\epsilon }\parallel {\overline{v}}_{0}^{\epsilon }-v{\parallel }_{{H}^{{s}^{\prime }}\left({\mathbb{T}}^{d}\right)},\frac{1}{{\epsilon }^{2}}\parallel {\overline{\rho }}^{\epsilon }-1{\parallel }_{{H}^{{s}^{\prime }}\left({\mathbb{T}}^{d}\right)}\right\}\le {C}_{1}.& & \end{array}$
• H2 The sequence ${f}_{0}^{\epsilon }$  satisfies ${H}_{\epsilon }\left(0\right)\le {C}_{2}{\epsilon }^{2}$  for some ${C}_{2}>0$  .
H0, H1, H2 imply that
• 1. There exists ${\stackrel{~}{C}}_{0}$  such that  $\begin{array}{ccc}\parallel v{\parallel }_{{L}^{\infty }\left(\left[0,T\right],{W}^{2,\infty }\left({\mathbb{T}}^{d}\right)\right)}\le {\stackrel{~}{C}}_{0}.& & \end{array}$ (77)
• 2. From Theorem  5.1 , there exists ${\stackrel{~}{C}}_{1}$  such that  $\begin{array}{ccc}{H}_{\epsilon }\left(t\right)\le {\stackrel{~}{C}}_{1}{\epsilon }^{2}\text{for}t\in \left[0,T\right].& & \end{array}$ (78)
• 3. From Theorem  5.4 and Sobolev imbeddings, there exists ${\stackrel{~}{C}}_{2}$  such that  $\begin{array}{ccc}{sup}_{\epsilon <{\epsilon }_{0}}\left\{{‖\frac{{\overline{v}}^{\epsilon }-v}{\epsilon },\frac{{\overline{\rho }}^{\epsilon }-1}{{\epsilon }^{2}}‖}_{{L}^{\infty }\left(\left[0,T\right],{W}^{2,\infty }\left({\mathbb{T}}^{d}\right)\right)}\right\}\le {\stackrel{~}{C}}_{2}.& & \end{array}$ (79)
We are now ready to prove the following result:
Theorem 5.4 Let ${f}_{0}^{\epsilon },{\overline{v}}_{0}^{\epsilon },{\overline{\rho }}_{0}^{\epsilon },v,T$  be as above, satisfying assumptions H0, H1, H2. Define
 $\begin{array}{ccc}{G}_{\epsilon }\left(t\right)=\frac{1}{2}\int {f}^{\epsilon }\left(t,x,\xi \right)|\xi -{\overline{v}}^{\epsilon }\left(x\right){|}^{2}dxd\xi +\frac{1}{2}\int |\nabla {\phi }^{\epsilon }-\nabla {\overline{\phi }}^{\epsilon }{|}^{2}dx.& & \end{array}$
Then there exists $C>0$  such that
 $\begin{array}{ccc}{G}_{\epsilon }\left(t\right)\le Cexp\left(Ct\right)\left({G}_{\epsilon }\left(0\right)+{\epsilon }^{3}\right),\forall t\in \left[0,T\right]& & \end{array}$
where C depends on ${s}^{\prime },{C}_{0},{C}_{1},{C}_{2},T$  .
Remark: the theorem shows that the distance between solutions of the $\left(EP\right)$  system and the $VM{A}_{p}$  system measured with ${G}_{\epsilon }$  is like $O\left({\epsilon }^{3}\right)$  whereas Theorem  5.1 showed that the distance between the solution of the Euler equation and the $VM{A}_{p}$  system was like $O\left({\epsilon }^{2}\right)$  . Note also that ${G}_{\epsilon }$  and ${H}_{\epsilon }$  can both be interpreted as the square of a distance.
Proof of Theorem  5.4 : For notational simplicity, we drop most $\epsilon$  ’s. Proceeding as in ( 68 ) and noticing that:
 $\begin{array}{ccc}\frac{d}{dt}{\int }_{{\mathbb{T}}^{d}}|\nabla \overline{\phi }{|}^{2}=\frac{1}{\epsilon }{\int }_{{\mathbb{T}}^{d}}\overline{\rho }\overline{v}\cdot \nabla \overline{\phi }& & \end{array}$
we obtain the following identity:
 $\begin{array}{ccc}\frac{d}{dt}{G}_{\epsilon }=& & -\int f\left(t,x,\xi \right)\left(\xi -\overline{v}\right)\nabla \overline{v}\left(\xi -\overline{v}\right)\end{array}$
 $\begin{array}{ccc}& & +\int f\left(t,x,\xi \right)\frac{1}{\epsilon }\overline{v}\cdot \nabla \stackrel{~}{\phi }-\int f\left(t,x,\xi \right)\frac{1}{\epsilon }\overline{v}\cdot \nabla \overline{\phi }\end{array}$
 $\begin{array}{ccc}& & +\int f\left(t,x,\xi \right)\frac{1}{\epsilon }\xi \cdot \nabla \overline{\phi }+\int \frac{1}{\epsilon }\overline{\rho }\overline{v}\cdot \nabla \overline{\phi }\end{array}$
 $\begin{array}{ccc}& & -\frac{d}{dt}\int \nabla \overline{\phi }\cdot \nabla \phi .\end{array}$ (80)
Then we notice
 $\begin{array}{ccc}\int f\left(t,x,\xi \right)\frac{1}{\epsilon }\xi \cdot \nabla \overline{\phi }=\frac{d}{dt}\left(\int \frac{1}{\epsilon }\rho \overline{\phi }\right)-\frac{1}{\epsilon }\int \rho {\partial }_{t}\overline{\phi }.& & \end{array}$
Next, we have the following lemma:
Lemma 5.5 Define for any $\theta \in {C}^{2}\left({\mathbb{T}}^{d}\right)$
 $\begin{array}{ccc}& & <\nabla \theta >\left(x\right)={\int }_{0}^{1}\nabla \theta \left(x+s\epsilon \nabla \phi \left(x\right)\right)ds,\end{array}$
 $\begin{array}{ccc}& & <{\nabla }^{2}\theta >\left(x\right)={\int }_{0}^{1}\left(1-s\right){\nabla }^{2}\theta \left(x+s\epsilon \nabla \phi \left(x\right)\right)ds.\end{array}$
Then
 $\begin{array}{ccc}\int \rho \theta & =& \int \theta +\epsilon \int <\nabla \theta >\cdot \nabla \phi \end{array}$
 $\begin{array}{ccc}& =& \int \theta +\epsilon \int \nabla \theta \cdot \nabla \phi +{\epsilon }^{2}\int <{\nabla }^{2}\theta >\nabla \phi \nabla \phi .\end{array}$
Proof: The proof just uses the Taylor expansion and the identity $\int \rho \theta =\int \theta \left(x+\epsilon \nabla \phi \right)$  .
$\square$  Using Lemma  5.5 , we get
 $\begin{array}{ccc}& & \frac{1}{\epsilon }\int \rho {\partial }_{t}\overline{\phi }\end{array}$
 $\begin{array}{ccc}=& & \frac{1}{\epsilon }\int {\partial }_{t}\overline{\phi }+\int {\partial }_{t}\nabla \overline{\phi }\cdot \nabla \phi +\epsilon \int <{\partial }_{t}{\nabla }^{2}\overline{\phi }>\nabla \phi \nabla \phi .\end{array}$
We claim that, under our assumptions, we have $\parallel {\partial }_{t}{\nabla }^{2}\overline{\phi }{\parallel }_{{L}^{\infty }\left(\left[0,{T}^{\prime }\right]×{\mathbb{T}}^{d}\right)}\le C.$  Proof: from ( 75 ), we have
 $\begin{array}{ccc}{\partial }_{t}\overline{\rho }=-\overline{\rho }\nabla \cdot \overline{v}-\overline{v}\cdot \nabla \overline{\rho }.& & \end{array}$
Using ( 79 ), we obtain that $\parallel {\partial }_{t}\overline{\rho }{\parallel }_{{H}^{{s}^{\prime }-1}}\le C\epsilon$  . Since ${H}^{{s}^{\prime }}\left({\mathbb{T}}^{d}\right)$  is continuously embedded in ${W}^{2,\infty }\left({\mathbb{T}}^{d}\right)$  , ${H}^{{s}^{\prime }-1}\left({\mathbb{T}}^{d}\right)$  is continuously embedded in ${L}^{\infty }\left({\mathbb{T}}^{d}\right)$  .
Then, using ( 76 ) and classical elliptic regularity, we have $\epsilon \parallel {\partial }_{t}{\nabla }^{2}\overline{\phi }{\parallel }_{{H}^{{s}^{\prime }-1}}\le C\parallel {\partial }_{t}\overline{\rho }{\parallel }_{{H}^{{s}^{\prime }-1}},$  and the desired result follows.
$\square$  This implies, using ( 78 ), that
 $\begin{array}{ccc}|\epsilon \int <{\partial }_{t}{\nabla }^{2}\overline{\phi }>\nabla \phi \nabla \phi |\le C{\epsilon }^{3}.& & \end{array}$
Next,
 $\begin{array}{ccc}& & \int {\partial }_{t}\nabla \overline{\phi }\cdot \nabla \phi =-\int {\partial }_{t}\Delta \overline{\phi }\phi \end{array}$
 $\begin{array}{ccc}& & =\frac{1}{\epsilon }\int {\partial }_{t}\overline{\rho }\phi =\frac{1}{\epsilon }\int \overline{\rho }\overline{v}\cdot \nabla \phi .\end{array}$
Using again Lemma  5.5 , we get
 $\begin{array}{ccc}& & \frac{d}{dt}\int \nabla \overline{\phi }\cdot \nabla \phi \end{array}$
 $\begin{array}{ccc}=& & \frac{1}{\epsilon }\frac{d}{dt}\left(\int \rho \overline{\phi }-{\epsilon }^{2}\int <{\nabla }^{2}\overline{\phi }>\nabla \phi \nabla \phi \right)\end{array}$
and for the same reasons we have $\parallel {\nabla }^{2}\overline{\phi }{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)\right)}\le C\epsilon$  . This yields $Q\left(t\right)=\epsilon \int <{\nabla }^{2}\overline{\phi }>\nabla \phi \nabla \phi =O\left({\epsilon }^{4}\right).$  Moreover, it does not cost to set $\int \overline{\phi }\equiv 0$  and deduce
 $\begin{array}{ccc}\int f\left(t,x,\xi \right)\frac{1}{\epsilon }\xi \cdot \nabla \overline{\phi }-\frac{d}{dt}\int \nabla \overline{\phi }\cdot \nabla \phi =-\frac{1}{\epsilon }\int \overline{\rho }\overline{v}\cdot \nabla \phi +O\left({\epsilon }^{3}\right)+\frac{d}{dt}Q.& & \end{array}$
Thus the remaining terms are
 $\begin{array}{ccc}R=\frac{1}{\epsilon }\int \left[\rho \nabla \stackrel{~}{\phi }-\rho \nabla \overline{\phi }+\overline{\rho }\nabla \overline{\phi }-\overline{\rho }\nabla \phi \right]\cdot \overline{v}.& & \end{array}$
Calculations that we postpone to the end of the proof show that
 $\begin{array}{ccc}R\le & & \int \left(\nabla \phi -\nabla \overline{\phi }\right)\nabla \overline{v}\left(\nabla \phi -\nabla \overline{\phi }\right)+C\int |\nabla \phi -\nabla \overline{\phi }{|}^{2}\end{array}$
 $\begin{array}{ccc}& & -\frac{1}{2}\int \nabla \cdot \overline{v}\left(|\nabla \overline{\phi }{|}^{2}-2\nabla \phi \cdot \nabla \overline{\phi }\right)+C{\epsilon }^{3}.\end{array}$ (81)
with $C$  depending on $\parallel {\nabla }^{2}\overline{v}{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)}$  and ${\epsilon }^{-1}\parallel {\nabla }^{3}\overline{\phi }{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)}$  , therefore uniformly bounded thanks to ( 79 ). Finally we obtain
 $\begin{array}{ccc}\frac{d}{dt}{G}_{\epsilon }\le & & -\int f\left(t,x,\xi \right)\left(\xi -\overline{v}\right)\nabla \overline{v}\left(\xi -\overline{v}\right)+\left(\nabla \phi -\nabla \overline{\phi }\right)\nabla \overline{v}\left(\nabla \phi -\nabla \overline{\phi }\right)\end{array}$
 $\begin{array}{ccc}& & -\frac{1}{2}\int \left(\nabla \cdot \overline{v}\right)\left(|\nabla \overline{\phi }{|}^{2}-2\nabla \overline{\phi }\cdot \nabla \phi \right)+C\int |\nabla \phi -\nabla \overline{\phi }{|}^{2}\end{array}$
 $\begin{array}{ccc}& & +C{\epsilon }^{3}+\frac{d}{dt}Q\end{array}$
with $|Q\left(t\right)|\le C{\epsilon }^{4}$  for $t\in \left[0,T\right]$  . From ( 79 ) we have $\parallel \nabla \cdot v{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)}\le C\epsilon$  and $\parallel \nabla \overline{\phi }{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)}\le C\epsilon$  , whereas ( 78 ) yields $\int {|\nabla \phi |}^{2}\le C{\epsilon }^{2}$  . Note that we also have
 $\begin{array}{ccc}& & -\int f\left(t,x,\xi \right)\left(\xi -\overline{v}\right)\nabla \overline{v}\left(\xi -\overline{v}\right)+\left(\nabla \phi -\nabla \overline{\phi }\right)\nabla \overline{v}\left(\nabla \phi -\nabla \overline{\phi }\right)\end{array}$
 $\begin{array}{ccc}& & +C\int |\nabla \phi -\nabla \overline{\phi }{|}^{2}\le C{G}_{\epsilon }.\end{array}$
We conclude that
 $\begin{array}{ccc}\frac{d}{dt}\left({G}_{\epsilon }-Q\right)\le C\left(\left({G}_{\epsilon }-Q\right)+{\epsilon }^{3}\right),& & \end{array}$
and the conclusion of Theorem  5.4 follows by Gronwall’s lemma.
$\square$  Proof of identity ( 81 ): Here we have to compute:
 $\begin{array}{ccc}R=\frac{1}{\epsilon }\int \overline{v}\left(x+\epsilon \nabla \phi \right)\cdot \nabla \phi -\left(\overline{v}\nabla \overline{\phi }\right)\left(x+\epsilon \nabla \phi \right)+\left(1-\epsilon \Delta \overline{\phi }\right)\left(\overline{v}\cdot \nabla \overline{\phi }-\overline{v}\cdot \nabla \phi \right)& & \end{array}$
Using Lemma  5.5 we have:
 $\begin{array}{ccc}R=& & \frac{1}{\epsilon }\int \overline{v}\cdot \nabla \phi -\overline{v}\cdot \nabla \overline{\phi }+\overline{v}\cdot \nabla \overline{\phi }-\overline{v}\cdot \nabla \phi \end{array}$
 $\begin{array}{ccc}& & +\int \nabla \overline{v}\cdot \nabla \phi \nabla \phi -\nabla \left(\overline{v}\nabla \overline{\phi }\right)\nabla \phi -\overline{v}\nabla \overline{\phi }\Delta \overline{\phi }+\overline{v}\nabla \phi \Delta \overline{\phi }\end{array}$
 $\begin{array}{ccc}& & +\int \left(<\nabla \overline{v}>-\nabla \overline{v}\right)\nabla \phi \nabla \phi -\epsilon <{\nabla }^{2}\left(\overline{v}\nabla \overline{\phi }\right)>\nabla \phi \nabla \phi .\end{array}$
We see that the first line cancels. Then we show that the last line is bounded by $C{\epsilon }^{3}$  .
This is obvious for the last term since from ( 77 ,  78 ) we have $\parallel \overline{v}{\parallel }_{{W}^{2,\infty }}\le C$  , and $\parallel \nabla \overline{\phi }{\parallel }_{{W}^{2,\infty }}\le C\epsilon$  .
Then for the first term we have the following lemma:
Lemma 5.6 We define $\Delta =\int \left(<\nabla \overline{v}>\left(x\right)-\nabla \overline{v}\left(x\right)\right)\nabla \phi \nabla \phi dx,$  then one has:
$|\Delta |\le C{\epsilon }^{10/3}+C\int |\nabla \phi -\nabla \overline{\phi }{|}^{2}dx.$
Proof: First we show that if $\Theta \left(R\right)={\int }_{\left\{|\nabla \phi |\ge R\right\}}{|\nabla \phi |}^{2}$  , $\Theta \left(R\right)\le C\int {|\nabla \phi -\nabla \overline{\phi }|}^{2}+\frac{C{\epsilon }^{4}}{{R}^{2}}.$  Proof: $\int {|\nabla \phi |}^{2}\le C{\epsilon }^{2}$  , implies that $\text{meas}\left\{|\nabla \phi |\ge R\right\}\le C\left(\frac{\epsilon }{R}{\right)}^{2}.$  Since $|\nabla \overline{\phi }\left(t,x\right)|\le \epsilon$  for $\left(t,x\right)\in \left[0,{T}^{\prime }x\right]×{\mathbb{T}}^{d}$
 $\begin{array}{ccc}& & \Theta \left(R\right)\le {\int }_{\left\{|\nabla \phi |\ge R\right\}}{|\nabla \overline{\phi }|}^{2}+{\int }_{\left\{|\nabla \phi |\ge R\right\}}{|\nabla \phi -\nabla \overline{\phi }|}^{2}\end{array}$
 $\begin{array}{ccc}& & \le \frac{C{\epsilon }^{4}}{{R}^{2}}+\int {|\nabla \phi -\nabla \overline{\phi }|}^{2}.\end{array}$
Then we have
 $\begin{array}{ccc}& & \Delta \le C\Theta \left(R\right)+{\int }_{|\nabla \phi |\le R}|<\nabla \overline{v}>\left(x\right)-\nabla \overline{v}\left(x\right)|\nabla \phi \nabla \phi \end{array}$
 $\begin{array}{ccc}\text{with}& & |<\nabla \overline{v}>\left(x\right)-\nabla \overline{v}\left(x\right)|\le C\epsilon |\nabla \phi |\end{array}$
 $\begin{array}{ccc}\text{thus}& & \Delta \le C\epsilon {\int }_{|\nabla \phi |\le R}{|\nabla \phi |}^{3}+C\Theta \left(R\right)\end{array}$
 $\begin{array}{ccc}& & \le C\left(\epsilon R\int {|\nabla \phi |}^{2}+\frac{{\epsilon }^{4}}{{R}^{2}}+\int {|\nabla \phi -\nabla \overline{\phi }|}^{2}\right)\end{array}$
 $\begin{array}{ccc}& & \le C\left({\epsilon }^{3}R++\frac{{\epsilon }^{4}}{{R}^{2}}+\int {|\nabla \phi -\nabla \overline{\phi }|}^{2}\right)\end{array}$
for all R, so for $R={\epsilon }^{\left(1/3\right)}$  one obtains:
 $\begin{array}{ccc}& & \Delta \le C{\epsilon }^{10/3}+C\int {|\nabla \phi -\nabla \overline{\phi }|}^{2}.\end{array}$
$\square$  Thus we have shown that $R=S+O\left({\epsilon }^{3}\right)$  , and $S={\Sigma }_{k=1}^{6}{T}_{k}$  where each ${T}_{k}$  is given by:
 $\begin{array}{ccc}& & {T}_{1}={\partial }_{j}{\overline{v}}_{i}{\partial }_{j}\phi {\partial }_{i}\phi \end{array}$
 $\begin{array}{ccc}& & {T}_{2}=-{\partial }_{j}{\overline{v}}_{i}{\partial }_{j}\phi {\partial }_{i}\overline{\phi }\end{array}$
 $\begin{array}{ccc}& & {T}_{3}=-{\overline{v}}_{i}{\partial }_{ij}\overline{\phi }{\partial }_{j}\phi \end{array}$
 $\begin{array}{ccc}& & {T}_{4}={\partial }_{j}{\overline{v}}_{i}{\partial }_{j}\overline{\phi }{\partial }_{i}\overline{\phi }\end{array}$
 $\begin{array}{ccc}& & {T}_{5}={\overline{v}}_{i}{\partial }_{ij}\overline{\phi }{\partial }_{j}\overline{\phi }\end{array}$
 $\begin{array}{ccc}& & {T}_{6}={\overline{v}}_{i}{\partial }_{jj}\overline{\phi }{\partial }_{i}\phi \end{array}$
 $\begin{array}{ccc}& & \end{array}$
where we have used Einstein’s convention for repeated indices. First we have
 $\begin{array}{ccc}{T}_{5}=-\frac{1}{2}\int \left(\nabla \cdot \overline{v}\right){|\nabla \overline{\phi }|}^{2}& & \end{array}$
 $\begin{array}{ccc}{T}_{1}+{T}_{2}+{T}_{4}=\int {\partial }_{j}{\overline{v}}_{i}\left({\partial }_{j}\phi -{\partial }_{j}\overline{\phi }\right)\left({\partial }_{i}\phi -{\partial }_{i}\overline{\phi }\right)+{T}_{7}& & \end{array}$
with ${T}_{7}=\int {\partial }_{j}{\overline{v}}_{i}{\partial }_{j}\overline{\phi }{\partial }_{i}\phi$  .
 $\begin{array}{ccc}{T}_{6}=-\int {\partial }_{i}{\overline{v}}_{i}{\partial }_{jj}\overline{\phi }\phi +{\overline{v}}_{i}{\partial }_{ijj}\overline{\phi }\phi & & \end{array}$
and
 $\begin{array}{ccc}-\int {\overline{v}}_{i}{\partial }_{ijj}\overline{\phi }\phi =\int {\partial }_{j}{\overline{v}}_{i}{\partial }_{ij}\overline{\phi }\phi +{\overline{v}}_{i}{\partial }_{ij}\overline{\phi }{\partial }_{j}\phi & & \end{array}$
thus
 $\begin{array}{ccc}{T}_{6}=\int -\left(\nabla \cdot \overline{v}\right)\Delta \overline{\phi }\phi +{T}_{8}-{T}_{3}& & \end{array}$
with ${T}_{8}=\int {\partial }_{j}{\overline{v}}_{i}{\partial }_{ij}\overline{\phi }\phi$  . Then
 $\begin{array}{ccc}{T}_{8}=& & -\int {\partial }_{j}{\overline{v}}_{i}{\partial }_{j}\overline{\phi }{\partial }_{i}\phi +{\partial }_{ij}{\overline{v}}_{i}{\partial }_{j}\overline{\phi }\phi \end{array}$
 $\begin{array}{ccc}& & =-{T}_{7}+\int \nabla \cdot \overline{v}\left(\Delta \overline{\phi }\phi +\nabla \overline{\phi }\nabla \phi \right)\end{array}$
and finally we obtain
 $\begin{array}{ccc}S\left(t\right)=& & \int \nabla \overline{v}\left(\nabla \overline{\phi }-\nabla \phi \right)\left(\nabla \overline{\phi }-\nabla \phi \right)-\frac{1}{2}\left(\nabla \cdot \overline{v}\right){|\nabla \overline{\phi }-\nabla \phi |}^{2}\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\int \left(\nabla \cdot \overline{v}\right){|\nabla \phi |}^{2}\end{array}$
and the identity ( 81 ) is proved.
$\square$

6 Appendix

Existence and uniqueness of solutions to second order ODE’s with BV field In this section we prove existence and a.e. uniqueness for ordinary differential equations of the form:
 $\begin{array}{ccc}\frac{d}{dt}\left(\begin{array}{c}X\\ V\end{array}\right)=\left(\begin{array}{c}V\\ E\left(t,X\right)\end{array}\right)& & \end{array}$ (82)
for $X\in {\mathbb{T}}^{d},Y\in {\mathbb{T}}^{d}$  , and where the field $E$  belongs to ${L}^{\infty }\left(\right]0,T\left[×{\mathbb{T}}^{d}\right)\cap {L}^{1}\left(\right]0,T\left[,BV\left({\mathbb{T}}^{d}\right)\right)$  . We work in the flat torus for simplicity, but our results are still valid in an open subset of ${\mathbb{R}}^{d}$  . This result is an adaptation of the proof of [? that uses the result of [? on renormalized solutions of transport equations.
Remark: After this proof was written, the authors learned of a result by L. Ambrosio ([?) that extends the results of [? to transport equations when the vector field is in $BV$  with bounded divergence.
Renormalized solutions for Vlasov equations with BV field Theorem 3.4 in [? adapted to the periodic case sates that if $f\in {L}^{\infty }\left(\right]0,T\left[×{\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)$  satisfies:
 $\begin{array}{ccc}\frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+{\nabla }_{\xi }\cdot \left(E\left(t,x\right)f\right)=0,& & \end{array}$ (83)
with $E\left(t,x\right)\in {L}^{1}\left(\right]0,T\left[×{\mathbb{T}}^{d}\right)\cap {L}^{1}\left(\right]0,T\left[,BV\left({\mathbb{T}}^{d}\right)\right)$  , then for all $g$  Lipschitz continuous we have
 $\begin{array}{ccc}\frac{\partial g\left(f\right)}{\partial t}+{\nabla }_{x}\cdot \left(\xi g\left(f\right)\right)+{\nabla }_{\xi }\cdot \left(E\left(t,x\right)g\left(f\right)\right)=0.& & \end{array}$
The property of renormalization implies that
• solutions to ( 83 ) with initial data in ${L}_{loc}^{\infty }\left({\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)$  belong to $C\left(\right]0,T\left[,{L}_{loc}^{p}\left({\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)\right)$  for any $1\le p<\infty$  ,
• solutions to ( 83 ) with prescribed initial data in ${L}^{\infty }\left({\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)$  are a.e.
unique,
• if ${E}_{n}$  converges to $E$  in ${L}^{1}\left(\right]0,T\left[×{\mathbb{T}}^{d}\right)$  then the solutions of ( 83 ) with ${E}_{n}$  instead of $E$  converge to the solution of ( 83 ).
We notice that equation ( 54 ) satisfies the assumptions of the Theorem, and thus will have the renormalization property. This renormalization property was used in [? to obtain a.e. uniqueness for solutions of the corresponding ODE’s. Indeed, the ODE
 $\begin{array}{ccc}& & {\partial }_{t}X\left(t,s,x\right)=b\left(t,X\right)\end{array}$
 $\begin{array}{ccc}& & X\left(s,s,x\right)=x\end{array}$
is associated to the transport equation:
 $\begin{array}{ccc}{\partial }_{t}u+b\left(t,x\right).\nabla u=0& & \end{array}$
whose solutions satisfy for all $\left(t,s\right)\in \right]0,T\left[$
 $\begin{array}{ccc}u\left(t,X\left(t,s,x\right)\right)=u\left(s,x\right).& & \end{array}$
We extend this consequence to the case of second order equations, with $BV$  acceleration field. To the kinetic equation
 $\begin{array}{ccc}{\partial }_{t}f+\xi \cdot {\nabla }_{x}f+E\left(t,x\right)\cdot {\nabla }_{\xi }f=0& & \end{array}$ (84)
we associate the second order ODE ( 82 ) which can rewritten as ${\partial }_{tt}X=E\left(t,X\right)$  . The result is then the following:
Theorem 6.1 Let $E\left(t,x\right)\in {L}^{\infty }\left(\right]0,T\left[×{\mathbb{T}}^{d}\right)\cap {L}^{1}\left(\right]0,T\left[,BV\left({\mathbb{T}}^{d}\right)\right)$  , then the ODE
 $\begin{array}{ccc}& & {\partial }_{tt}X\left(t,s,x,\xi \right)=E\left(t,X\right)\end{array}$ (85)
 $\begin{array}{ccc}& & \left(X\left(s,s,x,\xi \right),{\partial }_{t}X\left(s,s,x,\xi \right)\right)=\left(x,\xi \right)\end{array}$ (86)
Remark: Here almost everywhere must be understood for the Lebesgue measure of ${\mathbb{R}}^{6}$  .
Proof of Theorem  6.1 : We know that through equation ( 82 ) equation ( 85 ) can be considered as a first order differential equation. Let us first consider the case where $E$  is smooth. Note $Y\in {\mathbb{T}}^{d}×{\mathbb{R}}^{d}$  (resp. $y$  ) for $\left(X,V\right)$  (resp.
for $\left(x,\xi \right)$  ) and $B\in {\mathbb{R}}^{d}×{\mathbb{R}}^{d}$  for $\left(\xi ,E\right)$  . Then for all $s\in \right]0,T\left[$  , $Y$  solves:
 $\begin{array}{ccc}& & {\partial }_{t}Y\left(t,s,y\right)=B\left(t,Y\left(t,s,y\right)\right)\end{array}$ (87)
 $\begin{array}{ccc}& & Y\left(s,s,y\right)=y\end{array}$ (88)
Then for all $t,{t}_{1},{t}_{2},{t}_{3}\in \right]0,T\left[$  we have the following:
 $\begin{array}{ccc}& & Y\left({t}_{3},{t}_{2},Y\left({t}_{2},{t}_{1},y\right)\right)=Y\left({t}_{3},{t}_{1},y\right)\end{array}$
 $\begin{array}{ccc}& & Y\left(t,t,y\right)=y\end{array}$
 $\begin{array}{ccc}& & Y\left({t}_{1},{t}_{2},Y\left({t}_{2},{t}_{1},y\right)\right)=y.\end{array}$
Differentiating the last equation with respect to ${t}_{2}$  yields:
 $\begin{array}{ccc}& & {\partial }_{s}Y\left(t,s,y\right)+{\nabla }_{y}Y\left(t,s,y\right)\cdot B\left(s,y\right)=0\end{array}$ (89)
 $\begin{array}{ccc}& & Y\left(t,t,y\right)=y.\end{array}$ (90)
${Y}_{t}\left(s,y\right)=Y\left(t,s,y\right)$  thus solves a transport equation which is nothing but equation ( 84 ). Using Theorem  3.5 we know that for all $g:{\mathbb{R}}^{2d}\to \mathbb{R}$  Lipschitz continuous, $g\left(t,s,y\right)={g}_{0}\left(Y\left(t,s,y\right)\right)$  is the unique solution of
 $\begin{array}{ccc}& & {\partial }_{s}g\left(t,s,y\right)+{\nabla }_{y}g\left(t,s,y\right)\cdot B\left(s,y\right)=0\end{array}$ (91)
 $\begin{array}{ccc}& & g\left(t,t,y\right)={g}_{0}\left(y\right).\end{array}$ (92)
Now we show existence and uniqueness for solutions of ( 87 , 88 ). Let $t$  and $s$  be fixed. Let us consider a regularization ${E}_{n}$  of the the field $E$  and set ${B}_{n}=\left(\xi ,{E}_{n}\right)$  . We have
• $t\to {Y}_{1,n}\left(t,s,y\right)$  that satisfies ( 87 , 88 )
• $s\to {Y}_{2,n}\left(t,s,y\right)$  that satisfies ( 89 , 90 ).
From the stability Theorem 2.4 in [? we know that the whole sequence $t\to {Y}_{2,n}\left(t,s,.\right)$  converges in $C\left(\right]0,T\left[,{L}_{loc}^{p}\left({\mathbb{R}}^{d}×{\mathbb{T}}^{d}\right)\right)$  to $t\to {Y}_{2}\left(t,s,.\right)$  , the unique renormalized solution of ( 89 , 90 ). Thus for fixed $t$  the whole sequence ${Y}_{2,n}\left(t,s,.\right)$  converges strongly in ${L}_{loc}^{p}\left({\mathbb{R}}^{d}×{\mathbb{T}}^{d}\right)$  . Now since for every $n$  we have ${Y}_{1,n}\left(t,s,y\right)={Y}_{2,n}\left(t,s,y\right)$  the same property holds for ${Y}_{1,n}\left(s,t,.\right).$  Now we can pass to the limit in the term ${B}_{n}\left(t,{Y}_{1,n}\left(t,s,y\right)\right)$  . Indeed, by density of ${C}_{c}^{\infty }$  functions in ${L}^{1},$  if we have ${E}_{s}\in {C}_{c}^{\infty }$  approximating $E$  then
 $\begin{array}{ccc}& & \parallel B\left(t,{Y}_{n}\left(t,s,y\right)\right)-B\left(t,Y\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
 $\begin{array}{ccc}& \le & \parallel B\left(t,{Y}_{n}\left(t,s,y\right)\right)-{B}_{s}\left(t,{Y}_{n}\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
 $\begin{array}{ccc}& +& \parallel {B}_{s}\left(t,{Y}_{n}\left(t,s,y\right)\right)-{B}_{s}\left(t,Y\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
 $\begin{array}{ccc}& +& \parallel B\left(t,Y\left(t,s,y\right)\right)-{B}_{s}\left(t,Y\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
The second term goes to 0 because of the strong convergence of ${Y}_{n}$  , the first and the third go to 0 because $Y$  and ${Y}_{n}$  are measure preserving mappings, and so for example $\parallel B\left(t,Y\left(t,s,y\right)\right)-{B}_{s}\left(t,Y\left(t,s,y\right)\right){\parallel }_{{L}^{1}}=\parallel B\left(t,y\right)-{B}_{s}\left(t,y\right){\parallel }_{{L}^{1}}$  .
So finally we have
 $\begin{array}{ccc}& & \parallel {B}_{n}\left(t,{Y}_{n}\left(t,s,y\right)\right)-B\left(t,Y\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
 $\begin{array}{ccc}\le & & \parallel {B}_{n}\left(t,{Y}_{n}\left(t,s,y\right)\right)-B\left(t,{Y}_{n}\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
 $\begin{array}{ccc}& & +\parallel B\left(t,{Y}_{n}\left(t,s,y\right)\right)-B\left(t,Y\left(t,s,y\right)\right){\parallel }_{{L}^{1}}\end{array}$
that goes to 0 and we can pass to the limit in the equation ( 87 , 88 ) and the existence of a solution to ( 87 , 88 ) is proved.
To obtain uniqueness, we argue as in [?. Any function of the form ${g}_{0}\left(Y\left(t,s,y\right)\right)$  is a solution of ( 91 , 92 ), thus by uniqueness of the solution of the transport equation we obtain uniqueness of the ODE. $\square$  A remark on ODE’s of second order In this section, we want to solve the Cauchy problem for:
 $\begin{array}{ccc}& & {\partial }_{tt}X\left(t,x\right)=E\left(t,X\right)\end{array}$
 $\begin{array}{ccc}& & \left(X\left(0,x\right),{\partial }_{t}X\left(0,x\right)\right)=\left(x,v\left(x\right)\right)\end{array}$
with $E$  as above. We are thus interested in monokinetic initial data.
Theorem 6.2 for all ${v}^{0}\left(x\right)$  vector field on ${\mathbb{T}}^{d}$  , and for Lebesgue almost every $\delta v\in {\mathbb{R}}^{d}$  , there exists an a.e. unique solution to
 $\begin{array}{ccc}& & {\partial }_{tt}X\left(t,x\right)=E\left(t,X\left(t,x\right)\right)\end{array}$
 $\begin{array}{ccc}& & \left(X\left(0,x\right),{\partial }_{t}X\left(0,x\right)\right)=\left(x,{v}^{0}\left(x\right)+\delta v\right)\end{array}$
Proof: Let $g\left(x,\xi \right)$  be the indicator function of the set of those $\left(x;\xi \right)$  such that the trajectory coming from x is not well defined. We just have to prove that for a.e. $\delta v\in {\mathbb{R}}^{d}$  we have $\int g\left(x,{v}^{0}\left(x\right)+\delta v\right)dx=0,$  which is true because $\int g\left(x,{v}^{0}\left(x\right)+\xi \right)dxd\xi =\int g\left(x,\xi \right)dxd\xi =0.$

Stability

Using the fact that for ${E}_{n}$  converging to $E$  in ${L}^{1}$  with $E\in {L}^{1}\left(\right]0,T\left[,BV\left({\mathbb{T}}^{d}\right)\right)$  , we have ${X}_{n}\left(t,x,v\right)\to X\left(t,x,v\right)$  in $C\left(\left[0,T\right],{L}^{p}\right)$  , we have then, for all t, for almost every $\delta v,$  ${X}_{n}\left(t,x,{v}^{0}\left(x\right)+\delta v\right)\to X\left(t,x,{v}^{0}\left(x\right)+\delta v\right)$  in ${L}^{p}$  . Thus we have
Theorem 6.3 If ${E}_{n}$  converges to $E$  in ${L}^{1}$  let ${X}_{n}$  be solution of
 $\begin{array}{ccc}& & {\partial }_{tt}{X}_{n}\left(t,x\right)={E}_{n}\left(t,{X}_{n}\left(t,x\right)\right)\end{array}$
 $\begin{array}{ccc}& & \left({X}_{n}\left(0,x\right),{\partial }_{t}{X}_{n}\left(0,x\right)\right)=\left(x,{v}^{0}\left(x\right)+\delta v\right)\end{array}$
then for all t, for almost every $\delta v,$  ${X}_{n}$  converges in ${L}^{p}\left({\mathbb{R}}^{3}\right)-s$  to a solution (unique for almost every $\delta v$  ) of
 $\begin{array}{ccc}& & {\partial }_{tt}X\left(t,x\right)=E\left(t,X\right)\end{array}$
 $\begin{array}{ccc}& & \left(X\left(0,x\right),{\partial }_{t}X\left(0,x\right)\right)=\left(x,{v}^{0}\left(x\right)+\delta v\right).\end{array}$
Control of macroscopic density in kinetic equations We prove here Lemma  4.7 :
Lemma 6.4 Let $f\in {L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}×{\mathbb{R}}^{d}\right)$  satisfy
 $\begin{array}{ccc}& & \frac{\partial f}{\partial t}+{\nabla }_{x}\cdot \left(\xi f\right)+{\nabla }_{\xi }\cdot \left(E\left(t,x\right)f\right)=0\end{array}$ (93)
 $\begin{array}{ccc}& & f\left(0,.,.\right)={f}^{0}\end{array}$ (94)
with $E\in {L}^{1}\left(\left[0,T\right];BV\left({\mathbb{T}}^{d}\right)\right)$  and
 $\begin{array}{ccc}\parallel E{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{T}}^{d}\right)}\le F.& & \end{array}$ (95)
Let the initial condition ${f}_{0}$  be such that:
$a\left(x,\xi \right)\le f\left(0,x,\xi \right)\le b\left(x,\xi \right),$  with ${\rho }_{a}\left(x\right)=\int a\left(x,\xi \right)d\xi \ge m>0$  and ${\rho }_{b}\left(x\right)=\int b\left(x,\xi \right)d\xi \le M<\infty$  and $a,b$  satisfying
 $\begin{array}{ccc}|{\nabla }_{x,\xi }a,b|\le \frac{c}{1+|\xi {|}^{d+2}}.& & \end{array}$ (96)
Then there exists a constant $R>0$  such that
 $\begin{array}{ccc}\left({\rho }_{a}\left(x\right)-Rt\right)\le \rho \left(t,x\right)\le \left({\rho }_{b}\left(x\right)+Rt\right).& & \end{array}$
Proof: First suppose that the force field and the initial data are smooth.
For equation ( 93 , 94 ) we can exhibit characteristics $\left(x,\xi \right)\left(t;{t}_{0},{x}_{0},{\xi }_{0}\right),$  giving the evolution of the particles in the phase space. We have $f\left(t,x,\xi \right)=f\left({t}_{0},{x}_{0},{\xi }_{0}\right).$  Since the initial data is compactly supported and the force field is bounded in the ${L}^{\infty }$  norm, we have $|\xi -{\xi }_{0}|\le F|t-{t}_{0}|,$  $|x-{x}_{0}|\le \left(|{\xi }_{0}|+\frac{F}{2}|t-{t}_{0}|\right)|t-{t}_{0}|.$  If for $t=0$  we have $a\left(x,\xi \right)\le f\left(0,x,\xi \right)\le b\left(x,\xi \right)$  then
 $\begin{array}{ccc}& & \underline{A}\left(t,x,\xi \right)\le f\left(t,x,\xi \right)\le \overline{B}\left(t,x,\xi \right)\end{array}$
 $\begin{array}{ccc}& & \underline{A}\left(t,x,\xi \right)={inf}_{|{\sigma }_{1}|,|{\sigma }_{2}|\le 1}a\left(x+|t-{t}_{0}|\left(\xi +\frac{F}{2}|t-{t}_{0}|\right){\sigma }_{1},\xi +F|t-{t}_{0}|{\sigma }_{2}\right)\end{array}$
 $\begin{array}{ccc}& & \overline{B}\left(t,x,\xi \right)={sup}_{|{\sigma }_{1}|,|{\sigma }_{2}|\le 1}b\left(x+|t-{t}_{0}|\left(\xi +\frac{F}{2}|t-{t}_{0}|\right){\sigma }_{1},\xi +F|t-{t}_{0}|{\sigma }_{2}\right).\end{array}$
Using ( 96 ) and integrating in $\xi$  we find thus a constant $R=R\left(F,C,d\right)$  such that for $t-{t}_{0}\le 1$  we have:
 $\begin{array}{ccc}{\rho }_{a}\left(x\right)-R|t-{t}_{0}|\le \rho \left(t,x\right)\le {\rho }_{b}\left(x\right)+R|t-{t}_{0}|.& & \end{array}$
Next we need to show that the solution of the regularized equation converges to the solution we are studying: this result comes from the uniqueness of the solution to ( 93 , 94 ) which is a consequence of the renormalization property. Indeed since $E$  is bounded in $BV$  the system ( 93 , 94 ) admits a unique renormalized solution and the sequence of approximate solutions converge in $C\left(\left[0,T\right],{L}_{x,\xi }^{p}\right)$  for $1\le p<\infty$  thus the bounds obtained above are preserved.
$\square$  Regularity of the polar factorization on the flat torus Here we deduce from [?, [? and [?, [?, [? the Theorem  4.4 .
Theorem 6.5 If $\rho \in {C}^{\alpha }\left({\mathbb{T}}^{d}\right)$  with $0  is a probability measure on ${\mathbb{T}}^{d}$  then $\Psi =\Psi \left[\rho \right]$  (see Definition  4.3 ) is a classical solution of
 $\begin{array}{ccc}det{D}^{2}\Psi =\rho & & \end{array}$ (97)
and satisfies:
 $\begin{array}{ccc}& & \parallel \nabla \Psi \left(x\right)-x{\parallel }_{{L}^{\infty }}\le C\left(d\right)=\sqrt{d}/2\end{array}$ (98)
 $\begin{array}{ccc}& & \parallel {D}^{2}\Psi {\parallel }_{{C}^{\alpha }}\le K\left(m,M,\parallel \rho {\parallel }_{{C}^{\alpha }}\right)\end{array}$ (99)
Proof of Theorem  6.5 : Consider $\rho$  a ${\mathbb{Z}}^{d}$  periodic probability measure, satisfying
 $\begin{array}{ccc}0 (100)
and $\Phi \left[\rho \right]$  as in Definition  4.3 . First it is shown in [? that
 $\begin{array}{ccc}|\nabla \Phi \left[\rho \right]\left(x\right)-x|\le C\left(d\right).& & \end{array}$ (101)
It follows that the strict convexity argument of [? applies: indeed if $\Phi =\Phi \left[\rho \right]$  is not strictly convex its graph contains a line and this contradicts ( 101 ).
Moreover since $\Phi -|x{|}^{2}/2$  is globally Lipschitz and periodic there exists $N\left(d\right)$  such that $\parallel \Phi -|x{|}^{2}/2{\parallel }_{{L}^{\infty }}\le N\left(d\right)$  . It follows then that there exists $0  and $M\left(d\right)$  such that
 $\begin{array}{ccc}B\left(r\left(d\right)\right)\subset \left\{\Phi -\Phi \left(0\right)\le M\left(d\right)\right\}\subset B\left(R\left(d\right)\right)& & \end{array}$ (102)
It remains to show that our solution is a solution in the Aleksandrov sense of the Monge-Ampere equation
 $\begin{array}{ccc}m\le det{D}^{2}\Phi \le M.& & \end{array}$
This is a direct consequence with minor changes (to adapt to the periodic case) of Lemma 2 of [?. Then, normalizing $\Phi$  to $\stackrel{~}{\Phi }=\Phi -\Phi \left(0\right)-M\left(d\right)$  it follows that $\stackrel{~}{\Phi }$  is a solution of
 $\begin{array}{ccc}& & \rho \left(\nabla \stackrel{~}{\Phi }\right)det{D}^{2}\stackrel{~}{\Phi }=1\end{array}$
 $\begin{array}{ccc}& & \stackrel{~}{\Phi }=0\text{on}\partial \Omega \end{array}$
 $\begin{array}{ccc}& & B\left(r\left(d\right)\right)\subset \Omega \subset B\left(R\left(d\right)\right)\end{array}$
Thus the interior regularity results of [? apply uniformly to all $\Phi \left[\rho \right]$  with $\rho$  satisfying ( 100 ) and $\parallel \rho {\parallel }_{{C}^{\alpha }\left({\mathbb{T}}^{d}\right)}$  bounded and Theorem  4.4 follows.
$\square$  Acknowledgments The authors acknowledge the support of the European IHP network ”HYKE” HPRN-CT-2002-00282 and the LRC CEA-Cadarache/UNSA.