Showing 1–50 of 71 results for author: Bahé, Y

Euclid: Early Release Observations -- The intracluster light and intracluster globular clusters of the Perseus cluster

Authors: M. Kluge, N. A. Hatch, M. Montes, J. B. Golden-Marx, A. H. Gonzalez, J. -C. Cuillandre, M. Bolzonella, A. Lançon, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, A. Boselli, M. Cantiello, J. G. Sorce, F. R. Marleau, P. -A. Duc, E. Sola, M. Urbano, S. L. Ahad, Y. M. Bahé, S. P. Bamford, C. Bellhouse, F. Buitrago, P. Dimauro , et al. (163 additional authors not shown)

Abstract: We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus clu… ▽ More We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus cluster hosts 70000$\pm$2800 GCs and $1.6\times10^{12}$ L$_\odot$ of diffuse light from the BCG+ICL in the near-infrared H$_E$. This accounts for 37$\pm$6% of the cluster's total stellar luminosity within this radius. The ICL and ICGCs share a coherent spatial distribution, suggesting a common origin or that a common potential governs their distribution. Their contours on the largest scales (>200 kpc) are offset from the BCG's core westwards by 60 kpc towards several luminous cluster galaxies. This offset is opposite to the displacement observed in the gaseous intracluster medium. The radial surface brightness profile of the BCG+ICL is best described by a double Sérsic model, with 68$\pm$4% of the H$_E$ light in the extended, outer component. The transition between these components occurs at ~50 kpc, beyond which the isophotes become increasingly elliptical and off-centred. The radial ICGC number density profile closely follows the BCG+ICL profile only beyond this 50 kpc radius, where we find an average of 60 GCs per $10^9$ M$_\odot$ of diffuse stellar mass. The BCG+ICL colour becomes increasingly blue with radius, consistent with the stellar populations in the ICL having subsolar metallicities [Fe/H]~-0.6. The colour of the ICL, and the specific frequency and luminosity function of the ICGCs suggest that the ICL+ICGCs were tidally stripped from the outskirts of massive satellites with masses of a few $\times10^{10}$ M$_\odot$, with an increasing contribution from dwarf galaxies at large radii. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations. 24 pages, 17 figures, 5 tables. Submitted to A&A

Resolved UV and optical color gradients reveal environmental influence on galaxy evolution at redshift z$\sim$1.6

Authors: William J. Cramer, A. G. Noble, G. Rudnick, A. Pigarelli, G. Wilson, Y. M. Bahé, M. C. Cooper, R. Demarco, J. Matharu, T. B. Miller, A. Muzzin, J. Nantais, W. Sportsman, E. van Kampen, T. M. A. Webb, H. K. C. Yee

Abstract: The changes in colors across a galaxy are intimately connected to the galaxy's formation, growth, quenching history, and dust content. A particularly important epoch in the growth of galaxies is near $z \sim 2$ often referred to as 'cosmic noon', where galaxies on average reach the peak of their star formation. We study a population of 125 cluster galaxies at $z \sim 1.6$ in three Hubble Space Tel… ▽ More The changes in colors across a galaxy are intimately connected to the galaxy's formation, growth, quenching history, and dust content. A particularly important epoch in the growth of galaxies is near $z \sim 2$ often referred to as 'cosmic noon', where galaxies on average reach the peak of their star formation. We study a population of 125 cluster galaxies at $z \sim 1.6$ in three Hubble Space Telescope (HST) filters, F475W, F625W, and F160W, roughly corresponding to the rest-frame FUV, NUV, and r band, respectively. By comparing to a control sample of 200 field galaxies at similar redshift, we reveal clear, statistically significant differences in the overall spatially resolved colors and color gradients in galaxies across these two different environments. On average, cluster galaxies have redder UV colors in both the inner and outer regions bounded by $r_{\mathrm{50}}$, as well as an overall wider dispersion of outside-in color gradients. The presence of these observed differences, along with evidence from ancillary data from previous studies, strongly suggests that the environment drives these population-level color differences, by affecting the stellar populations and/or dust content. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

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

The first quenched galaxies, when and how?

Authors: Lizhi Xie, Gabriella De Lucia, Fabio Fontanot, Michaela Hirschmann, Yannick M Bahé, Michael L. Balogh, Adam Muzzin, Benedetta Vulcani, Devontae C. Baxter, Ben Forrest, Gillian Wilson, Gregory H. Rudnick, M. C. Cooper, Umberto Rescigno

Abstract: Many quiescent galaxies discovered in the early Universe by \textit{JWST} raise fundamental questions on when and how these galaxies became and stayed quenched. Making use of the latest version of the semi-analytic model GAEA that provides good agreement with the observed quenched fractions up to $z\sim 3$, we make predictions for the expected fractions of quiescent galaxies up to $z\sim 7$ and an… ▽ More Many quiescent galaxies discovered in the early Universe by \textit{JWST} raise fundamental questions on when and how these galaxies became and stayed quenched. Making use of the latest version of the semi-analytic model GAEA that provides good agreement with the observed quenched fractions up to $z\sim 3$, we make predictions for the expected fractions of quiescent galaxies up to $z\sim 7$ and analyze the main quenching mechanism. We find that in a simulated box of $685~{\rm Mpc}$ on a side, the first quenched massive ($M_{\star} \sim 10^{11} {\rm M}_{\odot}$), Milky Way mass, and low mass ($M_{\star} \sim 10^{9.5} {\rm M}_{\odot}$ ) galaxies appear at $z\sim 4.5$, $z\sim 6.2$, and before $z = 7$. Most quenched galaxies identified at early redshifts remain quenched for more than 1 Gyr. Independently of galaxy stellar mass, the dominant quenching mechanism at high redshift is accretion disk feedback (quasar winds) from a central massive black hole, which is triggered by mergers in massive and MW-mass galaxies, and by disk instabilities in low-mass galaxies. Environmental stripping becomes increasingly more important at lower redshift. △ Less

Submitted 1 April, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 8 pages, 3 figures, Accepted publication in APJL

$Spitzer$-selected $z > 1.3$ protocluster candidates in the LSST Deep Drilling Fields

Authors: Harry Gully, Nina Hatch, Yannick Bahé, Michael Balogh, Micol Bolzonella, M. C. Cooper, Adam Muzzin, Lucia Pozzetti, Gregory Rudnick, Benedetta Vulcani, Gillian Wilson

Abstract: We have identified 189 candidate $z > 1.3$ protoclusters and clusters in the LSST Deep Drilling Fields. This sample will enable the measurement of the metal enrichment and star formation history of clusters during their early assembly period through the direct measurement of the rate of supernovae identified through the LSST. The protocluster sample was selected from galaxy overdensities in a… ▽ More We have identified 189 candidate $z > 1.3$ protoclusters and clusters in the LSST Deep Drilling Fields. This sample will enable the measurement of the metal enrichment and star formation history of clusters during their early assembly period through the direct measurement of the rate of supernovae identified through the LSST. The protocluster sample was selected from galaxy overdensities in a $Spitzer$/IRAC colour-selected sample using criteria that were optimised for protocluster purity using a realistic lightcone. Our tests reveal that $60-80\%$ of the identified candidates are likely to be genuine protoclusters or clusters, which is corroborated by a $\sim4σ$ stacked X-ray signal from these structures. We provide photometric redshift estimates for 47 candidates which exhibit strong peaks in the photo-$z$ distribution of their candidate members. However, the lack of a photo-$z$ peak does not mean a candidate is not genuine, since we find a stacked X-ray signal of similar significance from both the candidates that exhibit photo-$z$ peaks and those that do not. Tests on the lightcone reveal that our pursuit of a pure sample of protoclusters results in that sample being highly incomplete ($\sim4\%$) and heavily biased towards larger, richer, more massive, and more centrally concentrated protoclusters than the total protocluster population. Most ($\sim75\%$) of the selected protoclusters are likely to have a maximum collapsed halo mass of between $10^{13}-10^{14}$ M$_{\odot}$, with only $\sim25\%$ likely to be collapsed clusters above $10^{14}$ M$_{\odot}$. However, the aforementioned bias ensures our sample is $\sim50\%$ complete for structures that have already collapsed into clusters more massive than $10^{14}$ M$_{\odot}$. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 18 pages, 14 figures

