Open Access
Issue
A&A
Volume 696, April 2025
Article Number A236
Number of page(s) 24
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202450785
Published online 28 April 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

The epoch known as the “cosmic noon” marks the peak of star formation and AGN activity in the Universe (Madau & Dickinson 2014). During this era, galaxy formation is largely driven by cold gas inflows from the surrounding structures (Dekel et al. 2009), and the star formation is regulated by such gas feeding as well as by gas outflows (feedback) due to supernova explosions and/or active galactic nucleus (AGN) activity (Genzel et al. 2014). On top of this, in high-density regions, there are more complexities introduced by the environmentally driven physical processes. According to current cosmological models, the channeling of cold gas streams towards the center of galaxies, the effect of dynamical friction, and the onset of environmental effects contribute to compressing and altering the gas distribution (Tacchella et al. 2016) of protoclusters members, possibly triggering central starburst (Gómez-Guijarro et al. 2019) and supporting their inside out mass growth (van Dokkum et al. 2015). Young forming high-z protoclusters are observed as strong overdensities of gas-rich, often dusty, star-forming and starbursting galaxies (Dannerbauer et al. 2014, 2017; Casey 2016; Casey et al. 2017; Calvi et al. 2023; Toshikawa et al. 2024) as simulations also predict (e.g., Chiang et al. 2017; Remus et al. 2023). This galaxy population is believed to be the predecessors of the red ellipticals that start dominating the cores of the clusters by z = 1 (e.g., Ivison et al. 2013; Smail et al. 2014). Thus, the rapid growth of galaxies within protoclusters is driven by the efficient transformation of the cold gas reservoir into stars until the sudden shutting down of star formation. However, the exact mechanisms terminating this phase of accelerated galaxy evolution remain unsettled, with both environmental effects and AGN feedback likely playing a role in the emergence of the first early-type galaxies in overdense environments (e.g., Krishnan et al. 2017; Polletta et al. 2021; Mei et al. 2023).

The last two decades have seen a growing number of works trying to unveil the interplay between the protocluster environment and the accelerated evolution of the galaxy populations therein in terms of star-formation and AGN activity, metal enrichment, and molecular gas content (see Overzier 2016; Alberts & Noble 2022 for a review). However, a unified picture depicting how the cluster assembly affects the early evolution of galaxies is still missing, with contradicting results between protoclusters at similar cosmic epochs probably due to the diversity of dynamical states that characterize these large-scale structures at the cosmic noon. In terms of star formation, several authors reported a higher activity within the most overdense regions at z > 2 thus suggesting the reversal of the so-called star formation density relation mediated by higher gas accretion and merging rates (e.g., Alberts et al. 2014; Shimakawa et al. 2018a; Lemaux et al. 2022; Monson et al. 2021; Shi et al. 2024a: Pérez-Martínez et al. 2024a; Laishram et al. 2024; Taamoli et al. 2024). However, others report no such differences (e.g., Toshikawa et al. 2014; Cucciati et al. 2014; Shi et al. 2021; Sattari et al. 2021; Pérez-Martínez et al. 2023, 2024a) in coeval protoclusters. A similar situation is found regarding the metal enrichment of galaxies with some authors claiming that the early development of an intracluster medium (ICM) in massive protoclusters would contribute to detaching infalling galaxies from the cosmic web, progressively hampering their gas accretion and suppressing their outflows, thus enriching their ISM (e.g., Kulas et al. 2013; Shimakawa et al. 2015; Pérez-Martínez et al. 2023). On the other hand, various levels of metallicity deficiency (e.g., Valentino et al. 2015; Chartab et al. 2021; Sattari et al. 2021; Pérez-Martínez et al. 2024a) or no significant environmental dependence at all (e.g., Tran et al. 2015; Kacprzak et al. 2015; Namiki et al. 2019) have been reported in other overdense structures at similar redshifts.

Furthermore, the slowdown of cold inflows in massive haloes is expected to take place over long periods (Zolotov et al. 2015), contradicting observational evidence showing that massive galaxies at z > 1.5 quench on much shorter timescales (Barro et al. 2016). Therefore, the reduction and eventual shutting down of star formation in massive galaxies is unlikely to occur solely because accretion is fading away. Instead, AGN feedback could provide an effective and complementary route to quenching massive protocluster galaxies after their starburst phase by expelling out a significant fraction of their gas reservoir via vigorous galactic winds and by directly heating the interstellar medium. Interestingly, several works have suggested that the AGN fraction is enhanced in protoclusters at 2 < z < 4 (see Lehmer et al. 2009, 2013; Krishnan et al. 2017; Vito et al. 2020, 2024; Polletta et al. 2021; Monson et al. 2021; Travascio et al. 2025). Furthermore, recent deep X-ray observations hinted that the massive and partially virialized Spiderweb protocluster at z = 2.16 may host up to 6 times higher AGN fraction than the coeval field (∼25% at M* > 1010.5 M, Tozzi et al. 2022a). By contrast, less evolved protoclusters in terms of their stellar-mass build-up (e.g. USS1558 at z = 2.53, Shimakawa et al. 2018a; Pérez-Martínez et al. 2024a; Daikuhara et al. 2024) do not show such enhancement of X-ray sources at similar depth (Macuga et al. 2019). This suggests that the transformative processes triggering AGN activity (e.g., angular momentum loss of the gas and/or merging) have yet to occur in such young protoclusters thus reflecting the diversity among the evolutionary stage of these large-scale structures. Furthermore, z ≲ 1 galaxies residing in virialized galaxy clusters display lower AGN activity than their coeval field counterparts (Haines et al. 2012; Ehlert et al. 2014; Koulouridis & Bartalucci 2019; Koulouridis et al. 2024), thus suggesting that the processes triggering the rapid supermassive black hole growth in cosmic noon protoclusters are short-lived and lose efficiency as these structures grow in mass and achieve virialization.

This puzzling situation has proven that the environmental effects present in protoclusters, if any, are mildly reflected on the properties of the ionized gas phase, and are easily washed out by the intrinsic scatter of the studied scaling relations due to relatively poor number statistics (e.g., Torrey et al. 2019; Huang et al. 2023). In order to get deeper insights we should move to investigate the physical quantities that give birth to the scaling relations responsible for the mass growth and metal enrichment of galaxies: the cold molecular gas reservoir. The molecular gas properties of protocluster galaxies at z > 2 are thus key to understanding the gas feeding and consumption processes fueling star formation and their changing efficiency as a function of both redshift, cluster assembly stage, and nuclear activity. Over the last years, there have been a growing number of works studying the molecular gas phase of star-forming protocluster galaxies using different CO transitions and dust continuum, although often over relatively small areas, sample sizes and focused on relatively bright sources (Tadaki et al. 2014, 2019; Dannerbauer et al. 2014, 2017; Emonts et al. 2016, 2018; Rudnick et al. 2017; Lee et al. 2017; Coogan et al. 2018; Wang et al. 2018; Zavala et al. 2019; Gómez-Guijarro et al. 2019; Champagne et al. 2021; Jin et al. 2021; Aoyama et al. 2022; Zhang et al. 2022; Polletta et al. 2022; Ikeda et al. 2022; Pensabene et al. 2024). Thus, it is important to expand these works to investigate the relatively mainstream star-forming population at this redshift range (i.e., main sequence galaxies) and cover a wider area within protoclusters to disentangle possible early environmental effects from other internal self-regulating ISM processes. However, very few fields have met these requirements up to now.

Over the last decade, a growing number of studies have identified star-forming (or star-bursting) protoclusters in the early universe with various techniques. The Spiderweb protocluster at z = 2.16 is a spectacular example of such high-z structures linked to a central radio galaxy MRC 1138–262 (so-called “Spiderweb”; Miley et al. 2006). This structure was originally discovered as an overdensity of Lyα emitters (LAEs, Pentericci et al. 2000; Kurk et al. 2000), and has gathered multiwavelength follow-up observations exceeding > 1200 hours over the last 25 years including X-ray and UgRIzJHKs photometry, as well as IR-to-millimetre mapping with Spitzer, APEX-LABOCA, Herschel, AzTEC, ALMA, and ATCA, thus confirming the overdensity through a large variety of galaxy populations (e.g., Kurk et al. 2004; Croft et al. 2005; Kodama et al. 2007; Zirm et al. 2008; Tanaka et al. 2010; Doherty et al. 2010; Kuiper et al. 2011; Hatch et al. 2011; Rigby et al. 2014; Shimakawa et al. 2014, 2015, 2018a; Shimakawa et al. 2024b; Dannerbauer et al. 2014, 2017; Zeballos et al. 2018; Tadaki et al. 2019; Emonts et al. 2016, 2018; Jin et al. 2021; De Breuck et al. 2022; Carilli et al. 2022; Tozzi et al. 2022a,b; Anderson et al. 2022; Pérez-Martínez et al. 2023; Chen et al. 2024; Naufal et al. 2023; Lepore et al. 2024; Daikuhara et al. 2024). In particular, a deep and wide-field (100 cMpc2) narrow-band search within this field, as part of the MAHALO-Subaru project (Kodama et al. 2013) yielded a huge (∼10 Mpc scale) overdensity traced by Hα emitters (HAEs, Koyama et al. 2013) whose membership has been confirmed by recent spectroscopic follow-ups (Shimakawa et al. 2015; Pérez-Martínez et al. 2023). In addition, recent results from a narrow-band JWST/NIRCam cycle 1 campaign have unveiled a new population of Paβ emitters that overlaps with the spatial distribution of HAEs and whose dust attenuation has been explored (Shimakawa et al. 2024a; Pérez-Martínez et al. 2024b). Furthermore, as part of the COALAS project, Jin et al. (2021) revealed an overdensity of CO(1−0) emitters at z = 2.1 − 2.2 hosting a large number of extended molecular gas reservoirs as shown by Chen et al. (2024). In addition, a large overdensity of bright dust-continuum emitting sources has been confirmed by Zhang et al. (2024) through ALMA observations. As a result, the dynamical mass of this protocluster exceeds Mcl ∼ 2 × 1014 M (Shimakawa et al. 2014). Recently, Tozzi et al. (2022a) revealed the presence of hot, diffuse baryons on a scale of ∼100 kpc from the Spiderweb galaxy using X-ray observations and Di Mascolo et al. (2023) reported signs of a nascent ICM via Sunyaev-Zeldovich effect implying this structure will soon become a bonafide galaxy cluster. Moreover, a handful of massive passive objects already populate the core of the protocluster (Naufal et al. 2024), adding evidence to the mature stage of this large-scale structure. Thus, the Spiderweb protocluster stands out as the best-known example of a confirmed galaxy cluster in formation at the cosmic noon and as an ideal test site to search for the first environmental effects from a multiwavelength perspective.

In this work, we search for environmental effects acting over the star formation activities and molecular gas reservoirs of Hα emitters belonging to the Spiderweb protocluster at z = 2.16. We structured this manuscript in the following way: Sect. 2 describes the observations and datasets used available within this field and part of our work. Sect. 3 outlines the methods used to analyze the physical properties of our targets and the environmental parameters measured. Sects. 4 and 5 respectively present our main results and their physical interpretation in the context of galaxy evolution in overdense environments at the cosmic noon. Finally, Sect. 6 outlines the major conclusions of this study. Throughout this article, we adopt a flat cosmology with ΩΛ = 0.7, Ωm = 0.3, and H0 = 70 km s−1 Mpc−1 and we assume a Chabrier (2003) initial mass function (IMF). Finally, all magnitudes quoted in this paper are in the AB system (Oke & Gunn 1983).

2. Observations and datasets

This section provides an overview of the datasets employed in our analyses, including their origins and key attributes. Our work relies on two main collections of data: the COALAS project1 whose focus is to identify CO(1−0) emitters across the Spiderweb protocluster structure and investigate their properties and dependencies with the environment; and the MAHALO project (Kodama et al. 2013), which aims at selecting narrow-band Hα emitters in high-z overdensities and characterize them through subsequent spectroscopic follow-ups. Additionally, we leverage the extensive repository of archival multiwavelength broad-band photometric data available for the Spiderweb protocluster field.

2.1. COALAS

The COALAS project is a large program (ID: C3181, PI: H. Dannerbauer) conducted with the Australian Telescope Compact Array (ATCA). Its primary objective is to explore into the influence of cosmic environments on the molecular gas reservoirs of galaxies in different environments during the cosmic noon, specifically through the ground-state transition of CO (i.e., J = 1). Furthermore, the CO(1−0) transition is the most robust tracer of the cold molecular gas and has been extensively used in the past for this purpose even at z > 1.5 (e.g., Dannerbauer et al. 2009; Carilli et al. 2010; Ivison et al. 2011; Emonts et al. 2013). The observations used in this work (i.e., approximately half of the COALAS project) took place from April 2017 to March 2020 targetting the central regions of the Spiderweb protocluster with 13 overlapping pointings, thus forming a mosaic of approximately 25.8 arcmin2 (Fig. 1). Throughout this program, approximately 475 hours of integration time were accumulated in this field with a central observed frequency of 36.5 GHz corresponding to the CO(1−0) emission line at z = 2.162, and bandwidth of 2 GHz spanning ±7000 km/s in velocity space across the line of sight of the protocluster in channels of ±40 km/s each and with a pixel size of 1.5″. Due to the variety of array configurations and integration times between pointings, the measured rms and beam size significantly vary across the mosaic. The typical rms level per beam for individual pointings ranges at 0.13 − 0.29 mJy per channel, while the maximum and minimum beam sizes found in the mosaic are 4.8″ × 3.5″ and 13.8″ × 13.1″, which corresponds to roughly 30 kpc and 100 kpc across respectively. Because of this low spatial resolution, the peak flux of a source is representative of the total flux of such an object. We refer to Table 1 in Jin et al. (2021) for further details on the array configurations for each individual pointing and data reduction details. In total, Jin et al. (2021) identified 46 CO(1−0) emitters at S/N > 4 forming an overdensity that spans a redshift range of z = 2.12 − 2.21, thus suggesting the presence of a large-scale filamentary structure surrounding the Spiderweb protocluster or a galaxy superprotocluster similar to the recently discovered Hyperion at z = 2.4 (Cucciati et al. 2018) or Elentári at z = 3.3 (Forrest et al. 2023). The probability of false CO(1−0) line detection for these sources as defined by Jin et al. (2021) is typically < 10%. In this work, we will use the COALAS mosaic to explore the evolution of molecular gas reservoirs in star-forming galaxies residing in the Spiderweb protocluster.

thumbnail Fig. 1.

Spiderweb protocluster field at z ≈ 2.16. Blue circles display the full sample of candidate HAEs from Koyama et al. (2013) and Shimakawa et al. (2018b). Red circles show those HAEs with measured spectroscopic redshift (Shimakawa et al. 2018b; Pérez-Martínez et al. 2023). Orange crosses and contours respectively show the CO(1−0) emitters reported by Jin et al. (2021) and the limits of the COALAS ATCA footprint. Empty squares depict the subsample of SMGs reported by Dannerbauer et al. (2014) in this field. Black stars show the LAEs from Pentericci et al. (2000), while X-ray sources from Tozzi et al. (2022a) are shown by black crosses. The Spiderweb galaxy is marked with a yellow star. We display a bar with 1 Mpc physical scale at z = 2.16 in the lower-left corner of the diagram for reference.

2.2. MAHALO-Subaru

The Spiderweb protocluster field also has extensive multiwavelength spectrophotometric coverage from the MAHALO-Subaru project, which aims at investigating the evolution of the star-forming activities of galaxies residing in overdensities up to z = 2.5 (Kodama et al. 2013). In particular, the available datasets include Subaru/Suprime-Cam B and z′ bands, Subaru/MOIRCS J and Ks-band and Subaru/MOIRCS narrow band observations with the filter NB2071 to identify HAEs (Koyama et al. 2013; Shimakawa et al. 2018b) across the Spiderweb protocluster structure. The MAHALO-Subaru project also conducted a series of spectroscopic follow-ups over the existing population of narrow-band selected HAEs yielding > 50 spectroscopic redshift through past programs with SUBARU/MOIRCS (Shimakawa et al. 2014, 2015) and VLT/KMOS (Pérez-Martínez et al. 2023). In this work, we will use these spectroscopically confirmed protocluster HAEs as our parent sample. In addition, we will crossmatch the position of these sources with the ATCA mosaic from the COALAS project (Jin et al. 2021) to obtain CO(1−0) molecular gas masses for some sources and upper limits for the rest. In total, the Spiderweb protocluster field has 54 spectroscopically confirmed narrow-band selected HAEs (Pérez-Martínez et al. 2023; Shimakawa et al. 2024b) and the ATCA mosaic from the COALAS project overlaps with the positions of 49 of them as shown in Fig. 1.

Given its distinct nature, we exclude the Spiderweb Galaxy complex from our analyses. This is often referred to as the ∼100 kpc region within the giant Ly-alpha nebula that surrounds the radio galaxy MRC 1138−262 and its satellite galaxies (“flies”, Kuiper et al. 2011). Through hierarchical merging, this region is likely destined to evolve into a single Brightest Cluster Galaxy (BCG) at low redshift (Pentericci et al. 1997; Hatch et al. 2008, 2009). In addition, we remove another five HAEs. Two of them (IDs 121 and 1185, Koyama et al. 2013) lie in the noisy region close to the edge of the ATCA mosaic. Consequently, they are unsuitable for accurate measurements (see Jin et al. 2021). Another two HAEs (IDs 469 and 1106) were spectroscopically confirmed only by their CO(1−0) emission. However, these redshifts would position their Hα emission at the edge of NB2071, where the transmission is less than 10%, rendering any attempts to measure Hα with the narrow-band filter unreliable. One additional source has a redshift mismatch of approximately Δz ≈ 0.04 between its Hα and CO(1−0) emissions, which complicates its redshift determination and leads to its exclusion from the analysis. Thus, the final sample of HAEs within the ATCA mosaic analyzed in this work consists of 43 sources. Figure 2 presents the redshift distributions for the spectroscopic HAEs (Pérez-Martínez et al. 2023; Shimakawa et al. 2024b), the CO(1−0) sources (Jin et al. 2021), and the specific subsample used in this study.

