in general relativity , kinetic theory has been used relatively sparsely to model phenomenological matter in comparison to fluid models , although interest has increased in recent years . from a mathematical point of view
there are fundamental advantages to using a kinetic description . in non - curved spacetimes kinetic theory has been studied intensively as a mathematical subject during several decades , and it has also played an important role from an engineering point of view . the main purpose of this review paper is to discuss mathematical results for the einstein - vlasov system .
however , in the first part of this introduction , we review kinetic theory in non - curved spacetimes and focus on the special - relativistic case , although some results in the non - relativistic case will also be mentioned .
the reason that we focus on the relativistic case is not only that it is more closely related to the main theme in this review , but also that the literature on relativistic kinetic theory is very sparse in comparison to the non - relativistic case , in particular concerning the relativistic and non - relativistic boltzmann equation .
we believe that a good understanding of kinetic theory in non - curved spacetimes is fundamental to good comprehension of kinetic theory in general relativity . moreover , it is often the case that mathematical methods used to treat the einstein - vlasov system are carried over from methods developed in the special relativistic or non - relativistic case . the purpose of kinetic theory is to model the time evolution of a collection of particles .
for instance , the particles are atoms and molecules in a neutral gas or electrons and ions in a plasma . in astrophysics the particles are stars , galaxies or even clusters of galaxies .
a characteristic feature of kinetic theory is that its models are statistical and the particle systems are described by density functions f = f(t , x , p ) , which represent the density of particles with given spacetime position ( t , x )     and momentum p. a density function contains a wealth of information , and macroscopic quantities are easily calculated from this function . in a fluid model
the quantities that describe the system do not depend on the momentum p but only on the spacetime point ( t , x ) .
a choice of model is usually made with regard to the physical properties of interest for the system or with regard to numerical considerations .
it should be mentioned that a too naive fluid model may give rise to shell - crossing singularities , which are unphysical . in a kinetic description
the time evolution of the system is determined by the interactions between the particles , which depend on the physical situation .
for instance , the driving mechanism for the time evolution of a neutral gas is the collision between particles ( the boltzmann equation ) . for a plasma
the interaction is through the electromagnetic field produced by the charges ( the vlasov - maxwell system ) , and in astrophysics the interaction is gravitational ( the vlasov - poisson system and the einstein - vlasov system ) .
of course , combinations of interaction processes are also considered but in many situations one of them is strongly dominating and the weaker processes are neglected . consider a collection of neutral particles in minkowski spacetime .
let the signature of the metric be (  , + , + , + ) .
in this section we assume that all the particles have rest mass m = 1 , and we normalize the speed of light c to one .
we point out that in section 2 on the einstein - vlasov system , the dependence on the rest mass and the speed of light will be included in the formulation of the system
. the four - momentum of a particle is denoted by p , a = 0 , 1 , 2 , 3 .
since all particles have equal rest mass , the four - momentum for each particle is restricted to the mass shell , ppa = m = 1 .
thus , by denoting the three - momentum by p   , p may be written p = ( p , p ) , where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${p^0 } = \sqrt { 1 + \vert p{\vert ^2}}$\end{document } is the energy of a particle with three - momentum p , and p is the usual euclidean length of p. the relativistic velocity of a particle with momentum p is denoted by p and is given by 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\hat p = { p \over { \sqrt { 1 + \vert p\vert ^{2}}}}.$$\end{document } note that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\vert \hat p\vert < 1 = c$\end{document}. the relativistic boltzmann equation models the spacetime behavior of the one - particle distribution function f = f(t , x , p ) , and it has the form 2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\left({{\partial _ t } + { p \over { { p^0 } } } \cdot { \nabla _ x } } \right)f = q(f , f),$$\end{document } where the relativistic collision operator q(f , g ) is defined by 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(f , g ) = \int\nolimits_{{{\mathbb r}^2 } } { \int\nolimits_{{{\mathbb s}^2 } } { k(p , q,\omega)[f(p + a(p , q,\omega)\omega)g(q - a(p , q,\omega)\omega ) - f(p)g(q)]d\omega \;dp.}}$$\end{document } note that g = f in equation ( 2 ) . here
d is the element of surface area on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${{\rm{{\mathbb s}}}^2}$\end{document } and k(p , q ,  ) is the scattering kernel , which depends on the differential cross - section in the interaction process .
we refer to ,   and   for examples of differential cross - sections in the relativistic case .
if two particles , with momentum p and q respectively , collide elastically with scattering angle \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\omega \in { { \rm{{\mathbb s}}}^2}$\end{document } , their momenta will change , i.e. , p  p and q  q. the relation between p , q and p , q is given by 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p^{\prime } = p + a(p , q,\omega)\omega , \quad q^{\prime } = q - a(p , q,\omega)\omega,$$\end{document } where 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(p , q,\omega ) = { { 2({p^0 } + { q^0}){p^0}{q^0}(\omega \cdot ( \hat q - \hat p ) ) } \over { { { ( { p^0 } + { q^0})}^2 } - { { ( \omega \cdot ( p + q))}^2}}}.$$\end{document } this relation is a consequence of four - momentum conservation , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p^a } + { q^a } = { p^{a^{\prime } } } + { q^{a{\prime}}},$$\end{document } or equivalently 6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p^0 } + { q^0 } = { p^{0{\prime } } } + { q^{0{\prime}}},$$\end{document } 
 7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p + q = { p^{\prime } } + { q^{\prime}}.$$\end{document } these are the conservation equations for relativistic particle dynamics . in the classical case the corresponding conservation equations read 8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\vert p\vert ^{2 } + \vert q\vert ^{2 } = \vert p^{\prime}\vert ^{2 } + \vert q^{\prime}\vert ^{2},$$\end{document } 
 9\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p + q = { p^{\prime } } + { q^{\prime}}.$$\end{document } the function ( p , q ,  ) gives the distance between p and p ( q and q ) in momentum space , and the analogue function in the non - relativistic , newtonian , classical case has the form 10\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${a_{{\rm{cl}}}}(p , q,\omega ) = \omega \cdot ( q - p).$$\end{document } by inserting ac1 in place of a in equation ( 3 ) we obtain the classical boltzmann collision operator ( disregarding the scattering kernel , which is also different ) .
we point out that there are other representations of the collision operator ( 3 ) , cf . . in   and  
classical solutions to the relativistic boltzmann equations are studied as c   , and it is proven that the limit as c   of these solutions satisfies the classical boltzmann equation .
the former work is more general since general initial data is considered , whereas the latter is concerned with data near vacuum .
the latter result is stronger in the sense that the limit , as c   , is shown to be uniform in time .
the main result concerning the existence of solutions to the classical boltzmann equation is a theorem by diperna and lions   that proves existence , but not uniqueness , of renormalized solutions .
an analogous result holds in the relativistic case , as was shown by dudyski and ekiel - jeewska , cf . also .
regarding classical solutions , illner and shinbrot   have shown global existence of solutions to the non - relativistic boltzmann equation for initial data close to vacuum .
glassey showed global existence for data near vacuum in the relativistic case in a technical work .
he only requires decay and integrability conditions on the differential cross - section , although these are not fully satisfactory from a physics point of view . by imposing more restrictive cut - off assumptions on the differential cross - section ,
strain   gives a different proof , which is more related to the proof in the non - relativistic case   than   is .
for the homogeneous relativistic boltzmann equation , global existence for small initial data has been shown in   under the assumption of a bounded differential cross - section .
for initial data close to equilibrium , global existence of classical solutions has been proven by glassey and strauss   using assumptions on the differential cross - section , which fall into the regime  hard potentials  , whereas strain   has shown existence in the case of soft potentials .
in the non - relativistic case , we refer to [ 189 , 172 , 119 ] for analogous results .
the collision operator q(f , g ) may be written in an obvious way as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(f , g ) = { q^ + } ( f , g ) - { q^ -}(f , g),$$\end{document } where q and q are called the gain and loss term , respectively . if the loss term is deleted the gain - term - only boltzmann equation is obtained .
it is interesting to note that the methods of proof for the small data results mentioned above concentrate on gain - term - only equations , and once that is solved it is easy to include the loss term . in   it is shown that the gain - term - only classical and relativistic boltzmann equations blow up for initial data not restricted to a small neighborhood of trivial data .
thus , if a global existence proof of classical solutions for unrestricted data will be given , it will necessarily use the full collision operator .
the gain term has a nice regularizing property in the momentum variable . in   it
is proven that given f  l( ) and g  l( ) with f , g  0 , then 11\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\vert \vert { q^ + } ( f , g)\vert \vert _ { { { h^1}({\mathbb r}_p^3 ) } } \leq c\vert \vert f\vert \vert _ { { { l^2}({\mathbb r}_p^3)}}\vert \vert g\vert \vert _ { { { l^1}({\mathbb r}_p^3)}},$$\end{document } under some technical requirements on the scattering kernel . here h is the usual sobolev space .
the proof relies on the theory of fourier integral operators and on the method of stationary phase , and requires a careful analysis of the collision geometry , which is very different in the relativistic case .
simplified proofs in the classical and relativistic case are given in [ 193 , 194 ] .
more precisely , lions used the regularizing theorem to prove that solutions to the classical boltzmann equation , with periodic boundary conditions , converge in l to a global maxwellian , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$m = { e^{- \alpha \vert p\vert ^{2 } + \beta \cdot p + \gamma}}\quad { \rm{with}}\;\alpha , \gamma \in r,\quad \alpha > 0,\quad \beta \in { { \mathbb r}^3},$$\end{document } as time goes to infinity .
it should be pointed out that the convergence takes place through a sequence of times tending to infinity and it is not known whether the limit is unique or depends on the sequence . in the relativistic situation , the analogous question of convergence to a relativistic maxwellian , or a jttner equilibrium solution , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$j = { e^{- \alpha \sqrt { 1 + \vert p\vert ^{2 } } + \beta \cdot p + \gamma}},\quad \alpha , \beta , \;{\rm{and}}\;\gamma \;{\rm{as}}\;{\rm{above}},\;{\rm{with}}\;\alpha >
\vert \beta \vert,$$\end{document } was studied by glassey and strauss [ 87 , 88 ] . in the periodic case , they proved convergence in a variety of function spaces for initial data close to a jttner solution .
having obtained the regularizing theorem for the relativistic gain term , it is a straightforward task to follow the method of lions and prove convergence to a global jttner solution for arbitrary initial data ( satisfying the natural bounds of finite energy and entropy ) , which are periodic in the space variables , cf . .
we also mention that in the non - relativistic case desvillettes and villani   have studied the convergence rate to equilibrium in detail .
a similar study in the relativistic case has not yet been achieved . for more information on the relativistic boltzmann equation on minkowski space
we refer to [ 54 , 68 , 181 , 79 ] and in the non - relativistic case we refer to [ 190 , 79 , 53 ] .
let us consider a collision - less plasma , which is a collection of particles for which collisions are relatively rare and the interaction is through their charges . for simplicity
we assume that the plasma consists of one type of particle , although the results below hold for plasmas with several particle species .
the particle rest mass and the particle charge are normalized to one . in the kinetic framework ,
the most general set of equations for modeling a collision - less plasma is the relativistic vlasov - maxwell system : 12\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ t}f + \hat v \cdot { \nabla _ x}f + ( e(t , x ) + \hat v \times b(t , x ) ) \cdot { \nabla _ v}f = 0$$\end{document } 
 13\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _
\quad \nabla \cdot e = \rho,$$\end{document } 
 14\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _
t}b = - c\nabla \times e,\quad \quad \nabla \cdot b = 0.$$\end{document } the notation follows the one already introduced with the exception that the momenta are now denoted by v instead of p. this has become a standard notation in this field .
e and b are the electric and magnetic fields , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\hat \upsilon}$\end{document } is the relativistic velocity , 15\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\hat v = { v \over { \sqrt { 1 + \vert v\vert ^{2}/{c^2}}}},$$\end{document } where is the speed of light .
the charge density  and current j are given by 16\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho = \int\nolimits_{{{\mathbb r}^3 } } { fdv,\quad j = } \int\nolimits_{{{\mathbb r}^3 } } { \hat vfdv.}$$\end{document } equation ( 12 ) is the relativistic vlasov equation and equations ( 13 , 14 ) are the maxwell equations .
a special case in three dimensions is obtained by considering spherically - symmetric initial data .
for such data it can be shown that the solution will also be spherically symmetric , and that the magnetic field has to be constant .
the maxwell equation   e = tb then implies that the electric field is the gradient of a potential . hence , in the spherically - symmetric case the relativistic vlasov - maxwell system takes the form 17\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ t}f + \hat v \cdot { \nabla _
x}f + \beta e(t , x ) \cdot { \nabla _ v}f = 0,$$\end{document } 
 18\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$e = \nabla \phi , \quad \nabla \phi = \rho.$$\end{document } here  = 1 , and the constant magnetic field has been set to zero , since a constant field has no significance in this discussion .
