We study the fermionic King model which may provide a relevant model of dark matter halos. The exclusion constraint can be due to quantum mechanics (for fermions such as massive neutrinos) or to Lynden-Bell's statistics (for collisionless systems undergoing violent relaxation). This model has a finite mass. Furthermore, a statistical equilibrium state exists for all accessible values of energy. Dwarf and intermediate size halos are degenerate quantum objects stabilized against gravitational collapse by the Pauli exclusion principle. Large halos at sufficiently high energies are in a gaseous phase where quantum effects are negligible. They are stabilized by thermal motion. Below a critical energy they undergo gravitational collapse (gravothermal catastrophe). This may lead to the formation of a central black hole that does not affect the structure of the halo. This may also lead to the formation of a compact degenerate object surrounded by a hot massive atmosphere extending at large distances. We argue that large dark matter halos should not contain a degenerate nucleus (fermion ball) because these nucleus-halo structures are thermodynamically unstable. We compare the rotation curves of the classical King model to observations of large dark matter halos (Burkert profile). Because of collisions and evaporation, the central density increases while the slope of the halo density profile decreases until an instability takes place. We find that the observations are compatible with a King profile at, or close to, the point of marginal stability in the microcanonical ensemble. At that point, the King profile can be fitted by the modified Hubble profile that has a flat core and a halo in which the density decreases as $r^{-3}$. This is qualitatively similar to the Burkert profile.

We discuss the nature of phase transitions in the fermionic King model which describes tidally truncated quantum self-gravitating systems. This distribution function takes into account the escape of high energy particles and has a finite mass. On the other hand, the Pauli exclusion principle puts an upper bound on the phase space density of the system and stabilizes it against gravitational collapse. As a result, there exists a statistical equilibrium state for any accessible values of energy and temperature. We plot the caloric curves and investigate the nature of phase transitions as a function of the degeneracy parameter in both microcanonical and canonical ensembles. We consider stable and metastable states and emphasize the importance of the latter for systems with long-range interactions. Phase transitions can take place between a "gaseous" phase unaffected by quantum mechanics and a "condensed" phase dominated by quantum mechanics. The phase diagram exhibits two critical points, one in each ensemble, beyond which the phase transitions disappear. There also exist a region of negative specific heats and a situation of ensemble inequivalence for sufficiently large systems. We apply the fermionic King model to the case of dark matter halos made of massive neutrinos. The gaseous phase describes large halos and the condensed phase describes dwarf halos. Partially degenerate configurations describe intermediate size halos. We argue that large dark matter halos cannot harbor a fermion ball because these nucleus-halo configurations are thermodynamically unstable (saddle points of entropy). Large dark matter halos may rather contain a central black hole resulting from a dynamical instability of relativistic origin occurring during the gravothermal catastrophe.

We consider the possibility that dark matter halos are described by the Fermi-Dirac distribution at finite temperature. This is the case if dark matter is a self-gravitating quantum gas made of massive neutrinos at statistical equilibrium. This is also the case if dark matter can be treated as a self-gravitating collisionless gas experiencing Lynden-Bell's type of violent relaxation. In order to avoid the infinite mass problem and carry out a rigorous stability analysis, we consider the fermionic King model. In this paper, we study the non-degenerate limit leading to the classical King model. This model was initially introduced to describe globular clusters. We propose to apply it also to large dark matter halos where quantum effects are negligible. We determine the caloric curve and study the thermodynamical stability of the different configurations. Equilibrium states exist only above a critical energy $E_c$ in the microcanonical ensemble and only above a critical temperature $T_c$ in the canonical ensemble. For $E<E_c$, the system undergoes a gravothermal catastrophe and, for $T<T_c$, it undergoes an isothermal collapse. We compute the profiles of density, circular velocity, and velocity dispersion. We compare the prediction of the classical King model to the observations of large dark matter halos. Because of collisions and evaporation, the central density increases while the slope of the halo density profile decreases until an instability takes place. We show that large dark matter halos are relatively well-described by the King model at, or close to, the point of marginal microcanonical stability. At that point, the King model generates a density profile that can be approximated by the modified Hubble profile. This profile has a flat core and decreases as $r^{-3}$ at large distances, like the observational Burkert profile. Less steep halos are unstable.

