Showing 1–50 of 484 results for author: schaye, J

The MUSE eXtremely Deep Field: Detections of circumgalactic SiII* emission at z>~2

Authors: Haruka Kusakabe, Valentin Mauerhofer, Anne Verhamme, Thibault Garel, Jeremy Blaizot, Lutz Wisotzki, Johan Richard, Leindert A. Boogaard, Floriane Leclercq, Yucheng Guo, Adelaide Claeyssens, Thierry Contini, Edmund Christian Herenz, Josephine Kerutt, Michael V. Maseda, Leo Michel-Dansac, Themiya Nanayakkara, Masami Ouchi, Ismael Pessa, Joop Schaye

Abstract: The circumgalactic medium (CGM) serves as a baryon reservoir that connects galaxies to the intergalactic medium and fuels star formation. The spatial distribution of the metal-enriched cool CGM has not yet been directly revealed at cosmic noon (z~2-4), as bright emission lines at these redshifts are not covered by optical integral field units. To remedy this situation, we aim for the first-ever de… ▽ More The circumgalactic medium (CGM) serves as a baryon reservoir that connects galaxies to the intergalactic medium and fuels star formation. The spatial distribution of the metal-enriched cool CGM has not yet been directly revealed at cosmic noon (z~2-4), as bright emission lines at these redshifts are not covered by optical integral field units. To remedy this situation, we aim for the first-ever detections and exploration of extended SiII* emission (low-ionization state, LIS), referred to as ``SiII* halos'', at redshifts ranging from z=2 to 4 as a means to trace the metal-enriched cool CGM. We use a sample of 39 galaxies with systemic redshifts of z=2.1-3.9 measured with the [CIII] doublet in the MUSE Hubble Ultra Deep Field catalog, which contains integration times spanning from ~30 to 140 hours. We search for extended SiII*1265, 1309, 1533 emission (fluorescent lines) around individual galaxies. We also stack a subsample of 14 UV-bright galaxies. We report five individual detections of SiII*1533 halos. We also confirm the presence of SiII*1533 halos in stacks for the subsample containing UV-bright sources. The other lines do not show secure detections of extended emission in either individual or stacking analyses. These detections may imply that the presence of metal-enriched CGM is a common characteristic for UV-bright galaxies. To investigate whether the origin of SiII* is continuum pumping as suggested in previous studies, we check the consistency of the equivalent width (EW) of SiII* emission and the EW of SiII absorption for the individual halo object with the most reliable detection. We confirm the equivalence, suggesting that photon conservation works for this object and pointing toward continuum pumping as the source of SiII*. We also investigate SiII* lines in a RAMSES-RT zoom-in simulation including continuum pumping and find ubiquitous presence of extended halos. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 13 pages. 8 figures. Submitted to A&A. Comments are welcome!

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

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

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

Submitted 5 June, 2024; originally announced June 2024.

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

Euclid. V. The Flagship galaxy mock catalogue: a comprehensive simulation for the Euclid mission

Authors: Euclid Collaboration, F. J. Castander, P. Fosalba, J. Stadel, D. Potter, J. Carretero, P. Tallada-Crespí, L. Pozzetti, M. Bolzonella, G. A. Mamon, L. Blot, K. Hoffmann, M. Huertas-Company, P. Monaco, E. J. Gonzalez, G. De Lucia, C. Scarlata, M. -A. Breton, L. Linke, C. Viglione, S. -S. Li, Z. Zhai, Z. Baghkhani, K. Pardede, C. Neissner , et al. (344 additional authors not shown)

Abstract: We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from… ▽ More We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from the combination of weak gravitational lensing and galaxy clustering data. The breath of Euclid's data will also foster a wide variety of scientific analyses. The Flagship simulation was developed to provide a realistic approximation to the galaxies that will be observed by Euclid and used in its scientific analyses. We ran a state-of-the-art N-body simulation with four trillion particles, producing a lightcone on the fly. From the dark matter particles, we produced a catalogue of 16 billion haloes in one octant of the sky in the lightcone up to redshift z=3. We then populated these haloes with mock galaxies using a halo occupation distribution and abundance matching approach, calibrating the free parameters of the galaxy mock against observed correlations and other basic galaxy properties. Modelled galaxy properties include luminosity and flux in several bands, redshifts, positions and velocities, spectral energy distributions, shapes and sizes, stellar masses, star formation rates, metallicities, emission line fluxes, and lensing properties. We selected a final sample of 3.4 billion galaxies with a magnitude cut of H_E<26, where we are complete. We have performed a comprehensive set of validation tests to check the similarity to observational data and theoretical models. In particular, our catalogue is able to closely reproduce the main characteristics of the weak lensing and galaxy clustering samples to be used in the mission's main cosmological analysis. (abridged) △ 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

Euclid. III. The NISP Instrument

Authors: Euclid Collaboration, K. Jahnke, W. Gillard, M. Schirmer, A. Ealet, T. Maciaszek, E. Prieto, R. Barbier, C. Bonoli, L. Corcione, S. Dusini, F. Grupp, F. Hormuth, S. Ligori, L. Martin, G. Morgante, C. Padilla, R. Toledo-Moreo, M. Trifoglio, L. Valenziano, R. Bender, F. J. Castander, B. Garilli, P. B. Lilje, H. -W. Rix , et al. (412 additional authors not shown)

Abstract: The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the proc… ▽ More The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated, and its technical potentials and limitations. Links to articles providing more details and technical background are included. NISP's 16 HAWAII-2RG (H2RG) detectors with a plate scale of 0.3" pix^-1 deliver a field-of-view of 0.57deg^2. In photo mode, NISP reaches a limiting magnitude of ~24.5AB mag in three photometric exposures of about 100s exposure time, for point sources and with a signal-to-noise ratio (SNR) of 5. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux ~2x10^-16erg/s/cm^2 integrated over two resolution elements of 13.4A, in 3x560s grism exposures at 1.6 mu (redshifted Ha). Our calibration includes on-ground and in-flight characterisation and monitoring of detector baseline, dark current, non-linearity, and sensitivity, to guarantee a relative photometric accuracy of better than 1.5%, and relative spectrophotometry to better than 0.7%. The wavelength calibration must be better than 5A. NISP is the state-of-the-art instrument in the NIR for all science beyond small areas available from HST and JWST - and an enormous advance due to its combination of field size and high throughput of telescope and instrument. During Euclid's 6-year survey covering 14000 deg^2 of extragalactic sky, NISP will be the backbone for determining distances of more than a billion galaxies. Its NIR data will become a rich reference imaging and spectroscopy data set for the coming decades. △ 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

Euclid. II. The VIS Instrument

Authors: Euclid Collaboration, M. Cropper, A. Al-Bahlawan, J. Amiaux, S. Awan, R. Azzollini, K. Benson, M. Berthe, J. Boucher, E. Bozzo, C. Brockley-Blatt, G. P. Candini, C. Cara, R. A. Chaudery, R. E. Cole, P. Danto, J. Denniston, A. M. Di Giorgio, B. Dryer, J. Endicott, J. -P. Dubois, M. Farina, E. Galli, L. Genolet, J. P. D. Gow , et al. (403 additional authors not shown)

Abstract: This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift ran… ▽ More This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information. △ 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

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Dark photon constraints from CMB temperature anisotropies

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

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

Submitted 8 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

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

MAGAZ3NE: Massive, Extremely Dusty Galaxies at $z\sim2$ Lead to Photometric Overestimation of Number Densities of the Most Massive Galaxies at $3<z<4$

