Search

Most of the visible matter in the Universe is in a gaseous state, subject to hydrodynamic forces and galaxy formation processes that are much more complex to model than gravity. These baryonic effects can potentially bias the analyses of several cosmological probes, such as weak gravitational lensing. In this work, we study the gas density and velocity fields of the FLAMINGO simulations and compare them with their gravity-only predictions. We find that, while the gas velocities do not differ from those of dark matter on large scales, the gas mass power spectrum is suppressed by up to $\approx 8\%$ relative to matter, even on gigaparsec scales. This is a consequence of star formation depleting gas in the densest and most clustered regions of the universe. On smaller scales, $k>0.1 \, h / \rm Mpc$, the power suppression for both gas densities and velocities is more significant and correlates with the strength of the active galactic nucleus (AGN) feedback. The impact of feedback can be understood in terms of outflows, identified as gas bubbles with positive radial velocities ejected from the central galaxy. With increasing feedback strength, the outflowing gas has higher velocities, and it reaches scales as large as $10$ times the virial radius of the halo, redistributing the gas and slowing its average infall velocity. Interestingly, different implementations of AGN feedback leave distinct features in these outflows in terms of their radial and angular profiles and their dependence on halo mass. In the future, such differences could be measured in observations using, for example, the kinetic Sunyaev-Zeldovich effect.

Energetic feedback processes associated with accreting supermassive black holes can expel gas from massive haloes and significantly alter various measures of clustering on ~Mpc scales, potentially biasing the values of cosmological parameters inferred from analyses of large-scale structure (LSS) if not modelled accurately. Here we use the state-of-the-art FLAMINGO suite of cosmological hydrodynamical simulations to gauge the impact of feedback on large-scale structure by comparing to Planck + ACT stacking measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect of SDSS BOSS galaxies. We make careful like-with-like comparisons to the observations, aided by high precision KiDS and DES galaxy-galaxy lensing measurements of the BOSS galaxies to inform the selection of the simulated galaxies. In qualitative agreement with several recent studies using dark matter only simulations corrected for baryonic effects, we find that the kSZ effect measurements prefer stronger feedback than predicted by simulations which have been calibrated to reproduce the gas fractions of low redshift X-ray-selected groups and clusters. We find that the increased feedback can help to reduce the so-called S8 tension between the observed and CMB-predicted clustering on small scales as probed by cosmic shear (although at the expense of agreement with the X-ray group measurements). However, the increased feedback is only marginally effective at reducing the reported offsets between the predicted and observed clustering as probed by the thermal SZ (tSZ) effect power spectrum and tSZ effect--weak lensing cross-spectrum, both of which are sensitive to higher halo masses than cosmic shear., Comment: 20 pages, 10 figures, submitted to MNRAS

Recent large-scale structure (LSS) surveys have revealed a persistent tension in the value of $S_8$ compared to predictions from the standard cosmological model. This tension may suggest the need for new physics beyond the standard model, but an accurate characterisation of baryonic effects is essential to avoid biases. Although some studies indicate that baryonic effects are too small to resolve this tension, others propose that more aggressive feedback mechanisms could reconcile differences between cosmic microwave background (CMB) measurements and low-redshift LSS observations. In this paper, we investigate the role of baryonic effects in alleviating the $S_8$ tension. We extend the SP(k) model (Salcido et al. 2023), which was trained on hundreds of cosmological hydrodynamical simulations to map the suppression of the matter power spectrum to the baryon fraction in groups and clusters, to predict the required baryon fraction for a given $P(k)$ suppression. We then compare predictions from recent cosmic shear (weak lensing) analyses with the latest baryon budget measurements from X-ray and weak gravitational lensing studies. Our findings show that studies marginalising over baryonic effects while fixing cosmological parameters to a Planck-like cosmology predict strong $P(k)$ suppression and baryon fractions that are much lower than existing low-redshift baryon budget estimates of galaxy groups and clusters. Conversely, most studies that marginalise over both cosmological parameters and baryonic effects imply baryon fractions that are consistent with observations but lower values of $S_8$ than inferred from the CMB. Unless the observed baryon fractions are biased high by a factor of several, these results suggest that a mechanism beyond baryonic physics alone is required to modify or slow down the growth of structure in the universe in order to resolve the $S_8$ tension., Comment: Submitted to MNRAS. 12 pages, 7 figures