In this work, we study the orbital stability of steady states and the existence of blow-up self-similar solutions to the so-called Vlasov-Manev (VM) system. This system is a kinetic model which has a similar Vlasov structure as the classical Vlasov-Poisson system, but is coupled to a potential in $-1/r- 1/r^2$ (Manev potential) instead of the usual gravitational potential in $-1/r$, and in particular the potential field does not satisfy a Poisson equation but a fractional-Laplacian equation. We first prove the orbital stability of the ground states type solutions which are constructed as minimizers of the Hamiltonian, following the classical strategy: compactness of the minimizing sequences and the rigidity of the flow. However, in driving this analysis, there are two mathematical obstacles: the first one is related to the possible blow-up of solutions to the VM system, which we overcome by imposing a sub-critical condition on the constraints of the variational problem. The second difficulty (and the most important) is related to the nature of the Euler-Lagrange equations (fractional-Laplacian equations) to which classical results for the Poisson equation do not extend. We overcome this difficulty by proving the uniqueness of the minimizer under equimeasurabilty constraints, using only the regularity of the potential and not the fractional-Laplacian Euler-Lagrange equations itself. In the second part of this work, we prove the existence of exact self-similar blow-up solutions to the Vlasov-Manev equation, with initial data arbitrarily close to ground states. This construction is based on a suitable variational problem with equimeasurability constraint.

We consider the three dimensional gravitational Vlasov Poisson system which is a canonical model in astrophysics to describe the dynamics of galactic clusters. A well known conjecture is the stability of spherical models which are nonincreasing radially symmetric steady states solutions. This conjecture was proved at the linear level by several authors in the continuation of the breakthrough work by Antonov in 1961. In a previous work (arXiv:0904.2443), we derived the stability of anisotropic models under {\it spherically symmetric perturbations} using fundamental monotonicity properties of the Hamiltonian under suitable generalized symmetric rearrangements first observed in the physics litterature. In this work, we show how this approach combined with a {\it new generalized} Antonov type coercivity property implies the orbital stability of spherical models under general perturbations.

We study the evaporation of stars from globular clusters using the simplified Chandrasekhar model. This is an analytically tractable model giving reasonable agreement with more sophisticated models that require complicated numerical integrations. In the Chandrasekhar model: (i) the stellar system is assumed to be infinite and homogeneous (ii) the evolution of the velocity distribution of stars f(v,t) is governed by a Fokker-Planck equation, the so-called Kramers-Chandrasekhar equation (iii) the velocities |v| that are above a threshold value R>0 (escape velocity) are not counted in the statistical distribution of the system. In fact, high velocity stars leave the system, due to free evaporation or to the attraction of a neighboring galaxy (tidal effects). Accordingly, the total mass and energy of the system decrease in time. If the star dynamics is described by the Kramers-Chandrasekhar equation, the mass decreases to zero exponentially rapidly. Our goal is to obtain non-perturbative analytical results that complement the seminal studies of Chandrasekhar, Michie and King valid for large times $t\to+\infty$ and large escape velocities $R\to +\infty$. In particular, we obtain an exact semi-explicit solution of the Kramers-Chandrasekhar equation with the absorbing boundary condition f(R,t)=0. We use it to obtain an explicit expression of the mass loss at any time t when $R\to +\infty$. We also derive an exact integral equation giving the exponential evaporation rate $\lambda(R)$, and the corresponding eigenfunction $f_{\lambda}(v)$, when $t\to +\infty$ for any sufficiently large value of the escape velocity R. For $R\to +\infty$, we obtain an explicit expression of the evaporation rate that refines the Chandrasekhar results.

We consider the three dimensional gravitational Vlasov Poisson system which describes the mechanical state of a stellar system subject to its own gravity. A well-known conjecture in astrophysics is that the steady state solutions which are nonincreasing functions of their microscopic energy are nonlinearly stable by the flow. This was proved at the linear level by several authors based on the pioneering work by Antonov in 1961. Since then, standard variational techniques based on concentration compactness methods as introduced by P.-L. Lions in 1983 have led to the nonlinear stability of subclasses of stationary solutions of ground state type. In this paper, inspired by pioneering works from the physics litterature (Lynden-Bell 94, Wiechen-Ziegler-Schindler MNRAS 88, Aly MNRAS 89), we use the monotonicity of the Hamiltonian under generalized symmetric rearrangement transformations to prove that non increasing steady solutions are local minimizer of the Hamiltonian under equimeasurable constraints, and extract compactness from suitable minimizing sequences. This implies the nonlinear stability of nonincreasing anisotropic steady states under radially symmetric perturbations.

We consider the three dimensional gravitational Vlasov-Poisson (GVP) system in both classical and relativistic cases. The classical problem is subcritical in the natural energy space and the stability of a large class of ground states has been derived by various authors. The relativistic problem is critical and displays finite time blow up solutions. Using standard concentration compactness techniques, we however show that the breaking of the scaling symmetry allows the existence of stable relativistic ground states. A new feature in our analysis which applies both to the classical and relativistic problem is that the orbital stability of the ground states does not rely as usual on an argument of uniqueness of suitable minimizers --which is mostly unknown-- but on strong rigidity properties of the transport flow, and this extends the class of minimizers for which orbital stability is now proved.