this system makes sense for any initial data , without symmetry constraints , and is called the relativistic vlasov - poisson system .
another special case of interest is the classical limit , obtained by letting c   in equations ( 12 , 13 , 14 ) , yielding : 19\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _
t}f + v \cdot { \nabla _ x}f + \beta e(t , x ) \cdot { \nabla _ v}f = 0,$$\end{document } 
 20\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$e = \nabla \phi , \quad \nabla \phi = \rho,$$\end{document } where  = 1 .
this is the vlasov - poisson system , and  = 1 corresponds to repulsive forces ( the plasma case ) .
taking  =  1 means attractive forces and the vlasov - poisson system is then a model for a newtonian self - gravitating system .
one of the fundamental problems in kinetic theory is to find out whether or not spontaneous shock formations will develop in a collision - less gas , i.e. , whether solutions to any of the equations above will remain smooth for all time , given smooth initial data .
if the initial data are small this problem has an affirmative solution in all cases considered above [ 81 , 86 , 32 , 33 ] . for initial data unrestricted in size
the picture is more involved . in order to obtain smooth solutions globally in time ,
the main issue is to control the support of the momenta 21\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t ) : = \sup \{\vert v\vert:\exists ( s , x ) \in [ 0,t ] \times { { \mathbb r}^3}\;{\rm{such}}\;{\rm{that}}\;f(s , x , v ) \neq 0\},$$\end{document } i.e. , to bound q(t ) by a continuous function so that q(t ) will not blow up in finite time . that such a control is sufficient for obtaining global existence of smooth solutions follows from well - known results in the different cases , cf .
[ 85 , 104 , 39 , 96 , 34 , 81 ] . for the full three - dimensional relativistic vlasov - maxwell system ,
the problem of establishing whether or not solutions will remain smooth for all time is open .
a different sufficient criterion for global existence in this case is given by pallard in , and he also shows a new bound for the electromagnetic field in terms of q(t ) in . in two space and three momentum dimensions , glassey and schaeffer [
82 , 83 ] have shown that q(t ) can be controlled for the relativistic vlasov - maxwell system , which thus yields global existence of smooth solutions in that case .
the relativistic and non - relativistic vlasov - poisson equations are very similar in form .
however , the mathematical results concerning the two systems are very different . in the non - relativistic case ,
batt   gave an affirmative solution in 1977 in the case of spherically - symmetric data .
pfaffelmoser   was the first one to give a proof for general smooth data . a simplified version of the proof is given by schaeffer in .
pfaffelmoser obtained the bound \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t ) \leq c{(1 + t)^{(51 + \delta)/11}},$$\end{document } where  > 0 can be taken as arbitrarily small .
the sharpest bound valid for  = 1 and  = 1 has been given by horst   and reads \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t ) \leq c(1 + t)\log
( 2 + t).$$\end{document } in the case of repulsive forces (  = 1 ) rein   has found a better estimate by using a new identity for the vlasov - poisson system , discovered independently by illner and rein   and by perthame .
rein s estimate reads \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t ) \leq c{(1 + t)^{2/3}}.$$\end{document } independently , and at about the same time as pfaffelmoser gave his proof , lions and perthame   used a different method for proving global existence .
their method is more generally applicable , and the two studies   and   are examples of problems in related systems , where their method has been successful . on the other hand , their method does not give such strong growth estimates on q(t ) as described above .
for the relativistic vlasov - poisson equation , glassey and schaeffer   showed in the case  = 1 that if the data are spherically symmetric , q(t ) can be controlled , which is analogous to the result by batt mentioned above .
also in the case of cylindrical symmetry they are able to control q(t ) ; see . if  = 1 it was shown in   that blow - up occurs in finite time for spherically - symmetric data with negative total energy .
more recently , lemou et al .   have investigated the structure of the blow - up solution .
they show that the blow - up is determined by the self - similar solution of the ultra - relativistic gravitational vlasov - poisson system .
it should be pointed out that the relativistic vlasov - poisson system is unphysical since it lacks the lorentz invariance ; it is a hybrid of a classical galilei invariant field equation and a relativistic transport equation ( 17 ) , cf . .
in particular , in the case  = 1 , it is not a special case of the einstein - vlasov system . only for spherically - symmetric data , in the case
the situation for weak solutions is different , in particular the existence of weak solutions to the relativistic vlasov - maxwell system is known [ 70 , 139 ] .
we also mention that models , which take into account both collisions and the electric and magnetic fields generated by the particles have been investigated .
classical solutions near a maxwellian for the vlasov - maxwell - boltzmann system are constructed by guo in .
a similar result for the vlasov - maxwell - landau system near a jttner solution is shown by guo and strain in .
we refer to the book by glassey   and the review article by rein   for more information on the relativistic vlasov - maxwell system and the vlasov - poisson system . before turning to the main theme of this review , i.e. , the einstein - vlasov system
nordstrm gravity   is an alternative theory of gravity introduced in 1913 . by coupling this model to a kinetic description of matter the nordstrm - vlasov system results . in nordstrm gravity
the nordstrm - vlasov system reads 22\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial _ t^2\phi - { \delta _ x}\phi = - { e^{4\phi}}\int\nolimits_{{{\mathbb r}^3 } } { { { { \mathfrak f}dp } \over { \sqrt { 1 + \vert p\vert ^{2}}}},}$$\end{document } 
 23\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ t}{\mathfrak f } + \hat{p } \cdot { \nabla _ x}{\mathfrak f } - \left [ { ( { \partial
_ t}\phi + \hat{p } \cdot { \nabla _ x}\phi)p + { { ( 1 + \vert p\vert ^{2})}^{- 1/2}}{\nabla _ x}\phi } \right ] \cdot { \nabla _ p}{\mathfrak f } = 0.$$\end{document } here \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\hat{p } = { p \over { \sqrt { 1 + \vert p\vert ^{2}}}},$$\end{document } denotes the relativistic velocity of a particle with momentum p. the mass of each particle , the gravitational constant , and the speed of light are all normalized to one . a solution ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$({\mathfrak f},\phi)$\end{document } ,  ) of this system is interpreted as follows .
the spacetime is a lorentzian manifold with a conformally - flat metric \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g_{\mu \nu } = e^{2\phi}{\rm diag}(- 1,1,1,1).$$\end{document}. the particle distribution f defined on the mass shell in this metric is given by 24\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f(t , x , p ) = { \mathfrak f}(t , x,{e^\phi}p).$$\end{document } the first mathematical study of this system was initiated by calogero in , where the existence of static solutions is established .
although the nordstrom - vlasov model of gravity does not describe physics correctly , the system approaches the vlasov - poisson system in the classical limit .
indeed , it is shown in   that solutions of the nordstrm - vlasov system tend to solutions of the vlasov - poisson system as the speed of light goes to infinity .
the cauchy problem was studied by several authors [ 51 , 50 , 15 , 108 , 131 ] and the question of global existence of classical solutions for general initial data was open for some time .
another interesting result for the nordstrm - vlasov system is given in , where a radiation formula , similar to the dipole formula in electrodynamics , is rigorously derived .
in this section we consider a self - gravitating collision - less gas in the framework of general relativity and we present the einstein - vlasov system . it is most often the case in the mathematics literature that the speed of light c and the gravitational constant g are normalized to one , but we keep these constants in the formulas in this section since in some problems they do play an important role .
however , in most of the problems discussed in the forthcoming sections these constants will be normalized to one .
let m be a four - dimensional manifold and let gab be a metric with lorentz signature (  , + , + , + ) so that ( m , gab ) is a spacetime .
the metric is assumed to be time - orientable so that there is a distinction between future and past directed vectors .
the possible values of the four - momentum p of a particle with rest mass m belong to the mass shell pm  tm , defined by 25\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p_m } : = \{({x^a},{p^a } ) \in tm:{g_{ab}}({x^a}){p^a}{p^b } = - { m^2}{c^2},{p^a}\;{\rm{is}}\;{\rm{future}}\;{\rm{directed}}\}.$$\end{document } hence , if m > 0 , pm(x ) is the set of all future - directed time - like vectors with length cm , and if m = 0 it is the set of all future - directed null vectors . on pm
we take ( x , p ) , a = 0 , 1 , 2 , 3 and j = 1 , 2 , 3 ( letters in the beginning of the alphabet always take values 0 , 1 , 2 , 3 and letters in the middle take 1 , 2 , 3 ) as local coordinates , and p is expressed in terms of p and the metric in view of equation ( 25 ) .
below we drop the index m on fm and simply write f. since we are considering a collisionless gas , the particles follow the geodesics in spacetime .
the geodesics are projections onto spacetime of the curves in pm defined in local coordinates by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { { d{x^a } } \over { ds } } = { p^a}},\\ { { { d{p^j } } \over { ds } } = - \gamma _ { bc}^j{p^b}{p^c}}.\\ \end{array}$$\end{document } here \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\gamma _ { bc}^a$\end{document } are the christoffel symbols . along a geodesic
the density function \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$f = f({x^a},{p^j})$\end{document } is invariant so that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${d \over { ds}}f({x^a}(s),{p^j}(s ) ) = 0,$$\end{document } which implies that 26\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p^a}{{\partial f } \over { \partial { x^a } } } - \gamma _ { ab}^j{p^a}{p^b}{{\partial f } \over { \partial { p^j } } } = 0.$$\end{document } this is accordingly the vlasov equation .
we point out that sometimes the density function is considered as a function on the entire tangent bundle tm rather than on the mass shell pm  tm . the vlasov equation for f = f(x , p ) then takes the form 27\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p^a}{{\partial f } \over { \partial { x^a } } } - \gamma _ { bc}^a{p^b}{p^c}{{\partial f } \over { \partial { p^a } } } = 0.$$\end{document } this equation follows from ( 26 ) if we take the mass shell condition gabpp = mc into account . indeed , by abuse of notation ,
we have \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { { \partial f } \over { \partial { x^a } } } = { { \partial f } \over { \partial { x^a } } } + { { \partial f } \over { \partial { p^0}}}{{\partial { p^0 } } \over { \partial { x^a}}},}\\ { { { \partial f } \over { \partial { p^j } } } = { { \partial f } \over { \partial { p^j } } } + { { \partial f } \over { \partial { p^0}}}{{\partial { p^0 } } \over { \partial { p^j}}}.}\\ \end{array}$$\end{document } here f is considered as a function on pm in the left - hand side , and on tm in the right - hand side . from the mass shell condition gabpp = mc we derive \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				
\usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { { \partial { p^0 } } \over { \partial { x^a } } } = - { 1 \over { { p_0}}}{p^b}{p_c}\gamma _ { ab}^c,}\\ { { { \partial { p^0 } } \over { \partial { p^j } } } = - { { { p_j } } \over { { p_0}}}.}\\ \end{array}$$\end{document } inserting these relations into ( 26 ) we obtain ( 27 ) .
if we let t = x , and divide the vlasov equation ( 26 ) by p we obtain the most common form in the literature of the vlasov equation 28\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\partial f } \over { \partial t } } + { { { p^j } } \over { { p^0}}}{{\partial f } \over { \partial { x^j } } } - { 1 \over { { p^0}}}\gamma _ { ab}^j{p^a}{p^b}{{\partial f } \over { \partial { p^j } } } = 0.$$\end{document } in a fixed spacetime the vlasov equation ( 28 ) is a linear hyperbolic equation for f and we can solve it by solving the characteristic system , 29\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{d{x^i } } \over { ds } } = { { { p^i } } \over { { p^0}}},$$\end{document } 
 30\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{d{p^i } } \over { ds } } = - \gamma _ { ab}^i{{{p^a}{p^b } } \over { { p^0}}}.$$\end{document } in terms of initial data f0 the solution of the vlasov equation can be written as 31\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f({x^a},{p^i } ) = { f_0}({x^i}(0,{x^a},{p^i}),{p^i}(0,{x^a},{p^i})),$$\end{document } where x(s , x , p ) and p(s , x , p ) solve equations ( 29 , 30 ) , and where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${x^i}(t,{x^a},{p^i } ) = { x^i}\;{\rm{and}}\;{p^i}(t,{x^a},{p^i } ) = { p^i}.$$\end{document } in order to write down the einstein - vlasov system we need to know the energy - momentum tensor \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t_{ab}^m$\end{document } in terms of f and gab .
