Search

To fully take advantage of the data provided by large-scale structure surveys, we need to quantify the potential impact of baryonic effects, such as feedback from active galactic nuclei (AGN) and star formation, on cosmological observables. In simulations, feedback processes originate on scales that remain unresolved. Therefore, they need to be sourced via subgrid models that contain free parameters. We use machine learning to calibrate the AGN and stellar feedback models for the FLAMINGO cosmological hydrodynamical simulations. Using Gaussian process emulators trained on Latin hypercubes of 32 smaller-volume simulations, we model how the galaxy stellar mass function and cluster gas fractions change as a function of the subgrid parameters. The emulators are then fit to observational data, allowing for the inclusion of potential observational biases. We apply our method to the three different FLAMINGO resolutions, spanning a factor of 64 in particle mass, recovering the observed relations within the respective resolved mass ranges. We also use the emulators, which link changes in subgrid parameters to changes in observables, to find models that skirt or exceed the observationally allowed range for cluster gas fractions and the stellar mass function. Our method enables us to define model variations in terms of the data that they are calibrated to rather than the values of specific subgrid parameters. This approach is useful, because subgrid parameters are typically not directly linked to particular observables, and predictions for a specific observable are influenced by multiple subgrid parameters., Comment: 24 pages, 10 figures (Including the appendix). Submitted to MNRAS. For visualisations, see the FLAMINGO website at https://flamingo.strw.leidenuniv.nl/

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 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 and baryonic particle masses of $1\times10^8$ and $1\times10^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\times10^{11}$ particles, making it the largest ever hydrodynamical simulation run to $z=0$. Lightcone 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 $\approx20$ per cent., Comment: 43 pages, 23 figures. Submitted to MNRAS. For visualisations, see the FLAMINGO website at https://flamingo.strw.leidenuniv.nl/

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations., Comment: The paper contains 23 pages, 13 figures. Accepted to MNRAS (2011, June 28)

We compare the results of two techniques used to calculate the evolution of cooling gas during galaxy formation: Smooth Particle Hydrodynamics (SPH) simulations and semi-analytic modelling. We improve upon the earlier statistical comparison of Benson et al. by taking halo merger histories from the dark matter component of the SPH simulation, which allows us to compare the evolution of galaxies on an object-by-object basis in the two treatments. We use a ``stripped-down'' version of the semi-analytic model described by Helly et al. which includes only shock heating and radiative cooling of gas and which is adjusted to mimic the resolution and other parameters of a comparison SPH simulation as closely as possible. We compare the total mass of gas that cools in halos of different mass as a function of redshift as well as the masses and spatial distribution of individual ``galaxies.'' At redshift z=0, the cooled gas mass in well-resolved halos agrees remarkably well (to better than ~20%) in the SPH simulation and stripped-down semi-analytic model. At high redshift, resolution effects in the simulation become increasingly important and, as a result, more gas tends to cool in low mass halos in the SPH simulation than in the semi-analytic model. The cold gas mass function of individual galaxies in the two treatments at z=0 also agrees very well and, when the effects of mergers are accounted for, the masses of individual galaxies and their 2-point correlation functions are also in excellent agreement in the two treatments. Thus, our comparison confirms and extends the earlier conclusion of Benson et al. that SPH simulations and semi-analytic models give consistent results for the evolution of cooling galactic gas., Comment: 14 pages, 9 figures, accepted by MNRAS. This paper was split at the referees request - the semi-analytic methods used here are now described in a separate paper, astro-ph/0210141