Source author record

P. H. Chavanis

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

13 published item(s)

preprint2014arXiv

Random transitions described by the stochastic Smoluchowski-Poisson system and by the stochastic Keller-Segel model

We study random transitions between two metastable states that appear below a critical temperature in a one dimensional self-gravitating Brownian gas with a modified Poisson equation experiencing a second order phase transition from a homogeneous phase to an inhomogeneous phase. We numerically solve the $N$-body Langevin equations and the stochastic Smoluchowski-Poisson system which takes fluctuations (finite $N$ effects) into account. The system switches back and forth between the two metastable states (bistability) and the particles accumulate successively at the center or at the boundary of the domain. We show that these random transitions exhibit the phenomenology of the ordinary Kramers problem for a Brownian particle in a double-well potential. The distribution of the residence time is Poissonian and the average lifetime of a metastable state is given by the Arrhenius law, i.e. it is proportional to the exponential of the barrier of free energy $ΔF$ divided by $k_B T$. Since the free energy is proportional to the number of particles $N$ for a system with long-range interactions, the lifetime of metastable states scales as $e^N$ and is considerable for $N\gg 1$. As a result, in many applications, metastable states of systems with long-range interactions can be considered as stable states. However, for moderate values of $N$, or close to a critical point $T_c$, the lifetime of the metastable states is reduced since the barrier of free energy decreases. In that case, the fluctuations become important and the mean field approximation is no more valid. This is the situation considered in this paper. By an appropriate change of notations, our results also apply to bacterial populations experiencing chemotaxis in biology. Their dynamics can be described by a stochastic Keller-Segel model that takes fluctuations into account and goes beyond the usual mean field approximation.

preprint2012arXiv

Caloric curves fitted by polytropic distributions in the HMF model

We perform direct numerical simulations of the HMF model starting from non-magnetized initial conditions with a velocity distribution that is (i) gaussian, (ii) semi-elliptical, and (iii) waterbag. Below a critical energy E_c, depending on the initial condition, this distribution is Vlasov dynamically unstable. The system undergoes a process of violent relaxation and quickly reaches a quasi-stationary state (QSS). We find that the distribution function of this QSS can be conveniently fitted by a polytrope with index (i) n=2, (ii) n=1, and (iii) n=1/2. Using the values of these indices, we are able to determine the physical caloric curve T_{kin}(E) and explain the negative kinetic specific heat region C_{kin}=dE/dT_{kin}<0 observed in the numerical simulations. At low energies, we find that the system takes a "core-halo" structure. The core corresponds to the pure polytrope discussed above but it is now surrounded by a halo of particles. We also consider unsteady initial conditions with magnetization M_0=1 and isotropic waterbag distribution and report the complex dynamics of the system creating phase space holes and dense filaments. We show that the kinetic caloric curve is approximately constant, corresponding to a polytrope with index n_0= 3.56. Finally, we consider the collisional evolution of an initially Vlasov stable distribution, and show that the time-evolving distribution function f(v,t) can be fitted by a sequence of polytropic distributions with a time-dependent index n(t) both in the non-magnetized and magnetized regimes. These numerical results show that polytropic distributions (also called Tsallis distributions) provide in many cases a good fit of the QSSs. However, in order to moderate our message, we also report a case where the Lynden-Bell theory provides an excellent prediction of an inhomogeneous QSS.

preprint2011arXiv

Bose-Einstein Condensate general relativistic stars

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

preprint2011arXiv

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

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

preprint2011arXiv

Relaxation equations for two-dimensional turbulent flows with a prior vorticity distribution

Using a Maximum Entropy Production Principle (MEPP), we derive a new type of relaxation equations for two-dimensional turbulent flows in the case where a prior vorticity distribution is prescribed instead of the Casimir constraints [Ellis, Haven, Turkington, Nonlin., 15, 239 (2002)]. The particular case of a Gaussian prior is specifically treated in connection to minimum enstrophy states and Fofonoff flows. These relaxation equations are compared with other relaxation equations proposed by Robert and Sommeria [Phys. Rev. Lett. 69, 2776 (1992)] and Chavanis [Physica D, 237, 1998 (2008)]. They can provide a small-scale parametrization of 2D turbulence or serve as numerical algorithms to compute maximum entropy states with appropriate constraints. We perform numerical simulations of these relaxation equations in order to illustrate geometry induced phase transitions in geophysical flows.

preprint2011arXiv

Statistical mechanics of Fofonoff flows in an oceanic basin

