ABSTRACT
The halo model formalism is widely adopted in cosmological studies for predicting the growth of large-scale structure in the Universe. However, to date there have been relatively few direct ...comparisons of the halo model with more accurate (but much more computationally expensive) cosmological simulations. We test the accuracy of the halo model in reproducing the non-linear matter power spectrum, P(k), when the main inputs of the halo model (specifically the matter density profiles, halo mass function, and linear bias) are taken directly from the BAHAMAS simulations, and we assess how well the halo model reproduces P(k) from the same simulations. We show that the halo model generally reproduces P(k) in the deep non-linear regime (1-halo) to typically a few per cent accuracy, but struggles to reproduce (approx. 15 per cent error) P(k) at intermediate scales of $0.1 \lesssim k \ h/{\rm Mpc} \lesssim 3$ at $z$ = 0, marking the transition between the 1-halo and 2-halo terms. We show that the magnitude of this error is a strong function of the halo mass definition (through its effects on radial extent of haloes) and of redshift. Furthermore, we test the accuracy of the halo model in recovering the relative impact of baryons on P(k). We show that the systematic errors in recovering the absolute P(k) largely cancel when considering the relative impact of baryons. This suggests that the halo model can make precise predictions for the baryonic suppression, offering a fast and accurate way to adjust collisionless matter power spectra for the presence of baryons and associated processes.
We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser ...Interferometer Space Antenna (eLISA), utilizing the reference cosmological hydrodynamical simulation from the eagle suite. These simulations assume a Lambda cold dark matter cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in eagle with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The eagle models predict ∼2 detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between z ∼ 2 and z ∼ 1. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger time-scale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.
Abstract
We investigate the effect of the accelerated expansion of the Universe due to a cosmological constant, Λ, on the cosmic star formation rate. We utilize hydrodynamical simulations from the ...eagle suite, comparing a ΛCDM (cold dark matter) Universe to an Einstein–de Sitter model with Λ = 0. Despite the differences in the rate of growth of structure, we find that dark energy, at its observed value, has negligible impact on star formation in the Universe. We study these effects beyond the present day by allowing the simulations to run forward into the future (t > 13.8 Gyr). We show that the impact of Λ becomes significant only when the Universe has already produced most of its stellar mass, only decreasing the total comoving density of stars ever formed by ${\approx }15\hbox{ per cent}$. We develop a simple analytic model for the cosmic star formation rate that captures the suppression due to a cosmological constant. The main reason for the similarity between the models is that feedback from accreting black holes dramatically reduces the cosmic star formation at late times. Interestingly, simulations without feedback from accreting black holes predict an upturn in the cosmic star formation rate for t > 15 Gyr due to the rejuvenation of massive (>1011 M⊙) galaxies. We briefly discuss the implication of the weak dependence of the cosmic star formation on Λ in the context of the anthropic principle.
How feedback shapes galaxies: an analytic model Salcido, Jaime; Bower, Richard G; Theuns, Tom
Monthly Notices of the Royal Astronomical Society,
02/2020, Letnik:
491, Številka:
4
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We introduce a simple analytic model of galaxy formation that links the growth of dark matter haloes in a cosmological background to the build-up of stellar mass within them. The model aims ...to identify the physical processes that drive the galaxy-halo co-evolution through cosmic time. The model restricts the role of baryonic astrophysics to setting the relation between galaxies and their haloes. Using this approach, galaxy properties can be directly predicted from the growth of their host dark matter haloes. We explore models in which the effective star formation efficiency within haloes is a function of mass (or virial temperature) and independent of time. Despite its simplicity, the model reproduces self-consistently the shape and evolution of the cosmic star formation rate density, the specific star formation rate of galaxies, and the galaxy stellar mass function, both at the present time and at high redshifts. By systematically varying the effective star formation efficiency in the model, we explore the emergence of the characteristic shape of the galaxy stellar mass function. The origin of the observed double Schechter function at low redshifts is naturally explained by two efficiency regimes in the stellar to halo mass relation, namely, a stellar feedback regulated stage, and a supermassive black hole regulated stage. By providing a set of analytic differential equations, the model can be easily extended and inverted, allowing the roles and impact of astrophysics and cosmology to be explored and understood.
ABSTRACT
Upcoming large-scale structure surveys will measure the matter power spectrum to approximately per cent level accuracy with the aim of searching for evidence for new physics beyond the ...standard model of cosmology. In order to avoid biasing our conclusions, the theoretical predictions need to be at least as accurate as the measurements for a given choice of cosmological parameters. However, recent theoretical work has shown that complex physical processes associated with galaxy formation (particularly energetic feedback processes associated with stars and especially supermassive black holes) can alter the predictions by many times larger than the required accuracy. Here we present SP(k), a model for the effects of baryon physics on the non-linear matter power spectrum based on a new large suite of hydrodynamical simulations. Specifically, the ANTILLES suite consists of 400 simulations spanning a very wide range of the ‘feedback landscape’ and show that the effects of baryons on the matter power spectrum can be understood at approaching the per cent level in terms of the mean baryon fraction of haloes, at scales of up to $k \, {\lesssim } \, 10\, h$ Mpc−1 and redshifts up to z = 3. For the range of scales and redshifts that will be probed by forthcoming cosmic shear measurements, most of the effects are driven by galaxy group mass haloes (M ∼ 1013–14 M⊙). We present a simple python implementation of our model, available at https://github.com/jemme07/pyspk, which can be used to incorporate baryon effects in standard gravity-only predictions, allowing for marginalization over baryon physics within cosmological pipelines.
ABSTRACT
We introduce the Virgo Consortium’s FLAMINGO suite of hydrodynamical simulations for cosmology and galaxy cluster physics.
To ensure the simulations are sufficiently realistic for studies of ...large-scale structure, the subgrid prescriptions for stellar and AGN feedback are calibrated to the observed low-redshift galaxy stellar mass function and cluster gas fractions. The calibration is performed using machine learning, separately for each of FLAMINGO’s three resolutions. This approach enables specification of the model by the observables to which they are calibrated. The calibration accounts for a number of potential observational biases and for random errors in the observed stellar masses. The two most demanding simulations have box sizes of 1.0 and 2.8 Gpc on a side and baryonic particle masses of 1 × 108 and $1\times 10^9\, \text{M}_\odot$, respectively. For the latter resolution, the suite includes 12 model variations in a 1 Gpc box. There are 8 variations at fixed cosmology, including shifts in the stellar mass function and/or the cluster gas fractions to which we calibrate, and two alternative implementations of AGN feedback (thermal or jets). The remaining 4 variations use the unmodified calibration data but different cosmologies, including different neutrino masses. The 2.8 Gpc simulation follows 3 × 1011 particles, making it the largest ever hydrodynamical simulation run to z = 0. Light-cone output is produced on-the-fly for up to 8 different observers. We investigate numerical convergence, show that the simulations reproduce the calibration data, and compare with a number of galaxy, cluster, and large-scale structure observations, finding very good agreement with the data for converged predictions. Finally, by comparing hydrodynamical and ‘dark-matter-only’ simulations, we confirm that baryonic effects can suppress the halo mass function and the matter power spectrum by up to ≈20 per cent.
ABSTRACT
The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational ...methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev–Zel’dovich (SZ) observations of the hot intracluster gas. Here, we use the FLAMINGO simulations to investigate whether an extremum gradient in a similar position to the splashback radius is predicted to occur in the cluster gas profiles. We find that the minimum in the gradient of the stacked 3D gas density and pressure profiles, and the maximum in the gradient of the entropy profile, broadly align with the splashback feature though there are significant differences. While the dark matter splashback radius varies with specific mass accretion rate, in agreement with previous work, the radial position of the deepest minimum in the log-slope of the gas density is more sensitive to halo mass. In addition, we show that a similar minimum is also present in projected 2D pseudo-observable profiles: emission measure (X-ray), Compton-y (SZ), and surface mass density (weak lensing). We find that the latter traces the dark matter results reasonably well albeit the minimum occurs at a slightly smaller radius. While results for the gas profiles are largely insensitive to accretion rate and various observable proxies for dynamical state, they do depend on the strength of the feedback processes.
Abstract Galaxy clusters are important probes for both cosmology and galaxy formation physics. We test the cosmological, hydrodynamical FLAMINGO simulations by comparing to observations of the ...gaseous properties of clusters measured from X-ray observations. FLAMINGO contains unprecedented numbers of massive galaxy groups (>106) and clusters (>105) and includes variations in both cosmology and galaxy formation physics. We predict the evolution of cluster scaling relations as well as radial profiles of the temperature, density, pressure, entropy, and metallicity for different masses and redshifts. We show that the differences between volume-, and X-ray-weighting of particles in the simulations, and between cool-core non cool-core samples, are similar in size as the differences between simulations for which the stellar and AGN feedback has been calibrated to produce significantly different gas fractions. Compared to thermally-driven AGN feedback, kinetic jet feedback calibrated to produce the same gas fraction at R500c yields a hotter core with higher entropies and lower densities, which translates into a smaller fraction of cool-core clusters. Stronger feedback, calibrated to produce lower gas fractions and hence lower gas densities, results in higher temperatures, entropies, and metallicities, but lower pressures. The scaling relations and thermodynamic profiles show almost no evolution with respect to self-similar expectations, except for the metallicity decreasing with redshift. We find that the temperature, density, pressure, and entropy profiles of clusters in the fiducial FLAMINGO simulation are in excellent agreement with observations, while the metallicities in the core are too high.
ABSTRACT
Cosmology is poised to measure the neutrino mass sum Mν and has identified several smaller-scale observables sensitive to neutrinos, necessitating accurate predictions of neutrino clustering ...over a wide range of length scales. The FlowsForTheMasses non-linear perturbation theory for the the massive neutrino power spectrum, $\Delta ^2_\nu (k)$, agrees with its companion N-body simulation at the $10~{{\ \rm per\ cent}}-15~{{\ \rm per\ cent}}$ level for k ≤ 1 h Mpc−1. Building upon the Mira-Titan IV emulator for the cold matter, we use FlowsForTheMasses to construct an emulator for $\Delta ^2_\nu (k)$, Cosmic-Eν, which covers a large range of cosmological parameters and neutrino fractions Ων, 0h2 ≤ 0.01 (Mν ≤ 0.93 eV). Consistent with FlowsForTheMasses at the 3.5 per cent level, it returns a power spectrum in milliseconds. Ranking the neutrinos by initial momenta, we also emulate the power spectra of momentum deciles, providing information about their perturbed distribution function. Comparing a Mν = 0.15 eV model to a wide range of N-body simulation methods, we find agreement to 3 per cent for k ≤ 3kFS = 0.17 h Mpc−1 and to 19 per cent for k ≤ 0.4 h Mpc−1. We find that the enhancement factor, the ratio of $\Delta ^2_\nu (k)$ to its linear-response equivalent, is most strongly correlated with Ων, 0h2, and also with the clustering amplitude σ8. Furthermore, non-linearities enhance the free-streaming-limit scaling $\partial \log (\Delta ^2_\nu /\Delta ^2_{\rm m}) / \partial \log (M_\nu)$ beyond its linear value of 4, increasing the Mν-sensitivity of the small-scale neutrino density.
ABSTRACT
A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological ...parameters fit to the cosmic microwave background (CMB). The origin of this ‘S8 tension’ remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in the observational measurements and/or uncertainties in the role that baryons play. Here, we carefully examine the latter possibility using the new FLAMINGO suite of large-volume cosmological hydrodynamical simulations. We project the simulations onto observable harmonic space and compare with observational measurements of the power and cross-power spectra of cosmic shear, CMB lensing, and the thermal Sunyaev-Zel’dovich (tSZ) effect. We explore the dependence of the predictions on box size and resolution and cosmological parameters, including the neutrino mass, and the efficiency and nature of baryonic ‘feedback’. Despite the wide range of astrophysical behaviours simulated, we find that baryonic effects are not sufficiently large to remove the S8 tension. Consistent with recent studies, we find the CMB lensing power spectrum is in excellent agreement with the standard model, while the cosmic shear power spectrum, tSZ effect power spectrum, and the cross-spectra between shear, CMB lensing, and the tSZ effect are all in varying degrees of tension with the CMB-specified standard model. These results suggest that some mechanism is required to slow the growth of fluctuations at late times and/or on non-linear scales, but that it is unlikely that baryon physics is driving this modification.