Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHow to Zoom: Bias, Contamination, and Lagrange Volumes in Multimass Cosmological Simulations
225
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We perform a suite of multimass cosmological zoom simulations of individual dark matter halos and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halos Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that halos living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting halos to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose halos with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing ones sample.
rate research
Read More
We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved H{alpha} and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in H{alpha} but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.
The advent of a new generation of large-scale galaxy surveys is pushing cosmological numerical simulations in an uncharted territory. The simultaneous requirements of high resolution and very large volume pose serious technical challenges, due to their computational and data storage demand. In this paper, we present a novel approach dubbed Dynamic Zoom Simulations -- or DZS -- developed to tackle these issues. Our method is tailored to the production of lightcone outputs from N-body numerical simulations, which allow for a more efficient storage and post-processing compared to standard comoving snapshots, and more directly mimic the format of survey data. In DZS, the resolution of the simulation is dynamically decreased outside the lightcone surface, reducing the computational work load, while simultaneously preserving the accuracy inside the lightcone and the large-scale gravitational field. We show that our approach can achieve virtually identical results to traditional simulations at half of the computational cost for our largest box. We also forecast this speedup to increase up to a factor of 5 for larger and/or higher-resolution simulations. We assess the accuracy of the numerical integration by comparing pairs of identical simulations run with and without DZS. Deviations in the lightcone halo mass function, in the sky-projected lightcone, and in the 3D matter lightcone always remain below 0.1%. In summary, our results indicate that the DZS technique may provide a highly-valuable tool to address the technical challenges that will characterise the next generation of large-scale cosmological simulations.
We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.
Cosmological N-body simulations have been a major tool of theorists for decades, yet many of the numerical issues that these simulations face are still unexplored. This paper measures numerical biases in these large, dark matter-only simulations that affect the properties of their dark matter haloes. We compare many simulation suites in order to provide several tools for simulators and analysts which help mitigate these biases. We summarise our comparisons with practical `convergence limits that can be applied to a wide range of halo properties, including halo properties which are traditionally overlooked by the testing literature. We also find that the halo properties predicted by different simulations can diverge from one another at unexpectedly high resolutions. We demonstrate that many halo properties depend strongly on force softening scale and that this dependence leads to much of the measured divergence between simulations. We offer an empirical model to estimate the impact of such effects on the rotation curves of a halo population. This model can serve as a template for future empirical models of the biases in other halo properties.
We present cosmological zoom-in hydro-dynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain size distribution, with the MUPPI star formation and feedback sub-resolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio and small-to-big grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