thumbnail Fig. 2.

Redshift distribution of sources. The complete sample of spectroscopic HAEs (Shimakawa et al. 2018b, 2024b; Pérez-Martínez et al. 2023) is shown in red. The sample of ATCA CO(1−0) emitters (Jin et al. 2021) is shown in yellow. The sample of HAEs within the ATCA footprint (see Fig. 1) that is used in this work is depicted by the blue contour histogram. The black solid line shows the shape and redshift coverage of the Subaru/MOIRCS narrow-band filter used to select the HAEs.

2.3. Multiwavelength broad-band photometry

In addition to datasets coming from the previous standalone projects, this field also counts with further NIR broadband deep imaging taken by VLT/HAWK-I in Y, H, Ks (PI: J. Kurk, program IDs 088.A-0754, 091.A-0106, 094.A-0104, see Dannerbauer et al. 2017; Shimakawa et al. 2018b). Furthermore, part of the field overlaps with broadband imaging at 3.6 and 4.5 μm made by Spitzer/IRAC (PI: D. Stern, campaign IDs 736 and 793, see Seymour et al. 2007) whose Post-BCD (PBCD) products can be retrieved from the Spitzer Heritage Archive (SHA). In addition, the Hubble Legacy Archive provides a small mosaic of reduced Hubble Space Telescope (HST) ACS/WFC data for filters F475W and F814W (PI: H. Ford, proposal ID 10327, see Miley et al. 2006). Finally, recent deep X-ray images with Chandra (PI: P. Tozzi, proposal ID: 20700463, see Tozzi et al. 2022a) surveyed both the COALAS and MAHALO-Subaru field in the Spiderweb protocluster thus providing a list of X-ray emitting sources whose positions have been crossmatched with the known existing HAEs to pinpoint the location AGN candidates within the protocluster structure. We refer to the main publications accompanying these datasets for further details on their instrument configurations and depths.

3. Methods

3.1. Molecular gas masses from CO(1–0)

In this section, we use the COALAS datacube mosaic in the Spiderweb protocluster to estimate the molecular gas masses for the spectroscopically confirmed protocluster sample of 43 HAEs. We gather the CO(1−0) fluxes at S/N ≳ 4 reported by Jin et al. (2021) for their sample of 14 dual Hα + CO(1 − 0) emitters. We remove the Spiderweb radio galaxy due to its distinct nature and three other sources from this list following Sect. 2.2. This results in a final sample of ten dual Hα + CO(1 − 0) emitters from Jin et al. (2021). We need to recover new and reliable CO(1−0) flux upper limits for the remaining 33 spectroscopic HAEs across the datacube footprint. Jin et al. (2021) based his CO(1−0) S/N measurements on the rms noise at various locations in the datacube. However, because the rms noise varies substantially between the pointings in the mosaic, and because we want to rule out that low-level instrumental effects across the ATCA band may affect our data (see Emonts et al. 2011), we will adopt a two-step process to derive our upper limits based on the known sky coordinates and redshifts of these 43 HAEs (Shimakawa et al. 2018b; Pérez-Martínez et al. 2023).

First, we take the sky coordinates of the spectroscopic HAEs (Koyama et al. 2013; Shimakawa et al. 2018a; Pérez-Martínez et al. 2023). The ATCA configuration during the COALAS project yielded a relatively wide but also variable beam size across the surveyed field. However, the typical distance between the Hα and CO(1−0) peak emission of dual emitters is less than twice the spaxel size (i.e., < 3″). Thus, we take grid of 3 × 3 spaxels around each Hα emitter’s position as the tentative locations to measure its flux upper limits (i.e., peak flux). However, to decide which spaxel contains the highest peak flux we need to define a spectral window to integrate the CO(1−0) flux upper limit, and then repeat the process across the grid in an iterative way.

Thus, we use the spectroscopic redshifts measured for each HAE as a reference point in the spectral dimension to search for CO(1−0) emission. Then, we define a spectral window of 1000 km/s width around the reference redshift of each source. This number has been chosen by taking into account that the median full width at zero intensity (FWZI) of the COALAS sample is approximately 600 ± 400 km/s (Jin et al. 2021) and thus, our spectral window encompass the flux contained up to 1σ beyond the median FWZI value. Clustering effects (i.e., signal overlapping) are only relevant in the vicinity of the Spiderweb complex given its bright and extended CO(1−0) emission. However, this part of the mosaic also has the best spatial resolution. Only two objects ID 1440 and 1501 would be affected by this situation and are considered as upper limits (see Table A.1), albeit they are separated from the systemic velocity of the Spiderweb galaxy by ≳2 × FWHMSW, CO(1 − 0) (Jin et al. 2021), implying that strong contamination is unlikely.

Finally, our approach integrates only the positive flux channels within the spectral window, thus guaranteeing that we are not underestimating the measured flux for our upper limits. Last, we repeat this process for all the spaxels within the 3 × 3 spatial grid centered on the peak position of the Hα emission and select the final upper limit for each source as the highest recovered flux between these 9 spaxels. In Fig. 3, we show two examples of this method applied to a source with CO(1−0) detection at S/N > 4 in Jin et al. (2021) and to a source below this threshold (i.e., not reported) in the same work. We achieve a median S/N ≈ 2.3 for our upper limits following the prescription given by Eq. (1) in Jin et al. (2021) and assuming FWZI = 1000 km/s, which represents the width of our spectral window. Upper limits at S/N < 1 are discarded and we take the rms value itself as the final flux estimate. We summarize our upper flux limits measurements among other properties of our sample in the appendix (Table A.1). Following our method, there are four cases (IDs 229, 999, 1139, 1385) out of the 33 spectroscopic HAEs in this field for which we recover CO(1−0) upper limits at S/N > 4 following our method. These objects were unaccounted for by Jin et al. (2021). However, as stated above, our method only integrates the positive flux within a spectral window, which can result in a flux overestimation. Consequently, we still consider these four CO(1−0) measurements as upper limits, which is reflected in Table A.1. These sources may have fallen near but below the S/N threshold imposed by Jin et al. (2021) in their blind emission line search, thus explaining why they were not selected by them.

thumbnail Fig. 3.

Two examples of our upper limit determination method. The upper panels display the HAE ID 790 for which Jin et al. (2021) reported CO(1−0) at S/N > 4, while the lower panel displays an undetected CO(1−0) source (HAE ID 1181) for which we compute upper limits. The first column displays the moment zero of the ATCA datacube after collapsing it at the redshift of our source ±500 km/s. The black cross shows the exact position of the HAE under scrutiny. The red empty square shows the position of the spaxel containing the peak flux after inspecting a 3 × 3 grid around the HAEs spaxel. The beam size is shown in the lower left corner. The second column depicts the spectra extracted from the peak flux spaxel. The vertical red dashed line shows the systemic velocity (i.e., redshift) of the galaxy from Hα while the solid red lines display the limits (±500 km/s) for flux integration (blue area). Velocities are relative to z = 2.1612 as in Jin et al. (2021).

To test the accuracy and reliability of our method, Fig. 4 compares the CO(1−0) fluxes of the sources as reported by Jin et al. (2021) S/N > 4 with the measurements obtained from those same sources following our method. There is a good agreement between both approaches within the error bars. The average signal-to-noise ratio of these ten sources following our method is S/N = 4.5 ± 1.5, while Jin et al. (2021) measurements yields S/N = 4.6 ± 1.0. This suggests that both our spatial spaxel grid and spectral window are reasonable choices to compute upper limits over known HAEs. Thus, we compute the molecular gas mass for each member of our sample based on the results published by Jin et al. (2021) on the ten dual Hα + CO(1 − 0) emitters and the new upper limits that we obtained for the 33 remaining HAEs in our sample. The CO(1−0) luminosity (LCO(1 − 0)) can be converted into a molecular gas mass estimate in a simple way assuming a conversion factor (αCO(1 − 0)) that traces the amount of molecular gas from optically thick virialized clouds (Dickman et al. 1986; Solomon et al. 1987):

M mol = α CO ( 1 0 ) × L CO ( 1 0 ) . $$ \begin{aligned} M_{\mathrm{mol}}=\alpha _{\rm CO(1{-}0)}\times L\prime _{\rm CO(1{-}0)}. \end{aligned} $$(1)

thumbnail Fig. 4.

Comparison of the LCO(1 − 0) obtained by applying our methodology and the source detection codes outlined in Jin et al. (2021) for HAEs with CO(1−0) detection at S/N > 4.

The αCO(1 − 0) factor is known to be approximately constant for star-forming galaxies close to the solar gas-phase metallicity value in the local Universe (αCO(1 − 0)= 4.36 ± 0.90 M/(K km s−1 pc2); Bolatto et al. 2013) albeit it depends on additional factors such as temperature, cloud density and metallicity (e.g., Narayanan et al. 2012). For example, lower metallicity values as those expected in the early universe may yield higher conversion factor values due to the contraction of the CO-emitting surface relative to the area where the gas is H2 for a fixed cloud size (Genzel et al. 2015; Tacconi et al. 2018).

Some sources within our sample have gas-phase metallicity estimates based on the [NII]/Hα ratio (Pérez-Martínez et al. 2023). However, the authors showed that these sources are representative of the metal-rich end of the mass-metallicity distribution and that the scatter of such relation for individual sources is ∼0.3 dex, which would introduce high uncertainties in the (αCO(1 − 0)) determination. Taking as a reference the stacking metallicity results of Pérez-Martínez et al. (2023) for HAEs in the Spiderweb protocluster we find that 12 + log(O/H) = 8.29, 8.55 and 8.64 for mass bins at 8.9 < log M*/M < 10.0, 10.0 < log M*/M < 10.8 and log M*/M > 10.8 respectively. If we take Eq. (2) in Tacconi et al. (2018) to compute the correction factors to αCO(1 − 0) based on these results we find that fcorr = 3.3, 1.15, and 1.0 respectively. This implies that galaxies in the intermediate and massive bin could increase their molecular gas mass by < 0.1 dex on average, but their effect in fgas would be negligible. The impact becomes stronger in the low-mass and low-metallicity end. Overall, the change of αCO(1 − 0) would increase the molecular gas mass of log M*/M < 10.0 sources by ∼0.5 dex, also increasing their gas fractions. However, fgas ≳ 0.9 for most of these (see Fig. 6) thus minimizing its impact on this property. Nevertheless, we note this correction diverges rapidly at low-metallicities and it must be used with caution in this regime. We refer to Tacconi et al. (2018) and references therein for further details. Therefore, we resort to using a single constant conversion factor value, αCO(1 − 0) = 4.36 M/(K km s−1 pc2), for our entire sample.

3.2. Star-formation activities

Our sample of protocluster HAEs for which we have obtained CO(1−0) flux measurements and upper limits has deep narrow-band and spectroscopic follow-ups (e.g., Shimakawa et al. 2015, 2018b; Pérez-Martínez et al. 2023) targeting the Hα emission line. To compute the SFR of our HAEs, we first retrieve Hα fluxes from VLT/KMOS datacubes as published in Pérez-Martínez et al. (2023), which contains 32 out of the 43 HAEs analyzed in this work. We note that the spectral resolution of KMOS allows us to separate the Hα and [NII] emission lines in our sources (see Pérez-Martínez et al. 2023) for further details). For the remaining 11 HAEs, we resort to the available narrow-band imaging to extract their Hα fluxes following Shimakawa et al. (2018b). We correct for filter transmission and [N II] contamination by using the mass-metallicity relation from Pérez-Martínez et al. (2023) and our SED-based stellar masses (see Sect. 3.3). Then, we use the measurements of Hα from these past works as a proxy for star formation within our sample and apply the Kennicutt (1998) calibration modified for a Chabrier IMF to compute their star-formation rates (SFR). This calibration has proven to be one of the most reliable both at local and high redshift (e.g. Moustakas et al. 2006; Wisnioski et al. 2019):

SFR ( H α ) = 4.65 × 10 42 × L ( H α ) $$ \begin{aligned} {\mathrm{SFR}(\mathrm{H}\alpha )}=4.65\times 10^{-42}\times L({\mathrm{H}\alpha }) \end{aligned} $$(2)

where L(Hα) is the luminosity of the Hα emission-line. We assume a Calzetti’s extinction law (Rv = 4.05, Calzetti et al. 2000) to account for the dust attenuation and follow the extinction correction calibration presented in Wuyts et al. (2013) to obtain our final SFR values. In this calibration, the total extinction over the Hα line (A(Hα)) is computed as a sum of both the continuum and the nebular contributions produced in active star-forming regions yielding A(Hα) = Acont + Aextra. We use the extinction Av values obtained from SED fitting to account for the diffuse dust attenuation in the galaxy’s continuum following Acont = 0.82Av, SED, while A extra = 0.9 A cont 0.15 A cont 2 $ {A_{\mathrm{extra}}=0.9A_{\mathrm{cont}}-0.15A_{\mathrm{cont}}^{2}} $ (Wuyts et al. 2013; Wisnioski et al. 2019). This calibration is in good agreement with the previous extinction estimates made by Calzetti et al. (2000) in the local universe. Our sample of 43 HAEs includes 9 X-ray emitters (Tozzi et al. 2022a), which are likely indicative of AGN activity. Two of them are identified as broad-line AGN (ID 647 and 911) by Pérez-Martínez et al. (2023) and their SFR is computed from their decomposed narrow-line component. Nevertheless, we are unable to quantify possible contamination from AGN activity or shocks into the narrow-line component of any of our sources and assume that it originates from star-formation activities in the following sections. The SFR of our protocluster sample can be found at the end of this work in Table A.1.

3.3. Stellar masses and dust attenuation

The field of the Spiderweb protocluster is covered with extensive and deep photometry comprising the rest-frame X-ray to NIR wavelength range for galaxies at z = 2.16 (see Sect. 2). We make use of SED fitting over the photometry retrieved from these datasets to derive stellar masses and dust extinctions for our sample of HAEs following the same procedure outlined in Pérez-Martínez et al. (2023) and using the same photometric bands. In the next paragraphs, we provide a brief summary of the main steps carried out during this process. First, we take the Subaru/MOIRCS narrow-band image as the reference in terms of seeing (FWHM = 0.7″) and pixel size (0.167″). We carried out PSF matching and pixel scale resampling over the remaining imaging except for the Subaru/Suprime-Cam B-band and the Spitzer/IRAC bands whose PSF is significantly larger (i.e., FWHM ≈ 1.1″ and ∼2″ respectively), than that of the narrow-band images.

Then, we perform dual-image photometry with SExtractor (Bertin & Arnouts 1996) using the NB2071 image for source detection while measuring the observed magnitudes (MAG_AUTO) on the rest of the images. In addition, we carry out independent single-band photometry for the Subaru/Suprime-Cam B-band and the Spitzer/IRAC 3.6 μm and 4.5 μm bands. As a result, we construct a 12-band multiwavelength catalog for our sources encompassing the rest frame 1400 − 14200 Å wavelength range for protocluster members at z ≈ 2.16. Each of our sources is covered by two to five bands with central wavelengths beyond 4000 Å in the rest frame (i.e., HAWKI/VLT H and Ks bands, Subaru/MOIRCS Ks band, and IRAC 3.6 and 4.5 μm). These bands are particularly important for estimating the stellar mass and dust attenuation of our sources. Spitzer/IRAC fluxes in the 3.6 and 4.5 μm bands were measured for approximately 60% of our sample. Despite the larger PSF size of Spitzer/IRAC compared to other bands in our dataset, blending issues do not significantly impact most of our sources.

We run the SED fitting code CIGALE (Code Investigating GALaxy Emission, Boquien et al. 2019) over this multiband catalog to derive the stellar mass and dust extinction for the sources within our sample. CIGALE is a widely used SED fitting routine that combines the modeling of composite stellar populations with nebular emission and dust attenuation while conserving the energy balance between the UV and the IR.

We create a grid of SED models based on CIGALE (Burgarella et al. 2005; Noll et al. 2009; Boquien et al. 2019) modules and parameters. Our grid is based on the stellar population synthesis models of Bruzual & Charlot (Bruzual & Charlot 2003) using the Chabrier IMF (Chabrier 2003) and subsolar metallicity (i.e. Z = 0.004). We assume an exponentially delayed star formation history with possible e-folding times between 0.1 and 8 Gyr and stellar population ages constraint to be younger than the age of the Universe at z = 2.16 (i.e., ∼3 Gyr). We allow for the presence of a recent (< 300 Myr old) but small (< 10% in mass) burst of star formation to account for possible recent interactions within the overdense protocluster environment. We account for the effects of nebular emission by setting the ionization parameter (U) range at −2.4 < log(U) <  − 2.8, which describes typical values for star-forming galaxies at z ∼ 2 (Cullen et al. 2016). Then, we apply Calzetti’s attenuation law (Calzetti et al. 2000) with extinction values ranging E(B − V)s = 0 − 2 mag in steps of 0.1 mag. Table 1 summarizes our CIGALE configuration. The parameters not included in this table are set to their default values.

Table 1.

Summary of CIGALE SED fitting modules and parameter configuration.