We study the minimization of potential enstrophy at fixed circulation and energy in an oceanic basin with arbitrary topography. For illustration, we consider a rectangular basin and a linear topography h=by which represents either a real bottom topography or the beta-effect appropriate to oceanic situations. Our minimum enstrophy principle is motivated by different arguments of statistical mechanics reviewed in the article. It leads to steady states of the quasigeostrophic (QG) equations characterized by a linear relationship between potential vorticity q and stream function psi. For low values of the energy, we recover Fofonoff flows [J. Mar. Res. 13, 254 (1954)] that display a strong westward jet. For large values of the energy, we obtain geometry induced phase transitions between monopoles and dipoles similar to those found by Chavanis and Sommeria [J. Fluid Mech. 314, 267 (1996)] in the absence of topography. In the presence of topography, we recover and confirm the results obtained by Venaille and Bouchet [Phys. Rev. Lett. 102, 104501 (2009)] using a different formalism. In addition, we introduce relaxation equations towards minimum potential enstrophy states and perform numerical simulations to illustrate the phase transitions in a rectangular oceanic basin with linear topography (or beta-effect).

preprint2011arXiv

Statistical mechanics of two-dimensional Euler flows and minimum enstrophy states

A simplified thermodynamic approach of the incompressible 2D Euler equation is considered based on the conservation of energy, circulation and microscopic enstrophy. Statistical equilibrium states are obtained by maximizing the Miller-Robert-Sommeria (MRS) entropy under these sole constraints. The vorticity fluctuations are Gaussian while the mean flow is characterized by a linear $\barω-ψ$ relationship. Furthermore, the maximization of entropy at fixed energy, circulation and microscopic enstrophy is equivalent to the minimization of macroscopic enstrophy at fixed energy and circulation. This provides a justification of the minimum enstrophy principle from statistical mechanics when only the microscopic enstrophy is conserved among the infinite class of Casimir constraints. A new class of relaxation equations towards the statistical equilibrium state is derived. These equations can provide an effective description of the dynamics towards equilibrium or serve as numerical algorithms to determine maximum entropy or minimum enstrophy states. We use these relaxation equations to study geometry induced phase transitions in rectangular domains. In particular, we illustrate with the relaxation equations the transition between monopoles and dipoles predicted by Chavanis and Sommeria [J. Fluid. Mech. 314, 267 (1996)]. We take into account stable as well as metastable states and show that metastable states are robust and have negative specific heats. This is the first evidence of negative specific heats in that context. We also argue that saddle points of entropy can be long-lived and play a role in the dynamics because the system may not spontaneously generate the perturbations that destabilize them.

preprint2010arXiv

Can one measure negative specific heat in the canonical statistical ensemble?

According to thermodynamics, the specific heat of Boltzmannian short-range interacting systems is a positive quantity. Less intuitive properties are instead displayed by systems characterized by long-range interactions. In that case, the sign of specific heat depends on the considered statistical ensemble: negative specific heat can be found in isolated systems, which are studied in the framework of the microcanonical ensemble; on the other hand, it is generally recognized that a positive specific heat should always be measured in systems in contact with a thermal bath, for which the canonical ensemble is the appropriate one. We demonstrate that the latter assumption is not generally true: one can in principle measure negative specific heat also in the canonical ensemble if the system under scrutiny is non-Boltzmannian and/or out-of-equilibrium.

preprint2010arXiv

Out-of-equilibrium phase transitions in the HMF model: a closer look

We provide a detailed discussion of out-of-equilibrium phase transitions in the Hamiltonian Mean Field (HMF) model in the framework of Lynden-Bell's statistical theory of the Vlasov equation. For two-levels initial conditions, the caloric curve $β(E)$ only depends on the initial value $f_0$ of the distribution function. We evidence different regions in the parameter space where the nature of phase transitions between magnetized and non-magnetized states changes: (i) for $f_0>0.10965$, the system displays a second order phase transition; (ii) for $0.109497<f_0<0.10965$, the system displays a second order phase transition and a first order phase transition; (iii) for $0.10947<f_0<0.109497$, the system displays two second order phase transitions; (iv) for $f_0<0.10947$, there is no phase transition. The passage from a first order to a second order phase transition corresponds to a tricritical point. The sudden appearance of two second order phase transitions from nothing corresponds to a second order azeotropy. This is associated with a phenomenon of phase reentrance. When metastable states are taken into account, the problem becomes even richer. In particular, we find a new situation of phase reentrance. We consider both microcanonical and canonical ensembles and report the existence of a tiny region of ensembles inequivalence. We also explain why the use of the initial magnetization $M_0$ as an external parameter, instead of the phase level $f_0$, may lead to inconsistencies in the thermodynamical analysis.

