BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Pierre-Henri Chavanis

Pierre-Henri Chavanis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechastro-ph.COgr-qcastro-ph.GAphysics.flu-dynastro-ph.HEastro-ph.SRphysics.ao-phQuantitative Methodsastro-phBiological Physicshep-phhep-thphysics.plasm-ph

Catalog footprint

What is connected

77works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

77 published item(s)

preprint2022arXiv

Generalized Logotropic Models and their Cosmological Constraints

We propose a new class of cosmological unified dark sector models called "{\em Generalized Logotropic Models}". They depend on a free parameter $n$. The original logotropic model [P.H. Chavanis, Eur. Phys. J. Plus {\bf 130}, 130 (2015)] is a special case of our generalized model corresponding to $n=1$. In our scenario, the Universe is filled with a single fluid, a generalized logotropic dark fluid (GLDF), whose pressure $P$ includes higher order logarithmic terms of the rest-mass density $ρ_m$. The total energy density $ε$ is the sum of the rest-mass energy density $ρ_m c^2$ and the internal energy density $u$ which play the role of dark matter energy density $ε_m$ and dark energy density $ε_{de}$, respectively. We investigate the cosmological behavior of the generalized logotropic models by focusing on the evolution of the energy density, scale factor, equation of state parameter, decceleration parameter and squared speed of sound. Low values of $n\le 3$ are favored. We also study the asymptotic behavior of the generalized logotropic models. In particular, we show that the model presents a phantom behavior and has three distinct ways of evolution depending on the value of $n$. For $n\le 2$, it leads to a little rip and for $n>2$ to a big rip. We predict the value of the big rip time as a function of $n$ without any free (undetermined) parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Orbital alignment and mass segregation in galactic nuclei via vector resonant relaxation

Supermassive black holes dominate the gravitational potential in galactic nuclei. In these dense environments, stars follow nearly Keplerian orbits and see their orbital planes relax through the potential fluctuations generated by the stellar cluster itself. For typical astrophysical galactic nuclei, the most likely outcome of this vector resonant relaxation (VRR) is that the orbital planes of the most massive stars spontaneously self-align within a narrow disc. We present a maximum entropy method to systematically determine this long-term distribution of orientations and use it for a wide range of stellar orbital parameters and initial conditions. The heaviest stellar objects are found to live within a thin equatorial disk. The thickness of this disk depends on the stars' initial mass function, and on the geometry of the initial cluster. This work highlights a possible (indirect) novel method to constrain the distribution of intermediate mass black holes in galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predictions from the logotropic model: the universal surface density of dark matter halos and the present proportion of dark matter and dark energy

The logotropic model [P.H. Chavanis, Eur. Phys. J. Plus {\bf 130}, 130 (2015)] may be an interesting alternative to the $Λ$CDM model. It is able to account for the present accelerating expansion of the universe while solving at the same time the core-cusp problem of the CDM model. In the logotropic model, there is a single dark fluid. Its rest-mass plays the role of dark matter and its internal energy plays the role of dark energy. We highlight two remarkable predictions of the logotropic model. It yields cored dark matter halos with a universal surface density equal to $Σ_0^{\rm th}=0.01955 c\sqrtΛ/G=133\, M_{\odot}/{\rm pc}^2$ without free parameter in very good agreement with the observational value $Σ_0^{\rm obs}=141_{-52}^{+83}\, M_{\odot}/{\rm pc}^2$. It also predicts the present ratio of dark energy and dark matter to be the pure number $Ω_{\rm de,0}^{\rm th}/Ω_{\rm dm,0}^{\rm th}=e=2.71828...$ in very good agreement with the observations giving $Ω_{\rm de,0}^{\rm obs}/Ω_{\rm dm,0}^{\rm obs}=2.669\pm 0.08$. Using the measured present proportion of baryonic matter $Ω_{\rm b,0}^{\rm obs}=0.0486\pm 0.0010$, we find that the values of the present proportion of dark matter and dark energy are $Ω_{\rm dm,0}^{\rm th}=\frac{1}{1+e}(1-Ω_{\rm b,0})=0.2559$ and $Ω_{\rm de,0}^{\rm th}=\frac{e}{1+e}(1-Ω_{\rm b,0})=0.6955$ in very good agreement with the observational values $Ω_{\rm dm,0}^{\rm obs}=0.2589\pm 0.0057$ and $Ω_{\rm de,0}^{\rm obs}=0.6911\pm 0.0062$ within the error bars. These theoretical predictions are obtained by advocating a mysterious strong cosmic coincidence (dubbed "dark magic") implying that our epoch plays a particular role in the history of the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Self-gravitating clusters of Bose-Einstein gas with planar, cylindrical, or spherical symmetry: gaseous density profiles and onset of condensation

We calculate density profiles for self-gravitating clusters of an ideal Bose-Einstein gas with nonrelativistic energy-momentum relation and macroscopic mass at thermal equilibrium. Our study includes clusters with planar symmetry in dimensions $\mathcal{D}=1,2,3$, clusters with cylindrical symmetry in $\mathcal{D}=2,3$, and clusters with spherical symmetry in $\mathcal{D}=3$. Wall confinement is imposed where needed to prevent escape. The length scale and energy scale in use for the gaseous phase render density profiles for gaseous macrostates independent of total mass. Density profiles for mixed-phase macrostates have a condensed core surrounded by a gaseous halo. The spatial extension of the core is negligibly small on the length scale tailored for the halo. The mechanical stability conditions as evident in caloric curves permit multiple macrostates to coexist. Their status regarding thermal equilibrium is examined by a comparison of free energies. The onset of condensation takes place at a nonzero temperature in all cases. The critical singularities and the nature of the phase transition vary with the symmetry of the cluster and the dimensionality of the space.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A heuristic wave equation parameterizing BEC dark matter halos with a quantum core and an isothermal atmosphere

The Gross-Pitaevskii-Poisson equations that govern the evolution of self-gravitating Bose-Einstein condensates, possibly representing dark matter halos, experience a process of gravitational cooling and violent relaxation. We propose a heuristic parametrization of this complicated process in the spirit of Lynden-Bell's theory of violent relaxation for collisionless stellar systems. We derive a generalized wave equation (that was introduced phenomenologically in [P.H. Chavanis, Eur. Phys. J. Plus {\bf 132}, 248 (2017)]) involving a logarithmic nonlinearity associated with an effective temperature $T_{\rm eff}$ and a damping term associated with a friction $ξ$. These terms can be obtained from a maximum entropy production principle and are linked by a form of Einstein relation expressing the fluctuation-dissipation theorem. The wave equation satisfies an $H$-theorem for the Lynden-Bell entropy and relaxes towards a stable equilibrium state which is a maximum of entropy at fixed mass and energy. This equilibrium state represents the most probable state of a Bose-Einstein condensate dark matter halo. It generically has a core-halo structure. The quantum core prevents gravitational collapse and may solve the core-cusp problem. The isothermal halo leads to flat rotation curves in agreement with the observations. These results are consistent with the phenomenology of dark matter halos. Furthermore, as shown in a previous paper [P.H. Chavanis, Phys. Rev. D {\bf 100}, 123506 (2019)], the maximization of entropy with respect to the core mass at fixed total mass and total energy determines a core mass--halo mass relation which agrees with the relation obtained in direct numerical simulations. We stress the importance of using a microcanonical description instead of a canonical one. We also explain how our formalism can be applied to the case of fermionic dark matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Cosmological models based on a complex scalar field with a power-law potential associated with a polytropic equation of state