we define 32\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{ab}^m = c\sqrt { \vert { g_{ab}}\vert } \int\nolimits_{{{\mathbb r}^3 } } { f{p_a}{p_b}{{d{p^1}d{p^2}d{p^3 } } \over { - { p_0}}}},$$\end{document } where , as usual , pa = gabp , and  gab denotes the absolute value of the determinant of gab .
we remark that the measure \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\mu : = { { \sqrt { \vert { g_{ab}}\vert } } \over { - { p_0}}}d{p^1}d{p^2}d{p^3},$$\end{document } is the induced metric of the submanifold \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${p_m}({x^a } ) \subset { t_{{x^a}}}m$\end{document } , and that  is invariant under lorentz transformations of the tangent space , and it is often the case in the literature that is written as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{ab}^m = c\int\nolimits_{{p_m}({x^a } ) } { f{p_a}{p_b}\mu.}$$\end{document } let us now consider a collisionless gas consisting of particles with different rest masses m1 , m2 ,  , mn , described by n density functions \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${f_{{m_j}}}$\end{document } , j = 1 ,  , n. then the vlasov equations for the different density functions \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${f_{{m_j}}}$\end{document } , together with the einstein equations , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${r_{ab } } - { 1 \over 2}r{g_{ab } } + \lambda { g_{ab } } = { { 8\pi g } \over { { c^4}}}\sum\limits_{k = 1}^n { t_{ab}^{{m_k}},}$$\end{document } form the einstein - vlasov system for the collision - less gas . here
rab is the ricci tensor , r is the scalar curvature and  is the cosmological constant .
henceforth , we always assume that there is only one species of particles in the gas and we write tab for its energy momentum tensor .
moreover , in what follows , we normalize the rest mass m of the particles , the speed of light c , and the gravitational constant g , to one , if not otherwise explicitly stated that this is not the case .
we define the particle current density \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${n^a } = - \int\nolimits_{{{\mathbb r}^3 } } { f{p^a } } \sqrt { \vert { g_{ab}}\vert } { { d{p^1}d{p^2}d{p^3 } } \over { { p_0}}}.$$\end{document } using normal coordinates based at a given point and assuming that f is compactly supported , it is not hard to see that tab is divergence - free , which is a necessary compatibility condition since the left - hand side of ( 2 ) is divergence - free by the bianchi identities .
a computation in normal coordinates also shows that n is divergence - free , which expresses the fact that the number of particles is conserved .
if v is a future - directed time - like or null vector then we have nav  0 with equality if and only if f = 0 at the given point .
hence , n is always future - directed time - like , if there are particles at that point . moreover , if v and w are future - directed time - like vectors then tabvw  0 , which is the dominant energy condition .
this is called the non - negative pressure condition , and it implies that the strong energy condition holds as well . that the energy conditions hold for vlasov matter is one reason that the vlasov equation defines a well - behaved matter model in general relativity .
another reason is the well - posedness theorem by choquet - bruhat   for the einstein - vlasov system that we state below . before stating that theorem we first discuss the conditions imposed on the initial data .
the initial data in the cauchy problem for the einstein - vlasov system consist of a 3-dimensional manifold s , a riemannian metric gij on s , a symmetric tensor kij on s , and a non - negative scalar function f0 on the tangent bundle ts of s. the relationship between a given initial data set ( gij , kij ) on s and the metric gab on the spacetime manifold is , that there exists an embedding  of s into the spacetime such that the induced metric and second fundamental form of (s ) coincide with the result of transporting ( gij , kij ) with . for the relation of the distribution functions f and f0 we have to note that f is defined on the mass shell .
the initial condition imposed is that the restriction of f to the part of the mass shell over (s ) should be equal to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${f_0 } \circ ( { \psi ^{- 1}},d{(\psi)^{- 1 } } ) \circ \phi$\end{document } , where  sends each point of the mass shell over (s ) to its orthogonal projection onto the tangent space to (s ) .
an initial data set for the einstein - vlasov system must satisfy the constraint equations , which read 33\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$r - { k_{ij}}{k^{ij } } + { ( { \rm{tr}}\;k)^2 } = 16\pi \rho,$$\end{document } 
 34\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ i}k_l^i - { \nabla _ l}({\rm{tr}}\;k ) = 8\pi j_{l}.$$\end{document } here  = tabnn 

j = htbcn , where n is the future directed unit normal vector to the initial hypersurface , and h = g + nn is the orthogonal projection onto the tangent space to the initial hypersurface . in terms of f0
we can express  and j by ( j satisfies naj = 0 , so it can naturally be identified with a vector intrinsic to s ) \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { \rho = - \int\nolimits_{{{\mathbb r}^3 } } { { f_0}{p^a}{p_a } } \sqrt { \vert { g_{ij}}\vert } { { d{p^1}d{p^2}d{p^3 } } \over { 1 + { p_j}{p^j}}},}\\ { { j_l } = \int\nolimits_{{{\mathbb r}^3 } } { { f_0}{p_l } } \sqrt { \vert { g_{ij}}\vert } d{p^1}d{p^2}d{p^3}.}\\ \end{array}$$\end{document } we can now state the local existence theorem by choquet - bruhat , for the einstein - vlasov system . theorem 1 
 let s be a 3-dimensional manifold , gij 
 a smooth riemannian metric on s , kij 
 a smooth symmetric tensor on s and f0 
 a smooth non - negative function of compact support on the tangent bundle ts of s. suppose that these objects satisfy the constraint equations ( 33 , 34 ) . then there exists a smooth spacetime ( m , gab ) ,
a smooth distribution function f on the mass shell of this spacetime , and a smooth embedding  into m , which induces the given initial data on s such that gab 
 and f satisfy the einstein - vlasov system and (s ) is a cauchy surface .
moreover , given any other spacetime ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$({m^\prime},g_{ab}^\prime)$\end{document } ) , distribution function f and embedding  satisfying these conditions , there exists a diffeomorphism  from an open neighborhood of if (s ) in m to an open neighborhood of (s ) in m , which satisfies    and carries gab 
 and f to 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$g_{ab}^\prime$\end{document } 
 and f , respectively .
the above formulation is in the case of smooth initial data ; for information on the regularity needed on the initial data we refer to   and . in this context we also mention that local existence has been proven for the yang - mills - vlasov system in , and that this problem for the einstein - maxwell - boltzmann system is treated in .
however , this result is not complete , as the non - negativity of f is left unanswered . also , the hypotheses on the scattering kernel in this work leave some room for further investigation . the local existence problem for physically reasonable assumptions on the scattering kernel
does not seem well understood in the context of the einstein - boltzmann system , and a careful study of this problem would be desirable .
the mathematical study of the einstein - boltzmann system has been very sparse in the last few decades , although there has been some activity in recent years . since most questions on the global properties are completely open let us only very briefly mention some of these works .
global existence for the homogeneous einstein - boltzmann system in robertson - walker spacetimes is proven in , and a generalization to bianchi type i symmetry is established in . in the following sections we present results on the global properties of solutions of the einstein - vlasov system , which have been obtained during the last two decades . before ending this section
we mention a few other sources for more background on the einstein - vlasov system , cf .
in this section , we discuss results on global existence and on the asymptotic structure of solutions of the cauchy problem in the asymptotically - flat case . in general relativity
two classes of initial data are distinguished in the study of the cauchy problem : asymptotically - flat initial data and cosmological initial data .
the initial hypersurface is topologically  and appropriate fall - off conditions are imposed to ensure that far away from the body spacetime is approximately flat .
spacetimes , which possess a compact cauchy hypersurface , are called cosmological spacetimes , and data are accordingly given on a compact 3-manifold . in this case , the whole universe is modeled rather than an isolated body .
one can also consider a case , which is un - physical , in which spacetime is asymptotically flat except in one direction , namely cylindrically - symmetric spacetimes , cf .
the majority of the work so far has been devoted to the spherically - symmetric case but recently a result on static axisymmetric solutions has been obtained .
in contrast to the asymptotically - flat case , cosmological spacetimes admit a large number of symmetry classes .
this provides the possibility to study many special cases for which the difficulties of the full einstein equations are reduced .
the main goal of these studies is to provide an answer to the weak and strong cosmic censorship conjectures , cf .
the study of the global properties of solutions to the spherically - symmetric einstein - vlasov system was initiated two decades ago by rein and rendall , cf . also [ 135 , 156 ] .
they chose to work in coordinates where the metric takes the form \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - { e^{2\mu ( t , r)}}d{t^2 } + { e^{2\lambda ( t , r)}}d{r^2 } + { r^2}(d{\theta ^2 } + { \sin ^2}\theta d{\varphi ^2}),$$\end{document } where t   , r  0 ,   [ 0, ] ,   [ 0 , 2 ] . these are called schwarzschild coordinates .
asymptotic flatness is expressed by the boundary conditions \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\underset{r \rightarrow \infty}{\lim } \lambda ( t , r ) = \underset{r \rightarrow \infty}{\lim } \mu ( t , r ) = 0,\quad \forall t \geq 0.$$\end{document } a regular center is also required and is guaranteed by the boundary condition \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\lambda ( t,0 ) = 0\quad \forall t \geq 0.$$\end{document } the coordinates ( r ,  ,  ) give rise to difficulties at r = 0 and it is advantageous to use cartesian coordinates . with \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				
\usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$x = ( r\sin \phi \cos \theta , r\sin \phi \sin \theta , r\cos \phi)$$\end{document } as spatial coordinates and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v^j } = { p^j } + ( { e^\lambda } - 1){{x \cdot p } \over r}{{{x^j } } \over r}$$\end{document } as momentum coordinates , the einstein - vlasov system reads 35\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _
t}f + { e^{\mu - \lambda}}{v \over { \sqrt { 1 + \vert v\vert ^{2 } } } } \cdot { \nabla _ x}f - \left({{\lambda
_ t}{{x \cdot v } \over r } + { e^{\mu - \lambda}}{\mu _ r}\sqrt { 1 + \vert v\vert ^{2 } } } \right){x \over r } \cdot { \nabla _ v}f = 0,$$\end{document } 
 36\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${e^{- 2\lambda}}(2r{\lambda _ r } - 1 ) + 1 = 8\pi { r^2}\rho,$$\end{document } 
 37\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${e^{- 2\lambda}}(2r{\mu _ r } + 1 ) - 1 = 8\pi { r^2}p,$$\end{document } 
 38\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\lambda _ t } = - 4\pi r{e^{\lambda + \mu}}j,$$\end{document } 
 39\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${e^{- 2\lambda}}({\mu _
{ rr } } + ( { \mu _ r } - { \lambda _ r})({\mu _ r } + { 1 \over r } ) ) - { e^{- 2\mu}}({\lambda _ { tt } } + { \lambda _ t}({\lambda _ t } - { \mu _ t } ) ) = 8\pi { p_t}.$$\end{document } the matter quantities are defined by 40\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho ( t , x ) = \int\nolimits_{{{\mathbb r}^3 } } { \sqrt { 1 + \vert v\vert ^{2 } } f(t , x , v)dv,}$$\end{document } 
 41\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p(t , x ) = \int\nolimits_{{{\mathbb r}^3 } } { { { \left({{{x \cdot v } \over r } } \right)}^2}f(t , x , v){{dv } \over { \sqrt { 1 + \vert v\vert ^{2}}}},}$$\end{document } 
 42\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$j(t , x ) = \int\nolimits_{{{\mathbb r}^3 } } { { { x\cdot v } \over r}f(t , x , v)dv,}$$\end{document } 
 43\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p_t}(t , x ) = { 1 \over 2}\int\nolimits_{{{\mathbb r}^3 } } { \vert { { x \times v } \over r}\vert ^{2}f(t , x , v)dv.}$$\end{document } here  is the energy density , j the current , p the radial pressure , and pt the tangential pressure .
