Search

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

The effect of baryon physics associated with galaxy formation onto the large-scale matter distribution of the Universe is a key uncertainty in the theoretical modelling required for the interpretation of Stage IV cosmology surveys. We use the FLAMINGO suite of simulations to study the baryon response due to galaxy formation of the total matter power spectrum. We find that it is only well converged for simulation volumes in excess of $(200~Mpc)^3$. We report results for simulations of varying feedback intensity, which either match the X-ray inferred gas fractions in clusters and the z=0 stellar mass function, or shifted versions of the data, as well as for different implementations of AGN feedback. We package our results in the form of a Gaussian process emulator which can rapidly reproduce all the simulations' predictions to better than 1% up to the comoving wavenumber $k = 10~h/Mpc$ and up to z=2 for all the feedback models present in the FLAMINGO suite. We find that the response becomes stronger, the range of scales affected increases, and the position of the minimum of the response moves to smaller scales as the redshift decreases. We find that lower gas fractions in groups and clusters lead to a stronger response and that the use of collimated jets instead of thermally driven winds for AGN feedback enhances the effect. Lowering the stellar masses at fixed cluster gas fractions also increases the magnitude of the response. We find only a small (1% at $k<10~h/Mpc$) dependence of our results on the background cosmology, but a wider range of cosmology variations will be needed to confirm this result. The response we obtain for our strongest feedback models is compatible with some of the recent analyses combining weak lensing with external data. Such a response is, however, in strong tension with the X-ray inferred gas fractions in clusters used to calibrate the FLAMINGO model., Comment: 14 pages, 13 figures, submitted to MNRAS, emulator will be released publicly with the final version of the paper, comments welcome

The masses of galaxy clusters are commonly measured from X-ray observations under the assumption of hydrostatic equilibrium (HSE). This technique is known to underestimate the true mass systematically. The fiducial FLAMINGO cosmological hydrodynamical simulation predicts the median hydrostatic mass bias to increase from $b_\text{HSE} \equiv (M_\text{HSE,500c}-M_\text{500c})/M_\text{500c} \approx -0.1$ to -0.2 when the true mass increases from group to cluster mass scales. However, the bias is nearly independent of the hydrostatic mass. The scatter at fixed true mass is minimum for $M_\text{500c}\sim 10^{14}~\text{M}_\odot$, where $\sigma(b_\text{HSE})\approx 0.1$, but increases rapidly towards lower and higher masses. At a fixed true mass, the hydrostatic masses increase (decrease) with redshift on group (cluster) scales, and the scatter increases. The bias is insensitive to the choice of analytic functions assumed to represent the density and temperature profiles, but it is sensitive to the goodness of fit, with poorer fits corresponding to a stronger median bias and a larger scatter. The bias is also sensitive to the strength of stellar and AGN feedback. Models predicting lower gas fractions yield more (less) biased masses for groups (clusters). The scatter in the bias at fixed true mass is due to differences in the pressure gradients rather than in the temperature at $R_\text{500c}$. The total kinetic energies within $r_\text{500c}$ in low- and high-mass clusters are sub- and super-virial, respectively, though all become sub-virial when external pressure is accounted for. Analyses of the terms in the virial and Euler equations suggest that non-thermal motions, including rotation, account for most of the hydrostatic mass bias. However, we find that the mass bias estimated from X-ray luminosity weighted profiles strongly overestimates the deviations from hydrostatic equilibrium., Comment: 19 pages, 14 figures (including appendix). Submitted to MNRAS

Galaxy cluster counts have historically been important for the measurement of cosmological parameters and upcoming surveys will greatly reduce the statistical errors. To exploit the potential of current and future cluster surveys, theoretical uncertainties on the predicted abundance must be smaller than the statistical errors. Models used to predict cluster counts typically combine a model for the dark matter only (DMO) halo mass function (HMF) with an observable - mass relation that is assumed to be a power-law with lognormal scatter. We use the FLAMINGO suite of cosmological hydrodynamical simulations to quantify the biases in the cluster counts and cosmological parameters resulting from the different ingredients of conventional models. For the observable mass proxy we focus on the Compton-Y parameter quantifying the thermal Sunyaev-Zel'dovich effect, which is expected to result in cluster samples that are relatively close to mass-selected samples. We construct three mock surveys based on existing (Planck and SPT) and upcoming (Simons Observatory) surveys. We ignore measurement uncertainties and compare the biases in the counts and inferred cosmological parameters to each survey's Poisson errors. We find that widely used models for the DMO HMF differ significantly from each other and from the DMO version of FLAMINGO, leading to significant biases for all three surveys. For upcoming surveys, dramatic improvements are needed for all additional model ingredients, i.e. the functional forms of the fits to the observable-mass scaling relation and the associated scatter, the priors on the scaling relation and the prior on baryonic effects associated with feedback processes on the HMF., Comment: 17 pages, 8 figures. Submitted to MNRAS