CIGALE follows a Bayesian-like approach by evaluating and weighing the models within a given grid depending on their proximity to the best fit, which has the heaviest weight. Based on the distribution of weights across the grid of models, the final physical properties and their uncertainties are estimated. This approach takes into account the uncertainties of the observations while also including the effect of intrinsic degeneracies between physical parameters (Boquien et al. 2019). Finally, the physical properties and their uncertainties are estimated as the likelihood-weighted means and standard deviations. The stellar mass and V-band attenuation (AV) values obtained for our sample carry uncertainties of the order of ±0.1 dex and ±0.2 mag respectively. We checked the effect of excluding the Spitzer/IRAC bands from our SED fitting analysis finding that the resulting stellar masses of our sources on average agree within their typical errors. Finally, in Fig. 5, we put our sample in context by comparisulting stellar massesng it with the so-called “Main Sequence” of star formation (Speagle et al. 2014) at z ≈ 2.16 similarly to what was shown in Pérez-Martínez et al. (2023) as we share 32 sources in common out of the 43 that comprise our HAE sample. Galaxies with CO(1−0) detection as well as those with upper limits scatter around the same locus and agree within uncertainties with the expectations from the coeval Main Sequence except for a few cases.

thumbnail Fig. 5.

Star-forming main sequence diagram. The black solid line and shaded grey region depict the Main Sequence of star formation at z = 2.16 from Speagle et al. (2014). Red small circles the spectroscopically confirmed HAEs of Shimakawa et al. (2018b) and Pérez-Martínez et al. (2023) within the ATCA datacube footprint. Blue circles show the location of the CO(1−0) emitters at S/N > 4 from Jin et al. (2021) which are part of our sample of HAEs. For context, we also add empty squares to show the overlap with the sample of SMGs in this field from Dannerbauer et al. (2014) and with the X-ray emitters reported by Tozzi et al. (2022a).

3.4. Environmental proxies

In this section, we aim to characterize the environment where our sources reside by inspecting both their immediate surroundings (i.e., local density) as well as their association with the larger scale structure of the Spiderweb protocluster (i.e., phase space). However, the quantification of the environment within overdense regions at high-z remains a challenge for most observational works due to three main reasons: projection effects, sample incompleteness, and structure definition. We can mitigate the first one thanks to the large number of spectroscopic redshifts (> 100) and narrow-band detections (HAEs and LAEs) available across the Spiderweb protocluster field (see Sect. 1). Nonetheless, limited projection effects could still be present due to the intrinsic proper motions of galaxies within protoclusters (Overzier 2016). The second is only partially alleviated thanks to the diversity of galaxy populations studied in this field over the last two decades, confirming the membership of a large number of both dust-obscured and unobscured star-forming galaxies and AGN throughout different methods (see Sect. 1), albeit the existence of a passive population and its distribution remains unexplored up to now. Finally, protoclusters are by definition large-scale structures studied amid their assembly process, and as such, measurements of their current size and total mass are highly uncertain (Muldrew et al. 2015). The discovery of a nascent ICM in the Spiderweb protocluster via Sunyaev-Zeldovich (SZ) effect confirms that this structure will evolve into a galaxy cluster by z = 0 (Di Mascolo et al. 2023) albeit the inner core might not yet be fully virialized, with current dynamical and SZ total mass estimates differing by a factor 2−4 (see Saro et al. 2009; Shimakawa et al. 2014; Di Mascolo et al. 2023). Nevertheless, we attempt to quantify the global environment by assuming the virialization of the protocluster core and inspecting the caustic profiles (η, Haines et al. 2012; Noble et al. 2013) in the phase-space diagram similarly to Wang et al. (2018) and Pérez-Martínez et al. (2023):

η = ( R proj / R 200 ) × ( | Δ v | / σ ) $$ \begin{aligned} \eta =(R_{\mathrm{proj}}/R_{200})\times (\left| \Delta v \right|/\sigma ) \end{aligned} $$(3)

where Rproj/R200 is the projected clustercentric distance relative to R200. This is the radius enclosing a volume with a mass density 200 times the critical density of the universe at the cluster’s redshift. In our analysis, we assume the Spiderweb Galaxy as the center of this protocluster. By comparison, the centroid of the X-ray diffuse emission (Tozzi et al. 2022a; Lepore et al. 2024) and the Sunyaev-Zeldovich signal (Di Mascolo et al. 2023) display offsets of 2.3 ± 0.6 and 6.2 ± 1.3 arcsec to the Spiderweb galaxy respectively. For context, the size of the Spiderweb complex (i.e., the radio galaxy and its surrounding interacting debris) encompasses at least 6−8 arcsec (Miley et al. 2006; Hatch et al. 2008; Kuiper et al. 2011) while our sources are distributed up to 3 arcmin away from the Spiderweb galaxy. Thus, we conclude that such small variations in the definition of the protocluster center have no significant impact on our environmental analysis. In addition, Δv is the line-of-sight velocity relative to the systemic velocity of the cluster defined as |Δv| = |(zzcl) c/(1+zcl)| with zcl = 2.156 being the systemic redshift of the Spiderweb galaxy (Kurk et al. 2000). Finally, σ is the dynamically estimated inner core’s velocity dispersion (σ = 683 km s−1) given in Shimakawa et al. (2014).

Alternatively, we can quantify the environment as a function of the local surface density of objects. To this end, we gather a sample of 145 unique protocluster members including spectroscopically confirmed LAEs (Pentericci et al. 2000), narrow-band selected and spectroscopically confirmed HAEs (Koyama et al. 2013; Shimakawa et al. 2024b), and spectroscopically confirmed CO(1−0) emitters (Jin et al. 2021). Narrow-band selected HAEs in this protocluster have been shown to have a membership success rate of > 90% (Pérez-Martínez et al. 2023), which is in line with the contamination rate from fore- and background sources reported by Sobral et al. (2013) in the field at z ∼ 2 using the same technique. Based on these results, we take narrow-band selected HAEs lacking spectroscopic confirmation as protocluster members. In addition, the survey area of our selected samples approximately overlaps with the full area of the COALAS footprint shown in Fig. 1, thus minimizing biases caused by uneven depths or coverage across this field. After these considerations, we compute the local surface density of objects as:

Σ N = N π R N 1 2 $$ \begin{aligned} \Sigma _N=\frac{N}{\pi R^2_{N-1}} \end{aligned} $$(4)

where N is the number of sources enclosed within a radius RN − 1, which is the distance to the Nth−1 nearest neighbor (i.e., not counting the galaxy of origin). The distance between sources is measured by their sky angular separation and transformed to the physical scale using a fixed cosmology described in Sect. 1. In this work, we compute local densities enclosing two, five, and ten neighboring galaxies (i.e. Σ2, Σ5, Σ10). This choice allows us to trace the properties of galaxies residing in local density peaks such as close companions or pairs (Σ2), compact groups (Σ5), and larger overdensities (Σ10).

4. Results

In this section, we present the molecular gas properties of our sample of 43 HAEs for which we computed CO(1−0) estimates within the Spiderweb protocluster. We recover ten sources from Jin et al. (2021) with CO(1−0) detections at S/N > 4 which represents ∼23 ± 12% of our spectroscopic sample and ∼20 ± 11% of all the detected narrow-band HAEs within the ATCA footprint (see Fig. 1 and references therein) assuming poissonian errors. This is consistent with the CO(1−0) detection ratio found by Tadaki et al. (2014) over HAEs within the USS1558−003 protocluster at z = 2.53 at a similar depth to ours (i.e., 0.13 − 0.29 mJy per beam in 40 km s−1 channels). However, they could only confirm 3 detections within 20 HAEs in the relatively small area surveyed by the VLA (∼1.5 arcmin2) in that protocluster. In our case, we also obtain CO(1−0) upper limits for the remaining 33 spectroscopically confirmed HAEs. This represents the largest sample of HAEs with available CO(1−0) measurements within a single cluster in formation at z > 2 up to date. To illustrate this further and to provide a broader context for the achievements of the COALAS project, we provide a collection of published CO(1−0) surveys carried out in clusters and protoclusters at the cosmic noon and several coeval field samples in the appendix (Table A.2). In the following sections, we will discuss the influence of the protocluster environment over the galaxies’ molecular gas reservoir and search for signs of accelerated galaxy evolution in overdense environments.

4.1. Gas fraction of protocluster HAEs

In Fig. 6, we inspect the gas fraction (i.e., fgas = Mmol/(M* + Mmol)) of our sources as a function of stellar mass and in comparison with other CO(1−0) coeval protocluster (Wang et al. 2018) and field samples (Kaasinen et al. 2019; Riechers et al. 2020; Frias Castillo et al. 2023) as well as field scaling relations (Popping et al. 2012; Sargent et al. 2014; Tacconi et al. 2018). We also present our median 1σFgas detection limit across the COALAS datacube as a function of stellar mass with a shaded grey area in Fig. 6. We note, however, that the rms across the ATCA mosaic is quite heterogeneous (see Sect. 2.1 and Jin et al. 2021). Consequently, our upper limits, which in some instances correspond to 1σ measurements (Sect. 3.1), naturally scatter around the boundary of this region, with variations depending on their position and associated local rms within the datacube. This constrains our ability to examine potential gas-depleted systems at the low stellar mass range (i.e., log M*/M ≲ 10.0), albeit it is sufficient to investigate the process of gas depletion at the massive end.

thumbnail Fig. 6.

Left: Total gas fraction versus stellar mass diagram. The solid green line and shaded region represent the field scaling relations for Main Sequence (MS) galaxies from Tacconi et al. (2018). We overplot our sample in blue color (both for measurements and upper limits) and compare it with other protocluster samples at the cosmic noon such as Wang et al. (2018) in orange, and coeval field samples (Kaasinen et al. 2019 in pink; Riechers et al. 2020 in grey; Frias Castillo et al. 2023 in green). In addition, we marked those galaxies with signs of AGN activity with red crosses (X-ray emission, Tozzi et al. 2022a) and red triangles (high [NII]/Hα, Pérez-Martínez et al. 2023) and those identified as SMGs by Dannerbauer et al. (2014) as black empty squares. The solid red line and red area depict the fit and uncertainty of the PKS1138 and CLJ1001 protocluster samples using a logistic function, as proposed by Popping et al. (2012). The grey shaded area displays the 1σ detection limit given by the median rms of the entire ATCA mosaic. Right: Total gas fraction versus stellar mass diagram after stacking. Our sample is divided into three bins: Low-mass galaxies (log M*/M < 10.5); massive galaxies (log M*/M > 10.5) excluding AGN candidates; and massive AGN candidates. In addition, median values and standard deviations for the comparison samples are displayed using the same symbol and color scheme as in the left-hand panel. The empty symbols depict the median values of the Kaasinen et al. (2019) and Frias Castillo et al. (2023) field samples after applying αCO ≈ 1, typical of high-z starbursts and SMGs (see Sects. 4.1 and 4.2). We apply a 0.02 dex cosmetic shift in M* to these two data points to improve the visibility of their error bars.

In our sample, HAEs above and below (upper limits) S/N = 4 in CO(1−0) display a similar behavior across the entire stellar mass range. In particular, at low and intermediate stellar masses (i.e., log M*/M < 10.5) we detect galaxies in CO(1−0) that display gas fractions close to unity (Fgas ≳ 0.8), something we do not see for galaxies with higher stellar mass (i.e., log M*/M ≥ 11.0). This suggests that these galaxies are still relatively young and will form the bulk of their stellar mass from their existing molecular gas reservoir. On the other hand, the most massive end displays a significant fraction of galaxies with rather low gas fraction values (fgas ≲ 0.4), indicating that if new fresh gas accretion is unable to replenish their reservoirs, they may soon start their quenching process due to the lack of cold gas to keep forming stars. Although our individual CO(1−0) data points cannot trace the lowest stellar mass regime, this trend nevertheless qualitatively agrees with the predictions of scaling relations in the field and galaxy evolution simulations.

The transition between these two regimes happens within a very narrow stellar mass range (log M*/M = 10.5 − 11.0) suggesting the presence of physical processes able to deplete the cold molecular gas reservoir of galaxies in a relatively efficient and quick way. This behavior is also reproduced by the z = 2.5 protocluster sample CLJ1001 of Wang et al. (2018) and other cosmic noon field samples (Kaasinen et al. 2019; Riechers et al. 2020; Frias Castillo et al. 2023) after renormalizing their Mgas values to the conversion factor used in this work (αCO = 4.36 M/(K km s−1 pc2)). This suggests that the physical origin of this trend might not depend on the environment. It is worth noting however that while the Riechers et al. (2020) sample overlaps well with our work and that of Wang et al. (2018) at 2 < z < 3, the field samples of Kaasinen et al. (2019) and Frias Castillo et al. (2023) appear to follow parallel sequences offset towards higher stellar masses. This offset is likely caused by sample selection effects as these two studies are predominately focused on massive (> 2 × 1010 M) IR-bright sources (e.g., SMGs selected by Herschel or ALMA), while our work traces stellar masses down to ≳109 M and our parent sample is dominated by typical main sequence HAEs at the cosmic noon (see Fig. 2, but also Shimakawa et al. 2018b; Pérez-Martínez et al. 2023). This also applies to the field sample of Riechers et al. (2020), which largely overlaps with the expected properties of Main Sequence galaxies as shown by Aravena et al. (2019) and Boogaard et al. (2019, 2020). We note that if instead we apply a conversion factor typical of high-z SMGs (i.e., αCO ≈ 1 M/(K km s−1 pc2), Hodge et al. 2012; Calistro Rivera et al. 2018; Birkin et al. 2021; Riechers et al. 2021; Frias Castillo et al. 2023; Liao et al. 2024) the median gas fraction of these two field samples would drop to fgas ≈ 0.25 − 0.30, consistent with the scatter of our sample at the massive end. We present some of the basic properties of these CO(1−0) samples in Table A.2, thus exposing additional differences in terms of median LCO(1 − 0) and surveyed area. Finally, we have excluded from our comparisons other literature protocluster samples whose gas masses are obtained using higher CO transitions (e.g., Tadaki et al. 2019) or dust continuum (e.g., Zavala et al. 2019; Aoyama et al. 2022) to prevent the propagation of additional uncertainties associated with their gas conversions. Nevertheless, these samples qualitatively agree with our trend albeit with a larger scatter.

Interestingly, the sharp change of properties at log M*/M = 10.5 − 11.0 seen in our sample was already reported by Popping et al. (2012) who attempted to model it by fitting a logistic function with the form:

M mol M mol + M = 1 exp ( B × ( log M A ) ) + 1 $$ \begin{aligned} \frac{{M_{\rm mol}}}{{M_{\rm mol}+M_*}}=\frac{1}{{\exp {(-B\times (\log M_*-A))}+1}} \end{aligned} $$(5)

where A and B are fitting parameters that depend on redshift in Popping et al. (2012). This function is shown in Fig. 6 by a dashed black line. Given that our sample displays a similar behavior, we use this same functional form to fit our protocluster sample and that of Wang et al. (2018) simultaneously obtaining that A = 10.93 ± 0.05 and B = −2.15 ± 0.33 (red solid line and shaded red area in Fig. 6). Removing Wang et al. (2018) from the fit would slightly smoothen the fit, thus increasing the values of A and B by 1σ and 0.5σ sigma respectively. Furthermore, based on previous works, we identified additional subpopulations such as submillimeter galaxies (SMGs, Dannerbauer et al. 2014) and X-ray emitters (Tozzi et al. 2022a) and emission line diagnostic AGN candidates (Pérez-Martínez et al. 2023) within our sample of 43 HAEs. Intriguingly, we find that while the SMGs are distributed across the whole range covered by our sample in the fgas − M* plane, the AGN candidates appear only at log M*/M ≥ 10.5 and dominate the low gas fraction regime in our sample with 5 out of 6 sources at fgas ≤ 0.4. Tozzi et al. (2022a) reported the AGN fraction is six times above the coeval field in this protocluster based on their X-ray Chandra imaging. This enhancement is driven however by relatively X-ray bright (L2 − 10 keV > 4 × 1043 erg s−1) and massive (log M*/M ≥ 10.5) galaxies in agreement with our findings. Nevertheless, Tozzi et al. (2022a) only find marginal detection corresponding to luminosities < 1041 erg s−1 when stacking both the spectroscopically confirmed and narrow-band selected protocluster members in the X-ray images below these thresholds, which indicates no significant AGN activity at lower masses and X-ray luminosities.

4.2. Gas fraction: Stacking analysis

To gain further insights into the typical gas fraction of protocluster HAEs across different mass regimes, we divide our sample into two bins and resort to stacking analysis regarding their CO(1−0) emission. Low-mass HAEs encompass 23 objects with log M*/M < 10.5 M while another 20 sources display log M*/M > 10.5 M. This latter group is divided into two bins, one representing massive HAEs without signs of AGN activity containing 8 galaxies, and another for the remaining 12 massive sources that are considered AGN candidates based on their X-ray emission or line diagnostics. The stacking analysis is performed using a similar approach than in Pérez-Martínez et al. (2024a). We stack the spectra of sources within each bin weighted by their individual noise values:

F stack = i n F i ( v ) σ i 2 / i n 1 σ i 2 $$ \begin{aligned} {F_{\rm stack}}=\sum _{i}^{n}\frac{F_{i}(v)}{\sigma _{i}^2}\Bigg /\sum _{i}^{n} \frac{1}{\sigma _{i}^2} \end{aligned} $$(6)

where Fi(v) is the peak flux density as a function of velocity for a given source, and σi is the average noise as described by Jin et al. (2021) for the location of each source in the COALAS datacube. We use the spectroscopic redshift of our sources to shift their expected CO(1−0) emission to v = 0 km/s before carrying out the stacking procedure. Fig. 7 displays the resulting stacked spectrum for each bin defined above. We measure the CO(1−0) line flux following the same procedure described in Sect. 3.1 for our individual sources, and compute its uncertainty by measuring the standard deviation of two regions 2500 Å wide to the left and right of the measured emission line velocity window (red solid lines in Fig. 7) and separated by ±500 km/s from it.

thumbnail Fig. 7.