Galaxy clusters provide an avenue to expand our knowledge of cosmology and galaxy evolution. Because it is difficult to accurately measure the total mass of a large number of individual clusters, cluster samples are typically selected using an observable proxy for mass. Selection effects are therefore a key problem in understanding galaxy cluster statistics. We make use of the $(2.8~\rm{Gpc})^3$ FLAMINGO hydrodynamical simulation to investigate how selection based on X-ray luminosity, thermal Sunyaev-Zeldovich effect or galaxy richness influences the halo mass distribution. We define our selection cuts based on the median value of the observable at a fixed mass and compare the resulting samples to a mass-selected sample. We find that all samples are skewed towards lower mass haloes. For X-ray luminosity and richness cuts below a critical value, scatter dominates over the trend with mass and the median mass becomes biased increasingly low with respect to a mass-selected sample. At $z\leq0.5$, observable cuts corresponding to median halo masses between $M_\text{500c}=10^{14}$ and $10^{15}~\rm{M_{\odot}}$ give nearly unbiased median masses for all selection methods, but X-ray selection results in biased medians for higher masses. For cuts corresponding to median masses $<10^{14}$ at $z\leq0.5$ and for all masses at $z\geq1$, only Compton-Y selection yields nearly unbiased median masses. Importantly, even when the median mass is unbiased, the scatter is not because for each selection the sample is skewed towards lower masses than a mass-selected sample. Each selection leads to a different bias in secondary quantities like cool-core fraction, temperature and gas fraction., Comment: 19 pages, 12 figures (Including the appendix). Accepted for publication in MNRAS. Minor changes w.r.t. version 1

We present several machine learning (ML) models developed to efficiently separate stars formed in-situ in Milky Way-type galaxies from those that were formed externally and later accreted. These models, which include examples from artificial neural networks, decision trees and dimensionality reduction techniques, are trained on a sample of disc-like, Milky Way-mass galaxies drawn from the ARTEMIS cosmological hydrodynamical zoom-in simulations. We find that the input parameters which provide an optimal performance for these models consist of a combination of stellar positions, kinematics, chemical abundances ([Fe/H] and [$\alpha$/Fe]) and photometric properties. Models from all categories perform similarly well, with area under the precision-recall curve (PR-AUC) scores of $\simeq 0.6$. Beyond a galactocentric radius of $5$~kpc, models retrieve $>90\%$ of accreted stars, with a sample purity close to $60\%$, however the purity can be increased by adjusting the classification threshold. For one model, we also include host galaxy-specific properties in the training, to account for the variability of accretion histories of the hosts, however this does not lead to an improvement in performance. The ML models can identify accreted stars even in regions heavily dominated by the in-situ component (e.g., in the disc), and perform well on an unseen suite of simulations (the Auriga simulations). The general applicability bodes well for application of such methods on observational data to identify accreted substructures in the Milky Way without the need to resort to selection cuts for minimising the contamination from in-situ stars.

Large-scale structure surveys have reported measurements of the density of matter, $\Omega_\mathrm{m}$, and the amplitude of clustering, $\sigma_8$, that are in tension with the values inferred from observations of the cosmic microwave background. While this may be a sign of new physics that slows the growth of structure at late times, strong astrophysical feedback processes could also be responsible. In this work, we argue that astrophysical processes are not independent of cosmology and that their coupling naturally leads to stronger baryonic feedback in cosmological models with suppressed structure formation or when combined with a mechanism that removes dark matter from halos. We illustrate this with two well-motivated extensions of the Standard Model known to suppress structure formation: massive neutrinos and decaying dark matter. Our results, based on the FLAMINGO suite of hydrodynamical simulations, show that the combined effect of baryonic and non-baryonic suppression mechanisms is greater than the sum of its parts, particularly for decaying dark matter. We also show that the dependence of baryonic feedback on cosmology can be modelled as a function of the ratio $f_\mathrm{b}/c^2_\mathrm{v}\sim f_\mathrm{b}/(\Omega_\mathrm{m}\sigma_8)^{1/4}$ of the universal baryon fraction, $f_\mathrm{b}$, to a velocity-based definition of halo concentration, $c^2_\mathrm{v}$, giving an accurate fitting formula for the baryonic suppression of the matter power spectrum. Although the combination of baryonic and non-baryonic suppression mechanisms can resolve the tension, the models with neutrinos and decaying dark matter are challenged by constraints on the expansion history., Comment: 18 pages, 14 figures, MNRAS, accepted for publication