let us point out that these equations are not independent , e.g. , equations ( 38 ) and ( 39 ) follow from ( 35)(37 ) . as initial data we take a spherically - symmetric , non - negative , and continuously diferentiable function f0 with compact support that satisfies 44\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\int\nolimits_{\vert y\vert < r } { \int\nolimits_{{{\mathbb r}^3 } } { \sqrt { 1 + \vert v{\vert ^2 } } { f_0}(y , v)dv\;dy\ ; < { r \over 2}.}}$$\end{document } this condition guarantees that no trapped surfaces are present initially .
the schwarzschild coordinates have the advantage that the resulting system of equations can be written in a quite condensed form .
moreover , for most initial data , solutions are expected to exist globally in schwarzschild time , which sometimes is called the polar time gauge .
let us point out here that there are initial data leading to spacetime singularities , cf .
hence , the question of global existence for general initial data is only relevant if the time slicing of the spacetime is expected to be singularity avoiding , which is the case for schwarzschild time .
however , one disadvantage is that these coordinates only cover a relatively small part of the spacetime , in particular trapped surfaces are not admitted .
hence , to analyze the black - hole region of a solution these coordinates are not appropriate . here
we only mention the other coordinates and time gauges that have been considered in the study of the spherically symmetric einstein - vlasov system .
the einstein - vlasov system is investigated in double null coordinates in [ 64 , 63 ] .
maximal - areal coordinates and eddington - finkelstein coordinates are used in [ 21 , 17 ] , and in   respectively . in   it
is shown that for initial data satisfying ( 44 ) there exists a unique , continuously - differentiable solution f with f(0 ) = f0 on some right maximal interval [ 0 , t ) .
if the solution blows up in finite time , i.e. , if t <  , then (t ) becomes unbounded as t  t. moreover , a continuation criterion is shown that says that a local solution can be extended to a global one , provided q(t ) can be bounded on [ 0 , t ) , where 45\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t):\sup \{\vert v\vert : \exists ( s , x ) \in [ 0,t ] \times { { \mathbb r}^3}\;{\rm{such}}\;{\rm{that}}\;f(s , x , v ) \neq 0\}.$$\end{document } this is analogous to the situation for the vlasov - maxwell system .
a control of the v - support immediately implies that  and p are bounded in view of equations ( 40 , 41 ) . in the vlasov - maxwell case the field equations have a regularizing effect in the sense that derivatives can be expressed through spatial integrals , and it follows   that the derivatives of f can also be bounded if the v - support is bounded . for the einstein - vlasov system such
a regularization is less clear , since , e.g. , r depends on p in a point - wise manner .
however , in view of equation ( 39 ) certain combinations of second and first order derivatives of the metric components can be expressed in terms of the matter component pt , which is a consequence of the geodesic deviation equation .
this fact turns out to also be sufficient for obtaining bounds on the derivatives of f , cf . [ 142 , 135 , 156 ] for details . the local existence result discussed above holds for compactly - supported initial data .
the compact support condition in the momentum variables is in   replaced by the fall - off condition 46\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\underset{(x , v ) \in { { \mathbb r}^6}}{\sup}{(1 + \vert v\vert)^5}\vert \overset{\circ}{f}(x , v)\vert < \infty.$$\end{document } we also refer to   where a subclass of non - compactly - supported data is treated . local existence of solutions in double null coordinates and in eddington - finkelstein coordinates is established in , and   respectively . in  
the authors also consider the problem of global existence in schwarzschild coordinates for small initial data for massive particles .
they show that for such data the v - support is bounded on [ 0 , t ) .
hence , the continuation criterion implies that t = . the resulting spacetime in   is geodesically complete , and the components of the energy - momentum tensor as well as the metric quantities decay with certain algebraic rates in t. the mathematical method used by rein and rendall is inspired by the analogous small data result for the vlasov - poisson equation by bardos and degond
. this should not be too surprising since for small data the gravitational fields are expected to be small and a newtonian spacetime should be a fair approximation .
in this context we point out that in   it is proven that the vlasov - poisson system is indeed the non - relativistic limit of the spherically - symmetric einstein - vlasov system , i.e. , the limit when the speed of light c  . in   this result is shown without symmetry assumptions .
as mentioned above the local and global existence problem has been studied using other time gauges , in particular rendall has shown global existence for small initial data in maximal - isotropic coordinates in .
the previous results refer to massive particles but they do not immediately carry over to massless particles .
this case is treated by dafermos in   where global existence for small initial data is shown in double null coordinates .
the spacetimes obtained in the studies [ 142 , 156 , 63 ] are all causally geodesically complete and appropriate decay rates of the metric and the matter quantities are given .
in the case of small initial data the resulting spacetime is geodesically complete and no singularities form .
a different scenario , which leads to a future geodesically complete spacetime , is to consider initial data where the particles are moving rapidly outwards .
if the particles move sufficiently fast the matter disperses and the gravitational attraction is not strong enough to reverse the velocities of the particles to create a collapsing system .
it is shown that the scenario described above can be realized , and that global existence holds . in section 3.7
we discuss results on the formation of black holes and trapped surfaces ; in particular , the results in   will be presented
. a corollary of the main result in   concerns the issue of global existence and thus we mention it here .
it is shown that a particular class of initial data , which lead to formation of black holes , have the property that the solutions exist for all schwarzschild time .
the initial data consist of two parts : an inner part , which is a static solution of the einstein - vlasov system , and an outer part with matter moving inwards .
the set - up is shown to preserve the direction of the momenta of the outer part of the matter , and it is also shown that in schwarzschild time the inner part and the outer part of the matter never interact in schwarzschild time . as was mentioned at the end of section 3.1 ,
the issue of global existence for general initial data is only relevant in certain time gauges since there are initial data leading to singular spacetimes .
however , it is reasonable to believe that global existence for general data may hold in a polar time gauge or a maximal time gauge , cf . , and it is often conjectured in the literature that these time slicings are singularity avoiding .
however , there is no proof of this statement for any matter model and it would be very satisfying to provide an answer to this conjecture for the einstein - vlasov system
. a proof of global existence in these time coordinates would also be of great importance due to its relation to the weak cosmic censorship conjecture , cf .
[ 61 , 62 , 65 ] . the methods of proofs in the cases described in sections 3.3 and 3.4 , where global existence has been shown , are all tailored to treat special classes of initial data and they will likely not apply in more general situations . in this section
these results are all conditional in the sense that assumptions are made on the solutions , and not only on the initial data .
the first study on global existence for general initial data is , which is carried out in schwarzschild coordinates .
the authors introduce the following variables in the momentum space adapted to spherical symmetry , 47\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$l : = \vert x{\vert ^2}\vert v{\vert ^2 } - { ( x \cdot v)^2},\;w = { { x \cdot v } \over r},$$\end{document } where l is the square of the angular momentum and w is the radial component of the momenta .
a consequence of spherical symmetry is that angular momentum is conserved along the characteristics . in these variables
the vlasov equation for f = f ( t , r , w , l ) becomes 48\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _
t}w + { e^{\mu - \lambda}}{\mu _ r}e - { e^{\mu - \lambda}}{l \over { { r^3}e } } } \right){\partial _
w}f = 0,$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$e = e(r , w , l ) = \sqrt { 1 + { w^2 } + l/{r^2}}.$$\end{document } the main result in   shows that as long as there is no matter in the ball \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\{x \in { { \mathbb r}^3}:\vert x\vert \leq \epsilon \},$$\end{document } the estimate 49\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t ) \leq { e^{\log q(0){e^{c(\epsilon)t}}}},$$\end{document } holds .
here c( ) is a constant , which depends on . thus , in view of the continuation criterion this can be viewed as a global existence result outside the center of symmetry for initial data with compact support .
this result rules out shell - crossing singularities , which are present when , e.g. , dust is used as a matter model .
the bound of q is obtained by estimating each term individually in the characteristic equation associated with the vlasov equation ( 48 ) for the radial momentum .
the einstein equations imply that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mu _ r } = { m \over { { r^2}}}{e^{2\lambda } } + 4\pi rp{e^{2\lambda}}$$\end{document } where 50\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$m(t,\,r ) = 4\pi \int\nolimits_0^r { { \eta ^2}\rho ( t,\eta ) d\eta,}$$\end{document } is the quasi local mass .
thus , using ( 38 ) the characteristic equation consists of the two terms t1 = 4re , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${t_2 } = { e^{\mu + \lambda}}{m \over { { r^2}}}$\end{document } , together with a term , which is independent of the matter quantities .
there is a distinct difference between the terms t1 and t2 due to the fact that m can be regarded as an average , since it is given as a space integral of the energy density  , whereas j and p are point - wise terms .
the method in   makes use of a cancellation property of the radial momenta in t1 so that outside the center this term is manageable but in general it seems very unpleasant to have to treat point - wise terms of this kind . in   rendall
the bound on q(t ) is again obtained by estimating each term in the characteristic equation . in this case
. however , the terms are , in analogy with the schwarzschild case , strongly singular at the center .
in particular , the method avoids the point - wise terms by using the fact that the characteristic system can be written in a form such that green s formula in the plane can be applied .
this results in a combination of terms involving second - order derivatives , which can be substituted for by one of the einstein equations .
this method was first introduced in   but the set - up is different in   and the application of green s formula becomes very natural .
in addition , the bound of q is improved compared to ( 49 ) and reads \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q(t ) \leq q(0){e^{c(1 + t)/\epsilon t}}.$$\end{document } this bound is sufficient to conclude that global existence outside the center also holds for non - compact initial data .
in addition to the global existence result outside the centre , it is shown in   that as long as 3m(t , r )  r and j  0 , singularities can not form .
note that in schwarzschild coordinates 2m(t , r )  r always , and that there are closed null geodesics if 3 m = r in the schwarzschild static spacetime .
the method in   also applies to the case of maximal - isotropic coordinates studied in .
there is an improvement concerning the regularity of the terms that need to be estimated to obtain global existence in the general case .
a consequence of   is accordingly that the quite different proofs in   and in   are put on the same footing .
we point out that the method can also be applied to the case of maximal - areal coordinates .
the results discussed above concern time gauges , which are expected to be singularity avoiding so that the issue of global existence makes sense .
an interpretation of these results is that  first singularities  ( where the notion of  first  is tied to the causal structure ) , in the non - trapped region , must emanate from the center and that this case has also been shown in double null - coordinates by dafermos and rendall in . the main motivation for studying the system in these coordinates has its origin from the method of proof of the cosmic - censorship conjecture for the einstein - scalar field system by christodoulou .
an essential part of his method is based on the understanding of the formation of trapped surfaces . in   it
is shown that a single trapped surface or marginally - trapped surface in the maximal development implies that weak cosmic censorship holds the theorem holds true for any spherically - symmetric matter spacetime if the matter model is such that  first  singularities necessarily emanate from the center .
the results in   and in   are not sufficient for concluding that the hypothesis of the matter needed in the theorem in   is satisfied , since they concern a portion of the maximal development covered by particular coordinates
. therefore , dafermos and rendall   choose double - null coordinates , which cover the maximal development , and they show that the mentioned hypothesis is satisfied for vlasov matter .
the main reason that the question of global existence in certain time coordinates discussed in the previous section 3.5 is of great importance is its relation to the cosmic censorship conjectures .
now there is , in fact , no theorem in the literature , which guarantees that weak cosmic censorship follows from such a global existence result , but there are strong reasons to believe that this is the case , cf .   and .
hence , if initial data can be constructed , which lead to naked singularities , then either the conjecture that global existence holds generally is false or the viewpoint that global existence implies the absence of naked singularities is wrong . in view of a recent result by rendall and velazquez   on self similar dust - like solutions for the massless einstein - vlasov system , this issue has much current interest .
let us mention here that there is a previous study on self - similar solutions in the massless case by martn - garca and gundlach . however , this result is based on a scaling of the density function itself and therefore makes the result less related to the cauchy problem .
also , their proof is , in part , based on numerics , which makes it harder to judge the relevance of the result .
the main aim of the work   is to establish self - similar solutions of the massive einstein - vlasov system and the present result can be viewed as a first step to achieving this . in the set - up , two simplifications are made .
first , the authors study the massless case in order to find a scaling group , which leaves the system invariant .
more precisely , the massless system is invariant under the scaling \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$r \rightarrow \theta r,\;\;t \rightarrow \theta t,\;\;w \rightarrow { 1 \over { \sqrt \theta}}w,\;\;l \rightarrow \theta l.$$\end{document } the massless assumption seems not very restrictive since , if a singularity forms , the momenta will be large and therefore the influence of the rest mass of the particles will be negligible , so that asymptotically the solution can be self - similar also in the massive case , cf .
the second simplification is that the possible radial momenta are restricted to two values , which means that the density function is a distribution in this variable .
thus , the solutions can be thought of as intermediate between smooth solutions of the einstein - vlasov system and dust . for this simplified system
it turns out that the existence question of self - similar solutions can be reduced to that of the existence of a certain type of solution of a four - dimensional system of ordinary differential equations depending on two parameters .
the proof is based on a shooting argument and involves relating the dynamics of solutions of the four - dimensional system to that of solutions of certain two- and three - dimensional systems obtained from it by limiting processes .
the reason that an ode system is obtained is due to the assumption on the radial momenta , and if regular initial data is considered , an ode system is not sufficient and a system of partial differential equations results .
the solution is not asymptotically flat but there are ideas outlined in   of how this can be overcome .
it should be pointed out here that a similar problem occurs in the work by christodoulou   for a scalar field , where the naked singularity solutions are obtained by truncating self - similar data .
the singularity of the self - similar solution by rendall and velazquez is real in the sense that the kretschmann scalar curvature blows up . the asymptotic structure of the solution is striking in view of the conditional global existence result in .
indeed , the self similar solution is such that j  0 , and 3m(t , r )  r asymptotically , but for any t , 3m(t , r ) > r for some t > t. in   global existence follows if j  0 and if 3m(t , r )  r for all t. it is also the case that if m / r is close to 1/2 , then global existence holds in certain situations , cf . .
hence , the asymptotic structure of the self - similar solution has properties , which have been shown to be difficult to treat in the search for a proof of global existence .
we have previously mentioned that there exist initial data for the spherically - symmetric einstein - vlasov system , which lead to formation of black holes . the first result in this direction
he shows that there exist initial data for the spherically - symmetric einstein - vlasov system such that a trapped surface forms in the evolution .
the occurrence of a trapped surface signals the formation of an event horizon . as mentioned above , dafermos   has proven that , if a spherically - symmetric spacetime contains a trapped surface and the matter model satisfies certain hypotheses , then weak cosmic censorship holds true . in   it
hence , by combining these results it follows that initial data exist , which lead to gravitational collapse and for which weak cosmic censorship holds . however , the proof in   rests on a continuity argument , and it is not possible to tell whether or not a given initial data set will give rise to a black hole .
this is in contrast to the result in , where explicit conditions on the initial data are given , which guarantee the formation of trapped surfaces in the evolution .
the analysis is carried out in eddington - finkelstein coordinates and a central result in   is to control the life span of the solution to ensure that there is sufficient time to form a trapped surface before the solution may break down . in particular , weak cosmic censorship holds for these initial data . in   the formation of the event horizon in gravitational collapse
one inner part and one outer part , in which all particles move inward initially .
the reason for the inner part is that it is possible to choose the parameters for the data such that the particles of the outer matter part continue to move inward for all schwarzschild time as long as the particles do not interact with the inner part .
this fact simplifies the analysis since the dynamics is much restricted when the particles keep the direction of their radial momenta .
the main result is that explicit conditions on the initial data with adm mass m are given such that there is a family of outgoing null geodesics for which the area radius r along each geodesic is bounded by 2 m .
it is furthermore shown that if \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t \geq 0,\;{\rm{and}}\;r \geq 2 m + \alpha { e^{- \beta t}},$$\end{document } where  and  are positive constants , then f(t , r ,  ,  ) = 0 , and the metric equals the schwarzschild metric 51\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - ( 1 - { { 2 m } \over r})d{t^2 } + { ( 1 - { { 2 m } \over r})^{- 1}}d{r^2 } + { r^2}(d{\theta ^2 } + { \sin ^2}\theta d{\phi ^2}),$$\end{document } representing a black hole with mass m. hence , spacetime converges asymptotically to the schwarzschild metric . the latter result does not reveal whether or not all the matter crosses r = 2 m or simply piles up at the event horizon . in   it is shown that for initial data , which are closely related to those in , but such that the radial momenta are unbounded , all the matter do cross the event horizon asymptotically in schwarzschild time .
this is in contrast to what happens to freely - falling observers in a static schwarzschild spacetime , since they will never reach the event horizon .
the result in   is reconsidered in , where an additional argument is given to match the definition of weak cosmic censorship given in .
it is natural to relate the results of [ 20 , 24 ] to those of christodoulou on the spherically - symmetric einstein - scalar field system   and . in   it
is shown that if the final bondi mass m is different from zero , the region exterior to the sphere r = 2 m tends to the schwarzschild metric with mass m similar to the result in . in   explicit conditions on the initial data
this paper played a crucial role in christodoulou s proof   of the weak and strong cosmic censorship conjectures .
the conditions on the initial data in   allow the ratio of the hawking mass and the area radius to cover the full range , i.e. , 2m / r  ( 0 , 1 ) , whereas the conditions in   require 2m / r to be close to one .
hence , it would be desirable to improve the conditions on the initial data in , although the conditions by christodoulou for a scalar field are not expected to be sufficient in the case of vlasov matter . in   a numerical study on critical collapse for the einstein - vlasov system
it has been shown by rein and rodewis   that the numerical scheme has desirable convergence properties .
( in the vlasov - poisson case , convergence was proven in , see also ) .
the speculation discussed above that there may be no naked singularities formed for any regular initial data is in part based on the fact that the naked singularities that occur in scalar field collapse appear to be associated with the existence of type ii critical collapse , while vlasov matter is of type i. the primary goal in   was indeed to decide whether vlasov matter is type i or type ii .
given small initial data , no black holes form and matter will disperse . for large data ,
black holes will form and consequently there is a transition regime separating dispersion of matter and formation of black holes .
if we introduce a parameter a on the initial data such that for small a dispersion occurs and for large a a black hole is formed , we get a critical value ac separating these regions .
if we take a > ac and denote by mb(a ) the mass of the black hole , then if mb(a )  0 as a  ac we have type ii matter , whereas for type i matter this limit is positive and there is a mass gap . for more information on critical collapse
the conclusion of   is that vlasov matter is of type i. there are two other independent numerical simulations on critical collapse for vlasov matter [ 128 , 21 ] . in these simulations ,
we end this section with a discussion of the spherically - symmetric einstein - vlasov - maxwell system , i.e. , the case considered above with charged particles . whereas the constraint equations in the uncharged case , written in schwarzschild coordinates , do not involve solving any difficulties once the distribution function is given , the charged case is more challenging .
however , in   it is shown that solutions to the constraint equations do exist for the einstein - vlasov - maxwell system . in   local existence
is shown together with a continuation criterion , and global existence for small initial data is shown in .
in this section we discuss the einstein - vlasov system for cosmological spacetimes , i.e. , spacetimes that possess a compact cauchy surface .
the  particles  in the kinetic description are in this case galaxies or even clusters of galaxies .
the main goal is to determine the global properties of the solutions to the einstein - vlasov system for initial data given on a compact 3-manifold . in order to do so , a global time
coordinate t must be found and the asymptotic behavior of the solutions when t tends to its limiting values has to be analyzed .
this might correspond to approaching a singularity , e.g. , the big bang singularity , or to a phase of unending expansion .
presently , the aim of most of the studies of the cosmological cauchy problem has been to show existence for unrestricted initial data and the results that have been obtained are in cases with symmetry ( see , however , , where to some extent global properties are shown in the case without symmetry ) .
a recent and very extensive work by ringstrm has , on the other hand , a different aim , i.e. , to show stability of homogeneous cosmological models , and concerns the general case without symmetry .
the size of the cauchy data is in this case very restricted but , since ringstrm allows general perturbations , there are no symmetries available to reduce the complexity of the einstein - vlasov system .
the only spatially - homogeneous spacetimes admitting a compact cauchy surface are the bianchi types i , ix and the kantowski - sachs model ; to allow for cosmological solutions with more general symmetry types , it is enough to replace the condition that the spacetime is spatially homogeneous , with the condition that the universal covering of spacetime is spatially homogeneous .
spacetimes with this property are called locally spatially homogeneous and these include , in addition , the bianchi types ii , iii , v , vi0 , vii0 , and viii .
one of the first studies on the einstein - vlasov system for spatially - homogeneous spacetimes is the work   by rendall .
he chooses a gaussian time coordinate and investigates the maximal range of this time coordinate for solutions evolving from homogeneous data . for bianchi ix and for kantowski - sachs spacetimes
he finds that the range is finite and that there is a curvature singularity in both the past and the future time directions . for the other bianchi types
there is a curvature singularity in the past , and to the future spacetime is causally geodesically complete .
although the questions on curvature singularities and geodesic completeness are very important , it is also desirable to have more detailed information on the asymptotic behavior of the solutions , and , in particular , to understand in which situations the choice of matter model is essential for the asymptotics . in recent years
several studies on the einstein - vlasov system for spatially locally homogeneous spacetimes have been carried out with the goal to obtain a deeper understanding of the asymptotic structure of the solutions .
roughly , these investigations can be divided into two cases : ( i ) studies on non - locally rotationally symmetric ( non - lrs ) bianchi i models and ( ii ) studies of lrs bianchi models . in case ( i )
rendall shows in   that solutions converge to dust solutions for late times . under the additional assumption of small initial data
this result is extended by nungesser , who gives the rate of convergence of the involved quantities . in   rendall
also raises the question of the existance of solutions with complicated oscillatory behavior towards the initial singularity may exist for vlasov matter , in contrast to perfect fluid matter .
note that for a perfect fluid the pressure is isotropic , whereas for vlasov matter the pressure may be anisotropic , and this fact could be sufficient to drastically change the dynamics .
this question is answered in , where the existence of a heteroclinic network is established as a possible asymptotic state .
this implies a complicated oscillating behavior , which differs from the dynamics of perfect fluid solutions .
the results in   were then put in a more general context by calogero and heinzle , where quite general anisotropic matter models are considered . in case ( ii ) the asymptotic behaviour of solutions has been analyzed in [ 159 , 160 , 48 , 47 ] . in ,
the case of massless particles is considered , whereas the massive case is studied in .
both the nature of the initial singularity and the phase of unlimited expansion are analyzed .
the main concern in these two works is the behavior of bianchi models i , ii , and iii .
a general conclusion is that the choice of matter model is very important since , for all symmetry classes studied , there are differences between the collision - less model and a perfect fluid model , both regarding the initial singularity and the expanding phase .
the most striking example is for the bianchi ii models , where they find persistent oscillatory behavior near the singularity , which is quite different from the known behavior of bianchi type ii perfect fluid models . in   it
is also shown that solutions for massive particles are asymptotic to solutions with massless particles near the initial singularity . for bianchi i and
ii , it is also proven that solutions with massive particles are asymptotic to dust solutions at late times .
the investigation   concerns a large class of anisotropic matter models , and , in particular , it is shown that solutions of the einstein - vlasov system with massless particles oscillate in the limit towards the past singularity for bianchi ix models .
considers the homogeneous spacetimes with a cosmological constant for all bianchi models except bianchi type ix .
she also obtains the time decay of the components of the energy momentum tensor as t   , and she shows that spacetime is asymptotically dust - like .
anguige   studies the conformal einstein - vlasov system for massless particles , which admit an isotropic singularity .
he shows that the cauchy problem is well posed with data consisting of the limiting density function at the singularity . in the spatially homogeneous case
the metric can be written in a form that is independent of the spatial variables and this leads to an enormous simplification .
another class of spacetimes that are highly symmetric but require the metric to be spatially dependent are those that admit a group of isometries acting on two - dimensional spacelike orbits , at least after passing to a covering manifold .
the group may be two - dimensional ( local u(1 )  u(1 ) or t symmetry ) or three - dimensional ( spherical , plane , or hyperbolic symmetry ) .
in all these cases , the quotient of spacetime by the symmetry group has the structure of a two - dimensional lorentzian manifold q. the orbits of the group action ( or appropriate quotients in the case of a local symmetry ) are called surfaces of symmetry .
thus , there is a one - to - one correspondence between surfaces of symmetry and points of q. there is a major difference between the cases where the symmetry group is two- or three - dimensional . in the three - dimensional case
no gravitational waves are admitted , in contrast to the two - dimensional case where the evolution part of the einstein equations are non - linear wave equations .
three types of time coordinates that have been studied in the inhomogeneous case are cmc , areal , and conformal coordinates .
a cmc time coordinate t is one where each hypersurface of constant time has constant mean curvature and on each hypersurface of this kind the value of t is the mean curvature of that slice . in the case of areal coordinates , the time coordinate is a function of the area of the surfaces of symmetry , e.g. , proportional to the area or proportional to the square root of the area . in the case of conformal coordinates ,
the advantage with a cmc approach is that the definition of a cmc hypersurface does not depend on any symmetry assumptions and it is possible that cmc foliations will exist for general spacetimes .
the areal coordinate foliation , on the other hand , is adapted to the symmetry of spacetime but it has analytical advantages and detailed information about the asymptotics can be derived .
the conformal coordinates have mainly served as a useful framework for the analysis to obtain geometrically - based time foliations .
let us now consider spacetimes ( m , g ) admitting a three - dimensional group of isometries .
the topology of m is assumed to be   s  f , with f a compact two - dimensional manifold .
the universal covering \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\hat f}$\end{document } of f induces a spacetime \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(\hat m,\hat g)$\end{document } by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\hat m = { \mathbb r } \times { s^1 } \times \hat f$\end{document } and  = p*g , where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				
a three - dimensional group g of isometries is assumed to act on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(\hat m,\hat g)$\end{document}. if f = s and g = so(3 ) , then ( m , g ) is called spherically symmetric ; if f = t and g = e2 ( euclidean group ) , then ( m , g ) is called plane symmetric ; and if f has genus greater than one and the connected component of the symmetry group g of the hyperbolic plane h acts isometrically on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\hat f = { h^2}$\end{document } , then ( m , g ) is said to have hyperbolic symmetry .
in the case of spherical symmetry the existence of one compact cmc hypersurface implies that the whole spacetime can be covered by a cmc time coordinate that takes all real values [ 151 , 42 ] .
the existence of one compact cmc hypersurface in this case was proven by henkel   using the concept of prescribed mean curvature ( pmc ) foliation .
accordingly , this gives a complete picture in the spherically symmetric case regarding cmc foliations . in the case of areal coordinates ,
rein   has shown , under a size restriction on the initial data , that the past of an initial hyper - surface can be covered , and that the kretschmann scalar blows up .
hence , the initial singularity for the restricted data is both a crushing and a curvature singularity . in the future direction
it is shown that areal coordinates break down in finite time . in the case of plane and hyperbolic symmetry , global existence to the past was shown by rendall   in cmc time .
this implies that the past singularity is a crushing singularity since the mean curvature blows up at the singularity .
also in these cases rein showed   under a size restriction on the initial data , that global existence to the past in areal time and blow up of the kretschmann scalar curvature as the singularity is approached .
hence , the singularity is both a crushing and a curvature singularity in these cases too . in both of these works the question of global existence to the future was left open .
this gap was closed in , and global existence to the future was established in both cmc and areal time coordinates .
the global existence result for cmc time is a consequence of the global existence theorem in areal coordinates , together with a theorem by henkel   which shows that there exists at least one hypersurface with ( negative ) constant mean curvature .
in addition , the past direction is analyzed in   using areal coordinates , and global existence is shown without a size restriction on the data .
it is not concluded if the past singularity , without the smallness condition on the data , is a curvature singularity as well .
the issues discussed above have also been studied in the presence of a cosmological constant , cf .
[ 182 , 184 ] . in particular it is shown that in the spherically - symmetric case , if  < 0 , global existence to the future holds in areal time for some special classes of initial data , which is in contrast to the case with  = 0 . in this context
we also mention that surface symmetric spacetimes with vlasov matter and with a maxwell field have been investigated in . an interesting question , which essentially was left open in the studies mentioned above , is whether the areal time coordinate , which is positive by definition , takes all values in the range ( 0 ,  ) or only in ( r0 ,  ) for some positive r0 . it should here be pointed out that there is an example for vacuum spacetimes with t symmetry ( which includes the plane symmetric case ) where r0 > 0 .
this question was first resolved by weaver   for t symmetric spacetimes with vlasov matter .
her result shows that if spacetime contains vlasov matter ( f  0 ) then r0 = 0 .
smulevici   has recently shown , under the same assumption , that r0 = 0 also in the hyperbolic case .
smulevici also includes a cosmological constant  and shows that both the results , for plane ( or t ) symmetry and hyperbolic symmetry , are valid for   0 .
the important question of strong cosmic censorship for surface - symmetric spacetimes has recently been investigated by neat methods by dafermos and rendall [ 67 , 65 ] .
the standard strategy to show cosmic censorship is to either show causal geodesic completeness in case there are no singularities , or to show that some curvature invariant blows up along any incomplete causal geodesic . in both cases no causal geodesic can leave the maximal cauchy development in any extension if we assume that the extension is c. in [ 67 , 65 ] two alternative approaches are investigated .
the first method is independent of the matter model and exploits a rigidity property of cauchy horizons inherited from the killing fields .
the areal time described above is defined in terms of the killing fields and a consequence of the method by dafermos and rendall is that the killing fields extend continuously to a cauchy horizon , if one exists .
now , since global existence has been shown in areal time it follows that there can not be an extension of the maximal hyperbolic development to the future .
the second method is dependent on vlasov matter and the idea is to follow the trajectory of a particle , which crosses the cauchy horizon and shows that the conservation laws for the particle motion associated with the symmetries of the spacetime , such as the angular momentum , lead to a contradiction . in most of the cases considered in  
there is an assumption on the initial data for the vlasov equation , which implies that the data have non - compact support in the momentum space
the results of the studies [ 67 , 65 ] can be summarized as follows . for plane and hyperbolic symmetry strong cosmic censorship
the restriction that matter has non - compact support in the momentum space is here imposed except in the plane case with  = 0 . in the spherically - symmetric case cosmic censorship
the difficulties to show cosmic censorship in this case are related to possible formation of extremal schwarzschild - de - sitter - type black holes .
cosmic censorship in the past direction is also shown for all symmetry classes , and for all values of  , for a special class of anti - trapped initial data .
although the methods developed in [ 67 , 65 ] provide a lot of information on the asymptotic structure of the solutions , questions on geodesic completeness and curvature blow up are not answered . in a few cases , information on these issues has been obtained .
as mentioned above , blow up of the kretschmann scalar curvature has been shown for restricted initial data . in the case of hyperbolic symmetry causal future geodesic completeness has been established by rein   when the initial data are small .
the authors show global existence to the future in areal time , and in particular they show that the spacetimes are future geodesically complete .
it is shown that the de sitter solution acts as a model for the dynamics of the solutions by proving that the generalized kasner exponents tend to 1/3 as t   , which in the plane case is the de sitter solution .
the first study of spacetimes admitting a two - dimensional isometry group was carried out by rendall   in the case of local t symmetry . for a discussion of the possible topologies of these spacetimes we refer to the original paper . in the model case
the spacetime is topologically of the form   t , and to simplify our discussion later on we write down the metric in areal coordinates for this type of spacetime : 52\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { g = { e^{2(\eta - u)}}(- \alpha d{t^2 } + d{\theta ^2 } ) + { e^{- 2u}}{t^2}{{[dy + h\;d\theta + m\;dt]}^2}}\\ { \quad \quad + { e^{2u}}{{[dx + a\;dy + ( g + ah)\;d\theta + ( l + am)dt]}^2}.}\\ \end{array}$$\end{document } here the metric coefficients  , u ,  , a , h , l , and m depend on t and  and  , x , y  s. in   cmc coordinates are in fact considered rather than areal coordinates . under the hypothesis that there exists at least one cmc hypersurface , rendall proves for general initial data that the past of the given cmc hypersurface can be globally foliated by cmc hypersurfaces and that the mean curvature of these hypersurfaces blows up at the past singularity .
the result in   holds for vlasov matter and for matter described by a wave map . that the choice of matter model is important was shown in , where a non - global existence result for dust is given , which leads to examples of spacetimes   that are not covered by a cmc foliation
the plane case , where the symmetry group is three - dimensional , is one subcase and the form of the metric in areal coordinates is obtained by letting a = g = h = l = m = 0 and u = logt/2 in equation ( 52 ) .
another subcase , which still admits only two killing fields ( and which includes plane symmetry as a special case ) , is gowdy symmetry .
it is obtained by letting g = h = l = m = 0 in equation ( 52 ) . in   gowdy symmetric spacetimes with vlasov matter
are considered , and it is proven that the entire maximal globally hyperbolic spacetime can be foliated by constant areal time slices for general initial data .
the areal coordinates are used in a direct way for showing global existence to the future , whereas the analysis for the past direction is carried out in conformal coordinates .
these coordinates are not fixed to the geometry of spacetime and it is not clear that the entire past has been covered .
a chain of geometrical arguments then shows that areal coordinates indeed cover the entire spacetime .
the method in   was in turn inspired by the work   for vacuum spacetimes , where the idea of using conformal coordinates in the past direction was introduced .
as pointed out in , the result by henkel   guarantees the existence of one cmc hypersurface in the gowdy case and , together with the global areal foliation in , it follows that gowdy spacetimes with vlasov matter can be globally covered by cmc hypersurfaces as well .
the more general case of t symmetry was considered in , where global cmc and areal time foliations were established for general initial data . in these results ,
the question whether or not the areal time coordinate takes values in ( 0 ,  ) or in ( r0 ,  ) , r0 > 0 , was left open .
as we pointed out in section 4.2.1 , this issue was solved by weaver   for t symmetric spacetimes with the conclusion that r0 = 0 , if the density function f is not identically zero initially . in the case of t
symmetric spacetimes , with a positive cosmological constant , smulevici   has shown global existence in areal time with the property that t  ( 0 ,  ) .
the issue of strong cosmic censorship for t symmetric spacetimes has been studied by dafermos and rendall using the methods , which were developed in the surface symmetric case described above . in   strong cosmic censorship
is shown under the same restriction on the initial data that was imposed in the surface symmetric case , which implies that the data have non - compact support in the momentum variable .
their result has been extended to the case with a positive cosmological constant by smulevici .
the present cosmological observations indicate that the expansion of the universe is accelerating , and this has influenced theoretical studies in the field during the last decade .
one way to produce models with accelerated expansion is to choose a positive cosmological constant .
another way is to include a non - linear scalar field among the matter fields , and in this section we review the results for the einstein - vlasov system , where a linear or non - linear scalar field have been included into the model .
lee considers in   the case where a non - linear scalar field is coupled to vlasov matter .
the form of the energy momentum tensor then reads 53\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${t_{\alpha \beta } } = t_{\alpha \beta}^{{\rm{vlasov } } } + { \nabla _ \alpha}\phi { \nabla _ \beta}\phi - \left({{1 \over 2}{\nabla ^\gamma}\phi { \nabla _ \gamma}\phi + v(\phi ) } \right){g_{\alpha \beta}}.$$\end{document } here  is the scalar field and v is a potential , and the bianchi identities lead to the following equation for the scalar field : 54\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla ^\gamma}{\nabla _
\gamma}\phi = v^{\prime}(\phi){.}$$\end{document } under the assumption that v is a non - negative c function , global existence to the future is obtained , and if the potential is restricted to the form \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$v(\phi ) = { v_0}{e^{- c\phi}},$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0 < c < 4\sqrt \pi$\end{document } , then future geodesic completeness is proven . in   the einstein - vlasov system with a linear scalar field is analyzed in the case of plane , spherical , and hyperbolic symmetry . here
. a local existence theorem and a continuation criterion , involving bounds on derivatives of the scalar field in addition to a bound on the support of one of the moment variables , is proven . for the einstein scalar field system ,
i.e. , when f = 0 , the continuation criterion is shown to be satisfied in the future direction , and global existence follows in that case .
the work   extends the result in the plane and hyperbolic case to a global result in the future direction . in the plane case
it is also shown that the singularity is crushing and that the kretschmann scalar diverges uniformly as the singularity is approached . in standard cosmology
this is a strong assumption leading to severe restrictions of the possible geometries as well as of the topologies of the universe .
thus , it is natural to ask if small perturbations of an initial data set , which corresponds to an expanding model of the standard type , give rise to solutions that are similar globally to the future ? in a recent work , ringstrm   considers the einstein - vlasov system and he gives an affirmative answer to the stability question for some of the standard cosmologies .
the standard model of the universe is spatially homogeneous and isotropic , has flat spatial hypersurfaces of homogeneity , a positive cosmological constant and the matter content consists of a radiation fluid and dust .
hence , to investigate the question on stability it is natural to consider cosmological solutions with perfect fluid matter and a positive cosmological constant .
however , as is shown by ringstrm , the standard model can be well approximated by a solution of the einstein - vlasov system with a positive cosmological constant .
approximating dust with vlasov matter is straightforward , whereas approximating a radiation fluid is not . by choosing the initial support of the distribution function suitably , ringstrm shows that vlasov matter can be made to mimic a radiation fluid for a prescribed amount of time ; sooner or later the matter will behave like dust , but the time at which the approximation breaks down can be chosen to be large enough that the radiation is irrelevant to the future of that time in the standard picture .
the main results in   are stability of expanding , spatially compact , spatially locally homogeneous solutions to the einstein - vlasov system with a positive cosmological constant as well as a construction of solutions with arbitrary compact spatial topology . in other words
, the assumption of almost spatial homogeneity and isotropy does not seem to impose a restriction on the allowed spatial topologies .
let us mention here some related works although these do not concern the einstein - vlasov system .
ringstrm considers the case where the matter model is a non - linear scalar field in   and .
the background solutions , which ringstrm perturb and which are shown to be stable , have accelerated expansion . in   the expansion
the corresponding problem for a fluid has been treated in   and , and the newtonian case is investigated in   and   for vlasov and fluid matter respectively .
equilibrium states in galactic dynamics can be described as static or stationary solutions of the einstein - vlasov system , or of the vlasov - poisson system in the newtonian case .
here we consider the relativistic case and we refer to the excellent review paper   for the newtonian case .
first , we discuss spherically - symmetric solutions for which the structure is quite well understood . on the other hand ,
almost nothing is known about the stability of the spherically - symmetric static solutions of the einstein - vlasov system , which is in sharp contrast to the situation for the vlasov - poisson system . at the end of this section
let the metric have the form \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - { e^{2\mu ( r)}}d{t^2 } + { e^{2\lambda ( r)}}d{r^2 } + { r^2}(d{\theta ^2 } + { \sin ^2}\theta d{\varphi ^2}),$$\end{document } where r  0 ,   [ 0 ,  ] ,   [ 0 , 2 ] . as before ,
asymptotic flatness is expressed by the boundary conditions \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\underset{r \rightarrow \infty}{\lim }
\lambda ( r ) = \underset{r \rightarrow \infty}{\lim } \mu ( r ) = 0,$$\end{document } and a regular center requires \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\lambda ( 0 ) = 0.$$\end{document } following the notation in section 3.1 , the time - independent einstein - vlasov system reads 55\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${e^{\mu - \lambda}}{v \over { \sqrt { 1 + \vert v{\vert ^2 } } } } \cdot { \nabla _ x}f - \sqrt { 1 + \vert v{\vert ^2 } } { e^{\mu - \lambda}}{\mu _ r}{x \over r } \cdot { \nabla _ v}f = 0,$$\end{document } 
 56\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${e^{- 2\lambda}}(2r{\lambda _ r } - 1 ) + 1 = 8\pi { r^2}\rho,$$\end{document } 
 57\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${e^{- 2\lambda}}(2r{\mu _ r } + 1 ) - 1 = 8\pi { r^2}p.$$\end{document } recall that there is an additional equation ( 39 ) of second order , which contains the tangential pressure pt , but we leave it out since it follows from the equations above .
