Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExploring the non-linear density field in the Millennium simulations with tessellations - I. The probability distribution function
102
0
0.0
(
0
)
Added by
Biswajit Pandey B.P.
Publication date
2013
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We use the Delaunay Tessellation Field Estimator (DTFE) to study the one-point density distribution functions of the Millennium (MS) and Millennium-II (MS-II) simulations. The DTFE technique is based directly on the particle positions, without requiring any type of smoothing or analysis grid, thereby providing high sensitivity to all non-linear structures resolved by the simulations. In order to identify the detailed origin of the shape of the one-point density probability distribution function (PDF), we decompose the simulation particles according to the mass of their host FoF halos, and examine the contributions of different halo mass ranges to the global density PDF. We model the one-point distribution of the FoF halos in each halo mass bin with a set of Monte Carlo realizations of idealized NFW dark matter halos, finding that this reproduces the measurements from the N-body simulations reasonably well, except for a small excess present in simulation results. This excess increases with increasing halo mass. We show that its origin lies in substructure, which becomes progressively more abundant and better resolved in more massive dark matter halos. We demonstrate that the high density tail of the one-point distribution function in less massive halos is severely affected by the gravitational softening length and the mass resolution. In particular, we find these two parameters to be more important for an accurate measurement of the density PDF than the simulated volume. Combining our results from individual halo mass bins we find that the part of the one-point density PDF originating from collapsed halos can nevertheless be quite well described by a simple superposition of a set of NFW halos with the expected cosmological abundance over the resolved mass range. The transition region to the low-density unbound material is however not well captured by such an analytic halo model.
rate research
Read More
We introduce the position-dependent probability distribution function (PDF) of the smoothed matter field as a cosmological observable. In comparison to the PDF itself, the spatial variation of the position-dependent PDF is simpler to model and has distinct dependence on cosmological parameters. We demonstrate that the position-dependent PDF is characterized by variations in the local mean density, and we compute the linear response of the PDF to the local density using separate universe N-body simulations. The linear response of the PDF to the local density field can be thought of as the linear bias of regions of the matter field selected based on density. We provide a model for the linear response, which accurately predicts our simulation measurements. We also validate our results and test the separate universe consistency relation for the local PDF using global universe simulations. We find excellent agreement between the two, and we demonstrate that the separate universe method gives a lower variance determination of the linear response.
In the context of tomographic cosmic shear surveys, a theoretical model for the one-point statistics of the aperture mass (Map) is developed. This formalism is based on the application of the large deviation principle to the projected matter density field and more specifically to the angular aperture masses. The latter holds the advantage of being an observable that can be directly extracted from the observed shear field and to be, by construction, independent from the long wave modes. Furthermore we show that, with the help of a nulling procedure based on the so-called BNT transform, it is possible to build observables that depend only on a finite range of redshifts making them also independent from the small-scale modes. This procedure makes predictions for the shape of the one-point Probability Distribution Function of such an observable very accurate, comparable to what had been previously obtained for 3D observables. Comparisons with specific simulations reveal however inconsistent results showing that synthetic lensing maps were not accurate enough for such refined observables. It points to the need for more precise dedicated numerical developments whose performances could be benchmarked with such observables. We furthermore review the possible systematics that could affect such a formalism in future weak-lensing surveys like Euclid, notably the impact of shape noise as well as leading corrections coming from lens-lens couplings, geodesic deviation, reduced shear and magnification bias.
The probability distribution functions (PDFs) for atomic, molecular, and total gas surface densities of M33 are determined at a resolution of about 50~pc over regions that share coherent morphological properties to unveil fingerprints of self-gravity across the star-forming disk. Most of the total gas PDFs from the central region to the edge of the star-forming disk are well-fitted by log-normal functions whose width decreases radially outwards. Because the HI velocity dispersion is approximately constant across the disk, the decrease of the PDF width is consistent with a lower Mach number for the turbulent ISM at large galactocentric radii where a higher fraction of HI is in the warm phase. The atomic gas is found mostly at face-on column densities below N$_{H}^{lim}$=2.5 10$^{21}$~cm$^{-2}$, with small radial variations of N$_{H}^{lim}$. The molecular gas PDFs do not show strong deviations from log-normal functions in the central region where molecular fractions are high. Here the high pressure and rate of star formation shapes the PDF as a log-normal function dispersing self-gravitating complexes with intense feedback at all column densities that are spatially resolved. Power law PDFs for the molecules are found near and above N$_H^{lim}$, in the well defined southern spiral arm and in a continuous dense filament extending at larger galactocentric radii; this is evident in cloud samples at different evolutionary stages along the star formation cycle. In the filament nearly half of the molecular gas departs from a log-normal PDF and power laws are also observed in pre-star forming molecular complexes. The slope of the power law is between -1 and -2. This slope, combined with maps showing where the different parts of the power law PDFs come from, suggest a power-law stratification of density within molecular cloud complexes, which is consistent with the dominance of self-gravity.
Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $zsim10$ galaxies. We characterize such lines of sight with the etendue $sigma_mu$, the area in the source plane magnified over some threshold $mu$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $sigma_mu$ beams on the sky, their properties, and efficient selection criteria. We define the best beams as having $sigma_{mu>3} >2000$ arcsec$^2$, for a $zsim10$ source plane, and predict $477 pm 21$ such beams on the sky. The total mass in the beam and $sigma_{mu>3}$ are strongly correlated. After controlling for total mass, we find a significant residual correlation between $sigma_{mu>3}$ and the number of cluster-scale halos ($>10^{14} M_odot h^{-1}$) in the beam. Beams with $sigma_{mu>3} >2000$ arcsec$^2$, which should be best at lensing $zsim10$ galaxies, are ten times more likely to contain multiple cluster-scale halos than a single cluster-scale halo. Beams containing an Abell 1689-like massive cluster halo often have additional structures along the line of sight, including at least one additional cluster-scale ($M_{200}>10^{14}M_odot h^{-1}$) halo 28% of the time. Selecting beams with multiple, massive structures will lead to enhanced detection of the most distant and intrinsically faint galaxies.
We use the spherical evolution approximation to investigate nonlinear evolution from the non-Gaussian initial conditions characteristic of the local f_nl model. We provide an analytic formula for the nonlinearly evolved probability distribution function of the dark matter which shows that the under-dense tail of the nonlinear PDF in the f_nl model should differ significantly from that for Gaussian initial conditions. Measurements of the under-dense tail in numerical simulations may be affected by discreteness effects, and we use a Poisson counting model to describe this effect. Once this has been accounted for, our model is in good quantitative agreement with the simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