We present the new ARTEMIS Emulator suite of high resolution (baryon mass of $2.23 \times 10^{4}$ $h^{-1}$M$_{\odot}$) zoom-in simulations of Milky Way mass systems. Here, three haloes from the original ARTEMIS sample have been rerun multiple times, systematically varying parameters for the stellar feedback model, the density threshold for star formation, the reionisation redshift and the assumed warm dark matter (WDM) particle mass (assuming a thermal relic). From these simulations emulators are trained for a wide range of statistics that allow for fast predictions at combinations of parameters not originally sampled, running in $\sim 1$ms (a factor of $\sim 10^{11}$ faster than the simulations). In this paper we explore the dependence of the central haloes' stellar mass on the varied parameters, finding the stellar feedback parameters to be the most important. When constraining the parameters to match the present-day stellar mass halo mass relation inferred from abundance matching we find that there is a strong degeneracy in the stellar feedback parameters, corresponding to a freedom in formation time of the stellar component for a fixed halo assembly history. We additionally explore the dependence of the satellite stellar mass function, where it is found that variations in stellar feedback, the reionisation redshift and the WDM mass all have a significant effect. The presented emulators are a powerful tool which allows for fundamentally new ways of analysing and interpreting cosmological hydrodynamic simulations. Crucially, allowing their free (subgrid) parameters to be varied and marginalised, leading to more robust constraints and predictions., Comment: 18 pages, 9 figures, accepted to MNRAS

We introduce GODMAX (Gas thermODynamics and Matter distribution using jAX), a novel code designed to calculate correlations between the cosmological matter distribution and various gas thermodynamic quantities. Utilizing the extensive ANTILLES suite of 200 hydrodynamical simulations with a diverse range of baryonic feedback strengths, we jointly fit the 3D profiles of total matter distribution, electron density, and electron pressure across various halo masses and redshifts. By accommodating significant variations in gas profiles expected due to baryonic feedback, solving exact hydrostatic equilibrium equation and offering flexible modeling of non-thermal pressure support, GODMAX has the capability to jointly fit all these profiles within the measurement uncertainties. This advancement enables, for the first time, robust joint analyses of multiple cosmic probes, including the kinetic and thermal Sunyaev-Zel'dovich effect, weak lensing, and X-ray observations. Furthermore, the model accurately captures correlations between the total matter power suppression due to baryonic feedback and local average thermodynamic quantities, such as the baryon fraction and integrated tSZ effect, in high-mass halos, aligning with observations from hydrodynamical simulations. Looking ahead, we forecast the expected constraints on cosmological and baryonic parameters from upcoming weak lensing catalogs from the LSST and tSZ maps from the Simons Observatory. This analysis underscores the importance of cross-correlations between weak lensing and tSZ in enhancing parameter constraints by resolving major systematic uncertainties due to baryonic physics. The GODMAX code leverages the JAX library, resulting in a fully differentiable halo model with native GPU compilation support., Comment: 17 pages, 9 figures

Weak gravitational lensing convergence peaks, the local maxima in weak lensing convergence maps, have been shown to contain valuable cosmological information complementary to commonly used two-point statistics. To exploit the full power of weak lensing for cosmology, we must model baryonic feedback processes because these reshape the matter distribution on non-linear and mildly non-linear scales. We study the impact of baryonic physics on the number density of weak lensing peaks using the FLAMINGO cosmological hydrodynamical simulation suite. We generate ray-traced full-sky convergence maps mimicking the characteristics of a Stage IV weak lensing survey. We compare the number densities of peaks in simulations that have been calibrated to reproduce the observed galaxy mass function and cluster gas fraction or to match a shifted version of these, and that use either thermally driven or jet AGN feedback. We show that the differences induced by realistic baryonic feedback prescriptions (typically $5 - 30\%$ for $\kappa = 0.1 - 0.4$) are smaller than those induced by reasonable variations in cosmological parameters ($20 - 60\%$ for $\kappa = 0.1 - 0.4$) but must be modeled carefully to obtain unbiased results. The reasons behind these differences can be understood by considering the impact of feedback on halo masses, or by considering the impact of different cosmological parameters on the halo mass function. Our analysis demonstrates that, for the range of models we investigated, the baryonic suppression is insensitive to changes in cosmology up to $\kappa \approx 0.4$ and that the higher $\kappa$ regime is dominated by Poisson noise and cosmic variance., Comment: 18 pages, 12 figures (including the appendix), Accepted for publication in MNRAS

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 ($>10^6$) and clusters ($>10^5$) 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 $R_{\rm 500c}$ 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., Comment: 21 pages, 14 figures (including appendix). Accepted for publication in MNRAS. Moved Section 3 to Appendix A, new Section 3 compares scaling relations for different physics models. Added lower mass bin for physics model comparison of thermodynamic profiles. Conclusions unchanged