Journal ref: MNRAS, Volume 527, Issue 4, February 2024, Pages 10680-10696

The stellar mass function of quiescent galaxies in 2 < z < 2.5 protoclusters

Authors: Adit H. Edward, Michael L. Balogh, Yannick M. Bahe, Michael C. Cooper, Nina A. Hatch, Justin Marchioni, Adam Muzzin, Allison Noble, Gregory H. Rednick, Benedetta Vulcani, Gillian Wilson, Gabriella De Lucia, Ricardo Demarco, Ben Forrest, Michaela Hirschmann, Gianluca Castignani, Pierluigi Cerulo, Rose A. Finn, Guillaume Hewitt, Pascale Jablonka, Yadayuki Kodama, Sophie Maurogordato, Julie Nantais, Lizhi Xie

Abstract: We present an analysis of the galaxy stellar mass function (SMF) of 14 known protoclusters between $2.0 < z < 2.5$ in the COSMOS field, down to a mass limit of $10^{9.5}$ M$_{\odot}$. We use existing photometric redshifts with a statistical background subtraction, and consider star-forming and quiescent galaxies identified from $(NUV - r)$ and $(r - J)$ colours separately. Our fiducial sample incl… ▽ More We present an analysis of the galaxy stellar mass function (SMF) of 14 known protoclusters between $2.0 < z < 2.5$ in the COSMOS field, down to a mass limit of $10^{9.5}$ M$_{\odot}$. We use existing photometric redshifts with a statistical background subtraction, and consider star-forming and quiescent galaxies identified from $(NUV - r)$ and $(r - J)$ colours separately. Our fiducial sample includes galaxies within 1 Mpc of the cluster centres. The shape of the protocluster SMF of star-forming galaxies is indistinguishable from that of the general field at this redshift. Quiescent galaxies, however, show a flatter SMF than in the field, with an upturn at low mass, though this is only significant at $\sim 2σ$. There is no strong evidence for a dominant population of quiescent galaxies at any mass, with a fraction of $< 15\%$ at $1σ$ confidence for galaxies with log$M_{\ast}/M_{\odot} < 10.5$. We compare our results with a sample of galaxies groups at $1 < z < 1.5$, and demonstrate that a significant amount of environmental quenching must take place between these epochs, increasing the relative abundance of high-mass ($\rm M > 10^{10.5} M_{\odot}$) quiescent galaxies by a factor of $\gtrsim$ 2. However, we find that at lower masses ($\rm M < 10^{10.5} M_{\odot}$), no additional environmental quenching is required. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 23 pages, 22 figures. Accepted for publication in MNRAS

Preparing for low surface brightness science with the Vera C. Rubin Observatory: A Comparison of Observable and Simulated Intracluster Light Fractions

Authors: Sarah Brough, Syeda Lammim Ahad, Yannick M. Bahe, Amaël Ellien, Anthony H. Gonzalez, Yolanda Jiménez-Teja, Lucas C. Kimmig, Garreth Martin, Cristina Martínez-Lombilla, Mireia Montes, Annalisa Pillepich, Rossella Ragusa, Rhea-Silvia Remus, Chris A. Collins, Johan H. Knapen, J. Chris Mihos