the matter quantities are defined as before : 58\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho ( x ) = \int\nolimits_{{{\mathbb r}^3 } } { \sqrt { 1 + \vert v{\vert ^2 } } f(x , v)dv,}$$\end{document } 
 59\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p(x ) = \int\nolimits_{{{\mathbb r}^3 } } { { { \left({{{x \cdot v } \over r } } \right)}^2}f(x , v){{dv } \over { \sqrt { 1 + \vert v{\vert ^2}}}}.}$$\end{document } the quantities \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$e : = { e^{\mu ( r)}}\sqrt { 1 + \vert v{\vert ^2 } } , \quad l = \vert x{\vert ^2}\vert v{\vert ^2 } - { ( x \cdot v)^2 } = \vert x \times v{\vert ^2},$$\end{document } are conserved along characteristics .
if we let 60\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f(x , v ) = \phi ( e , l),$$\end{document } for some function  , the vlasov equation is automatically satisfied .
a common assumption in the literature is to restrict the form of  to 61\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\phi ( e , l ) = \phi ( e)(l - { l_0 } ) _
+ ^l,$$\end{document } where l > 1/2 , l0  0 and x+ = max{x , 0}. if we furthermore assume that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\phi ( e ) = ( e - { e_0 } ) _ + ^k,\;k > - 1,$$\end{document } for some positive constant e0 , then we obtain the polytropic ansatz .
the case of isotropic pressure is obtained by letting l = 0 so that  only depends on e. in passing , we mention that for the vlasov - poisson system it has been shown   that every static spherically - symmetric solution must have the form f = (e , l ) . this is referred to as jeans theorem .
it was an open question for some time whether or not this was also true for the einstein - vlasov system .
this was settled in 1999 by schaeffer , who found solutions that do not have this particular form globally on phase space , and consequently , jeans theorem is not valid in the relativistic case .
however , almost all results on static solutions are based on this ansatz . by inserting the ansatz f(x ,
v ) = (e , l ) in the matter quantities  and p , a non - linear system for  and  is obtained , in which \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { e^{- 2\lambda}}(2r{\lambda _ r } - 1 ) + 1 = 8\pi { r^2}{g_\phi}(r,\mu),}\\ { { e^{- 2\lambda}}(2r{\mu _ r } + 1 ) - 1 = 8\pi { r^2}{h_\phi}(r,\mu),}\\ \end{array}$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { g_\phi}(r,\mu ) = { { 2\pi } \over { { r^2}}}\int\nolimits_1^\infty { \int\nolimits_0^{{r^2}({\epsilon ^2 } - 1 ) } { \phi ( { e^{\mu ( r)}}\epsilon , l){{{\epsilon ^2 } } \over { \sqrt { { \epsilon ^2 } - 1 - l/{r^2}}}}dld\epsilon}}}\\ { { h_\phi}(r,\mu ) = { { 2\pi } \over { { r^2}}}\int\nolimits_1^\infty { \int\nolimits_0^{{r^2}({\epsilon ^2 } - 1 ) } { \phi ( { e^{\mu ( r)}}\epsilon , l)\sqrt { { \epsilon ^2 } - 1 - l/{r^2 } } dld\epsilon.}}}\\ \end{array}$$\end{document } existence of solutions to this system was first proven in the case of isotropic pressure in , and extended to anisotropic pressure in . the main difficulty is to show that the solutions have finite adm mass and compact support .
the argument in these works to obtain a solution of compact support is to perturb a steady state of the vlasov - poisson system , which is known to have compact support .
two different types of solutions are constructed , those with a regular centre [ 144 , 134 ] , and those with a schwarzschild singularity in the centre . in ,
rein and rendall go beyond the polytropic ansatz and obtain steady states with compact support and finite mass under the assumption that  satisfies \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\phi ( e , l ) = c(e - { e_0 } ) _
+ ^{\delta + k}){l^l}\;{\rm{as}}\;e \rightarrow { e_0},$$\end{document } where k > 1 , l > 1/2 , k + l + 1/2 > 0 , k < l + 3/2 .
this result is obtained in a more direct way and is not based on the perturbation argument used in [ 144 , 134 ] .
their method is inspired by a work on stellar models by makino , in which he considers steady states of the euler - einstein system . in  
there is also a discussion about steady states that appear in the astrophysics literature , and it is shown that their result applies to most of these steady states . an alternative method to obtain steady states with finite radius and finite mass , which is based on a dynamical system analysis ,
all solutions described so far have the property that the support of  contains a ball about the centre . in   rein
showed that steady states also exist whose support is a finite , spherically - symmetric shell with a vacuum region in the center . in   it
a systematic study of the structure of spherically - symmetric static solutions was carried out mainly by numerical means in   and we now present the conclusions of this investigation . by prescribing the value (0 ) , the equations can be solved , but the resulting solution will in general not satisfy the boundary condition ( ) = 0 , but it will have some finite limit .
it is then possible to shift both the cut - off energy e0 and the solution by this limit to obtain a solution , which satisfies ( ) = 0 . a convenient way to handle the problem that e0 and (0 )
can not both be treated as free parameters is to use the ansatz \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f(x , v ) = ( 1 - e/{e_0 } ) _
+ ^l , k \geq 0,\;l > - 1/2,k < 3l + 7/2,$$\end{document } as in .
this gives an equation for  , which can be rewritten in terms of the function \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$y(r ) = { { { e^{\mu ( r ) } } } \over { { e_0}}}.$$\end{document } in this way the cut - of energy disappears as a free parameter of the problem and we thus have the four free parameters k , l , l0 and y(0 ) .
the structure of the static solutions obtained in   is as follows : if l0 = 0 the energy density can be strictly positive or vanish at r = 0 ( depending on l ) but it is always strictly positive sufficiently close to r = 0 . hence , the support of the matter is an interval [ 0 , r1 ] with r1 > 0 , and we call such states ball configurations .
if l0 > 0 the support is in an interval [ r0 , r1 ] , r0 > 0 , and we call such steady states for shells
. the value y(0 ) determines how compact or relativistic the steady state is , and the smaller values the more relativistic . for large values , recall y(0 )  1/e0 , a pure shell or a pure ball configuration is obtained , cf .
note that we depict the behavior of  but we remark that the pressure terms behave similarly but the amplitudes of p and pt can be very different , i.e. , the steady states can be highly anisotropic . 