Context. Mock galaxy catalogues are essential for correctly interpreting current and future generations of galaxy surveys. Despite their significance in galaxy formation and cosmology, little to no work has been done to validate the predictions of these mocks for high-order clustering statistics. Aims. We compare the predicting power of the latest generation of empirical models used in the creation of mock galaxy catalogues: a 13-parameter Halo Occupation Distribution (HOD) and an extension of the SubHalo Abundance Matching technique (SHAMe). Methods. We build GalaxyEmu-Planck, an emulator that makes precise predictions for the two-point correlation function, galaxy-galaxy lensing (restricted to distances greater than 1 $h^{-1} {\rm Mpc}$ to avoid baryonic effects), and other high-order statistics resulting from the evaluation of SHAMe and HOD models. Results. We evaluate the precision of GalaxyEmu-Planck using two galaxy samples extracted from the FLAMINGO hydrodynamical simulation that mimic the properties of DESI-BGS and BOSS galaxies, finding that the emulator reproduces all the predicted statistics precisely. The HOD showed comparable performance when fitting galaxy clustering and galaxy-galaxy lensing. In contrast, the SHAMe model showed better predictions for higher-order statistics, especially regarding the galaxy assembly bias. We also tested the performance of the models after removing some of their extensions, finding that we can withdraw two of the HOD parameters without a loss of performance. Conclusions. The results of this paper validate the current generation of empirical models as a way to reproduce galaxy clustering, galaxy-galaxy lensing and other high-order statistics. The excellent performance of the SHAMe model with a small number of free parameters suggests that it is a valid method to extract cosmological constraints from galaxy clustering., Comment: 22 pages, 12 figures, submitted to A&A. Comments welcomed!

This study concerns an inventory of the gravitational force and tidal field induced by filaments, walls, cluster nodes and voids on Megaparsec scales and how they assemble and shape the Cosmic Web. The study is based on a N$_{\rm Part}=512^3$ $\Lambda$CDM dark matter only N-body simulation in a (300$h^{-1}$~Mpc)$^3$ box at $z=0$. We invoke the density field NEXUS+ multiscale morphological procedure to assign the appropriate morphological feature to each location. We then determine the contribution by each of the cosmic web components to the local gravitational and tidal forces. We find that filaments are, by far, the dominant dynamical component in the interior of filaments, in the majority of underdense void regions and in all wall regions. The gravitational influence of cluster nodes is limited, and they are only dominant in their immediate vicinity. The force field induced by voids is marked by divergent outflowing patterns, yielding the impression of a segmented volume in which voids push matter towards their boundaries. Voids manifest themselves strongly in the tidal field as a cellular tapestry that is closely linked to the multiscale cosmic web. However, even within the interior of voids, the dynamical influence of the surrounding filaments is stronger than the outward push by voids. Therefore, the dynamics of voids cannot be understood without taking into account the influence of the environment. We conclude that filaments constitute the overpowering gravitational agent of the cosmic web, while voids are responsible for the cosmic web's spatial organisation and hence of its spatial connectivity., Comment: 33 pages, 21 figures (Including the appendix). Submitted to MNRAS

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