Stacked spectra around the CO(1−0) emission line for the three bins defined in Sect. 4.2: Low-mass galaxies (log M*/M < 10.5), massive galaxies (log M*/M > 10.5) excluding AGN candidates, and massive AGN candidates. We achieve S/N ≈ 8, 10 and 13 for our stacked measurements respectively.

The results of this process are displayed in Table 2. This table gathers the median values and standard deviation of the main properties of our sample and the comparison samples introduced in the previous section. The right-hand panel of Fig. 6 shows the median Fgas and M* for all protocluster and field samples discussed assuming αCO(1 − 0) = 4.36 (filled symbols). For reference, and as discussed in Sect. 4.1, we compute the median values of the field samples of Kaasinen et al. (2019) and Frias Castillo et al. (2023) when assuming αCO(1 − 0) = 1 (empty pink circle and green triangle). This is a more appropriate conversion factor for high-z starbursts and SMGs (e.g., Hodge et al. 2012) and shifts these two samples from a parallel sequence to Eq. (5) as described in Sect. 4.1 to the massive and gas depleted end of our sequence (red solid line and shaded area in Fig. 6). In such case and given the high SFR of those samples, they may deplete their gas reservoir and become passive within a very short period. The three bins of our sample in PKS1138 define a sequence where galaxies at low and intermediate stellar mass display high gas fractions (Fgas ≈ 0.86). However, the massive end display diminished gas fraction values, with Fgas ≈ 0.55 for massive HAEs without signs of AGN activity and Fgas ≈ 0.38 for massive HAEs labeled as AGN candidates. This supports the scenario where both stellar-mass growth and AGN activity contribute to the gradual exhaustion of the gas reservoir in our sample, particularly for objects with a stellar mass beyond log M*/M > 10.5 M.

Table 2.

Median physical properties of the samples used in Sect. 4.

4.3. Depletion times

Similarly, we inspect the depletion times (τdep = Mmol/SFR) of our sample in Fig. 8 using dashed gray lines over the Mmol − SFR plane. Most of our individual sources lie within the lines of τdep = 1 − 3 Gyr indicating that, in the absence of further inflows and outflows of cold gas, the HAEs populating the Spiderweb protocluster at z = 2.16 will deplete their gas reservoirs and thus shut down their star formation by z = 1.0 − 1.6 concurring with the surge and dominance of the red sequence in the cores of the most massive galaxy clusters known at this cosmic epoch (e.g., XMMU J2235−2557 at z ≈ 1.4, Rosati et al. 2009; but also see Grützbauch et al. 2012; Nantais et al. 2016; Beifiori et al. 2017). A few objects within our sample display even lower depletion times τdep ≤ 0.5 Gyr with SFRs exceeding 100 M/yr. These objects are possibly experiencing their last episode of star formation while rapidly consuming their gas reservoirs. Furthermore, three out of the four sources within this category are labeled as AGN candidates, hinting at the co-evolution of star formation and supermassive black hole growth.

thumbnail Fig. 8.

Molecular gas mass versus star formation rate diagrams. Left: Individual sources from the samples displayed in Fig. 6. Right: Stacked (PKS1138) and median results of each sample displayed in the left-hand panel. Colors and symbols are the same as in Fig. 6. Dashed grey lines display depletion times between 0.1 and 10 Gyr.

Nevertheless, the current ICM conditions within the Spiderweb protocluster (Tozzi et al. 2022b; Di Mascolo et al. 2023) likely still allow the channeling of cold gas streams to some of its members albeit with diminishing efficiency over time as the protocluster virializes, thus extending the lifespan of their star formation activities beyond the reported depletion times. At the same time, additional processes such as mergers (Mei et al. 2023), gas stripping (Boselli et al. 2022), or AGN feedback (Heckman & Best 2014; Davé et al. 2020) may contribute to either extending or shortening these depletion times. We will discuss the possible interpretations of these results in Sect. 5 and their implications in the context of galaxy evolution in overdense environments at the cosmic noon.

The stacked and median measurements for the Spiderweb and the comparison samples are shown in the right-hand panel of Fig. 8. We observe a good agreement between Wang et al. (2018), Riechers et al. (2020), and our Spiderweb sample. On the other hand, the samples of Kaasinen et al. (2019) and Frias Castillo et al. (2023) are clearly shifted towards higher SFR and lower depletion times when assuming αCO(1 − 0) ≈ 1, consistent with the expectations for starburst galaxies.

4.4. Molecular gas properties and environmental effects

In Fig. 9 we examine the relation of the molecular-to-stellar mass ratio (i.e., μgas = Mmol/M*) with the star formation activities within protocluster galaxies (i.e., SFR and SFE) and with the environment (η) as defined by Eq. (3) (see also Haines et al. 2012; Noble et al. 2013). First, we show the phase-space distribution of our targets (left panel in Fig. 9) assuming the virialization of its core (Shimakawa et al. 2014). This diagram includes both the solid S/N > 4 detection of CO(1−0) over HAEs (N = 10) and the CO(1−0) upper limits following the method of Sect. 3.1. Both samples have been color-coded to reflect their offsets with respect to the μgas scaling relations of Tacconi et al. (2018). In addition, we also divide the phase-space into three areas to separate galaxies in the outskirts (η ≥ 2), infalling (0.4 ≤ η ≤ 2) or central regions (η ≤ 0.4). Most galaxies in the outskirts display μgas values more than 0.2 dex above the scaling relation of Tacconi et al. (2018), while a lower number of them can be found in the infalling and core regions, which hints at the presence of some weak environmental segregation. In the right-hand panel of Fig. 9, we find a similar picture over the log μgas − η plane with galaxies with higher absolute μgas values populating the outskirts of the protocluster (i.e., high η) and a declining trend in μgas towards the cluster core (i.e., low η). Similar trends were also reported by Wang et al. (2018) in CLJ 1001+0220 at z = 2.506 as a function of η, albeit most of their sources are concentrated in the protocluster inner region (R < 150 kpc) and their total surveyed area is relatively small (5.3 arcmin2) while this work examines a large mosaic of 25.8 arcmin2 thus probing a larger range of environments. On the other hand, we also examine the distribution of PKS1138 massive and low-mass galaxies, as well as AGN candidates in this diagram by using the three bins (pentagon and hexagons) defined in Table 2. Our results show that while the massive end and AGN candidates show a lower median μgas (similar to Fig. 6), they tend to lie in similar environments in terms of η although with relatively high dispersion.

thumbnail Fig. 9.

Left: Phase space diagram. Objects are color-coded by their molecular to stellar mass ratio μgas offsets with respect to the scaling relation of Tacconi et al. (2018). Solid grey lines divide the phase space into three regimes with galaxies being in the outskirts (η ≥ 2), infalling (0.4 ≤ η ≤ 2) or central regions (η ≤ 0.4) similarly to Pérez-Martínez et al. (2023). CO(1−0) measurements at S/N > 4 are depicted as solid circles while inverted triangles show upper limits. Right: Molecular to stellar mass ratio (μgas) offsets from the scaling relation as a function of the environmental parameter η. In addition, the pentagon and two hexagons depict the three stellar mass and AGN activity bins presented in Table 2. The remaining symbols and colors remain the same as in the left-hand diagram.

In addition, Fig. 10 examines the relation between the star formation activities of our HAEs and their molecular gas reservoirs. The left panel of Fig. 10 shows a flat trend between the offsets from the main sequence of star formation and η. This means that regardless of the distance to the center of the cluster, most of the surveyed galaxies remain part of the Main Sequence even when reaching the cluster core, thus suggesting that the environment is not actively promoting or suppressing star formation over those objects. However, the color code shows again that galaxies in the outskirts of the protocluster have larger μgas than their core counterparts, implying that the cold gas is gradually removed, heated up, or consumed as the galaxies approach the central regions. This can be seen more clearly in the right-hand panel of Fig. 10 where we display the offsets of the star formation efficiency (i.e., SFE = SFR/Mmol = 1/τdep) with respect to the scaling relation of Tacconi et al. (2018) as a function of η. Most objects display higher log(SFE/SFET18) values at lower η while this situation is inverted as we go further away from the protocluster core. However, our sample of HAEs is comprised of Main Sequence galaxies (Fig. 5 and the left panel in Fig. 10) regardless of their position across the structure protocluster. This suggests that the observed increase in SFE towards the protocluster center is predominantly driven by the depletion of the gas reservoir in the central regions of the protocluster, in agreement with our previous results and those of Wang et al. (2018). Nevertheless, it is common to find more massive objects in overdensities compared to the coeval field as well as towards the center of a given overdensity compared to its outskirts. This stellar mass segregation (e.g., Contini et al. 2012; van den Bosch et al. 2016) is usually linked with the earlier formation times of these objects compared to those in less dense environments. Thus, stellar mass build-up may also play a role in the observed depletion of molecular gas trends towards the protocluster center on top of which environmental effects may be acting as discussed in Sect. 5.

thumbnail Fig. 10.

Left: Star formation rate offsets from the Main Sequence of Tacconi et al. (2018) as a function of the environmental parameter η. Right: Star formation efficiency offsets with respect to the scaling relation of Tacconi et al. (2018) as a function of η. The symbols and colors remain the same as in Fig. 9.

Finally, Fig. 11 explore possible correlations between the gas fraction and the projected local surface density computed using the minimum area to enclose two, five, and ten galaxies (see Sect. 3.4). This approach allows us to identify projected local density peaks (e.g., close companions or compact groups) across the structure of the protocluster. Our results show no clear correlation between Fgas and any of these three proxies, with galaxies displaying gas fractions above and below Fgas = 0.5 spreading across most of the local density range probed by these proxies. Furthermore, we split our sample into three bins depending on stellar mass and AGN activity as in Fig. 6 but find no signs of strong environmental segregation in terms of these physical properties. However, we acknowledge that the statistical limitations of our sample, the large number of upper limits, and the intrinsic bias of this proxy due to projection effects may contribute to erasing potential underlying correlations between the local density and other physical properties.

thumbnail Fig. 11.

Gas fraction as a function of local density defined by the area enclosing two, five, and ten neighboring galaxies (Σ2, Σ5, Σ10. The symbols and colors remain the same as in Fig. 9.

5. Discussion

This work has investigated the evolution of the molecular gas reservoir of a sample of main-sequence HAEs within the Spiderweb protocluster at z = 2.16. The sharp evolution of the total molecular gas fraction (fgas) at log M*/M = 10.5 − 11.0 shown in Fig. 6 yields an evolutionary scenario composed of two main stages connected by a rapid transition at a given characteristic mass. In the low mass regime (log M*/M ≤ 10.5), protocluster galaxies display molecular gas fractions close to unity despite having star formation activities already compatible with those of the main sequence (Fig. 5). This implies that, while the bulk of their stellar mass remains to be formed in the future, the past and present stellar mass growth (i.e., SFR) has not depleted their gas reservoirs, suggesting that these objects keep replenishing their gas reservoirs through cold gas inflows from the cosmic web (Tadaki et al. 2019). On the other hand, the high mass end of our HAEs displays fgas ≤ 0.4 with some objects going down to fgas ≤ 0.2 at log M*/M = 11.0 − 11.5. These objects represent the final stage of the in-situ stellar mass growth in galaxies, showing still active star formation albeit with clear signs of gas reservoir depletion, hinting at the formation of passive objects by the time this process is completed (e.g., Falkendal et al. 2019; Williams et al. 2021; Whitaker et al. 2021; Zanella et al. 2023; Blánquez-Sesé et al. 2023; D’Eugenio et al. 2023). The gap between these two phases is bridged by a population of rapidly transitioning galaxies at (log M*/M = 10.5 − 11.0) displaying 0.3 < fgas < 0.9 and with a high fraction of AGN (12/20) at log M*/M > 10.5. Thus, the physical mechanism behind this transition must act on a relatively short timescale and predominantly act at and beyond the aforementioned stellar mass threshold, for which AGN feedback is a good candidate.

Furthermore, Fig. 8 shows that most HAEs in our sample have depletion times ranging from 1 to 3 Gyr, with a few of them showing even shorter time scales. This means that, in the absence of inflows and/or outflows, the star-forming population that dominates the Spiderweb protocluster would become passive by 1 < z < 1.6. Observational evidence suggests that galaxy protoclusters start the surge of their red sequence by z ≈ 2. In fact, past works in the Spiderweb protocluster already confirmed the presence of a handful of passively evolving galaxies near its core (Kodama et al. 2007; Zirm et al. 2008; Tanaka et al. 2010; Doherty et al. 2010; Naufal et al. 2024) in line with other such studies in protoclusters (e.g., Strazzullo et al. 2016; Noordeh et al. 2021). However, most observational evidence suggests that the red sequence experiences significant growth in massive clusters at 1.5 < z ≲ 2 (Grützbauch et al. 2012; Fassbender et al. 2014; Nantais et al. 2016; Beifiori et al. 2017) until it becomes the dominant population within the cluster cores by z = 1 in agreement with our depletion times estimates.

On the other hand, we have examined possible environmental impacts over our sample in Figs. 9 and 10 through the molecular-to-stellar mass ratio (μgas), the star formation rate (SFR), and the star formation efficiency (SFE) as a function of the phase-space environmental parameter (η, Eq. 3). Our findings suggest that the Spiderweb protoclusters host a large population of molecular gas-rich galaxies with lower star formation efficiency in its outskirts (η > 2) while this trend is reversed when approaching the protocluster core (i.e., lower μgas and higher SFE at η < 0.4). However, our sample displays star formation activities in line with the main sequence regardless of their location across the cluster structure. This was previously reported by Pérez-Martínez et al. (2023) and suggests that if environmental effects are at play, they do not have a significant impact on the star formation activities of the protocluster members for now and thus, any change in the molecular-to-stellar mass ratio μgas is primarily driven by changes in the cold gas reservoir while star formation proceeds unaltered. We propose three different mechanisms to explain this process: changes in gas accretion, gas removal through environmental effects, and AGN activity.

5.1. Changes between accretion regimes

Environmental effects and gas accretion from the cosmic web (Dekel & Birnboim 2006; Daddi et al. 2022a) could respectively hasten or delay the depletion of the gas reservoir. Current observational evidence can not rule out the presence of cold gas accretion within the Spiderweb protocluster. However, the recent discovery of a nascent ICM within its inner regions using Sunyaev-Zeldovich effect observations (Di Mascolo et al. 2023) suggests that if such accretion is still in place it may not proceed homogeneously across the protocluster structure. Thus, we would expect some galaxies to be fed through the so-called “cold streams in hot media” (Dekel & Birnboim 2006; Dekel et al. 2009) within the protocluster core. This process would progressively lose efficiency as the Spiderweb grows in mass and achieves virialization becoming a bona fide galaxy cluster.

Nonetheless, should cold accretion still be in place, we can roughly estimate its relevance in sustaining the star formation activities of the Spiderweb protocluster member galaxies through a series of simple assumptions. Several authors have provided practical approximations to estimate the expected baryonic accretion rate (BAR) within the ΛCDM cosmology for a given redshift and halo mass (e.g., Genel et al. 2008; Goerdt et al. 2010; Dekel et al. 2013) based on hydrodynamical simulations. Recently, Daddi et al. (2022a) investigated the transition between cold and hot accretion in a series of groups at the cosmic noon following Eq. (5) in Goerdt et al. (2010), which we reproduce here:

BAR 137 × ( M DM 10 12 M ) 1.15 ( 1 + z 4 ) 2.25 M yr 1 . $$ \begin{aligned} {\mathrm{BAR}\approx 137\times \left(\frac{M_{\rm DM}}{10^{12}\,\mathrm{M}_\odot }\right)^{1.15}\left(\frac{1+z}{4}\right)^{2.25}\,\mathrm{M}_{\odot }\,\mathrm{yr}^{-1}}. \end{aligned} $$(7)

Following the theory outlined by Dekel & Birnboim (2006), dark matter haloes below a given mass threshold (MDMMshock = 6 × 1011 M) accrete all their baryons throughout cold streams. Thus Eq. (7) approximately describes their total but also cold baryonic accretion rate. Above this limit, the cooling times are longer than dynamical times, and shocks can efficiently heat part of the gas. However, numerical simulations predict that cold accretion continues to penetrate at high-z (z > zcrit) in the form of cold streams, surviving the shocks at the virial radii of massive halos up to MDM = Mstream with log Mstream ≈ log Mshock + 1.11 × (z − 1.4) and zcrit = 1.4 (Dekel & Birnboim 2006). However, recent observational works claim a flatter slope and lower zcrit for this boundary (e.g., Daddi et al. 2022b). Dark matter haloes beyond this mass threshold (MDM ≳ Mstream) experience a smooth transition with decreasing cold accretion efficiency as their mass grows. Daddi et al. (2022a) modeled such transition as BARcold ≈ BAR × (Mstream/MDM)α with α being the modulation factor of such efficiency. Fig. 12 displays the Spiderweb protocluster accretion regime for the halo mass estimates using the Sunyaev-Zeldovich effect (M500 = 3.5 × 1013 M yr−1, Di Mascolo et al. 2023) and dynamical estimates (M200 = 1.7 × 1014 M yr−1, Shimakawa et al. 2014). By comparison, we have also included the sample of SZ-detected clusters at z < 1.5 (Bleem et al. 2015), a sample of protoclusters a 1.5 < z < 6.5 (Overzier 2016), and individual Spiderweb protocluster HAEs studied in this work after estimating their dark matter halo masses using the stellar-to-halo mass relation (SHMR) for satellite galaxies from Behroozi et al. (2019). At z = 2.16, this relation is essentially flat for M > 1010.8 M . We assign a maximum dark matter halo mass of MDM ≈ 1012.9 M to the most massive HAEs in our sample, thus preventing unrealistically large variations of MDM with small changes in stellar mass.

thumbnail Fig. 12.

Accretion regime diagram. Colored stars display the dark matter halo mass of the Spiderweb protocluster derived using the SZ effect and dynamical estimates by Di Mascolo et al. (2023) and Shimakawa et al. (2014) respectively. Grey squares depict the position of our sample of HAES after applying the stellar-to-halo mass relation in Behroozi et al. (2019) for satellite galaxies (see Sect. 5.1). For comparison, we use green big circles to show a sample of protoclusters from Overzier (2016) while the small red circles represent a sample of SZ-detected z < 1.5 clusters from Bleem et al. (2015).

Given these conditions and considering individual satellite dark matter haloes, a typical M = 1010 M HAE will lie within a halo of MDM ≈ 1012 M. This value is well below Mstream at z = 2.16 and thus, its cold accretion is equal to a BAR of roughly 80 M yr−1. This level of cold gas accretion is comparatively higher than the expected gas consumption via star formation for a main sequence galaxy (Speagle et al. 2014) of this stellar mass, with SFR ≈ 30 M yr−1. Thus, the gas reservoir of galaxies within this stellar mass regime would be replenished as they grow in mass through star formation, keeping their gas fraction close to unity even if the efficiency of converting the cold gas into stars is only 50%, as described by Dekel et al. (2009) models (see also Daddi et al. 2022a). On the other hand, a typical M = 1011 M galaxy will inhabit a dark matter halo of MDM ≈ 1013 M just above Mstream at this redshift and thus lying within the so-called hot accretion regime (Dekel & Birnboim 2006). Here, the cold BAR depends on the modulation factor α with recent studies suggesting that α ≈ 1.0 ± 0.2 for haloes within approximately ±1.5 dex of Mstream (Daddi et al. 2022a,b), which would be the case of our most massive HAEs. Thus, their cold BAR is roughly 475 77 + 90 M yr 1 $ {475^{+90}_{-77}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}} $ which would be enough to feed a M = 1011 M main sequence star-forming galaxy displaying SFR ≈ 180 M yr−1 when assuming the same cold gas-to-stellar mass conversion factor (i.e., ∼50%). This implies that in the absence of additional mechanisms such as gas removal or AGN feedback (see Sects. 5.2 and 5.3), the gas that is consumed by star formation in these objects would be replenished by cold gas accretion at this redshift. Thus, gas deprivation (i.e., so-called starvation) due to dark matter halo growth of individual objects could not explain alone the low gas fractions reported in Sect. 4 for the massive end of our sample of HAEs.