figure 1a pure shell for moderate values of y(0 ) the solutions have a distinct inner peak and a tail - like outer peak , and by making y(0 ) smaller more peaks appear , cf .
figure 2multi - peaks of ball configurations , l0 = 0 multi - peaks of ball configurations , l0 = 0 in the case of shells there is a similar structure but in this case the peaks can either be separated by vacuum regions or by thin atmospheric regions as in the case of ball configurations .
an example with multi - peaks , where some of the peaks are separated by vacuum regions , is given in figure 3 . 

figure 3multi - peaks of a shell multi - peaks of a shell a different feature of the structure of static solutions is the issue of spirals . for a fixed ansatz of the density function
f , there is a one - parameter family of static solutions , which are parameterized by y(0 ) .
a natural question to ask is how the adm mass m and the radius of the support r change along such a family . by plotting for each y(0 ) the resulting values for r and m a curve
is obtained , which reflects how radius and mass are related along such a one - parameter family of steady states .
figure 4 . it is shown in   that in the isotropic case , where l = l0 = 0 the radius - mass curves always have a spiral form . 
 figure 4(r , m ) spiral for k = 0 , l = 10.5 , lo = 0 , 0.01  y(0 )  0.99 ( r , m ) spiral for k = 0 , l = 10.5 , lo = 0 , 0.01  y(0 )  0.99 another aspect of the structure of steady states investigated numerically in   concerns the buchdahl inequality .
