ABSTRACT
We introduce a new self-consistent model of galaxy evolution and reionization, astraeus (seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr ...simUlationS), which couples a state-of-the-art N-body simulation with the semi-analytical galaxy evolution delphi and the seminumerical reionization scheme cifog. astraeus includes all the key processes of galaxy formation and evolution (including accretion, mergers, supernova, and radiative feedback) and follows the time and spatial evolution of the ionized regions in the intergalactic medium (IGM). Importantly, it explores different radiative feedback models that cover the physically plausible parameter space, ranging from a weak and delayed to a strong and immediate reduction of gas mass available for star formation. From our simulation suite that covers the different radiative feedback prescriptions and ionization topologies, we find that radiative feedback continuously reduces star formation in galaxies with $M_\mathrm{ h}\lesssim 10^{9.5}\, {\rm \rm M_\odot }$ upon local reionization; larger mass haloes are unaffected even for the strongest and immediate radiative feedback cases during reionization. For this reason, the ionization topologies of different radiative feedback scenarios differ only on scales smaller than 1–2 comoving Mpc, and significant deviations are found only when physical parameters (e.g. the escape fraction of ionizing photons) are altered based on galactic properties. Finally, we find that observables (the ultraviolet luminosity function, stellar mass function, reionization histories and ionization topologies) are hardly affected by the choice of the used stellar population synthesis models that model either single stars or binaries.
Predicting structural properties of dark matter haloes is one of the fundamental goals of modern cosmology. We use the suite of MultiDark cosmological simulations to study the evolution of dark ...matter halo density profiles, concentrations, and velocity anisotropies. We find that in order to understand the structure of dark matter haloes and to make 1–2 per cent accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the ratio of the virial radius to the core radius in the Navarro–Frenk–White (NFW) profile. For massive haloes, the average density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter α in the Einasto approximation. We show that haloes progress through three stages of evolution. They start as rare density peaks and experience fast and nearly radial infall that brings mass closer to the centre, producing a highly concentrated halo. Here, the halo concentration increases with increasing halo mass and the concentration is defined by the α parameter with a nearly constant core radius. Later haloes slide into the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. At this stage, the concentration does not depend on halo mass. Once the rate of accretion and merging slows down, haloes move into the domain of declining concentration–mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. Accurate analytical fits are provided.
TURNING AROUND ALONG THE COSMIC WEB Lee, Jounghun; Yepes, Gustavo
The Astrophysical journal,
12/2016, Letnik:
832, Številka:
2
Journal Article
Recenzirano
Odprti dostop
ABSTRACT A bound violation designates a case in which the turnaround radius of a bound object exceeds the upper limit imposed by the spherical collapse model based on the standard ΛCDM paradigm. ...Given that the turnaround radius of a bound object is a stochastic quantity and that the spherical model overly simplifies the true gravitational collapse, which actually proceeds anisotropically along the cosmic web, the rarity of the occurrence of a bound violation may depend on the web environment. Assuming a Planck cosmology, we numerically construct the bound-zone peculiar velocity profiles along the cosmic web (filaments and sheets) around the isolated groups with virial mass identified in the Small MultiDark Planck simulations and determine the radial distances at which their peculiar velocities equal the Hubble expansion speed as the turnaround radii of the groups. It is found that although the average turnaround radii of the isolated groups are well below the spherical bound limit on all mass scales, the bound violations are not forbidden for individual groups, and the cosmic web has an effect of reducing the rarity of the occurrence of a bound violation. Explaining that the spherical bound limit on the turnaround radius in fact represents the threshold distance up to which the intervention of the external gravitational field in the bound-zone peculiar velocity profiles around the nonisolated groups stays negligible, we discuss the possibility of using the threshold distance scale to constrain locally the equation of state of dark energy.
Abstract
The halo shape plays a central role in determining important observational properties of the haloes such as mass, concentration and lensing cross-sections. The triaxiality of lensing galaxy ...clusters has a substantial impact on the distribution of the largest Einstein radii, while weak lensing techniques are sensitive to the intrinsic halo ellipticity. In this work, we provide scaling relations for the shapes of dark matter haloes as a function of mass (peak height) and redshift over more than four orders of magnitude in halo masses, namely from 1011.5 to 1015.8 h
−1 M⊙. We have analysed four dark matter only simulations from the MultiDark cosmological simulation suite with more than 56 billion particles within boxes of 4.0, 2.5, 1.0 and 0.4 h
−1 Gpc size assuming Planck cosmology. The dark matter haloes have been identified in the simulations using the rockstar halo finder, which also determines the axis ratios in terms of the diagonalization of the inertia tensor. In order to infer the shape for a hypothetical halo of a given mass at a given redshift, we provide fitting functions to the minor-to-major and intermediate-to-major axis ratios as a function of the peak height.
ABSTRACT
Cosmic Dawn II (CoDa II) is a new, fully coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift z < 6. With 40963 ...particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above 108 M⊙. Using the same hybrid CPU–GPU code RAMSES–CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star formation algorithm, making reionization end earlier, at z ≳ 6, thereby better matching the observations of intergalactic Lyman α opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at z = 6. As in CoDa I, reionization feedback suppresses star formation in haloes below ∼2 × 109 M⊙, though suppression here is less severe, a possible consequence of modifying the star formation algorithm. Suppression is environment dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of lambda cold dark matter constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, or by nearby groups, or by each other.
ABSTRACT
The James Webb Space Telescope (JWST) is expected to observe galaxies at z > 10 that are presently inaccessible. Here, we use a self-consistent empirical model, the universemachine, to ...generate mock galaxy catalogues and light-cones over the redshift range z = 0−15. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy–halo relationships, and galaxy–galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at z < 10. We predict that Cycle 1 JWST surveys will very likely detect galaxies with M* > 107 M⊙ and/or M1500 < −17 out to at least z ∼ 13.5. Number density uncertainties at z > 12 expand dramatically, so efforts to detect z > 12 galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits (M* < 107 M⊙ or M1500 > −17). For reionization models, extrapolating luminosity functions with a constant faint-end slope from M1500 = −17 down to M1500 = −12 gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to z ∼ 12. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and light-cones for common fields are available online.
ABSTRACT
In this work, we have implemented a detailed physical model of galaxy chemical enrichment into the Astraeus (seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in ...N-body dark matter simUlationS) framework which couples galaxy formation and reionization in the first billion years. Simulating galaxies spanning over 2.5 orders of magnitude in halo mass with $M_{\rm h} \sim 10^{8.9}{-}10^{11.5}\,{\rm M_\odot}$ ($M_{\rm h} \sim 10^{8.9}{-}10^{12.8}\rm M_\odot$) at z ∼ 10 (5), we find: (i) smooth accretion of metal-poor gas from the intergalactic medium (IGM) plays a key role in diluting the interstellar medium interstellar medium metallicity which is effectively restored due to self-enrichment from star formation; (ii) a redshift averaged gas-mass loading factor that depends on the stellar mass as $\eta _{\rm g} \approx 1.38 ({M_*}/{10^{10}\, {\rm \rm M_\odot }})^{-0.43}$; (iii) the mass–metallicity relation is already in place at z ∼ 10 and shows effectively no redshift evolution down to z ∼ 5; (iv) for a given stellar mass, the metallicity decreases with an increase in the star formation rate (SFR); (v) the key properties of the gas-phase metallicity (in units of 12 + log(O/H), stellar mass, SFR and redshift are linked through a high-redshift fundamental plane of metallicity (HFPZ) for which we provide a functional form; (vi) the mass–metallicity–SFR relations are effectively independent of the reionization radiative feedback model for $M_* {\,\, \buildrel\gt \over \sim \,\,}10^{6.5}\rm M_\odot$ galaxies; (vii) while low-mass galaxies ($M_{\rm h} {\,\, \buildrel\lt \over \sim \,\,}10^9\,\rm M_\odot$) are the key contributors to the metal budget of the IGM at early times, higher mass haloes provide about 50 per cent of the metal budget at lower redshifts.
Abstract
Radio relics in galaxy clusters are believed to trace merger shock fronts. If cosmological structure formation determines the luminosity, size and shape distributions of radio relics, then ...merger shocks need to be lighted up in a homogeneous way. We investigate if a mock relic sample, obtained from zoomed galaxy cluster simulations, is able to match the properties of relics measured in the NRAO VLA Sky Survey (NVSS). We compile a list of all radio relics known to date and homogeneously measure their parameters in all NVSS images and apply the same procedure to relics in our simulations. Number counts in the mock relic sample increase more steeply towards lower relic flux densities, suggesting an incompleteness of NVSS in this regime. Overall, we find that NVSS and mock samples show similar properties. However, large simulated relics tend to be somewhat smaller and closer to the cluster centre than observed ones. Besides this, the mock sample reproduces very well-known correlations for radio relics, in particular those relating the radio luminosity with the largest linear size and the X-ray luminosity. We show that these correlations are largely governed by the sensitivity of the NVSS observations. Mock relics show a similar orientation with respect to the direction to the cluster centre as the NVSS sample. Moreover, we find that their maximum radio luminosity roughly correlates with cluster mass, although displaying a large scatter. The overall good agreement between NVSS and the mock sample suggests that properties of radio relics are indeed governed by merger shock fronts, emitting in a homogeneous fashion. Our study demonstrates that the combination of mock observations and data from upcoming radio surveys will allow us to shed light on both the origin of radio relics and the nature of the intracluster medium.
We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat Delta *LCDM cosmology. Halos are identified using the spherical overdensity ...algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta *D times the background. We calibrate fitting functions for the large-scale bias that are adaptable to any value of Delta *D we examine. We find a ~6% scatter about our best-fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et al. such that the bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth et al. and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta *D = 200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.
Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization itself. Star formation suppression, for example, may explain the observed underabundance ...of Local Group dwarfs relative to N-body predictions for cold dark matter. Reionization modelling requires simulating volumes large enough ~ (100 Mpc) super( 3) to sample reionization 'patchiness', while resolving millions of galaxy sources above ~10 super( 8) M... combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the Local Universe today. Cosmic Dawn ('CoDa') is the first such fully coupled, radiation-hydrodynamics simulation of reionization of the Local Universe. Our new hybrid CPU-GPU code, ramses-cudaton, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc) super( 3) with 4096 super( 3) particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10 super( 9) M..., decreasing the abundance of faint galaxies around M sub( AB1600) = -10, -12. For most of reionization, star formation was dominated by haloes between 10 super( 10)-10 super( 11) M..., so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. Intergalactic filaments display sheathed structures, with hot envelopes surrounding cooler cores, but do not self-shield, unlike regions denser than 100 ... (ProQuest: ... denotes formulae/symbols omitted.)