Showing 1–50 of 156 results for author: Schaller, M

Validating the clustering predictions of empirical models with the FLAMINGO simulations

Authors: Sergio Contreras, Raul E. Angulo, Jonás Chaves-Montero, Roi Kugel, Matthieu Schaller, Joop Schaye

Abstract: 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 creat… ▽ More 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. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: 22 pages, 12 figures, submitted to A&A. Comments welcomed!

REMIX SPH -- improving mixing in smoothed particle hydrodynamics simulations using a generalised, material-independent approach

Authors: Thomas D. Sandnes, Vincent R. Eke, Jacob A. Kegerreis, Richard J. Massey, Sergio Ruiz-Bonilla, Matthieu Schaller, Luis F. A. Teodoro

Abstract: We present REMIX, a smoothed particle hydrodynamics (SPH) scheme designed to alleviate effects that typically suppress mixing and instability growth at density discontinuities in SPH simulations. We approach this problem by directly targeting sources of kernel smoothing error and discretisation error, resulting in a generalised, material-independent formulation that improves the treatment both of… ▽ More We present REMIX, a smoothed particle hydrodynamics (SPH) scheme designed to alleviate effects that typically suppress mixing and instability growth at density discontinuities in SPH simulations. We approach this problem by directly targeting sources of kernel smoothing error and discretisation error, resulting in a generalised, material-independent formulation that improves the treatment both of discontinuities within a single material, for example in an ideal gas, and of interfaces between dissimilar materials. This approach also leads to improvements in capturing hydrodynamic behaviour unrelated to mixing, such as in shocks. We demonstrate marked improvements in three-dimensional test scenarios, focusing on more challenging cases with particles of equal mass across the simulation. This validates our methods for use-cases relevant across applications spanning astrophysics and engineering, where particles are free to evolve over a large range of density scales, or where emergent and evolving density discontinuities cannot easily be corrected by choosing bespoke particle masses in the initial conditions. We achieve these improvements while maintaining sharp discontinuities; without introducing additional equation of state dependence in, for example, particle volume elements; and without contrived or targeted corrections. Our methods build upon a fully compressible and thermodynamically consistent core-SPH construction, retaining Galilean invariance as well as conservation of mass, momentum, and energy. REMIX is integrated in the open-source, state-of-the-art \swift code and is designed with computational efficiency in mind, which means that its improved hydrodynamic treatment can be used for high-resolution simulations without significant cost to run-speed. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: 51 pages, 33 figures, submitted to Journal of Computational Physics

The FLAMINGO Project: A comparison of galaxy cluster samples selected on mass, X-ray luminosity, Compton-Y parameter, or galaxy richness

Authors: Roi Kugel, Joop Schaye, Matthieu Schaller, Ian G. McCarthy, Joey Braspenning, John C. Helly, Victor J. Forouhar Moreno, Robert J. McGibbon

Abstract: 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$ F… ▽ More 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. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 18 pages, 10 figures (Including the appendix). Submitted to MNRAS

Dark photon constraints from CMB temperature anisotropies

Authors: Andres Aramburo-Garcia, Kyrylo Bondarenko, Alexey Boyarsky, Pavlo Kashko, Josef Pradler, Anastasia Sokolenko, Roi Kugel, Matthieu Schaller, Joop Schaye

Abstract: 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 t… ▽ More 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. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

Report number: UWThPh 2024-9, FERMILAB-PUB-24-0219-T, NORDITA 2024-013

Relativistic SZ temperatures and hydrostatic mass bias for massive clusters in the FLAMINGO simulations

Authors: Scott T. Kay, Joey Braspenning, Jens Chluba, John C. Helly, Roi Kugel, Matthieu Schaller, Joop Schaye

