Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRelaxation in self-gravitating systems
102
0
0.0
(
0
)
Added by
Jean-Baptiste Fouvry
Publication date
2017
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
A new discrete model for energy relaxation of a quantum particle is described via a projection operator, causing the wave function collapse. Power laws for the evolution of the particle coordinate and momentum dispersions are derived. A new dissipative Schrodinger equation is proposed and solved for particular cases. A new dissipative Liouville equation is heuristically constructed.
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.
We show, using the N-body code GADGET-2, that stellar scattering by massive clumps can produce exponential discs, and the effectiveness of the process depends on the mass of scattering centres, as well as the stability of the galactic disc. Heavy, dense scattering centres in a less stable disc generate an exponential profile quickly, with a timescale shorter than 1 Gyr. The profile evolution due to scattering can make a near-exponential disc under various initial stellar distributions. This result supports analytic theories that predict the scattering processes always favour the zero entropy gradient solution to the Jeans/Poisson equations, whose profile is a near-exponential. Profile changes are accompanied by disc thickening, and a power-law increase in stellar velocity dispersion in both vertical and radial directions is also observed through the evolution. Close encounters between stars and clumps can produce abrupt changes in stellar orbits and shift stars radially. These events can make trajectories more eccentric, but many leave eccentricities little changed. On average, orbital eccentricities of stars increase moderately with time.
When an open system of classical point particles interacting by Newtonian gravity collapses and relaxes violently, an arbitrary amount of energy may in principle be carried away by particles which escape to infinity. We investigate here, using numerical simulations, how this released energy and other related quantities (notably the binding energy and size of the virialized structure) depends on the initial conditions, for the one parameter family of starting configurations given by randomly distributing N cold particles in a spherical volume. Previous studies have established that the minimal size reached by the system scales approximately as N^{-1/3}, a behaviour which follows trivially when the growth of perturbations (which regularize the singularity of the cold collapse in the infinite N limit) are assumed to be unaffected by the boundaries. Our study shows that the energy ejected grows approximately in proportion to N^{1/3}, while the fraction of the initial mass ejected grows only very slowly with N, approximately logarithmically, in the range of N simulated. We examine in detail the mechanism of this mass and energy ejection, showing explicitly that it arises from the interplay of the growth of perturbations with the finite size of the system. A net lag of particles compared to their uniform spherical collapse trajectories develops first at the boundaries and then propagates into the volume during the collapse. Particles in the outer shells are then ejected as they scatter through the time dependent potential of an already re-expanding central core. Using modified initial configurations we explore the importance of fluctuations at different scales, and discreteness (i.e. non-Vlasov) effects in the dynamics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