if a steady state has support in [ r0 , r1 ] , then the adm mass m of the configuration is m = m(r1 ) , where the quasi local mass m(r ) is given by equation ( 50 ) in schwarzschild coordinates . in view of the schwarzschild metric ( 51 )
, schwarzschild asked already in 1916 the question : how large can 2m / r possibly be ?
he gave the answer   2m / r  8/9 in the special case of the schwarzschild interior solution , which has constant energy density and isotropic pressure . in 1959 buchdahl   extended his result to isotropic solutions for which the energy density is non - increasing outwards and he showed that also in this case 62\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{2 m } \over r } \leq { 8 \over 9}.$$\end{document } this is sometimes called the buchdahl inequality .
let us remark that the buchdahl inequality can obviously be written as m / r  4/9 , but since it is the quantity 2m / r , which appears in the schwarzschild metric ( 51 ) , it is common to keep the form of equation ( 62 ) .
a bound on 2m / r has an immediate observational consequence since it limits the possible gravitational red shift of a spherically - symmetric static object .
in particular , the overwhelming number of the steady states of the einstein - vlasov system have neither an isotropic pressure nor a non - increasing energy density , but nevertheless 2m / r is always found to be less than 8/9 in the numerical study . also for other matter models
as pointed out by guven and  murchadha , neither of the buchdahl assumptions hold in a simple soap bubble and they do not approximate any known topologically stable field configuration .
in addition , there are also several astrophysical models of stars , which are anisotropic .
lemaitre   proposed a model of an anisotropic star already in 1933 , and binney and tremaine   explicitly allow for an anisotropy coefficient .
hence , it is an important question to investigate bounds on 2m / r under less restrictive assumptions . in   it
is shown that for any static solution of the spherically - symmetric einstein equation , not necessarily of the einstein - vlasov system , for which p  0 , and 63\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p + 2{p_t } \leq \omega \rho,$$\end{document } the following inequality holds 64\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{2m(r ) } \over r } \leq { { { { ( 1 + 2\omega)}^2 } - 1 } \over { { { ( 1 + 2\omega)}^2}}}.$$\end{document } moreover , the inequality is sharp and sharpness is obtained uniquely by an infinitely thin shell solution .
note in particular that for vlasov matter  = 1 and that the right - hand side then equals 8/9 as in the buchdahl inequality .
an alternative proof was given in   and their method applies to a larger class of conditions on  , p and pt than the one given in equation ( 63 ) . on the other hand ,
the result in   is weaker than the result in   in the sense that the latter method implies that the steady state that saturates the inequality is unique ; it is an infinitely thin shell .
the studies [ 10 , 103 ] are of general character and in particular it is not shown that solutions exist to the coupled einstein - matter system , which can saturate the inequality .
for instance , it is natural to ask if there are solutions of the einstein - vlasov system , which have 2m / r arbitrarily close to 8/9 .
this question is given an affirmative answer in , where in particular it is shown that arbitrarily thin shells exist , which are regular solutions of the spherically - symmetric einstein - vlasov system .
using the strategy in   it follows that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\underset{r}{\sup}{{2m(r ) } \over r } \rightarrow { 8 \over 9},$$\end{document } in the limit when the shells become infinitely thin . the question of finding an upper bound on 2m / r can be extended to charged objects and to the case with a positive cosmological constant .
the spacetime outside a spherically - symmetric charged object is given by the reissner - nordstrm metric \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - ( 1 - { { 2 m } \over r } + { q \over { { r^2}}})d{t^2 } + { ( 1 - { { 2 m } \over r } + { q \over { { r^2}}})^{- 1}}d{r^2 } + { r^2}(d{\theta ^2 } + { \sin ^2}\theta d{\phi ^2}),$$\end{document } where q is the total charge of the object . the quantity \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$1 - { { 2 m } \over r } + { q \over { { r^2}}}$\end{document } is zero when \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${r _ \pm } = m \pm \sqrt { { m^2 } - { q^2}}$\end{document } , and r is called the inner and outer horizon respectively of a reissner - nordstrm black hole .
a buchdahl type inequality gives a lower bound of the area radius of a static object and this radius is thus often called the critical stability radius .
it is shown in   that a spherically - symmetric static solution of the einstein - maxwell system for which p  0 , p + 2pt   , and q
< m satisfy 65\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sqrt m \leq { { \sqrt r } \over 3 } + \sqrt { { r \over 9 } + { { { q^2 } } \over { 3r}}}.$$\end{document } note , in particular , that the inequality holds for solutions of the einstein - vlasov - maxwell system , since the conditions above are always satisfied in this case . this inequality ( 65 ) implies that the stability radius is outside the outer horizon of a reissner - nordstrm black hole . in   the relevance of an inequality of this kind on aspects in black - hole physics
an infinitely thin shell solution does saturate the inequality ( 65 ) , but numerical evidence is given in   that there is also another type of solution , which saturates the inequality for which the inner and outer horizon coincide .
the study in   is concerned with the non - charged situation when a positive cosmological constant  is included .
the following inequality is derived \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${m \over r } \leq { 2 \over 9 } - { { \lambda { r^2 } } \over 3 } + { 2 \over 9}\sqrt { 1 + 3\lambda { r^2}},$$\end{document } for solutions for which p  0 , p + 2pt   , and 0  r  1 . in this situation ,
an infinitely thin shell solution does not generally saturate the inequality but does so in the two degenerate situations r = 0 and r = 1 . in the latter case
there is a constant density solution , and the exterior spacetime is the nariai solution , which saturates the inequality and the saturating solution is thus non - unique . in this case , the cosmological horizon and the black hole horizon coincide , which is in analogy with the charged situation described above where the inner and outer horizons coincide when uniqueness is likely lost .
an important problem is the question of the stability of spherically - symmetric steady states . at present , there are almost no theoretical results on the stability of the steady states of the einstein - vlasov system .
wolansky   has applied the energy - casimir method and obtained some insights , but the theory is much less developed than in the vlasov - poisson case and the stability problem is essentially open .
the situation is very different for the vlasov - poisson system , and we refer to   for a review on the results in this case .
however , there are numerical studies [ 21 , 100 , 171 ] on the stability of spherically - symmetric steady states for the einstein - vlasov system .
the latter two studies concern isotropic steady states , whereas the first , in addition , treats anisotropic steady states . here
we present the conclusions of , emphasizing that these agree with the conclusions in [ 171 , 100 ] for isotropic states . to allow for trapped surfaces ,
maximal - areal coordinates are used , i.e. , the metric is written in the following form in   \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - ( { \alpha ^2 } - { a^2}{\beta ^2})d{t^2 } + 2{a^2}\beta dtdr + { a^2}d{r^2 } + { r^2}(d{\theta ^2 } + { \sin ^2}\theta \,d{\phi ^2}).$$\end{document } here the metric coefficients  ,  , and a depend on t   and r  0 ,  and a are positive , and the polar angles   [ 0 ,  ] and   [ 0 , 2 ] parameterize the unit sphere .
a maximal gauge condition is then imposed , which means that each hypersurface of constant t has vanishing mean curvature .
the boundary conditions , which guarantee asymptotic flatness and a regular center , are given by 66\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(t,0 ) = a(t,\infty ) = \alpha ( t,\infty ) = 1.$$\end{document } steady states are numerically constructed , and these are then perturbed in order to investigate the stability . more precisely , to construct the steady states the polytropic ansatz is used , cf .
section 5.1 , 67\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f(r , w , l ) = \phi ( e , l ) = ( { e_0 } - e ) _
+ ^l.$$\end{document } by specifying values on e0 , l0 and (0 ) steady states are numerically constructed .
the distribution function fs of the steady state is then multiplied by an amplitude a , so that a new , perturbed distribution function is obtained .
we remark that also other types of perturbations are analyzed in . for k and l
fixed each steady state is characterized by its central red shift zc and its fractional binding energy eb , which are defined by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${z_c } = { 1 \over { \alpha ( 0 ) } } - 1,\;{e_b } = { { { e_b } } \over { { m_0}}},\;{\rm{where}}\;{e_b } = { m_0 } - m.$$\end{document } here \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${m_0 } = 4{\pi ^2}\int\nolimits_0^\infty { \int\nolimits_{- \infty}^\infty { \int\nolimits_0^\infty { a(t , r)f(t , r , w , l ) } dl\;dw\;dr}}$$\end{document } is the total number of particles , which , since all particles have rest mass one , equals the rest mass of the system .
m is the adm mass given by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$m = \int\nolimits_0^\infty { \left({4\pi \rho ( t , r ) + { 3 \over 2}{\kappa ^2}(t , r ) } \right){r^2}dr,}$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$k = \beta / r\alpha$\end{document } both m0 and m are conserved quantities . the central redshift is the redshift of a photon emitted from the center and received at infinity , and the binding energy eb is the difference of the rest mass and the adm mass . in figure 5 and
figure 6 the relation between the fractional binding energy and the central redshift is given for two different cases . 