Abstract: Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-S… ▽ More Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-SB, Non-Parametric Measure-NP, Composite Models-CM, Multi-Galaxy Fitting-MGF) and new approaches under development (Wavelet Decomposition-WD) applied to mock images of 61 galaxy clusters (14<log10 M_200c/M_solar <14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions respectively. The Ahad (CM), MGF and WD algorithms are best set up to process larger samples, however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory's Legacy Survey of Space and Time will provide. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Submitted for publication in MNRAS, posted to arXiv after responding to two positive rounds of referee comments. Key results in Figs 3, 5, 6 and 11

Prospects for detecting the circum- and intergalactic medium in X-ray absorption using the extended intracluster medium as a backlight

Authors: Lýdia Štofanová, Aurora Simionescu, Nastasha A. Wijers, Joop Schaye, Jelle S. Kaastra, Yannick M. Bahé, Andrés Arámburo-García

Abstract: The warm-hot plasma in cosmic web filaments is thought to comprise a large fraction of the gas in the local Universe. So far, the search for this gas has focused on mapping its emission, or detecting its absorption signatures against bright, point-like sources. Future, non-dispersive, high spectral resolution X-ray detectors will, for the first time, enable absorption studies against extended obje… ▽ More The warm-hot plasma in cosmic web filaments is thought to comprise a large fraction of the gas in the local Universe. So far, the search for this gas has focused on mapping its emission, or detecting its absorption signatures against bright, point-like sources. Future, non-dispersive, high spectral resolution X-ray detectors will, for the first time, enable absorption studies against extended objects. Here, we use the Hydrangea cosmological hydrodynamical simulations to predict the expected properties of intergalactic gas in and around massive galaxy clusters, and investigate the prospects of detecting it in absorption against the bright cores of nearby, massive, relaxed galaxy clusters. We probe a total of $138$ projections from the simulation volumes, finding $16$ directions with a total column density $N_{O VII} > 10^{14.5}$ cm$^{-2}$. The strongest absorbers are typically shifted by $\pm 1000$ km/s with respect to the rest frame of the cluster they are nearest to. Realistic mock observations with future micro-calorimeters, such as the Athena X-ray Integral Field Unit or the proposed Line Emission Mapper (LEM) X-ray probe, show that the detection of cosmic web filaments in O VII and O VIII absorption against galaxy cluster cores will be feasible. An O VII detection with a $5σ$ significance can be achieved in $10-250$ ks with Athena for most of the galaxy clusters considered. The O VIII detection becomes feasible only with a spectral resolution of around $1$ eV, comparable to that envisioned for LEM. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 17 pages, 10 figures, 4 tables, accepted for publication in MNRAS

A case study of early galaxy cluster with the Athena X-IFU

Authors: F. Castellani, N. Clerc, E. Pointecouteau, Y. M. Bahé, J. Schaye, F. Pajot

Abstract: Context: Observations of the hot gas in distant clusters of galaxies, though challenging, are key to understand the role of intense galaxy activity, super-massive black hole feedback and chemical enrichment in the process of massive halos assembly. Aims: We assess the feasibility to retrieve, using X-ray hyperspectral data only, the thermodymamical hot gas properties and chemical abundances of a… ▽ More Context: Observations of the hot gas in distant clusters of galaxies, though challenging, are key to understand the role of intense galaxy activity, super-massive black hole feedback and chemical enrichment in the process of massive halos assembly. Aims: We assess the feasibility to retrieve, using X-ray hyperspectral data only, the thermodymamical hot gas properties and chemical abundances of a $z=2$ galaxy cluster of mass M500=7 x $10^{13} M_{\odot}$, extracted from the Hydrangea hydrodynamical simulation. Methods: We create mock X-ray observations of the future X-ray Integral Field Unit (X-IFU) onboard the Athena mission. By forward-modeling the measured 0.4-1 keV surface brightness, the projected gas temperature and abundance profiles, we reconstruct the three-dimensional distribution for the gas density, pressure, temperature and entropy. Results: Thanks to its large field-of-view, high throughput and exquisite spectral resolution, one X-IFU exposure lasting 100ks enables reconstructing density and pressure profiles with 20% precision out to a characteristic radius of R500, accounting for each quantity's intrinsic dispersion in the Hydrangea simulations. Reconstruction of abundance profiles requires both higher signal-to-noise ratios and specific binning schemes. We assess the enhancement brought by longer exposures and by observing the same object at later evolutionary stages ($z=1-1.5$). Conclusions: Our analysis highlights the importance of scatter in the radially binned gas properties, which induces significant effects on the observed projected quantities. The fidelity of the reconstruction of gas profiles is sensitive to the degree of gas components mixing along the line-of-sight. Future analyses should aim at involving dedicated hyper-spectral models and fitting methods that are able to grasp the complexity of such three-dimensional, multi-phase, diffuse gas structures. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 15 pages, 11 figures, 3 tables. Accepted for publication in A&A

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

VERTICO VII: Environmental quenching caused by suppression of molecular gas content and star formation efficiency in Virgo Cluster galaxies

Authors: Toby Brown, Ian D. Roberts, Mallory Thorp, Sara L. Ellison, Nikki Zabel, Christine D. Wilson, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Alessandro Boselli, Aeree Chung, Luca Cortese, Barbara Catinella, Timothy A. Davis, María J. Jiménez-Donaire, Claudia D. P. Lagos, Bumhyun Lee, Laura C. Parker, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Vicente Villanueva, Adam B. Watts

Abstract: We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with resp… ▽ More We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

An environment-dependent halo mass function as a driver for the early quenching of $z\geq1.5$ cluster galaxies

Authors: Syeda Lammim Ahad, Adam Muzzin, Yannick M. Bahé, Henk Hoekstra

Abstract: Many $z=1.5$ galaxies with a stellar mass ($M_{\star}$) $\geq 10^{10}\,\mathrm{M}_\odot$ are already quenched in both galaxy clusters ($>50$ per cent) and the field ($>20$ per cent), with clusters having a higher quenched fraction at all stellar masses compared to the field. A puzzling issue is that these massive quenched galaxies have stellar populations of similar age in both clusters and the fi… ▽ More Many $z=1.5$ galaxies with a stellar mass ($M_{\star}$) $\geq 10^{10}\,\mathrm{M}_\odot$ are already quenched in both galaxy clusters ($>50$ per cent) and the field ($>20$ per cent), with clusters having a higher quenched fraction at all stellar masses compared to the field. A puzzling issue is that these massive quenched galaxies have stellar populations of similar age in both clusters and the field. This suggests that, despite the higher quenched fraction in clusters, the dominant quenching mechanism for massive galaxies is similar in both environments. In this work, we use data from the cosmological hydrodynamic simulations Hydrangea and EAGLE to test whether the excess quenched fraction of massive galaxies in $z = 1.5$ clusters results from fundamental differences in their halo properties compared to the field. We find that (i) at $10^{10} \leq$ $M_{\star}/\,\mathrm{M}_\odot\leq 10^{11}$, quenched fractions in the redshift range $1.5 < z < 3.5$ are consistently higher for galaxies with higher peak maximum circular velocity of the dark matter halo ($v_{\mathrm{max, peak}}$), and (ii) the distribution of $v_{\mathrm{max, peak}}$ is strongly biased towards higher values for cluster satellites compared to the field centrals. Due to this difference in the halo properties of cluster and field galaxies, secular processes alone may account for (most of) the environmental excess of massive quenched galaxies in high-redshift (proto) clusters. Taken at face value, our results challenge a fundamental assumption of popular quenching models, namely that clusters are assembled from an unbiased subset of infalling field galaxies. If confirmed, this would imply that such models must necessarily fail at high redshift, as indicated by recent observations. △ Less

Submitted 5 February, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 11 pages, 4+2 figures (including appendices), Accepted for publication in MNRAS

Roman-Cosmic Noon: A Legacy Spectroscopic Survey of Massive Field and Protocluster Galaxies at $2<z<3$

Authors: Gregory Rudnick, Yannick Bahé, Michael Balogh, Mike Cooper, Nina Hatch, Benedetta Vulcani, Gillian Wilson, Gianluca Castignani, Pierluigi Cerulo, Gabriella De Lucia, Ricardo DeMarco, Benjamin Forrest, Pascale Jablonka

Abstract: Protoclusters are the densest regions in the distant universe ($z>2$) and are the progenitors of massive galaxy clusters ($M_{halo}>10^{14}{\rm M}_\odot$) in the local universe. They undoubtedly play a key role in early massive galaxy evolution and they may host the earliest sites of galaxy quenching or even induce extreme states of star formation. Studying protoclusters therefore not only gives u… ▽ More Protoclusters are the densest regions in the distant universe ($z>2$) and are the progenitors of massive galaxy clusters ($M_{halo}>10^{14}{\rm M}_\odot$) in the local universe. They undoubtedly play a key role in early massive galaxy evolution and they may host the earliest sites of galaxy quenching or even induce extreme states of star formation. Studying protoclusters therefore not only gives us a window into distant galaxy formation but also provides an important link in our understanding of how dense structures grow over time and modify the galaxies within them. Current protocluster samples are completely unable to address these points because they are small and selected in a heterogeneous way. We propose the Roman-Cosmic Noon survey, whose centerpiece is an extremely deep (30ksec) and wide area (10 deg$^2$) prism slitless spectroscopy survey to identify the full range of galaxy structures at $2<z<3$. This survey will include 1500 uniformly selected protoclusters, their surrounding cosmic web environments, and at least 15,000 protocluster galaxies with $M_\star>10^{10.5} {\rm M}_\odot$ across the full range of star formation histories as well as many more lower mass star-forming galaxies. The survey will also contain field galaxies to much lower masses than in the High Latitude Wide Area Survey, but over an area dwarfing any current or planned deep spectroscopy probe at $z>2$. With the prism spectroscopy and some modest additional imaging this survey will measure precise stellar mass functions, quenched fractions, galaxy and protocluster morphologies, stellar ages, emission-line based SFRs, and metallicities. It will have extensive legacy value well beyond the key protocluster science goals. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: Submitted to the 2023 call for White Papers regarding the Roman Mission Core Community Surveys, 6 pages, 3 figures

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

Intracluster light in the core of z~2 galaxy proto-clusters

Authors: S. V. Werner, N. A. Hatch, J. Matharu, A. H. Gonzalez, Y. M. Bahé, S. Mei, G. Noirot, D. Wylezalek

Abstract: Intracluster light is thought to originate from stars that were ripped away from their parent galaxies by gravitational tides and galaxy interactions during the build up of the cluster. The stars from such interactions will accumulate over time, so semi-analytic models suggest that the abundance of intracluster stars is negligible in young proto-clusters at z~2 and grows to around a quarter of the… ▽ More Intracluster light is thought to originate from stars that were ripped away from their parent galaxies by gravitational tides and galaxy interactions during the build up of the cluster. The stars from such interactions will accumulate over time, so semi-analytic models suggest that the abundance of intracluster stars is negligible in young proto-clusters at z~2 and grows to around a quarter of the stellar mass in the oldest, most mature clusters. In contrast to these theoretical expectations, we report on the detection of intracluster light within two proto-clusters at z=2 using deep HST images. We use the colour of the intracluster light to estimate its mass-to-light ratio in annuli around the brightest cluster galaxies (BCG), up to a radius of 100 kpc. We find that $54\pm5$% and $71\pm3$% of the stellar mass in these regions is located more than 10 kpc away from the BCGs in the two proto-clusters. This low concentration is similar to BCGs in lower redshift clusters, and distinct from other massive proto-cluster galaxies. This suggests that intracluster stars are already present within the core 100 kpc of proto-clusters. We compare these observations to the Hydrangea hydrodynamical galaxy cluster simulations and find that intracluster stars are predicted to be a generic feature of group-sized halos at z=2. These intracluster stars will gradually move further away from the BCG as the proto-cluster assembles into a cluster. △ Less

Submitted 12 May, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: 14 pages, 9 figures, accepted for publication in MNRAS

VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Toby Brown, Christine D. Wilson, Nikki Zabel, Ian D. Roberts, Timothy A. Davis, Mallory Thorp, Aeree Chung, Adam R. H. Stevens, Sara L. Ellison, Kristine Spekkens, Laura C. Parker, Yannick M. Bahé, Vicente Villanueva, María Jiménez-Donaire, Dhruv Bisaria, Alessandro Boselli, Alberto D. Bolatto, Bumhyun Lee

Abstract: The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virg… ▽ More The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($Σ_{\star}$), HI ($Σ_\mathrm{HI}$), and molecular gas ($Σ_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average $Σ_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower $Σ_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical $Σ_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $Σ_\mathrm{HI}$ or $Σ_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of $Σ_\mathrm{mol}$-to-$Σ_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

Cosmic gas highways in C-EAGLE simulations

Authors: I. Vurm, J. Nevalainen, S. E. Hong, Y. M. Bahé, C. Dalla Vecchia, P. Heinämäki

Abstract: A substantial fraction of the cosmic baryons is expected to hide in the form of diffuse warm-hot intergalactic medium (WHIM), the majority of which resides in the filaments of the Cosmic Web and has proven very difficult to detect due to its low density. Close to galaxy clusters, the filament gas is affected by the cluster's gravitational potential and attains substantial infall velocities, eventu… ▽ More A substantial fraction of the cosmic baryons is expected to hide in the form of diffuse warm-hot intergalactic medium (WHIM), the majority of which resides in the filaments of the Cosmic Web and has proven very difficult to detect due to its low density. Close to galaxy clusters, the filament gas is affected by the cluster's gravitational potential and attains substantial infall velocities, eventually undergoing a termination shock which may boost its X-ray signal. We aim to identify the optimal locations of the enhanced X-ray emission and absorption arising from cluster-filament interactions, as well as improve our understanding of the various physical processes affecting the WHIM as it approaches the cluster. We applied the DisPerSE filament finder to the galaxy distribution in the surroundings of a Coma-like ($M_{200} \sim 10^{15.4}~{\rm M}_{\odot}$) simulated C-EAGLE galaxy cluster. We characterised the thermodynamic properties of the gas in such filaments as well as their dependence on distance from the cluster, and provided a physical interpretation for the results. The identified filaments account for $\sim 50$% of the hot WHIM ($T > 10^{5.5}$ K) in the cluster vicinity. The filament gas is in approximate free-fall all the way down to $\sim 2 \ r_{200}$ from the cluster, at which stage it begins to slow down due to the increasing pressure of the ambient gas. The deceleration is accompanied by the conversion of gas bulk kinetic energy into heat and the increase of density and temperature from the general Cosmic Web level of $ρ\sim 10ρ_{\rm av}$ and $T = 10^5-10^6$ K towards $ρ\sim 100ρ_{\rm av}$ and $T = 10^7-10^8$ K near the cluster boundary. We conclude that the detection of the cosmic filaments of galaxies around clusters may provide a practical observational avenue for locating the densest and hottest phase of the missing baryons. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 20 pages, 16 figures, accepted for publication in A&A

Star formation concentration as a tracer of environmental quenching in action: a study of the Eagle and C-Eagle simulations

Authors: Di Wang, Claudia D. P. Lagos, Scott M. Croom, Ruby J. Wright, Yannick M. Bahé, Julia J. Bryant, Jesse van de Sande, Sam P. Vaughan

Abstract: We study environmental quenching in the Eagle}/C-Eagle cosmological hydrodynamic simulations over the last 11 Gyr (i.e. $z=0-2$). The simulations are compared with observations from the SAMI Galaxy Survey at $z=0$. We focus on satellite galaxies in galaxy groups and clusters ($10^{12}\,\rm M_{\odot}$ $\lesssim$ $M_{200}$ < $3 \times 10^{15}\, \rm M_{\odot}$). A star-formation concentration index [… ▽ More We study environmental quenching in the Eagle}/C-Eagle cosmological hydrodynamic simulations over the last 11 Gyr (i.e. $z=0-2$). The simulations are compared with observations from the SAMI Galaxy Survey at $z=0$. We focus on satellite galaxies in galaxy groups and clusters ($10^{12}\,\rm M_{\odot}$ $\lesssim$ $M_{200}$ < $3 \times 10^{15}\, \rm M_{\odot}$). A star-formation concentration index [$C$-index $= \log_{10}(r_\mathrm{50,SFR} / r_\mathrm{50,rband})$] is defined, which measures how concentrated star formation is relative to the stellar distribution. Both Eagle/C-Eagle and SAMI show a higher fraction of galaxies with low $C$-index in denser environments at $z=0-0.5$. Low $C$-index galaxies are found below the SFR-$M_{\star}$ main sequence (MS), and display a declining specific star formation rate (sSFR) with increasing radii, consistent with ``outside-in'' environmental quenching. Additionally, we show that $C$-index can be used as a proxy for how long galaxies have been satellites. These trends become weaker at increasing redshift and are absent by $z=1-2$. We define a quenching timescale $t_{\rm quench}$ as how long it takes satellites to transition from the MS to the quenched population. We find that simulated galaxies experiencing ``outside-in'' environmental quenching at low redshift ($z=0\sim0.5$) have a long quenching timescale (median $t_{\rm quench}$ > 2 Gyr). The simulated galaxies at higher redshift ($z=0.7\sim2$) experience faster quenching (median $t_{\rm quench}$ < 2Gyr). At $z\gtrsim 1-2$ galaxies undergoing environmental quenching have decreased sSFR across the entire galaxy with no ``outside-in'' quenching signatures and a narrow range of $C$-index, showing that on average environmental quenching acts differently than at $z\lesssim 1$. △ Less

Submitted 18 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 21 pages, 17 figures,

Early results from GLASS-JWST. XIV: A spectroscopically confirmed protocluster 650 million years after the Big Bang

Authors: Takahiro Morishita, Guido Roberts-Borsani, Tommaso Treu, Gabriel Brammer, Charlotte A. Mason, Michele Trenti, Benedetta Vulcani, Xin Wang, Ana Acebron, Yannick Bahé, Pietro Bergamini, Kristan Boyett, Marusa Bradac, Antonello Calabrò, Marco Castellano, Wenlei Chen, Gabriella De Lucia, Alexei V. Filippenko, Adriano Fontana, Karl Glazebrook, Claudio Grillo, Alaina Henry, Tucker Jones, Patrick L. Kelly, Anton M. Koekemoer , et al. (17 additional authors not shown)

Abstract: We present the spectroscopic confirmation of a protocluster at $z=7.88$ behind the galaxy cluster Abell2744 (hereafter A2744-z7p9OD). Using JWST NIRSpec, we find seven galaxies within a projected radius of 60kpc. Although the galaxies reside in an overdensity around $>20\times$ greater than a random volume, they do not show strong Lyman-alpha emission. We place 2-$σ$ upper limits on the rest-frame… ▽ More We present the spectroscopic confirmation of a protocluster at $z=7.88$ behind the galaxy cluster Abell2744 (hereafter A2744-z7p9OD). Using JWST NIRSpec, we find seven galaxies within a projected radius of 60kpc. Although the galaxies reside in an overdensity around $>20\times$ greater than a random volume, they do not show strong Lyman-alpha emission. We place 2-$σ$ upper limits on the rest-frame equivalent width $<16$-$28$AA. Based on the tight upper limits to the Lyman-alpha emission, we constrain the volume-averaged neutral fraction of hydrogen in the intergalactic medium to be $x_{\rm HI} > 0.45$ (68% CI). Using an empirical $M_{\rm UV}$-$M_{\rm halo}$ relation for individual galaxies, we estimate that the total halo mass of the system is $\gtrsim 4\times10^{11}\,M_\odot$. Likewise, the line of sight velocity dispersion is estimated to be $1100 \pm 200$km/s. Using an empirical relation, we estimate the present-day halo mass of A2744-z7p9OD to be $\sim2\times10^{15}\,M_\odot$, comparable to the Coma cluster. A2744-z7p9OD is the highest redshift spectroscopically confirmed protocluster to date, demonstrating the power of JWST to investigate the connection between dark-matter halo assembly and galaxy formation at very early times with medium-deep observations at $<20$hrs total exposure time. Follow-up spectroscopy of the remaining photometric candidates of the overdensity will further refine the features of this system and help characterize the role of such overdensities in cosmic reionization. △ Less

Submitted 24 April, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: The title has been updated to reflect the published numbering; a minor change has been made to Figure 1 with regard to the MSA shutters on the rgb stamp image. NASA press release article can be found at: https://www.nasa.gov/feature/goddard/2023/webb-reveals-early-universe-prequel-to-huge-galaxy-cluster

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

How to Interpret Measurements of Diffuse Light in Stacked Observations of Groups and Clusters of Galaxies

Authors: Syeda Lammim Ahad, Yannick M. Bahé, Henk Hoekstra

Abstract: The diffuse light within galaxy groups and clusters provides valuable insight into the growth of massive cosmic structures. Groups are particularly interesting in this context, because they represent the link between galactic haloes and massive clusters. However, low surface brightness makes their diffuse light extremely challenging to detect individually. Stacking many groups is a promising alter… ▽ More The diffuse light within galaxy groups and clusters provides valuable insight into the growth of massive cosmic structures. Groups are particularly interesting in this context, because they represent the link between galactic haloes and massive clusters. However, low surface brightness makes their diffuse light extremely challenging to detect individually. Stacking many groups is a promising alternative, but its physical interpretation is complicated by possible systematic variations of diffuse light profiles with other group properties. Another issue is the often ambiguous choice of group centre. We explore these challenges using mock observations for 497 galaxy groups and clusters with halo masses from $~ 10^{12} \textrm{M}_{\odot}$ to $1.5 \times 10^{15}\textrm{M}_{\odot}$ at redshift $0.1$ from the Hydrangea cosmological hydrodynamic simulations. In 18 per cent of groups with at least five galaxies above $10^{9} \textrm{M}_{\odot}$ in stellar mass, the $r$-band brightest galaxy is not the one at the centre of the gravitational potential; line-of-sight projections account for half of these cases. Miscentring does not significantly affect the ensemble average mass density profile or the surface brightness profile for our sample: even within ambiguously centred haloes, different centring choices lead to only a 1 per cent change in the total fraction of diffuse intra-group light, $f_{\textrm{IGL}}$. We find strong correlations of $f_{\textrm{IGL}}$ with the luminosity of the central group galaxy and halo mass. Stacking groups in narrow bins of central galaxy luminosity will therefore make the physical interpretation of the signal more straightforward than combining systems across a wide range of mass. △ Less

Submitted 18 November, 2022; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: 18 pages, 12 figures, Accepted for publication in MNRAS

GOGREEN: a critical assessment of environmental trends in cosmological hydrodynamical simulations at z ~ 1

Authors: Egidijus Kukstas, Michael L. Balogh, Ian G. McCarthy, Yannick M. Bahe, Gabriella De Lucia, Pascale Jablonka, Benedetta Vulcani, Devontae C. Baxter, Andrea Biviano, Pierluigi Cerulo, Jeffrey C. Chan, M. C. Cooper, Ricardo Demarco, Alexis Finoguenov, Andreea S. Font, Chris Lidman, Justin Marchioni, Sean McGee, Adam Muzzin, Julie Nantais, Lyndsay Old, Irene Pintos-Castro, Bianca Poggianti, Andrew M. M. Reeves, Gregory Rudnick , et al. (6 additional authors not shown)

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

Submitted 19 October, 2022; originally announced October 2022.

Comments: 20 pages, 10 figures, resubmitted to MNRAS following minor revisions

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

X-IFU/Athena view of the most distant galaxy clusters in the Universe

Authors: Florent Castellani, N. Clerc, Etienne Pointecouteau, Yannick Bahé, F. Pajot

Abstract: The X-ray Integral Field Unit (X-IFU) on-board the second large ESA mission ''Athena'' will be a high spatial (5'') and spectral (2.5eV) resolution X-ray imaging spectrometer, operating in the 0.2-12 keV energy band. It will address the science question of the assembly and evolution through cosmic time of the largest halos of matter in the Universe, groups and clusters of galaxies. To this end, we… ▽ More The X-ray Integral Field Unit (X-IFU) on-board the second large ESA mission ''Athena'' will be a high spatial (5'') and spectral (2.5eV) resolution X-ray imaging spectrometer, operating in the 0.2-12 keV energy band. It will address the science question of the assembly and evolution through cosmic time of the largest halos of matter in the Universe, groups and clusters of galaxies. To this end, we present an on-going feasibility study to demonstrate the X-IFU capabilities to unveil the physics of massive halos at their epoch of formation. Starting from a distant (z=2) group of galaxies ($M_{500} = 7\cdot 10^{13} M_\odot/h$) extracted from the HYDRANGEA cosmological and hydrodynamical numerical simulations, we perform an end-to-end simulation of X-IFU observations. From the reconstruction of the global, 1D and 2D quantities, we plan to investigate the various X-IFU science cases for clusters of galaxies, such as the chemical enrichment of the intra-cluster medium (ICM), the dynamical assembly of groups and clusters and the impact of feedback from galaxy and super-massive black hole evolution. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: Journées 2021 de la Société Française d'Astronomie \& d'Astrophysique, Jun 2021, CarbonFreeConf, France

The luminosity of cluster galaxies in the Cluster-EAGLE simulations

Authors: Andrea Negri, Claudio Dalla Vecchia, J. Alfonso L. Aguerri, Yannick Bahé

Abstract: We computed the luminosity of simulated galaxies of the C-EAGLE project, a suite of 30 high-resolution zoom-in simulations of galaxy clusters based on the EAGLE simulation. The AB magnitudes are derived for different spectral bands, from ultraviolet to infrared, using the simple stellar population modeling based on the E-MILES stellar spectra library. We take into account obscuration due to dust i… ▽ More We computed the luminosity of simulated galaxies of the C-EAGLE project, a suite of 30 high-resolution zoom-in simulations of galaxy clusters based on the EAGLE simulation. The AB magnitudes are derived for different spectral bands, from ultraviolet to infrared, using the simple stellar population modeling based on the E-MILES stellar spectra library. We take into account obscuration due to dust in star forming regions and diffuse interstellar medium. The $g-r$ colour-stellar mass diagram, at z=0.1, presents a defined red sequence, reaching $g-r \simeq 0.8$, 0.05 dex redder than EAGLE at high masses, and a well populated blue cloud, when field galaxies are included. The clusters' inner regions are dominated by red-sequence galaxies at all masses, although a non-negligible amount of blue galaxies are still present. We adopt Bayesian inference to compute the clusters LFs, testing for statistical significance of both single and double Schechter functions. The multicolour LFs at z=0 show a knee luminosity that peaks in the infrared and increases with the cluster's mass. The faint-end is weakly dependent on colour and mass and shows an upturn in the optical, bounded between -1.25 and -1.39, just moderately steeper than the field. The simulations reproduce, within the observational errors, the spectroscopic LFs of the Hercules and Abell 85 clusters, including their faint end upturn. C-EAGLE LFs are in broad agreement with observed LFs taken from SDSS and XXL surveys, up to z=0.67, showing a rather flat faint end when the observational constrains are taken into account. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: 18 pages, 10 figures, 7 tables. Accepted for publication in MNRAS

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

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

VERTICO: The Virgo Environment Traced In CO Survey

Authors: Toby Brown, Christine D. Wilson, Nikki Zabel, Timothy A. Davis, Alessandro Boselli, Aeree Chung, Sara L. Ellison, Claudia D. P. Lagos, Adam R. H. Stevens, Luca Cortese, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Claire R. Cashmore, Barbara Catinella, Ryan Chown, Benedikt Diemer, Pascal J. Elahi, Maan H. Hani, María J. Jiménez-Donaire, Bumhyun Lee, Katya Leidig, Angus Mok, Karen Pardos Olsen, Laura C. Parker , et al. (11 additional authors not shown)

Abstract: We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star forma… ▽ More We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

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

Simulating Groups and the IntraGroup Medium: The Surprisingly Complex and Rich Middle Ground Between Clusters and Galaxies

Authors: Benjamin D. Oppenheimer, Arif Babul, Yannick Bahé, Iryna S. Butsky, Ian G. McCarthy

Abstract: Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories.… ▽ More Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories. Widely used cosmological simulations with ~(100 Mpc)^3 volumes contain statistical samples of groups that provide important tests of galaxy evolution influenced by environmental processes. Larger volumes capable of reproducing LSS while following the redistribution of baryons by cooling and feedback are essential tools necessary to constrain cosmological parameters. Higher resolution simulations can currently model satellite interactions, the processing of cool (T~10^4 K) multi-phase gas, and non-thermal physics including turbulence, magnetic fields, and cosmic ray transport. We review simulation results regarding the gas and stellar contents of groups, cooling flows and the relation to the central galaxy, the formation and processing of multi-phase gas, satellite interactions with the intragroup medium, and the impact of groups for cosmological parameter estimation. Cosmological simulations provide evolutionarily consistent predictions of these observationally difficult-to-define objects, and have untapped potential to accurately model their gaseous, stellar, and dark matter distributions. △ Less

Submitted 24 June, 2021; originally announced June 2021.

Comments: 61 pages, 18 figures, and 2 tables. This review article is part of the special issue "The Physical Properties of the Groups of Galaxies", edited by L. Lovisari and S. Ettori. Published in MDPI - Universe: https://www.mdpi.com/journal/universe/special_issues/PPGG

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

Strongly lensed cluster substructures are not in tension with $Λ$CDM

Authors: Yannick M. Bahé

Abstract: Strong gravitational lensing observations can test structure formation models by constraining the masses and concentrations of subhaloes in massive galaxy clusters. Recent work has concluded that cluster subhaloes are more abundant and/or concentrated than predicted by $Λ$CDM simulations; this finding has been interpreted as arising from unidentified issues with simulations or an incorrect underst… ▽ More Strong gravitational lensing observations can test structure formation models by constraining the masses and concentrations of subhaloes in massive galaxy clusters. Recent work has concluded that cluster subhaloes are more abundant and/or concentrated than predicted by $Λ$CDM simulations; this finding has been interpreted as arising from unidentified issues with simulations or an incorrect understanding of the nature of dark matter. We test these hypotheses by comparing observed subhalo masses and maximum circular velocities $v_\mathrm{max}$ to predictions from the high resolution Hydrangea galaxy cluster simulation suite, which is based on the successful EAGLE galaxy formation model. The simulated subhalo mass distribution and mass-$v_\mathrm{max}$ relation agrees well with observations, due to the presence of baryons during tidal stripping. Similar agreement is found for the lower-resolution Illustris-TNG300 simulation. In combination, our results suggest that the abundance and concentration of cluster substructures are not in tension with $Λ$CDM, but may provide useful constraints for the refinement of baryon physics models in simulations. △ Less

Submitted 8 March, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 6 pages, 4 figures; minor changes following first referee report (conclusions unchanged)

Too dense to go through: The importance of low-mass clusters for satellite quenching

Authors: Diego Pallero, Facundo A. Gómez, Nelson D. Padilla, Yannick M. Bahé, Cristian A. Vega-Martínez, S. Torres-Flores

Abstract: We study the evolution of satellite galaxies in clusters of the \textsc{c-eagle} simulations, a suite of 30 high-resolution cosmological hydrodynamical zoom-in simulations based on the \textsc{eagle} code. We find that the majority of galaxies that are quenched at $z=0$ ($\gtrsim$ 80$\%$) reached this state in a dense environment (log$_{10}$M$_{200}$[M$_{\odot}$]$\geq$13.5). At low redshift, regar… ▽ More We study the evolution of satellite galaxies in clusters of the \textsc{c-eagle} simulations, a suite of 30 high-resolution cosmological hydrodynamical zoom-in simulations based on the \textsc{eagle} code. We find that the majority of galaxies that are quenched at $z=0$ ($\gtrsim$ 80$\%$) reached this state in a dense environment (log$_{10}$M$_{200}$[M$_{\odot}$]$\geq$13.5). At low redshift, regardless of the final cluster mass, galaxies appear to reach their quenching state in low-mass clusters. Moreover, galaxies quenched inside the cluster that they reside in at $z=0$ are the dominant population in low-mass clusters, while galaxies quenched in a different halo dominate in the most massive clusters. When looking at clusters at $z>0.5$, their in situ quenched population dominates at all cluster masses. This suggests that galaxies are quenched inside the first cluster they fall into. After galaxies cross the cluster's $r_{200}$ they rapidly become quenched ($\lesssim$ 1Gyr). Just a small fraction of galaxies ($\lesssim 15\%$) is capable of retaining their gas for a longer period of time, but after 4Gyr, almost all galaxies are quenched. This phenomenon is related to ram pressure stripping and is produced when the density of the intracluster medium reaches a threshold of $ρ_{\rm ICM}$ $\sim 3 \times 10 ^{-5}$ n$_{\rm H}$ (cm$^{-3}$). These results suggest that galaxies start a rapid-quenching phase shortly after their first infall inside $r_{200}$ and that, by the time they reach $r_{500}$, most of them are already quenched. △ Less

Submitted 3 October, 2023; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 14 pages, 8 figures, 2 appendices

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 511, Issue 3, pp.3210-3227, April 2022

The mass-size relation of galaxy clusters

Authors: O. Contigiani, Y. M. Bahé, H. Hoekstra

Abstract: The outskirts of accreting dark matter haloes exhibit a sudden drop in density delimiting their multi-stream region. Due to the dynamics of accretion, the location of this physically motivated edge strongly correlates with the halo growth rate. Using hydrodynamical zoom-in simulations of high-mass clusters, we explore this definition in realistic simulations and find an explicit connection between… ▽ More The outskirts of accreting dark matter haloes exhibit a sudden drop in density delimiting their multi-stream region. Due to the dynamics of accretion, the location of this physically motivated edge strongly correlates with the halo growth rate. Using hydrodynamical zoom-in simulations of high-mass clusters, we explore this definition in realistic simulations and find an explicit connection between this feature in the dark matter and galaxy profiles. We also show that the depth of the splashback feature correlates well with the direction of filaments and, surprisingly, the orientation of the brightest cluster galaxy. Our findings suggest that galaxy profiles and weak-lensing masses can define an observationally viable mass-size scaling relation for galaxy clusters, which can be used to extract cosmological information. △ Less

Submitted 17 May, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: 9 pages, 9 figures. Fig. 6 and 9 updated with correct labels

Dynamical masses of brightest cluster galaxies II: constraints on the stellar IMF

Authors: S. I. Loubser, H. Hoekstra, A. Babul, Y. M. Bahé, M. Donahue

Abstract: We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30, to place constraints on their stellar Initial Mass Function (IMF). We measure the spatially-resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($Υ_{\star \rm POP}$). We find… ▽ More We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30, to place constraints on their stellar Initial Mass Function (IMF). We measure the spatially-resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($Υ_{\star \rm POP}$). We find young stellar populations ($<$200 Myr) in the centres of 22 per cent of the sample, and constant $Υ_{\star \rm POP}$ within 15 kpc for 60 per cent of the sample. We further use the stellar mass-to-light ratio from the dynamical mass profiles of the BCGs ($Υ_{\star \rm DYN}$), modelled using a Multi-Gaussian Expansion (MGE) and Jeans Anisotropic Method (JAM), with the dark matter contribution explicitly constrained from weak gravitational lensing measurements. We directly compare the stellar mass-to-light ratios derived from the two independent methods, $Υ_{\star \rm POP}$ (assuming some IMF) to $Υ_{\star \rm DYN}$ for the subsample of BCGs with no young stellar populations and constant $Υ_{\star \rm POP}$. We find that for the majority of these BCGs, a Salpeter (or even more bottom-heavy) IMF is needed to reconcile the stellar population and dynamical modelling results although for a small number of BCGs, a Kroupa (or even lighter) IMF is preferred. For those BCGs better fit with a Salpeter IMF, we find that the mass-excess factor against velocity dispersion falls on an extrapolation (towards higher masses) of known literature correlations. We conclude that there is substantial scatter in the IMF amongst the highest-mass galaxies. △ Less

Submitted 10 November, 2020; originally announced November 2020.

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

The stellar mass function and evolution of the density profile of galaxy clusters from the Hydrangea simulations at $0<z<1.5$

Authors: Syeda Lammim Ahad, Yannick M. Bahé, Henk Hoekstra, Remco F. J. van der Burg, Adam Muzzin

Abstract: Galaxy clusters are excellent probes to study the effect of environment on galaxy formation and evolution. Along with high-quality observational data, accurate cosmological simulations are required to improve our understanding of galaxy evolution in these systems. In this work, we compare state-of-the-art observational data of massive galaxy clusters ($>10^{14} \textrm{M}_{\odot}$) at different re… ▽ More Galaxy clusters are excellent probes to study the effect of environment on galaxy formation and evolution. Along with high-quality observational data, accurate cosmological simulations are required to improve our understanding of galaxy evolution in these systems. In this work, we compare state-of-the-art observational data of massive galaxy clusters ($>10^{14} \textrm{M}_{\odot}$) at different redshifts ($0<z<1.5$) with predictions from the Hydrangea suite of cosmological hydrodynamic simulations of 24 massive galaxy clusters ($>10^{14} \textrm{M}_{\odot}$ at $z=0$). We compare three fundamental observables of galaxy clusters: the total stellar mass to halo mass ratio, the stellar mass function (SMF), and the radial mass density profile of the cluster galaxies. In the first two of these, the simulations agree well with the observations, albeit with a slightly too high abundance of $M_\star \lesssim 10^{10} \textrm{M}_{\odot}$ galaxies at $z \gtrsim 1$. The NFW concentrations of cluster galaxies increase with redshift, in contrast to the decreasing dark matter halo concentrations. This previously observed behaviour is therefore due to a qualitatively different assembly of the smooth DM halo compared to the satellite population. Quantitatively, we however find a discrepancy in that the simulations predict higher stellar concentrations than observed at lower redshifts ($z<0.3$), by a factor of $\approx$2. This may be due to selection bias in the simulations, or stem from shortcomings in the build-up and stripping of their inner satellite halo. △ Less

Submitted 12 April, 2021; v1 submitted 30 October, 2020; originally announced October 2020.

Comments: 14 pages, 9 figures (excluding appendices), Accepted for publication in MNRAS

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

A homogeneous measurement of the delay between the onsets of gas stripping and star formation quenching in satellite galaxies of groups and clusters

Authors: Kyle A. Oman, Yannick M. Bahé, Julia Healy, Kelley M. Hess, Michael J. Hudson, Marc A. W. Verheijen

Abstract: We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster environment on star formation in satellite galaxies. We also consider a model for gas stripping, using the same input supplemented with HI fluxes from the ALFALFA survey. The models are motivated by and tested on the… ▽ More We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster environment on star formation in satellite galaxies. We also consider a model for gas stripping, using the same input supplemented with HI fluxes from the ALFALFA survey. The models are motivated by and tested on the Hydrangea cosmological hydrodynamical simulation suite. We recover the characteristic times when satellite galaxies are stripped and quenched. Stripping in massive ($M_\mathrm{ vir}\sim 10^{14.5}\,\mathrm{M}_\odot$) clusters typically occurs at or just before the first pericentric passage. Lower mass ($\sim10^{13.5}\,\mathrm{M}_\odot$) groups strip their satellites on a significantly longer (by $\sim3\,\mathrm{Gyr}$) timescale. Quenching occurs later: Balmer emission lines typically fade $\sim3.5\,\mathrm{Gyr}$ ($5.5\,\mathrm{Gyr}$) after first pericentre in clusters (groups), followed a few hundred $\mathrm{Myr}$ later by reddenning in $(g-r)$ colour. These `delay timescales' are remarkably constant across the entire satellite stellar mass range probed ($\sim10^{9.5}-10^{11}\,\mathrm{M}_\odot$), a feature closely tied to our treatment of `group pre-processing'. The lowest mass groups in our sample ($\sim10^{12.5}\,\mathrm{M}_\odot$) strip and quench their satellites extremely inefficiently: typical timescales may approach the age of the Universe. Our measurements are qualitatively consistent with the `delayed-then-rapid' quenching scenario advocated for by several other studies, but we find significantly longer delay times. Our combination of a homogeneous analysis and input catalogues yields new insight into the sequence of events leading to quenching across wide intervals in host and satellite mass. △ Less

Submitted 7 June, 2024; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: MNRAS accepted version. 25 pages, 11 figures. Appendices included in the source tarball (Download: Other formats: source)

Dynamical masses of brightest cluster galaxies I: stellar velocity anisotropy and mass-to-light ratios

Authors: S. I. Loubser, A. Babul, H. Hoekstra, Y. M. Bahé, E. O'Sullivan, M. Donahue

Abstract: We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for… ▽ More We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for $h_{4}$ for all the BCGs. We derive the stellar mass-to-light ratio ($Υ_{\star \rm DYN}$), and velocity anisotropy ($β$) based on a Multi-Gaussian Expansion (MGE) and axisymmetric Jeans Anisotropic Methods (JAM, cylindrically- and spherically-aligned). We explicitly include a dark matter halo mass component, which is constrained by weak gravitational lensing measurements for these clusters. We find a strong correlation between anisotropy and velocity dispersion profile slope, with rising velocity dispersion profiles corresponding to tangential anisotropy and decreasing velocity dispersion profiles corresponding to radial anisotropy. The rising velocity dispersion profiles can also indicate a significant contribution from the intracluster light (ICL) to the total light (in projection) in the centre of the galaxy. For a small number of BCGs with rising velocity dispersion profiles, a variable stellar mass-to-light ratio can also account for the profile shape, instead of tangential anisotropy or a significant ICL contribution. We note that, for some BCGs, a variable $β_{z}(r)$ (from radial to tangential anisotropy) can improve the model fit to the observed kinematic profiles. The observed diversity in these properties illustrates that BCGs are not the homogeneous class of objects they are often assumed to be. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: 25 pages, accepted for publication in MNRAS

Redshift evolution of the hot intracluster gas metallicity in the C-EAGLE cluster simulations

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Yannick M. Bahe, Richard G. Bower

Abstract: The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investi… ▽ More The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, submitted to MNRAS

The intra-cluster light as a tracer of the total matter density distribution: a view from simulations

Authors: Isaac Alonso Asensio, Claudio Dalla Vecchia, Yannick M. Bahé, David J. Barnes, Scott T. Kay

Abstract: By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-den… ▽ More By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-density contours using the procedure of Montes \& Trujillo (2019), and find that the shape of the stellar mass distribution follows that of the total matter even more closely than observed, although their radial profiles differ substantially. The ratio between stellar and total matter density profiles in circular apertures, shows a slope close to -1, with a small dependence on the cluster's total mass. We propose an indirect method to calculate the halo mass and mass density profile from the radial profile of the intra-cluster stellar mass density. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 7 pages, 5 figures, submitted to MNRAS

CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

Authors: Ricardo Herbonnet, Cristóbal Sifón, Henk Hoekstra, Yannick Bahé, Remco F. J. van der Burg, Jean-Baptiste Melin, Anja von der Linden, David Sand, Scott Kay, David Barnes

Abstract: Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and r… ▽ More Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$σ$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 20 pages, 9 figures, 4 tables, submitted to MNRAS

Constraining the inner density slope of massive galaxy clusters

Authors: Qiuhan He, Hongyu Li, Ran Li, Carlos S. Frenk, Matthieu Schaller, David Barnes, Yannick Bahé, Scott T. Kay, Liang Gao, Claudio Dalla Vecchia

Abstract: We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from… ▽ More We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, $γ_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of $γ_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of $γ_{\rm gNFW}$. △ Less

Submitted 21 June, 2020; v1 submitted 2 July, 2019; originally announced July 2019.

Comments: 19 pages, 16 figures, accepted by MNRAS

Disruption of satellite galaxies in simulated groups and clusters: the roles of accretion time, baryons, and pre-processing

Authors: Yannick M. Bahé, Joop Schaye, David J. Barnes, Claudio Dalla Vecchia, Scott T. Kay, Richard G. Bower, Henk Hoekstra, Sean L. McGee, Tom Theuns

Abstract: We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters t… ▽ More We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes. △ Less

Submitted 1 February, 2019; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 26 pages, 24 figures (16 pages, 14 figures without appendices); accepted by MNRAS (a few extra references added)

Galaxies with monstrous black holes in galaxy cluster environments

Authors: Lieke A. C. van Son, Christopher Barber, Yannick M. Bahe, Joop Schaye, David J. Barnes, Robert A. Crain, Scott T. Kay, Tom Theuns, Claudio Dalla Vecchia

Abstract: Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations… ▽ More Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations, designed to enable the study of high-mass galaxy formation and evolution in cluster environments. We find 69 $M_{\rm BH}(M_{\star})$ outliers at $z=0$, defined as those with $ \rm M_{BH} >10^{7} M_{\odot}$ and $\rm M_{BH}/\rm M_{\star}> 0.01$. This paper focusses on a sample of 5 extreme outliers, that have been selected based on their $\rm M_{BH}$ and $\rm M_{\star}$ values, which are comparable to the most recent estimates of observed positive outliers. This sample of 5 outliers, classified as `Black hole monster galaxies' (BMGs), was traced back in time to study their origin and evolution. In agreement with the results of previous simulations for lower-mass $\rm M_{BH}(\rm M_{\star})$ outliers, we find that these galaxies became outliers due to a combination of their early formation times and tidal stripping. For BMGs with $\rm M_{BH} > 10^9 M_{\odot}$, major mergers (with a stellar mass ratio of $μ> 0.25$) at early times ($z>2$) precede the rapid growth of their supermassive BHs. Furthermore, the scatter in the relation between $\rm M_{BH}$ and stellar velocity dispersion, $σ$, correlates positively with the scatter in [Mg/Fe]($σ$). This indicates that the alpha enhancement of these galaxies, which is closely related to their star formation history, is related to the growth of their central BHs. △ Less

Submitted 10 January, 2019; originally announced January 2019.

Journal ref: Monthly Notices of the Royal Astronomical Society 485 (2019) 396 - 407

Galaxy And Mass Assembly (GAMA): Environmental Quenching of Centrals and Satellites in Groups

Authors: L. J. M. Davies, A. S. G. Robotham, C. del P. Lagos, S. P. Driver, A. R. H. Stevens, Y. M. Bahé, M. Alpaslan, M. N. Bremer, M. J. I. Brown, S. Brough, J. Bland-Hawthorn, L. Cortese, P. Elahi, M. W. Grootes, B. W. Holwerda, A. D. Ludlow, S. McGee, M. Owers, S. Phillipps

Abstract: Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as str… ▽ More Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method. △ Less

Submitted 6 January, 2019; originally announced January 2019.

Comments: 15 pages, 8 figures, MNRAS accepted

The SAMI Galaxy Survey: comparing 3D spectroscopic observations with galaxies from cosmological hydrodynamical simulations

Authors: Jesse van de Sande, Claudia D. P. Lagos, Charlotte Welker, Joss Bland-Hawthorn, Felix Schulze, Rhea-Silvia Remus, Yannick Bahe, Sarah Brough, Julia J. Bryant, Luca Cortese, Scott M. Croom, Julien Devriendt, Yohan Dubois, Michael Goodwin, Iraklis S. Konstantopoulos, Jon S. Lawrence, Anne M. Medling, Christophe Pichon, Samuel N. Richards, Sebastian F. Sanchez, Nicholas Scott, Sarah M. Sweet

Abstract: Cosmological hydrodynamical simulations are rich tools to understand the build-up of stellar mass and angular momentum in galaxies, but require some level of calibration to observations. We compare predictions at $z\sim0$ from the Eagle, Hydrangea, Horizon-AGN, and Magneticum simulations with integral field spectroscopic (IFS) data from the SAMI Galaxy Survey, ATLAS3D, CALIFA and MASSIVE surveys.… ▽ More Cosmological hydrodynamical simulations are rich tools to understand the build-up of stellar mass and angular momentum in galaxies, but require some level of calibration to observations. We compare predictions at $z\sim0$ from the Eagle, Hydrangea, Horizon-AGN, and Magneticum simulations with integral field spectroscopic (IFS) data from the SAMI Galaxy Survey, ATLAS3D, CALIFA and MASSIVE surveys. The main goal of this work is to simultaneously compare structural, dynamical, and stellar population measurements in order to identify key areas of success and tension. We have taken great care to ensure that our simulated measurement methods match the observational methods as closely as possible. We find that the Eagle and Hydrangea simulations reproduce many galaxy relations but with some offsets at high stellar masses. There are moderate mismatches in $R_e$ (+), $ε$ (-), $σ_e$ (-), and mean stellar age (+), where a plus sign indicates that quantities are too high on average, and minus sign too low. The Horizon-AGN simulations qualitatively reproduce several galaxy relations, but there are a number of properties where we find a quantitative offset to observations. Massive galaxies are better matched to observations than galaxies at low and intermediate masses. Overall, we find mismatches in $R_e$ (+), $ε$ (-), $σ_e$ (-) and $(V/σ)_e$ (-). Magneticum matches observations well: this is the only simulation where we find ellipticities typical for disk galaxies, but there are moderate differences in $σ_e$ (-), $(V/σ)_e$ (-) and mean stellar age (+). Our comparison between simulations and observational data has highlighted several areas for improvement, such as the need for improved modelling resulting in a better vertical disk structure, yet our results demonstrate the vast improvement of cosmological simulations in recent years. △ Less

Submitted 12 July, 2019; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 25 pages and 14 figures. Published in the Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, 2019, Volume 484, Issue 1, p.869-891

The signal of decaying dark matter with hydrodynamical simulations

Authors: Mark R. Lovell, David Barnes, Yannick Bahé, Joop Schaye, Matthieu Schaller, Tom Theuns, Sownak Bose, Robert A. Crain, Claudio dalla Vecchia, Carlos S. Frenk, Wojciech Hellwing, Scott T. Kay, Aaron D. Ludlow, Richard G. Bower

Abstract: Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate tha… ▽ More Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 20 pages, 16 Figures. Highlights: Figs 4, 13, 15. To be submitted to MNRAS, comments welcome. Contact: lovell@hi.is