## Uniqueness of the solution to the Vlasov-Poisson system with bounded density

### November 27, 2006

Abstract
In this note, we show uniqueness of weak solutions to the Vlasov-Poisson system on the only condition that the macroscopic density $\rho$  defined by $\rho \left(t,x\right)={\int }_{{\mathbb{R}}^{d}}f\left(t,x,\xi \right)d\xi$  is bounded in ${L}^{\infty }$  . Our proof is based on optimal transportation.

1 Introduction

The Vlasov-Poisson system (herafter (VP)) describes the evolution of cloud of electrons or gravitational matter through the equations
 $\begin{array}{ccc}& & {\partial }_{t}f+\xi \cdot {\nabla }_{x}f+\nabla \Psi \cdot {\nabla }_{\xi }f=0,\end{array}$ (1)
 $\begin{array}{ccc}& & \Delta \Psi =\epsilon \rho ,\end{array}$ (2)
where $\rho \left(t,x\right)=\int f\left(t,x,\xi \right)d\xi$  , and $\epsilon >0$  in the electrostatic (repulsive) case, $\epsilon <0$  in the gravitational (attractive) case. Here $f\left(t,x,\xi \right)\ge 0$  , denotes the density of electrons (or matter) at time $t\in {\mathbb{R}}^{+}$  , position $x\in {\mathbb{R}}^{3}$  , velocity $\xi \in {\mathbb{R}}^{3}$  . Equation ( 2 )means
 $\begin{array}{ccc}\Psi \left(t,x\right)=-\epsilon {\int }_{{\mathbb{R}}^{3}}\rho \left(t,y\right)\frac{1}{4\pi |x-y|}dy.& & \end{array}$ (3)
We denote by $\mathcal{ℳ}\left({\mathbb{R}}^{6}\right)$  (resp. ${\mathcal{ℳ}}^{+}\left({\mathbb{R}}^{6}\right)$  ) the set of bounded (resp. bounded and positive) measures on ${\mathbb{R}}^{6}$  . Given an initial datum ${f}^{0}\in {\mathcal{ℳ}}^{+}\left({\mathbb{R}}^{6}\right)$  , we look for solutions to ( 1 ,  2 ) such that
 $\begin{array}{ccc}f{|}_{t=0}={f}^{0}.& & \end{array}$ (4)