figure 5k = 0 , l = 0 , l0 = 0.1 
 figure 6k = 0 , l = 3/2 , l0 = 0.1 k = 0 , l = 0 , l0 = 0.1 k = 0 , l = 3/2 , l0 = 0.1 the relevance of these concepts for the stability properties of steady states was first discussed by zeldovich and podurets , who argued that it should be possible to diagnose the stability from binding energy considerations . zeldovich and novikov   then conjectured that the binding energy maximum along a steady state sequence signals the onset of instability .
varying the parameters k , l and l0 give rise to essentially the same tables , cf . . 

table 1k = 0 and l = 1/2 . 
 z 
 c 

 e 
 b 
 a <
1a > 10.210.032stablestable0.340.040stablestable0.390.040stablestable0.420.041stableunstable0.460.040stableunstable0.560.036stableunstable0.650.029stableunstable0.820.008stableunstable0.950.015unstableunstable1.200.078unstableunstable 
  if we first consider perturbations with a > 1 , it is found that steady states with small values on zc ( less than approximately 0.40 in this case ) are stable , i.e. , the perturbed solutions stay in a neighbourhood of the static solution .
a careful investigation of the perturbed solutions indicates that they oscillate in a periodic way . for larger values of zc
hence , for perturbations with a > 1 the value of zc alone seems to determine the stability features of the steady states .
plotting eb versus zc with higher resolution , cf . , gives support to the conjecture by novikov and zeldovich mentioned above that the maximum of eb along a sequence of steady states signals the onset of instability .
consider a steady state with eb > 0 and a perturbation with a < 1 but close to 1 so that the fractional binding energy remains positive .
the perturbed solution then drifts outwards , turns back and reimplodes , and comes close to its initial state , and then continues to expand and reimplode and thus oscillates , cf .
figure 7zc = 0.47 , eb = 0.04 , a = 0.99 , t = 90.0 zc = 0.47 , eb = 0.04 , a = 0.99 , t = 90.0 in   it is stated ( without proof ) that if eb > 0 the solution must ultimately reimplode and the simulations in   support that it is true . for negative values of eb , the solutions with a < 1 disperse to infinity .
a simple analytic argument is given in , which relates the question , whether a solution disperses or not .
it is shown that if a shell solution has an expanding vacuum region of radius r(t ) at the center with r(t )   for t   , i.e. , the solution disperses in a strong sense , then necessarily m0  m , i.e. , eb  0 . as we have seen above , a broad variety of static solutions of the einstein - vlasov system has been established , all of which share the restriction that they are spherically symmetric .
the recent investigation   removes this restriction and proves the existence of static solutions of the einstein - vlasov system , which are axially symmetric but not spherically symmetric . from the applications point of view
this symmetry assumption is more  realistic  than spherical symmetry , and from the mathematics point of view
the complexity of the einstein field equations increases drastically if one gives up spherical symmetry . before discussing this result ,
let us mention that similar results have been obtained for two other matter models . in the case of a perfect fluid , heilig showed the existence of axisymmetric stationary solutions in .
these solutions have non - zero angular momentum since static solutions are necessarily spherically symmetric . in this respect
the existence of static axisymmetric solutions of elastic matter , which are not spherically symmetric , was proven in .
let us now briefly discuss the method of proof in , which relies on an application of the implicit function theorem . also , the proofs in [ 92 , 1 , 2 ] make use of the implicit function theorem , but apart from this fact the methods are quite different .
the set - up of the problem in   follows the work of bardeen , where the metric is written in the form 68\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - { c^2}{e^{2\nu/{c^2}}}d{t^2 } + { e^{2\mu}}d{\rho ^2 } + { e^{2\mu}}d{z^2 } + { \rho ^2}{b^2}{e^{- 2\nu/{c^2}}}d{\varphi ^2}$$\end{document } for functions  , b ,  depending on  and z , where t   ,   [ 0 ,  [ , z   and   [ 0 , 2 ] . the killing vector fields t and  correspond to the stationarity and axial symmetry of the spacetime .
solutions are obtained by perturbing off spherically symmetric steady states of the vlasov - poisson system via the implicit function theorem and the reason for writing /c in the metric , instead of  , is that  converges to the newtonian potential un of the steady state in the limit c  . asymptotic flatness is expressed by the boundary conditions 69\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\underset{\vert ( \rho , z)\vert \rightarrow \infty}{\lim}\nu ( \rho , z ) = \underset{\vert ( \rho , z)\vert \rightarrow \infty}{\lim } \mu ( \rho , z ) = 0,\;\;\underset{\vert ( \rho , z)\vert \rightarrow \infty}{\lim } b(\rho , z ) = 1.$$\end{document } in addition the solutions are required to be locally flat at the axis of symmetry , which implies the condition 70\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\nu ( 0,z)/{c^2 } + \mu ( 0,z ) = \ln b(0,z),\;z \in { \mathbb r}.$$\end{document } let us now recall from section 5.1 the strategy to construct static solutions by using an ansatz of the form \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f(x , v ) = \phi ( e , l),$$\end{document } where e and l are conserved quantities along characteristics . due to
the symmetries of the metric ( 68 ) the following quantities are constant along geodesics : 71\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { e : = - g(\partial/\partial t,{p^a } ) = { c^2}{e^{2\nu/{c^2}}}{p^0}}\\ { \quad \ ; = { c^2}{e^{\nu/{c^2}}}\sqrt { 1 + { c^{- 2}}({e^{2\mu}}{{({p^1})}^2 } + { e^{2\mu}}{{({p^2})}^2 } + { \rho ^2}{b^2}{e^{- 2\nu/{c^2}}}{{({p^3})}^2})},}\\ \end{array}$$\end{document } 
 72\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$l : = - g(\partial/\partial \varphi , { p^a } ) = { \rho ^2}{b^2}{e^{- 2\nu/{c^2}}}{p^3}.$$\end{document } here p are the canonical momenta .
e can be thought of as a local or particle energy and l is the angular momentum of a particle with respect to the axis of symmetry . for a sufficiently regular 
the ansatz function f satisfies the vlasov equation and upon insertion of this ansatz into the definition of the energy momentum tensor ( 32 ) the latter becomes a functional ta = t( , b ,  ) of the unknown metric functions  , b , . it then remains to solve the einstein equations with this energy momentum tensor as right - hand side .
the newtonian limit of the einstein - vlasov system is the vlasov - poisson system and the strategy in   is to perturb off spherically symmetric steady states of the vlasov - poisson system via the implicit function theorem to obtain axisymmetric solutions .
indeed , the main result of   specifies conditions on the ansatz function  such that a two parameter (  and  ) family of axially - symmetric solutions of the einstein - vlasov system passes through the corresponding spherically symmetric , newtonian steady state , whose ansatz function we denote by . the parameter  = 1/c turns on general relativity and the parameter  turns on the dependence on l. since l is not invariant under arbitrary rotations about the origin the solution is not spherically symmetric if f depends on l. it should also be mentioned that although  is a priori small , which means that c is large , the scaling symmetry of the einstein - vlasov system can be used to obtain solutions corresponding to the physically correct value of c. the most striking condition on the ansatz function  , or rather on the ansatz function  of the corresponding vlasov - poisson system , needed to carry out the proof is that it must satisfy \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$6 + 4\pi { r^2}{a_n}(r ) > 0,\quad r \in [ 0,\infty [ , $ $ \end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${a_n}(r ) : = \int\nolimits_{{{\mathbb r}^3 } } { \phi ^{\prime } } \left({{1 \over 2}\vert v{\vert ^2 } + { u_n}(r ) } \right)dv.$$\end{document } an important argument in the proof is indeed to justify that there are steady states of the vlasov - poisson system satisfying this condition .
it is of course desirable to extend the result in   to stationary solutions with rotation .
moreover , the deviation from spherically symmetry of the solutions in   is small and an interesting open question is the existence of disk - like models of galaxies . in the vlasov - poisson case this has been shown in .
 the main purpose of this article is to provide a guide to theorems on global properties of solutions to the einstein - vlasov system . 
 this system couples einstein s equations to a kinetic matter model . 
 kinetic theory has been an important field of research during several decades in which the main focus has been on non - relativistic and special relativistic physics , i.e. , to model the dynamics of neutral gases , plasmas , and newtonian self - gravitating systems . in 1990 , rendall and rein initiated a mathematical study of the einstein - vlasov system . since then many theorems on global properties of solutions to this system have been established . 
 this paper gives introductions to kinetic theory in non - curved spacetimes and then the einstein - vlasov system is introduced . 
 we believe that a good understanding of kinetic theory in non - curved spacetimes is fundamental to a good comprehension of kinetic theory in general relativity . 
Introduction to Kinetic Theory
The Einstein-Vlasov System
The Asymptotically-Flat Cauchy Problem: Spherically-Symmetric Solutions
The Cosmological Cauchy Problem
Stationary Asymptotically-Flat Solutions