We construct cosmological models based on a complex scalar field with a power-law potential $V=\frac{K}{γ-1}(\frac{m}{\hbar})^{2γ}|φ|^{2γ}$ associated with a polytropic equation of state $P=Kρ^γ$ (the potential associated with an isothermal equation of state $P=ρk_B T/m$ is $V=\frac{m k_B T}{\hbar^2}|φ|^2 [\ln(m^2|φ|^2/ρ_*\hbar^2)-1]$ and the potential associated with a logotropic equation of state $P=A\ln(ρ/ρ_P)$ is $V=-A[\ln(m^2|φ|^2/\hbar^2ρ_P)+1]$). We consider a fast oscillation regime of ``spintessence'' where the equations of the problem can be simplified. We study all possible cases with arbitrary (positive and negative) values of the polytropic constant and polytropic index. The $Λ$CDM model, the Chaplygin gas model and the Bose-Einstein condensate model are recovered as particular cases of our study corresponding to a constant potential ($γ=0$), an inverse square-law potential ($γ=-1$), and a quartic potential ($γ=2$). We also derive the two-fluid representation of the Chaplygin gas model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Predictive model of fermionic dark matter halos with a quantum core and an isothermal atmosphere

We develop a thermodynamical model of fermionic dark matter halos at finite temperature. Statistical equilibrium states may be justified by a process of violent collisionless relaxation in the sense of Lynden-Bell or from a collisional relaxation of nongravitational origin if the fermions are self-interacting. The most probable state (maximum entropy state) generically has a "core-halo" structure with a quantum core (fermion ball) surrounded by an isothermal atmosphere. The quantum core is equivalent to a polytrope of index $n=3/2$. The Pauli exclusion principle creates a quantum pressure that prevents gravitational collapse and solves the core-cusp problem of the cold dark matter model. The isothermal atmosphere (which is similar to the NFW profile of cold dark matter) accounts for the flat rotation curves of the galaxies at large distances. We numerically solve the equation of hydrostatic equilibrium with the Fermi-Dirac equation of state and determine the density profiles and rotation curves of fermionic dark matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Self-gravitating clusters of Fermi-Dirac gas with planar, cylindrical, or spherical symmetry: evolution of density profiles with temperature

We calculate density profiles for self-gravitating clusters of an ideal Fermi-Dirac gas with nonrelativistic energy-momentum relation and macroscopic mass at thermal equilibrium. Our study includes clusters with planar symmetry in dimensions $\mathcal{D}=1,2,3$, clusters with cylindrical symmetry in $\mathcal{D}=2,3$, and clusters with spherical symmetry in $\mathcal{D}=3$. Wall confinement is imposed where needed for stability against escape. The length scale and energy scale in use render all results independent of total mass and prove adequate at all temperatures. We present exact analytic expressions for (fully degenerate) $T=0$ density profiles in four of the six combinations of symmetry and dimensionality. Our numerical results for $T>0$ describe the emergence, upon quasistatic cooling, of a core with incipient degeneracy surrounded by a more dilute halo. The equilibrium macrostates are found to depend more strongly on the cluster symmetry than on the space dimensionality. We demonstrate the mechanical and thermal stability of spherical clusters with coexisting phases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Landau equation for self-gravitating classical and quantum particles: Application to dark matter

We develop the kinetic theory of classical and quantum particles (fermions and bosons) in gravitational interaction. The kinetic theory of quantum particles may have applications in the context of dark matter. For simplicity, we consider an infinite and spatially homogeneous system (or make a local approximation) and neglect collective effects. This leads to the quantum Landau equation derived heuristically in [Chavanis, Physica A 332, 89 (2004)]. We establish its main properties: conservation laws, $H$-theorem, equilibrium state, relaxation time, quantum diffusion and friction coefficients, quantum Rosenbluth potentials, self-consistent evolution, (thermal) bath approximation, quantum Fokker-Planck equation, quantum King model... For bosonic particles, the Landau equation can describe the process of Bose-Einstein condensation. We discuss the relation of our study with the works of [Levkov et al., Phys. Rev. Lett. 121, 151301 (2018); Bar-Or et al., Astrophys. J. 871, 28 (2019)] on fuzzy dark matter halos and the formation of Bose stars and solitons.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Caloric curves of classical self-gravitating systems in general relativity

We determine the caloric curves of classical self-gravitating systems at statistical equilibrium in general relativity. In the classical limit, the caloric curves of a self-gravitating gas depend on a unique parameter $ν=GNm/Rc^2$, called the compactness parameter, where $N$ is the particle number and $R$ the system's size. Typically, the caloric curves have the form of a double spiral. The "cold spiral", corresponding to weakly relativistic configurations, is a generalization of the caloric curve of nonrelativistic classical self-gravitating systems. The "hot spiral'", corresponding to strongly relativistic configurations, is similar (but not identical) to the caloric curve of the ultrarelativistic self-gravitating black-body radiation. We introduce two types of normalization of energy and temperature in order to obtain asymptotic caloric curves describing respectively the cold and the hot spirals in the limit $ν\rightarrow 0$. As the number of particles increases, the cold and the hot spirals approach each other, merge at $ν'_S=0.128$, form a loop above $ν_S=0.1415$, reduce to a point at $ν_{\rm max}=0.1764$, and finally disappear. Therefore, the double spiral shrinks when the compactness parameter $ν$ increases, implying that general relativistic effects render the system more unstable. We discuss the nature of the gravitational collapse at low and high energies with respect to a dynamical (fast) or a thermodynamical (slow) instability.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Core mass -- halo mass relation of bosonic and fermionic dark matter halos harbouring a supermassive black hole

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Gravitational phase transition of self-gravitating systems of fermions in General Relativity

The Thomas-Fermi model is extended at finite temperature, to describe the gravitational phase transition occurring in massive fermionic systems in a general-relativistic framework. It is shown that, when a nondegenerate fermionic gas (for $N < N_{OV}$, where $N_{OV}$ is the Oppenheimer-Volkoff limit) is cooled down below a critical temperature, a condensed phase emerges and the gravitational collapse is prevented by quantum mechanics. If $N > N_{OV}$, by contrast, the system is destined to collapse towards a Black Hole because no equilibrium states exist.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Gravitational phase transitions and instabilities of self-gravitating fermions in general relativity

We discuss the occurrence of gravitational phase transitions and instabilities in a gas of self-gravitating fermions within the framework of general relativity. In the classical (nondegenerate) limit, the system undergoes a gravitational collapse at low energies $E<E_c$ and low temperatures $T<T_c$. This is called "gravothermal catastrophe" in the microcanonical ensemble and "isothermal collapse" in the canonical ensemble. When quantum mechanics is taken into account and when the particle number is below the Oppenheimer-Volkoff limit ($N<N_{\rm OV}$), complete gravitational collapse is prevented by the Pauli exclusion principle. In that case, the Fermi gas undergoes a gravitational phase transition from a gaseous phase to a condensed phase. The condensed phase represents a compact object like a white dwarf, a neutron star, or a dark matter fermion ball. When $N>N_{\rm OV}$, there can be a subsequent gravitational collapse below a lower critical energy $E<E''_c$ or a lower critical temperature $T<T'_c$ leading presumably to the formation of a black hole. The evolution of the system is different in the microcanonical and canonical ensembles. In the microcanonical ensemble, the system takes a "core-halo" structure. The core consists in a compact quantum object or a black hole while the hot halo is expelled at large distances. This is reminiscent of the red giant structure of low-mass stars or the implosion-explosion of massive stars (supernova). In the canonical ensemble, the system collapses as a whole towards a compact object or a black hole. This is reminiscent of the implosion of supermassive stars (hypernova).

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Implosion-explosion in supernovae

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é 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Relativistic Gravitational Phase Transitions and Instabilities of the Fermi Gas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Supernova implosion-explosion in the light of catastrophe theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy Balance of a Bose Gas in Curved Spacetime