Authors: Ben Forrest, M. C. Cooper, Adam Muzzin, Gillian Wilson, Danilo Marchesini, Ian McConachie, Percy Gomez, Marianna Annunziatella, Z. Cemile Marsan, Joey Braspenning, Wenjun Chang, Gabriella de Lucia, Fabio Fontanot, Michaela Hirschmann, Dylan Nelson, Annalisa Pillepich, Joop Schaye, Stephanie M. Urbano Stawinski, Mauro Stefanon, Lizhi Xie

Abstract: We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distribution… ▽ More We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distributions (SEDs) from the UltraVISTA and VIDEO surveys. In contrast to previous spectroscopic observations of blue star-forming and post-starburst ultramassive galaxies, candidates in this sample have very red SEDs implying significant dust attenuation, old stellar ages, and/or active galactic nuclei (AGN). Of these galaxies, eight are revealed to be heavily dust-obscured $2.0<z<2.7$ galaxies with strong emission lines, some showing broad features indicative of AGN, three are Type I AGN hosts at $z>3$, one is a $z\sim1.2$ dusty galaxy, and four galaxies do not have a confirmed spectroscopic redshift. In fact, none of the sample has |$z_{\rm spec}-z_{\rm phot}$|$<0.5$, suggesting difficulties for photometric redshift programs in fitting similarly red SEDs. The prevalence of these red interloper galaxies suggests that the number densities of high-mass galaxies are overestimated at $z\gtrsim3$ in large photometric surveys, helping to resolve the `impossibly early galaxy problem' and leading to much better agreement with cosmological galaxy simulations. A more complete spectroscopic survey of ultramassive galaxies is required to pin down the uncertainties on their number densities in the early universe. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ. Comments welcome!

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

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

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

Submitted 12 April, 2024; originally announced April 2024.

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

MusE GAs FLOw and Wind (MEGAFLOW) XI. Scaling relations between outflows and host galaxy properties

Authors: Ilane Schroetter, Nicolas F. Bouché, Johannes Zabl, Martin Wendt, Maxime Cherrey, Ivanna Langan, Joop Schaye, Thierry Contini

Abstract: Absorption line spectroscopy using background quasars can provide strong constraints on galactic outflows. In this paper, we investigate possible scaling relations between outflow properties, namely outflow velocity \Vout, the mass ejection rate $\dot M_{\rm out}$, and the mass loading factor $η$ and the host galaxy properties, such as star formation rate (SFR), SFR surface density, redshift, and… ▽ More Absorption line spectroscopy using background quasars can provide strong constraints on galactic outflows. In this paper, we investigate possible scaling relations between outflow properties, namely outflow velocity \Vout, the mass ejection rate $\dot M_{\rm out}$, and the mass loading factor $η$ and the host galaxy properties, such as star formation rate (SFR), SFR surface density, redshift, and stellar mass using galactic outflows probed by background quasars from MEGAFLOW and other surveys. We find that $V_{\rm out}$ ($η$) is (anti-)correlated with SFR and SFR surface density. We extend the formalism of momentum-driven outflows of Heckman et al. to show that it applies not only to down the barrel studies but also to winds probed by background quasars, suggesting a possible universal wind formalism. Under this formalism, we find a clear distinction between ``strong'' and ``weak'' outflows where ``strong'' outflows seem to have tighter correlations with galaxy properties (SFR or galaxy stellar mass) than ``weak'' outflows. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 13 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 687, A39 (2024)

Discovery of optically emitting circumgalactic nebulae around the majority of UV-luminous quasars at intermediate redshift

Authors: Sean D. Johnson, Zhuoqi Will Liu, Jennifer I. Li, Joop Schaye, Jenny E. Greene, Sebastiano Cantalupo, Gwen C. Rudie, Zhijie Qu, Hsiao-Wen Chen, Marc Rafelski, Sowgat Muzahid, Mandy C. Chen, Thierry Contini, Wolfram Kollatschny, Nishant Mishra, Michael Rauch, Patrick Petitjean, Fakhri S. Zahedy

Abstract: We report the discovery of large ionized, [O II] emitting circumgalactic nebulae around the majority of thirty UV luminous quasars at $z=0.4-1.4$ observed with deep, wide-field integral field spectroscopy (IFS) with the Multi-Unit Spectroscopy Explorer (MUSE) by the Cosmic Ultraviolet Baryon Survey (CUBS) and MUSE Quasar Blind Emitters Survey (MUSEQuBES). Among the 30 quasars, seven (23%) exhibit… ▽ More We report the discovery of large ionized, [O II] emitting circumgalactic nebulae around the majority of thirty UV luminous quasars at $z=0.4-1.4$ observed with deep, wide-field integral field spectroscopy (IFS) with the Multi-Unit Spectroscopy Explorer (MUSE) by the Cosmic Ultraviolet Baryon Survey (CUBS) and MUSE Quasar Blind Emitters Survey (MUSEQuBES). Among the 30 quasars, seven (23%) exhibit [O II] emitting nebulae with major axis sizes greater than 100 kpc, twenty greater than 50 kpc (67%), and 27 (90%) greater than 20 kpc. Such large, optically emitting nebulae indicate that cool, dense, and metal-enriched circumgalactic gas is common in the halos of luminous quasars at intermediate redshift. Several of the largest nebulae exhibit morphologies that suggest interaction-related origins. We detect no correlation between the sizes and cosmological dimming corrected surface brightnesses of the nebulae and quasar redshift, luminosity, black hole mass, or radio-loudness, but find a tentative correlation between the nebulae and rest-frame [O II] equivalent width in the quasar spectra. This potential trend suggests a relationship between ISM content and gas reservoirs on CGM scales. The [O II]-emitting nebulae around the $z\approx1$ quasars are smaller and less common than Ly$α$ nebulae around $z\approx3$ quasars. These smaller sizes can be explained if the outer regions of the Ly$α$ halos arise from scattering in more neutral gas, by evolution in the cool CGM content of quasar host halos, by lower-than-expected metallicities on $\gtrsim50$ kpc scales around $z\approx1$ quasars, or by changes in quasar episodic lifetimes between $z=3$ and $1$. △ Less

Submitted 3 April, 2024; v1 submitted 29 March, 2024; originally announced April 2024.

Comments: 18 pages, 5 figures, 2 tables. Accepted for publication in the Astrophysical Journal

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

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

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

Submitted 19 March, 2024; originally announced March 2024.

Comments: 16 pages, 9 figures, submitted to MNRAS

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

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

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

Submitted 18 March, 2024; originally announced March 2024.

Comments: Submitted to MNRAS; 15 pages, 4 figures

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

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

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

Submitted 14 March, 2024; originally announced March 2024.

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

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

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

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

Submitted 12 March, 2024; originally announced March 2024.

Comments: Submitted to AAS Journals

The Cosmic Ultraviolet Baryon Survey (CUBS) VIII: Group Environment of the Most Luminous Quasars at $z\approx1$

Authors: Jennifer I. Li, Sean D. Johnson, Erin Boettcher, Sebastiano Cantalupo, Hsiao-Wen Chen, Mandy C. Chen, David R. DePalma, Zhuoqi, Liu, Nishant Mishra, Patrick Petitjean, Zhijie Qu, Gwen C. Rudie, Joop Schaye, Fakhri S. Zahedy

Abstract: We investigate the group-scale environment of 15 luminous quasars (luminosity $L_{\rm 3000}>10^{46}$ erg s$^{-1}$) from the Cosmic Ultraviolet Baryon Survey (CUBS) at redshift $z\approx1$. Using the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph on the Very Large Telescope (VLT), we conduct a deep galaxy redshift survey in the CUBS quasar fields to identify group members and… ▽ More We investigate the group-scale environment of 15 luminous quasars (luminosity $L_{\rm 3000}>10^{46}$ erg s$^{-1}$) from the Cosmic Ultraviolet Baryon Survey (CUBS) at redshift $z\approx1$. Using the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph on the Very Large Telescope (VLT), we conduct a deep galaxy redshift survey in the CUBS quasar fields to identify group members and measure the physical properties of individual galaxies and galaxy groups. We find that the CUBS quasars reside in diverse environments. The majority (11 out of 15) of the CUBS quasars reside in overdense environments with typical halo masses exceeding $10^{13}{\rm M}_{\odot}$, while the remaining quasars reside in moderate-size galaxy groups. No correlation is observed between overdensity and redshift, black hole (BH) mass, or luminosity. Radio-loud quasars (5 out of 15 CUBS quasars) are more likely to be in overdense environments than their radio-quiet counterparts in the sample, consistent with the mean trends from previous statistical observations and clustering analyses. Nonetheless, we also observe radio-loud quasars in moderate groups and radio-quiet quasars in overdense environments, indicating a large scatter in the connection between radio properties and environment. We find that the most UV luminous quasars might be outliers in the stellar mass-to-halo mass relations or may represent departures from the standard single-epoch BH relations. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 18 pages, 8 figures, accepted for publication in ApJ

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

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

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

Submitted 27 February, 2024; originally announced February 2024.

Comments: 21 pages, 15 figures, comments are welcome

The Cosmic Ultraviolet Baryon Survey (CUBS) VII: on the warm-hot circumgalactic medium probed by O VI and Ne VIII at 0.4 $\lesssim$ z $\lesssim$ 0.7

Authors: Zhijie Qu, Hsiao-Wen Chen, Sean D. Johnson, Gwen C. Rudie, Fakhri S. Zahedy, David DePalma, Joop Schaye, Erin T. Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Claude-André Faucher-Giguère, Jennifer I-Hsiu Li, John S. Mulchaey, Patrick Petitjean, Marc Rafelski

Abstract: This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at… ▽ More This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at $z_{\rm QSO}\gtrsim 0.8$. A total of 30 galaxies or galaxy groups exhibit associated O VI $λλ$ 1031, 1037 doublet absorption within a line-of-sight velocity interval of $\pm250$ km/s, while the rest show no trace of O VI to a detection limit of $\log N_{\rm OVI}/{\rm cm^{-2}}\approx13.7$. Meanwhile, only five galaxies or galaxy groups exhibit the Ne VIII $λλ$ 770,780 doublet absorption, down to a limiting column density of $\log N_{\rm NeVIII}/{\rm cm^{-2}}\approx14.0$. These O VI- and Ne VIII-bearing halos reside in different galaxy environments with stellar masses ranging from $\log M_{\rm star}/M_\odot \approx 8$ to $\approx11.5$. The warm-hot CGM around galaxies of different stellar masses and star formation rates exhibits different spatial profiles and kinematics. In particular, star-forming galaxies with $\log M_{\rm star}/M_\odot\approx9-11$ show a significant concentration of metal-enriched warm-hot CGM within the virial radius, while massive quiescent galaxies exhibit flatter radial profiles of both column densities and covering fractions. In addition, the velocity dispersion of O VI absorption is broad with $σ_v > 40$ km/s for galaxies of $\log M_{\rm star}/M_\odot>9$ within the virial radius, suggesting a more dynamic warm-hot halo around these galaxies. Finally, the warm-hot CGM probed by O VI and Ne VIII is suggested to be the dominant phase in sub-$L^*$ galaxies with $\log M_{\rm star}/M_\odot\approx9-10$ based on their high ionization fractions in the CGM. △ Less

Submitted 14 February, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ after addressing the referee's comments; 28 pages, 16 figures

Lyman Continuum Leaker Candidates at $z\sim3-4$ in the HDUV Based on a Spectroscopic Sample of MUSE LAEs

Authors: J. Kerutt, P. A. Oesch, L. Wisotzki, A. Verhamme, H. Atek, E. C. Herenz, G. D. Illingworth, H. Kusakabe, J. Matthee, V. Mauerhofer, M. Montes, R. P. Naidu, E. Nelson, N. Reddy, J. Schaye, C. Simmonds, T. Urrutia, E. Vitte

Abstract: In recent years, a number of Lyman continuum (LyC) leaker candidates at intermediate redshifts have been found, providing insight into how the Universe was reionised at early cosmic times. Here we identify new LyC leaker candidates at $z\approx 3-4.5$ and compare them to objects from the literature to get an overview of the different observed escape fractions and their relation to the properties o… ▽ More In recent years, a number of Lyman continuum (LyC) leaker candidates at intermediate redshifts have been found, providing insight into how the Universe was reionised at early cosmic times. Here we identify new LyC leaker candidates at $z\approx 3-4.5$ and compare them to objects from the literature to get an overview of the different observed escape fractions and their relation to the properties of the Lyman $α$ (Ly$α$) emission line. The aim of this work is to test indicators for LyC leakage and to improve our understanding of the kind of galaxies from which LyC radiation can escape. We use data from the Hubble Deep Ultraviolet (HDUV) legacy survey to search for LyC emission based on a sample of $\approx 2000$ Ly$α$ emitters (LAEs) detected previously in two surveys with the Multi-Unit Spectroscopic Explorer (MUSE), MUSE-Deep and MUSE-Wide. Based on their known redshifts and positions, we look for potential LyC leakage in the WFC3/UVIS F336W band of the HDUV. The escape fractions are measured and compared based on spectral energy distribution (SED) fitting performed using the CIGALE software. We add twelve objects to the sample of known LyC leaker candidates, one of which was previously known, and compare their Ly$α$ properties to their escape fractions. We find escape fractions between $\sim 20\%$ and $\sim 90\%$, assuming a high transmission in the intergalactic medium (IGM). We show a method to use the number of LyC leaker candidates we find to infer the underlying average escape fraction of galaxies, which is $\approx 12\%$. Based on their Ly$α$ properties, we conclude that LyC leakers are not very different from other high-z LAEs and suggest that most LAEs could be leaking LyC even if this can not always be detected due to the direction of emission and the transmission properties of the IGM. △ Less

Submitted 14 December, 2023; originally announced December 2023.

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

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

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

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

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

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

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

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

Abstract: Galaxy clusters are important probes for both cosmology and galaxy formation physics. We test the cosmological, hydrodynamical FLAMINGO simulations by comparing to observations of the gaseous properties of clusters measured from X-ray observations. FLAMINGO contains unprecedented numbers of massive galaxy groups ($>10^6$) and clusters ($>10^5$) and includes variations in both cosmology and galaxy… ▽ More Galaxy clusters are important probes for both cosmology and galaxy formation physics. We test the cosmological, hydrodynamical FLAMINGO simulations by comparing to observations of the gaseous properties of clusters measured from X-ray observations. FLAMINGO contains unprecedented numbers of massive galaxy groups ($>10^6$) and clusters ($>10^5$) and includes variations in both cosmology and galaxy formation physics. We predict the evolution of cluster scaling relations as well as radial profiles of the temperature, density, pressure, entropy, and metallicity for different masses and redshifts. We show that the differences between volume-, and X-ray-weighting of particles in the simulations, and between cool-core non cool-core samples, are similar in size as the differences between simulations for which the stellar and AGN feedback has been calibrated to produce significantly different gas fractions. Compared to thermally-driven AGN feedback, kinetic jet feedback calibrated to produce the same gas fraction at $R_{\rm 500c}$ yields a hotter core with higher entropies and lower densities, which translates into a smaller fraction of cool-core clusters. Stronger feedback, calibrated to produce lower gas fractions and hence lower gas densities, results in higher temperatures, entropies, and metallicities, but lower pressures. The scaling relations and thermodynamic profiles show almost no evolution with respect to self-similar expectations, except for the metallicity decreasing with redshift. We find that the temperature, density, pressure, and entropy profiles of clusters in the fiducial FLAMINGO simulation are in excellent agreement with observations, while the metallicities in the core are too high. △ Less

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

Comments: 21 pages, 14 figures (including appendix). Accepted for publication in MNRAS. Moved Section 3 to Appendix A, new Section 3 compares scaling relations for different physics models. Added lower mass bin for physics model comparison of thermodynamic profiles. Conclusions unchanged

Bipolar Outflows out to 10~kpc for Massive Galaxies at Redshift $z\approx 1$

Authors: Yucheng Guo, Roland Bacon, Nicolas F. Bouché, Lutz Wisotzki, Joop Schaye, Jérémy Blaizot, Anne Verhamme, Sebastiano Cantalupo, Leindert A. Boogaard, Jarle Brinchmann, Maxime Cherrey, Haruka Kusakabe, Ivanna Langan, Floriane Leclercq, Jorryt Matthee, Léo Michel-Dansac, Ilane Schroetter, Martin Wendt

Abstract: Galactic outflows are believed to play a critical role in the evolution of galaxies by regulating their mass build-up and star formation. Theoretical models assumes bipolar shapes for the outflows that extends well into the circumgalctic medium (CGM), up to tens of kpc perpendicular to the galaxies. They have been directly observed in the local Universe in several individual galaxies, e.g., around… ▽ More Galactic outflows are believed to play a critical role in the evolution of galaxies by regulating their mass build-up and star formation. Theoretical models assumes bipolar shapes for the outflows that extends well into the circumgalctic medium (CGM), up to tens of kpc perpendicular to the galaxies. They have been directly observed in the local Universe in several individual galaxies, e.g., around the Milky Way and M82. At higher redshifts, cosmological simulations of galaxy formation predict an increase in the frequency and efficiency of galactic outflows due to the increasing star formation activity. Outflows are responsible for removing potential fuel for star formation from the galaxy, while at the same enriching the CGM and the intergalactic medium. These feedback processes, although incorporated as key elements of cosmological simulations, are still poorly constrained on CGM scales. Here we present an ultra-deep MUSE image of the mean MgII emission surrounding a sample of galaxies at z~1 that strongly suggests the presence of outflowing gas on physical scales of more than 10kpc. We find a strong dependence of the detected signal on the inclination of the central galaxy, with edge-on galaxies clearly showing enhanced MgII emission along the minor axis, while face-on galaxies display much weaker and more isotropic emission. We interpret these findings as supporting the idea that outflows typically have a bipolar cone geometry perpendicular to the galactic disk. We demonstrate that the signal is not dominated by a few outliers. After dividing the galaxy sample in subsamples by mass, the bipolar emission is only detected in galaxies with stellar mass $\mathrm{M_* \gtrsim 10^{9.5} M_\odot}$. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: published in Nature. 25 pages, 10 figures

Journal ref: 2023, Nature, 624, 53

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

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

Abstract: The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev-Zel'dovich (SZ) observations… ▽ More The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev-Zel'dovich (SZ) observations of the hot intracluster gas. Here, we use the FLAMINGO simulations to investigate whether an extremum gradient in a similar position to the splashback radius is predicted to occur in the cluster gas profiles. We find that the minimum in the gradient of the stacked 3D gas density and pressure profiles, and the maximum in the gradient of the entropy profile, broadly align with the splashback feature though there are significant differences. While the dark matter splashback radius varies with specific mass accretion rate, in agreement with previous work, the radial position of the deepest minimum in the log-slope of the gas density is more sensitive to halo mass. In addition, we show that a similar minimum is also present in projected 2D pseudo-observable profiles: emission measure (X-ray); Compton-$y$ (SZ) and surface mass density (weak lensing). We find that the latter traces the dark matter results reasonably well albeit the minimum occurs at a slightly smaller radius. While results for the gas profiles are largely insensitive to accretion rate and various observable proxies for dynamical state, they do depend on the strength of the feedback processes. △ Less

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

Comments: Accepted by MNRAS

MusE GAs FLOw and Wind (MEGAFLOW) X. The cool gas and covering fraction of MgII in galaxy groups

Authors: Maxime Cherrey, Nicolas Bouché, Johannes Zabl, Ilane Schroetter, Martin Wendt, Ivanna Langan, Johan Richard, Joop Schaye, Wilfried Mercier, Benoit Epinat, Thierry Contini

Abstract: We present a study of the cool gas ($\approx 10^4$ K) traced by MgII absorptions around groups of galaxies in the MEGAFLOW survey. Using a combination of two algorithms we blindly identify 32 groups of more than 5 galaxies at $0.3 < z < 1.5$ with $10.7 < \log_{10}(M/\rm M_{\odot}) < 13.7$. Among them 26 can be used to study potential counterpart MgII absorptions. We report that 21 out of the total… ▽ More We present a study of the cool gas ($\approx 10^4$ K) traced by MgII absorptions around groups of galaxies in the MEGAFLOW survey. Using a combination of two algorithms we blindly identify 32 groups of more than 5 galaxies at $0.3 < z < 1.5$ with $10.7 < \log_{10}(M/\rm M_{\odot}) < 13.7$. Among them 26 can be used to study potential counterpart MgII absorptions. We report that 21 out of the total 120 MgII absorption systems present in MEGAFLOW are associated with groups. We observe that the MgII rest-frame equivalent width ($W^{2796}_r$) drops at an impact parameter of $\approx 150$ projected kpc from the closest galaxy and $\approx$ one virial radius from the identified group center indicating that MgII halos scale with the mass of the groups.The impact parameter where the covering fraction exceeds $50\%$ is $\log_{10}(b/\rm kpc) = 2.17 \pm 0.47$ $(2 σ)$ and $(b/R_{\rm vir}) = 1.67 \pm 0.98$, which is $\approx 3$ times larger than for field galaxies ($\log_{10}(b/\rm kpc)=1.67\pm0.15$). Finally, we estimate the cool gas column density profile in groups (from the $W^{2796}_r$) and show that its shape follows closely the typical dark matter column density profile for halos at similar redshift and masses. △ Less

Submitted 4 December, 2023; originally announced December 2023.

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

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

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

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

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

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

Euclid preparation. XXXIX. The effect of baryons on the Halo Mass Function

Authors: Euclid Collaboration, T. Castro, S. Borgani, M. Costanzi, J. Dakin, K. Dolag, A. Fumagalli, A. Ragagnin, A. Saro, A. M. C. Le Brun, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero , et al. (198 additional authors not shown)

Abstract: The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being sub-dominant to dark matter and dark energy, the baryonic component in our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model to precisely quantify the… ▽ More The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being sub-dominant to dark matter and dark energy, the baryonic component in our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model to precisely quantify the impact of baryons on galaxy cluster virial halo masses, using the baryon fraction within a cluster as proxy for their effect. Constructed on the premise of quasi-adiabaticity, the model includes two parameters calibrated using non-radiative cosmological hydrodynamical simulations and a single large-scale simulation from the Magneticum set, which includes the physical processes driving galaxy formation. As a main result of our analysis, we demonstrate that this model delivers a remarkable one percent relative accuracy in determining the virial dark matter-only equivalent mass of galaxy clusters, starting from the corresponding total cluster mass and baryon fraction measured in hydrodynamical simulations. Furthermore, we demonstrate that this result is robust against changes in cosmological parameters and against varying the numerical implementation of the sub-resolution physical processes included in the simulations. Our work substantiates previous claims about the impact of baryons on cluster cosmology studies. In particular, we show how neglecting these effects would lead to biased cosmological constraints for a Euclid-like cluster abundance analysis. Importantly, we demonstrate that uncertainties associated with our model, arising from baryonic corrections to cluster masses, are sub-dominant when compared to the precision with which mass-observable relations will be calibrated using Euclid, as well as our current understanding of the baryon fraction within galaxy clusters. △ Less

Submitted 16 April, 2024; v1 submitted 25 October, 2023; originally announced November 2023.

Comments: 18 pages, 10 figures, 4 tables, 1 appendix, abstract abridged for arXiv submission; v2 matches published version

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

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

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

Submitted 31 October, 2023; originally announced November 2023.

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

An ensemble study of turbulence in extended QSO nebulae at $z\approx0.5$--1

Authors: Mandy C. Chen, Hsiao-Wen Chen, Michael Rauch, Zhijie Qu, Sean D. Johnson, Joop Schaye, Gwen C. Rudie, Jennifer I-Hsiu Li, Zhuoqi, Liu, Fakhri S. Zahedy, Sebastiano Cantalupo, Erin Boettcher

Abstract: Turbulent motions in the circumgalactic medium (CGM) play a critical role in regulating the evolution of galaxies, yet their detailed characterization remains elusive. Using two-dimensional velocity maps constructed from spatially-extended [OII] and [OIII] emission, Chen et al. (2023b) measured the velocity structure functions (VSFs) of four quasar nebulae at $z\approx\!0.5$--1.1. One of these exh… ▽ More Turbulent motions in the circumgalactic medium (CGM) play a critical role in regulating the evolution of galaxies, yet their detailed characterization remains elusive. Using two-dimensional velocity maps constructed from spatially-extended [OII] and [OIII] emission, Chen et al. (2023b) measured the velocity structure functions (VSFs) of four quasar nebulae at $z\approx\!0.5$--1.1. One of these exhibits a spectacular Kolmogorov relation. Here we carry out an ensemble study using an expanded sample incorporating four new nebulae from three additional QSO fields. The VSFs measured for all eight nebulae are best explained by subsonic turbulence revealed by the line-emitting gas, which in turn strongly suggests that the cool gas ($T\!\sim\!10^4$ K) is dynamically coupled to the hot ambient medium. Previous work demonstrates that the largest nebulae in our sample reside in group environments with clear signs of tidal interactions, suggesting that environmental effects are vital in seeding and enhancing turbulence within the gaseous halos, ultimately promoting the formation of the extended nebulae. No discernible differences are observed in the VSF properties between radio-loud and radio-quiet QSO fields. We estimate the turbulent heating rate per unit volume, $Q_{\rm turb}$, in the QSO nebulae to be $\sim 10^{-26}$--$10^{-22}$ erg cm$^{-3}$ s$^{-1}$ for the cool phase and $\sim 10^{-28}$--$10^{-25}$ erg cm$^{-3}$ s$^{-1}$ for the hot phase. This range aligns with measurements in the intracluster medium and star-forming molecular clouds but is $\sim10^3$ times higher than the $Q_{\rm turb}$ observed inside cool gas clumps on scales $\lesssim1$ kpc using absorption-line techniques. We discuss the prospect of bridging the gap between emission and absorption studies by pushing the emission-based VSF measurements to below $\approx\!10$ kpc. △ Less

Submitted 12 January, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 23 pages; 7 figures, and 4 tables in main text; 9 figures in Appendix; accepted by ApJ. Comments welcome

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

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

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

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

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

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

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

An efficient and robust method to estimate halo concentration

Authors: Kai Wang, H. J. Mo, Yangyao Chen, Joop Schaye

Abstract: We propose an efficient and robust method to estimate the halo concentration based on the first moment of the density distribution, which is $R_1\equiv \int_0^{r_{\rm vir}}4πr^3ρ(r)dr/M_{\rm vir}/r_{\rm vir}$. We find that $R_1$ has a monotonic relation with the concentration parameter of the NFW profile, and that a cubic polynomial function can fit the relation with an error $\lesssim 3\%$. Tests… ▽ More We propose an efficient and robust method to estimate the halo concentration based on the first moment of the density distribution, which is $R_1\equiv \int_0^{r_{\rm vir}}4πr^3ρ(r)dr/M_{\rm vir}/r_{\rm vir}$. We find that $R_1$ has a monotonic relation with the concentration parameter of the NFW profile, and that a cubic polynomial function can fit the relation with an error $\lesssim 3\%$. Tests on ideal NFW halos show that the conventional NFW profile fitting method and the $V_{\rm max}/V_{\rm vir}$ method produce biased halo concentration estimation by $\approx 10\%$ and $\approx 30\%$, respectively, for halos with 100 particles. In contrast, the systematic error for our $R_1$ method is smaller than $0.5\%$ even for halos containing only 100 particles. Convergence tests on realistic halos in $N$-body simulations show that the NFW profile fitting method underestimates the concentration parameter for halos with $\lesssim 300$ particles by $\gtrsim 20\%$, while the error for the $R_1$ method is $\lesssim 8\%$. We also show other applications of $R_1$, including estimating $V_{\rm max}$ and the Einasto concentration $c_{\rm e}\equiv r_{\rm vir}/r_{-2}$. The calculation of $R_1$ is efficient and robust, and we recommend including it as one of the halo properties in halo catalogs of cosmological simulations. △ Less

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

Comments: 17 pages, 9 + 9 figures. Main figure: Fig. 3. MNRAS accepted

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

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

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

Submitted 24 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS, 23 pages, 20 figures

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

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

Abstract: A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological parameters fit to the cosmic microwave background (CMB). The origin of this '$S_8$ tension' remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in… ▽ More A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological parameters fit to the cosmic microwave background (CMB). The origin of this '$S_8$ tension' remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in the observational measurements and/or uncertainties in the role that baryons play. Here we carefully examine the latter possibility using the new FLAMINGO suite of large-volume cosmological hydrodynamical simulations. We project the simulations onto observable harmonic space and compare with observational measurements of the power and cross-power spectra of cosmic shear, CMB lensing, and the thermal Sunyaev-Zel'dovich (tSZ) effect. We explore the dependence of the predictions on box size and resolution, cosmological parameters including the neutrino mass, and the efficiency and nature of baryonic 'feedback'. Despite the wide range of astrophysical behaviours simulated, we find that baryonic effects are not sufficiently large to remove the $S_8$ tension. Consistent with recent studies, we find the CMB lensing power spectrum is in excellent agreement with the standard model, whilst the cosmic shear power spectrum, tSZ effect power spectrum, and the cross-spectra between shear, CMB lensing, and the tSZ effect are all in varying degrees of tension with the CMB-specified standard model. These results suggest that some mechanism is required to slow the growth of fluctuations at late times and/or on non-linear scales, but that it is unlikely that baryon physics is driving this modification. △ Less

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

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

Spatially-resolved Spectroscopic Analysis of Ly$α$ Haloes: Radial Evolution of the Ly$α$ Line Profile out to 60 kpc

Authors: Yucheng Guo, Roland Bacon, Lutz Wisotzki, Thibault Garel, Jérémy Blaizot, Joop Schaye, Jorryt Matthee, Floriane Leclercq, Leindert Boogaard, Johan Richard, Anne Verhamme, Jarle Brinchmann, Léo Michel-Dansac, Haruka Kusakabe

Abstract: Deep MUSE observations have unveiled pervasive Ly$α$ haloes (LAHs) surrounding high-redshift star-forming galaxies. However, the origin of the extended Ly$α$ emission is still a subject of debate. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Ly$α$ emitters (LAEs) at redshift $3<z<4$ in the MUSE Extremely Deep Field (MXDF). With… ▽ More Deep MUSE observations have unveiled pervasive Ly$α$ haloes (LAHs) surrounding high-redshift star-forming galaxies. However, the origin of the extended Ly$α$ emission is still a subject of debate. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Ly$α$ emitters (LAEs) at redshift $3<z<4$ in the MUSE Extremely Deep Field (MXDF). With respect to the Ly$α$ red peak of the target LAE, the Ly$α$ line peak becomes increasingly more blueshifted out to a projected distance of at least 60 kpc, where the velocity offset is $\approx$ 250 km/s. This signal is observed in both the mean and median stacks, and is thus a generic property of the LAE sample with typical Ly$α$ luminosity $\mathrm{\approx 10^{41.1} erg\,s^{-1}}$. We discuss multiple scenarios to explain the blueshift of the circumgalactic Ly$α$ line. The most plausible one is a combination of outflows and inflows. In the inner region of the LAH, the Ly$α$ photons are produced by the central star formation and then scattered within outflows. At larger radii, the infalling cool gas shapes the observed Ly$α$ blueshift. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: Submitted to A&A Letters

The Cosmic Ultraviolet Baryon Survey: Empirical Characterization of Turbulence in the Cool Circumgalactic Medium

Authors: Hsiao-Wen Chen, Zhijie Qu, Michael Rauch, Mandy C. Chen, Fakhri S. Zahedy, Sean D. Johnson, Joop Schaye, Gwen C. Rudie, Erin Boettcher, Sebastiano Cantalupo, Claude-Andre Faucher-Giguere, Jenny E. Greene, Sebastian Lopez, Robert A. Simcoe

Abstract: This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower tha… ▽ More This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower than nH = 0.001 cm^-3 are identified in galaxy halos. In addition, comparing the line widths between different elements for kinematically matched components provides robust empirical constraints on the thermal temperature T and the non-thermal motions bNT, independent of the ionization models. On average, bNT is found to increase with l_cl following bNT \propto l_cl^0.3 over three decades in spatial scale from l_cl~1 pc to l_cl~1 kpc. Attributing the observed bNT to turbulent motions internal to the clumps, the best-fit bNT-l_cl relation shows that the turbulence is consistent with Kolmogorov at <1 kpc with a roughly constant energy transfer rate per unit mass of epsilon~0.003 cm^2 s^-3 and a dissipation time scale of <~ 100 Myr. No significant difference is found between massive quiescent and star-forming halos in the sample on scales less than 1 kpc. While the inferred epsilon is comparable to what is found in CIV absorbers at high redshift, it is considerably smaller than observed in star-forming gas or in extended line-emitting nebulae around distant quasars. A brief discussion of possible sources to drive the observed turbulence in the cool CGM is presented. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ Letters; 10 pages, 4 figures

Median Surface Brightness Profiles of Lyman-$α$ Haloes in the MUSE Extremely Deep Field

Authors: Yucheng Guo, Roland Bacon, Lutz Wisotzki, Thibault Garel, Jérémy Blaizot, Joop Schaye, Johan Richard, Yohana Herrero Alonso, Floriane Leclercq, Leindert Boogaard, Haruka Kusakabe, John Pharo

Abstract: We present the median surface brightness profiles of diffuse Ly$α$ haloes (LAHs) around star-forming galaxies by stacking 155 spectroscopically confirmed Ly$α$ emitters (LAEs) at 3<z<4 in the MUSE Extremely Deep Field (MXDF), with median Ly$α$ luminosity $\mathrm{L_{Lyα} \approx 10^{41.1} erg\,s^{-1}}$. After correcting for a systematic surface brightness offset we identified in the datacube, we d… ▽ More We present the median surface brightness profiles of diffuse Ly$α$ haloes (LAHs) around star-forming galaxies by stacking 155 spectroscopically confirmed Ly$α$ emitters (LAEs) at 3<z<4 in the MUSE Extremely Deep Field (MXDF), with median Ly$α$ luminosity $\mathrm{L_{Lyα} \approx 10^{41.1} erg\,s^{-1}}$. After correcting for a systematic surface brightness offset we identified in the datacube, we detect extended Ly$α$ emission out to a distance of 270 kpc. The median Ly$α$ surface brightness profile shows a power-law decrease in the inner 20 kpc, and a possible flattening trend at larger distance. This shape is similar for LAEs with different Ly$α$ luminosities, but the normalisation of the surface brightness profile increases with luminosity. At distances larger than 50 kpc, we observe strong overlap of adjacent LAHs, and the Ly$α$ surface brightness is dominated by the LAHs of nearby LAEs. We find no clear evidence of redshift evolution of the observed Ly$α$ profiles when comparing with samples at 4<z<5 and 5<z<6. Our results are consistent with a scenario in which the inner 20 kpc of the LAH is powered by star formation in the central galaxy, while the LAH beyond a radius of 50 kpc is dominated by photons from surrounding galaxies. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Submitted to A&A

The first comprehensive study of a giant nebula around a radio-quiet quasar in the $z < 1$ Universe

Authors: Zhuoqi Will Liu, Sean D. Johnson, Jennifer I-Hsiu Li, Gwen C. Rudie, Joop Schaye, Hsiao-Wen Chen, Jarle Brinchmann, Sebastiano Cantalupo, Mandy C. Chen, Wolfram Kollatschny, Michael V. Maseda, Nishant Mishra, Sowgat Muzahid

Abstract: We present the first comprehensive study of a giant, $\approx \! \! 70$ kpc-scale nebula around a radio-quiet quasar at $z<1$. The analysis is based on deep integral field spectroscopy with MUSE of the field of HE$\,$0238$-$1904, a luminous quasar at $z=0.6282$. The nebula emits strongly in $\mathrm{[O \, II]}$, $\rm H β$, and $\mathrm{[O \, III]}$, and the quasar resides in an unusually overdense… ▽ More We present the first comprehensive study of a giant, $\approx \! \! 70$ kpc-scale nebula around a radio-quiet quasar at $z<1$. The analysis is based on deep integral field spectroscopy with MUSE of the field of HE$\,$0238$-$1904, a luminous quasar at $z=0.6282$. The nebula emits strongly in $\mathrm{[O \, II]}$, $\rm H β$, and $\mathrm{[O \, III]}$, and the quasar resides in an unusually overdense environment for a radio-quiet system. The environment likely consists of two groups which may be merging, and in total have an estimated dynamical mass of $M_{\rm dyn}\approx 4\times 10^{13}$ to $10^{14}\ {\rm M_\odot}$. The nebula exhibits largely quiescent kinematics and irregular morphology. The nebula may arise primarily through interaction-related stripping of circumgalactic and interstellar medium (CGM/ISM) of group members, with some potential contributions from quasar outflows. The simultaneous presence of the giant nebula and a radio-quiet quasar in a rich environment suggests a correlation between such circum-quasar nebulae and environmental effects. This possibility can be tested with larger samples. The upper limits on the electron number density implied by the $\mathrm{[O \, II]}$ doublet ratio range from $\log(n_{\rm e, \, [O \, II]} / \mathrm{cm^{-3}}) < 1.2$ to $2.8$. However, assuming a constant quasar luminosity and negligible projection effects, the densities implied from the measured line ratios between different ions (e.g., $\mathrm{[O\,II]}$, $\mathrm{[O\,III]}$, and $\mathrm{[Ne\,V]}$) and photoionization simulations are often $10{-}400$ times larger. This large discrepancy can be explained by quasar variability on a timescale of $\approx 10^4{-}10^5$ years. △ Less

Submitted 30 October, 2023; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: 19 pages, 9 figures, 3 tables; Accepted by MNRAS

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

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

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

Submitted 18 August, 2023; originally announced August 2023.

Comments: 14 pages, 21 figures, 1 table

The contribution of faint Lyman-$α$ emitters to extended Lyman-$α$ halos constrained by MUSE clustering measurements

Authors: Yohana Herrero Alonso, L. Wisotzki, T. Miyaji, J. Schaye, J. Pharo, M. Krumpe

Abstract: Detections of extended Ly$α$ halos (LAHs) around Ly$α$ emitters (LAEs) have lately been reported on a regular basis, but their origin is still under investigation. Simulation studies predict that the outer regions of the extended LAHs contain a major contribution from the Ly$α$ emission of faint, individually undetected LAEs. To address this matter from an observational angle, we use halo occupati… ▽ More Detections of extended Ly$α$ halos (LAHs) around Ly$α$ emitters (LAEs) have lately been reported on a regular basis, but their origin is still under investigation. Simulation studies predict that the outer regions of the extended LAHs contain a major contribution from the Ly$α$ emission of faint, individually undetected LAEs. To address this matter from an observational angle, we use halo occupation distribution (HOD) modeling to reproduce the clustering of a spectroscopic sample of 1265 LAEs at $3<z<5$ from the MUSE-Wide survey. We integrate the Ly$α$ luminosity function (LF) to estimate the background surface brightness due to discrete faint LAEs. We then extend the HOD statistics inwards towards small separations and compute the factor by which the measured Ly$α$ surface brightness (SB) is enhanced by undetected close physical neighbors. We consider various clustering scenarios for the undetected sources and compare the corresponding radial profiles. The resulting inferred Ly$α$ SB of faint LAEs ranges between $(0.4-2)\times10^{20}\;\rm{erg}\;\rm{s}^{-1}~\rm{cm}^{-2}~\rm{arcsec}^{-2}$, with a very slow radial decline outwards. Our results suggest that the outer regions of observed LAHs ($R\gtrsim50$~pkpc) could indeed contain a strong component from external (but physically associated) LAEs, possibly even be dominated by them. Only for a relatively shallow faint-end slope of the Ly$α$ LF would this contribution from clustered LAEs become unimportant. We also confirm that the observed emission from the inner regions ($R\le20-30$~pkpc) is too bright to be significantly affected by clustering. We compare our findings with predicted profiles from simulations and find good overall agreement. We outline possible future measurements to further constrain the impact of discrete undetected LAEs on observed extended LAHs. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 13 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 677, A125 (2023)

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

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

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

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

Comments: Accepted for publication in MNRAS

Mapping the imprints of stellar and AGN feedback in the circumgalactic medium with X-ray microcalorimeters

Authors: Gerrit Schellenberger, Ákos Bogdán, John A. ZuHone, Benjamin D. Oppenheimer, Nhut Truong, Ildar Khabibullin, Fred Jennings, Annalisa Pillepich, Joseph Burchett, Christopher Carr, Priyanka Chakraborty, Robert Crain, William Forman, Christine Jones, Caroline A. Kilbourne, Ralph P. Kraft, Maxim Markevitch, Daisuke Nagai, Dylan Nelson, Anna Ogorzalek, Scott Randall, Arnab Sarkar, Joop Schaye, Sylvain Veilleux, Mark Vogelsberger , et al. (2 additional authors not shown)

Abstract: The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future, large grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate… ▽ More The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future, large grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate a wide range of potential measurements, which will address the open questions in galaxy formation and evolution. By including all background and foreground components in our mock observations, we show why it is impossible to perform these measurements with current instruments, such as X-ray CCDs, and only microcalorimeters will allow us to distinguish the faint CGM emission from the bright Milky Way (MW) foreground emission lines. We find that individual halos of MW mass can, on average and depending on star formation rate, be traced out to large radii, around R500, and for larger galaxies even out to R200, using prominent emission lines, such as OVII, or OVIII. Furthermore, we show that emission line ratios for individual halos can reveal the radial temperature structure. Substructure measurements show that it will be possible to relate azimuthal variations to the feedback mode of the galaxy. We demonstrate the ability to construct temperature, velocity, and abundance ratio maps from spectral fitting for individual galaxy halos, which reveal rotation features, AGN outbursts, and enrichment. △ Less

Submitted 29 April, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 41 pages, 19 figures, accepted for publication in ApJ

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

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

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

Submitted 26 June, 2023; originally announced June 2023.

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

The Cosmic Ultraviolet Baryon Survey (CUBS) VI: Connecting Physical Properties of the Cool Circumgalactic Medium to Galaxies at $z\approx 1$

Authors: Zhijie Qu, Hsiao-Wen Chen, Gwen C. Rudie, Sean D. Johnson, Fakhri S. Zahedy, David DePalma, Erin Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jennifer I-Hsiu Li, Sebastian Lopez, Joop Schaye, Robert A. Simcoe

Abstract: This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$σ$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental… ▽ More This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$σ$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental properties of the galaxies, including galaxy overdensities, the total stellar mass and gravitational potential summed over all nearby neighbors, and the presence of local ionizing sources, are found to have a significant impact on the observed CGM absorption properties. Specifically, massive galaxies and galaxies in overdense regions exhibit a higher rate of incidence of absorption. At the same time, the observed CGM absorption properties in galaxy groups appear to be driven by the galaxy closest to the QSO sightline, rather than by the most massive galaxy or by mass-weighted properties. We introduce a total projected gravitational potential $ψ$, defined as $-ψ/G =\sum M_{\rm halo}/d_{\rm proj}$ summed over all group members, to characterize the overall galaxy environment. This projected gravitational potential correlates linearly with the maximum density detected in each sightline, consistent with higher-pressure gas being confined in deeper gravitational potential wells. In addition, we find that the radial profile of cool gas density exhibits a general decline from the inner regions to the outskirts, being in pressure balance with the hot halo. Finally, we note that the ionizing flux from nearby galaxies can generate an elevated $N$(HI)/$N$(HeI) ratio, which in turn provides a unique diagnostic of possible local sources contributing to the ionizing radiation field. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 17 pages, 11 figures and 1 table in the main text; 40 figures and 1 table in the appendix; accepted for publication in MNRAS

Spurious heating of stellar motions by dark matter particles in cosmological simulations of galaxy formation

Authors: Aaron D. Ludlow, S. Michael Fall, Matthew J. Wilkinson, Joop Schaye, Danail Obreschkow

Abstract: We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. Many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution… ▽ More We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. Many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution of their gaseous baryons. The distribution and kinematics of stellar and DM particles, however, are affected. Below a resolution-dependent virial mass, $M_{200}^{\rm spur}$, galaxies have higher characteristic velocities, larger sizes, and more angular momentum in the simulation with lower DM mass resolution; haloes have higher central densities and lower velocity dispersions. Above $M_{200}^{\rm spur}$, galaxies and haloes have similar properties in both runs. The differences arise due to spurious heating, which transfers energy from DM to stellar particles, causing galaxies to heat up and haloes to cool down. The value of $M_{200}^{\rm spur}$ can be derived from an empirical disc heating model, and coincides with the mass below which the predicted {\em spurious} velocity dispersion exceeds the {\em measured} velocity dispersion of simulated galaxies. We predict that galaxies in the $100^3\, {\rm Mpc}^3$ \eagle\, run and IllustrisTNG-100 are robust to spurious collisional effects at their half-mass radii provided $M_{200}^{\rm spur}\approx 10^{11.7}{\rm M_\odot}$; for the $25^3\, {\rm Mpc}^3$ \eagle\, run and IllustrisTNG-50, we predict $M_{200}^{\rm spur}\approx 10^{11}{\rm M_\odot}$. Suppressing spurious heating at smaller/larger radii, or for older/younger stellar populations, requires haloes to be resolved with more/fewer DM particles. △ Less

Submitted 4 September, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 17 pages, 15 figures. MNRAS published. A basic python script for calculating spurious disc heating rates is available at https://github.com/AaronDLudlow/hot-disc.git

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

Large-scale excess HI absorption around $z\approx4$ galaxies detected in a background galaxy spectrum in the MUSE eXtremely Deep Field

Authors: Jorryt Matthee, Christopher Golling, Ruari Mackenzie, Gabriele Pezzulli, Simon Lilly, Joop Schaye, Roland Bacon, Haruka Kusakabe, Tanya Urrutia, Leindert Boogaard, Jarle Brinchmann, Michael V. Maseda, Thibault Garel, Nicolas F. Bouché, Lutz Wisotzki

Abstract: Observationally mapping the relation between galaxies and the intergalactic medium (IGM) is of key interest for studies of cosmic reionization. Diffuse hydrogen gas has typically been observed in HI Lyman-$α$ (Ly$α$) absorption in the spectra of bright background quasars. However, it is important to extend these measurements to background galaxies as quasars become increasingly rare at high redshi… ▽ More Observationally mapping the relation between galaxies and the intergalactic medium (IGM) is of key interest for studies of cosmic reionization. Diffuse hydrogen gas has typically been observed in HI Lyman-$α$ (Ly$α$) absorption in the spectra of bright background quasars. However, it is important to extend these measurements to background galaxies as quasars become increasingly rare at high redshift and rarely probe closely separated sight-lines. Here we use deep integral field spectroscopy in the MUSE eXtremely Deep Field to demonstrate the measurement of the Ly$α$ transmission at $z\approx4$ in absorption to a background galaxy at $z=4.77$. The HI transmission is consistent with independent quasar sight-lines at similar redshifts. Exploiting the high number of spectroscopic redshifts of faint galaxies (500 between $z=4.0-4.7$ within a radius of 8 arcmin) that are tracers of the density field, we show that Ly$α$ transmission is inversely correlated with galaxy density, i.e. transparent regions in the Ly$α$ forest mark under-dense regions at $z\approx4$. Due to large-scale clustering, galaxies are surrounded by excess HI absorption over the cosmic mean out to 4 cMpc/h. We also find that redshifts from the peak of the Ly$α$ line are typically offset from the systemic redshift by +170 km/s. This work extends results from $z\approx 2 - 3$ to higher redshifts and demonstrates the power of deep integral field spectroscopy to simultaneously measure the ionization structure of the IGM and the large-scale density field in the early Universe. △ Less

Submitted 5 March, 2024; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS. Main text 10 pages, 9 figures. Key results in Fig 4 (Lya forest transmission in the MXDF field) and Fig 9 (transmission - galaxy distance cross-correlation)

JWST/NIRSpec Measurements of Extremely Low Metallicities in High Equivalent Width Lyman-$α$ Emitters

Authors: Michael V. Maseda, Zach Lewis, Jorryt Matthee, Joseph F. Hennawi, Leindert Boogaard, Anna Feltre, Themiya Nanayakkara, Roland Bacon, Amy Barger, Jarle Brinchmann, Marijn Franx, Takuya Hashimoto, Hanae Inami, Haruka Kusakabe, Floriane Leclercq, Lucie Rowland, Anthony J. Taylor, Christy Tremonti, Tanya Urrutia, Joop Schaye, Charlotte Simmonds, Eloïse Vitte

Abstract: Deep VLT/MUSE optical integral field spectroscopy has recently revealed an abundant population of ultra-faint galaxies ($M_{UV} \approx -15$; 0.01 $L_{\star}$) at $z=$2.9$-$6.7 due to their strong Lyman-$α$ emission with no detectable continuum. The implied Lyman-$α$ equivalent widths can be in excess of 100-200 Angstrom, challenging existing models of normal star formation and indicating extremel… ▽ More Deep VLT/MUSE optical integral field spectroscopy has recently revealed an abundant population of ultra-faint galaxies ($M_{UV} \approx -15$; 0.01 $L_{\star}$) at $z=$2.9$-$6.7 due to their strong Lyman-$α$ emission with no detectable continuum. The implied Lyman-$α$ equivalent widths can be in excess of 100-200 Angstrom, challenging existing models of normal star formation and indicating extremely young ages, small stellar masses, and a very low amount of metal enrichment. We use JWST/NIRSpec's microshutter array to follow-up 45 of these galaxies (11h in G235M/F170LP and 7h in G395M/F290LP), as well as 45 lower-equivalent width Lyman-$α$ emitters. Our spectroscopy covers the range 1.7$-$5.1 micron in order to target strong optical emission lines: H$α$, [N II], [O III], and H$β$. Individual measurements as well as stacks reveal line ratios consistent with a metal poor nature (2$-$40% $Z_{\odot}$, depending on the calibration). The galaxies with the highest equivalent widths of Lyman-$α$, in excess of 90 Angstrom, have lower [N II]/H$α$ (1.9-$σ$) and [O III]/H$β$ (2.2-$σ$) ratios than those with lower equivalent widths, implying lower gas-phase metallicities at a combined significance of 2.4-$σ$. This implies a selection based on Lyman-$α$ equivalent width is an efficient technique for identifying younger, less chemically enriched systems. △ Less

Submitted 30 August, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 14 pages, 4 appendices; ApJ Accepted

MUSEQuBES: The relation between Ly$α$ emitters and C IV absorbers at $z \approx 3.3$

Authors: Eshita Banerjee, Sowgat Muzahid, Joop Schaye, Sean D. Johnson, Sebastiano Cantalupo

Abstract: We present a detailed study of the column density and covering fraction profiles of C IV absorption around 86 redshift $z \approx 3.3$ Ly$α$ emitters (LAEs) detected in 8 Multi-Unit Spectroscopic Explorer (MUSE) fields of $1'\times 1'$ centered on 8 bright background quasars as part of the MUSEQuBES survey. Using Voigt profile fitting of all the C IV absorbers detected along these 8 sightlines, we… ▽ More We present a detailed study of the column density and covering fraction profiles of C IV absorption around 86 redshift $z \approx 3.3$ Ly$α$ emitters (LAEs) detected in 8 Multi-Unit Spectroscopic Explorer (MUSE) fields of $1'\times 1'$ centered on 8 bright background quasars as part of the MUSEQuBES survey. Using Voigt profile fitting of all the C IV absorbers detected along these 8 sightlines, we generated a ``blind'' absorbers' catalog consisting of 489 C IV absorption components. We cross-matched this blind C IV catalog with the MUSE-detected LAE catalog and found a significant enhancement of C IV components within $\approx \pm$400 $\rm km\, s^{-1}$ of the systemic redshifts of the LAEs. Neither the C IV column density ($N$) nor the Doppler parameter ($b$) of individual C IV components shows any significant anti-correlation with impact parameter ($ρ$) of the LAEs in the 68 percentile range of $90\leq ρ\leq 230$ physical kpc (pkpc). We find a covering fraction of $\approx 60\%$ for a threshold $N$(C IV) of $10^{12.5}\, \rm \rm cm^{-2}$, which is roughly twice as high as in random regions. The C IV covering fraction remains constant at $\approx50\%$ for impact parameters in the range 150--250~pkpc ($\approx 3-6 R_{200}$). Using the covering fraction profile, we constrained the LAE--C IV absorber two-point correlation function, and obtained $r_0 = 3.4$ comoving Mpc (cMpc) and $γ= 1.2$ for a threshold $N$(C IV) of $10^{13.0}\, \rm cm^{-2}$. The C IV covering fraction is found to be enhanced for the LAEs that are part of a ``pair/group'' compared to the isolated ones. △ Less

Submitted 8 July, 2023; v1 submitted 10 April, 2023; originally announced April 2023.

Comments: 19 (+15 in the appendix) pages, 15 (+2 in the appendix) figures. Accepted for publication in MNRAS