Nonetheless, the virial mass of the main dark matter halo associated with the Spiderweb protocluster is of the order of ∼ 1014 M with some of the satellite HAEs lying well within the virial radius (R200 = 0.53 Mpc, Shimakawa et al. 2014) and thus the accretion modulation factor, α, could be larger than unity thus diminishing the efficiency of cold accretion. To test this possibility, we inspect the gas fraction of the 12 HAEs (excluding the Spiderweb galaxy) lying within R200 in our sample ignoring projection effects. We find that seven of them are relatively low-mass galaxies and display fgas ≳ 0.7 despite their proximity to the Spiderweb galaxy, which we take as the approximate center for the halo of this protocluster. The other five objects, on the other hand, have M > 1010.5 M and 0.1 < fgas < 0.6. However, we find another five objects in our sample with the same stellar mass and gas fraction ranges outside R200. These findings suggest that the changes between accretion regimes that satellite galaxies experience as they infall towards the core of the Spiderweb protocluster are unlikely to be the main mechanism behind the gas depletion we observe at the massive end of Fig. 6.

5.2. Gas removal through environmental effects

Gas removal is a common effect discussed in the context of environmental effects and can obey different causes. Its most studied version in clusters at z < 1 is ram pressure stripping (RPS, Gunn & Gott 1972). Ram pressure stripping originates due to the dynamical friction between the hot dense gas of the ICM and the outer layers of the galaxies’ ISM (Boselli et al. 2022). The effectiveness of RPS is proportional to the square of the satellite’s infalling velocity and to the density of the medium it goes through and thus, while RPS is very effective in quenching low-mass satellite galaxies in low-z clusters (e.g., Laganá et al. 2013; Roberts et al. 2019; Mao et al. 2022), it is unlikely it has a significant impact in protoclusters due to the low-velocity proper motions and the relatively low density and likely heterogeneous distribution of a nascent ICM during the protocluster assembling phase. Nevertheless, recent works have measured an electron density of the order of ne ≈ 10−1 − 10−2 cm−3 and temperatures of kT ≈ 0.6 − 2 keV for a thermal component within the inner 10 − 100 kpc of the Spiderweb protocluster (Tozzi et al. 2022b; Lepore et al. 2024) implying the presence of a cool-core for the first time at z > 2. Despite these values being similar to those found in low-z galaxy clusters (e.g., Böhringer et al. 2004, 2007; Sanderson et al. 2006, 2009), the small extent of the X-ray emission, the steepness of its electron density profile, and the asymmetries found in the X-ray and SZ emission suggest that we are witnessing the first stages of a nascent ICM, still far from virialization, as suggested by Di Mascolo et al. (2023) and Lepore et al. (2024). Thus, we would expect RPS to be a suitable mechanism to remove significant amounts of gas only within close proximity of the Spiderweb galaxy (e.g., Kuiper et al. 2011) and in low mass galaxies (i.e., shallower gravitational potentials).

We explore this possibility by inspecting the four HAEs within our sample at a distance < 200 kpc from the Spiderweb galaxy. Two of them (IDs 1420 and 1498) are relatively low-mass galaxies (log M* ≤ 10.1) displaying high gas fractions (fgas = 0.7 − 0.9) arguing against the presence of RPS. In contrast, the other two objects (IDs 880 and 1501) are massive galaxies (log M* ≥ 10.9) displaying lower gas fractions (fgas < 0.3 − 0.5) and with confirmed AGN activity via X-ray emission, pointing towards different mechanism at play to explain the depletion of their gas reservoir. Furthermore, local universe studies in massive galaxy clusters have demonstrated that while loosely bound gas lying in the outskirts of low-mass satellite galaxies could be easily removed in large quantities via RPS (Fumagalli et al. 2014; Boselli et al. 2016, 2022) only extreme cases would have a similar effect over the denser and typically more centrally concentrated cold molecular gas (e.g., Stevens et al. 2021; Roberts et al. 2023). In that sense, the detection of extended (≥40 kpc) molecular gas reservoirs in satellite galaxies across the protocluster structure (Dannerbauer et al. 2017; Chen et al. 2024) also argues against RPS being an efficient gas removal mechanism at this stage. We have checked the abundance of extended molecular gas reservoirs as described by Chen et al. (2024) within our sample of HAEs finding that only four HAEs out of the ten CO(1−0) emitters at S/N > 4 (i.e., ∼40%) within our sample fulfill this requisite either as robust (3) or tentative (1) extended molecular gas reservoirs. One of these robust detections (ID 902 in Table A.1) hosts a remarkable extended molecular gas reservoir previously studied by Dannerbauer et al. (2017) and Chen et al. (2024). However, this ratio is similar but slightly inferior to the ∼46% of extended molecular gas reservoir candidates (robust + tentative = 21) found by Chen et al. (2024) out of the full COALAS sample of 46 CO(1−0) emitters. This suggests that Hα emitters do not preferentially host large molecular gas reservoirs compared to the COALAS sample.

Alongside RPS, simulations have demonstrated that significant gas outflows can also be mediated by stellar feedback (e.g., Feldmann & Mayer 2015; Bahé & McCarthy 2015; Tollet et al. 2019). The efficiency of this mechanism on the removal of cold gas and thus star formation quenching depends on the mass of the host galaxy (i.e., the depth of the gravitational potential) and the magnitude of the ongoing star-formation episode (Larson 1974). Therefore, low-mass galaxies would be more likely to lose part of their gas reservoirs than their most massive counterparts (e.g., El-Badry et al. 2016; Christensen et al. 2016; Keller et al. 2016; Pontzen et al. 2017; Bassini et al. 2023). The stronger gravitational pull of massive galaxies effectively confines the outflowing gas, causing it to fall back into the galaxy and form galactic fountains (Shapiro & Field 1976). During the early Universe, however, starburst episodes were more frequent thus resulting in higher levels of supernova feedback which in extreme cases can eject large amounts of metal-enriched gas beyond the virial radius (Kimmig et al. 2025), thus temporarily exhausting their star formation activities (Remus & Kimmig 2023). However, these same works show that stellar feedback alone is not enough to completely shut down star formation or prevent the inflow of new gas after a short period unless it couples with AGN activity and other environmentally driven processes (see also Sect. 5.3). The results presented in Sect. 4 do not allow us to discard that stellar feedback may be playing a role in individual cases, albeit the lack of a starbursting population (Fig. 5) and the high gas fractions reported in the low-mass end of our sample (Fig. 6) argue against this mechanism being dominant on the evolution of the observed populations at this stage.

Finally, tidal interactions between infalling satellites and the global cluster potential, or frequent high-speed galaxy encounters (i.e., harassment, Moore et al. 1996, 1998) are often referred to as one of the mechanisms behind the morphological transformation of galaxies in overdense environments (Moore et al. 1999; Moss & Whittle 2000). In extreme cases, such processes may cause the partial removal of the gas reservoir of a galaxy through tidal stripping (e.g., Spilker et al. 2022). The efficiency of harassment mainly depends on the geometry of the infalling satellite’s orbit, the closest distance to the center of the potential, the cumulative number of fly-bys, and the intrinsic properties of the interacting galaxies (e.g., Gnedin 2003; Mastropietro et al. 2005; Smith et al. 2010; Chang et al. 2013; Bialas et al. 2015). Figs. 6 and 9 (left panel) show that the low gas fraction population in this work is dominated by massive galaxies (log M*/M ≥ 10.5), spanning a large variety of projected clustercentric distances 0 ≲ R/R200 ≲ 3 with a high fraction of AGN candidates (12/20). While harassment is capable of inducing instabilities in the disk funneling gas towards the nucleus thus triggering AGN activity, the widespread distribution of the X-ray sources across the protocluster global and local density regimes (see also Pérez-Martínez et al. 2023) argue against tidal stripping or harassment being an efficient mechanism to deplete the gas reservoir of the studied sample.

5.3. Effects of AGN feedback

Numerous observational and numerical studies have shown that gas depletion and subsequent quenching can also be a consequence of AGN feedback, expelling to the IGM part of the gas near the galaxy center that would otherwise be available for star formation while heating the remaining one (see Heckman & Best 2014 and references therein). Simulations predict that violent AGN outflows would be particularly efficient in removing gas from low-mass galaxies in very short timescales (< 100 Myr), while this process would be inefficient in their massive counterparts, impacting only the galaxy’s inner few kpc (Guo et al. 2022). Furthermore, Oppenheimer & Davé (2008) suggested that in this latter case, the large gravitational potential of the massive objects would be able to rebound (part of) the outflowing gas which will return to the galaxy after cooling down to be eventually recycled. Several authors have argued that a nascent ICM such as the one reported by Di Mascolo et al. (2023) in the Spiderweb protocluster would contribute to this scenario by impeding the outflowing gas leaving the galaxies, which would rain down over their disks contributing to enhancing their gas-phase metallicity (e.g., Kulas et al. 2013; Shimakawa et al. 2018a; Pérez-Martínez et al. 2023). Gas depletion throughout violent outflows is thus unlikely to be the only (or the main) driver acting on the depletion of the molecular gas reservoir in our subsample of HAEs hosting AGN, as they largely trace the intermediate-to-high mass regime (log M*/M > 10.5) for fgas < 0.7 (see Fig. 6).

Some hydrodynamical simulations suggest that the gradual depletion of cold gas is caused by the cumulative energy release from AGN feedback which, acting over sufficiently long timescales, end up heating the medium, suppressing cold accretion, and shutting down star formation (Ma et al. 2022; Ghodsi et al. 2024). In this scenario, black hole mass (MBH) would be an accurate predictor for gas depletion due to AGN feedback as it traces the black hole’s growth history (Piotrowska et al. 2022) while current activity indicators such as the Eddington ratio would play a secondary role except for extreme growth episodes (e.g., Guo et al. 2022). Interestingly, Zinger et al. (2020) found that a black hole mass scale of MBH = 107.5 − 108.5 is linked to drastic changes in several key physical parameters (e.g., cooling times, gas fraction and star formation) driving the transformation between low-mass, gas-rich, star-forming galaxies to massive gas-depleted quiescent using the IllustrisTNG project suite of simulations (Springel 2010; Pillepich et al. 2018). A similar MBH scale was found by Ma et al. (2022) and Ghodsi et al. (2024) when examining the depletion of the molecular gas reservoir and quenching of star formation for both IllustrisTNG and SIMBA (Davé et al. 2016, 2019) or by Bower et al. (2017) using the EAGLE simulation (Schaye et al. 2015). This MBH roughly corresponds with a stellar mass of M ≈ 1010.5 M assuming an almost invariant stellar to black hole mass relation at 0 < z < 2.0 (Reines & Volonteri 2015; Suh et al. 2020; Tanaka et al. 2024). This stellar mass threshold and sharp transition of properties resembles our Fig. 6 results where the molecular gas fraction estimates of our HAEs drops from fgas ≈ 0.9 at M ≲ 1010.5 M to fgas ≈ 0.3 at the massive end. Furthermore, 12 out of the 20 objects at M ≳ 1010.5 M show signs of current AGN activity throughout X-ray emission or emission line diagnostics, reinforcing the idea that the abrupt change of gas properties may be linked to AGN activity. This possibility is also discussed by Kolwa et al. (2023) for a sample of field massive (M ≳ 1011 M) and radio-loud AGN host at z = 2.9 − 4.3. These authors find high levels of gas depletion (fgas < 0.2) through the [C I](1−0) molecular gas tracer in contrast with star-forming or SMG coeval field samples. Nevertheless, this comparison between observational and simulated results might be limited by several factors. For example, the various specific implementations of AGN feedback between theoretical works (Hopkins 2015; Weinberger et al. 2017), which are tuned to be the main quenching mechanism of massive objects in SIMBA and the TNG simulations. In addition, several of the aforementioned works focus on the evolution of central galaxies at z = 0, thus minimizing the impact of additional environmental effects and neglecting intrinsic differences between the star-forming population of the local Universe and the Cosmic Noon.

On the other hand, Kimmig et al. (2025) proposes a somewhat different scenario using the Magneticum Pathfinder simulation suite. In this case, the best predictor for the quenching of massive high redshift galaxies would depend on the timescale scrutinized. On short timescales (< 500 Myr) galaxies may be quenched by the combined feedback of a violent starburst followed by a powerful episode of AGN activity. These events would be mediated either by strong accretion or gravitational interactions, able to trigger star formation and funnel gas to the central engine simultaneously. In this case, the best predictor for quenching is the black hole to stellar mass ratio while the correlation with the environment is subdominant. Thus, the population of recently quenched or soon-to-be-quenched massive satellite galaxies identified in z > 2 protoclusters during the last years (Kubo et al. 2021, 2022; McConachie et al. 2022; Ito et al. 2023; Jin et al. 2024; Shi et al. 2024b) may have predefined their main quenching mechanism independently of the environment they reside in at the time of their discovery. However, on quenching timescales longer than 1 Gyr the environment becomes the best predictor for quenching as the suppression of cold gas accretion eventually shuts down star formation (overconsumption, McGee et al. 2014) and prevents future episodes of rejuvenation as discussed by Remus & Kimmig (2023) using the same suit of simulations. Interestingly, this environmentally driven lack of accretion has different origins for massive central and satellite galaxies. Massive central quenched galaxies end up residing in underdense regions and are unable to accrete enough fuel to overcome the heating of their haloes and rekindle. However, the gas accretion of massive satellite galaxies is stopped due to the growth of the nascent ICM within the overdensity they reside in, and due to the onset of classic environmental effects such as RPS.

Our results in the Spiderweb protocluster do not show a synchronous decrease in SFR as the molecular gas reservoir is depleted (Figures 5 and 6) in contrast with the expectations from Zinger et al. (2020) and Ma et al. (2022). At the same time, we do not detect ongoing starburst activity in our massive galaxies based on the available Hα-based SFR, albeit we do see signs of AGN activity for a significant number of them. One possibility is that the ISM conditions for star formation may remain relatively stable while the gas reservoir is exhausted over time. In this scenario, quenching would be the consequence of overconsumption and thus felt by the galaxy only after the gas fraction reaches a critically low value. During this process, gas heating through previous and ongoing stellar and AGN feedback would prevent the incoming gas from cooling (Li et al. 2015, 2017; Dubois et al. 2013), thus ensuring the gradual depletion of the molecular gas reservoir. Taking as a reference that the transition phase in fgas happens for galaxies at 10.5 < log M*/M < 11 and that the typical SFR for a main sequence galaxy at log M*/M = 10.5 is SFR ≈ 30 − 100 M/yr we obtain that the time required to bridge this mass growth by in-situ star formation oscillates between 0.7 and 2.3 Gyr, which is shorter but not far from the typical depletion times measured in our sample (Fig. 8, τdep = 1 − 3 Gyr). Naturally, further environmental processes and their interrelations may extend or shorten this timescale. For example, as the Spiderweb protocluster continues growing its nascent ICM (Di Mascolo et al. 2023), RPS may become a more relevant and widespread phenomenon, which coupled with AGN and stellar feedback episodes could rapidly remove the gas reservoir of some satellite galaxies (e.g., Bahé & McCarthy 2015; Noble et al. 2019). However, its efficiency would depend on the stellar mass of the galaxy, geometrical factors, and the local physical conditions on its orbit through the ICM (see Boselli et al. 2022 for a review). At the same time, mergers could both accelerate the mass growth by bringing external material or inhibit it by directly expelling gas or triggering AGN activity depending on the specific gas and stellar mass properties of the satellite galaxies involved (Pontzen et al. 2017; Sparre & Springel 2017; Lagos et al. 2018).

