Abstract
We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum depending on the six cosmological parameters ...ωb, ωm, ns, h, $w$0, and σ8. It is constructed using the uncertainty quantification software UQLab using a spectral decomposition method called polynomial chaos expansion. All steps in its construction have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a simulation from the Euclid Flagship campaign and demonstrated to have converged up to wavenumbers $k\approx 5\, h\, {\rm Mpc}^{-1}$ for redshifts $z$ ≤ 5. The emulator is based on 100 input cosmologies simulated in boxes of (1250 Mpc/h)3 using 20483 particles. We show that by creating mock emulators it is possible to successfully predict and optimize the performance of the final emulator prior to performing any N-body simulations. The absolute accuracy of the final nonlinear power spectrum is as good as one obtained with N-body simulations, conservatively, ${\sim } 1$ per cent for $k\lesssim 1\, h\, {\rm Mpc}^{-1}$ and $z$ ≲ 1. This enables efficient forward modelling in the nonlinear regime, allowing for estimation of cosmological parameters using Markov Chain Monte Carlo methods. EuclidEmulator has been compared to HALOFIT, CosmicEmu, and NGenHalofit, and shown to be more accurate than these other approaches. This work paves a new way for optimal construction of future emulators that also consider other cosmological observables, use higher resolution input simulations, and investigate higher dimensional cosmological parameter spaces.
ABSTRACT We present initial results from the "Ponos" zoom-in numerical simulations of dark matter substructures in massive ellipticals. Two very highly resolved dark matter halos with Mvir = 1.2 × ...1013 and Mvir = 6.5 × 1012 and different ("violent" versus "quiescent") assembly histories have been simulated down to z = 0 in a ΛCDM cosmology with a total of 921,651,914 and 408,377,544 particles, respectively. Within the virial radius, the total mass fraction in self-bound Msub > 106 subhalos at the present epoch is 15% for the violent host and 16.5% for the quiescent one. At z = 0.7, these fractions increase to 19% and 33%, respectively, as more recently accreted satellites are less prone to tidal destruction. In projection, the average fraction of surface mass density in substructure at a distance of R/Rvir = 0.02 (∼5-10 kpc) from the two halo centers ranges from 0.6% to 2%, significantly higher than that measured in simulations of Milky Way-sized halos. The contribution of subhalos with Msub < 109 to the projected mass fraction is between one-fifth and one-third of the total, with the smallest share found in the quiescent host. We assess the impact of baryonic effects via twin, lower-resolution hydrodynamical simulations that include metallicity-dependent gas cooling, star formation, and a delayed-radiative-cooling scheme for supernova feedback. Baryonic contraction produces a super-isothermal total density profile and increases the number of massive subhalos in the inner regions of the main host. The host density profiles and projected subhalo mass fractions appear to be broadly consistent with observations of gravitational lenses.
We use a series of cosmological N-body simulations for a flat Λ cold dark matter (ΛCDM) cosmology to investigate the structural properties of dark matter haloes, at redshift zero, in the mass range 3 ...× 109 h−1 ≲Mvir ≲ 3 × 1013 h−1 M⊙. These properties include the concentration parameter, c, the spin parameter, λ, and the mean axis ratio, . For the concentration–mass relation we find c∝ in agreement with the model proposed by Bullock et al., but inconsistent with the alternative model of Eke et al. The normalization of the concentration–mass relation, however, is 15 per cent lower than suggested by Bullock et al. The results for λ and are in good agreement with previous studies, when extrapolated to the lower halo masses probed here, while c and λ are anticorrelated, in that high-spin haloes have, on average, lower concentrations. In an attempt to remove unrelaxed haloes from the sample, we compute for each halo the offset parameter, xoff, defined as the distance between the most bound particle and the centre of mass, in units of the virial radius. Removing haloes with large xoff increases the mean concentration by ∼10 per cent, lowers the mean spin parameter by ∼15 per cent, and removes the most prolate haloes. In addition, it largely removes the anticorrelation between c and λ, though not entirely. We also investigate the relation between halo properties and their large-scale environment density. For low-mass haloes we find that more concentrated haloes live in denser environments than their less concentrated counterparts of the same mass, consistent with recent correlation function analyses. Note, however, that the trend is weak compared to the scatter. For the halo spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in slightly denser environments. Finally, using a simple model for disc galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased subset of haloes. Not only do these haloes have spin parameters that are larger than average, they also have concentration parameters that are ∼15 per cent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from low surface brightness rotation curves, and those expected for a ΛCDM cosmology.
We have carried out a comparison study of hydrodynamical codes by investigating their performance in modelling interacting multiphase fluids. The two commonly used techniques of grid and smoothed ...particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin–Helmholtz or Rayleigh–Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of an SPH smoothing kernel radius over which interactions are severely damped.
Haloes gone MAD: The Halo-Finder Comparison Project Knebe, Alexander; Knollmann, Steffen R.; Muldrew, Stuart I. ...
Monthly notices of the Royal Astronomical Society,
08/2011, Letnik:
415, Številka:
3
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including ...friends‐of‐friends, spherical‐overdensity and phase‐space‐based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allow halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large‐scale structure of the universe.
All the halo‐finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo, and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Through a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30–40 particles. However, also here the phase‐space finders excelled by resolving substructure down to 10–20 particles. By comparing the halo finders using a high‐resolution cosmological volume, we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity and peak of the rotation curve).
We further suggest to utilize the peak of the rotation curve, vmax, as a proxy for mass, given the arbitrariness in defining a proper halo edge.
Merging haloes with similar masses (i.e. major mergers) pose significant challenges for halo finders. We compare five halo-finding algorithms’ (ahf, hbt, rockstar, subfind, and velociraptor) recovery ...of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z > 1). This raises significant uncertainties for theoretical models that require major mergers for, e.g. galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.