Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

We study the core mass -- halo mass relation of bosonic dark matter halos, in the form of self-gravitating Bose-Einstein condensates, harbouring a supermassive black hole. We use the ``velocity dispersion tracing'' relation according to which the velocity dispersion in the core $v_c^2\sim GM_c/R_c$ is of the same order as the velocity dispersion in the halo $v_h^2\sim GM_h/r_h$ (this relation can be justified from thermodynamical arguments) and the approximate analytical mass-radius relation of the quantum core in the presence of a central black hole obtained in our previous paper [P.H. Chavanis, Eur. Phys. J. Plus 134, 352 (2019)]. For a given minimum halo mass $(M_h)_{\rm min}\sim 10^8\, M_{\odot}$ determined by the observations, the only free parameter of our model is the scattering length $a_s$ of the bosons (their mass $m$ is then determined by the characteristics of the minimum halo). For noninteracting bosons and for bosons with a repulsive self-interaction, we find that the core mass $M_c$ increases with the halo mass $M_h$ and achieves a maximum value $(M_c)_{\rm max}$ at some halo mass $(M_h)_{*}$ before decreasing. The whole series of equilibria is stable. For bosons with an attractive self-interaction, we find that the core mass achieves a maximum value $(M_c)_{\rm max}$ at some halo mass $(M_h)_{*}$ before decreasing. The series of equilibria becomes unstable above a maximum halo mass $(M_h)_{\rm max}\ge (M_h)_{*}$. In the absence of black hole $(M_h)_{\rm max}=(M_h)_{*}$. At that point, the quantum core (similar to a dilute axion star) collapses. We perform a similar study for fermionic dark matter halos. We find that they behave similarly to bosonic dark matter halos with a repulsive self-interaction, the Pauli principle for fermions playing the role of the repulsive self-interaction for bosons.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.

We consider the possibility that dark matter halos are made of quantum particles such as fermions or bosons in the form of Bose-Einstein condensates. In that case, they generically have a "core-halo" structure with a quantum core that depends on the type of particle considered and a halo that is relatively independent of the dark matter particle and that is similar to the NFW profile of cold dark matter. We model the halo by an isothermal gas with an effective temperature $T$. We then derive the core mass -- halo mass relation $M_c(M_v)$ of dark matter halos from an effective thermodynamical model by extremizing the free energy $F(M_c)$ with respect to the core mass $M_c$. We obtain a general relation that is equivalent to the "velocity dispersion tracing" relation according to which the velocity dispersion in the core $v_c^2\sim GM_c/R_c$ is of the same order as the velocity dispersion in the halo $v_v^2\sim GM_v/r_v$. We provide therefore a justification of this relation from thermodynamical arguments. In the case of fermions, we obtain a relation $M_c\propto M_v^{1/2}$ that agrees with the relation found numerically by Ruffini {\it et al.}. In the case of noninteracting bosons, we obtain a relation $M_c\propto M_v^{1/3}$ that agrees with the relation found numerically by Schive {\it et al.}. In the case of bosons with a repulsive self-interaction in the Thomas-Fermi limit, we predict a relation $M_c\propto M_v^{2/3}$ that still has to be confirmed numerically. We also obtain a general approximate core mass -- halo mass relation that is valid for bosons with arbitrary repulsive or attractive self-interaction. For an attractive self-interaction, we determine the maximum halo mass that can harbor a stable quantum core (dilute axion "star").

Supernovae explosions of massive stars are nowadays believed to result from a two-step process, with an initial gravitational core collapse followed by an expansion of matter after a bouncing on the core. This scenario meets several difficulties. We show that it is not the only possible one: a simple model based on fluid mechanics and stability properties of the equilibrium state shows that one can have also a simultaneous inward/outward motion in the early stage of the instability of the supernova. This shows up in the slow sweeping across a saddle-center bifurcation found when considering equilibrium states associated to the constraint of energy conservation. We first discuss the weakly nonlinear regime in terms of a Painlev\'e I equation. We then show that the strongly nonlinear regime displays a self-similar behavior of the core collapse. Finally, the expansion of the remnants is revisited as an isentropic process leading to shocks formation.

The present understanding of supernova explosion of massive stars as a two-step process, with an initial gravitational collapse toward the center of the star followed by an expansion of matter after a bouncing on the core, meets several difficulties. We show that it is not the only possible one: a simple model based on fluid mechanics, catastrophe theory, and stability properties of the equilibrium state shows that one can have also a simultaneous inward/outward motion in the early stage of the instability of the supernova described by a dynamical saddle-center bifurcation. The existence of this simultaneous inward/outward motion is sensitive to the model in such systems with long-range interactions. If a constant temperature is assumed (canonical ensemble), an overall inward motion occurs, but if one imposes with the same equation of state the constraint of energy conservation (microcanonical ensemble) there is an inward velocity field near the center of the star together with an outward velocity field in the rest of the star. We discuss the expansion stage of the remnants away from the collapsed core, and propose a new explanation for the formation of shock waves in the ejecta which differs from the usual Sedov-Taylor self-similar description.