Redshift space distortions are an important probe of the growth of large-scale structure and for constraining cosmological parameters in general. As galaxy redshift surveys approach percent level precision in their observations of the two point clustering statistics, it is timely to review what effects baryons and associated processes such as feedback may have on small-scale clustering in redshift space. Contrary to previous studies in the literature, we show using the large-volume BAHAMAS hydrodynamic simulations that the effect of baryons can be as much as 1% in the k ~ 0.1 h/Mpc range for the monopole and 5% for quadrupole, and that this could rise to as much as 10% at k~10 h/Mpc in both measurements. For the halo power spectra, this difference can be as much 3-4% in the monopole on scales of 0.05 < k < 0.3 h/Mpc for 10^{13} M_sun/h haloes. We find that these deviations can be mitigated to the sub-percent level in the both the monopole and quadrupole up to k ~ 0.3 h/Mpc if the baryon corrected halo masses are used to calculate the redshift space power spectra. Finally, we use the cosmo-OWLS simulation suite to explore the changes in the redshift space power spectra with different feedback prescriptions, finding that there is a maximum of 15-20% difference between the redshift space monopole and quadrupole with and without baryons at k ~1-2 h/Mpc within these models., Comment: 16 pages, 11 figues, replaced with accepted version, fixed typos

Cosmology is poised to measure the neutrino mass sum $M_\nu$ 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 massive neutrino power spectrum, $\Delta^2_\nu(k)$, agrees with its companion N-body simulation at the $10\%-15\%$ level for $k \leq 1~h/$Mpc. Building upon the Mira-Titan IV emulator for the cold matter, we use FlowsForTheMasses to construct an emulator for $\Delta^2_\nu(k)$ covering a large range of cosmological parameters and neutrino fractions $\Omega_{\nu,0} h^2 \leq 0.01$, which corresponds to $M_\nu \leq 0.93$~eV. Consistent with FlowsForTheMasses at the $3.5\%$ 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_\nu=0.15$~eV model to a wide range of N-body simulation methods, we find agreement to $3\%$ for $k \leq 3 k_\mathrm{FS} = 0.17~h/$Mpc and to $19\%$ for $k \leq 0.4~h/$Mpc. We find that the enhancement factor, the ratio of $\Delta^2_\nu(k)$ to its linear-response equivalent, is most strongly correlated with $\Omega_{\nu,0} h^2$, and also with the clustering amplitude $\sigma_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_\nu$-sensitivity of the small-scale neutrino density., Comment: 17 pages, 14 figures, 3 tables. Emulator code available at: https://github.com/upadhye/Cosmic-Enu

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 '$S_8$ 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, 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 $S_8$ tension. Consistent with recent studies, we find the CMB lensing power spectrum is in excellent agreement with the standard model, whilst 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., Comment: 26 pages, 12 figures, MNRAS, accepted with minor revisions

Weak lensing of the cosmic microwave background is rapidly emerging as a powerful probe of neutrinos, dark energy, and new physics. We present a fast computation of the non-linear CMB lensing power spectrum which combines non-linear perturbation theory at early times with power spectrum emulation using cosmological simulations at late times. Comparing our calculation with lightcones from the FLAMINGO 5.6 Gpc cube dark-matter-only simulation, we confirm its accuracy to 1% (2%) up to multipoles L = 3000 (L = 5000) for a nuLambdaCDM cosmology consistent with current data. Clustering suppression due to small-scale baryonic phenomena such as feedback from active galactic nuclei can reduce the lensing power by of order 10%. To our perturbation theory and emulator-based calculation we add SP(k), a new fitting function for this suppression, and confirm its accuracy compared to the FLAMINGO hydrodynamic simulations to 4% at L = 5000, with similar accuracy for massive neutrino models. We further demonstrate that scale-dependent suppression due to neutrinos and baryons approximately factorize, implying that a careful treatment of baryonic feedback can limit biasing neutrino mass constraints., Comment: 14 pages, 21 figures, 1 table

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). Accepted for publication in MNRAS. Minor changes w.r.t. version 1. 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: 44 pages, 23 figures. Accepted for publication in MNRAS. V3 includes changes made in published version: jet simulations were redone to fix a bug, but the differences are nearly invisible. For visualizations, see the FLAMINGO website at https://flamingo.strw.leidenuniv.nl/

