No Arabic abstract
In the mean field limit, isolated gravitational systems often evolve towards a steady state through a violent relaxation phase. One question is to understand the nature of this relaxation phase, in particular the role of radial instabilities in the establishment/destruction of the steady profile. Here, through a detailed phase-space analysis based both on a spherical Vlasov solver, a shell code and a $N$-body code, we revisit the evolution of collisionless self-gravitating spherical systems with initial power-law density profiles $rho(r) propto r^n$, $0 leq n leq -1.5$, and Gaussian velocity dispersion. Two sub-classes of models are considered, with initial virial ratios $eta=0.5$ (warm) and $eta=0.1$ (cool). Thanks to the numerical techniques used and the high resolution of the simulations, our numerical analyses are able, for the first time, to show the clear separation between two or three well known dynamical phases: (i) the establishment of a spherical quasi-steady state through a violent relaxation phase during which the phase-space density displays a smooth spiral structure presenting a morphology consistent with predictions from self-similar dynamics, (ii) a quasi-steady state phase during which radial instabilities can take place at small scales and destroy the spiral structure but do not change quantitatively the properties of the phase-space distribution at the coarse grained level and (iii) relaxation to non spherical state due to radial orbit instabilities for $n leq -1$ in the cool case.
The long timescale evolution of a self-gravitating system is generically driven by two-body encounters. In many cases, the motion of the particles is primarily governed by the mean field potential. When this potential is integrable, particles move on nearly fixed orbits, which can be described in terms of angle-action variables. The mean field potential drives fast orbital motions (angles) whose associated orbits (actions) are adiabatically conserved on short dynamical timescales. The long-term stochastic evolution of the actions is driven by the potential fluctuations around the mean field and in particular by resonant two-body encounters, for which the angular frequencies of two particles are in resonance. We show that the stochastic gravitational fluctuations acting on the particles can generically be described by a correlated Gaussian noise. Using this approach, the so-called $eta$-formalism, we derive a diffusion equation for the actions in the test particle limit. We show that in the appropriate limits, this diffusion equation is equivalent to the inhomogeneous Balescu-Lenard and Landau equations. This approach provides a new view of the resonant diffusion processes associated with long-term orbital distortions. Finally, by investigating the example of the Hamiltonian Mean Field Model, we show how the present method generically allows for alternative calculations of the long-term diffusion coefficients in inhomogeneous systems.
We conduct numerical experiments to determine the density probability distribution function (PDF) produced in supersonic, isothermal, self-gravitating turbulence of the sort that is ubiquitous in star-forming molecular clouds. Our experiments cover a wide range of turbulent Mach number and virial parameter, allowing us for the first time to determine how the PDF responds as these parameters vary, and we introduce a new diagnostic, the dimensionless star formation efficiency versus density ($epsilon_{rm ff}(s)$) curve, which provides a sensitive diagnostic of the PDF shape and dynamics. We show that the PDF follows a universal functional form consisting of a log-normal at low density with two distinct power law tails at higher density; the first of these represents the onset of self-gravitation, and the second reflects the onset of rotational support. Once the star formation efficiency reaches a few percent, the PDF becomes statistically steady, with no evidence for secular time-evolution at star formation efficiencies from about five to 20 percent. We show that both the Mach number and the virial parameter influence the characteristic densities at which the log-normal gives way to the first power-law, and the first to the second, and we extend (for the former) and develop (for the latter) simple theoretical models for the relationship between these density thresholds and the global properties of the turbulent medium.
Dynamics and collapse of collisionless self-gravitating systems is described by the coupled collisionless Boltzmann and Poisson equations derived from $f(R)$-gravity in the weak field approximation. Specifically, we describe a system at equilibrium by a time-independent distribution function $f_0(x,v)$ and two potentials $Phi_0(x)$ and $Psi_0(x)$ solutions of the modified Poisson and collisionless Boltzmann equations. Considering a small perturbation from the equilibrium and linearizing the field equations, it can be obtained a dispersion relation. A dispersion equation is achieved for neutral dust-particle systems where a generalized Jeans wave-number is obtained. This analysis gives rise to unstable modes not present in the standard Jeans analysis (derived assuming Newtonian gravity as weak filed limit of $f(R)=R$). In this perspective, we discuss several self-gravitating astrophysical systems whose dynamics could be fully addressed in the framework of $f(R)$-gravity.
The short--wave asymptotics (WKB) of spiral density waves in self-gravitating stellar discs is well suited for the study of the dynamics of tightly--wound wavepackets. But the textbook WKB theory is not well adapted to the study of the linear eigenmodes in a collisionless self-gravitating disc because of the transcendental nature of the dispersion relation. We present a modified WKB of spiral density waves, for collisionless discs in the epicyclic limit, in which the perturbed gravitational potential is related to the perturbed surface density by the Poisson integral in Kalnajs logarithmic spiral form. An integral equation is obtained for the surface density perturbation, which is seen to also reduce to the standard WKB dispersion relation. We specialize to a low mass (or Keplerian) self-gravitating disc around a massive black hole, and derive an integral equation governing the eigenspectra and eigenfunctions of slow precessional modes. For a prograde disc, the integral kernel turns out be real and symmetric, implying that all slow modes are stable. We apply the slow mode integral equation to two unperturbed disc profiles, the Jalali--Tremaine annular discs, and the Kuzmin disc. We determine eigenvalues and eigenfunctions for both $m = 1$ and $m = 2$ slow modes for these profiles and discuss their properties. Our results compare well with those of Jalali--Tremaine.
We study the statistical mechanics of binary systems under gravitational interaction of the Modified Newtonian Dynamics (MOND) in three-dimensional space. Considering the binary systems, in the microcanonical and canonical ensembles, we show that in the microcanonical systems, unlike the Newtonian gravity, there is a sharp phase transition, with a high-temperature homogeneous phase and a low temperature clumped binary one. Defining an order parameter in the canonical systems, we find a smoother phase transition and identify the corresponding critical temperature in terms of physical parameters of the binary system.