We describe microcanonical phase transitions and instabilities of the ideal Fermi gas in general relativity at nonzero temperature confined in the interior of a spherical shell. The thermodynamic behaviour is governed by the compactness of rest mass, namely of the total rest mass over radius of the system. For a fixed value of rest mass compactness, we study the caloric curves as a function of the size of the spherical box. At low compactness values, low energies and for sufficiently big systems the system is subject to a gravothermal catastrophe, which cannot be halted by quantum degeneracy pressure, and the system collapses towards the formation of a black hole. For small systems, there appears no instability at low energies. For intermediate sizes, between two marginal values, gravothermal catastrophe is halted and a microcanonical phase transition occurs from a gaseous phase to a condensed phase with a nearly degenerate core. The system is subject to a relativistic instability at low energy, when the core gets sufficiently condensed above the Oppenheimer-Volkoff limit. For sufficiently high values of rest mass compactness the microcanonical phase transitions are suppressed. They are replaced either by an Antonov type gravothermal catastrophe for sufficiently big systems or by stable equilibria for small systems. At high energies the system is subject to the `relativistic gravothermal instability', identified by Roupas in [1], for all values of compactness and any size.

In the vicinity of a massive black hole, stars move on precessing Keplerian orbits. The mutual stochastic gravitational torques between the stellar orbits drive a rapid reorientation of their orbital planes, through a process called vector resonant relaxation. We derive, from first principles, the correlation of the potential fluctuations in such a system, and the statistical properties of random walks undergone by the stellar orbital orientations. We compare this new analytical approach with effective $N$-body simulations. We also provide a simple scheme to generate the random walk of a test star's orbital orientation using a stochastic equation of motion. We finally present quantitative estimations of this process for a nuclear stellar cluster such as the one of the Milky Way.

We develop a model of Bose-Einstein condensate dark matter halos with a solitonic core and an isothermal atmosphere based on a generalized Gross-Pitaevskii-Poisson equation [P.H. Chavanis, Eur. Phys. J. Plus 132, 248 (2017)]. This equation provides a heuristic coarse-grained parametrization of the ordinary Gross-Pitaevskii-Poisson equation accounting for violent relaxation and gravitational cooling. It involves a cubic nonlinearity taking into account the self-interaction of the bosons, a logarithmic nonlinearity associated with an effective temperature, and a source of dissipation. It leads to superfluid dark matter halos with a core-halo structure. The quantum potential or the self-interaction of the bosons generates a solitonic core that solves the cusp problem of the cold dark matter model. The logarithmic nonlinearity generates an isothermal atmosphere accounting for the flat rotation curves of the galaxies. The dissipation ensures that the system relaxes towards an equilibrium configuration. In the Thomas-Fermi approximation, the dark matter halo is equivalent to a barotropic gas with an equation of state $P=2\pi a_s\hbar^2\rho^2/m^3+\rho k_B T/m$, where $a_s$ is the scattering length of the bosons and $m$ is their individual mass. We numerically solve the equation of hydrostatic equilibrium and determine the corresponding density profiles and rotation curves. Our model has no free parameter so it is completely predictive. Extension of this model to noninteracting bosons and fermions will be presented in forthcoming papers.

The Schr\"odinger-Poisson equations describe the behavior of a superfluid Bose-Einstein condensate under self-gravity with a 3D wave function. As $\hbar/m\to 0$, $m$ being the boson mass, the equations have been postulated to approximate the collisionless Vlasov-Poisson equations also known as the collisionless Boltzmann-Poisson equations. The latter describe collisionless matter with a 6D classical distribution function. We investigate the nature of this correspondence with a suite of numerical test problems in 1D, 2D, and 3D along with analytic treatments when possible. We demonstrate that, while the density field of the superfluid always shows order unity oscillations as $\hbar/m\to 0$ due to interference and the uncertainty principle, the potential field converges to the classical answer as $(\hbar/m)^{2}$. Thus, any dynamics coupled to the superfluid potential is expected to recover the classical collisionless limit as $\hbar/m\to 0$. The quantum superfluid is able to capture rich phenomena such as multiple phase-sheets, shell-crossings, and warm distributions. Additionally, the quantum pressure tensor acts as a regularizer of caustics and singularities in classical solutions. This suggests the exciting prospect of using the Schr\"odinger-Poisson equations as a low-memory method for approximating the high-dimensional evolution of the Vlasov-Poisson equations. As a particular example we consider dark matter composed of ultra-light axions, which in the classical limit ($\hbar/m\to 0$) is expected to manifest itself as collisionless cold dark matter.

We present a derivation of the kinetic equation describing the secular evolution of spatially inhomogeneous systems with long-range interactions, the so-called inhomogeneous Landau equation, by relying on a functional integral formalism. We start from the BBGKY hierarchy derived from the Liouville equation. At the order ${1/N}$, where $N$ is the number of particles, the evolution of the system is characterised by its 1-body distribution function and its 2-body correlation function. Introducing associated auxiliary fields, the evolution of these quantities may be rewritten as a traditional functional integral. By functionally integrating over the 2-body autocorrelation, one obtains a new constraint connecting the 1-body DF and the auxiliary fields. When inverted, this constraint allows us to obtain the closed non-linear kinetic equation satisfied by the 1-body distribution function. This derivation provides an alternative to previous methods, either based on the direct resolution of the truncated BBGKY hierarchy or on the Klimontovich equation. It may turn out to be fruitful to derive more accurate kinetic equations, e.g., accounting for collective effects, or higher order correlation terms.

We study the nature of phase transitions between dilute and dense axion stars interpreted as self-gravitating Bose-Einstein condensates. We develop a Newtonian model based on the Gross-Pitaevskii-Poisson equations for a complex scalar field with a self-interaction potential $V(|\psi|^2)$ involving an attractive $|\psi|^4$ term and a repulsive $|\psi|^6$ term. Using a Gaussian ansatz for the wave function, we analytically obtain the mass-radius relation of dilute and dense axion stars for arbitrary values of the self-interaction parameter $\lambda\le 0$. We show the existence of a critical point $|\lambda|_c\sim (m/M_P)^2$ above which a first order phase transition takes place. We qualitatively estimate general relativistic corrections on the mass-radius relation of axion stars. For weak self-interactions $|\lambda|<|\lambda|_c$, a system of self-gravitating axions forms a stable dilute axion star below a general relativistic maximum mass $M_{\rm max,GR}^{\rm dilute}\sim M_P^2/m$ and collapses into a black hole above that mass. For strong self-interactions $|\lambda|>|\lambda|_c$, a system of self-gravitating axions forms a stable dilute axion star below a Newtonian maximum mass $M_{\rm max,N}^{\rm dilute}=5.073 M_P/\sqrt{|\lambda|}$, collapses into a dense axion star above that mass, and collapses into a black hole above a general relativistic maximum mass $M_{\rm max,GR}^{\rm dense}\sim \sqrt{|\lambda|}M_P^3/m^2$. Dense axion stars explode below a Newtonian minimum mass $M_{\rm min,N}^{\rm dense}\sim m/\sqrt{|\lambda|}$ and form dilute axion stars of large size or disperse away. We determine the phase diagram of self-gravitating axions and show the existence of a triple point $(|\lambda|_*,M_*/(M_P^2/m))$ separating dilute axion stars, dense axion stars, and black holes. We make numerical applications for QCD axions and ultralight axions.

The secular thickening of a discrete self-gravitating galactic disc is investigated using the inhomogeneous multi-component Balescu-Lenard equation. The thick WKB limit for the diffusion and drift coefficients is found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields a simple double quadrature for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc induced by the effects of a finite number of particles. When applied to a tepid stable tapered disc, the Balescu-Lenard formalism predicts the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but over-estimates the timescale involved in their appearance. Swing amplication is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. The joint evolution of a population of giant molecular clouds within the disc may accelerate the secular disc's thickening induced by finite${-N}$ effects, but the observed number of clouds in the Milky Way does not seem to be sufficient to explain its thick disc.

The secular thickening of a self-gravitating stellar galactic disc is investigated using the dressed collisionless Fokker-Planck equation and the inhomogeneous multicomponent Balescu-Lenard equation. The thick WKB limits for the diffusion fluxes are found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields simple quadratures for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc, induced by fluctuations either external or due to the finite number of particles. These thick limits also offer a consistent derivation of a thick disc Toomre parameter, which is shown to be exponentially boosted by the ratio of the vertical to radial scale heights. Dressed potential fluctuations within the disc statistically induce a vertical bending of a subset of resonant orbits, triggering the corresponding increase in vertical velocity dispersion. When applied to a tepid stable tapered disc perturbed by shot noise, these two frameworks reproduce qualitatively the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but overestimates the time-scale involved in their appearance. Swing amplification is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. Other sources of thickening are also investigated, such as fading sequences of slowing bars, or the joint evolution of a population of giant molecular clouds within the disc.

General self-consistent expressions for the coefficients of diffusion and dynamical friction in a stable, bound, multicomponent self-gravitating and inhomogeneous system are derived. They account for the detailed dynamics of the colliding particles and their self-consistent dressing by collective gravitational interactions. The associated Fokker-Planck equation is shown to be fully consistent with the corresponding inhomogeneous Balescu-Lenard equation and, in the weak self-gravitating limit, to the inhomogeneous Landau equation. Hence it provides an alternative derivation to both and demonstrates their equivalence. The corresponding stochastic Langevin equations are presented: they can be a practical alternative to numerically solving the inhomogeneous Fokker-Planck and Balescu-Lenard equations. The present formalism allows for a self-consistent description of the secular evolution of different populations covering a spectrum of masses, with a proper accounting of the induced secular mass segregation, which should be of interest to various astrophysical contexts, from galactic centers to protostellar discs.

We perform a detailed comparison between the Logotropic model [P.H. Chavanis, Eur. Phys. J. Plus 130 (2015) 130] and the $\Lambda$CDM model. These two models behave similarly at large (cosmological) scales up to the present. Differences will appear only in the far future, in about $25\, {\rm Gyrs}$, when the Logotropic Universe becomes phantom while the $\Lambda$CDM Universe enters in the de Sitter era. However, the Logotropic model differs from the $\Lambda$CDM model at small (galactic) scales, where the latter encounters serious problems. Having a nonvanishing pressure, the Logotropic model can solve the cusp problem and the missing satellite problem of the $\Lambda$CDM model. In addition, it leads to dark matter halos with a constant surface density $\Sigma_0=\rho_0 r_h$, and can explain its observed value $\Sigma_0=141 \, M_{\odot}/{\rm pc}^2$ without adjustable parameter. This makes the logotropic model rather unique among all the models attempting to unify dark matter and dark energy. In this paper, we compare the Logotropic and $\Lambda$CDM models at the cosmological scale where they are very close to each other in order to determine quantitatively how much they differ. This comparison is facilitated by the fact that these models depend on only two parameters, the Hubble constant $H_0$ and the present fraction of dark matter $\Omega_{\rm m0}$. Using the latest observational data from Planck 2015+Lensing+BAO+JLA+HST, we find that the best fit values of $H_0$ and $\Omega_{\rm m0}$ are $H_0=68.30\, {\rm km}\, {\rm s}^{-1}\,{\rm Mpc}^{-1}$ and $\Omega_{\rm m0}=0.3014$ for the Logotropic model, and $H_0=68.02\, {\rm km}\, {\rm s}^{-1}\,{\rm Mpc}^{-1}$ and $\Omega_{\rm m0}=0.3049$ for the $\Lambda$CDM model. The difference between the two models appears at the percent level.

The drift and diffusion coefficients of the inhomogeneous multi-mass degenerate Landau equation are computed to describe the self-induced resonant relaxation of a discrete self-gravitating quasi-Keplerian razor-thin axisymmetric disc orbiting a massive black hole while relying on Gauss' method. For a disc-like configuration in our Galactic centre, secular diffusion induces an adiabatic distortion of orbits. When considering a disc composed of multiple masses similarly distributed, the population of lighter stars will gain eccentricity, driving it closer to the central black hole provided the distribution function increases with angular momentum. The quenching of the diffusion of a test star in the vicinity of the black hole due to the divergence of the relativistic precessions (the "Schwarzschild barrier") is correctly recovered by the kinetic equation. The dual stochastic Langevin formulation yields consistent results and provides a versatile framework in which to incorporate other stochastic processes.

We study the collapse of a self-gravitating Bose-Einstein condensate with attractive self-interaction. Equilibrium states in which the gravitational attraction and the attraction due to the self-interaction are counterbalanced by the quantum pressure exist only below a maximum mass $M_{\rm max}=1.012\hbar/\sqrt{Gm|a_s|}$ where $a_s<0$ is the scattering length of the bosons and $m$ is their mass. For $M>M_{\rm max}$ the system is expected to collapse and form a black hole. We study the collapse dynamics by making a Gaussian ansatz for the wave function. We find that the collapse time scales as $(M/M_{\rm max}-1)^{-1/4}$ for $M\rightarrow M_{\rm max}^+$ and as $M^{-1/2}$ for $M\gg M_{\rm max}$. We apply our results to standard axions with mass $m=10^{-4}\, {\rm eV}/c^2$ and scattering length $a_s=-5.8\times 10^{-53}\, {\rm m}$ for which $M_{\rm max}=6.5\times 10^{-14}M_{\odot}$ and $R=3.3\times 10^{-4}\, R_{\odot}$. We confirm our previous claim that bosons with attractive self-interaction, such as standard axions, may form low mass stars but cannot form dark matter halos of relevant mass and size. These mini axions stars could be the constituents of dark matter. They can collapse into mini black holes of mass $\sim 10^{-14}\, M_{\odot}$ in a few hours. In that case, dark matter halos would be made of mini black holes. We also apply our results to ultralight axions with mass $m=1.93\times 10^{-20}\, {\rm eV}/c^2$ and scattering length $a_s=-8.29\times 10^{-60}\, {\rm fm}$ for which $M_{\rm max}=0.39\times 10^6\, M_{\odot}$ and $R=33\, {\rm pc}$. These ultralight axions could cluster into dark matter halos. Axionic dark matter halos with attractive self-interaction can collapse into supermassive black holes of mass $\sim 10^{6}\, M_{\odot}$ in about one million years.

Using a generalization of the Madelung transformation, we derive the hydrodynamic representation of the Klein-Gordon-Einstein equations in the weak field limit. We consider a complex self-interacting scalar field with a $\lambda|\varphi|^4$ potential. We study the evolution of the homogeneous background in the fluid representation and derive the linearized equations describing the evolution of small perturbations in a static and in an expanding universe. We compare the results with simplified models in which the gravitational potential is introduced by hand in the Klein-Gordon equation, and assumed to satisfy a (generalized) Poisson equation. We study the evolution of the perturbations in the matter era using the nonrelativistic limit of our formalism. Perturbations whose wavelength is below the Jeans length oscillate in time while pertubations whose wavelength is above the Jeans length grow linearly with the scale factor as in the cold dark matter model. The growth of perturbations in the scalar field model is substantially faster than in the cold dark matter model. When the wavelength of the pertubations approaches the cosmological horizon (Hubble length), a relativistic treatment is mandatory. In that case, we find that relativistic effects attenuate or even prevent the growth of pertubations. This paper exposes the general formalism and provides illustrations in simple cases. Other applications of our formalism will be considered in companion papers.

We develop a cosmological model based on a quadratic equation of state p/c^2=-(\alpha+1){\rho^2}/{\rho_P}+\alpha\rho-(\alpha+1)\rho_{\Lambda} (where \rho_P is the Planck density and \rho_{\Lambda} the cosmological density) "unifying" vacuum energy and dark energy in the spirit of a generalized Chaplygin gas model. For $\rho\rightarrow \rho_P$, it reduces to p=-\rho c^2 leading to a phase of early accelerated expansion (early inflation) with a constant density equal to the Planck density \rho_P (vacuum energy). For $\rho_{\Lambda}\ll\rho\ll \rho_P$, we recover the standard linear equation of state p=\alpha \rho c^2 describing radiation (\alpha=1/3) or pressureless matter (\alpha=0) and leading to an intermediate phase of decelerating expansion. For $\rho\rightarrow \rho_{\Lambda}$, we get p=-\rho c^2 leading to a phase of late accelerated expansion (late inflation) with a constant density equal to the cosmological density \rho_{\Lambda} (dark energy). We show a nice symmetry between the early universe (vacuum energy + \alpha-fluid) and the late universe (\alpha-fluid + dark energy). In our model, they are described by two polytropic equations of state with index n=+1 and n=-1 respectively. Furthermore, the Planck density \rho_P in the early universe plays a role similar to the cosmological density \rho_{\Lambda} in the late universe. They represent fundamental upper and lower density bounds differing by 122 orders of magnitude. This quadratic equation of state leads to a fully analytical model describing the evolution of the universe from the early inflation (Planck era) to the late accelerated expansion (de Sitter era). These two phases are bridged by a decelerating algebraic expansion (\alpha-era). This model does not present any singularity at t=0 and exists eternally in the past. It admits a scalar field interpretation based on a quintessence field or a tachyon field.

We propose a heuristic unification of dark matter and dark energy in terms of a single dark fluid with a logotropic equation of state $P=A\ln(\rho/\rho_P)$, where $\rho$ is the rest-mass density, $\rho_P$ is the Planck density, and $A$ is the logotropic temperature. The energy density $\epsilon$ is the sum of a rest-mass energy term $\rho c^2$ mimicking dark matter and an internal energy term $u(\rho)=-P(\rho)-A$ mimicking dark energy. The logotropic temperature is approximately given by $A \simeq \rho_{\Lambda}c^2/\ln(\rho_P/\rho_{\Lambda})\simeq\rho_{\Lambda}c^2/[123 \ln(10)]$, where $\rho_{\Lambda}$ is the cosmological density. More precisely, we obtain $A=2.13\times 10^{-9} \, {\rm g}\, {\rm m}^{-1}\, {\rm s}^{-2}$ that we interpret as a fundamental constant. At the cosmological scale, this model fullfills the same observational constraints as the $\Lambda$CDM model. However, it has a nonzero velocity of sound and a nonzero Jeans length which, at the beginning of the matter era, is about $\lambda_J=40.4\, {\rm pc}$, in agreement with the minimum size of the dark matter halos observed in the universe. At the galactic scale, the logotropic pressure balances gravitational attraction and solves the cusp problem and the missing satellite problem. The logotropic equation of state generates a universal rotation curve that agrees with the empirical Burkert profile of dark matter halos up to the halo radius. In addition, it implies that all the dark matter halos have the same surface density $\Sigma_0=\rho_0 r_h=141\, M_{\odot}/{\rm pc}^2$ and that the mass of dwarf galaxies enclosed within a sphere of fixed radius $r_{u}=300\, {\rm pc}$ has the same value $M_{300}=1.93\times 10^{7}\, M_{\odot}$, in remarkable agreement with the observations.

The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is investigated using the inhomogeneous Balescu-Lenard equation expressed in terms of angle-action variables. The matrix method is implemented numerically in order to model the induced gravitational polarization. Special care is taken to account for the amplification of potential fluctuations of mutually resonant orbits and the unwinding of the induced swing amplified transients. Quantitative comparisons with ${N-}$body simulations yield consistent scalings with the number of particles and with the self-gravity of the disc: the fewer particles and the colder the disc, the faster the secular evolution. Secular evolution is driven by resonances, but does not depend on the initial phases of the disc. For a Mestel disc with ${Q \sim 1.5}$, the polarization cloud around each star boosts up its secular effect by a factor of the order of a thousand or more, promoting accordingly the dynamical relevance of self-induced collisional secular evolution. The position and shape of the induced resonant ridge are found to be in very good agreement with the prediction of the Balescu-Lenard equation, which scales with the square of the susceptibility of the disc. In astrophysics, the inhomogeneous Balescu-Lenard equation may describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic centre. It could be used as a valuable check of the accuracy of ${N-}$body integrators over secular timescales.

We consider the possibility that the universe is made of a single dark fluid described by a logotropic equation of state $P=A\ln(\rho/\rho_*)$, where $\rho$ is the rest-mass density, $\rho_*$ is a reference density, and $A$ is the logotropic temperature. The energy density $\epsilon$ is the sum of two terms: a rest-mass energy term $\rho c^2$ that mimics dark matter and an internal energy term $u(\rho)=-P(\rho)-A$ that mimics dark energy. This decomposition leads to a natural, and physical, unification of dark matter and dark energy, and elucidates their mysterious nature. The logotropic model depends on a single parameter $B=A/\rho_{\Lambda}c^2$ where $\rho_{\Lambda}$ is the cosmological density. For $B=0$, we recover the $\Lambda$CDM model. Using cosmological constraints, we find that $0\le B\le 0.09425$. We consider the possibility that dark matter halos are described by the same logotropic equation of state. When $B>0$, pressure gradients prevent gravitational collapse and provide halo density cores instead of cuspy density profiles, in agreement with the observations. The universal rotation curve of logotropic dark matter halos is consistent with the observational Burkert profile up to the halo radius. Interestingly, if we assume that all the dark matter halos have the same logotropic temperature $B$, we find that their surface density $\Sigma=\rho_0 r_h$ is constant. This result is in agreement with the observations where it is found that $\Sigma_0=141\, M_{\odot}/{\rm pc}^2$ for dark matter halos differing by several orders of magnitude in size. Using this observational result, we obtain $B=3.53\times 10^{-3}$. Assuming that $\rho_*=\rho_P$, where $\rho_P$ is the Planck density, we predict $B=3.53\times 10^{-3}$, in perfect agreement with the value obtained from the observations.

The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is described using the inhomogeneous Balescu-Lenard equation. Assuming that only tightly wound transient spirals are present in the disc, a WKB approximation provides a simple and tractable quadrature for the corresponding drift and diffusion coefficients. It provides insight into the physical processes at work during the secular diffusion of a self-gravitating discrete disc and makes quantitative predictions on the initial variations of the distribution function in action space. When applied to the secular evolution of an isolated stationary self-gravitating Mestel disc, this formalism predicts initially the importance of the corotation resonance in the inner regions of the disc leading to a regime involving radial migration and heating. It predicts in particular the formation of a "ridge like" feature in action space, in agreement with simulations, but over-estimates the timescale involved in its appearance. Swing amplification is likely to resolve this discrepancy. In astrophysics, the inhomogeneous Balescu-Lenard equation and its WKB limit may also describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic center.

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.

The supernovae are typical complex phenomena in fluid mechanics with very different time scales. We describe them in the light of catastrophe theory, assuming that successive equilibria between pressure and gravity present a saddle-node bifurcation. In the early stage we show that the loss of equilibrium may be described by a generic equation of the Painlev\'e I form. In the final stage of the collapse, just before the divergence of the central density, we show that the existence of a self-similar collapsing solution compatible with the numerical observations imposes that the gravity forces are stronger than the pressure ones. This situation differs drastically in its principle from the one generally admitted where pressure and gravity forces are assumed to be of the same order. Our new self-similar solution (based on the hypothesis of dominant gravity forces) which matches the smooth solution of the outer core solution, agrees globally well with our numerical results. Whereas some differences with the earlier self-similar solutions are minor, others are very important. For example, we find that the velocity field becomes singular at the collapse time, diverging at the center, and decreasing slowly outside the core, whereas previous works described a finite velocity field in the core which tends to a supersonic constant value at large distances. This discrepancy should be important for explaining the emission of remnants in the post-collapse regime. Finally we describe the post-collapse dynamics, when mass begins to accumulate in the center, also within the hypothesis that gravity forces are dominant.

We review and complete the kinetic theory of spatially inhomogeneous stellar systems when collective effects (dressing of the stars by their polarization cloud) are neglected. We start from the BBGKY hierarchy issued from the Liouville equation and consider an expansion in powers of 1/N in a proper thermodynamic limit. For $N\rightarrow +\infty$, we obtain the Vlasov equation describing the evolution of collisionless stellar systems like elliptical galaxies. At the order 1/N, we obtain a kinetic equation describing the evolution of collisional stellar systems like globular clusters. This equation does not suffer logarithmic divergences at large scales since spatial inhomogeneity is explicitly taken into account. Making a local approximation, and introducing an upper cut-off at the Jeans length, it reduces to the Vlasov-Landau equation which is the standard kinetic equation of stellar systems. Our approach provides a simple and pedagogical derivation of these important equations from the BBGKY hierarchy which is more rigorous for systems with long-range interactions than the two-body encounters theory. Making an adiabatic approximation, we write the generalized Landau equation in angle-action variables and obtain a Landau-type kinetic equation that is valid for fully inhomogeneous stellar systems and is free of divergences at large scales. This equation is less general than the Lenard Balescu-type kinetic equation recently derived by Heyvaerts (2010) since it neglects collective effects, but it is substantially simpler and could be useful as a first step. We discuss the evolution of the system as a whole and the relaxation of a test star in a bath of field stars. We derive the corresponding Fokker-Planck equation in angle-action variables and provide expressions for the diffusion coefficient and friction force.

We construct a simple model of universe with a generalized equation of state $p=(\alpha +k\rho^{1/n})\rho c^2$ having a linear component $p=\alpha\rho c^2$ and a polytropic component $p=k\rho^{1+1/n}c^2$. For $\alpha=1/3$, $n=1$ and $k=-4/(3\rho_P)$, where $\rho_P=5.16 10^{99} g/m^3$ is the Planck density, this equation of state provides a model of the early universe without singularity describing the transition between the pre-radiation era and the radiation era. The universe starts from $t=-\infty$ but, when $t<0$, its size is less than the Planck length $l_P=1.62 10^{-35} m$. The universe undergoes an inflationary expansion that brings it to a size $a_1=2.61 10^{-6} m$ on a timescale of a few Planck times $t_P=5.39 10^{-44} s$. When $t\gg t_P$, the universe decelerates and enters in the radiation era. For $\alpha=0$, $n=-1$ and $k=-\rho_{\Lambda}$, where $\rho_{\Lambda}=7.02 10^{-24} g}/m^3$ is the cosmological density, this equation of state describes the transition from a decelerating universe dominated by baryonic and dark matter to an accelerating universe dominated by dark energy (second inflation). The transition takes place at a size $a_2=8.95 10^{25} m$ corresponding to a time of the order of the cosmological time $t_{\Lambda}=1.46 10^{18} s$. This polytropic model reveals a nice "symmetry" between the early and late evolution of the universe, the cosmological constant $\Lambda$ in the late universe playing a role similar to the Planck constant $\hbar$ in the early universe. We interpret the cosmological constant as a fundamental constant of nature describing the "cosmophysics" just like the Planck constant describes the microphysics. The Planck density and the cosmological density represent fundamental upper and lower bounds differing by ${122}$ orders of magnitude. The cosmological constant "problem" may be a false problem.

We construct models of universe with a generalized equation of state $p=(\alpha \rho+k\rho^{1+1/n})c^2$ having a linear component and a polytropic component. The linear equation of state $p=\alpha\rho c^2$ with $-1\le \alpha\le 1$ describes radiation ($\alpha=1/3$), pressureless matter ($\alpha=0$), stiff matter ($\alpha=1$), and vacuum energy ($\alpha=-1$). The polytropic equation of state $p=k\rho^{1+1/n} c^2$ may be due to Bose-Einstein condensates with repulsive ($k>0$) or attractive ($k<0$) self-interaction, or have another origin. In this paper, we consider the case where the density increases as the universe expands. This corresponds to a "phantom universe" for which $w=p/\rho c^2<-1$ (this requires $k<0$). We complete previous investigations on this problem and analyze in detail the different possibilities. We describe the singularities using the classification of [S. Nojiri, S.D. Odintsov, S. Tsujikawa, Phys. Rev. D {\bf 71}, 063004 (2005)]. We show that for $\alpha>-1$ there is no Big Rip singularity although $w\le -1$. For $n=-1$, we provide an analytical model of phantom bouncing universe "disappearing" at $t=0$. We also determine the potential of the phantom scalar field and phantom tachyon field corresponding to the generalized equation of state $p=(\alpha \rho+k\rho^{1+1/n})c^2$.

We construct models of universe with a generalized equation of state $p=(\alpha \rho+k\rho^{1+1/n})c^2$ having a linear component and a polytropic component. In this paper, we consider positive indices $n>0$. In that case, the polytropic component dominates in the early universe where the density is high. For $\alpha=1/3$, $n=1$ and $k=-4/(3\rho_P)$, we obtain a model of early universe describing the transition from a pre-radiation era to the radiation era. The universe exists at any time in the past and there is no singularity. However, for $t<0$, its size is less than the Planck length $l_P=1.62 10^{-35} m$. In this model, the universe undergoes an inflationary expansion with the Planck density $\rho_P=5.16 10^{99} g/m^3$ that brings it to a size $a_1=2.61 10^{-6} m$ at $t_1=1.25 10^{-42} s$ (about 20 Planck times $t_P$). For $\alpha=1/3$, $n=1$ and $k=4/(3\rho_P)$, we obtain a model of early universe with a new form of primordial singularity: The universe starts at t=0 with an infinite density and a finite radius $a=a_1$. Actually, this universe becomes physical at a time $t_i=8.32 10^{-45} s$ from which the velocity of sound is less than the speed of light. When $a\gg a_1$, the universe evolves like in the standard model. We describe the transition from the pre-radiation era to the radiation era by analogy with a second order phase transition where the Planck constant $\hbar$ plays the role of finite size effects (the standard Big Bang theory is recovered for $\hbar=0$).

We construct models of universe with a generalized equation of state $p=(\alpha \rho+k\rho^{1+1/n})c^2$ having a linear component and a polytropic component. In this paper, we consider negative indices $n<0$. In that case, the polytropic component dominates in the late universe where the density is low. For $\alpha=0$, $n=-1$ and $k=-\rho_{\Lambda}$, we obtain a model of late universe describing the transition from the matter era to the dark energy era. The universe exists eternally in the future and undergoes an inflationary expansion with the cosmological density $\rho_{\Lambda}=7.02 10^{-24} g/m^3$ on a timescale $t_{\Lambda}=1.46 10^{18} s$. For $\alpha=0$, $n=-1$ and $k=\rho_{\Lambda}$, we obtain a model of cyclic universe appearing and disappearing periodically. If we were living in this universe, it would disappear in about 2.38 billion years. We make the connection between the early and the late universe and propose a simple equation describing the whole evolution of the universe. This leads to a model of universe that is eternal in past and future without singularity (aioniotic universe). This model exhibits a nice "symmetry" between an early and late phase of inflation, the cosmological constant in the late universe playing the same role as the Planck constant in the early universe. The Planck density and the cosmological density represent fundamental upper and lower bounds differing by 122 orders of magnitude. The cosmological constant "problem" may be a false problem. We determine the potential of the scalar field (quintessence, tachyon field) corresponding to the generalized equation of state $p=(\alpha \rho+k\rho^{1+1/n})c^2$. We also propose a unification of pre-radiation, radiation and dark energy through the quadratic equation of state $p/c^2=-4\rho^2/3\rho_P+\rho/3-4\rho_{\Lambda}/3$.

We develop the suggestion that dark matter could be a Bose-Einstein condensate. We determine the mass-radius relation of a Newtonian self-gravitating Bose-Einstein condensate with short-range interactions described by the Gross-Pitaevskii-Poisson system. We numerically solve the equation of hydrostatic equilibrium describing the balance between the gravitational attraction and the pressure due to quantum effects (Heisenberg's uncertainty principle) and short-range interactions (scattering). We connect the non-interacting limit to the Thomas-Fermi limit. We also consider the case of attractive self-interaction. We compare the exact mass-radius relation obtained numerically with the approximate analytical relation obtained with a Gaussian ansatz. An overall good agreement is found.

We provide an approximate analytical expression of the mass-radius relation of a Newtonian self-gravitating Bose-Einstein condensate (BEC) with short-range interactions described by the Gross-Pitaevskii-Poisson system. These equations model astrophysical objects such as boson stars and, presumably, dark matter galactic halos. Our study connects the non-interacting case studied by Ruffini & Bonazzola (1969) to the Thomas-Fermi limit studied by B\"ohmer & Harko (2007). For repulsive short-range interactions (positive scattering lengths), there exists configurations of arbitrary mass but their radius is always larger than a minimum value. For attractive short-range interactions (negative scattering lengths), equilibrium configurations only exist below a maximum mass. Above that mass, the system is expected to collapse and form a black hole. We also study the radius versus scattering length relation for a given mass. We find that stable configurations only exist above a (negative) minimum scattering length. Our approximate analytical solution, based on a Gaussian ansatz, provides a very good agreement with the exact solution obtained by numerically solving a nonlinear differential equation representing hydrostatic equilibrium. Our treatment is, however, easier to handle and permits to study the stability problem, and derive an analytical expression of the pulsation period, by developing an analogy with a simple mechanical problem.

We analyze the possibility that due to their superfluid properties some compact astrophysical objects may contain a significant part of their matter in the form of a Bose-Einstein condensate. To study the condensate we use the Gross-Pitaevskii equation, with arbitrary non-linearity. By introducing the Madelung representation of the wave function, we formulate the dynamics of the system in terms of the continuity equation and of the hydrodynamic Euler equations. The non-relativistic and Newtonian Bose-Einstein gravitational condensate can be described as a gas, whose density and pressure are related by a barotropic equation of state. In the case of a condensate with quartic non-linearity, the equation of state is polytropic with index one. In the framework of the Thomas-Fermi approximation the structure of the Newtonian gravitational condensate is described by the Lane-Emden equation, which can be exactly solved. The case of the rotating condensate is also discussed. General relativistic configurations with quartic non-linearity are studied numerically with both non-relativistic and relativistic equations of state, and the maximum mass of the stable configuration is determined. Condensates with particle masses of the order of two neutron masses (Cooper pair) and scattering length of the order of 10-20 fm have maximum masses of the order of 2 M_sun, maximum central density of the order of 0.1-0.3 10^16 g/cm^3 and minimum radii in the range of 10-20 km. In this way we obtain a large class of stable astrophysical objects, whose basic astrophysical parameters (mass and radius) sensitively depend on the mass of the condensed particle, and on the scattering length. We also propose that the recently observed neutron stars with masses in the range of 2-2.4 M_sun are Bose-Einstein Condensate stars.

If the dark matter in the universe is a self-gravitating Bose-Einstein condensate (BEC) with quartic self-interaction described by the Gross-Pitaevskii-Poisson system, the adhesion model, the Burgers equation and the cosmological Kardar-Parisi-Zhang (KPZ) equation that have been introduced heuristically to solve the problems inherent to cold dark matter (CDM) models find a natural justification and an interesting generalization.

We study the growth of perturbations in an expanding Newtonian universe with Bose-Einstein condensate dark matter. We first ignore special relativistic effects and derive a differential equation governing the evolution of the density contrast in the linear regime taking into account quantum pressure and self-interaction. This equation can be solved analytically in several cases. We argue that an attractive self-interaction can enhance the Jeans instability and fasten the formation of structures. Then, we take into account pressure effects (coming from special relativity) in the evolution of the cosmic fluid and add the contribution of radiation, baryons and dark energy (cosmological constant). For a BEC dark matter with repulsive self-interaction (positive pressure) the scale factor increases more rapidly than in the standard \Lambda CDM model where dark matter is pressureless while for a BEC dark matter with attractive self-interaction (negative pressure) it increases less rapidly. We study the linear development of the perturbations in these two cases and show that the perturbations grow faster in a BEC dark matter than in a pressureless dark matter. This confirms a recent result of Harko (2011). Finally, we consider a "dark fluid" with a generalized equation of state p=(\alpha \rho + k \rho ^2)c^2 having a component p=k \rho ^2 c^2 similar to a BEC dark matter and a component p=\alpha \rho c^2 mimicking the effect of the cosmological constant (dark energy). We find optimal parameters that give a good agreement with the standard \Lambda CDM model assuming a finite cosmological constant.

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 review the out-of-equilibrium properties of a self-gravitating gas of particles in the presence of a strong friction and a random force (canonical gas). We assume a bare diffusion coefficient of the form $D(\rho)=T\rho^{1/n}$, where $\rho$ is the local particle density, so that the equation of state is $P(\rho)=D(\rho)\rho$. Depending on the spatial dimension $d$, the index $n$, the temperature $T$, and whether the system is confined to a finite box or not, the system can reach an equilibrium state, collapse or evaporate. This article focuses on the latter cases, presenting a complete dynamical phase diagram of the system.

We complete classical investigations concerning the dynamical stability of an infinite homogeneous gaseous medium described by the Euler-Poisson system or an infinite homogeneous stellar system described by the Vlasov-Poisson system (Jeans problem). To determine the stability of an infinite homogeneous stellar system with respect to a perturbation of wavenumber k, we apply the Nyquist method. We first consider the case of single-humped distributions and show that, for infinite homogeneous systems, the onset of instability is the same in a stellar system and in the corresponding barotropic gas, contrary to the case of inhomogeneous systems. We show that this result is true for any symmetric single-humped velocity distribution, not only for the Maxwellian. If we specialize on isothermal and polytropic distributions, analytical expressions for the growth rate, damping rate and pulsation period of the perturbation can be given. Then, we consider the Vlasov stability of symmetric and asymmetric double-humped distributions (two-stream stellar systems) and determine the stability diagrams depending on the degree of asymmetry. We compare these results with the Euler stability of two self-gravitating gaseous streams. Finally, we determine the corresponding stability diagrams in the case of plasmas and compare the results with self-gravitating systems.

We complete our previous investigation concerning the structure and the stability of "isothermal" spheres in general relativity. This concerns objects that are described by a linear equation of state $P=q\epsilon$ so that the pressure is proportional to the energy density. In the Newtonian limit $q\to 0$, this returns the classical isothermal equation of state. We consider specifically a self-gravitating radiation (q=1/3), the core of neutron stars (q=1/3) and a gas of baryons interacting through a vector meson field (q=1). We study how the thermodynamical parameters scale with the size of the object and find unusual behaviours due to the non-extensivity of the system. We compare these scaling laws with the area scaling of the black hole entropy. We also determine the domain of validity of these scaling laws by calculating the critical radius above which relativistic stars described by a linear equation of state become dynamically unstable. For photon stars, we show that the criteria of dynamical and thermodynamical stability coincide. Considering finite spheres, we find that the mass and entropy as a function of the central density present damped oscillations. We give the critical value of the central density, corresponding to the first mass peak, above which the series of equilibria becomes unstable. Finally, we extend our results to d-dimensional spheres. We show that the oscillations of mass versus central density disappear above a critical dimension d_{crit}(q). For Newtonian isothermal stars (q=0) we recover the critical dimension d_{crit}=10. For the stiffest stars (q=1) we find d_{crit}=9 and for a self-gravitating radiation (q=1/d) we find d_{crit}=9.96404372... very close to 10. Finally, we give analytical solutions of relativistic isothermal spheres in 2D gravity.

We discuss the distribution of the gravitational force created by a Poissonian distribution of field sources (stars, galaxies,...) in different dimensions of space d. In d=3, it is given by a Levy law called the Holtsmark distribution. It presents an algebraic tail for large fluctuations due to the contribution of the nearest neighbor. In d=2, it is given by a marginal Gaussian distribution intermediate between Gaussian and Levy laws. In d=1, it is exactly given by the Bernouilli distribution (for any particle number N) which becomes Gaussian for N>>1. Therefore, the dimension d=2 is critical regarding the statistics of the gravitational force. We generalize these results for inhomogeneous systems with arbitrary power-law density profile and arbitrary power-law force in a d-dimensional universe.

We determine the mass-radius relation of relativistic white dwarf stars (a self-gravitating degenerate Fermi gas at T=0) in a D-dimensional universe and study the influence of the dimension of space on the laws of physics when we combine quantum mechanics, special relativity and gravity. We exhibit characteristic dimensions D=1, D=2, D=3, D={1\over 2}(3+\sqrt{17}), D=4, D=2(1+\sqrt{2}) and show that quantum mechanics cannot balance gravitational collapse for D\ge 4, even in the non-relativistic regime. This is similar to a result found by Ehrenfest (1917) at the atomic level for Coulomb forces (in Bohr's model) and for the Kepler problem. This makes the dimension of our universe D=3 very particular. We discuss some historic aspects concerning the discovery of the Chandrasekhar (1931b) limiting mass in relation with previous investigations by Anderson (1929) and Stoner (1930). We also show the analogy between the Chandrasekhar limiting mass of white dwarf stars and the critical mass of biological populations in the two-dimensional chemotactic problem.

We study the influence of the dimension of space on the thermodynamics of the classical and quantum self-gravitating gas. We consider Hamiltonian systems of self-gravitating particles described by the microcanonical ensemble and self-gravitating Brownian particles described by the canonical ensemble. We present a gallery of caloric curves in different dimensions of space and discuss the nature of phase transitions as a function of the dimension d. We also provide the general form of the Virial theorem in d dimensions and discuss the particularity of the dimension d=4 for Hamiltonian systems and the dimension d=2 for Brownian systems.

We complete previous investigations on the dynamical stability of barotropic stars and collisionless stellar systems. A barotropic star that minimizes the energy functional at fixed mass is a nonlinearly dynamically stable stationary solution of the Euler-Poisson system. Formally, this minimization problem is similar to a condition of ``canonical stability'' in thermodynamics. A stellar system that maximizes an H-function at fixed mass and energy is a nonlinearly dynamically stable stationary solution of the Vlasov-Poisson system. Formally, this maximization problem is similar to a condition of ``microcanonical stability'' in thermodynamics. Using a thermodynamical analogy, we provide a derivation and an interpretation of the nonlinear Antonov first law in terms of ``ensembles inequivalence'': a spherical stellar system with f=f(epsilon) and f'(epsilon)<0 is nonlinearly dynamically stable with respect to the Vlasov-Poisson system if the corresponding barotropic star with the same equilibrium density distribution is nonlinearly dynamically stable with respect to the Euler-Poisson system. This is similar to the fact that ``canonical stability implies microcanonical stability'' in thermodynamics. The converse is wrong in case of ``ensembles inequivalence'' which is generic for systems with long-range interactions like gravity. We show that criteria of nonlinear dynamical stability can be obtained very simply from purely graphical constructions by using the method of series of equilibria and the turning point argument of Poincare, as in thermodynamics.

We develop the analogy between self-gravitating Brownian particles and bacterial populations. In the high friction limit, the self-gravitating Brownian gas is described by the Smoluchowski-Poisson system. These equations can develop a self-similar collapse leading to a finite time singularity. Coincidentally, the Smoluchowski-Poisson system corresponds to a simplified version of the Keller-Segel model of bacterial populations. In this biological context, it describes the chemotactic aggregation of the bacterial colonies. We extend these classical models by introducing a small-scale regularization. In the gravitational context, we consider a gas of self-gravitating Brownian fermions and in the biological context we consider finite size effects. In that case, the collapse stops when the system feels the influence of the small-scale regularization. A phenomenon of ``explosion'', reverse to the collapse, is also possible.

We investigate a model describing the dynamics of a gas of self-gravitating Brownian particles. This model can also have applications for the chemotaxis of bacterial populations. We focus here on the collapse phase obtained at sufficiently low temperature/energy and on the post-collapse regime following the singular time where the central density diverges. Several analytical results are illustrated by numerical simulations.

We discuss the physical basis of the statistical mechanics of self-gravitating systems. We show the correspondance between statistical mechanics methods based on the evaluation of the density of states and partition function and thermodynamical methods based on the maximization of a thermodynamical potential (entropy or free energy). We address the question of the thermodynamic limit of self-gravitating systems, the justification of the mean-field approximation, the validity of the saddle point approximation near the transition point, the lifetime of metastable states and the fluctuations in isothermal spheres. In particular, we emphasize the tremendously long lifetime of metastable states of self-gravitating systems which increases exponentially with the number of particles N except in the vicinity of the critical point. More specifically, using an adaptation of the Kramers formula justified by a kinetic theory, we show that the lifetime of a metastable state scales as e^{N Delta s} in microcanonical ensemble and e^{N Delta j} in canonical ensemble, where Delta s and Delta j are the barriers of entropy and free energy j=s-beta epsilon (per particle) respectively. The physical caloric curve must take these metastable states (local entropy maxima) into account. As a result, it becomes multi-valued and leads to microcanonical phase transitions and ``dinosaur's necks''(Chavanis 2002b, Chavanis & Rieutord 2003). The consideration of metastable states answers the critics raised by D.H.E. Gross [cond-mat/0307535; cond-mat/0403582].

We address the post-collapse dynamics of a self-gravitating gas of Brownian particles in D dimensions, in both canonical and microcanonical ensembles. In the canonical ensemble, the post-collapse evolution is marked by the formation of a Dirac peak with increasing mass. The density profile outside the peak evolves self-similarly with decreasing central density and increasing core radius. In the microcanonical ensemble, the post-collapse regime is marked by the formation of a ``binary''-like structure surrounded by an almost uniform halo with high temperature. These results are consistent with thermodynamical predictions.

We determine an asymptotic expression of the blow-up time t_coll for self-gravitating Brownian particles or bacterial populations (chemotaxis) close to the critical point. We show that t_coll=t_{*}(eta-eta_c)^{-1/2} with t_{*}=0.91767702..., where eta represents the inverse temperature (for Brownian particles) or the mass (for bacterial colonies), and eta_c is the critical value of eta above which the system blows up. This result is in perfect agreement with the numerical solution of the Smoluchowski-Poisson system. We also determine the asymptotic expression of the relaxation time close but above the critical temperature and derive a large time asymptotic expansion for the density profile exactly at the critical point.

We compute statistical equilibrium states of rotating self-gravitating systems enclosed within a box by maximizing the Fermi-Dirac entropy at fixed mass, energy and angular momentum. We increase the rotation up to the Keplerian limit and describe the flattening of the configuration until mass shedding occurs. At the maximum rotation, the system develops a cusp at the equator. We draw the equilibrium phase diagram of the rotating self-gravitating Fermi gas and discuss the structure of the caloric curve as a function of degeneracy parameter and angular velocity. We argue that systems described by the Fermi-Dirac distribution in phase space do not bifurcate to non-axisymmetric structures, in continuity with the case of polytropes with index n>0.808 (the Fermi gas at T=0 corresponds to n=3/2). This contrasts with the study of Votyakov et al. (2002) who consider a Fermi-Dirac distribution in configuration space and find ``double star'' structures (their model at T=0 corresponds to n=0). We also discuss the influence of rotation on the onset of the gravothermal catastrophe for globular clusters. On general grounds, we complete previous investigations concerning the nature of phase transitions in self-gravitating systems. We emphasize the inequivalence of statistical ensembles regarding the formation of binaries (or low-mass condensates) in the microcanonical ensemble and Dirac peaks (or massive condensates) in the canonical ensemble. We also describe an hysteretic cycle between the gaseous phase and the condensed phase that are connected by a ``collapse'' or an ``explosion''. This notion of hysteresis in self-gravitating systems is new.

We address the generalized thermodynamics and the collapse of a system of self-gravitating Langevin particles exhibiting anomalous diffusion in a space of dimension D. The equilibrium states correspond to polytropic distributions. The index n of the polytrope is related to the exponent of anomalous diffusion. We consider a high-friction limit and reduce the problem to the study of the nonlinear Smoluchowski-Poisson system. We show that the associated Lyapunov functional is the Tsallis free energy. We discuss in detail the equilibrium phase diagram of self-gravitating polytropes as a function of D and n and determine their stability by using turning points arguments and analytical methods. When no equilibrium state exists, we investigate self-similar solutions describing the collapse. These results can be relevant for astrophysical systems, two-dimensional vortices and for the chemotaxis of bacterial populations. Above all, this model constitutes a prototypical dynamical model of systems with long-range interactions which possesses a rich structure and which can be studied in great detail.

We complete previous investigations on the thermodynamics of self-gravitating systems by studying the grand canonical, grand microcanonical and isobaric ensembles. We also discuss the stability of polytropic spheres in the light of a generalized thermodynamics proposed by Tsallis. We determine in each case the onset of gravitational instability by analytical methods and graphical constructions in the Milne plane. We also discuss the relation between dynamical and thermodynamical stability of stellar systems and gaseous spheres. Our study provides an aesthetic and simple approach to this otherwise complicated subject.

The formation of large-scale vortices is an intriguing phenomenon in two-dimensional turbulence. Such organization is observed in large-scale oceanic or atmospheric flows, and can be reproduced in laboratory experiments and numerical simulations. A general explanation of this organization was first proposed by Onsager (1949) by considering the statistical mechanics for a set of point vortices in two-dimensional hydrodynamics. Similarly, the structure and the organization of stellar systems (globular clusters, elliptical galaxies,...) in astrophysics can be understood by developing a statistical mechanics for a system of particles in gravitational interaction as initiated by Chandrasekhar (1942). These statistical mechanics turn out to be relatively similar and present the same difficulties due to the unshielded long-range nature of the interaction. This analogy concerns not only the equilibrium states, i.e. the formation of large-scale structures, but also the relaxation towards equilibrium and the statistics of fluctuations. We will discuss these analogies in detail and also point out the specificities of each system.

We discuss the statistical mechanics of violent relaxation in stellar systems following the pioneering work of Lynden-Bell (1967). The solutions of the gravitational Vlasov-Poisson system develop finer and finer filaments so that a statistical description is appropriate to smooth out the small-scales and describe the ``coarse-grained'' dynamics. In a coarse-grained sense, the system is expected to reach an equilibrium state of a Fermi-Dirac type within a few dynamical times. We describe in detail the equilibrium phase diagram and the nature of phase transitions which occur in self-gravitating systems. Then, we introduce a small-scale parametrization of the Vlasov equation and propose a set of relaxation equations for the coarse-grained dynamics. These relaxation equations, of a generalized Fokker-Planck type, are derived from a Maximum Entropy Production Principle (MEPP). We make a link with the quasilinear theory of the Vlasov-Poisson system and derive a truncated model appropriate to collisionless systems subject to tidal forces. With the aid of this kinetic theory, we qualitatively discuss the concept of ``incomplete relaxation'' and the limitations of Lynden-Bell's theory.

Self-gravitating systems are expected to reach a statistical equilibrium state either through collisional relaxation or violent collisionless relaxation. However, a maximum entropy state does not always exist and the system may undergo a ``gravothermal catastrophe'': it can achieve ever increasing values of entropy by developing a dense and hot ``core'' surrounded by a low density ``halo''. In this paper, we study the phase transition between ``equilibrium'' states and ``collapsed'' states with the aid of a simple relaxation equation [Chavanis, Sommeria and Robert, Astrophys. J. 471, 385 (1996)] constructed so as to increase entropy with an optimal rate while conserving mass and energy. With this numerical algorithm, we can cover the whole bifurcation diagram in parameter space and check, by an independent method, the stability limits of Katz [Mon. Not. R. astr. Soc. 183, 765 (1978)] and Padmanabhan [Astrophys. J. Supp. 71, 651 (1989)]. When no equilibrium state exists, our relaxation equation develops a self-similar collapse leading to a finite time singularity.

We discuss the nature of phase transitions in the self-gravitating Fermi gas at non-zero temperature. This study can be relevant for massive neutrinos in Dark Matter models and for collisionless self-gravitating systems experiencing a ``violent relaxation''. Below a critical energy (in the microcanonical ensemble) or below a critical temperature (in the canonical ensemble), the system undergoes a gravitational collapse leading to a compact object with a massive degenerate nucleus (fermion ball).

Phase diagram of microcanonical ensembles of self-attracting particles is studied for two types of short-range potential regularizations: self-gravitating fermions and classical particles interacting via attractive soft $-(r^2+r_0^2)^{-1/2}$ Coulomb potential. When the range of regularization is sufficiently short, the self-attracting systems exhibit gravitational or collapse-like transition. As the fermionic degeneracy or the softness radius increases, the gravitational phase transition crosses over to a normal first-order phase transition, becomes second-order at a critical point, and finally disappears. Applicability of a commonly used saddle-point or mean-field approximation and importance of metastable states is discussed.

We address the thermodynamics (equilibrium density profiles, phase diagram, instability analysis...) and the collapse of a self-gravitating gas of Brownian particles in D dimensions, in both canonical and microcanonical ensembles. In the canonical ensemble, we derive the analytic form of the density scaling profile which decays as f(x)=x^{-\alpha}, with alpha=2. In the microcanonical ensemble, we show that f decays as f(x)=x^{-\alpha_{max}}, where \alpha_{max} is a non-trivial exponent. We derive exact expansions for alpha_{max} and f in the limit of large D. Finally, we solve the problem in D=2, which displays rather rich and peculiar features.

We discuss the statistical mechanics of rotating self-gravitating systems by allowing properly for the conservation of angular momentum. We study analytically the case of slowly rotating isothermal spheres by expanding the solutions of the Boltzmann-Poisson equation in a series of Legendre polynomials, adapting the procedure introduced by Chandrasekhar (1933) for distorted polytropes. We show how the classical spiral of Lynden-Bell & Wood (1967) in the temperature-energy plane is deformed by rotation. We find that gravitational instability occurs sooner in the microcanonical ensemble and later in the canonical ensemble. According to standard turning point arguments, the onset of the collapse coincides with the minimum energy or minimum temperature state in the series of equilibria. Interestingly, it happens to be close to the point of maximum flattening. We determine analytically the generalization of the singular isothermal solution to the case of a slowly rotating configuration. We also consider slowly rotating configurations of the self-gravitating Fermi gas at non zero temperature.

We discuss the nature of phase transitions in self-gravitating systems both in the microcanonical and in the canonical ensemble. We avoid the divergence of the gravitational potential at short distances by considering the case of self-gravitating fermions and hard spheres models. Three kinds of phase transitions (of zeroth, first and second order) are evidenced. They separate a ``gaseous'' phase with a smoothly varying distribution of matter from a ``condensed'' phase with a core-halo structure. We propose a simple analytical model to describe these phase transitions. We determine the value of energy (in the microcanonical ensemble) and temperature (in the canonical ensemble) at the transition point and we study their dependance with the degeneracy parameter (for fermions) or with the size of the particles (for a hard spheres gas). Scaling laws are obtained analytically in the asymptotic limit of a small short distance cut-off. Our analytical model captures the essential physics of the problem and compares remarkably well with the full numerical solutions.

We complete the existing literature on the structure and stability of polytropic gas spheres reported in the classical monograph of Chandrasekhar (1942). For isolated polytropes with index $1<n<5$, we provide a new, alternative, proof that the onset of instability occurs for $n=3$ and we express the perturbation profiles of density and velocity at the point of marginal stability in terms of the Milne variables. Then, we consider the case of polytropes confined within a box of radius $R$ (an extension of the Antonov problem for isothermal gas spheres). For $n\ge 3$, the mass-density relation presents some damped oscillations and there exists a limiting mass above which no hydrostatic equilibrium is possible. Like for isothermal gas spheres, the onset of instability occurs precisely at the point of maximum mass. Analytical results are obtained for the particular index $n=5$. We also discuss the relation of our study with extended thermodynamics (Tsallis entropy) recently investigated by Taruya & Sakagami (cond-mat/0107494).

We investigate the effects of relativity on the gravitational instability of finite isothermal gaseous spheres. In the first part of the paper, we treat the gravitational field within the framework of Newtonian mechanics but we use a relativistic equation of state in the condition of hydrostatic equilibrium. In the second part of the paper, we study the full general relativistic problem for a gas described by an equation of state $p=q\epsilon$ such that the pressure is proportional to the energy density (``isothermal'' distribution). For $q=1/3$, this equation of state describes the core of neutron stars. The mass-density diagram displays some damped oscillations and there exists a critical value of mass-energy above which no equilibrium state is possible. We show analytically that the mass peaks are associated with new modes of instability. These results are strikingly similar to those obtained by Antonov [Vest. Leningr. Gos. Univ. 7, 135 (1962)] and Lynden-Bell & Wood (1968) for a classical isothermal gas. Our study completes the analogy between isothermal spheres and neutron stars investigated by Yabushita [MNRAS 167, 95 (1974)].

We investigate the stability of bounded self-gravitating systems in the canonical ensemble by using a thermodynamical approach. Our study extends the earlier work of Padmanabhan [Astrophys. J. Supp. 71, 651 (1989)] in the microcanonical ensemble. By studying the second variations of the free energy, we find that instability sets in precisely at the point of minimum temperature in agreement with the theorem of Katz [ Mon. Not. R. astr. Soc. 183, 765 (1978)]. The perturbation that induces instability at this point is calculated explicitly; it has not a ``core-halo'' structure contrary to what happens in the microcanonical ensemble. We also study Jeans type gravitational instability of isothermal gaseous spheres described by Navier-Stokes equations. We show analytically the equivalence between dynamical stability and thermodynamical stability and the fact that the stability of isothermal gas spheres does not depend on the viscosity. This confirms the findings of Semelin et al. [astro-ph/9908073] who used numerical methods. We also give a simpler derivation of the geometric hierarchy of scales inducing instability discovered by these authors. The density profiles that trigger these instabilities are calculated analytically; they present more and more oscillations that also follow a geometric progression. This suggests that the system will fragmentate in a series of `clumps' and that these `clumps' will themselves fragmentate in substructures. The fact that both the domain sizes leading to instability and the `clumps' sizes within a box follow a geometric progression with the same ratio suggests a fractal-like behaviour. This gives further support to the interpretation of de Vega et al. [Nature, 383, 56 (1996)].

We develop the idea proposed by Barge & Sommeria (1995) and Tanga et al. (1996) that large-scale vortices present in the solar nebula can concentrate dust particles and facilitate the formation of planetesimals and planets. We introduce an exact vortex solution of the incompressible 2D Euler equation and study the motion of dust particles in that vortex. In particular, we derive analytical expressions for the capture time and the mass capture rate as a function of the friction parameter. Then, we study how small-scale turbulent fluctuations affect the motion of the particles in the vortex and determine their rate of escape by solving a problem of quantum mechanics. We apply these results to the solar nebula and find that the capture is optimum near Jupiter's orbit (as noticed already by Barge & Sommeria 1995) but also in the Earth region. This second optimum corresponds to the transition between the Epstein and the Stokes regime which takes place, for relevant particles, at the separation between telluric and giant planets (i.e near the asteroid belt). At these locations, the particles are efficiently captured and concentrated by the vortices and can undergo gravitational collapse to form the planetesimals.

In the problem of planetary formation one seeks a mechanism to gather small solid particles together into larger accumulations of solid matter. Here we describe a scenario in which turbulence mediates this process by aggregating particles into anticyclonic regions. If, as our simulations suggest, anticyclonic vortices form as long-lived coherent structures, the process becomes more powerful because such vortices trap particles effectively. Even if the turbulence is decaying, following the upheaval that formed the disk, there is enough time to make the dust distribution quite lumpy.