Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA new framework for numerical simulations of structure formation
113
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The diversity of structures in the Universe (from the smallest galaxies to the largest superclusters) has formed under the pull of gravity from the tiny primordial perturbations that we see imprinted in the cosmic microwave background. A quantitative description of this process would require description of motion of zillions of dark matter particles. This impossible task is usually circumvented by coarse-graining the problem: one either considers a Newtonian dynamics of particles with macroscopically large masses or approximates the dark matter distribution with a continuous density field. There is no closed system of equations for the evolution of the matter density field alone and instead it should still be discretized at each timestep. In this work we describe a method of solving the full 6-dimensional Vlasov-Poisson equation via a system of auxiliary Schroedinger-like equations. The complexity of the problem gets shifted into the choice of the number and shape of the initial wavefunctions that should only be specified at the beginning of the computation (we stress that these wavefunctions have nothing to do with quantum nature of the actual dark matter particles). We discuss different prescriptions to generate the initial wave functions from the initial conditions and demonstrate the validity of the technique on two simple test cases. This new simulation algorithm can in principle be used on an arbitrary distribution function, enabling the simulation of warm and hot dark matter structure formation scenarios.
rate research
Read More
Dust emission is the main foreground for cosmic microwave background (CMB) polarization. Its statistical characterization must be derived from the analysis of observational data because the precision required for a reliable component separation is far greater than what is currently achievable with physical models of the turbulent magnetized interstellar medium. This letter takes a significant step toward this goal by proposing a method that retrieves non-Gaussian statistical characteristics of dust emission from noisy Planck polarization observations at 353 GHz. We devised a statistical denoising method based on wavelet phase harmonics (WPH) statistics, which characterize the coherent structures in non-Gaussian random fields and define a generative model of the data. The method was validated on mock data combining a dust map from a magnetohydrodynamic simulation and Planck noise maps. The denoised map reproduces the true power spectrum down to scales where the noise power is an order of magnitude larger than that of the signal. It remains highly correlated to the true emission and retrieves some of its non-Gaussian properties. Applied to Planck data, the method provides a new approach to building a generative model of dust polarization that will characterize the full complexity of the dust emission. We also release PyWPH, a public Python package, to perform GPU-accelerated WPH analyses on images.
We employ cosmological hydrodynamical simulations to study the growth of massive black holes (BHs) at high redshifts subject to BH merger recoils from gravitational wave emission. We select the most massive dark matter halo at z=6 from the Millennium simulation, and resimulate its formation at much higher resolution including gas physics and a model for BH seeding, growth and feedback. Assuming that the initial BH seeds are relatively massive, of the order of 10^5 Msun, and that seeding occurs around z~15 in dark matter haloes of mass 10^9-10^10 Msun, we find that it is possible to build up supermassive BHs (SMBHs) by z=6 that assemble most of their mass during extended Eddington-limited accretion periods. The properties of the simulated SMBHs are consistent with observations of z=6 quasars in terms of the estimated BH masses and bolometric luminosities, the amount of star formation occurring within the host halo, and the presence of highly enriched gas in the innermost regions of the host galaxy. After a peak in the BH accretion rate at z=6, the most massive BH has become sufficiently massive for the growth to enter into a much slower phase of feedback-regulated accretion. We explore the full range of expected recoils and radiative efficiencies, and also consider models with spinning BHs. In the most `pessimistic case where BH spins are initially high, we find that the growth of the SMBHs can be potentially hampered if they grow mostly in isolation and experience only a small number of mergers. Whereas BH kicks can expel a substantial fraction of low mass BHs, they do not significantly affect the build up of the SMBHs. On the contrary, a large number of BH mergers has beneficial consequences for the growth of the SMBHs by considerably reducing their spin. [Abridged]
We review the numerical techniques for ideal and non-ideal magneto-hydrodynamics (MHD) used in the context of star formation simulations. We outline the specific challenges offered by modeling star forming environments, which are dominated by supersonic and super-Alfvenic turbulence in a radiative, self-gravitating fluid. These conditions are rather unique in physics and engineering and pose particularly severe restrictions on the robustness and accuracy of numerical codes. One striking aspect is the formation of collapsing fluid elements leading to the formation of singularities that represent point-like objects, namely the proto-stars. Although a few studies have attempted to resolve the formation of the first and second Larson cores, resolution limitations force us to use sink particle techniques, with sub-grid models to compute the accretion rates of mass, momentum and energy, as well as their ejection rate due to radiation and jets from the proto-stars. We discuss the most popular discretisation techniques used in the community, namely smoothed particle hydrodynamics, finite difference and finite volume methods, stressing the importance to maintain a divergence-free magnetic field. We discuss how to estimate the truncation error of a given numerical scheme, and its importance in setting the magnitude of the numerical diffusion. This can have a strong impact on the outcome of these MHD simulations, where both viscosity and resistivity are implemented at the grid scale. We then present various numerical techniques to model non-ideal MHD effects, such as Ohmic and ambipolar diffusion, as well as the Hall effect. These important physical ingredients are posing strong challenges in term of resolution and time stepping. For the latter, several strategies are discussed to overcome the limitations due to prohibitively small time steps (abridged).
We present the curation and verification of a new combined optical and near infrared dataset for cosmology and astrophysics, derived from the combination of $ugri$-band imaging from the Kilo Degree Survey (KiDS) and $ZY!J!H!K_{rm s}$-band imaging from the VISTA Kilo degree Infrared Galaxy (VIKING) survey. This dataset is unrivaled in cosmological imaging surveys due to its combination of area ($458$ deg$^2$ before masking), depth ($rle25$), and wavelength coverage ($ugriZY!J!H!K_{rm s}$). The combination of survey depth, area, and (most importantly) wavelength coverage allows significant reductions in systematic uncertainties (i.e. reductions of between 10 and 60% in bias, outlier rate, and scatter) in photometric-to-spectroscopic redshift comparisons, compared to the optical-only case at photo-$z$ above $0.7$. The complementarity between our optical and NIR surveys means that over $80%$ of our sources, across all photo-$z$, have significant detections (i.e. not upper limits) in our $8$ reddest bands. We derive photometry, photo-$z$, and stellar masses for all sources in the survey, and verify these data products against existing spectroscopic galaxy samples. We demonstrate the fidelity of our higher-level data products by constructing the survey stellar mass functions in 8 volume-complete redshift bins. We find that these photometrically derived mass functions provide excellent agreement with previous mass evolution studies derived using spectroscopic surveys. The primary data products presented in this paper are publicly available at http://kids.strw.leidenuniv.nl/.
We present a novel population-based Bayesian inference approach to model the average and population variance of spatial distribution of a set of observables from ensemble analysis of low signal-to-noise ratio measurements. The method consists of (1) inferring the average profile using Gaussian Processes and (2) computing the covariance of the profile observables given a set of independent variables. Our model is computationally efficient and capable of inferring average profiles of a large population size from noisy measurements, without stacking and binning data nor parameterizing the shape of the mean profile. We demonstrate the performance of our method using dark matter, gas and stellar profiles extracted from hydrodynamical cosmological simulations of galaxy formation. Population Profile Estimator (PoPE) is publicly available in a GitHub repository. Our new method should be useful for measuring the spatial distribution and internal structure of a variety of astrophysical systems using large astronomical surveys.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