Nevertheless, our SFRs rely on narrow-band and spectroscopic Hα flux measurements carried out by Shimakawa et al. (2018a) and Pérez-Martínez et al. (2023) respectively. This tracer is unable to capture heavily obscured star formation episodes common in SMGs. However, only 5 objects within our sample are identified as such (Dannerbauer et al. 2014) and they reside in the protocluster outskirts. In addition, these previous works assumed a star-forming origin as the main ionization source for Hα. This assumption is somehow reinforced by the low fraction of broad-line emitters identified by Pérez-Martínez et al. (2023), with only two broad-line AGN out of 39 spectroscopically analyzed HAEs. However, the X-ray identifications reported in Tozzi et al. (2022a) and the multi-wavelength SED analysis carried out by Shimakawa et al. (2024b) over these same sources suggest that a fraction of the Hα flux of these sources may be contaminated by the ionization from the central engine which could result in the overestimation of our SFR for this sources, thus placing them up to ∼1 dex below the main sequence of star formation and in qualitative agreement with the predictions from the numerical simulations (Zinger et al. 2020).

Regardless of AGN feedback being (or not) the dominant factor for the cold gas depletion and quenching of star formation at the massive end of our sample, it remains unclear if the cumulative growth of the supermassive black hole or the triggering event for its current activity are linked to environmental factors such as gas angular momentum loss caused by enhanced galaxy-galaxy interactions in protoclusters (e.g., Lotz et al. 2013; Hine et al. 2016; Watson et al. 2019; Liu et al. 2023; Naufal et al. 2023; Shibuya et al. 2024) or it rather merely responds to secular processes (i.e., enhanced accretion) acting over the assembly history of massive galaxies which are statistically amplified by the fact that protoclusters tend to host larger number counts of massive objects than the coeval field (e.g., Nantais et al. 2016; Shimakawa et al. 2018b; Hill et al. 2022). During the last years, several teams hinted at diverse quenching mechanisms for passively evolving galaxies residing within protoclusters albeit current observational evidence does not allow to pinpoint a purely environmental origin (e.g., Kalita et al. 2021; Kubo et al. 2021; McConachie et al. 2022; Ito et al. 2023; Shi et al. 2024b; Jin et al. 2024). In the Spiderweb protocluster, AGNs are widely distributed in terms of phase-space (Fig. 9) as well as projected local densities (Fig. 11, see also Pérez-Martínez et al. 2023) making it difficult to establish a causal link between their local environment and their ongoing AGN activity. Furthermore, current broad-band photometry in the NIR lacks the resolution to identify tidal tails or heavily disturbed morphologies that could act as the smoking gun for recent gravitational interactions (Jin et al. 2024).

6. Conclusions

In this work, we have analyzed the molecular gas properties of a sample of HAEs within the Spiderweb protocluster using CO(1−0) measurements and upper limits from the COALAS survey as a proxy to weigh their cold gas reservoir. We have examined fgas,  τdep, SFR, and SFE in our sample and put in compare our results with coeval scaling relations and observational works both in the field and protoclusters (Sect. 4). Finally, we have discussed the implications of these results on different scenarios of galaxy evolution during the early stages of protocluster assembly (Sect. 5). In the following paragraphs, we summarize the main conclusions of this work:

  1. We have identified 43 spectroscopically confirmed HAEs within the COALAS ATCA CO(1−0) footprint and obtained CO(1−0) detection for ten objects in agreement with Jin et al. (2021). In addition, we have obtained 26 CO(1−0) upper flux limits for the rest of our sample.

  2. Our protocluster sample displays a very sharp transition from high to low molecular gas fractions (fgas) values at stellar masses log M*/M = 10.5 − 11.0, in agreement with previous coeval observational CO(1−0) works in protoclusters (Wang et al. 2018) and the field (e.g., Riechers et al. 2020). The distribution of these sources in the fgas − M* plane can be fitted by a logistic function (Eq. 5) as suggested by Popping et al. (2012). This is in contrast with the smoother functions reported by Sargent et al. (2014) and Tacconi et al. (2018).

  3. More than half (12/20) of the HAEs with log M*/M ≥ 10.5 host X-ray (Tozzi et al. 2022a) or emission line diagnostics indicating the presence of an AGN (Pérez-Martínez et al. 2023). Furthermore, the stacking of these sources reveals a lower gas fraction value (Fgas = 0.38 ± 0.23) than that of equally massive sources but without signs of AGN activity (Fgas = 0.55 ± 0.12). This hints that nuclear activity contributes to molecular gas depletion beyond this mass limit.

  4. The measured depletion times (τdep) show that, in the absence of inflows, most of these HAEs will consume their gas reservoir via star formation in 1 − 3 Gyr (i.e., by 1 < z < 1.6), concurrently with the observational evidence suggesting that massive galaxy clusters establish their red sequence at 1.5 < z ≲ 2 which will end up dominating their cores by z = 1. Environmental effects during this period could further shorten the depletion time while cold gas accretion could delay it.

  5. Our environmental analysis suggests that the fraction of HAEs with higher molecular-to-stellar mass ratios (μgas) increases towards the outskirts of the cluster (Fig. 9). The opposite trend is found in terms of SFE, similarly to Wang et al. (2018). However, SFR shows no clear correlation as a function of phase space (η) in this protocluster (Fig. 10, see also Pérez-Martínez et al. 2023). This implies that the changes in SFE are predominately driven by decreasing molecular gas reservoirs while star formation activity proceeds unaffected for most galaxies regardless of their location within the protocluster structure.

  6. We propose three possible mechanisms to explain the depletion of the molecular gas reservoir: changes in gas accretion due to the halo growth, gas removal through environmental effects, and AGN activity. The sharp change in gas fraction at log M*/M = 10.5 − 11.0 suggests that the dominant mechanism predominately acts in the high stellar mass regime. Thus, we propose AGN gas heating coupled with overconsumption as the likely mechanism behind this behavior. As the nascent ICM grows, the efficiency of gas accretion diminishes. On top of this, AGN heating contributes to the confinement of the host galaxy while allowing its star formation activities to continue in a relatively unaltered way. As a result, the galaxy will end up depleting its gas reservoir and quenching. This scenario is supported by the high AGN fraction between massive galaxies in this protocluster. On the other hand gas removal would have a secondary role in this stellar mass range, as it is typically more efficient on lower-mass galaxies (e.g., ram pressure or tidal stripping) at later cosmic epochs when clusters have fully virialized and developed their ICM.

Overall, this work has conducted one of the most comprehensive studies of star formation and gas consumption in a protocluster during the cosmic noon despite the spatial resolution constraints of both ATCA submillimetre and NIR ground-based observations. The recent advent of ALMA Band-1 in Cycle 10 provides new opportunities to examine and replicate our findings using CO(1−0) at higher spatial resolution and across protoclusters at different evolutionary stages during this epoch. Furthermore, the integration of ALMA with the exceptional NIR and MIR capabilities of the James Webb Space Telescope (JWST) will enable a deeper understanding of some of the remaining challenges for galaxy evolution in early overdense environments in the near future. These include determining when, where, and how intensely star-forming activities occur within protocluster members, and what roles the environment plays in enhancing or supressing them during the early stages of cluster assembly.


Acknowledgments

