Abstract
Recent studies have presented evidence for tension between the constraints on Ωm and σ8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This ...tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully calibrated cosmological hydrodynamical simulations, which thus capture the backreaction of the baryons on the total matter distribution. We have extended the BAryons and HAloes of MAssive Sysmtes simulations to include a treatment of massive neutrinos, which currently represents the best-motivated extension to the standard model. We make synthetic thermal Sunyaev–Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: (i) in general, there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; (ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and (iii) if one accounts for internal tensions in the Planck CMB data set (by allowing the lensing amplitude, ALens, to vary), invoking a non-minimal neutrino mass, typically of 0.2–0.4 eV, can resolve the tension. This solution is fully consistent with separate constraints from the primary CMB and baryon acoustic oscillations.
ABSTRACT
The recent advent of deep observational surveys of local Milky Way ‘analogues’ and their satellite populations allows us to place the Milky Way in a broader cosmological context and to test ...models of galaxy formation on small scales. In this study, we use the Lambda cold dark matter (ΛCDM)-based ARTEMIS suite of cosmological hydrodynamical simulations containing 45 Milky Way analogue host haloes to make comparisons to the observed satellite luminosity functions, radial distribution functions, and abundance scaling relations from the recent Local Volume and SAGA observational surveys, in addition to the Milky Way and M31. We find that, contrary to some previous claims, ΛCDM-based simulations can successfully and simultaneously capture the mean trends and the diversity in both the observed luminosity and radial distribution functions of Milky Way analogues once important observational selection criteria are factored in. Furthermore, we show that, at fixed halo mass, the concentration of the simulated satellite radial distribution is partly set by that of the underlying smooth dark matter halo, although stochasticity due to the finite number of satellites is the dominant driver of scatter in the radial distribution of satellites at fixed halo mass.
ABSTRACT
We introduce the Assembly of high-ResoluTion Eagle-simulations of MIlky Way-type galaxieS (artemis) simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of $\approx ...2\times 10^4 \, {\rm M}_{\odot }\, h^{-1}$), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the eagle galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour, and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ≈−3, outer slopes of ≈−4, and break radii that are typically ≈20–40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour, and age profiles show mild large-scale gradients, particularly when spherically averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes.
ABSTRACT
Using the artemis set of 45 high-resolution cosmological simulations, we investigate a range of merger-induced dynamical transformations of Milky Way-like galaxies. We first identify ...populations of accreted stars on highly radial orbits, similar to the ‘Gaia Sausage’ in the Milky Way. We show that ≈1/3 of the artemis galaxies contain a similar feature, and confirm that they usually comprise stellar debris from the most massive accreted satellite. Selecting 15 galaxies with discs at the present-day, we study their changes around the times of the GS-like mergers. Dark matter haloes of many of these exhibit global changes in shape and orientation, with almost half becoming significantly more spherical when the mergers occur. Focusing on the galaxies themselves, we find that 4/15 have stellar discs which experience large changes in the orientation of their angular momentum (AM) axes, at rates of up to ∼60 degrees Gyr−1. By calculating the orbital angular momentum axes of the satellites before they are accreted, we show that there is a tendency for the disc’s AM to become more aligned with this axis after the merger. We also investigate the origin of in situ retrograde stars, analogous to the ‘Splash’ in the Milky Way. Tracing them back to earlier snapshots, we demonstrate that they were often disrupted on to their extreme orbits by multiple early mergers. We also find that the total mass of these stars outside the central regions positively correlates with the total accreted stellar mass.
Abstract
We quantify the correlations between gas-phase and stellar metallicities and global properties of galaxies, such as stellar mass, halo mass, age and gas fraction, in the Evolution and ...Assembly of GaLaxies and their Environments suite of cosmological hydrodynamical simulations. The slope of the correlation between stellar mass and metallicity of star-forming (SF) gas (M
*–Z
SF,gas relation) depends somewhat on resolution, with the higher resolution run reproducing a steeper slope. This simulation predicts a non-zero metallicity evolution, increasing by ≈0.5 dex at ∼109 M⊙ since z = 3. The simulated relation between stellar mass, metallicity and star formation rate at z ≲ 5 agrees remarkably well with the observed fundamental metallicity relation. At M
* ≲ 1010.3 M⊙ and fixed stellar mass, higher metallicities are associated with lower specific star formation rates, lower gas fractions and older stellar populations. On the other hand, at higher M
*, there is a hint of an inversion of the dependence of metallicity on these parameters. The fundamental parameter that best correlates with the metal content, in the simulations, is the gas fraction. The simulated gas fraction–metallicity relation exhibits small scatter and does not evolve significantly since z = 3. In order to better understand the origin of these correlations, we analyse a set of lower resolution simulations in which feedback parameters are varied. We find that the slope of the simulated M
*–Z
SF,gas relation is mostly determined by stellar feedback at low stellar masses (M
* ≲ 1010 M⊙), and at high masses (M
* ≳ 1010 M⊙) by the feedback from active galactic nuclei.
In the local Universe, galaxies in groups and clusters contain less gas and are less likely to be forming stars than their field counterparts. This effect is not limited to the central group/cluster ...regions, but is shown by recent observations to persist out to several virial radii. To gain insight into the extent and cause of this large-scale environmental influence, we use a suite of high-resolution cosmological hydrodynamic simulations to analyse galaxies around simulated groups and clusters of a wide range of mass (log 10
M
host/M = 13.0, 15.2). In qualitative agreement with the observations, we find a systematic depletion of both hot and cold gas and a decline in the star-forming fraction of galaxies as far out as ∼5r
200 from the host centre. A substantial fraction of these galaxies are on highly elliptical orbits and are not infalling for the first time (∼50 per cent at 2r
200, independent of host mass) or are affected by 'pre-processing' (less than 10 per cent of galaxies around groups, increasing to ∼50 per cent near a massive cluster), even a combination of these indirect mechanisms does not fully account for the environmental influence, particularly in the case of the hot gas content. Direct ram pressure interaction with an extended gas 'halo' surrounding groups and clusters is shown to be sufficiently strong to strip the hot gas atmospheres of infalling galaxies out to ∼5r
200. We show that this influence is highly anisotropic, with ram pressure along filaments enhanced by up to a factor of 100 despite significant co-flow of gas and galaxies.
ABSTRACT
We investigate the formation (spin-up) of galactic discs in the artemis simulations of Milky Way (MW)-mass galaxies. In almost all galaxies, discs spin up at higher Fe/H than the MW. Those ...galaxies that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which formed their discs ∼8–11 Gyr ago. Five of these experience a merger similar to the GSE. The spin-up times correlate with the halo masses at early times: galaxies with early spin-up have larger virial masses at a lookback time tL = 12 Gyr. The fraction of stars accreted from outside the host galaxy is smaller in galaxies with earlier spin-ups. Accreted fractions small enough to be comparable to the MW are only found in galaxies with the earliest disc formation and large initial virial masses (M200c ≈ 2 × 1011 M⊙ at tL = 12 Gyr). We find that discs form when the halo’s virial mass reaches a threshold of M200c ≈ (6 ± 3) × 1011 M⊙, independent of the spin-up time. However, the failure to form a disc in other galaxies appears to be instead related to mergers at early times. We also find that discs form when the central potential is not particularly steep. Our results indicate that the MW assembled its mass and formed its disc earlier than the average galaxy of a similar mass.
ABSTRACT
We present redshift-zero synthetic dust-aware observations for the 45 Milky Way-mass simulated galaxies of the ARTEMIS project, calculated with the skirt radiative transfer code. The ...post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust. We produce and publicly release realistic high-resolution images for 50 commonly used broad-band filters from ultraviolet to submillimetre wavelengths and for 18 different viewing angles. We compare the simulated ARTEMIS galaxies to observed galaxies in the DustPedia data base with similar stellar mass and star formation rate, and to synthetic observations of the simulated galaxies of the Auriga project produced in previous work using a similar post-processing technique. In all cases, global galaxy properties are derived using SED fitting. We find that, similar to Auriga, the post-processed ARTEMIS galaxies generally reproduce the observed scaling relations for global fluxes and physical properties, although dust extinction at FUV/UV wavelengths is underestimated and representative dust temperatures are lower than observed. At a resolved scale, we compare multiwavelength non-parametric morphological properties of selected disc galaxies across the data sets. We find that the ARTEMIS galaxies largely reproduce the observed morphological trends as a function of wavelength, although they appear to be more clumpy and less symmetrical than observed. We note that the ARTEMIS and Auriga galaxies occupy adjacent regions in the specific star formation versus stellar mass plane, so that the synthetic observation data sets supplement each other.
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
The vast majority of low-mass satellite galaxies around the Milky Way and M31 appear virtually devoid of cool gas and show no signs of recent or ongoing star formation. Cosmological ...simulations demonstrate that such quenching is expected and is due to the harsh environmental conditions that satellites face when joining the Local Group (LG). However, recent observations of Milky Way analogues in the SAGA survey present a very different picture, showing the majority of observed satellites to be actively forming stars, calling into question the realism of current simulations and the typicality of the LG. Here, we use the ARTEMIS suite of high-resolution cosmological hydrodynamical simulations to carry out a careful comparison with observations of dwarf satellites in the LG, SAGA, and the Local Volume (LV) survey. We show that differences between SAGA and the LG and LV surveys, as well as between SAGA and the ARTEMIS simulations, can be strongly reduced by considering differences in the host mass distributions and (more importantly) observational selection effects, specifically that low-mass satellites which have only recently been accreted are more likely to be star forming, have a higher optical surface brightness, and are therefore more likely to be included in the SAGA survey. This picture is confirmed using data from the deeper LV survey, which shows pronounced quenching at low masses, in accordance with the predictions of ΛCDM-based simulations.