Abstract: 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 spectroscop… ▽ More 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, $σ_{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 $σ_{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. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 19 pages, 22 figures, submitted to MNRAS. Comments welcome

The FLAMINGO project: the coupling between baryonic feedback and cosmology in light of the $S_8$ tension

Authors: Willem Elbers, Carlos S. Frenk, Adrian Jenkins, Baojiu Li, John C. Helly, Roi Kugel, Matthieu Schaller, Joop Schaye, Joey Braspenning, Juliana Kwan, Ian G. McCarthy, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke, Silvia Pascoli

Abstract: Large-scale structure surveys have reported measurements of the density of matter, $Ω_\mathrm{m}$, and the amplitude of clustering, $σ_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 th… ▽ More Large-scale structure surveys have reported measurements of the density of matter, $Ω_\mathrm{m}$, and the amplitude of clustering, $σ_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}/(Ω_\mathrm{m}σ_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. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 16 pages, 9 figures, submitted to MNRAS

A unified model for the clustering of quasars and galaxies at $z\approx6$

Authors: Elia Pizzati, Joseph F. Hennawi, Joop Schaye, Matthieu Schaller, Anna-Christina Eilers, Feige Wang, Carlos S. Frenk, Willem Elbers, John C. Helly, Ruari Mackenzie, Jorryt Matthee, Rongmon Bordoloi, Daichi Kashino, Rohan P. Naidu, Minghao Yue

Abstract: Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we… ▽ More Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx2-4$. We make use of a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx6$ quasars and galaxies (i.e., their clustering properties and luminosity functions), and infer key quantities such as their luminosity-halo mass relation, the mass function of their host halos, and their duty cycle/occupation fraction. Our key findings are: (i) quasars reside on average in $\approx10^{12.4}\,{\rm M}_\odot$ halos (corresponding to $\approx5σ$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\lesssim1\%$); (iii) galaxies (that are bright in [OIII]) live in much smaller halos ($\approx10^{10.9}\,{\rm M}_\odot$) and have a larger duty cycle (occupation fraction) of $\approx20\%$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx2-6$, and calls for new investigations of the role of quasar activity across cosmic time. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: Submitted to MNRAS; 15 pages, 4 figures

TangoSIDM Project: Is the Stellar Mass Tully-Fisher relation consistent with SIDM?

Authors: Camila Correa, Matthieu Schaller, Joop Schaye, Sylvia Ploeckinger, Josh Borrow, Yannick Bahe

Abstract: Self-interacting dark matter (SIDM) has the potential to significantly influence galaxy formation in comparison to the cold, collisionless dark matter paradigm (CDM), resulting in observable effects. This study aims to elucidate this influence and to demonstrate that the stellar mass Tully-Fisher relation imposes robust constraints on the parameter space of velocity-dependent SIDM models. We prese… ▽ More Self-interacting dark matter (SIDM) has the potential to significantly influence galaxy formation in comparison to the cold, collisionless dark matter paradigm (CDM), resulting in observable effects. This study aims to elucidate this influence and to demonstrate that the stellar mass Tully-Fisher relation imposes robust constraints on the parameter space of velocity-dependent SIDM models. We present a new set of cosmological hydrodynamical simulations that include the SIDM scheme from the TangoSIDM project and the SWIFT-EAGLE galaxy formation model. Two cosmological simulations suites were generated: one (Reference model) which yields good agreement with the observed $z=0$ galaxy stellar mass function, galaxy mass-size relation, and stellar-to-halo mass relation; and another (WeakStellarFB model) in which the stellar feedback is less efficient, particularly for Milky Way-like systems. Both galaxy formation models were simulated under four dark matter cosmologies: CDM, SIDM with two different velocity-dependent cross sections, and SIDM with a constant cross section. While SIDM does not modify global galaxy properties such as stellar masses and star formation rates, it does make the galaxies more extended. In Milky Way-like galaxies, where baryons dominate the central gravitational potential, SIDM thermalises, causing dark matter to accumulate in the central regions. This accumulation results in density profiles that are steeper than those produced in CDM from adiabatic contraction. The enhanced dark matter density in the central regions of galaxies causes a deviation in the slope of the Tully-Fisher relation, which significantly diverges from the observational data. In contrast, the Tully-Fisher relation derived from CDM models aligns well with observations. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 17 pages, 10 figures and 3 tables. Submitted to MNRAS. Comments welcome

EIGER VI. The Correlation Function, Host Halo Mass and Duty Cycle of Luminous Quasars at $z\gtrsim6$

Authors: Anna-Christina Eilers, Ruari Mackenzie, Elia Pizzati, Jorryt Matthee, Joseph F. Hennawi, Haowen Zhang, Rongmon Bordoloi, Daichi Kashino, Simon J. Lilly, Rohan P. Naidu, Robert A. Simcoe, Minghao Yue, Carlos S. Frenk, John C. Helly, Matthieu Schaller, Joop Schaye

Abstract: We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detecte… ▽ More We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detected [OIII]-emitting galaxies in the quasars' environments, ranging between a density enhancement of $δ>100$ within a $2$ cMpc radius - one of the largest proto-clusters during the Epoch of Reionization discovered to date - to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy auto-correlation function, we infer a correlation length of quasars at $\langle z\rangle=6.25$ of $r_0^{\rm QQ}=21.3^{+2.7}_{-2.6}~{\rm cMpc}\,h^{-1}$, while we obtain a correlation length of the [OIII]-emitting galaxies of $r_0^{\rm GG}=4.2\pm0.1~{\rm cMpc}\,h^{-1}$. By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars' host dark matter halos to be $\log_{10}(M_{\rm halo, min}/M_\odot)=12.30\pm0.14$ (and $\log_{10}(M_{\rm halo, min}^{\rm [OIII]}/M_\odot) = 10.72\pm0.03$ for the [OIII]-emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is $f_{\rm duty}\ll1$. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: Submitted to AAS Journals

The FLAMINGO simulation view of cluster progenitors observed in the epoch of reionization with JWST

Authors: Seunghwan Lim, Sandro Tacchella, Joop Schaye, Matthieu Schaller, Jakob M. Helton, Roi Kugel, Roberto Maiolino

Abstract: 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… ▽ More 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\,σ$ tensions with the number density of structures found in recent JWST observations. △ Less

Submitted 24 July, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

Comments: 22 pages, 15 figures, accepted for publication in MNRAS

The baryon cycle in modern cosmological hydrodynamical simulations

Authors: Ruby J. Wright, Rachel S. Somerville, Claudia del P. Lagos, Matthieu Schaller, Romeel Davé, Daniel Anglés-Alcázar, Shy Genel

Abstract: In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a like-for-like comparison between the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predict… ▽ More In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a like-for-like comparison between the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, M_{\odot}$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2-3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the ISM, recycle within the CGM (within $R_{\rm 200c}$). We find that AGN-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses up to $M_{\rm 200c}\approx10^{13.5}\, M_{\odot}$. In both TNG and EAGLE, AGN feedback can eject gas beyond $R_{\rm 200c}$ at this mass scale, but seldom beyond $2-3\times R_{\rm 200c}$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $R_{\rm 200c}$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform sub-grid feedback models in the next generation of simulations. △ Less

Submitted 9 July, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures. Accepted for publication in MNRAS 03/07/24

The FLAMINGO project: Baryonic impact on weak gravitational lensing convergence peak counts

Authors: Jeger C. Broxterman, Matthieu Schaller, Joop Schaye, Henk Hoekstra, Konrad Kuijken, John C. Helly, Roi Kugel, Joey Braspenning, Willem Elbers, Carlos S. Frenk, Juliana Kwan, Ian G. McCarthy, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke

Abstract: 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.… ▽ More 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 $κ= 0.1 - 0.4$) are smaller than those induced by reasonable variations in cosmological parameters ($20 - 60\%$ for $κ= 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 $κ\approx 0.4$ and that the higher $κ$ regime is dominated by Poisson noise and cosmic variance. △ Less

Submitted 5 March, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 18 pages, 12 figures (including the appendix), Accepted for publication in MNRAS

Journal ref: MNRAS, Volume 529, Issue 3, April 2024, Pages 2309-2326

The FLAMINGO Project: Galaxy clusters in comparison to X-ray observations

Authors: Joey Braspenning, Joop Schaye, Matthieu Schaller, Ian G. McCarthy, Scott T. Kay, John C. Helly, Roi Kugel, Willem Elbers, Carlos S. Frenk, Juliana Kwan, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke

Abstract: 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… ▽ More 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. △ Less

Submitted 7 June, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 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

Inferring the dark matter splashback radius from cluster gas and observable profiles in the FLAMINGO simulations

Authors: Imogen Towler, Scott T. Kay, Joop Schaye, Roi Kugel, Matthieu Schaller, Joey Braspenning, Willem Elbers, Carlos S. Frenk, Juliana Kwan, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke, Edoardo Altamura

Abstract: 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… ▽ More 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. △ Less

Submitted 8 March, 2024; v1 submitted 8 December, 2023; originally announced December 2023.

Comments: Accepted by MNRAS

Revisiting the Extreme Clustering of $z \approx 4$ Quasars with Large Volume Cosmological Simulations

Authors: Elia Pizzati, Joseph F. Hennawi, Joop Schaye, Matthieu Schaller

Abstract: Observations from wide-field quasar surveys indicate that the quasar auto-correlation length increases dramatically from $z\approx2.5$ to $z\approx4$. This large clustering amplitude at $z\approx4$ has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter halos residing in the most extreme environments at that redshift. In this work, we prese… ▽ More Observations from wide-field quasar surveys indicate that the quasar auto-correlation length increases dramatically from $z\approx2.5$ to $z\approx4$. This large clustering amplitude at $z\approx4$ has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter halos residing in the most extreme environments at that redshift. In this work, we present a model that simultaneously reproduces both the observed quasar auto-correlation and quasar luminosity functions. The spatial distribution of halos and their relative abundance are obtained via a novel method that computes the halo mass and halo cross-correlation functions by combining multiple large-volume dark-matter-only cosmological simulations with different box sizes and resolutions. Armed with these halo properties, our model exploits the conditional luminosity function framework to describe the stochastic relationship between quasar luminosity, $L$, and halo mass, $M$. Assuming a simple power-law relation $L\propto M^γ$ with log-normal scatter, $σ$, we are able to reproduce observations at $z\sim 4$ and find that: (a) the quasar luminosity-halo mass relation is highly non-linear ($γ\gtrsim2$), with very little scatter ($σ\lesssim 0.3$ dex); (b) luminous quasars ($\log_{10} L/{\rm erg s}^{-1} \gtrsim 46.5-47$) are hosted by halos with mass $\log_{10} M/{\rm M}_\odot\gtrsim 13-13.5$; and (c) the implied duty cycle for quasar activity approaches unity ($\varepsilon_{\rm DC}\approx10-60\%$). We also consider observations at $z\approx2.5$ and find that the quasar luminosity-halo mass relation evolves significantly with cosmic time, implying a rapid change in quasar host halo masses and duty cycles, which in turn suggests concurrent evolution in black hole scaling relations and/or accretion efficiency. △ Less

Submitted 1 February, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 19 pages, 6 figures; accepted for publication in MNRAS

Realistic simulated galaxies form [$α$/Fe]-[Fe/H] knees due to a sustained decline in their star formation rates

Authors: Andrew C. Mason, Robert A. Crain, Ricardo P. Schiavon, David H. Weinberg, Joel Pfeffer, Joop Schaye, Matthieu Schaller, Tom Theuns

Abstract: We examine the stellar [$α$/Fe]-[Fe/H] distribution of $\simeq1000$ present-day galaxies in a high-resolution EAGLE simulation. A slight majority of galaxies exhibit the canonical distribution, characterised by a sequence of low-metallicity stars with high [$α$/Fe] that transitions at a "knee" to a sequence of declining [$α$/Fe] with increasing metallicity. This population yields a knee metallicit… ▽ More We examine the stellar [$α$/Fe]-[Fe/H] distribution of $\simeq1000$ present-day galaxies in a high-resolution EAGLE simulation. A slight majority of galaxies exhibit the canonical distribution, characterised by a sequence of low-metallicity stars with high [$α$/Fe] that transitions at a "knee" to a sequence of declining [$α$/Fe] with increasing metallicity. This population yields a knee metallicity - galaxy-mass relation similar to that observed in Local Group galaxies, both in slope and scatter. However, many simulated galaxies lack a knee or exhibit more complicated distributions. Knees are found only in galaxies with star formation histories (SFHs) featuring a sustained decline from an early peak ($t\simeq7~{\rm Gyr}$), which enables enrichment by Type Ia supernovae to dominate that due to Type II supernovae (SN II), reducing [$α$/Fe] in the interstellar gas. The simulation thus indicates that, contrary to the common interpretation implied by analytic galactic chemical evolution (GCE) models, knee formation is not a consequence of the onset of enrichment by SN Ia. We use the SFH of a simulated galaxy exhibiting a knee as input to the VICE GCE model, finding it yields an $α$-rich plateau enriched only by SN II, but the plateau comprises little stellar mass and the galaxy forms few metal-poor ([Fe/H]$\lesssim - 1$) stars. This follows from the short, constant gas consumption timescale typically assumed by GCEs, which implies the presence of a readily-enriched, low-mass gas reservoir. When an initially longer, evolving consumption timescale is adopted, VICE reproduces the simulated galaxy's track through the [$α$/Fe]-[Fe/H] plane and its metallicity distribution function. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 17 pages, 12 figures, 1 table, submitted to MNRAS

Resolution criteria to avoid artificial clumping in Lagrangian hydrodynamic simulations with a multi-phase interstellar medium

Authors: Sylvia Ploeckinger, Folkert S. J. Nobels, Matthieu Schaller, Joop Schaye

Abstract: Large-scale cosmological galaxy formation simulations typically prevent gas in the interstellar medium (ISM) from cooling below $\approx 10^4$ K. This has been motivated by the inability to resolve the Jeans mass in molecular gas (>>$10^5\,\mathrm{M}_{\odot}$) which would result in undesired artificial clumping. We show that the classical Jeans criteria derived for Newtonian gravity are not applic… ▽ More Large-scale cosmological galaxy formation simulations typically prevent gas in the interstellar medium (ISM) from cooling below $\approx 10^4$ K. This has been motivated by the inability to resolve the Jeans mass in molecular gas (>>$10^5\,\mathrm{M}_{\odot}$) which would result in undesired artificial clumping. We show that the classical Jeans criteria derived for Newtonian gravity are not applicable in the simulated ISM if the spacing of resolution elements representing the dense ISM is below the gravitational force softening length and gravity is therefore softened and not Newtonian. We re-derive the Jeans criteria for softened gravity in Lagrangian codes and use them to analyse gravitational instabilities at and below the hydrodynamical resolution limit for simulations with adaptive and constant gravitational softening lengths. In addition, we define criteria for which a numerical runaway collapse of dense gas clumps can occur caused by over-smoothing of the hydrodynamical properties relative to the gravitational force resolution. This effect is illustrated using simulations of isolated disk galaxies with the smoothed particle hydrodynamics code Swift. We also demonstrate how to avoid the formation of artificial clumps in gas and stars by adjusting the gravitational and hydrodynamical force resolutions. △ Less

Submitted 15 February, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 24 pages, 15 figures, accepted for publication in MNRAS, smaller updates to match published version

Journal ref: MNRAS, Volume 528, Issue 2, pp.2930-2951, 2024

On the anisotropic distribution of clusters in the local Universe

Authors: Matthieu Schaller

Abstract: In his 2021 lecture to the Canadian Association of Physicists Congress, P.J.E. Peebles pointed out that the brightest extra-galactic radio sources tend to be aligned with the plane of the de Vaucouleur Local Supercluster up to redshifts of $z=0.02$ ($d_{\rm MW}\approx 85~\rm{Mpc}$). He then asked whether such an alignment of clusters is anomalous in the standard $Λ$CDM framework. In this letter, w… ▽ More In his 2021 lecture to the Canadian Association of Physicists Congress, P.J.E. Peebles pointed out that the brightest extra-galactic radio sources tend to be aligned with the plane of the de Vaucouleur Local Supercluster up to redshifts of $z=0.02$ ($d_{\rm MW}\approx 85~\rm{Mpc}$). He then asked whether such an alignment of clusters is anomalous in the standard $Λ$CDM framework. In this letter, we employ an alternative, absolute orientation agnostic, measure of the anisotropy based on the inertia tensor axis ratio of these brightest sources and use a large cosmological simulation from the FLAMINGO suite to measure how common such an alignment of structures is. We find that only 3.5% of randomly selected regions display an anisotropy of their clusters more extreme than the one found in the local Universe's radio data. This sets the region around the Milky Way as a $1.85σ$ outlier. Varying the selection parameters of the objects in the catalogue, we find that the clusters in the local Universe are never more than $2σ$ away from the simulations' prediction for the same selection. We thus conclude that the reported anisotropy, whilst note-worthy, is not in tension with the $Λ$CDM paradigm. △ Less

Submitted 19 December, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 5 pages, 4 figures, Accepted for publication in MNRAS

Journal ref: MNRAS Letter, Volume 529, Issue 1, March 2024, Pages L23-L27

Tests of subgrid models for star formation using simulations of isolated disk galaxies

Authors: Folkert S. J. Nobels, Joop Schaye, Matthieu Schaller, Sylvia Ploeckinger, Evgenii Chaikin, Alexander J. Richings

Abstract: We use smoothed-particle hydrodynamics simulations of isolated Milky Way-mass disk galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Ke… ▽ More We use smoothed-particle hydrodynamics simulations of isolated Milky Way-mass disk galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt-Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially-resolved and azimuthally-averaged observed KS relations for neutral, atomic and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially-resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity. △ Less

Submitted 24 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS, 23 pages, 20 figures

The FLAMINGO project: revisiting the $S_8$ tension and the role of baryonic physics

Authors: Ian G. McCarthy, Jaime Salcido, Joop Schaye, Juliana Kwan, Willem Elbers, Roi Kugel, Matthieu Schaller, John C. Helly, Joey Braspenning, Carlos S. Frenk, Marcel P. van Daalen, Bert Vandenbroucke, Jonah T. Conley, Andreea S. Font, Amol Upadhye

Abstract: 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… ▽ More 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. △ Less

Submitted 9 October, 2023; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 26 pages, 12 figures, MNRAS, accepted with minor revisions

Non-linear CMB lensing with neutrinos and baryons: FLAMINGO simulations vs. fast approximations

Authors: Amol Upadhye, Juliana Kwan, Ian G. McCarthy, Jaime Salcido, John C. Helly, Roi Kugel, Matthieu Schaller, Joop Schaye, Joey Braspenning, Willem Elbers, Carlos S. Frenk, Marcel P. van Daalen, Bert Vandenbroucke, Jeger C. Broxterman

Abstract: 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 FLAMI… ▽ More 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. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 14 pages, 21 figures, 1 table

Winds versus jets: a comparison between black hole feedback modes in simulations of idealized galaxy groups and clusters

Authors: Filip Huško, Cedric G. Lacey, Joop Schaye, Folkert S. J. Nobels, Matthieu Schaller

Abstract: Using the SWIFT simulation code we study different forms of active galactic nuclei (AGN) feedback in idealized galaxy groups and clusters. We first present a physically motivated model of black hole (BH) spin evolution and a numerical implementation of thermal isotropic feedback (representing the effects of energy-driven winds) and collimated kinetic jets that they launch at different accretion ra… ▽ More Using the SWIFT simulation code we study different forms of active galactic nuclei (AGN) feedback in idealized galaxy groups and clusters. We first present a physically motivated model of black hole (BH) spin evolution and a numerical implementation of thermal isotropic feedback (representing the effects of energy-driven winds) and collimated kinetic jets that they launch at different accretion rates. We find that kinetic jet feedback is more efficient at quenching star formation in the brightest cluster galaxies (BCGs) than thermal isotropic feedback, while simultaneously yielding cooler cores in the intracluster medium (ICM). A hybrid model with both types of AGN feedback yields moderate star formation rates, while having the coolest cores. We then consider a simplified implementation of AGN feedback by fixing the feedback efficiencies and the jet direction, finding that the same general conclusions hold. We vary the feedback energetics (the kick velocity and the heating temperature), the fixed efficiencies and the type of energy (kinetic versus thermal) in both the isotropic and the jet case. The isotropic case is largely insensitive to these variations. In particular, we highlight that kinetic isotropic feedback (used e.g. in IllustrisTNG) is similar in its effects to its thermal counterpart (used e.g. in EAGLE). On the other hand, jet feedback must be kinetic in order to be efficient at quenching. We also find that it is much more sensitive to the choice of energy per feedback event (the jet velocity), as well as the efficiency. The former indicates that jet velocities need to be carefully chosen in cosmological simulations, while the latter motivates the use of BH spin evolution models. △ Less

Submitted 15 November, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Accepted for publication in MNRAS

Hydrodynamic simulations of the Disk of Gas Around Supermassive black holes (HDGAS) -I; Molecular Gas Dynamics

Authors: Mojtaba Raouf, Serena Viti, S. García-Burillo, Alexander J. Richings, Joop schaye, Ashley Bemis, Folkert S. J. Nobels, Matteo Guainazzi, Ko-Yun Huang, Matthieu Schaller, Violette Impellizzeri, Jon Holdship

Abstract: We present hydrodynamic simulations of the interstellar medium (ISM) within the circumnuclear disk (CND) of a typical AGN-dominated galaxy influenced by mechanical feedback from an active galactic nucleus(AGN). The simulations are coupled with the CHIMES non-equilibrium chemistry network to treat the radiative-cooling and AGN-heating. A focus is placed on the central 100 pc scale where AGN outflow… ▽ More We present hydrodynamic simulations of the interstellar medium (ISM) within the circumnuclear disk (CND) of a typical AGN-dominated galaxy influenced by mechanical feedback from an active galactic nucleus(AGN). The simulations are coupled with the CHIMES non-equilibrium chemistry network to treat the radiative-cooling and AGN-heating. A focus is placed on the central 100 pc scale where AGN outflows are coupled to the ISM and constrained by observational Seyfert-2 galaxies. AGN-feedback models are implemented with different wind-velocity and mass-loading factors. We post-process the simulation snapshots with a radiative-transfer code to obtain the molecular emission lines. We find that the inclusion of an AGN promotes the formation of CO in clumpy and dense regions surrounding supermassive-blackholes (SMBH). The CO(1-0) intensity maps ($<$6 Myr) in the CND seem to match well with observations of NGC 1068 with a best match for a model with 5000 $\rm km/s$ wind-velocity and a high mass-loading factor. We attempt to discern between competing explanations for the apparent counter-rotating gas disk in the NGC 1068 through an analysis of kinematic maps of the CO line emission. We suggest that mechanical AGN-feedback could explain the alignment-stability of position-angle across the different CND radii around the SMBH through momentum and energy loading of the wind. It is the wind-velocity that drives the disk out of alignment on a 100 pc scale for a long period of time. The position-velocity diagrams are in broad agreement with the predicted Keplerian rotation-curve in the model without-AGN, but the AGN models exhibit a larger degree of scatter, in better agreement with NGC 1068 observations. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

FLAMINGO: Calibrating large cosmological hydrodynamical simulations with machine learning

Authors: Roi Kugel, Joop Schaye, Matthieu Schaller, John C. Helly, Joey Braspenning, Willem Elbers, Carlos S. Frenk, Ian G. McCarthy, Juliana Kwan, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke, Yannick M. Bahé, Josh Borrow, Evgenii Chaikin, Filip Huško, Adrian Jenkins, Cedric G. Lacey, Folkert S. J. Nobels, Ian Vernon

Abstract: 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 paramete… ▽ More 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. △ Less

Submitted 23 October, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 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/

The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys

Authors: Joop Schaye, Roi Kugel, Matthieu Schaller, John C. Helly, Joey Braspenning, Willem Elbers, Ian G. McCarthy, Marcel P. van Daalen, Bert Vandenbroucke, Carlos S. Frenk, Juliana Kwan, Jaime Salcido, Yannick M. Bahé, Josh Borrow, Evgenii Chaikin, Oliver Hahn, Filip Huško, Adrian Jenkins, Cedric G. Lacey, Folkert S. J. Nobels

Abstract: 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 us… ▽ More 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. △ Less

Submitted 20 October, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 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/

Journal ref: 2023, MNRAS, 526, 4978

SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications

Authors: Matthieu Schaller, Josh Borrow, Peter W. Draper, Mladen Ivkovic, Stuart McAlpine, Bert Vandenbroucke, Yannick Bahé, Evgenii Chaikin, Aidan B. G. Chalk, Tsang Keung Chan, Camila Correa, Marcel van Daalen, Willem Elbers, Pedro Gonnet, Loïc Hausammann, John Helly, Filip Huško, Jacob A. Kegerreis, Folkert S. J. Nobels, Sylvia Ploeckinger, Yves Revaz, William J. Roper, Sergio Ruiz-Bonilla, Thomas D. Sandnes, Yolan Uyttenhove , et al. (2 additional authors not shown)

Abstract: Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gra… ▽ More Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code SWIFT. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. SWIFT also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarise the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with SWIFT. △ Less

Submitted 29 March, 2024; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 43 pages, 20 figures, accepted for publication in MNRAS. Code, documentation, and examples available at www.swiftsim.com

Journal ref: MNRAS, Volume 530, Issue 2, May 2024, Pages 2378-2419

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

KURVS: The outer rotation curve shapes and dark matter fractions of $z \sim 1.5 $ star-forming galaxies

Authors: Annagrazia Puglisi, Ugnė Dudzevičiūtė, Mark Swinbank, Steven Gillman, Alfred L. Tiley, Richard G. Bower, Michele Cirasuolo, Luca Cortese, Karl Glazebrook, Chris Harrison, Edo Ibar, Juan Molina, Danail Obreschkow, Kyle A. Oman, Matthieu Schaller, Francesco Shankar, Ray M. Sharples

Abstract: We present first results from the KMOS Ultra-deep Rotation Velocity Survey (KURVS), aimed at studying the outer rotation curves shape and dark matter content of 22 star-forming galaxies at $z\sim1.5$. These galaxies represent `typical' star-forming discs at $z \sim 1.5$, being located within the star-forming main sequence and stellar mass-size relation with stellar masses $9.5\leqslant$log… ▽ More We present first results from the KMOS Ultra-deep Rotation Velocity Survey (KURVS), aimed at studying the outer rotation curves shape and dark matter content of 22 star-forming galaxies at $z\sim1.5$. These galaxies represent `typical' star-forming discs at $z \sim 1.5$, being located within the star-forming main sequence and stellar mass-size relation with stellar masses $9.5\leqslant$log$(M_{\star}/\mathrm{M_{\odot}})\leqslant11.5$. We extract individual rotation curves out to 4 times the effective radius, on average, or $\sim 10-15$ kpc. Most rotation curves are flat or rising between three- and six-disc scale radii. Only three objects with dispersion-dominated dynamics ($v_{\rm rot}/σ_0\sim0.2$) have declining outer rotation curves at more than 5$σ$ significance. After accounting for seeing and pressure support, the nine rotation-dominated discs with $v_{\rm rot}/σ_0\geqslant1.5$ have average dark matter fractions of $50 \pm 20\%$ at the effective radius, similar to local discs. Together with previous observations of star-forming galaxies at cosmic noon, our measurements suggest a trend of declining dark matter fraction with increasing stellar mass and stellar mass surface density at the effective radius. Simulated EAGLE galaxies are in quantitative agreement with observations up to log$(M_{\star}R_{\rm eff}^{-2}/\mathrm{M_{\odot}kpc^{-2}}) \sim 9.2$, and over-predict the dark matter fraction of galaxies with higher mass surface densities by a factor of $\sim 3$. We conclude that the dynamics of typical rotationally-supported discs at $z \sim 1.5$ is dominated by dark matter from effective radius scales, in broad agreement with cosmological models. The tension with observations at high stellar mass surface density suggests that the prescriptions for baryonic processes occurring in the most massive galaxies (such as bulge growth and quenching) need to be reassessed. △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: 23 pages, 9 figures. Resubmitted to MNRAS after addressing the referee's comments. Abstract slightly modified to compile with the arXiv formatting

The impact of stochastic modeling on the predictive power of galaxy formation simulations

Authors: Josh Borrow, Matthieu Schaller, Yannick M. Bahe, Joop Schaye, Aaron D. Ludlow, Sylvia Ploeckinger, Folkert S. J. Nobels, Edoardo Altamura

Abstract: All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar… ▽ More All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long ($>$ Gyr) timescales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the SWIFT code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses ($M_{\rm BH} > 10^8$ M$_\odot$), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25\%, even far into the well-resolved regime, driven by bursty physics (black hole feedback) and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results. △ Less

Submitted 22 September, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

A thermal-kinetic subgrid model for supernova feedback in simulations of galaxy formation

Authors: Evgenii Chaikin, Joop Schaye, Matthieu Schaller, Alejandro Benítez-Llambay, Folkert S. J. Nobels, Sylvia Ploeckinger

Abstract: We present a subgrid model for supernova feedback designed for cosmological simulations of galaxy formation that may include a cold interstellar medium (ISM). The model uses thermal and kinetic channels of energy injection, which are built upon the stochastic kinetic and thermal models for stellar feedback used in the OWLS and EAGLE simulations, respectively. In the thermal channel, the energy is… ▽ More We present a subgrid model for supernova feedback designed for cosmological simulations of galaxy formation that may include a cold interstellar medium (ISM). The model uses thermal and kinetic channels of energy injection, which are built upon the stochastic kinetic and thermal models for stellar feedback used in the OWLS and EAGLE simulations, respectively. In the thermal channel, the energy is distributed statistically isotropically and injected stochastically in large amounts per event, which minimizes spurious radiative energy losses. In the kinetic channel, we inject the energy in small portions by kicking gas particles in pairs in opposite directions. The implementation of kinetic feedback is designed to conserve energy, linear and angular momentum, and is statistically isotropic. To test the model, we run simulations of isolated Milky Way-mass and dwarf galaxies, in which the gas is allowed to cool down to 10 K. Using the thermal and kinetic channels together, we obtain smooth star formation histories and powerful galactic winds with realistic mass loading factors. Furthermore, the model produces spatially resolved star formation rates (SFRs) and velocity dispersions that are in agreement with observations. We vary the numerical resolution by several orders of magnitude and find excellent convergence of the global SFRs and wind mass loading. We show that large thermal-energy injections generate a hot phase of the ISM and modulate the star formation by ejecting gas from the disc, while the low-energy kicks increase the turbulent velocity dispersion in the neutral ISM, which in turn helps suppress star formation. △ Less

Submitted 27 May, 2023; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: 23 pages, 17 figures (including appendix); Accepted for publication in MNRAS

EAGLE-like simulation models do not solve the entropy core problem in groups and clusters of galaxies

Authors: Edoardo Altamura, Scott T. Kay, Richard G. Bower, Matthieu Schaller, Yannick M. Bahé, Joop Schaye, Josh Borrow, Imogen Towler

Abstract: Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to inve… ▽ More Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to investigate the sensitivity of the central entropy and the shape of the profiles to changes in the sub-grid model applied to a suite of zoom-in cosmological simulations of a group of mass $M_{500} = 8.8 \times 10^{12}~{\rm M}_\odot$ and a cluster of mass $2.9 \times 10^{14}~{\rm M}_\odot$. Using our reference model, calibrated to match the stellar mass function of field galaxies, we confirm that our simulated groups and clusters contain hot gas with too high entropy in their cores. Additional simulations run without artificial conduction, metal cooling or AGN feedback produce lower entropy levels but still fail to reproduce observed profiles. Conversely, the two objects run without supernova feedback show a significant entropy increase which can be attributed to excessive cooling and star formation. Varying the AGN heating temperature does not greatly affect the profile shape, but only the overall normalisation. Finally, we compared runs with four AGN heating schemes and obtained similar profiles, with the exception of bipolar AGN heating, which produces a higher and more uniform entropy distribution. Our study leaves open the question of whether the entropy core problem in simulations, and particularly the lack of power-law cool-core profiles, arise from incorrect physical assumptions, missing physical processes, or insufficient numerical resolution. △ Less

Submitted 9 February, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS on 2023 February 06. The code used in the analysis is publicly available on the corresponding author's GitHub repository: https://github.com/edoaltamura/entropy-core-problem

Continuous Simulation Data Stream: A dynamical timescale-dependent output scheme for simulations

Authors: Loic Hausammann, Pedro Gonnet, Matthieu Schaller

Abstract: Exa-scale simulations are on the horizon but almost no new design for the output has been proposed in recent years. In simulations using individual time steps, the traditional snapshots are over resolving particles/cells with large time steps and are under resolving the particles/cells with short time steps. Therefore, they are unable to follow fast events and use efficiently the storage space. Th… ▽ More Exa-scale simulations are on the horizon but almost no new design for the output has been proposed in recent years. In simulations using individual time steps, the traditional snapshots are over resolving particles/cells with large time steps and are under resolving the particles/cells with short time steps. Therefore, they are unable to follow fast events and use efficiently the storage space. The Continuous Simulation Data Stream (CSDS) is designed to decrease this space while providing an accurate state of the simulation at any time. It takes advantage of the individual time step to ensure the same relative accuracy for all the particles. The outputs consist of a single file representing the full evolution of the simulation. Within this file, the particles are written independently and at their own frequency. Through the interpolation of the records, the state of the simulation can be recovered at any point in time. In this paper, we show that the CSDS can reduce the storage space by 2.76x for the same accuracy than snapshots or increase the accuracy by 67.8x for the same storage space whilst retaining an acceptable reading speed for analysis. By using interpolation between records, the CSDS provides the state of the simulation, with a high accuracy, at any time. This should largely improve the analysis of fast events such as supernovae and simplify the construction of light-cone outputs. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Comments: Accepted for publication in "Astronomy and Computing"

A common origin for the Fundamental Plane of quiescent and star-forming galaxies in the EAGLE simulations

Authors: Anna de Graaff, Marijn Franx, Eric F. Bell, Rachel Bezanson, Matthieu Schaller, Joop Schaye, Arjen van der Wel

Abstract: We use the EAGLE cosmological simulations to perform a comprehensive and systematic analysis of the $z=0.1$ Fundamental Plane (FP), the tight relation between galaxy size, mass and velocity dispersion. We first measure the total mass and velocity dispersion (including both random and rotational motions) within the effective radius to show that simulated galaxies obey a total mass FP that is very c… ▽ More We use the EAGLE cosmological simulations to perform a comprehensive and systematic analysis of the $z=0.1$ Fundamental Plane (FP), the tight relation between galaxy size, mass and velocity dispersion. We first measure the total mass and velocity dispersion (including both random and rotational motions) within the effective radius to show that simulated galaxies obey a total mass FP that is very close to the virial relation ($<10\%$ deviation), indicating that the effects of non-homology are weak. When we instead use the stellar mass, we find a strong deviation from the virial plane, which is driven by variations in the dark matter content. The dark matter fraction is a smooth function of the size and stellar mass, and thereby sets the coefficients of the stellar mass FP without substantially increasing the scatter. Hence, both star-forming and quiescent galaxies obey the same FP, with equally low scatter (0.02 dex). We employ simulations with a variable stellar initial mass function (IMF) to show that IMF variations have a modest additional effect on this FP. Moreover, when we use luminosity-weighted mock observations of the size and spatially-integrated velocity dispersion, the inferred FP changes only slightly. However, the scatter increases significantly, due to the luminosity-weighting and line-of-sight projection of the velocity dispersions, and measurement uncertainties on the half-light radii. Importantly, we find significant differences between the simulated FP and observations, which likely reflects a systematic difference in the stellar mass distributions. Therefore, we suggest the stellar mass FP offers a simple test for cosmological simulations, requiring minimal post-processing of simulation data. △ Less

Submitted 8 November, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: accepted for publication in MNRAS; resolution of some figures reduced for arxiv version

TangoSIDM: Tantalizing models of Self-Interacting Dark Matter

Authors: Camila A. Correa, Matthieu Schaller, Sylvia Ploeckinger, Noemi Anau Montel, Christoph Weniger, Shinichiro Ando

Abstract: We introduce the TangoSIDM project, a suite of cosmological simulations of structure formation in a $Λ$-Self-Interacting Dark Matter (SIDM) universe. TangoSIDM explores the impact of large dark matter (DM) scattering cross sections over dwarf galaxy scales. Motivated by DM interactions that follow a Yukawa potential, the cross section per unit mass, $σ/m_χ$, assumes a velocity dependent form that… ▽ More We introduce the TangoSIDM project, a suite of cosmological simulations of structure formation in a $Λ$-Self-Interacting Dark Matter (SIDM) universe. TangoSIDM explores the impact of large dark matter (DM) scattering cross sections over dwarf galaxy scales. Motivated by DM interactions that follow a Yukawa potential, the cross section per unit mass, $σ/m_χ$, assumes a velocity dependent form that avoids violations of current constraints on large scales. We demonstrate that our implementation accurately models not only core formation in haloes, but also gravothermal core collapse. For central haloes in cosmological volumes, frequent DM particle collisions isotropise the particles orbit, making them largely spherical. We show that the velocity-dependent $σ/m_χ$ models produce a large diversity in the circular velocities of satellites haloes, with the spread in velocities increasing as the cross sections reach 20, 60 and 100 cm$^2$/g in $10^9~\rm{M}_{\odot}$ haloes. The large variation in the haloes internal structure is driven by DM particles interactions, causing in some haloes the formation of extended cores, whereas in others gravothermal core collapse. We conclude that the SIDM models from the Tango project offer a promising explanation for the diversity in the density and velocity profiles of observed dwarf galaxies. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 19 pages, 14 figures, submitted to MNRAS

Spin-driven jet feedback in idealised simulations of galaxy groups and clusters

Authors: Filip Huško, Cedric G. Lacey, Joop Schaye, Matthieu Schaller, Folkert S. J. Nobels

Abstract: We implement a black hole spin evolution and jet feedback model into SWIFT, a smoothed particle hydrodynamics code. The jet power is determined self-consistently assuming Bondi accretion, using a realistic, spin-dependant efficiency. The jets are launched along the spin axis of the black hole, resulting in natural reorientation and precession. We apply the model to idealised simulations of galaxy… ▽ More We implement a black hole spin evolution and jet feedback model into SWIFT, a smoothed particle hydrodynamics code. The jet power is determined self-consistently assuming Bondi accretion, using a realistic, spin-dependant efficiency. The jets are launched along the spin axis of the black hole, resulting in natural reorientation and precession. We apply the model to idealised simulations of galaxy groups and clusters, finding that jet feedback successfully quenches gas cooling and star formation in all systems. Our group-size halo ($M_\mathrm{200}=10^{13}$ $\mathrm{M}_\odot$) is quenched by a strong jet episode triggered by a cooling flow, and it is kept quenched by a low-power jet fed from hot halo accretion. In more massive systems ($M_\mathrm{200}\geq 10^{14}$ $\mathrm{M}_\odot$), hot halo accretion is insufficient to quench the galaxies, or to keep them quenched after the first cooling episode. These galaxies experience multiple episodes of gas cooling, star formation and jet feedback. In the most massive galaxy cluster that we simulate ($M_\mathrm{200}=10^{15}$ $\mathrm{M}_\odot$), we find peak cold gas masses of $10^{10}$ $\mathrm{M}_\odot$ and peak star formation rates of a few times $100$ $\mathrm{M}_\odot\mathrm{yr}^{-1}$. These values are achieved during strong cooling flows, which also trigger the strongest jets with peak powers of $10^{47}$ $\mathrm{erg}\hspace{0.3mm}\mathrm{s}^{-1}$. These jets subsequently shut off the cooling flows and any associated star formation. Jet-inflated bubbles draw out low-entropy gas that subsequently forms dense cooling filaments in their wakes, as seen in observations. △ Less

Submitted 10 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS

A sparse regression approach for populating dark matter halos and subhalos with galaxies

Authors: M. Icaza-Lizaola, Richard G. Bower, Peder Norberg, Shaun Cole, Matthieu Schaller

Abstract: We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM… ▽ More We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM methods are simple and depend on a relatively small set of adjustable parameters. We collect data from 35,459 galaxies from the EAGLE simulation using 19 redshift slices between $z=0$ and $z=4$ to parameterize the mass evolution of the host halos. Using an appropriate formulation of input parameters, our methodology can model satellite and central halos using a single predictive model that achieves the same accuracy as when predicted separately. This allows us to remove the somewhat arbitrary distinction between those two galaxy types and model them based only on their halo growth history. Our models can accurately reproduce the total galaxy stellar mass function and the stellar mass-dependent galaxy correlation functions ($ξ(r)$) of EAGLE. We show that our SRM model predictions of $ξ(r)$ is competitive with those from sub-halo abundance matching and might be comparable to results from extremely randomized trees. We suggest SRM as an encouraging approach for populating the halos of dark matter only simulations with galaxies and for generating mock catalogues that can be used to explore galaxy evolution or analyse forthcoming large-scale structure surveys. △ Less

Submitted 8 November, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 19 pages, 12 figures

The Milky Way's plane of satellites: consistent with $Λ$CDM

Authors: Till Sawala, Marius Cautun, Carlos S. Frenk, John Helly, Jens Jasche, Adrian Jenkins, Peter H. Johansson, Guilhem Lavaux, Stuart McAlpine, Matthieu Schaller

Abstract: The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is… ▽ More The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is strongly contingent on a lopsided radial distribution, which earlier simulations have failed to reproduce, combined with the close but fleeting conjunction of the two most distant satellites, Leo I and Leo II. Using Gaia proper motions, we show that the orbital pole alignment is much more common than previously reported, and reveal the plane of satellites to be transient rather than rotationally supported. Comparing to new simulations, where such short-lived planes are common, we find the Milky Way satellites to be compatible with standard model expectations. △ Less

Submitted 5 May, 2022; originally announced May 2022.

The interplay between AGN feedback and precipitation of the intracluster medium in simulations of galaxy groups and clusters

Authors: Folkert S. J. Nobels, Joop Schaye, Matthieu Schaller, Yannick M. Bahé, Evgenii Chaikin

Abstract: Using high-resolution hydrodynamical simulations of galaxy clusters, we study the interaction between the brightest cluster galaxy, its supermassive black hole (BH) and the intracluster medium (ICM). We create initial conditions for which the ICM is in hydrostatic equilibrium within the gravitational potential from the galaxy and an NFW dark matter halo. Two free parameters associated with the the… ▽ More Using high-resolution hydrodynamical simulations of galaxy clusters, we study the interaction between the brightest cluster galaxy, its supermassive black hole (BH) and the intracluster medium (ICM). We create initial conditions for which the ICM is in hydrostatic equilibrium within the gravitational potential from the galaxy and an NFW dark matter halo. Two free parameters associated with the thermodynamic profiles determine the cluster gas fraction and the central temperature, where the latter can be used to create cool-core or non-cool-core systems. Our simulations include radiative cooling, star formation, BH accretion, and stellar and active galactic nucleus (AGN) feedback. Even though the energy of AGN feedback is injected thermally and isotropically, it leads to anisotropic outflows and buoyantly rising bubbles. We find that the BH accretion rate (BHAR) is highly variable and only correlates strongly with the star formation rate (SFR) and the ICM when it is averaged over more than $1~\rm Myr$. We generally find good agreement with the theoretical precipitation framework. In $10^{13}~\rm M_\odot$ haloes, AGN feedback quenches the central galaxy and converts cool-core systems into non-cool-core systems. In contrast, higher-mass, cool-core clusters evolve cyclically. Episodes of high BHAR raise the entropy of the ICM out to the radius where the ratio of the cooling time and the local dynamical time $t_{\rm cool}/t_{\rm dyn} > 10$, thus suppressing condensation and, after a delay, the BHAR. The corresponding reduction in AGN feedback allows the ICM to cool and become unstable to precipitation, thus initiating a new episode of high SFR and BHAR. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 22 pages, 15 figures (including appendix); submitted to MNRAS; Supplementary material is available on Youtube at https://youtu.be/HQhc_mytj0A and online on a single website page at https://home.strw.leidenuniv.nl/~nobels/supplementary_idealised_cluster_paper.php

Sensitivity of non-radiative cloud-wind interactions to the hydrodynamics solver

Authors: Joey Braspenning, Joop Schaye, Josh Borrow, Matthieu Schaller

Abstract: Cloud-wind interactions are common in the interstellar and circumgalactic media. Many studies have used simulations of such interactions to investigate the effect of particular physical processes, but the impact of the choice of hydrodynamics solver has largely been overlooked. Here we study the cloud-wind interaction, also known as the "blob test", using seven different hydrodynamics solvers: Thr… ▽ More Cloud-wind interactions are common in the interstellar and circumgalactic media. Many studies have used simulations of such interactions to investigate the effect of particular physical processes, but the impact of the choice of hydrodynamics solver has largely been overlooked. Here we study the cloud-wind interaction, also known as the "blob test", using seven different hydrodynamics solvers: Three flavours of SPH, a moving mesh, adaptive mesh refinement and two meshless schemes. The evolution of masses in dense gas and intermediate-temperature gas, as well as the covering fraction of intermediate-temperature gas, are systematically compared for initial density contrasts of 10 and 100, and five numerical resolutions. To isolate the differences due to the hydrodynamics solvers, we use idealised non-radiative simulations without physical conduction. We find large differences between these methods. SPH methods show slower dispersal of the cloud, particularly for the higher density contrast, but faster convergence, especially for the lower density contrast. Predictions for the intermediate-temperature gas differ particularly strongly, also between non-SPH codes, and converge most slowly. We conclude that the hydrodynamical interaction between a dense cloud and a supersonic wind remains an unsolved problem. Studies aiming to understand the physics or observational signatures of cloud-wind interactions should test the robustness of their results by comparing different hydrodynamics solvers. △ Less

Submitted 28 April, 2023; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: 12 pages, 8 figures; Accepted for publication in MNRAS. Added higher resolution simulations, fixed an issue with the GDF SPH run, clarified the discussion. Conclusions unchanged

The importance of the way in which supernova energy is distributed around young stellar populations in simulations of galaxies

Authors: Evgenii Chaikin, Joop Schaye, Matthieu Schaller, Yannick M. Bahé, Folkert S. J. Nobels, Sylvia Ploeckinger

Abstract: Supernova (SN) feedback plays a crucial role in simulations of galaxy formation. Because blastwaves from individual SNe occur on scales that remain unresolved in modern cosmological simulations, SN feedback must be implemented as a subgrid model. Differences in the manner in which SN energy is coupled to the local interstellar medium and in which excessive radiative losses are prevented have resul… ▽ More Supernova (SN) feedback plays a crucial role in simulations of galaxy formation. Because blastwaves from individual SNe occur on scales that remain unresolved in modern cosmological simulations, SN feedback must be implemented as a subgrid model. Differences in the manner in which SN energy is coupled to the local interstellar medium and in which excessive radiative losses are prevented have resulted in a zoo of models used by different groups. However, the importance of the selection of resolution elements around young stellar particles for SN feedback has largely been overlooked. In this work, we examine various selection methods using the smoothed particle hydrodynamics code SWIFT. We run a suite of isolated disk galaxy simulations of a Milky Way-mass galaxy and small cosmological volumes, all with the thermal stochastic SN feedback model used in the EAGLE simulations. We complement the original mass-weighted neighbour selection with a novel algorithm guaranteeing that the SN energy distribution is as close to isotropic as possible. Additionally, we consider algorithms where the energy is injected into the closest, least dense, or most dense neighbour. We show that different neighbour-selection strategies cause significant variations in star formation rates, gas densities, wind mass loading factors, and galaxy morphology. The isotropic method results in more efficient feedback than the conventional mass-weighted selection. We conclude that the manner in which the feedback energy is distributed among the resolution elements surrounding a feedback event is as important as changing the amount of energy by factors of a few. △ Less

Submitted 22 April, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 16 pages, 13 figures (including appendix); Accepted for publication in MNRAS

SIBELIUS-DARK: a galaxy catalogue of the Local Volume from a constrained realisation simulation

Authors: Stuart McAlpine, John C. Helly, Matthieu Schaller, Till Sawala, Guilhem Lavaux, Jens Jasche, Carlos S. Frenk, Adrian Jenkins, John R. Lucey, Peter H. Johansson

Abstract: We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population… ▽ More We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population calculated using the semi-analytic model of galaxy formation, GALFORM. We demonstrate that the large-scale structure that emerges from the SIBELIUS constrained initial conditions matches well the observational data. The inferred galaxy population of SIBELIUS-DARK also match well the observational data, both statistically for the whole volume and on an object-by-object basis for the most massive clusters. For example, the $K$-band number counts across the whole sky, and when divided between the northern and southern Galactic hemispheres, are well reproduced by SIBELIUS-DARK. We find that the local volume is somewhat unusual in the wider context of $Λ$CDM: it contains an abnormally high number of supermassive clusters, as well as an overall large-scale underdensity at the level of $\approx 5$\% relative to the cosmic mean. However, whilst rare, the extent of these peculiarities does not significantly challenge the $Λ$CDM model. SIBELIUS-DARK is the most comprehensive constrained realisation simulation of the local volume to date, and with this paper we publicly release the halo and galaxy catalogues at $z=0$, which we hope will be useful to the wider astronomy community. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 26 pages, 19 figures. Paper as been accepted for publication in MNRAS. Comments welcome

SEAGLE--III: Towards resolving the mismatch in the dark-matter fraction in early-type galaxies between simulations and observations

Authors: Sampath Mukherjee, Léon V. E. Koopmans, Crescenzo Tortora, Matthieu Schaller, R. Benton Metcalf, Joop Schaye, Georgios Vernardos

Abstract: The central dark-matter fraction of galaxies is sensitive to feedback processes during galaxy formation. Strong gravitational lensing has been effective in the precise measurement of the dark-matter fraction inside massive early-type galaxies. Here, we compare the projected dark-matter fraction of early-type galaxies inferred from the SLACS strong-lens survey, with those obtained from the EAGLE, I… ▽ More The central dark-matter fraction of galaxies is sensitive to feedback processes during galaxy formation. Strong gravitational lensing has been effective in the precise measurement of the dark-matter fraction inside massive early-type galaxies. Here, we compare the projected dark-matter fraction of early-type galaxies inferred from the SLACS strong-lens survey, with those obtained from the EAGLE, Illustris, and IllustrisTNG hydro-dynamical simulations. Previous comparisons with some simulations revealed a large discrepancy, with considerably higher inferred dark-matter fractions -- by factors 2-3 -- inside half of the effective radius in observed strong-lens galaxies as compared to simulated galaxies. Here, we report good agreement between EAGLE and SLACS for the dark-matter fractions inside both half of the effective radius and the effective radius as a function of the galaxy's stellar mass, effective radius, and total mass-density slope. However, for IllustrisTNG and Illustris, the dark-matter fractions are lower than observed. This work consistently assumes a Chabrier IMF, which suggests that a different IMF (although not excluded) is not necessary to resolve this mismatch. The differences in the stellar feedback model between EAGLE and Illustris and IllustrisTNG, are likely the dominant cause of the difference in their dark-matter fraction, and density slope. △ Less

Submitted 12 November, 2021; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: Matches with the published version. Fixed an error in Fig 2. Updated reference. Conclusions unchanged. 7 pages, 4 figures

Observed structural parameters of EAGLE galaxies: reconciling the mass-size relation in simulations with local observations

Authors: Anna de Graaff, James Trayford, Marijn Franx, Matthieu Schaller, Joop Schaye, Arjen van der Wel

Abstract: We use mock images of $z=0.1$ galaxies in the 100 Mpc EAGLE simulation to establish the differences between the sizes and morphologies inferred from the stellar mass distributions and the optical light distributions. The optical, $r$-band images used were constructed with a radiative transfer method to account for the effects of dust, and we measure galaxy structural parameters by fitting Sérsic m… ▽ More We use mock images of $z=0.1$ galaxies in the 100 Mpc EAGLE simulation to establish the differences between the sizes and morphologies inferred from the stellar mass distributions and the optical light distributions. The optical, $r$-band images used were constructed with a radiative transfer method to account for the effects of dust, and we measure galaxy structural parameters by fitting Sérsic models to the images with Galfit. We find that the derived $r$-band half-light radii differ systematically from the stellar half-mass radii, as the $r$-band sizes are typically 0.1 dex larger, and can deviate by as much as $\approx0.5$ dex. The magnitude of this size discrepancy depends strongly on the dust attenuation and star formation activity within the galaxy, as well as the measurement method used. Consequently, we demonstrate that the $r$-band sizes significantly improve the agreement between the simulated and observed stellar mass-size relation: star-forming and quiescent galaxies in EAGLE are typically only slightly larger than observed in the GAMA survey (by 0.1 dex), and the slope and scatter of the local mass-size relation are reproduced well for both populations. Finally, we also compare the obtained morphologies with measurements from GAMA, finding that too few EAGLE galaxies have light profiles that are similar to local early-type galaxies (Sérsic indices of $n\sim 4$). Despite the presence of a significant population of triaxial systems among the simulated galaxies, the surface brightness and stellar mass density profiles tend to be closer to exponential discs ($n\sim1-2$). Our results highlight the need to measure the sizes and morphologies of simulated galaxies using common observational methods in order to perform a meaningful comparison with observations. △ Less

Submitted 2 December, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 22 pages, 18 figures; accepted for publication in MNRAS. The Sersic model parameters are available through the public EAGLE database (see http://icc.dur.ac.uk/Eagle/database.php )

The importance of black hole repositioning for galaxy formation simulations

Authors: Yannick M. Bahé, Joop Schaye, Matthieu Schaller, Richard G. Bower, Josh Borrow, Evgenii Chaikin, Roi Kugel, Folkert Nobels, Sylvia Ploeckinger

Abstract: 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 SM… ▽ More 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. △ Less

Submitted 16 June, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures (including appendix); accepted by MNRAS. Minor updates to match accepted version

A machine learning approach to mapping baryons onto dark matter haloes using the EAGLE and C-EAGLE simulations

Authors: Christopher C. Lovell, Stephen M. Wilkins, Peter A. Thomas, Matthieu Schaller, Carlton M. Baugh, Giulio Fabbian, Yannick Bahé

Abstract: High-resolution cosmological hydrodynamic simulations are currently limited to relatively small volumes due to their computational expense. However, much larger volumes are required to probe rare, overdense environments, and measure clustering statistics of the large scale structure. Typically, zoom simulations of individual regions are used to study rare environments, and semi-analytic models and… ▽ More High-resolution cosmological hydrodynamic simulations are currently limited to relatively small volumes due to their computational expense. However, much larger volumes are required to probe rare, overdense environments, and measure clustering statistics of the large scale structure. Typically, zoom simulations of individual regions are used to study rare environments, and semi-analytic models and halo occupation models applied to dark matter only (DMO) simulations are used to study the Universe in the large-volume regime. We propose a new approach, using a machine learning framework to explore the halo-galaxy relationship in the periodic EAGLE simulations, and zoom C-EAGLE simulations of galaxy clusters. We train a tree based machine learning method to predict the baryonic properties of galaxies based on their host dark matter halo properties. The trained model successfully reproduces a number of key distribution functions for an infinitesimal fraction of the computational cost of a full hydrodynamic simulation. By training on both periodic simulations as well as zooms of overdense environments, we learn the bias of galaxy evolution in differing environments. This allows us to apply the trained model to a larger DMO volume than would be possible if we only trained on a periodic simulation. We demonstrate this application using the $(800 \; \mathrm{Mpc})^3$ P-Millennium simulation, and present predictions for key baryonic distribution functions and clustering statistics from the EAGLE model in this large volume. △ Less

Submitted 2 May, 2023; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 17 pages, 16 figures, Accepted to MNRAS. Updated Acknowledgements

A sparse regression approach to modeling the relation between galaxy stellar masses and their host halos

Authors: M. Icaza-Lizaola, Richard G. Bower, Peder Norberg, Shaun Cole, Matthieu Schaller, Stefan Egan

Abstract: Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 cen… ▽ More Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 central galaxies from the EAGLE hydrodynamic simulation. By matching the host halos to a DM-only simulation, we collect the halo mass and specific angular momentum at present time and for their main progenitors in 10 redshift bins from $z=0$ to $z=4$. The principal component of our governing equation is a third-order polynomial of the host halo mass, which models the stellar-mass halo-mass relation. The scatter about this relation is driven by the halo mass evolution and is captured by second and third-order correlations of the halo mass evolution with the present halo mass. An advantage of sparse regression approaches is that unnecessary terms are removed. Although we include information on halo specific angular momentum, these parameters are discarded by our methodology. This suggests that halo angular momentum has little connection to galaxy formation efficiency. Our model has a root mean square error (RMSE) of $0.167 \log_{10}(M^*/M_\odot)$, and accurately reproduces both the stellar mass function and central galaxy correlation function of EAGLE. The methodology appears to be an encouraging approach for populating the halos of DM-only simulations with galaxies, and we discuss the next steps that are required. △ Less

Submitted 13 August, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: 20 pages, 15 figures

Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations

Authors: Josh Borrow, Matthieu Schaller, Richard G. Bower, Joop Schaye

Abstract: Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was designed with three key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute an… ▽ More Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was designed with three key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute and memory), and to be Lagrangian. Sphenix uses a Density-Energy equation of motion, along with variable artificial viscosity and conduction, including limiters designed to work with common sub-grid models of galaxy formation. In particular, we present and test a novel limiter that prevents conduction across shocks, preventing spurious radiative losses in feedback events. Sphenix is shown to solve many difficult test problems for traditional SPH, including fluid mixing and vorticity conservation, and it is shown to produce convergent behaviour in all tests where this is appropriate. Crucially, we use the same parameters within Sphenix for the various switches throughout, to demonstrate the performance of the scheme as it would be used in production simulations. △ Less

Submitted 28 October, 2021; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: Accepted for publication in MNRAS

Inconsistencies arising from the coupling of galaxy formation sub-grid models to Pressure-Smoothed Particle Hydrodynamics

Authors: Josh Borrow, Matthieu Schaller, Richard G. Bower

Abstract: Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice of variable to smooth in SPH has been the topic of contention, with smoothed pressure (P-SPH) being introduced to reduce errors at contact discontinuities relative to smoothed density schemes. Smoothed pressure… ▽ More Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice of variable to smooth in SPH has been the topic of contention, with smoothed pressure (P-SPH) being introduced to reduce errors at contact discontinuities relative to smoothed density schemes. Smoothed pressure schemes produce excellent results in isolated hydrodynamics tests; in more complex situations however, especially when coupling to the `sub-grid' physics and multiple time-stepping used in many state-of-the-art astrophysics simulations, these schemes produce large force errors that can easily evade detection as they do not manifest as energy non-conservation. Here two scenarios are evaluated: the injection of energy into the fluid (common for stellar feedback) and radiative cooling. In the former scenario, force and energy conservation errors manifest (of the same order as the injected energy), and in the latter large force errors that change rapidly over a few timesteps lead to instability in the fluid (of the same order as the energy lost to cooling). Potential ways to remedy these issues are explored with solutions generally leading to large increases in computational cost. Schemes using a Density-based formulation do not create these instabilities and as such it is recommended that they are preferred over Pressure-based solutions when combined with an energy diffusion term to reduce errors at contact discontinuities. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: Submitted to MNRAS; comments to be sent to joshua.borrow@durham.ac.uk

papaya2: 2D Irreducible Minkowski Tensor computation

Authors: Fabian M. Schaller, Jenny Wagner, Sebastian C. Kapfer

Abstract: A common challenge in scientific and technical domains is the quantitative description of geometries and shapes, e.g. in the analysis of microscope imagery or astronomical observation data. Frequently, it is desirable to go beyond scalar shape metrics such as porosity and surface to volume ratios because the samples are anisotropic or because direction-dependent quantities such as conductances or… ▽ More A common challenge in scientific and technical domains is the quantitative description of geometries and shapes, e.g. in the analysis of microscope imagery or astronomical observation data. Frequently, it is desirable to go beyond scalar shape metrics such as porosity and surface to volume ratios because the samples are anisotropic or because direction-dependent quantities such as conductances or elasticity are of interest. Minkowski Tensors are a systematic family of versatile and robust higher-order shape descriptors that allow for shape characterization of arbitrary order and promise a path to systematic structure-function relationships for direction-dependent properties. Papaya2 is a software to calculate 2D higher-order shape metrics with a library interface, support for Irreducible Minkowski Tensors and interpolated marching squares. Extensions to Matlab, JavaScript and Python are provided as well. While the tensor of inertia is computed by many tools, we are not aware of other open-source software which provides higher-rank shape characterization in 2D. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 5 pages, 3 figures, published in the Journal of Open Source Software, code available at https://morphometry.org/software/papaya2/

Journal ref: Journal of Open Source Software, 5(54) (2020)

Stellar splashback: the edge of the intracluster light

Authors: Alis J. Deason, Kyle A. Oman, Azadeh Fattahi, Matthieu Schaller, Mathilde Jauzac, Yuanyuan Zhang, Mireia Montes, Yannick M. Bahé, Claudio Dalla Vecchia, Scott T. Kay, Tilly A. Evans

Abstract: We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interesti… ▽ More We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the `stellar splashback' feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (mu >> 32 mag/arcsec^2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid. △ Less

Submitted 29 November, 2020; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 12 pages, 9 figures. Accepted for publication in MNRAS