The Australia Telescope Compact Array is part of the Australia Telescope National Facility (grid.421683.a) which is funded by the Australian Government for operation as a National Facility managed by CSIRO. We acknowledge the Gomeroi people as the traditional owners of the Observatory site. This research is based in part on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 095.A-0500(B). This research is based in part on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan (NAOJ). We are honored and grateful for the opportunity of observing the Universe from Maunakea, which has cultural, historical, and natural significance in Hawaii. This work is based in part on observations made with the Spitzer Space Telescope, which was operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. This research is based in part on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. The National Radio Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This research made use of Astropy (http://www.astropy.org), a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018). JMPM acknowledges funding from the European Union’s Horizon-Europe research and innovation program under the Marie Skłodowska-Curie grant agreement No 101106626. HD, ZC, YZ, and JMPM acknowledge financial support from the Agencia Estatal de Investigación del Ministerio de Ciencia e Innovación (AEI-MCINN) under grant (La evolución de los cúmulos de galaxias desde el amanecer hasta el mediodía cósmico) with reference (PID2019-105776GB-I00/DOI:10.13039/501100011033). HD, YZ, and JMPM acknowledge financial support from the Ministerio de Ciencia, Innovación y Universidades (MCIU/AEI) under grant (Construcción de cúmulos de galaxias en formación a través de la formación estelar oscurecida por el polvo) and the European Regional Development Fund (ERDF) with reference (PID2022-143243NB-I00/DOI:10.13039/501100011033). C.D.E. acknowledges funding from the MCIN/AEI (Spain) and the “NextGenerationEU”/PRTR (European Union) through the Juan de la Cierva-Formación program (FJC2021-047307-I). SJ acknowledges the support from the European Union’s Horizon Europe research and innovation program under the Marie Skłodowska-Curie Action grant No. 101060888. The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant DNRF140. M.S.P. acknowledges the support of the Spanish Ministry of Science, Innovation and Universities through the project PID-2021-122544NB-C43. YZ acknowledges the support from the China Scholarship Council (202206340048). This article is based upon work from COST Action CA21126 – Carbon molecular nanostructures in space (NanoSpace), supported by COST (European Cooperation in Science and Technology).

References

  1. Alberts, S., & Noble, A. 2022, Universe, 8, 554 [Google Scholar]
  2. Alberts, S., Pope, A., Brodwin, M., et al. 2014, MNRAS, 437, 437 [NASA ADS] [CrossRef] [Google Scholar]
  3. Anderson, C. S., Carilli, C. L., Tozzi, P., et al. 2022, ApJ, 937, 45 [NASA ADS] [CrossRef] [Google Scholar]
  4. Aoyama, K., Kodama, T., Suzuki, T. L., et al. 2022, ApJ, 924, 74 [NASA ADS] [CrossRef] [Google Scholar]
  5. Aravena, M., Decarli, R., Gónzalez-López, J., et al. 2019, ApJ, 882, 136 [Google Scholar]
  6. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  8. Bahé, Y. M., & McCarthy, I. G. 2015, MNRAS, 447, 969 [Google Scholar]
  9. Barro, G., Faber, S. M., Dekel, A., et al. 2016, ApJ, 820, 120 [Google Scholar]
  10. Bassini, L., Feldmann, R., Gensior, J., et al. 2023, MNRAS, 525, 5388 [NASA ADS] [CrossRef] [Google Scholar]
  11. Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 2019, MNRAS, 488, 3143 [NASA ADS] [CrossRef] [Google Scholar]
  12. Beifiori, A., Mendel, J. T., Chan, J. C. C., et al. 2017, ApJ, 846, 120 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Bialas, D., Lisker, T., Olczak, C., Spurzem, R., & Kotulla, R. 2015, A&A, 576, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Birkin, J. E., Weiss, A., Wardlow, J. L., et al. 2021, MNRAS, 501, 3926 [Google Scholar]
  16. Blánquez-Sesé, D., Gómez-Guijarro, C., Magdis, G. E., et al. 2023, A&A, 674, A166 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Bleem, L. E., Stalder, B., de Haan, T., et al. 2015, ApJS, 216, 27 [Google Scholar]
  18. Böhringer, H., Schuecker, P., Guzzo, L., et al. 2004, A&A, 425, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Böhringer, H., Schuecker, P., Pratt, G. W., et al. 2007, A&A, 469, 363 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 [CrossRef] [Google Scholar]
  21. Boogaard, L. A., Decarli, R., González-López, J., et al. 2019, ApJ, 882, 140 [NASA ADS] [CrossRef] [Google Scholar]
  22. Boogaard, L. A., van der Werf, P., Weiss, A., et al. 2020, ApJ, 902, 109 [NASA ADS] [CrossRef] [Google Scholar]
  23. Boquien, M., Burgarella, D., Roehlly, Y., et al. 2019, A&A, 622, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Boselli, A., Cuillandre, J. C., Fossati, M., et al. 2016, A&A, 587, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Boselli, A., Fossati, M., & Sun, M. 2022, A&ARv, 30, 3 [NASA ADS] [CrossRef] [Google Scholar]
  26. Bower, R. G., Schaye, J., Frenk, C. S., et al. 2017, MNRAS, 465, 32 [NASA ADS] [CrossRef] [Google Scholar]
  27. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  28. Burgarella, D., Buat, V., & Iglesias-Páramo, J. 2005, MNRAS, 360, 1413 [NASA ADS] [CrossRef] [Google Scholar]
  29. Calistro Rivera, G., Hodge, J. A., Smail, I., et al. 2018, ApJ, 863, 56 [Google Scholar]
  30. Calvi, R., Castignani, G., & Dannerbauer, H. 2023, A&A, 678, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  32. Carilli, C. L., Daddi, E., Riechers, D., et al. 2010, ApJ, 714, 1407 [NASA ADS] [CrossRef] [Google Scholar]
  33. Carilli, C. L., Anderson, C. S., Tozzi, P., et al. 2022, ApJ, 928, 59 [NASA ADS] [CrossRef] [Google Scholar]
  34. Casey, C. M. 2016, ApJ, 824, 36 [CrossRef] [Google Scholar]
  35. Casey, C. M., Cooray, A., Killi, M., et al. 2017, ApJ, 840, 101 [NASA ADS] [CrossRef] [Google Scholar]
  36. Chabrier, G. 2003, PASP, 115, 763 [Google Scholar]
  37. Champagne, J. B., Casey, C. M., Zavala, J. A., et al. 2021, ApJ, 913, 110 [NASA ADS] [CrossRef] [Google Scholar]
  38. Chang, J., Macciò, A. V., & Kang, X. 2013, MNRAS, 431, 3533 [NASA ADS] [CrossRef] [Google Scholar]
  39. Chartab, N., Mobasher, B., Shapley, A. E., et al. 2021, ApJ, 908, 120 [CrossRef] [Google Scholar]
  40. Chen, Z., Dannerbauer, H., Lehnert, M. D., et al. 2024, MNRAS, 527, 8950 [Google Scholar]
  41. Chiang, Y.-K., Overzier, R. A., Gebhardt, K., & Henriques, B. 2017, ApJ, 844, L23 [Google Scholar]
  42. Christensen, C. R., Davé, R., Governato, F., et al. 2016, ApJ, 824, 57 [CrossRef] [Google Scholar]
  43. Contini, E., De Lucia, G., & Borgani, S. 2012, MNRAS, 420, 2978 [NASA ADS] [CrossRef] [Google Scholar]
  44. Coogan, R. T., Daddi, E., Sargent, M. T., et al. 2018, MNRAS, 479, 703 [NASA ADS] [Google Scholar]
  45. Croft, S., Kurk, J., van Breugel, W., et al. 2005, AJ, 130, 867 [NASA ADS] [CrossRef] [Google Scholar]
  46. Cucciati, O., Zamorani, G., Lemaux, B. C., et al. 2014, A&A, 570, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Cucciati, O., Lemaux, B. C., Zamorani, G., et al. 2018, A&A, 619, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Cullen, F., Cirasuolo, M., Kewley, L. J., et al. 2016, MNRAS, 460, 3002 [NASA ADS] [CrossRef] [Google Scholar]
  49. Daddi, E., Rich, R. M., Valentino, F., et al. 2022a, ApJ, 926, L21 [NASA ADS] [CrossRef] [Google Scholar]
  50. Daddi, E., Delvecchio, I., Dimauro, P., et al. 2022b, A&A, 661, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Daikuhara, K., Kodama, T., Pérez-Martínez, J. M., et al. 2024, MNRAS, 531, 2335 [NASA ADS] [CrossRef] [Google Scholar]
  52. Dannerbauer, H., Daddi, E., Riechers, D. A., et al. 2009, ApJ, 698, L178 [NASA ADS] [CrossRef] [Google Scholar]
  53. Dannerbauer, H., Kurk, J. D., De Breuck, C., et al. 2014, A&A, 570, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Dannerbauer, H., Lehnert, M. D., Emonts, B., et al. 2017, A&A, 608, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Davé, R., Thompson, R., & Hopkins, P. F. 2016, MNRAS, 462, 3265 [Google Scholar]
  56. Davé, R., Anglés-Alcázar, D., Narayanan, D., et al. 2019, MNRAS, 486, 2827 [Google Scholar]
  57. Davé, R., Crain, R. A., Stevens, A. R. H., et al. 2020, MNRAS, 497, 146 [Google Scholar]
  58. De Breuck, C., Lundgren, A., Emonts, B., et al. 2022, A&A, 658, L2 [CrossRef] [EDP Sciences] [Google Scholar]
  59. Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 [NASA ADS] [CrossRef] [Google Scholar]
  60. Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451 [Google Scholar]
  61. Dekel, A., Zolotov, A., Tweed, D., et al. 2013, MNRAS, 435, 999 [Google Scholar]
  62. D’Eugenio, C., Daddi, E., Liu, D., & Gobat, R. 2023, A&A, 678, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Dickman, R. L., Snell, R. L., & Schloerb, F. P. 1986, ApJ, 309, 326 [NASA ADS] [CrossRef] [Google Scholar]
  64. Di Mascolo, L., Saro, A., Mroczkowski, T., et al. 2023, Nature, 615, 809 [NASA ADS] [CrossRef] [Google Scholar]
  65. Doherty, M., Tanaka, M., De Breuck, C., et al. 2010, A&A, 509, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Dubois, Y., Gavazzi, R., Peirani, S., & Silk, J. 2013, MNRAS, 433, 3297 [CrossRef] [Google Scholar]
  67. Ehlert, S., von der Linden, A., Allen, S. W., et al. 2014, MNRAS, 437, 1942 [CrossRef] [Google Scholar]
  68. El-Badry, K., Wetzel, A., Geha, M., et al. 2016, ApJ, 820, 131 [Google Scholar]
  69. Emonts, B. H. C., Norris, R. P., Feain, I., et al. 2011, MNRAS, 415, 655 [Google Scholar]
  70. Emonts, B. H. C., Feain, I., Röttgering, H. J. A., et al. 2013, MNRAS, 430, 3465 [NASA ADS] [CrossRef] [Google Scholar]
  71. Emonts, B. H. C., Lehnert, M. D., Villar-Martín, M., et al. 2016, Science, 354, 1128 [NASA ADS] [CrossRef] [Google Scholar]
  72. Emonts, B. H. C., Lehnert, M. D., Dannerbauer, H., et al. 2018, MNRAS, 477, L60 [NASA ADS] [CrossRef] [Google Scholar]
  73. Falkendal, T., De Breuck, C., Lehnert, M. D., et al. 2019, A&A, 621, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Fassbender, R., Nastasi, A., Santos, J. S., et al. 2014, A&A, 568, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Feldmann, R., & Mayer, L. 2015, MNRAS, 446, 1939 [Google Scholar]
  76. Forrest, B., Lemaux, B. C., Shah, E., et al. 2023, MNRAS, 526, L56 [NASA ADS] [CrossRef] [Google Scholar]
  77. Frias Castillo, M., Hodge, J., Rybak, M., et al. 2023, ApJ, 945, 128 [NASA ADS] [CrossRef] [Google Scholar]
  78. Fumagalli, M., Fossati, M., Hau, G. K. T., et al. 2014, MNRAS, 445, 4335 [Google Scholar]
  79. Genel, S., Genzel, R., Bouché, N., et al. 2008, ApJ, 688, 789 [NASA ADS] [CrossRef] [Google Scholar]
  80. Genzel, R., Förster Schreiber, N. M., Rosario, D., et al. 2014, ApJ, 796, 7 [Google Scholar]
  81. Genzel, R., Tacconi, L. J., Lutz, D., et al. 2015, ApJ, 800, 20 [Google Scholar]
  82. Ghodsi, L., Man, A. W. S., Donevski, D., et al. 2024, MNRAS, 528, 4393 [NASA ADS] [CrossRef] [Google Scholar]
  83. Gnedin, O. Y. 2003, ApJ, 582, 141 [Google Scholar]
  84. Goerdt, T., Dekel, A., Sternberg, A., et al. 2010, MNRAS, 407, 613 [Google Scholar]
  85. Gómez-Guijarro, C., Magdis, G. E., Valentino, F., et al. 2019, ApJ, 886, 88 [Google Scholar]
  86. Grützbauch, R., Bauer, A. E., Jørgensen, I., & Varela, J. 2012, MNRAS, 423, 3652 [Google Scholar]
  87. Gunn, J. E., & Gott, J. R. 1972, ApJ, 176, 1 [NASA ADS] [CrossRef] [Google Scholar]
  88. Guo, H., Jones, M. G., & Wang, J. 2022, ApJ, 933, L12 [NASA ADS] [CrossRef] [Google Scholar]
  89. Haines, C. P., Pereira, M. J., Sanderson, A. J. R., et al. 2012, ApJ, 754, 97 [NASA ADS] [CrossRef] [Google Scholar]
  90. Hatch, N. A., Overzier, R. A., Röttgering, H. J. A., Kurk, J. D., & Miley, G. K. 2008, MNRAS, 383, 931 [Google Scholar]
  91. Hatch, N. A., Overzier, R. A., Kurk, J. D., et al. 2009, MNRAS, 395, 114 [NASA ADS] [CrossRef] [Google Scholar]
  92. Hatch, N. A., Kurk, J. D., Pentericci, L., et al. 2011, MNRAS, 415, 2993 [NASA ADS] [CrossRef] [Google Scholar]
  93. Heckman, T. M., & Best, P. N. 2014, ARA&A, 52, 589 [Google Scholar]
  94. Hill, R., Chapman, S., Phadke, K. A., et al. 2022, MNRAS, 512, 4352 [NASA ADS] [CrossRef] [Google Scholar]
  95. Hine, N. K., Geach, J. E., Alexander, D. M., et al. 2016, MNRAS, 455, 2363 [NASA ADS] [CrossRef] [Google Scholar]
  96. Hodge, J. A., Carilli, C. L., Walter, F., et al. 2012, ApJ, 760, 11 [Google Scholar]
  97. Hopkins, P. F. 2015, MNRAS, 450, 53 [Google Scholar]
  98. Huang, R., Battisti, A. J., Grasha, K., et al. 2023, MNRAS, 520, 446 [NASA ADS] [CrossRef] [Google Scholar]
  99. Ikeda, R., Tadaki, K.-I., Iono, D., et al. 2022, ApJ, 933, 11 [NASA ADS] [CrossRef] [Google Scholar]
  100. Ito, K., Tanaka, M., Valentino, F., et al. 2023, ApJ, 945, L9 [NASA ADS] [CrossRef] [Google Scholar]
  101. Ivison, R. J., Papadopoulos, P. P., Smail, I., et al. 2011, MNRAS, 412, 1913 [Google Scholar]
  102. Ivison, R. J., Swinbank, A. M., Smail, I., et al. 2013, ApJ, 772, 137 [NASA ADS] [CrossRef] [Google Scholar]
  103. Jin, S., Dannerbauer, H., Emonts, B., et al. 2021, A&A, 652, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Jin, S., Sillassen, N. B., Magdis, G. E., et al. 2024, A&A, 683, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Kaasinen, M., Scoville, N., Walter, F., et al. 2019, ApJ, 880, 15 [NASA ADS] [CrossRef] [Google Scholar]
  106. Kacprzak, G. G., Yuan, T., Nanayakkara, T., et al. 2015, ApJ, 802, L26 [NASA ADS] [CrossRef] [Google Scholar]
  107. Kalita, B. S., Daddi, E., D’Eugenio, C., et al. 2021, ApJ, 917, L17 [NASA ADS] [CrossRef] [Google Scholar]
  108. Keller, B. W., Wadsley, J., & Couchman, H. M. P. 2016, MNRAS, 463, 1431 [Google Scholar]
  109. Kennicutt, R. C. 1998, ARA&A, 36, 189 [NASA ADS] [CrossRef] [Google Scholar]
  110. Kimmig, L. C., Remus, R.-S., Seidel, B., et al. 2025, ApJ, 979, 15 [Google Scholar]
  111. Kodama, T., Tanaka, I., Kajisawa, M., et al. 2007, MNRAS, 377, 1717 [NASA ADS] [CrossRef] [Google Scholar]
  112. Kodama, T., Hayashi, M., Koyama, Y., et al. 2013, Proc. Int. Astron. Union, 295, 74 [Google Scholar]
  113. Kolwa, S., De Breuck, C., Vernet, J., et al. 2023, MNRAS, 525, 5831 [NASA ADS] [CrossRef] [Google Scholar]
  114. Koulouridis, E., & Bartalucci, I. 2019, A&A, 623, L10 [EDP Sciences] [Google Scholar]
  115. Koulouridis, E., Gkini, A., & Drigga, E. 2024, A&A, 684, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Koyama, Y., Kodama, T., Tadaki, K.-I., et al. 2013, MNRAS, 428, 1551 [NASA ADS] [CrossRef] [Google Scholar]
  117. Krishnan, C., Hatch, N. A., Almaini, O., et al. 2017, MNRAS, 470, 2170 [NASA ADS] [CrossRef] [Google Scholar]
  118. Kubo, M., Umehata, H., Matsuda, Y., et al. 2021, ApJ, 919, 6 [NASA ADS] [CrossRef] [Google Scholar]
  119. Kubo, M., Umehata, H., Matsuda, Y., et al. 2022, ApJ, 935, 89 [NASA ADS] [CrossRef] [Google Scholar]
  120. Kuiper, E., Hatch, N. A., Miley, G. K., et al. 2011, MNRAS, 415, 2245 [NASA ADS] [CrossRef] [Google Scholar]
  121. Kulas, K. R., McLean, I. S., Shapley, A. E., et al. 2013, ApJ, 774, 130 [NASA ADS] [CrossRef] [Google Scholar]
  122. Kurk, J. D., Röttgering, H. J. A., Pentericci, L., et al. 2000, A&A, 358, L1 [NASA ADS] [Google Scholar]
  123. Kurk, J. D., Pentericci, L., Overzier, R. A., Röttgering, H. J. A., & Miley, G. K. 2004, A&A, 428, 817 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Laganá, T. F., Martinet, N., Durret, F., et al. 2013, A&A, 555, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  125. Lagos, C. D. P., Stevens, A. R. H., Bower, R. G., et al. 2018, MNRAS, 473, 4956 [NASA ADS] [CrossRef] [Google Scholar]
  126. Laishram, R., Kodama, T., Morishita, T., et al. 2024, ApJ, 964, L33 [Google Scholar]
  127. Larson, R. B. 1974, MNRAS, 169, 229 [NASA ADS] [CrossRef] [Google Scholar]
  128. Lee, M. M., Tanaka, I., Kawabe, R., et al. 2017, ApJ, 842, 55 [CrossRef] [Google Scholar]
  129. Lehmer, B. D., Alexander, D. M., Geach, J. E., et al. 2009, ApJ, 691, 687 [Google Scholar]
  130. Lehmer, B. D., Lucy, A. B., Alexander, D. M., et al. 2013, ApJ, 765, 87 [NASA ADS] [CrossRef] [Google Scholar]
  131. Lemaux, B. C., Cucciati, O., Le Fèvre, O., et al. 2022, A&A, 662, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  132. Lepore, M., Di Mascolo, L., Tozzi, P., et al. 2024, A&A, 682, A186 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  133. Li, Y., Bryan, G. L., Ruszkowski, M., et al. 2015, ApJ, 811, 73 [Google Scholar]
  134. Li, Y., Ruszkowski, M., & Bryan, G. L. 2017, ApJ, 847, 106 [Google Scholar]
  135. Liao, C.-L., Chen, C.-C., Wang, W.-H., et al. 2024, ApJ, 961, 226 [NASA ADS] [CrossRef] [Google Scholar]
  136. Liu, S., Zheng, X. Z., Shi, D. D., et al. 2023, MNRAS, 523, 2422 [CrossRef] [Google Scholar]
  137. Lotz, J. M., Papovich, C., Faber, S. M., et al. 2013, ApJ, 773, 154 [NASA ADS] [CrossRef] [Google Scholar]
  138. Ma, W., Liu, K., Guo, H., et al. 2022, ApJ, 941, 205 [NASA ADS] [CrossRef] [Google Scholar]
  139. Macuga, M., Martini, P., Miller, E. D., et al. 2019, ApJ, 874, 54 [NASA ADS] [CrossRef] [Google Scholar]
  140. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [Google Scholar]
  141. Mao, Z., Kodama, T., Pérez-Martínez, J. M., et al. 2022, A&A, 666, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  142. Mastropietro, C., Moore, B., Mayer, L., et al. 2005, MNRAS, 364, 607 [Google Scholar]
  143. McConachie, I., Wilson, G., Forrest, B., et al. 2022, ApJ, 926, 37 [NASA ADS] [CrossRef] [Google Scholar]
  144. McGee, S. L., Bower, R. G., & Balogh, M. L. 2014, MNRAS, 442, L105 [Google Scholar]
  145. Mei, S., Hatch, N. A., Amodeo, S., et al. 2023, A&A, 670, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  146. Miley, G. K., Overzier, R. A., Zirm, A. W., et al. 2006, ApJ, 650, L29 [NASA ADS] [CrossRef] [Google Scholar]
  147. Monson, E. B., Lehmer, B. D., Doore, K., et al. 2021, ApJ, 919, 51 [NASA ADS] [CrossRef] [Google Scholar]
  148. Moore, B., Katz, N., Lake, G., Dressler, A., & Oemler, A. 1996, Nature, 379, 613 [Google Scholar]
  149. Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 [Google Scholar]
  150. Moore, B., Lake, G., Quinn, T., & Stadel, J. 1999, MNRAS, 304, 465 [NASA ADS] [CrossRef] [Google Scholar]
  151. Moss, C., & Whittle, M. 2000, MNRAS, 317, 667 [NASA ADS] [CrossRef] [Google Scholar]
  152. Moustakas, J., Kennicutt, R. C., & Tremonti, C. A. 2006, ApJ, 642, 775 [Google Scholar]
  153. Muldrew, S. I., Hatch, N. A., & Cooke, E. A. 2015, MNRAS, 452, 2528 [NASA ADS] [CrossRef] [Google Scholar]
  154. Namiki, S. V., Koyama, Y., Hayashi, M., et al. 2019, ApJ, 877, 118 [NASA ADS] [CrossRef] [Google Scholar]
  155. Nantais, J. B., van der Burg, R. F. J., Lidman, C., et al. 2016, A&A, 592, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  156. Narayanan, D., Bothwell, M., & Davé, R. 2012, MNRAS, 426, 1178 [Google Scholar]
  157. Naufal, A., Koyama, Y., Shimakawa, R., & Kodama, T. 2023, ApJ, 958, 170 [NASA ADS] [CrossRef] [Google Scholar]
  158. Naufal, A., Koyama, Y., D’Eugenio, C., et al. 2024, ApJ, 977, 58 [Google Scholar]
  159. Noble, A. G., Webb, T. M. A., Muzzin, A., et al. 2013, ApJ, 768, 118 [Google Scholar]
  160. Noble, A. G., Muzzin, A., McDonald, M., et al. 2019, ApJ, 870, 56 [Google Scholar]
  161. Noll, S., Burgarella, D., Giovannoli, E., et al. 2009, A&A, 507, 1793 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  162. Noordeh, E., Canning, R. E. A., Willis, J. P., et al. 2021, MNRAS, 507, 5272 [NASA ADS] [CrossRef] [Google Scholar]
  163. Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 [NASA ADS] [CrossRef] [Google Scholar]
  164. Oppenheimer, B. D., & Davé, R. 2008, MNRAS, 387, 577 [NASA ADS] [CrossRef] [Google Scholar]
  165. Overzier, R. A. 2016, A&ARv, 24, 14 [Google Scholar]
  166. Pavesi, R., Sharon, C. E., Riechers, D. A., et al. 2018, ApJ, 864, 49 [Google Scholar]
  167. Pensabene, A., Cantalupo, S., Cicone, C., et al. 2024, A&A, 684, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  168. Pentericci, L., Roettgering, H. J. A., Miley, G. K., Carilli, C. L., & McCarthy, P. 1997, A&A, 326, 580 [Google Scholar]
  169. Pentericci, L., Kurk, J. D., Röttgering, H. J. A., et al. 2000, A&A, 361, L25 [NASA ADS] [Google Scholar]
  170. Pérez-Martínez, J. M., Dannerbauer, H., Kodama, T., et al. 2023, MNRAS, 518, 1707 [Google Scholar]
  171. Pérez-Martínez, J. M., Kodama, T., Koyama, Y., et al. 2024a, MNRAS, 527, 10221 [Google Scholar]
  172. Pérez-Martínez, J. M., Dannerbauer, H., Koyama, Y., et al. 2024b, ApJ, 977, 74 [Google Scholar]
  173. Pillepich, A., Springel, V., Nelson, D., et al. 2018, MNRAS, 473, 4077 [Google Scholar]
  174. Piotrowska, J. M., Bluck, A. F. L., Maiolino, R., & Peng, Y. 2022, MNRAS, 512, 1052 [NASA ADS] [CrossRef] [Google Scholar]
  175. Polletta, M., Soucail, G., Dole, H., et al. 2021, A&A, 654, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  176. Polletta, M., Dole, H., Martinache, C., et al. 2022, A&A, 662, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  177. Pontzen, A., Tremmel, M., Roth, N., et al. 2017, MNRAS, 465, 547 [Google Scholar]
  178. Popping, G., Caputi, K. I., Somerville, R. S., & Trager, S. C. 2012, MNRAS, 425, 2386 [Google Scholar]
  179. Reines, A. E., & Volonteri, M. 2015, ApJ, 813, 82 [NASA ADS] [CrossRef] [Google Scholar]
  180. Remus, R. S., & Kimmig, L. C. 2023, ArXiv e-prints [arXiv:2310.16089] [Google Scholar]
  181. Remus, R.-S., Dolag, K., & Dannerbauer, H. 2023, ApJ, 950, 191 [NASA ADS] [CrossRef] [Google Scholar]
  182. Riechers, D. A., Boogaard, L. A., Decarli, R., et al. 2020, ApJ, 896, L21 [NASA ADS] [CrossRef] [Google Scholar]
  183. Riechers, D. A., Nayyeri, H., Burgarella, D., et al. 2021, ApJ, 907, 62 [NASA ADS] [CrossRef] [Google Scholar]
  184. Rigby, E. E., Hatch, N. A., Röttgering, H. J. A., et al. 2014, MNRAS, 437, 1882 [NASA ADS] [CrossRef] [Google Scholar]
  185. Roberts, I. D., Parker, L. C., Brown, T., et al. 2019, ApJ, 873, 42 [CrossRef] [Google Scholar]
  186. Roberts, I. D., Brown, T., Zabel, N., et al. 2023, A&A, 675, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  187. Rosati, P., Tozzi, P., Gobat, R., et al. 2009, A&A, 508, 583 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  188. Rudnick, G., Hodge, J., Walter, F., et al. 2017, ApJ, 849, 27 [NASA ADS] [CrossRef] [Google Scholar]
  189. Sanderson, A. J. R., Ponman, T. J., & O’Sullivan, E. 2006, MNRAS, 372, 1496 [NASA ADS] [CrossRef] [Google Scholar]
  190. Sanderson, A. J. R., O’Sullivan, E., & Ponman, T. J. 2009, MNRAS, 395, 764 [NASA ADS] [CrossRef] [Google Scholar]
  191. Sargent, M. T., Daddi, E., Béthermin, M., et al. 2014, ApJ, 793, 19 [NASA ADS] [CrossRef] [Google Scholar]
  192. Saro, A., Borgani, S., Tornatore, L., et al. 2009, MNRAS, 392, 795 [NASA ADS] [CrossRef] [Google Scholar]
  193. Sattari, Z., Mobasher, B., Chartab, N., et al. 2021, ApJ, 910, 57 [CrossRef] [Google Scholar]
  194. Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 [Google Scholar]
  195. Seymour, N., Stern, D., De Breuck, C., et al. 2007, ApJS, 171, 353 [CrossRef] [Google Scholar]
  196. Shapiro, P. R., & Field, G. B. 1976, ApJ, 205, 762 [NASA ADS] [CrossRef] [Google Scholar]
  197. Shi, K., Toshikawa, J., Lee, K.-S., et al. 2021, ApJ, 911, 46 [NASA ADS] [CrossRef] [Google Scholar]
  198. Shi, K., Malavasi, N., Toshikawa, J., & Zheng, X. 2024a, ApJ, 961, 39 [NASA ADS] [CrossRef] [Google Scholar]
  199. Shi, D. D., Wang, X., Zheng, X. Z., et al. 2024b, ApJ, 963, 21 [NASA ADS] [CrossRef] [Google Scholar]
  200. Shibuya, T., Ito, Y., Asai, K., et al. 2024, PASJ, 77, 21 [Google Scholar]
  201. Shimakawa, R., Kodama, T., Tadaki, K. I., et al. 2014, MNRAS, 441, L1 [NASA ADS] [CrossRef] [Google Scholar]
  202. Shimakawa, R., Kodama, T., Tadaki, K.-I., et al. 2015, MNRAS, 448, 666 [NASA ADS] [CrossRef] [Google Scholar]
  203. Shimakawa, R., Kodama, T., Hayashi, M., et al. 2018a, MNRAS, 473, 1977 [NASA ADS] [CrossRef] [Google Scholar]
  204. Shimakawa, R., Koyama, Y., Röttgering, H. J. A., et al. 2018b, MNRAS, 481, 5630 [NASA ADS] [CrossRef] [Google Scholar]
  205. Shimakawa, R., Pérez-Martínez, J. M., Dannerbauer, H., et al. 2024a, ApJ, 977, 73 [NASA ADS] [Google Scholar]
  206. Shimakawa, R., Pérez-Martínez, J. M., Koyama, Y., et al. 2024b, MNRAS, 528, 3679 [NASA ADS] [CrossRef] [Google Scholar]
  207. Smail, I., Geach, J. E., Swinbank, A. M., et al. 2014, ApJ, 782, 19 [NASA ADS] [CrossRef] [Google Scholar]
  208. Smith, R., Davies, J. I., & Nelson, A. H. 2010, MNRAS, 405, 1723 [NASA ADS] [Google Scholar]
  209. Sobral, D., Smail, I., Best, P. N., et al. 2013, MNRAS, 428, 1128 [NASA ADS] [CrossRef] [Google Scholar]
  210. Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 [Google Scholar]
  211. Sparre, M., & Springel, V. 2017, MNRAS, 470, 3946 [Google Scholar]
  212. Speagle, J. S., Steinhardt, C. L., Capak, P. L., & Silverman, J. D. 2014, ApJS, 214, 15 [Google Scholar]
  213. Spilker, J. S., Suess, K. A., Setton, D. J., et al. 2022, ApJ, 936, L11 [Google Scholar]
  214. Springel, V. 2010, MNRAS, 401, 791 [Google Scholar]
  215. Stevens, A. R. H., Lagos, C. D. P., Cortese, L., et al. 2021, MNRAS, 502, 3158 [NASA ADS] [CrossRef] [Google Scholar]
  216. Strazzullo, V., Daddi, E., Gobat, R., et al. 2016, ApJ, 833, L20 [NASA ADS] [CrossRef] [Google Scholar]
  217. Suh, H., Civano, F., Trakhtenbrot, B., et al. 2020, ApJ, 889, 32 [NASA ADS] [CrossRef] [Google Scholar]
  218. Taamoli, S., Mobasher, B., Chartab, N., et al. 2024, ApJ, 966, 18 [NASA ADS] [CrossRef] [Google Scholar]
  219. Tacchella, S., Dekel, A., Carollo, C. M., et al. 2016, MNRAS, 458, 242 [NASA ADS] [CrossRef] [Google Scholar]
  220. Tacconi, L. J., Genzel, R., Saintonge, A., et al. 2018, ApJ, 853, 179 [Google Scholar]
  221. Tadaki, K.-I., Kodama, T., Tamura, Y., et al. 2014, ApJ, 788, L23 [NASA ADS] [CrossRef] [Google Scholar]
  222. Tadaki, K.-I., Kodama, T., Hayashi, M., et al. 2019, PASJ, 71, 40 [NASA ADS] [CrossRef] [Google Scholar]
  223. Tanaka, M., De Breuck, C., Venemans, B., & Kurk, J. 2010, A&A, 518, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  224. Tanaka, T. S., Silverman, J. D., Ding, X., et al. 2024, ArXiv e-prints [arXiv:2401.13742] [Google Scholar]
  225. Tollet, É., Cattaneo, A., Macciò, A. V., Dutton, A. A., & Kang, X. 2019, MNRAS, 485, 2511 [NASA ADS] [CrossRef] [Google Scholar]
  226. Torrey, P., Vogelsberger, M., Marinacci, F., et al. 2019, MNRAS, 484, 5587 [NASA ADS] [Google Scholar]
  227. Toshikawa, J., Kashikawa, N., Overzier, R., et al. 2014, ApJ, 792, 15 [NASA ADS] [CrossRef] [Google Scholar]
  228. Toshikawa, J., Wuyts, S., Kashikawa, N., et al. 2024, MNRAS, 527, 6276 [Google Scholar]
  229. Tozzi, P., Pentericci, L., Gilli, R., et al. 2022a, A&A, 662, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  230. Tozzi, P., Gilli, R., Liu, A., et al. 2022b, A&A, 667, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  231. Tran, K.-V. H., Nanayakkara, T., Yuan, T., et al. 2015, ApJ, 811, 28 [NASA ADS] [CrossRef] [Google Scholar]
  232. Travascio, A., Cantalupo, S., Tozzi, P., et al. 2025, A&A, 694, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  233. Valentino, F., Daddi, E., Strazzullo, V., et al. 2015, ApJ, 801, 132 [NASA ADS] [CrossRef] [Google Scholar]
  234. van den Bosch, F. C., Jiang, F., Campbell, D., & Behroozi, P. 2016, MNRAS, 455, 158 [NASA ADS] [CrossRef] [Google Scholar]
  235. van Dokkum, P. G., Nelson, E. J., Franx, M., et al. 2015, ApJ, 813, 23 [NASA ADS] [CrossRef] [Google Scholar]
  236. Vito, F., Brandt, W. N., Lehmer, B. D., et al. 2020, A&A, 642, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  237. Vito, F., Brandt, W. N., Comastri, A., et al. 2024, A&A, 689, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  238. Wang, T., Elbaz, D., Daddi, E., et al. 2018, ApJ, 867, L29 [NASA ADS] [CrossRef] [Google Scholar]
  239. Watson, C., Tran, K.-V., Tomczak, A., et al. 2019, ApJ, 874, 63 [NASA ADS] [CrossRef] [Google Scholar]
  240. Weinberger, R., Springel, V., Hernquist, L., et al. 2017, MNRAS, 465, 3291 [Google Scholar]
  241. Whitaker, K. E., Williams, C. C., Mowla, L., et al. 2021, Nature, 597, 485 [NASA ADS] [CrossRef] [Google Scholar]
  242. Williams, C. C., Spilker, J. S., Whitaker, K. E., et al. 2021, ApJ, 908, 54 [NASA ADS] [CrossRef] [Google Scholar]
  243. Wisnioski, E., Förster Schreiber, N. M., Fossati, M., et al. 2019, ApJ, 886, 124 [Google Scholar]
  244. Wuyts, S., Förster Schreiber, N. M., Nelson, E. J., et al. 2013, ApJ, 779, 135 [Google Scholar]
  245. Zanella, A., Valentino, F., Gallazzi, A., et al. 2023, MNRAS, 524, 923 [CrossRef] [Google Scholar]
  246. Zavala, J. A., Casey, C. M., Scoville, N., et al. 2019, ApJ, 887, 183 [NASA ADS] [CrossRef] [Google Scholar]
  247. Zeballos, M., Aretxaga, I., Hughes, D. H., et al. 2018, MNRAS, 479, 4577 [Google Scholar]
  248. Zhang, Y., Zheng, X. Z., Shi, D. D., et al. 2022, MNRAS, 512, 4893 [NASA ADS] [CrossRef] [Google Scholar]
  249. Zhang, Y. H., Dannerbauer, H., Pérez-Martínez, J. M., et al. 2024, A&A, 692, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  250. Zinger, E., Pillepich, A., Nelson, D., et al. 2020, MNRAS, 499, 768 [CrossRef] [Google Scholar]
  251. Zirm, A. W., Stanford, S. A., Postman, M., et al. 2008, ApJ, 680, 224 [NASA ADS] [CrossRef] [Google Scholar]
  252. Zolotov, A., Dekel, A., Mandelker, N., et al. 2015, MNRAS, 450, 2327 [Google Scholar]

Appendix A: Summary Tables

Table A.1.

Physical properties of the Spiderweb protocluster HAES for which we have obtained ATCA CO(1-0) measurements.

Table A.2.

List of protocluster (PCL), protocluster core (PCL-core), cluster (CL), and field samples with CO(1-0) emitters from the literature at z > 1.5.

All Tables

Table 1.

Summary of CIGALE SED fitting modules and parameter configuration.

Table 2.

Median physical properties of the samples used in Sect. 4.

Table A.1.

Physical properties of the Spiderweb protocluster HAES for which we have obtained ATCA CO(1-0) measurements.

Table A.2.

List of protocluster (PCL), protocluster core (PCL-core), cluster (CL), and field samples with CO(1-0) emitters from the literature at z > 1.5.

All Figures

thumbnail Fig. 1.

Spiderweb protocluster field at z ≈ 2.16. Blue circles display the full sample of candidate HAEs from Koyama et al. (2013) and Shimakawa et al. (2018b). Red circles show those HAEs with measured spectroscopic redshift (Shimakawa et al. 2018b; Pérez-Martínez et al. 2023). Orange crosses and contours respectively show the CO(1−0) emitters reported by Jin et al. (2021) and the limits of the COALAS ATCA footprint. Empty squares depict the subsample of SMGs reported by Dannerbauer et al. (2014) in this field. Black stars show the LAEs from Pentericci et al. (2000), while X-ray sources from Tozzi et al. (2022a) are shown by black crosses. The Spiderweb galaxy is marked with a yellow star. We display a bar with 1 Mpc physical scale at z = 2.16 in the lower-left corner of the diagram for reference.

In the text
thumbnail Fig. 2.

Redshift distribution of sources. The complete sample of spectroscopic HAEs (Shimakawa et al. 2018b, 2024b; Pérez-Martínez et al. 2023) is shown in red. The sample of ATCA CO(1−0) emitters (Jin et al. 2021) is shown in yellow. The sample of HAEs within the ATCA footprint (see Fig. 1) that is used in this work is depicted by the blue contour histogram. The black solid line shows the shape and redshift coverage of the Subaru/MOIRCS narrow-band filter used to select the HAEs.

In the text
thumbnail Fig. 3.

Two examples of our upper limit determination method. The upper panels display the HAE ID 790 for which Jin et al. (2021) reported CO(1−0) at S/N > 4, while the lower panel displays an undetected CO(1−0) source (HAE ID 1181) for which we compute upper limits. The first column displays the moment zero of the ATCA datacube after collapsing it at the redshift of our source ±500 km/s. The black cross shows the exact position of the HAE under scrutiny. The red empty square shows the position of the spaxel containing the peak flux after inspecting a 3 × 3 grid around the HAEs spaxel. The beam size is shown in the lower left corner. The second column depicts the spectra extracted from the peak flux spaxel. The vertical red dashed line shows the systemic velocity (i.e., redshift) of the galaxy from Hα while the solid red lines display the limits (±500 km/s) for flux integration (blue area). Velocities are relative to z = 2.1612 as in Jin et al. (2021).

In the text
thumbnail Fig. 4.

Comparison of the LCO(1 − 0) obtained by applying our methodology and the source detection codes outlined in Jin et al. (2021) for HAEs with CO(1−0) detection at S/N > 4.

In the text
thumbnail Fig. 5.

Star-forming main sequence diagram. The black solid line and shaded grey region depict the Main Sequence of star formation at z = 2.16 from Speagle et al. (2014). Red small circles the spectroscopically confirmed HAEs of Shimakawa et al. (2018b) and Pérez-Martínez et al. (2023) within the ATCA datacube footprint. Blue circles show the location of the CO(1−0) emitters at S/N > 4 from Jin et al. (2021) which are part of our sample of HAEs. For context, we also add empty squares to show the overlap with the sample of SMGs in this field from Dannerbauer et al. (2014) and with the X-ray emitters reported by Tozzi et al. (2022a).

In the text
thumbnail Fig. 6.

Left: Total gas fraction versus stellar mass diagram. The solid green line and shaded region represent the field scaling relations for Main Sequence (MS) galaxies from Tacconi et al. (2018). We overplot our sample in blue color (both for measurements and upper limits) and compare it with other protocluster samples at the cosmic noon such as Wang et al. (2018) in orange, and coeval field samples (Kaasinen et al. 2019 in pink; Riechers et al. 2020 in grey; Frias Castillo et al. 2023 in green). In addition, we marked those galaxies with signs of AGN activity with red crosses (X-ray emission, Tozzi et al. 2022a) and red triangles (high [NII]/Hα, Pérez-Martínez et al. 2023) and those identified as SMGs by Dannerbauer et al. (2014) as black empty squares. The solid red line and red area depict the fit and uncertainty of the PKS1138 and CLJ1001 protocluster samples using a logistic function, as proposed by Popping et al. (2012). The grey shaded area displays the 1σ detection limit given by the median rms of the entire ATCA mosaic. Right: Total gas fraction versus stellar mass diagram after stacking. Our sample is divided into three bins: Low-mass galaxies (log M*/M < 10.5); massive galaxies (log M*/M > 10.5) excluding AGN candidates; and massive AGN candidates. In addition, median values and standard deviations for the comparison samples are displayed using the same symbol and color scheme as in the left-hand panel. The empty symbols depict the median values of the Kaasinen et al. (2019) and Frias Castillo et al. (2023) field samples after applying αCO ≈ 1, typical of high-z starbursts and SMGs (see Sects. 4.1 and 4.2). We apply a 0.02 dex cosmetic shift in M* to these two data points to improve the visibility of their error bars.

In the text
thumbnail Fig. 7.

Stacked spectra around the CO(1−0) emission line for the three bins defined in Sect. 4.2: Low-mass galaxies (log M*/M < 10.5), massive galaxies (log M*/M > 10.5) excluding AGN candidates, and massive AGN candidates. We achieve S/N ≈ 8, 10 and 13 for our stacked measurements respectively.

In the text
thumbnail Fig. 8.

Molecular gas mass versus star formation rate diagrams. Left: Individual sources from the samples displayed in Fig. 6. Right: Stacked (PKS1138) and median results of each sample displayed in the left-hand panel. Colors and symbols are the same as in Fig. 6. Dashed grey lines display depletion times between 0.1 and 10 Gyr.

In the text
thumbnail Fig. 9.

Left: Phase space diagram. Objects are color-coded by their molecular to stellar mass ratio μgas offsets with respect to the scaling relation of Tacconi et al. (2018). Solid grey lines divide the phase space into three regimes with galaxies being in the outskirts (η ≥ 2), infalling (0.4 ≤ η ≤ 2) or central regions (η ≤ 0.4) similarly to Pérez-Martínez et al. (2023). CO(1−0) measurements at S/N > 4 are depicted as solid circles while inverted triangles show upper limits. Right: Molecular to stellar mass ratio (μgas) offsets from the scaling relation as a function of the environmental parameter η. In addition, the pentagon and two hexagons depict the three stellar mass and AGN activity bins presented in Table 2. The remaining symbols and colors remain the same as in the left-hand diagram.

In the text
thumbnail Fig. 10.

Left: Star formation rate offsets from the Main Sequence of Tacconi et al. (2018) as a function of the environmental parameter η. Right: Star formation efficiency offsets with respect to the scaling relation of Tacconi et al. (2018) as a function of η. The symbols and colors remain the same as in Fig. 9.

In the text
thumbnail Fig. 11.

Gas fraction as a function of local density defined by the area enclosing two, five, and ten neighboring galaxies (Σ2, Σ5, Σ10. The symbols and colors remain the same as in Fig. 9.

In the text
thumbnail Fig. 12.

Accretion regime diagram. Colored stars display the dark matter halo mass of the Spiderweb protocluster derived using the SZ effect and dynamical estimates by Di Mascolo et al. (2023) and Shimakawa et al. (2014) respectively. Grey squares depict the position of our sample of HAES after applying the stellar-to-halo mass relation in Behroozi et al. (2019) for satellite galaxies (see Sect. 5.1). For comparison, we use green big circles to show a sample of protoclusters from Overzier (2016) while the small red circles represent a sample of SZ-detected z < 1.5 clusters from Bleem et al. (2015).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.