preprint2010arXiv

Phase transitions in self-gravitating systems and bacterial populations with a screened attractive potential

We consider a system of particles interacting via a screened Newtonian potential and study phase transitions between homogeneous and inhomogeneous states in the microcanonical and canonical ensembles. Like for other systems with long-range interactions, we obtain a great diversity of microcanonical and canonical phase transitions depending on the dimension of space and on the importance of the screening length. We also consider a system of particles in Newtonian interaction in the presence of a ``neutralizing background''. By a proper interpretation of the parameters, our study describes (i) self-gravitating systems in a cosmological setting, and (ii) chemotaxis of bacterial populations in the original Keller-Segel model.

preprint2010arXiv

Self-gravitating Brownian particles in two dimensions: the case of N=2 particles

We study the motion of N=2 overdamped Brownian particles in gravitational interaction in a space of dimension d=2. This is equivalent to the simplified motion of two biological entities interacting via chemotaxis when time delay and degradation of the chemical are ignored. This problem also bears some similarities with the stochastic motion of two point vortices in viscous hydrodynamics [Agullo & Verga, Phys. Rev. E, 63, 056304 (2001)]. We analytically obtain the density probability of finding the particles at a distance r from each other at time t. We also determine the probability that the particles have coalesced and formed a Dirac peak at time t (i.e. the probability that the reduced particle has reached r=0 at time t). Finally, we investigate the variance of the distribution <r^2> and discuss the proper form of the virial theorem for this system. The reduced particle has a normal diffusion behaviour for small times with a gravity-modified diffusion coefficient <r^2>=r_0^2+(4k_B/ξμ)(T-T_*)t, where k_BT_{*}=Gm_1m_2/2 is a critical temperature, and an anomalous diffusion for large times <r^2>~t^(1-T_*/T). As a by-product, our solution also describes the growth of the Dirac peak (condensate) that forms in the post-collapse regime of the Smoluchowski-Poisson system (or Keller-Segel model) for T<T_c=GMm/(4k_B). We find that the saturation of the mass of the condensate to the total mass is algebraic in an infinite domain and exponential in a bounded domain.

preprint2009arXiv

Out-of-equilibrium phase re-entrance(s) in long-range interacting systems

Systems with long-range interactions display a short-time relaxation towards Quasi Stationary States (QSSs) whose lifetime increases with system size. The application of Lynden-Bell's theory of "violent relaxation" to the Hamiltonian Mean Field model leads to the prediction of out-of-equilibrium first and second order phase transitions between homogeneous (zero magnetization) and inhomogeneous (non-zero magnetization) QSSs, as well as an interesting phenomenon of phase re-entrances. We compare these theoretical predictions with direct $N$-body numerical simulations. We confirm the existence of phase re-entrance in the typical parameter range predicted from Lynden-Bell's theory, but also show that the picture is more complicated than initially thought. In particular, we exhibit the existence of secondary re-entrant phases: we find un-magnetized states in the theoretically magnetized region as well as persisting magnetized states in the theoretically unmagnetized region.

preprint2008arXiv

Dynamical stability of systems with long-range interactions: application of the Nyquist method to the HMF model

We apply the Nyquist method to the Hamiltonian Mean Field (HMF) model in order to settle the linear dynamical stability of a spatially homogeneous distribution function with respect to the Vlasov equation. We consider the case of Maxwell (isothermal) and Tsallis (polytropic) distributions and show that the system is stable above a critical kinetic temperature T_c and unstable below it. Then, we consider a symmetric double-humped distribution, made of the superposition of two decentered Maxwellians, and show the existence of a re-entrant phase in the stability diagram. When we consider an asymmetric double-humped distribution, the re-entrant phase disappears above a critical value of the asymmetry factor Delta>1.09. We also consider the HMF model with a repulsive interaction. In that case, single-humped distributions are always stable. For asymmetric double-humped distributions, there is a re-entrant phase for 1<Delta<25.6, a double re-entrant phase for 25.6<Delta<43.9 and no re-entrant phase for Delta>43.9. Finally, we extend our results to arbitrary potentials of interaction and mention the connexion between the HMF model, Coulombian plasmas and gravitational systems. We discuss the relation between linear dynamical stability and formal nonlinear dynamical stability and show their equivalence for spatially homogeneous distributions. We also provide a criterion of dynamical stability for inhomogeneous systems.