Cosmology

Abstract:

Galactic scale tests have proven to be powerful tools in constraining fundamental physics in previously under-explored regions of parameter space. The astrophysical regime which they probe is inherently complicated, and the inference methods used to make these constraints should be robust to baryonic effects. Previous analyses have assumed simple empirical models for astrophysical noise without detailed calibration or justification. We outline a framework for assessing the reliability of such methods by constructing and testing more advanced baryonic models using cosmological hydrodynamical simulations. As a case study, we use the Horizon-AGN simulation to investigate warping of stellar disks and offsets between gas and stars within galaxies, which are powerful probes of screened fifth forces. We show that the degree of `U'-shaped warping of galaxies is well modelled by Gaussian random noise, but that the magnitude of the gas-star offset is correlated with the virial radius of the host halo. By incorporating this correlation we confirm recent results ruling out astrophysically relevant Hu-Sawicki $f(R)$ gravity, and identify a $\sim 30\%$ systematic uncertainty due to baryonic physics. Such an analysis must be performed case-by-case for future galactic tests of fundamental physics.

Abstract:

The anisotropic distribution of satellites around the central galaxy of their host halo is well-documented. However the relative impact of baryons and dark matter in shaping this distribution is still debated. Using the simulation Horizon-AGN, the angular distribution of satellite galaxies with respect to their central counterpart and halo is quantified. Below one Rvir, satellites cluster more strongly in the plane of the central, rather than merely tracing the shape of their host halo. This is due to the increased isotropy of inner haloes acquired through their inside-out assembly in vorticity-rich flows along the cosmic web. While the effect of centrals decreases with distance, halos' triaxiality increases, impacting more and more the satellite's distribution. Effects become comparable just outside one virial radius. Above this scale, the filamentary infall also impacts the satellites distribution, dominating above two virial radii. The central's morphology plays a governing role: the alignment w.r.t. the central plane is four times stronger in haloes hosting stellar discs than in spheroids. But the impact of the galactic plane decreases for lower satellite-to-central mass ratios, suggesting this might not hold for dwarf satellites of the Local group. The orientation of the Milky-Way's satellites traces their cosmic filament, their level of coplanarity is consistent with systems of similar mass and cosmic location in Horizon-AGN. However, the strong impact of galactic planes in massive groups and clusters bounds the likelihood of finding a relaxed region where satellites can be used to infer halo shape. The minor-to-major axis ratios for haloes with log(M0/Msun)>13.5 is underestimated by 10%. This error soars quickly to 30-40% for individual halo measurements.

Abstract:

Hydrodynamical cosmological simulations have recently made great advances in reproducing galaxy mass assembly over cosmic time - as often quantified from the comparison of their predicted stellar mass functions to observed stellar mass functions from data. In this paper we compare the clustering of galaxies from the hydrodynamical cosmological simulated lightcone Horizon-AGN, to clustering measurements from the VIDEO survey observations. Using mocks built from a VIDEO-like photometry, we first explore the bias introduced into clustering measurements by using stellar masses and redshifts derived from SED-fitting, rather than the intrinsic values. The propagation of redshift and mass statistical and systematic uncertainties in the clustering measurements causes us to underestimate the clustering amplitude. We find then that clustering and halo occupation distribution (HOD) modelling results are qualitatively similar in Horizon-AGN and VIDEO. However at low stellar masses Horizon-AGN underestimates the observed clustering by up to a factor of ~3, reflecting the known excess stellar mass to halo mass ratio for Horizon-AGN low mass haloes, already discussed in previous works. This reinforces the need for stronger regulation of star formation in low mass haloes in the simulation. Finally, the comparison of the stellar mass to halo mass ratio in the simulated catalogue, inferred from angular clustering, to that directly measured from the simulation, validates HOD modelling of clustering as a probe of the galaxy-halo connection.

Abstract:

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

Abstract:

We present a public suite of weak lensing mock data, extending the Scinet Light Cone Simulations (SLICS) to simulate cross-correlation analyses with different cosmological probes. These mocks include KiDS-450- and LSST-like lensing data, cosmic microwave background lensing maps and simulated spectroscopic surveys that emulate the GAMA, BOSS and 2dFLenS galaxy surveys. With 817 independent realisations, our mocks are optimised for combined-probe covariance estimation, which we illustrate for the case of a joint measurement involving cosmic shear, galaxy-galaxy lensing and galaxy clustering from KiDS-450 and BOSS data. With their high spatial resolution, the SLICS are also optimal for predicting the signal for novel lensing estimators, for the validation of analysis pipelines, and for testing a range of systematic effects such as the impact of neighbour-exclusion bias on the measured tomographic cosmic shear signal. For surveys like KiDS and DES, where the rejection of neighbouring galaxies occurs within ~2 arcseconds, we show that the measured cosmic shear signal will be biased low, but by less than a percent on the angular scales that are typically used in cosmic shear analyses. The amplitude of the neighbour-exclusion bias doubles in deeper, LSST-like data. The simulation products described in this paper are made available at http://slics.roe.ac.uk/.