Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userLagrangian cosmological perturbation theory at shell-crossing
74
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We consider the growth of primordial dark matter halos seeded by three crossed initial sine waves of various amplitudes. Using a Lagrangian treatment of cosmological gravitational dynamics, we examine the convergence properties of a high-order perturbative expansion in the vicinity of shell-crossing, by comparing the analytical results with state-of-the-art high resolution Vlasov-Poisson simulations. Based on a quantitative exploration of parameter space, we study explicitly for the first time the convergence speed of the perturbative series, and find, in agreement with intuition, that it slows down when going from quasi one-dimensional initial conditions (one sine wave dominating) to quasi triaxial symmetry (three sine waves with same amplitude). In most cases, the system structure at collapse time is, as expected, very similar to what is obtained with simple one-dimensional dynamics, except in the quasi-triaxial regime, where the phase-space sheet presents a velocity spike. In all cases, the perturbative series exhibits a generic convergence behavior as fast as an exponential of a power-law of the order of the expansion, allowing one to numerically extrapolate it to infinite order. The results of such an extrapolation agree remarkably well with the simulations, even at shell-crossing.
rate research
Read More
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 distribution 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 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 accurately 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.
Cosmological perturbation theory is crucial for our understanding of the universe. The linear theory has been well understood for some time, however developing and applying the theory beyond linear order is currently at the forefront of research in theoretical cosmology. This thesis studies the applications of perturbation theory to cosmology and, specifically, to the early universe. Starting with some background material introducing the well-tested standard model of cosmology, we move on to develop the formalism for perturbation theory up to second order giving evolution equations for all types of scalar, vector and tensor perturbations, both in gauge dependent and gauge invariant form. We then move on to the main result of the thesis, showing that, at second order in perturbation theory, vorticity is sourced by a coupling term quadratic in energy density and entropy perturbations. This source term implies a qualitative difference to linear order. Thus, while at linear order vorticity decays with the expansion of the universe, the same is not true at higher orders. This will have important implications on future measurements of the polarisation of the Cosmic Microwave Background, and could give rise to the generation of a primordial seed magnetic field. Having derived this qualitative result, we then estimate the scale dependence and magnitude of the vorticity power spectrum, finding, for simple power law inputs a small, blue spectrum. The final part of this thesis concerns higher order perturbation theory, deriving, for the first time, the metric tensor, gauge transformation rules and governing equations for fully general third order perturbations. We close with a discussion of natural extensions to this work and other possible ideas for off-shooting projects in this continually growing field.
We compare and contrast two different metric based formulations of non- linear cosmological perturbation theory: the MW2009 approach in [K. A. Malik and D. Wands, Phys. Rept. 475 (2009), 1.] following Bardeen and the recent approach of the paper KN2010 [K. Nakamura, Advances in Astronomy 2010 (2010), 576273]. We present each formulation separately. In the MW2009 approach, one considers the gauge transformations of perturbative quantities, choosing a gauge by requiring that certain quantities vanish, rendering all other variables gauge invariant. In the KN2010 formalism, one decomposes the metric tensor into a gauge variant and gauge invariant part from the outset. We compare the two approaches in both the longitudinal and uniform curvature gauges. In the longitudinal gauge, we find that Nakamuras gauge invariant variables correspond exactly to those in the longitudinal gauge (i.e., for scalar perturbations, to the Bardeen potentials), and in the uniform curvature gauge we obtain the usual relationship between gauge invariant variables in the flat and longitudinal gauge. Thus, we show that these two approaches are equivalent.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