For $T>0$  , we will call $f$  a solution to ( 1 ,  2 ,  4 ) in ${\mathcal{D}}^{\prime }\left(\left[0,T\left[×{\mathbb{R}}^{6}\right)$  , if
• - $f\in C\left(\left[0,T\left[,{\mathcal{ℳ}}^{+}\left({\mathbb{R}}^{6}\right)-w*\right)$  ,
• - $\forall \phi \in {C}_{c}^{\infty }\left(\left[0,T\left[×{R}^{6}\right),$   $\begin{array}{ccc}\int f\left({\partial }_{t}\phi +\xi \cdot {\nabla }_{x}\phi +{\nabla }_{x}\Psi \cdot {\nabla }_{\xi }\phi \right)=-{\int }_{{\mathbb{R}}^{6}}{f}^{0}\phi {|}_{t=0},& & \end{array}$ (5)
• - $\Psi$  solves ( 3 ).
We will not discuss the conditions needed on $\Psi ,f$  to give sense to the product $f{\nabla }_{x}\Psi$  or to the singular integral ( 3 ), since we will only consider the case where $\rho \in {L}^{\infty }$  . In this case, ${\nabla }_{x}\Psi$  will be continuous, and the product $f{\nabla }_{x}\Psi$  will be well defined for $f$  a bounded measure.
Our result is the following:
Theorem 1.1. Given ${f}^{0}$  in ${\mathcal{ℳ}}^{+}\left({\mathbb{R}}^{6}\right)$  , given $T>0$  , there exists at most one weak solution to ( 1 ,  2 ,  4 ) in ${\mathcal{D}}^{\prime }\left(\left[0,T\left[×{\mathbb{R}}^{6}\right)$  such that
 $\begin{array}{ccc}{sup}_{t\in \left[0,T\left[}\parallel \rho {\parallel }_{{L}^{\infty }\left({\mathbb{R}}^{d}\right)}<+\infty .& & \end{array}$ (6)
Remark 1. To establish the existence of a solution to (VP) satisfying the bound ( 6 ) requires much more assumptions on the initial datum than what we need here ! This question is treated in [?.
Remark 2. Note that we do not ask for any bound on the moments of $f$  , and also that we do not ask the energy to be finite.
A sufficient condition for uniqueness had been given by Lions and Perthame in [?, relying on Lipschitz bounds on the initial data ${f}^{0}$  , but they expected a uniqueness result under the weaker assumption of bounded density. The Lipschitz condition had indeed later been relaxed by Robert in [? down to $f\in {L}^{\infty }$  compactly supported in $x$  and $\xi$  for $t\in \left[0,T\right]$  . Here we relax the bound on the support of $f$  , and we do not ask either $f$  to be bounded in ${L}^{\infty }$  . We only need a ${L}^{\infty }\left(\left[0,T\left[×{\mathbb{R}}^{3}\right)$  bound on $\rho \left(t,x\right)$  . Hence our result applies also to monokinetic solutions of ( 1 ,  2 ). In that that case, we have $f\left(t,x,\xi \right)=\rho \left(t,x\right)\delta \left(\xi -v\left(t,x\right)\right)$  for some vector field $v$  , and this gives formally a solution to the Euler-Poisson system
 $\begin{array}{ccc}& & {\partial }_{t}\rho +\nabla \cdot \left(\rho v\right)=0,\end{array}$ (7)
 $\begin{array}{ccc}& & {\partial }_{t}\left(\rho v\right)+\nabla \cdot \left(\rho v\otimes v\right)=\rho \nabla \Psi ,\end{array}$ (8)
 $\begin{array}{ccc}& & \Delta \Psi =\epsilon \rho .\end{array}$ (9)
Our proof will rely on optimal transportation, and the next section is devoted to recall some facts concerning this subject. A complete reference on this topic is [?. The technique we will use adapts to many similar problems, where a transport equation and an elliptic equation are coupled. The velocity field is the gradient of a potential satisfying an elliptic equation whose right hand side depends smoothly on the density. This has been observed in [? in the case of the 2-d incompressible Euler equations and the semi-geostrophic equations. It is interesting to notice that our technique gives a new proof of the uniqueness part in Youdovich’s Theorem, while the technique used by Robert in [? was and adaptation of Youdovich’s original proof (see [?).

2 Preliminary results on optimal transportation and Wasserstein distances

Definition 2.1. Let ${\rho }_{1},{\rho }_{2}$  be two probability measures on ${\mathbb{R}}^{d}$  .
We define the Wasserstein distance between ${\rho }_{1}$  and ${\rho }_{2}$  , that we denote ${W}_{2}\left({\rho }_{1},{\rho }_{2}\right)$  , by
 $\begin{array}{ccc}{W}_{2}\left({\rho }_{1},{\rho }_{2}\right)={\left({inf}_{\pi }{\int }_{{\mathbb{R}}^{d}×{\mathbb{R}}^{d}}\pi \left(x,y\right)|x-y{|}^{2}\right)}^{\frac{1}{2}},& & \end{array}$
where the infimum runs over probability measures $\pi$  on ${\mathbb{R}}^{d}×{\mathbb{R}}^{d}$  with marginals ${\rho }_{1}$  and ${\rho }_{2}$  .
Then we gather several results of optimal tranpsortations in the following theorem. These results can be found in Benamou & Brenier [?, McCann [?, Gangbo & McCann [?.
Theorem 2.2 (Benamou, Brenier, Gangbo, McCann). Let ${\rho }_{1}$  , ${\rho }_{2}$  be two probability measures on ${\mathbb{R}}^{d}$  , such that ${\rho }_{1},{\rho }_{2}$  are absolutely continuous with respect to the Lebesgue measure. Assume that ${W}_{2}\left({\rho }_{1},{\rho }_{2}\right)<+\infty$  .
• 1. There exists a pair $\left({\rho }_{\theta },{v}_{\theta }\right),\theta \in \left[1,2\right]$  such that ${\rho }_{\theta }\in C\left(\left[1,2\right];\mathcal{ℳ}\left({\mathbb{R}}^{d}\right)-w*\right)$  and ${v}_{\theta }$  is a $d{\rho }_{\theta }$  measurable vector field, that satisfies  $\begin{array}{ccc}& & {\partial }_{\theta }{\rho }_{\theta }+\nabla \cdot \left({\rho }_{\theta }{v}_{\theta }\right)=0,\end{array}$
 $\begin{array}{ccc}& & \rho {|}_{\theta =1}={\rho }_{1},\rho {|}_{\theta =2}={\rho }_{2},\end{array}$
 $\begin{array}{ccc}& & {\int }_{{\mathbb{R}}^{d}}|{v}_{\theta }\left(x\right){|}^{2}d{\rho }_{\theta }\left(x\right)\equiv {W}_{2}^{2}\left({\rho }_{1},{\rho }_{2}\right)\forall \theta \in \left[1,2\right].\end{array}$
• 2. For this ’optimal’ path we have also, when both ${\rho }_{1}$  and ${\rho }_{2}$  have densities in ${L}^{\infty }$  with respect to the Lebesgue measure,  $\begin{array}{ccc}\forall \theta \in \left[1,2\right],\parallel {\rho }_{\theta }{\parallel }_{{L}^{\infty }}\le max\left\{\parallel {\rho }_{1}{\parallel }_{{L}^{\infty }},\parallel {\rho }_{2}{\parallel }_{{L}^{\infty }}\right\}.& & \end{array}$
Remark. The path $\left\{{\rho }_{\theta },\theta \in \left[1,2\right]\right\}$  is the geodesic linking ${\rho }_{1}$  to ${\rho }_{2}$  , with respect to the Wasserstein metric (see [?, and also [? where this notion was introduced).

3 Proof of Theorem  1.1

From now on, we assume for simplicity that ${\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)=1$  , $\epsilon =1$  , and the reader can check that this choice does not play any role in the proof. In particular, the result of the previous section adapt with minor changes to the case of two positive measures of equal total mass.

3.1 ${H}^{-1}$  estimates along geodesics

In this section, we show that the ${H}^{-1}\left({\mathbb{R}}^{3}\right)$  norm of the difference ${\rho }_{1}-{\rho }_{2}$  can be controlled by the Wasserstein distance between ${\rho }_{1}$  and ${\rho }_{2}$  , provided both measures have densities in ${L}^{\infty }$  with respect to the Lebesgue measure.
This is the crucial estimate at the core of our result.
Proposition 3.1. Let ${\rho }_{1},{\rho }_{2}$  be two probability measures on ${\mathbb{R}}^{d}$  with ${L}^{\infty }$  densities with respect to the Lebesgue measure . Let ${\Psi }_{i},i=1,2$  solve
 $\begin{array}{ccc}\Delta {\Psi }_{i}={\rho }_{i}& & \end{array}$
in the sense of ( 3 ). Then $\parallel \nabla {\Psi }_{1}-\nabla {\Psi }_{2}{\parallel }_{{L}^{2}\left({\mathbb{R}}^{d}\right)}\le {\left[max\left\{\parallel {\rho }_{1}{\parallel }_{{L}^{\infty }},\parallel {\rho }_{2}{\parallel }_{{L}^{\infty }}\right\}\right]}^{\frac{1}{2}}{W}_{2}\left({\rho }_{1},{\rho }_{2}\right),$  where ${W}_{2}\left({\rho }_{1},{\rho }_{2}\right)$  is the Wasserstein distance between ${\rho }_{1}$  and ${\rho }_{2}$  given in Definition  2.1 .
Proof of Proposition  3.1 . We suppose that ${W}_{2}\left({\rho }_{1},{\rho }_{2}\right)<+\infty$  otherwise there is nothing to prove. For $\theta \in \left[1,2\right]$  , we take ${\rho }_{\theta }$  that interpolates between ${\rho }_{1}$  and ${\rho }_{2}$  as in Theorem  2.2 , and ${v}_{\theta }$  the corresponding velocity field.
If we consider, for every $\theta \in \left[1,2\right]$  , ${\Psi }_{\theta }$  solution of
 $\begin{array}{ccc}\Delta {\Psi }_{\theta }={\rho }_{\theta },& & \end{array}$ (10)
then ${\Psi }_{\theta }$  interpolates between ${\Psi }_{1}$  and ${\Psi }_{2}$  . If we differentiate ( 10 ) with respect to $\theta$  , we obtain
 $\begin{array}{ccc}\Delta {\partial }_{\theta }{\Psi }_{\theta }& =& {\partial }_{\theta }{\rho }_{\theta }\end{array}$
 $\begin{array}{ccc}& =& -\nabla \cdot \left({\rho }_{\theta }{v}_{\theta }\right).\end{array}$
Note that since ${\rho }_{1},{\rho }_{2}$  are bounded, so will be ${\rho }_{\theta }$  from Theorem  2.2 , hence using that ${W}_{2}^{2}\left({\rho }_{1},{\rho }_{2}\right)=\int {\rho }_{\theta }|{v}_{\theta }{|}^{2}<+\infty$  , we have ${\rho }_{\theta }{v}_{\theta }$  bounded in ${L}^{2}$  . This implies ${\partial }_{\theta }{\Psi }_{\theta }\in {W}^{1,2}$  .
We integrate the above equation against ${\partial }_{\theta }{\Psi }_{\theta }$  to obtain
 $\begin{array}{ccc}\int |\nabla {\partial }_{\theta }{\Psi }_{\theta }{|}^{2}=-\int {\rho }_{\theta }{v}_{\theta }\cdot \nabla {\partial }_{\theta }{\Psi }_{\theta },& & \end{array}$
and this yields
 $\begin{array}{ccc}\parallel {\partial }_{\theta }\nabla {\Psi }_{\theta }{\parallel }_{{L}^{2}\left({\mathbb{R}}^{d}\right)}& \le & \parallel {\rho }_{\theta }{v}_{\theta }{\parallel }_{{L}^{2}\left({\mathbb{R}}^{d}\right)}\end{array}$
 $\begin{array}{ccc}& \le & \parallel {\rho }_{\theta }{\parallel }_{{L}^{\infty }\left({\mathbb{R}}^{d}\right)}^{1/2}{W}_{2}\left({\rho }_{1},{\rho }_{2}\right).\end{array}$
Using then the second point of Theorem  2.2 , and integrating over $\theta \in \left[1,2\right]$  , this proves Proposition  3.1 .
$\square$

3.2 Lagrangian formulation of the Vlasov-Poisson system

Given a solution of (VP) with bounded density $\rho$  on $\left[0,T\left[$  , we consider for $t\in \left[0,T\left[$  the characteristics of equation ( 1 ), that solve the ODE
 $\begin{array}{ccc}& & \stackrel{˙}{X}=\xi ,\end{array}$ (11)
 $\begin{array}{ccc}& & \stackrel{˙}{\xi }=\nabla \Psi \left(t,X\right).\end{array}$ (12)
Since we assume an ${L}^{\infty }\left(\left[0,T\left[×{\mathbb{R}}^{3}\right)$  bound on the density $\rho$  , the field $\nabla \Psi$  classically satisfies a log-Lispchitz condition:
 $\begin{array}{ccc}& & \forall t\in \left[0,T\left[,\forall \left(x,y\right)\in {\mathbb{R}}^{3}×{\mathbb{R}}^{3},|x-y|\le \frac{1}{2},\end{array}$
 $\begin{array}{ccc}& & |\nabla \Psi \left(t,x\right)-\nabla \Psi \left(t,y\right)|\le C|x-y|log\frac{1}{|x-y|},\end{array}$ (13)
where $C$  depends on $\parallel \rho {\parallel }_{{L}^{\infty }\left(\left[0,T\left[×{\mathbb{R}}^{3}\right)}$  . This condition is enough to define a Holder continuous flow
 $\begin{array}{ccc}\Xi \left(t,x,v\right)=\left(X,\xi \right)\left(t,x,v\right)& & \end{array}$
for the ODE ( 11 , 12 ), where $\left(X,\xi \right)$  is the pair (velocity, position) at time $t$  of the trajectory having (velocity, position) equal to $\left(x,v\right)$  at time $0$  .
Then we use the following Theorem, proved in [?:
Theorem 3.2. Let $u\left(t,x\right)$  be a vector field on ${\mathbb{R}}^{d}$  . Consider the ODE
 $\begin{array}{ccc}\stackrel{˙}{\gamma }\left(t\right)=u\left(t,\gamma \left(t\right)\right),& & \end{array}$
and the PDE
 $\begin{array}{ccc}{\partial }_{t}\mu \left(t,x\right)+\nabla \cdot \left(\mu \left(t,x\right)u\left(t,x\right)\right)=0.& & \end{array}$
Let $B\subset {\mathbb{R}}^{d}$  be a Borel set. The following are equivalent:
• (a) For all $x$  in $B$  , there exists a unique solution to the ODE starting at $x$  .
• (b) Non negative measure-valued solutions to the PDE with initial data ${\mu }^{0}$  concentrated in $B$  are unique.
From this result, we deduce the following corollary:
Corollary 3.3. The potential $\Psi$  being held fixed, and satisfying $\Delta \Psi \in {L}^{\infty }\left(\left[0,t\right]×{\mathbb{R}}^{3}\right)$  , for any ${f}^{0}\in {\mathcal{ℳ}}^{+}\left({\mathbb{R}}^{6}\right)$  there exists a unique weak solution to ( 1 ) (i.e. in the sense of ( 5 )) with initial datum ${f}^{0}$  which is given by
 $\begin{array}{ccc}f\left(t\right)=\Xi \left(t,\cdot ,\cdot {\right)}_{#}{f}^{0},& & \end{array}$ (14)
where $\Xi =\left(X,\xi \right)$  solves ( 11 ,  12 ). Note also that we will have
 $\begin{array}{ccc}\rho \left(t\right)=X\left(t,\cdot ,\cdot {\right)}_{#}{f}^{0}.& & \end{array}$ (15)
We remind the reader that the measure $f\left(t\right)=\Xi \left(t,\cdot ,\cdot {\right)}_{#}{f}^{0}$  is defined by $f\left(t\right)\left(B\right)={f}^{0}\left({\Xi }^{-1}\left(t\right)\left(B\right)\right)$  for all Borel subsets $B$  of ${\mathbb{R}}^{6}$  .
Remark. This corollary does not solve the uiqueness problem, but says only that if we suppress the coupling between $\Psi$  and $\rho$  , there is a unique weak measure-valued solution to the transport equation ( 1 ), that we can represent with the help of characteristics.

3.3 Final estimate

Given an initial distribution ${f}^{0}\left(x,\xi \right)\in {\mathcal{ℳ}}^{+}\left({\mathbb{R}}^{6}\right)$  with ${\int }_{{\mathbb{R}}^{6}}{f}^{0}=1$  , we take two solutions $\left({f}_{1},{f}_{2}\right)$  to (VP) with bounded density and initial datum ${f}^{0}$  . We have $\Delta \Psi ={\rho }_{i},i=1,2$  in the sense of ( 3 ). We then consider the associated characteristics ${\Xi }_{1}$  and ${\Xi }_{2}$  , where for $i=1,2$  , ${\Xi }_{i}=\left({X}_{i},{\xi }_{i}\right)\left(t,x,\xi \right)$  and ${X}_{i},{\xi }_{i}$  solve ( 11 ,  12 ) with force field $\nabla {\Psi }_{i}$  . Note that we will have ${f}_{i}\left(t\right)={\Xi }_{i}\left(t{\right)}_{#}{f}^{0},i=1,2$  . We then consider $Q\left(t\right)=\frac{1}{2}{\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|{\Xi }_{1}\left(t,x,\xi \right)-{\Xi }_{2}\left(t,x,\xi \right)|}^{2}.$  Remark. Notice that $\left({\Xi }_{1}\left(t\right),{\Xi }_{2}\left(t\right){\right)}_{#}{f}^{0}$  is a probability measure on ${\mathbb{R}}^{3}×{\mathbb{R}}^{3}$  , with marginals ${f}_{1}\left(t\right)$  and ${f}_{2}\left(t\right)$  , hence by Definition  2.1 ,
 $\begin{array}{ccc}Q\left(t\right)\ge \frac{1}{2}{W}_{2}^{2}\left({f}_{1}\left(t\right),{f}_{2}\left(t\right)\right).& & \end{array}$
(We will repeat this argument in Lemma  3.4 .) This implies in particular that $\left(Q=0\right)⟺\left({f}_{1}={f}_{2}\right)$  . Our proof will rely on an estimate on the Wasserstein distance between ${f}_{1}$  and ${f}_{2}$  , while the proof of [? was obtained by estimating the ${H}^{-1}$  norm of ${f}_{1}-{f}_{2}$  .
Of course $Q\left(0\right)=0$  , and
 $\begin{array}{ccc}\frac{d}{dt}Q\left(t\right)& =& {\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)\left({\Xi }_{1}\left(t,x,\xi \right)-{\Xi }_{2}\left(t,x,\xi \right)\right)\cdot {\partial }_{t}\left({\Xi }_{1}\left(t,x,\xi \right)-{\Xi }_{2}\left(t,x,\xi \right)\right)\end{array}$
 $\begin{array}{ccc}& =& {\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)\left[\left({X}_{1}\left(t,x,\xi \right)-{X}_{2}\left(t,x,\xi \right)\right)\cdot \left({\xi }_{1}\left(t,x,\xi \right)-{\xi }_{2}\left(t,x,\xi \right)\right)\right]\end{array}$
 $\begin{array}{ccc}& +& {\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)\left[\left({\xi }_{1}\left(t,x,\xi \right)-{\xi }_{2}\left(t,x,\xi \right)\right)\cdot \end{array}$
 $\begin{array}{ccc}& & \left(\nabla {\Psi }_{1}\left(t,{X}_{1}\left(t,x,\xi \right)\right)-\nabla {\Psi }_{2}\left(t,{X}_{2}\left(t,x,\xi \right)\right)\right)\right].\end{array}$
The second line is bounded by $Q\left(t\right)$  , and using Cauchy-Schwartz inequality, the third line is bounded by
 $\begin{array}{ccc}& & {Q}^{1/2}\left(t\right){\left({\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|\nabla {\Psi }_{1}\left(t,{X}_{1}\left(t,x,\xi \right)\right)-\nabla {\Psi }_{2}\left(t,{X}_{2}\left(t,x,\xi \right)\right)|}^{2}\right)}^{1/2}\end{array}$
 $\begin{array}{ccc}& \le & {Q}^{1/2}\left(t\right)\left({\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|\nabla {\Psi }_{2}\left(t,{X}_{1}\left(t,x,\xi \right)\right)-\nabla {\Psi }_{2}\left(t,{X}_{2}\left(t,x,\xi \right)\right)|}^{2}\end{array}$
 $\begin{array}{ccc}& & +{{\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|\nabla {\Psi }_{2}\left(t,{X}_{1}\left(t,x,\xi \right)\right)-\nabla {\Psi }_{1}\left(t,{X}_{1}\left(t,x,\xi \right)\right)|}^{2}\right)}^{1/2}\end{array}$
 $\begin{array}{ccc}& =& {Q}^{1/2}\left(t\right){\left({T}_{1}\left(t\right)+{T}_{2}\left(t\right)\right)}^{1/2}.\end{array}$
Hence we have
 $\begin{array}{ccc}\frac{d}{dt}Q\left(t\right)\le Q\left(t\right)+{Q}^{1/2}\left(t\right){\left({T}_{1}\left(t\right)+{T}_{2}\left(t\right)\right)}^{1/2},& & \end{array}$ (16)
and we will now estimate ${T}_{2}$  and then ${T}_{1}$  .
For ${T}_{2}$  we have, using ( 15 ) and Proposition  3.1 ,
 $\begin{array}{ccc}{T}_{2}\left(t\right)& =& {\int }_{{\mathbb{R}}^{3}}{\rho }_{1}\left(t,x\right){|\nabla {\Psi }_{1}\left(t,x\right)-\nabla {\Psi }_{2}\left(t,x\right)|}^{2}\end{array}$
 $\begin{array}{ccc}& \le & C{W}_{2}^{2}\left({\rho }_{1}\left(t\right),{\rho }_{2}\left(t\right)\right),\end{array}$
where $C$  depends on the ${L}^{\infty }$  norms of ${\rho }_{1},{\rho }_{2}$  . Moreover, from the very definition of the Wasserstein distance given in Definition  2.1 , we have the elementary lemma:
Lemma 3.4. Given ${X}_{1},{X}_{2},{\rho }_{1},{\rho }_{2}$  as above, we have
 $\begin{array}{ccc}{W}_{2}\left({\rho }_{1}\left(t\right),{\rho }_{2}\left(t\right)\right)\le {\left({\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|{X}_{1}\left(t,x,\xi \right)-{X}_{2}\left(t,x,\xi \right)|}^{2}\right)}^{\frac{1}{2}}.& & \end{array}$
Proof. The proof follows immediately from Definition  2.1 by noticing that $\pi =\left({X}_{1}\left(t\right),{X}_{2}\left(t\right){\right)}_{#}{f}^{0}$  is a probability measure on ${\mathbb{R}}^{3}×{\mathbb{R}}^{3}$  with marginals ${\rho }_{1}\left(t\right)$  and ${\rho }_{2}\left(t\right)$  .
$\square$  Hence
 $\begin{array}{ccc}{W}_{2}^{2}\left({\rho }_{1},{\rho }_{2}\right)& \le & {\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|{X}_{2}\left(t,x,\xi \right)-{X}_{1}\left(t,x,\xi \right)|}^{2}\end{array}$
 $\begin{array}{ccc}& \le & {\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|{\Xi }_{1}\left(t,x,\xi \right)-{\Xi }_{2}\left(t,x,\xi \right)|}^{2}=2Q\left(t\right),\end{array}$
and we conclude that ${T}_{2}\left(t\right)\le CQ\left(t\right)$  for some $C$  depending on the ${L}^{\infty }$  bounds on ${\rho }_{1},{\rho }_{2}$  .
Now, we evaluate ${T}_{1}$  by standard arguments, using the log-Lipschitz regularity of $\nabla {\Psi }_{2}$  : note first that since ${\rho }_{i},i=1,2$  are bounded in ${L}^{\infty }$  , $\nabla {\Psi }_{i},i=1,2$  are also bounded in ${L}^{\infty }$  , hence for any $C>0$  , we can take $T$  small enough such that $\parallel {\Xi }_{i}-{\Xi }_{2}{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{R}}^{6}\right)}\le C$  . Thus we have, for some other $C$  depending on $\parallel {\rho }_{i}{\parallel }_{{L}^{\infty }},i=1,2$  , and as long as $\parallel {\Xi }_{1}-{\Xi }_{2}{\parallel }_{{L}^{\infty }}\le \frac{1}{2}$  ,
 $\begin{array}{ccc}{T}_{1}& =& {\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right){|\nabla {\Psi }_{2}\left(t,{X}_{1}\left(t,x,\xi \right)\right)-\nabla {\Psi }_{2}\left(t,{X}_{2}\left(t,x,\xi \right)\right)|}^{2}\end{array}$
 $\begin{array}{ccc}& \le & C{\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)\left(|{X}_{1}-{X}_{2}{|}^{2}{log}^{2}\frac{1}{|{X}_{1}-{X}_{2}|}\right)\left(t,x,\xi \right)\end{array}$
 $\begin{array}{ccc}& =& \frac{C}{4}{\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)\left(|{X}_{1}-{X}_{2}{|}^{2}{log}^{2}\left(|{X}_{1}-{X}_{2}{|}^{2}\right)\right)\left(t,x,\xi \right).\end{array}$
Then we use that $x↦x{log}^{2}x$  is concave for $0\le x\le 1/e$  , and we can assume (taking $T$  small enough) that $\parallel {\Xi }_{1}-{\Xi }_{2}{\parallel }_{{L}^{\infty }\left(\left[0,T\right]×{\mathbb{R}}^{6}\right)}\le 1/e$  , therefore by Jensen’s inequality we have
 $\begin{array}{ccc}& & {T}_{1}\left(t\right)\end{array}$
 $\begin{array}{ccc}\le & & \frac{C}{4}\left[{\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)|{X}_{1}-{X}_{2}{|}^{2}\left(t,x,\xi \right)\right]{log}^{2}\left[{\int }_{{\mathbb{R}}^{6}}{f}^{0}\left(x,\xi \right)|{X}_{1}-{X}_{2}{|}^{2}\left(t,x,\xi \right)\right]\end{array}$
 $\begin{array}{ccc}\le & & \frac{C}{2}Q\left(t\right){log}^{2}\left(2Q\left(t\right)\right).\end{array}$
Combining all these bounds in ( 16 ), we obtain that
 $\begin{array}{ccc}\frac{d}{dt}Q\left(t\right)\le CQ\left(t\right)\left(1+log\frac{1}{Q\left(t\right)}\right),& & \end{array}$
and we conclude by standard arguments that if $Q\left(0\right)=0$  , $Q\equiv 0$  on $\left[0,T\left[$  .
This achieves the proof of Theorem  1.1 .
$\square$
G. Loeper EPFL, SB, IMA 10015 Lausanne e-mail: gregoire.loeper@epfl.ch