Upcoming large-scale structure surveys will measure the matter power spectrum to approximately percent 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 $\texttt{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 percent 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 \sim 10^{13-14}$ M$_\odot$). We present a simple Python implementation of our model, available at $\href{https://github.com/jemme07/pyspk}{\mathrm{https{:}//github.com/jemme07/pyspk}}$, which can be used to incorporate baryon effects in standard gravity-only predictions, allowing for marginalisation over baryon physics within cosmological pipelines., Comment: Refereed version accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS). 16 pages, 11 figures

Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ~ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest-mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses (Mstar <~ 1E10 Msun), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help to shed light on the low-mass tension., Comment: 20 pages, 10 figures, resubmitted to MNRAS following minor revisions

The radial acceleration relation (RAR) represents a tight empirical relation between the inferred total and baryonic centripetal accelerations, $g_{\rm{tot}}=GM_{\rm{tot}}(

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 broadband filters from ultraviolet to sub-millimetre wavelengths and for 18 different viewing angles. We compare the simulated ARTEMIS galaxies to observed galaxies in the DustPedia database 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 multi-wavelength 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., Comment: Accepted for publication by MNRAS

The power spectrum has been a workhorse for cosmological studies of large-scale structure. However, the present-day matter distribution is highly non-Gaussian and significant cosmological information is also contained in higher-order correlation functions. Meanwhile, baryon physics (particularly AGN feedback) has previously been shown to strongly affect the two-point statistics but there has been limited exploration of its effects on higher-order functions to date. Here we use the BAHAMAS suite of cosmological hydrodynamical simulations to explore the effects of baryon physics and massive neutrinos on the halo bispectrum. In contrast to matter clustering which is suppressed by baryon physics, we find that the halo clustering is typically enhanced. The strength of the effect and the scale over which it extends depends on how haloes are selected. On small scales (k > 1 $h$ Mpc$^{-1}$, dominated by satellites of groups/clusters), we find that the bispectrum is highly sensitive to the efficiency of star formation and feedback, making it an excellent testing ground for galaxy formation models. We show that the effects of feedback and the effects of massive neutrinos are largely separable (independent of each other) and that massive neutrinos strongly suppress the halo bispectrum on virtually all scales up to the free-streaming length (apart from the smallest scales, where baryon physics dominates). The strong sensitivity of the bispectrum to neutrinos on the largest scales and galaxy formation physics on the smallest scales bodes well for upcoming precision measurements from the next generation of wide-field surveys., Comment: 13 pages, 6+3 figures. Published in MNRAS

We present measurements of the intrinsic alignments (IAs) of the star-forming gas of galaxies in the EAGLE simulations. Radio continuum imaging of this gas enables cosmic shear measurements complementary to optical surveys. We measure the orientation of star-forming gas with respect to the direction to, and orientation of, neighbouring galaxies. Star-forming gas exhibits a preferentially radial orientation-direction alignment that is a decreasing function of galaxy pair separation, but remains significant to $\gtrsim 1$ Mpc at $z=0$. The alignment is qualitatively similar to that exhibited by the stars, but is weaker at fixed separation. Pairs of galaxies hosted by more massive subhaloes exhibit stronger alignment at fixed separation, but the strong alignment of close pairs is dominated by ${\sim}L^\star$ galaxies and their satellites. At fixed comoving separation, the radial alignment is stronger at higher redshift. The orientation-orientation alignment is consistent with random at all separations, despite subhaloes exhibiting preferential parallel minor axis alignment. The weaker IA of star-forming gas than for stars stems from the former's tendency to be less well aligned with the dark matter structure of galaxies than the latter, and implies that the systematic uncertainty due to IA may be less severe in radio continuum weak lensing surveys than in optical counterparts. Alignment models equating the orientation of star-forming gas discs to that of stellar discs or the DM structure of host subhaloes will therefore overestimate the impact of IAs on radio continuum cosmic shear measurements., Comment: 18 pages, 13 figures. Paper submitted to MNRAS