اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدVlasov-Poisson in 1D for initially cold systems: post-collapse Lagrangian perturbation theory
148
0
0.0
(
0
)
تأليف
Stephane Colombi
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We study analytically the collapse of an initially smooth, cold, self-gravitating collisionless system in one dimension. The system is described as a central S shape in phase-space surrounded by a nearly stationary halo acting locally like a harmonic background on the S. To resolve the dynamics of the S under its self-gravity and under the influence of the halo, we introduce a novel approach using post-collapse Lagrangian perturbation theory. This approach allows us to follow the evolution of the system between successive crossing times and to describe in an iterative way the interplay between the central S and the halo. Our theoretical predictions are checked against measurements in entropy conserving numerical simulations based on the waterbag method. While our post-collapse Lagrangian approach does not allow us to compute rigorously the long term behavior of the system, i.e. after many crossing times, it explains the close to power-law behavior of the projected density observed in numerical simulations. Pushing the model at late time suggests that the system could build at some point a very small flat core, but this is very speculative. This analysis shows that understanding the dynamics of initially cold systems requires a fine grained approach for a correct description of their very central part. The analyses performed here can certainly be extended to spherical symmetry.
قيم البحث
اقرأ أيضاً
We develop a new perturbation theory (PT) treatment that can describe gravitational dynamics of large-scale structure after shell-crossing in the one-dimensional cosmological case. Starting with cold initial conditions, the motion of matter distribut
ion follows at early stages the single-stream regime, which can, in one dimension, be described exactly by the first-order Lagrangian perturbation, i.e. the Zeldovich solution. However, the single-stream flow no longer holds after shell-crossing and a proper account of the multi-stream flow is essential for post-collapse dynamics. In this paper, extending previous work by Colombi (2015, MNRAS 446, 2902), we present a perturbative description for the multi-stream flow after shell-crossing in a cosmological setup. In addition, we introduce an adaptive smoothing scheme to deal with the bulk properties of phase-space structures. The filtering scales in this scheme are linked to the next-crossing time in the post-collapse region, estimated from our PT calculations. Our PT treatment combined with adaptive smoothing is illustrated in several cases. Predictions are compared to simulations and we find that post-collapse PT with adaptive smoothing reproduces the power spectrum and phase-space structures remarkably well even at small scales, where Zeldovich solution substantially deviates from simulations.
We propose a new semi-Lagrangian Vlasov-Poisson solver. It employs elements of metric to follow locally the flow and its deformation, allowing one to find quickly and accurately the initial phase-space position $Q(P)$ of any test particle $P$, by exp
anding at second order the geometry of the motion in the vicinity of the closest element. It is thus possible to reconstruct accurately the phase-space distribution function at any time $t$ and position $P$ by proper interpolation of initial conditions, following Liouville theorem. When distorsion of the elements of metric becomes too large, it is necessary to create new initial conditions along with isotropic elements and repeat the procedure again until next resampling. To speed up the process, interpolation of the phase-space distribution is performed at second order during the transport phase, while third order splines are used at the moments of remapping. We also show how to compute accurately the region of influence of each element of metric with the proper percolation scheme. The algorithm is tested here in the framework of one-dimensional gravitational dynamics but is implemented in such a way that it can be extended easily to four or six-dimensional phase-space. It can also be trivially generalised to plasmas.
We revisit in one dimension the waterbag method to solve numerically Vlasov-Poisson equations. In this approach, the phase-space distribution function $f(x,v)$ is initially sampled by an ensemble of patches, the waterbags, where $f$ is assumed to be
constant. As a consequence of Liouville theorem it is only needed to follow the evolution of the border of these waterbags, which can be done by employing an orientated, self-adaptive polygon tracing isocontours of $f$. This method, which is entropy conserving in essence, is very accurate and can trace very well non linear instabilities as illustrated by specific examples. As an application of the method, we generate an ensemble of single waterbag simulations with decreasing thickness, to perform a convergence study to the cold case. Our measurements show that the system relaxes to a steady state where the gravitational potential profile is a power-law of slowly varying index $beta$, with $beta$ close to $3/2$ as found in the literature. However, detailed analysis of the properties of the gravitational potential shows that at the center, $beta > 1.54$. Moreover, our measurements are consistent with the value $beta=8/5=1.6$ that can be analytically derived by assuming that the average of the phase-space density per energy level obtained at crossing times is conserved during the mixing phase. These results are incompatible with the logarithmic slope of the projected density profile $beta-2 simeq -0.47$ obtained recently by Schulz et al. (2013) using a $N$-body technique. This sheds again strong doubts on the capability of $N$-body techniques to converge to the correct steady state expected in the continuous limit.
We present the one-loop 2-point function of biased tracers in redshift space computed with Lagrangian perturbation theory, including a full resummation of both long-wavelength (infrared) displacements and associated velocities. The resulting model ac
curately predicts the power spectrum and correlation function of halos and mock galaxies from two different sets of N-body simulations at the percent level for quasi-linear scales, including the damping of the baryon acoustic oscillation signal due to the bulk motions of galaxies. We compare this full resummation with other, approximate, techniques including the moment expansion and Gaussian streaming model. We discuss infrared resummation in detail and compare our Lagrangian formulation with the Eulerian theory augmented by an infrared resummation based on splitting the input power spectrum into wiggle and no-wiggle components. We show that our model is able to recover unbiased cosmological parameters in mock data encompassing a volume much larger than what will be available to future galaxy surveys. We demonstrate how to efficiently compute the resulting expressions numerically, making available a fast Python code capable of rapidly computing these statistics in both configuration and Fourier space.
Self-gravitating quantum matter may exist in a wide range of cosmological and astrophysical settings from the very early universe through to present-day boson stars. Such quantum matter arises in a number of different theories, including the Peccei-Q
uinn axion and UltraLight (ULDM) or Fuzzy (FDM) dark matter scenarios. We consider the dynamical evolution of perturbations to the spherically symmetric soliton, the ground state solution to the Schr{o}dinger-Poisson system common to all these scenarios. We construct the eigenstates of the Schr{o}dinger equation, holding the gravitational potential fixed to its ground state value. We see that the eigenstates qualitatively capture the properties seen in full ULDM simulations, including the soliton breathing mode, the random walk of the soliton center, and quadrupolar distortions of the soliton. We then show that the time-evolution of the gravitational potential and its impact on the perturbations can be well described within the framework of time-dependent perturbation theory. Applying our formalism to a synthetic ULDM halo reveals considerable mixing of eigenstates, even though the overall density profile is relatively stable. Our results provide a new analytic approach to understanding the evolution of these systems as well as possibilities for faster approximate simulations.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