We derive a general energy balance equation for a self-interacting boson gas at vanishing temperature in a curved spacetime. This represents a first step towards a formulation of the first law of thermodynamics for a scalar field in general relativity. By using a $3+1$ foliation of the spacetime and performing a Madelung transformation, we rewrite the Klein-Gordon-Maxwell equations in a general curved spacetime into its hydrodynamic version where we can identify the different energy contributions of the system and separate them into kinetic, quantum, electromagnetic, and gravitational.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Euler-like modelling of dense granular flows: application to a rotating drum

General conservation equations are derived for 2D dense granular flows from the Euler equation within the Boussinesq approximation. In steady flows, the 2D fields of granular temperature, vorticity and stream function are shown to be encoded in two scalar functions only. We checked such prediction on steady surface flows in a rotating drum simulated through the Non-Smooth Contact Dynamics method. This result is non trivial because granular flows are dissipative and therefore not necessarily compatible with Euler equation. Finally, we briefly discuss some possible ways to predict theoretically these two functions using statistical mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A cosmological model describing the early inflation, the intermediate decelerating expansion, and the late accelerating expansion by a quadratic equation of state

We develop a cosmological model based on a quadratic equation of state p/c^2=-(α+1){ρ^2}/{ρ_P}+αρ-(α+1)ρ_Λ (where ρ_P is the Planck density and ρ_Λ the cosmological density) "unifying" vacuum energy and dark energy in the spirit of a generalized Chaplygin gas model. For $ρ\rightarrow ρ_P$, it reduces to p=-ρc^2 leading to a phase of early accelerated expansion (early inflation) with a constant density equal to the Planck density ρ_P (vacuum energy). For $ρ_Λ\llρ\ll ρ_P$, we recover the standard linear equation of state p=αρc^2 describing radiation (α=1/3) or pressureless matter (α=0) and leading to an intermediate phase of decelerating expansion. For $ρ\rightarrow ρ_Λ$, we get p=-ρc^2 leading to a phase of late accelerated expansion (late inflation) with a constant density equal to the cosmological density ρ_Λ (dark energy). We show a nice symmetry between the early universe (vacuum energy + α-fluid) and the late universe (α-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 ρ_P in the early universe plays a role similar to the cosmological density ρ_Λ 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 (α-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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A parametrization of two-dimensional turbulence based on a maximum entropy production principle with a local conservation of energy

In the context of two-dimensional (2D) turbulence, we apply the maximum entropy production principle (MEPP) by enforcing a local conservation of energy. This leads to an equation for the vorticity distribution that conserves all the Casimirs, the energy, and that increases monotonically the mixing entropy ($H$-theorem). Furthermore, the equation for the coarse-grained vorticity dissipates monotonically all the generalized enstrophies. These equations may provide a parametrization of 2D turbulence. They do not generally relax towards the maximum entropy state. The vorticity current vanishes for any steady state of the 2D Euler equation. Interestingly, the equation for the coarse-grained vorticity obtained from the MEPP turns out to coincide, after some algebraic manipulations, with the one obtained with the anticipated vorticity method. This shows a connection between these two approaches when the conservation of energy is treated locally. Furthermore, the newly derived equation, which incorporates a diffusion term and a drift term, has a nice physical interpretation in terms of a selective decay principle. This gives a new light to both the MEPP and the anticipated vorticity method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamic representation of the Klein-Gordon-Einstein equations in the weak field limit

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 an arbitrary potential of the form $V(|φ|^2)$. 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. Nonrelativistic hydrodynamic equations based on the Schrödinger-Poisson equations or on the Gross-Pitaevskii-Poisson equations are recovered in the limit $c\rightarrow +\infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamic representation of the Klein-Gordon-Einstein equations in the weak field limit: I. General formalism and perturbations analysis

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 $λ|φ|^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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Is the Universe logotropic?

We consider the possibility that the universe is made of a single dark fluid described by a logotropic equation of state $P=A\ln(ρ/ρ_*)$, where $ρ$ is the rest-mass density, $ρ_*$ is a reference density, and $A$ is the logotropic temperature. The energy density $ε$ is the sum of two terms: a rest-mass energy term $ρc^2$ that mimics dark matter and an internal energy term $u(ρ)=-P(ρ)-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/ρ_Λc^2$ where $ρ_Λ$ is the cosmological density. For $B=0$, we recover the $Λ$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 $Σ=ρ_0 r_h$ is constant. This result is in agreement with the observations where it is found that $Σ_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 $ρ_*=ρ_P$, where $ρ_P$ is the Planck density, we predict $B=3.53\times 10^{-3}$, in perfect agreement with the value obtained from the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Secular diffusion in discrete self-gravitating tepid discs I : analytic solution in the tightly wound limit

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Secular diffusion in discrete self-gravitating tepid discs II: accounting for swing amplification via the matrix method

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmology with a stiff matter era

We provide a simple analytical solution of the Friedmann equations for a universe made of stiff matter, dust matter, and dark energy. A stiff matter era is present in the cosmological model of Zel'dovich (1972) where the primordial universe is assumed to be made of a cold gas of baryons. It also occurs in certain cosmological models where dark matter is made of relativistic self-gravitating Bose-Einstein condensates (BECs). When the energy density of the stiff matter is positive, the primordial universe is singular. It starts from a state with a vanishing scale factor and an infinite density. We consider the possibility that the energy density of the stiff matter is negative (anti-stiff matter). This happens, for example, when the BECs have an attractive self-interaction. In that case, the primordial universe is non-singular. It starts from a state in which the scale factor is finite and the energy density is equal to zero. For the sake of generality, we consider a cosmological constant of arbitrary sign. When the cosmological constant is positive, the universe asymptotically reaches a de Sitter phase where the scale factor increases exponentially rapidly. This can account for the accelerating expansion of the universe that we observe at present. When the cosmological constant is negative (anti-de Sitter), the evolution of the universe is cyclic. Therefore, depending on the sign of the energy density of the stiff matter and of the dark energy, we obtain singular and non-singular expanding or cyclic universes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Explicit formula of energy-conserving Fokker-Planck type collision term for single species point vortex systems with weak mean flow

This paper derives a kinetic equation for a two-dimensional single species point vortex system. We consider a situation (different from the ones considered previously) of weak mean flow where the time scale of the macroscopic motion is longer than the decorrelation time so that the trajectory of the point vortices can be approximated by a straight line on the decorrelation time scale. This may be the case when the number $N$ of point vortices is not too large. Using a kinetic theory based on the Klimontovich formalism, we derive a collision term consisting of a diffusion term and a drift term, whose structure is similar to the Fokker-Planck equation. The collision term exhibits several important properties: (a) it includes a nonlocal effect; (b) it conserves the mean field energy; (c) it satisfies the H theorem; (d) its effect vanishes in each local equilibrium region with the same temperature. When the system reaches a global equilibrium state, the collision term completely converges to zero all over the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Models of dark matter halos based on statistical mechanics: I. The classical King model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Models of dark matter halos based on statistical mechanics: II. The fermionic King model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic self-gravitating Bose-Einstein condensates and cold baryons with a stiff equation of state

Because of their superfluid properties, some compact astrophysical objects such as neutron stars may contain a significant part of their matter in the form of a Bose-Einstein condensate (BEC). We consider a partially-relativistic model of self-gravitating BECs where the relation between the pressure and the rest-mass density is assumed to be quadratic (as in the case of classical BECs) but pressure effects are taken into account in the relation between the energy density and the rest-mass density. At high densities, we get a stiff equation of state similar to the one considered by Zel'dovich (1961) in the context of baryon stars in which the baryons interact through a vector meson field. We determine the maximum mass of general relativistic BEC stars described by this equation of state by using the formalism of Tooper (1965). This maximum mass is slightly larger than the maximum mass obtained by Chavanis and Harko (2012) using a fully-relativistic model. We also consider the possibility that dark matter is made of BECs and apply the partially-relativistic model of BECs to cosmology. In this model, we show that the universe experiences a stiff matter phase, followed by a dust matter phase, and finally by a dark energy phase (equivalent to a cosmological constant). The same evolution is obtained in Zel'dovich (1972) model which assumes that initially, near the cosmological singularity, the universe is filled with cold baryons. Interestingly, the Friedmann equations can be solved analytically in that case and provide a simple generalization of the $Λ$CDM model. We point out, however, the limitations of the partially-relativistic model for BECs and show the need for a fully-relativistic one.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical mechanics of two-dimensional point vortices: relaxation equations and strong mixing limit

We complete the literature on the statistical mechanics of point vortices in two-dimensional hydrodynamics. Using a maximum entropy principle, we determine the multi-species Boltzmann-Poisson equation and establish a form of virial theorem. Using a maximum entropy production principle (MEPP), we derive a set of relaxation equations towards statistical equilibrium. These relaxation equations can be used as a numerical algorithm to compute the maximum entropy state. We mention the analogies with the Fokker-Planck equations derived by Debye and Hückel for electrolytes. We then consider the limit of strong mixing (or low energy). To leading order, the relationship between the vorticity and the stream function at equilibrium is linear and the maximization of the entropy becomes equivalent to the minimization of the enstrophy. This expansion is similar to the Debye-Hückel approximation for electrolytes, except that the temperature is negative instead of positive so that the effective interaction between like-sign vortices is attractive instead of repulsive. This leads to an organization at large scales presenting geometry-induced phase transitions, instead of Debye shielding. We compare the results obtained with point vortices to those obtained in the context of the statistical mechanics of continuous vorticity fields described by the Miller-Robert-Sommeria (MRS) theory. At linear order, we get the same results but differences appear at the next order. In particular, the MRS theory predicts a transition between sinh and tanh-like ω-ψrelationships depending on the sign of Ku-3 (where Ku is the Kurtosis) while there is no such transition for point vortices which always show a sinh-like ω-ψrelationship. We derive the form of the relaxation equations in the strong mixing limit and show that the enstrophy plays the role of a Lyapunov functional.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The fermionic King model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational phase transitions with an exclusion constraint in position space

We discuss the statistical mechanics of a system of self-gravitating particles with an exclusion constraint in position space in a space of dimension $d$. The exclusion constraint puts an upper bound on the density of the system and can stabilize it against gravitational collapse. We plot the caloric curves giving the temperature as a function of the energy and investigate the nature of phase transitions as a function of the size of the system and of the dimension of space 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. For $d\le 2$, there is no phase transition. For $d>2$, phase transitions can take place between a "gaseous" phase unaffected by the exclusion constraint and a "condensed" phase dominated by this constraint. The condensed configurations have a core-halo structure made of a "rocky core" surrounded by an "atmosphere", similar to a giant gaseous planet. For large systems there exist microcanonical and canonical first order phase transitions. For intermediate systems, only canonical first order phase transitions are present. For small systems there is no phase transition at all. As a result, the phase diagram exhibits two critical points, one in each ensemble. There also exist a region of negative specific heats and a situation of ensemble inequivalence for sufficiently large systems. By a proper interpretation of the parameters, our results have application for the chemotaxis of bacterial populations in biology described by a generalized Keller-Segel model including an exclusion constraint in position space. They also describe colloids at a fluid interface driven by attractive capillary interactions when there is an excluded volume around the particles. Connexions with two-dimensional turbulence are also mentioned.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Initial value problem for the linearized mean field Kramers equation with long-range interactions

We solve the initial value problem for the linearized mean field Kramers equation describing Brownian particles with long-range interactions in the $N\rightarrow +\infty$ limit. We show that the dielectric function can be expressed in terms of incomplete Gamma functions. The dielectric functions associated with the linearized Vlasov equation and with the linearized mean field Smoluchowski equation are recovered as special cases corresponding to the no friction limit or to the strong friction limit respectively. Although the stability of the Maxwell-Boltzmann distribution is independent on the friction parameter, the evolution of the perturbation depends on it in a non-trivial manner. For illustration, we apply our results to self-gravitating systems, plasmas, and to the attractive and repulsive BMF models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic theory of spatially homogeneous systems with long-range interactions: II. Basic equations

We provide a short historic of the early development of kinetic theory in plasma physics and synthesize the basic kinetic equations describing the evolution of systems with long-range interactions derived in Paper I. We describe the evolution of the system as a whole and the relaxation of a test particle in a bath at equilibrium or out-of-equilibrium. We write these equations for an arbitrary long-range potential of interaction in a space of arbitrary dimension d. We discuss the scaling of the relaxation time with the number of particles for non-singular potentials. For always spatially homogeneous systems, the relaxation time of the system as a whole scales like N in d>1 and like N^2 (presumably) in d=1. For always spatially inhomogeneous systems, the relaxation time of the system as a whole scales like N in any dimension of space. For one dimensional systems undergoing a dynamical phase transition from a homogeneous to an inhomogeneous phase, we expect a relaxation time of the form N^δ with 1<δ<2 intermediate between the two previous cases. The relaxation time of a test particle in a bath always scales like N. We also briefly discuss the kinetic theory of systems with long-range interactions submitted to an external stochastic potential. This paper gathers basic equations that are applied to specific systems in Paper III.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic theory of spatially homogeneous systems with long-range interactions: III. Application to power-law potentials, plasmas, stellar systems, and to the HMF model

We apply the general results of the kinetic theory of systems with long-range interactions to particular systems of physical interest. We consider repulsive and attractive power-law potentials of interaction r^{-γ} with γ<d in a space of dimension d. For γ>γ_c= (d-1)/2, strong collisions must be taken into account and the evolution of the system is governed by the Boltzmann equation or by a modified Landau equation; for γ<γ_c, strong collisions are negligible and the evolution of the system is governed by the Lenard-Balescu equation. In the marginal case γ=γ_c, we can use the Landau equation (with appropriately justified cut-offs) as a relevant approximation of the Boltzmann and Lenard-Balescu equations. The divergence at small scales that appears in the original Landau equation is regularized by the effect of strong collisions. In the case of repulsive interactions with a neutralizing background (e.g. plasmas), the divergence at large scales that appears in the original Landau equation is regularized by collective effects accounting for Debye shielding. In the case of attractive interactions (e.g. gravity), it is regularized by the spatial inhomogeneity of the system and its finite extent. We provide explicit analytical expressions of the diffusion and friction coefficients, and of the relaxation time, depending on the value of the exponent γand on the dimension of space d. We treat in a unified framework the case of Coulombian plasmas and stellar systems in various dimensions of space, and the case of the attractive and repulsive HMF models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Brownian Mean Field model

We discuss the dynamics and thermodynamics of the Brownian Mean Field (BMF) model which is a system of N Brownian particles moving on a circle and interacting via a cosine potential. It can be viewed as the canonical version of the Hamiltonian Mean Field (HMF) model. We first complete the description of this system in the mean field approximation. Then, we take fluctuations into account and study the stochastic evolution of the magnetization both in the homogeneous phase and in the inhomogeneous phase. We discuss its behavior close to the critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A simple model of universe with a polytropic equation of state

We construct a simple model of universe with a generalized equation of state $p=(α+kρ^{1/n})ρc^2$ having a linear component $p=αρc^2$ and a polytropic component $p=kρ^{1+1/n}c^2$. For $α=1/3$, $n=1$ and $k=-4/(3ρ_P)$, where $ρ_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 $α=0$, $n=-1$ and $k=-ρ_Λ$, where $ρ_Λ=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_Λ=1.46 10^{18} s$. This polytropic model reveals a nice "symmetry" between the early and late evolution of the universe, the cosmological constant $Λ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinetic theory of spatially inhomogeneous stellar systems without collective effects

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Linear response theory for hydrodynamic and kinetic equations with long-range interactions

We apply the linear response theory to systems with long-range interactions described by hydrodynamic equations such as the Euler, Smoluchowski, and damped Euler equations. We analytically determine the response of the system submitted to a pulse and to a step function. We compare these results with those obtained for collisionless systems described by the Vlasov equation. We show that, in the linear regime, the evolution of a collisionless system (Vlasov) with the waterbag distribution is the same as the evolution of a collision-dominated gas without dissipation (Euler). In this analogy, the maximum velocity of the waterbag distribution plays the role of the velocity of sound in the corresponding barotropic gas. When submitted to a step function, these systems exhibit permanent oscillations. Other distributions exhibit Landau damping and relax towards a steady state. We illustrate this behaviour with the Cauchy distribution which can be studied analytically. We apply our results to the HMF model and obtain a generalized Curie-Weiss law for the magnetic susceptibility. Finally, we compare the linear response theory to the initial value problem for the linearized Vlasov equation and report a case of algebraic damping of the initial perturbation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Models of universe with a polytropic equation of state: I. The early universe

We construct models of universe with a generalized equation of state $p=(αρ+kρ^{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 $α=1/3$, $n=1$ and $k=-4/(3ρ_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 $ρ_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 $α=1/3$, $n=1$ and $k=4/(3ρ_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$).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Models of universe with a polytropic equation of state: II. The late universe

We construct models of universe with a generalized equation of state $p=(αρ+kρ^{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 $α=0$, $n=-1$ and $k=-ρ_Λ$, 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 $ρ_Λ=7.02 10^{-24} g/m^3$ on a timescale $t_Λ=1.46 10^{18} s$. For $α=0$, $n=-1$ and $k=ρ_Λ$, 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=(αρ+kρ^{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ρ^2/3ρ_P+ρ/3-4ρ_Λ/3$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Models of universe with a polytropic equation of state: III. The phantom universe

We construct models of universe with a generalized equation of state $p=(αρ+kρ^{1+1/n})c^2$ having a linear component and a polytropic component. The linear equation of state $p=αρc^2$ with $-1\le α\le 1$ describes radiation ($α=1/3$), pressureless matter ($α=0$), stiff matter ($α=1$), and vacuum energy ($α=-1$). The polytropic equation of state $p=kρ^{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/ρ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 $α>-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=(αρ+kρ^{1+1/n})c^2$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase transitions and marginal ensemble equivalence for freely evolving flows on a rotating sphere

The large-scale circulation of planetary atmospheres like that of the Earth is traditionally thought of in a dynamical framework. Here, we apply the statistical mechanics theory of turbulent flows to a simplified model of the global atmosphere, the quasi-geostrophic model, leading to non-trivial equilibria. Depending on a few global parameters, the structure of the flow may be either a solid-body rotation (zonal flow) or a dipole. A second order phase transition occurs between these two phases, with associated spontaneous symmetry-breaking in the dipole phase. This model allows us to go beyond the general theory of marginal ensemble equivalence through the notion of Goldstone modes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Statistical mechanics of quasi-geostrophic flows on a rotating sphere

Statistical mechanics provides an elegant explanation to the appearance of coherent structures in two-dimensional inviscid turbulence: while the fine-grained vorticity field, described by the Euler equation, becomes more and more filamented through time, its dynamical evolution is constrained by some global conservation laws (energy, Casimir invariants). As a consequence, the coarse-grained vorticity field can be predicted through standard statistical mechanics arguments (relying on the Hamiltonian structure of the two-dimensional Euler flow), for any given set of the integral constraints. It has been suggested that the theory applies equally well to geophysical turbulence; specifically in the case of the quasi-geostrophic equations, with potential vorticity playing the role of the advected quantity. In this study, we demonstrate analytically that the Miller-Robert-Sommeria theory leads to non-trivial statistical equilibria for quasi-geostrophic flows on a rotating sphere, with or without bottom topography. We first consider flows without bottom topography and with an infinite Rossby deformation radius, with and without conservation of angular momentum. When the conservation of angular momentum is taken into account, we report a case of second order phase transition associated with spontaneous symmetry breaking. In a second step, we treat the general case of a flow with an arbitrary bottom topography and a finite Rossby deformation radius. Previous studies were restricted to flows in a planar domain with fixed or periodic boundary conditions with a beta-effect. In these different cases, we are able to classify the statistical equilibria for the large-scale flow through their sole macroscopic features. We build the phase diagrams of the system and discuss the relations of the various statistical ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Virial theorem for Onsager vortices in two-dimensional hydrodynamics

We derive the virial theorem appropriate to two-dimensional point vortices at statistical equilibrium in the microcanonical and canonical ensembles. In an unbounded domain, it relates the angular velocity to the angular momentum and the temperature. Our expression is valid for an arbitrary number of point vortices of possibly different species. In the single-species case, and in the mean field approximation, it reduces to the relation empirically obtained by J.H. Williamson [J. Plasma Physics 17, 85 (1977)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

BEC dark matter, Zeldovich approximation and generalized Burgers equation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effective merging dynamics of two and three fluid vortices: Application to two-dimensional decaying turbulence

We present a kinetic theory of two-dimensional decaying turbulence in the context of two-body and three-body vortex merging processes. By introducing the equations of motion for two or three vortices in the effective noise due to all the other vortices, we demonstrate analytically that a two-body mechanism becomes inefficient at low vortex density $n\ll 1$. When the more efficient three-body vortex mergings are considered {(involving vortices of different signs)}, we show that $n\sim t^{-ξ}$, with $ξ=1$. We generalize this argument to three-dimensional geostrophic turbulence, finding $ξ=5/4$, in excellent agreement with direct Navier-Stokes simulations [J.\,C. McWilliams \emph{et al.}, J. Fluid Mech. {\bf 401}, 1 (1999)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exact analytical solution of the collapse of self-gravitating Brownian particles and bacterial populations at zero temperature

We provide an exact analytical solution of the collapse dynamics of self-gravitating Brownian particles and bacterial populations at zero temperature. These systems are described by the Smoluchowski-Poisson system or Keller-Segel model in which the diffusion term is neglected. As a result, the dynamics is purely deterministic. A cold system undergoes a gravitational collapse leading to a finite time singularity: the central density increases and becomes infinite in a finite time t_coll. The evolution continues in the post collapse regime. A Dirac peak emerges, grows and finally captures all the mass in a finite time t_end, while the central density excluding the Dirac peak progressively decreases. Close to the collapse time, the pre and post collapse evolution is self-similar. Interestingly, if one starts from a parabolic density profile, one obtains an exact analytical solution that describes the whole collapse dynamics, from the initial time to the end, and accounts for non self-similar corrections that were neglected in previous works. Our results have possible application in different areas including astrophysics, chemotaxis, colloids and nanoscience.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Growth of perturbations in an expanding universe with Bose-Einstein condensate dark matter

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 Λ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=(αρ+ k ρ^2)c^2 having a component p=k ρ^2 c^2 similar to a BEC dark matter and a component p=αρc^2 mimicking the effect of the cosmological constant (dark energy). We find optimal parameters that give a good agreement with the standard ΛCDM model assuming a finite cosmological constant.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Instability of a uniformly collapsing cloud of classical and quantum self-gravitating Brownian particles

We study the growth of perturbations in a uniformly collapsing cloud of self-gravitating Brownian particles. This problem shares analogies with the formation of large-scale structures in a universe experiencing a "big-crunch" or with the formation of stars in a molecular cloud experiencing gravitational collapse. Starting from the barotropic Smoluchowski-Poisson system, we derive a new equation describing the evolution of the density contrast in the comoving (collapsing) frame. This equation can serve as a prototype to study the process of self-organization in complex media with structureless initial conditions. We solve this equation analytically in the linear regime and compare the results with those obtained by using the "Jeans swindle" in a static medium. The stability criteria, as well as the laws for the time evolution of the perturbations, are different. The Jeans criterion is expressed in terms of a critical wavelength $λ_J$ while our criterion is expressed in terms of a critical polytropic index $γ_{4/3}$. We also study the fragmentation process in the nonlinear regime. We determine the growth of the skewness, the long-wavelength tail of the power spectrum and find a self-similar solution to the nonlinear equations valid for large times. Finally, we consider dissipative self-gravitating Bose-Einstein condensates with short-range interactions and show that, in a strong friction limit, the dissipative Gross-Pitaevskii-Poisson system is equivalent to the quantum barotropic Smoluchowski-Poisson system. This yields a new type of nonlinear mean field Fokker-Planck equations including quantum effects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic theory of long-range interacting systems with angle-action variables and collective effects

We develop a kinetic theory of systems with long-range interactions taking collective effects and spatial inhomogeneity into account. Starting from the Klimontovich equation and using a quasilinear approximation, we derive a Lenard-Balescu-type kinetic equation written in angle-action variables. We confirm the result obtained by Heyvaerts [Mon. Not. R. Astron. Soc. {\bf 407}, 355 (2010)] who started from the Liouville equation and used the BBGKY hierarchy truncated at the level of the two-body distribution function. When collective effects are neglected, we recover the Landau-type kinetic equation obtained in our previous papers [P.H. Chavanis, Physica A {\bf 377}, 469 (2007); J. Stat. Mech., P05019 (2010)]. We also consider the relaxation of a test particle in a bath of field particles. Its stochastic motion is described by a Fokker-Planck equation written in angle-action variables. We determine the diffusion tensor and the friction force by explicitly calculating the first and second order moments of the increment of action of the test particle from its equations of motion, taking collective effects into account. This generalizes the expressions obtained in our previous works. We discuss the scaling with $N$ of the relaxation time for the system as a whole and for a test particle in a bath.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic theory of Onsager's vortices in two-dimensional hydrodynamics

Starting from the Liouville equation, and using a BBGKY-like hierarchy, we derive a kinetic equation for the point vortex gas in two-dimensional (2D) hydrodynamics, taking two-body correlations and collective effects into account. This equation is valid at the order 1/N where N>>1 is the number of point vortices in the system (we assume that their individual circulation scales like γ~ 1/N). It gives the first correction, due to graininess and correlation effects, to the 2D Euler equation that is obtained for $N\rightarrow +\infty$. For axisymmetric distributions, this kinetic equation does not relax towards the Boltzmann distribution of statistical equilibrium. This implies either that (i) the "collisional" (correlational) relaxation time is larger than Nt_D, where t_D is the dynamical time, so that three-body, four-body... correlations must be taken into account in the kinetic theory, or (ii) that the point vortex gas is non-ergodic (or does not mix well) and will never attain statistical equilibrium. Non-axisymmetric distributions may relax towards the Boltzmann distribution on a timescale of the order Nt_D due to the existence of additional resonances, but this is hard to prove from the kinetic theory. On the other hand, 2D Euler unstable vortex distributions can experience a process of "collisionless" (correlationless) violent relaxation towards a non-Boltzmannian quasistationary state (QSS) on a very short timescale of the order of a few dynamical times. This QSS is possibly described by the Miller-Robert-Sommeria (MRS) statistical theory which is the counterpart, in the context of two-dimensional hydrodynamics, of the Lynden-Bell statistical theory of violent relaxation in stellar dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic theory of spatially homogeneous systems with long-range interactions: I. General results

We review and complete the existing literature on the kinetic theory of spatially homogeneous systems with long-range interactions taking collective effects into account. The evolution of the system as a whole is described by the Lenard-Balescu equation which is valid in a weak coupling approximation. When collective effects are neglected it reduces to the Landau equation and when collisions (correlations) are neglected it reduces to the Vlasov equation. The relaxation of a test particle in a bath is described by a Fokker-Planck equation involving a diffusion term and a friction term. For a thermal bath, the diffusion and friction coefficients are connected by an Einstein relation. General expressions of the diffusion and friction coefficients are given, depending on the potential of interaction and on the dimension of space. We also discuss the scaling with $N$ (number of particles) or with $Λ$ (plasma parameter) of the relaxation time towards statistical equilibrium. Finally, we consider the effect of an external stochastic forcing on the evolution of the system and compute the corresponding term in the kinetic equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic theory of two-dimensional point vortices with collective effects

We develop a kinetic theory of point vortices in two-dimensional hydrodynamics taking collective effects into account. We first recall the approach of Dubin & O'Neil [Phys. Rev. Lett. 60, 1286 (1988)] that leads to a Lenard-Balescu-type kinetic equation for axisymmetric flows. When collective effects are neglected, it reduces to the Landau-type kinetic equation obtained independently in our previous papers [P.H. Chavanis, Phys. Rev. E 64, 026309 (2001); Physica A 387, 1123 (2008)]. We also consider the relaxation of a test vortex in a "sea" (bath) of field vortices. Its stochastic motion is described in terms of a Fokker-Planck equation. We determine the diffusion coefficient and the drift term by explicitly calculating the first and second order moments of the radial displacement of the test vortex from its equations of motion, taking collective effects into account. This generalizes the expressions obtained in our previous works. We discuss the scaling with N of the relaxation time for the system as a whole and for a test vortex in a bath.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mass-radius relation of Newtonian self-gravitating Bose-Einstein condensates with short-range interactions: I. Analytical results

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öhmer & 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Brownian particles with long and short range interactions

We develop a kinetic theory of Brownian particles with long and short range interactions. We consider both overdamped and inertial models. In the overdamped limit, the evolution of the spatial density is governed by the generalized mean field Smoluchowski equation including a mean field potential due to long-range interactions and a generically nonlinear barotropic pressure due to short-range interactions. This equation describes various physical systems such as self-gravitating Brownian particles (Smoluchowski-Poisson system), bacterial populations experiencing chemotaxis (Keller-Segel model) and colloidal particles with capillary interactions. We also take into account the inertia of the particles and derive corresponding kinetic and hydrodynamic equations generalizing the usual Kramers, Jeans, Euler and Cattaneo equations. For each model, we provide the corresponding form of free energy and establish the H-theorem and the virial theorem. Finally, we show that the same hydrodynamic equations are obtained in the context of nonlinear mean field Fokker-Planck equations associated with generalized thermodynamics. However, in that case, the nonlinear pressure is due to the bias in the transition probabilities from one state to the other leading to non-Boltzmannian distributions while in the former case the distribution is Boltzmannian but the nonlinear pressure arises from the two-body correlation function induced by the short-range potential of interaction. As a whole, our paper develops connections between the topics of long-range interactions, short-range interactions, nonlinear mean field Fokker-Planck equations and generalized thermodynamics. It also justifies from a kinetic theory based on microscopic processes, the basic equations that were introduced phenomenologically in gravitational Brownian dynamics, chemotaxis and colloidal suspensions with attractive interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Collapse and evaporation of a canonical self-gravitating gas

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(ρ)=Tρ^{1/n}$, where $ρ$ is the local particle density, so that the equation of state is $P(ρ)=D(ρ)ρ$. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical and thermodynamical stability of isothermal distributions in the HMF model

We provide a new derivation of the conditions of dynamical and thermodynamical stability of homogeneous and inhomogeneous isothermal distributions in the Hamiltonian Mean Field (HMF) model. This proof completes the original thermodynamical approach of Inagaki [Prog. Theor. Phys. 90, 557 (1993)]. Our formalism, based on variational principles, is simple and the method can be applied to more general situations. For example, it can be used to settle the dynamical stability of polytropic distributions with respect to the Vlasov equation [Chavanis & Campa, arXiv:1001.2109]. For isothermal distributions, the calculations can be performed fully analytically, providing therefore a clear illustration of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical stability criterion for inhomogeneous quasi-stationary states in long-range systems

We derive a necessary and sufficient condition of linear dynamical stability for inhomogeneous Vlasov stationary states of the Hamiltonian Mean Field (HMF) model. The condition is expressed by an explicit disequality that has to be satisfied by the stationary state, and it generalizes the known disequality for homogeneous stationary states. In addition, we derive analogous disequalities that express necessary and sufficient conditions of formal stability for the stationary states. Their usefulness, from the point of view of linear dynamical stability, is that they are simpler, although they provide only sufficient criteria of linear stability. We show that for homogeneous stationary states the relations become equal, and therefore linear dynamical stability and formal stability become equivalent.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical stability of infinite homogeneous self-gravitating systems: application of the Nyquist method

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Escape of stars from gravitational clusters in the Chandrasekhar model

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 $λ(R)$, and the corresponding eigenfunction $f_λ(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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kinetic equations for systems with long-range interactions: a unified description

We complete the existing literature on the kinetic theory of systems with long-range interactions. Starting from the BBGKY hierarchy, or using projection operator technics or a quasilinear theory, a general kinetic equation can be derived when collective effects are neglected. This equation (which is not well-known) applies to possibly spatially inhomogeneous systems, which is specific to systems with long-range interactions. Interestingly, the structure of this kinetic equation bears a clear physical meaning in terms of generalized Kubo relations. Furthermore, this equation takes a very similar form for stellar systems and two-dimensional point vortices providing therefore a unified description of the kinetic theory of these systems. If we assume that the system is spatially homogeneous (or axisymmetric for point vortices), this equation can be simplified and reduces to the Landau equation (or its counterpart for point vortices). Our formalism thus offers a simple derivation of Landau-type equations. We also use this general formalism to derive a kinetic equation, written in angle-action variables, describing spatially inhomogeneous systems with long-range interactions. This new derivation solves the shortcomings of our previous derivation [P.H. Chavanis, Physica A 377, 469 (2007)]. Finally, we consider a test particle approach and derive general expressions for the diffusion and friction (or drift) coefficients of a test particle evolving in a bath of field particles. We make contact with the expressions previously obtained in the literature. As an application of the kinetic theory, we argue that the relaxation time is shorter for inhomogeneous (or high-dimensional) systems than for homogeneous (or low-dimensional) systems because there are potentially more resonances.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Noise-induced dynamical phase transitions in long-range systems

In the thermodynamic limit, the time evolution of isolated long-range interacting systems is properly described by the Vlasov equation. This equation admits non-equilibrium dynamically stable stationary solutions characterized by a zero order parameter. We show that the presence of external noise sources, like for instance a heat bath, can induce at a specific time a dynamical phase transition marked by a non-zero order parameter. This transition corresponds to a restoring of the full ergodic properties of the system and may be used as a distinctive experimental signature of the existence of non-equilibrium Vlasov-stable states. In particular, we evidence for the first time a regime characterized by an order parameter pulse. Our analytical results are corroborated by numerical simulations of a paradigmatic long-range model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical mechanics of Beltrami flows in axisymmetric geometry: Equilibria and bifurcations

We characterize the thermodynamical equilibrium states of axisymmetric Euler-Beltrami flows. They have the form of coherent structures presenting one or several cells. We find the relevant control parameters and derive the corresponding equations of state. We prove the coexistence of several equilibrium states for a given value of the control parameter like in 2D turbulence [Chavanis and Sommeria, J. Fluid Mech. 314, 267 (1996)]. We explore the stability of these equilibrium states and show that all states are saddle points of entropy and can, in principle, be destabilized by a perturbation with a larger wavenumber, resulting in a structure at the smallest available scale. This mechanism is therefore reminiscent of the 3D Richardson energy cascade towards smaller and smaller scales. Therefore, our system is truly intermediate between 2D turbulence (coherent structures) and 3D turbulence (energy cascade). We further explore numerically the robustness of the equilibrium states with respect to random perturbations using a relaxation algorithm in both canonical and microcanonical ensembles. We show that saddle points of entropy can be very robust and therefore play a role in the dynamics. We evidence differences in the robustness of the solutions in the canonical and microcanonical ensembles. A scenario of bifurcation between two different equilibria (with one or two cells) is proposed and discussed in connection with a recent observation of a turbulent bifurcation in a von Karman experiment [Ravelet et al., Phys. Rev. Lett. 93, 164501 (2004)].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical mechanics of Beltrami flows in axisymmetric geometry: Theory reexamined

A simplified thermodynamic approach of the incompressible axisymmetric Euler equations is considered based on the conservation of helicity, angular momentum and microscopic energy. Statistical equilibrium states are obtained by maximizing the Boltzmann entropy under these sole constraints. We assume that these constraints are selected by the properties of forcing and dissipation. The fluctuations are found to be Gaussian while the mean flow is in a Beltrami state. Furthermore, we show that the maximization of entropy at fixed helicity, angular momentum and microscopic energy is equivalent to the minimization of macroscopic energy at fixed helicity and angular momentum. This provides a justification of this selective decay principle from statistical mechanics. These theoretical predictions are in good agreement with experiments of a von Karman turbulent flow and provide a way to measure the temperature of turbulence and check Fluctuation-Dissipation Relations (FDR). Relaxation equations are derived that could provide an effective description of the dynamics towards the Beltrami state and the progressive emergence of a Gaussian distribution. They can also provide a numerical algorithm to determine maximum entropy states or minimum energy states.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The HMF model for fermions and bosons

We study the thermodynamics of quantum particles with long-range interactions at T=0. Specifically, we generalize the Hamiltonian Mean Field (HMF) model to the case of fermions and bosons. In the case of fermions, we consider the Thomas-Fermi approximation that becomes exact in a proper thermodynamic limit. The equilibrium configurations, described by the Fermi (or waterbag) distribution, are equivalent to polytropes with index n=1/2. In the case of bosons, we consider the Hartree approximation that becomes exact in a proper thermodynamic limit. The equilibrium configurations are solutions of the mean field Schrödinger equation with a cosine interaction. We show that the homogeneous phase, that is unstable in the classical regime, becomes stable in the quantum regime. This takes place through a first order phase transition for fermions and through a second order phase transition for bosons where the control parameter is the normalized Planck constant. In the case of fermions, the homogeneous phase is stabilized by the Pauli exclusion principle while for bosons the stabilization is due to the Heisenberg uncertainty principle. As a result, the thermodynamic limit is different for fermions and bosons. We point out analogies between the quantum HMF model and the concepts of fermion and boson stars in astrophysics. Finally, as a by-product of our analysis, we obtain new results concerning the Vlasov dynamical stability of the waterbag distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermodynamics of the HMF model with a magnetic field

We study the thermodynamics of the Hamiltonian Mean Field (HMF) model with an external potential playing the role of a "magnetic field". If we consider only fully stable states, this system does not present any phase transition. However, if we take into account metastable states (for a restricted class of perturbations), we find a very rich phenomenology. In particular, the system displays a region of negative specific heats in the microcanonical ensemble in which the temperature decreases as the energy increases. This leads to ensembles inequivalence and to zeroth order phase transitions similar to the "gravothermal catastrophe" and to the "isothermal collapse" of self-gravitating systems. In the present case, they correspond to the reorganization of the system from an "anti-aligned" phase (magnetization pointing in the direction opposite to the magnetic field) to an "aligned" phase (magnetization pointing in the same direction as the magnetic field). We also find that the magnetic susceptibility can be negative in the microcanonical ensemble so that the magnetization decreases as the magnetic field increases. The magnetic curves can take various shapes depending on the values of energy or temperature. We describe hysteretic cycles involving positive or negative susceptibilities. We also show that this model exhibits gaps in the magnetization at fixed energy, resulting in ergodicity breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Statistics of the gravitational force in various dimensions of space: from Gaussian to Levy laws

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Virial theorem for rotating self-gravitating Brownian particles and two-dimensional point vortices

We derive the proper form of Virial theorem for a system of rotating self-gravitating Brownian particles. We show that, in the two-dimensional case, it takes a very simple form that can be used to obtain general results about the dynamics of the system without being required to solve the Smoluchowski-Poisson system explicitly. We also develop the analogy between self-gravitating systems and two-dimensional point vortices and derive a Virial-like relation for the vortex system.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Generalized Keller-Segel models of chemotaxis. Analogy with nonlinear mean field Fokker-Planck equations

We consider a generalized class of Keller-Segel models describing the chemotaxis of biological populations (bacteria, amoebae, endothelial cells, social insects,...). We show the analogy with nonlinear mean field Fokker-Planck equations and generalized thermodynamics. As an illustration, we introduce a new model of chemotaxis incorporating both effects of anomalous diffusion and exclusion principle (volume filling). We also discuss the analogy between biological populations described by the Keller-Segel model and self-gravitating Brownian particles described by the Smoluchowski-Poisson system.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Kinetic theory with angle-action variables

We present a kinetic theory for inhomogeneous systems with weak long-range interactions. Starting from the Klimontovich equation and using a quasilinear theory, we obtain a closed kinetic equation describing the relaxation of the distribution function of the system as a whole due to resonances between different orbits. This equation is written in angle-action variables. It conserves mass and energy and increases the Boltzmann entropy (H-theorem). Using a thermal bath approximation, we derive a Fokker-Planck equation that describes the relaxation of a test particle towards the Boltzmann distribution under the combined effect of diffusion and friction terms. We mention some analogies with the kinetic theory of point vortices in two-dimensional hydrodynamics. We also stress the limitations of our approach and the connection with recent works.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Out of equilibrium phase transitions in mean field Hamiltonian dynamics

Systems with long-range interactions display a short-time relaxation towards Quasi-Stationary States (QSSs), whose lifetime increases with system size. With reference to the Hamiltonian Mean Field (HMF) model, we here review Lynden-Bell's theory of ``violent relaxation''. The latter results in a maximum entropy scheme for a water-bag initial profile which predicts the presence of out-of-equilibrium phase transitions} separating homogeneous (zero magnetization) from inhomogeneous (non-zero magnetization) QSSs. Two different parametric representations of the initial condition are analyzed and the features of the phase diagram are discussed. In both representations we find a second order and a first order line of phase transitions that merge at a tricritical point. Particular attention is payed to the condition of existence and stability of the homogenous phase.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Dynamics and thermodynamics of axisymmetric flows: I. Theory

We develop new variational principles to study stability and equilibrium of axisymmetric flows. We show that there is an infinite number of steady state solutions. We show that these steady states maximize a (non-universal) $H$-function. We derive relaxation equations which can be used as numerical algorithm to construct stable stationary solutions of axisymmetric flows. In a second part, we develop a thermodynamical approach to the equilibrium states at some fixed coarse-grained scale. We show that the resulting distribution can be divided in a universal part coming from the conservation of robust invariants and one non-universal determined by the initial conditions through the fragile invariants (for freely evolving systems) or by a prior distribution encoding non-ideal effects such as viscosity, small-scale forcing and dissipation (for forced systems). Finally, we derive a parameterization of inviscid mixing to describe the dynamics of the system at the coarse-grained scale.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Virial theorem and dynamical evolution of self-gravitating Brownian particles and bacterial populations in an unbounded domain

We derive the Virial theorem appropriate to the generalized Smoluchowski-Poisson system describing self-gravitating Brownian particles and bacterial populations (chemotaxis). We extend previous works by considering the case of an unbounded domain and an arbitrary equation of state. We use the Virial theorem to study the diffusion (evaporation) of an isothermal Brownian gas above the critical temperature T_c in dimension d=2 and show how the effective diffusion coefficient and the Einstein relation are modified by self-gravity or chemotactic attraction. We also study the collapse at T=T_c and show that the central density increases logarithmically with time instead of exponentially in a bounded domain. Finally, for d>2, we show that the evaporation of the system is essentially a pure diffusion slightly slowed-down by self-gravity. We also study the linear dynamical stability of stationary solutions of the generalized Smoluchowski-Poisson system representing isolated clusters of particles and investigate the influence of the equation of state and of the dimension of space on the dynamical stability of the system. Finally, we propose a general kinetic and hydrodynamic description of self-gravitating Brownian particles and biological populations and recover known models in some particular limits.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Generalized thermodynamics and kinetic equations: Boltzmann, Landau, Kramers and Smoluchowski

We propose a formal extension of thermodynamics and kinetic theories to a larger class of entropy functionals. Kinetic equations associated to Boltzmann, Fermi, Bose and Tsallis entropies are recovered as a special case. This formalism first provides a unifying description of classical and quantum kinetic theories. On the other hand, a generalized thermodynamical framework is justified to describe complex systems exhibiting anomalous diffusion. Finally, a notion of generalized thermodynamics emerges in the context of the the violent relaxation of collisionless stellar systems and two-dimensional vortices due to the existence of Casimir invariants and incomplete relaxation. A thermodynamical analogy can also be developed to analyze the nonlinear dynamical stability of stationary solutions of the Vlasov and 2D Euler-Poisson systems. On general grounds, we suggest that generalized entropies arise due to the existence of ``hidden constraints'' that modify the form of entropy that we would naively expect. Generalized kinetic equations are therefore ``effective'' equations that are introduced heuristically to describe complex systems.

0 citations0 reviewsNo code linkedOpen access