The resonant conversion, within the inter-galactic medium, of regular photons into dark photons amplifies the anisotropy observed in the CMB, thereby imposing stringent constraints on the existence of light dark photons. In this study, we investigate the impact of light dark photons, with masses in the range $3\times 10^{-15} ~\rm{eV} < m_{A'} < 3\times 10^{-12}~\rm{eV}$ on the power spectrum of temperature anisotropies within the cosmic microwave background (CMB) radiation utilizing the state-of-the-art large-volume FLAMINGO cosmological simulations. Our results show that using full Planck data, one can expect the existing constraints on the dark photon mixing parameter in this mass range to improve by an order of magnitude., Comment: 16 pages, 8 figures

The relativistic Sunyaev-Zel'dovich (SZ) effect can be used to measure intracluster gas temperatures independently of X-ray spectroscopy. Here, we use the large-volume FLAMINGO simulation suite to determine whether SZ $y$-weighted temperatures lead to more accurate hydrostatic mass estimates in massive ($M_{\rm 500c} > 7.5\times 10^{14}\,{\rm M}_{\odot}$) clusters than when using X-ray spectroscopic-like temperatures. We find this to be the case, on average. The median bias in the SZ mass at redshift zero is $\left< b \right> \equiv 1-\left< M_{\rm 500c,hse}/M_{\rm 500c,true} \right> = -0.05 \pm 0.01$, over 4 times smaller in magnitude than the X-ray spectroscopic-like case, $\left< b \right> = 0.22 \pm 0.01$. However, the scatter in the SZ bias, $\sigma_{b} \approx 0.2$, is around 40 per cent larger than for the X-ray case. We show that this difference is strongly affected by clusters with large pressure fluctuations, as expected from shocks in ongoing mergers. Selecting the clusters with the best-fitting generalized NFW pressure profiles, the median SZ bias almost vanishes, $\left< b \right> = -0.009 \pm 0.005$, and the scatter is halved to $\sigma_{b} \approx 0.1$. We study the origin of the SZ/X-ray difference and find that, at $R_{\rm 500c}$ and in the outskirts, SZ weighted gas better reflects the hot, hydrostatic atmosphere than the X-ray weighted gas. The SZ/X-ray temperature ratio increases with radius, a result we find to be insensitive to variations in baryonic physics, cosmology and numerical resolution., Comment: 19 pages, 22 figures, submitted to MNRAS. Comments welcome

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: 16 pages, 9 figures, submitted to MNRAS

Motivated by the recent JWST discovery of galaxy overdensities during the Epoch of Reionzation, we examine the physical properties of high-$z$ protoclusters and their evolution using the FLAMINGO simulation suite. We investigate the impact of the apertures used to define protoclusters, because the heterogeneous apertures used in the literature have limited our understanding of the population. Our results are insensitive to the uncertainties of the subgrid models at a given resolution, whereas further investigation into the dependence on numerical resolution is needed. When considering galaxies more massive than $M_\ast\,{\simeq}\,10^8\,{\rm M_\odot}$, the FLAMINGO simulations predict a dominant contribution from progenitors similar to those of the Coma cluster to the cosmic star-formation rate density during the reionization epoch. Our results indicate the onset of suppression of star formation in the protocluster environments as early as $z\,{\simeq}\,5$. The galaxy number density profiles are similar to NFW at $z\,{\lesssim}\,1$ while showing a steeper slope at earlier times before the formation of the core. Different from most previous simulations, the predicted star-formation history for individual protoclusters is in good agreement with observations. We demonstrate that, depending on the aperture, the integrated physical properties including the total (dark matter and baryonic) mass can be biased by a factor of 2 to 5 at $z\,{=}\,5.5$--$7$, and by an order of magnitude at $z\,{\lesssim}\,4$. This correction suffices to remove the ${\simeq}\,3\,\sigma$ tensions with the number density of structures found in recent JWST observations., Comment: 22 pages, 15 figures, accepted for publication in MNRAS

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

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., Comment: Accepted by MNRAS

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/

Active galactic nucleus (AGN) feedback from accreting supermassive black holes (SMBHs) is an essential ingredient of galaxy formation simulations. The orbital evolution of SMBHs is affected by dynamical friction that cannot be predicted self-consistently by contemporary simulations of galaxy formation in representative volumes. Instead, such simulations typically use a simple "repositioning" of SMBHs, but the effects of this approach on SMBH and galaxy properties have not yet been investigated systematically. Based on a suite of smoothed particle hydrodynamics simulations with the SWIFT code and a Bondi-Hoyle-Lyttleton subgrid gas accretion model, we investigate the impact of repositioning on SMBH growth and on other baryonic components through AGN feedback. Across at least a factor ~1000 in mass resolution, SMBH repositioning (or an equivalent approach) is a necessary prerequisite for AGN feedback; without it, black hole growth is negligible. Limiting the effective repositioning speed to $\lesssim$ 10 km/s delays the onset of AGN feedback and severely limits its impact on stellar mass growth in the centre of massive galaxies. Repositioning has three direct physical consequences. It promotes SMBH mergers and thus accelerates their initial growth. In addition, it raises the peak density of the ambient gas and reduces the SMBH velocity relative to it, giving a combined boost to the accretion rate that can reach many orders of magnitude. Our results suggest that a more sophisticated and/or better calibrated treatment of SMBH repositioning is a critical step towards more predictive galaxy formation simulations., Comment: 17 pages, 11 figures (including appendix); accepted by MNRAS. Minor updates to match accepted version