alphaFeedback

WorksPhysics and Astronomy

The HST-Hyperion Survey: X-ray AGN Enhancement in the Hyperion Proto-supercluster

AI contributions
Not recorded for this paper
View details

Abstract

Dense environments can influence both galaxy evolution and supermassive black hole growth, but whether they enhance or suppress (or both) active galactic nuclei (AGN) activity at high redshift (z≳2z\gtrsim2) remains uncertain. As one of the most massive structures identified in the early universe, the Hyperion proto-supercluster at z≈2.5z\approx2.5 provides a unique opportunity to investigate this connection during an epoch near the peak of both cosmic star formation rate and AGN activity density. To examine the impact of galaxy environment on AGN activity, we measure the fraction of X-ray-selected AGN in Hyperion as a whole, its highly overdense peaks, and its outskirts, comparing these measurements to a coeval field sample. We identify X-ray AGN using the deep Chandra COSMOS-Legacy Survey and combine optical/near-infrared (NIR) photometry and photometric redshift information from COSMOS2020 with extensive ground-based spectroscopy and HST-grism observations. Redshift information is incorporated statistically using a Monte Carlo framework. We find a higher X-ray AGN fraction in Hyperion (7.4−1.8+3.2%7.4^{+3.2}_{-1.8}\%) than in the coeval field (3.2−0.5+0.6%3.2^{+0.6}_{-0.5}\%), corresponding to a ∼2.3σ\sim2.3σ difference. Strikingly, this enhancement is concentrated in the overdense peaks, where the X-ray AGN fraction reaches 12.3−3.4+5.5%12.3^{+5.5}_{-3.4}\% (15.0−4.4+6.9%15.0^{+6.9}_{-4.4}\% in the most massive peaks versus 6.2−2.1+11.9%6.2^{+11.9}_{-2.1}\% in the less massive peaks), while the lower-density Hyperion outskirts (2.1−1.1+3.6%2.1^{+3.6}_{-1.1}\%) show no evidence of enhancement relative to the field. These results add to the limited but growing evidence for enhanced AGN activity in high-redshift protostructures, suggesting that the dense environments at these epochs may preferentially trigger AGN during the early stages of cluster assembly.

Introduction

Galaxy properties vary strongly with environment in the lower-redshift, local universe, from late-types typically found in low-density field regions to the prevalence of early-types in dense galaxy clusters (e.g., A. Dressler [17]). Gas accretion can regulate the supply of material available to galaxies (e.g., F. van de Voort et al. 2017), while environmental processes such as ram-pressure stripping can remove gas from galaxies (e.g., J. E. Gunn & J. R. Gott [24]; A. Boselli et al. [6]; K. Xu et al. 2026). At higher redshifts, galaxy interactions increasingly provide a mechanism for redistributing gas (e.g., E. A. Shah et al. [56]; M. Romano et al. [53]; E. A. Shah et al. [57]). Specifically, interactions and mergers can drive gas toward galaxy centers, potentially triggering active galactic nuclei (AGN; e.g., T. Storchi-Bergmann & A. Schnorr-Müller 2019). Once triggered, AGN can suppress star formation and contribute to galaxy quenching, working with environmental effects to enhance gas removal (e.g., A. Fabian [18]; B. R. McNamara & P. E. J. Nulsen [48]). AGN feedback may also augment environmental processes that remove or limit gas in dense environments (e.g., Y. M. Bahé & I. G. McCarthy [5]). Enhanced AGN activity in dense environments may therefore influence both host galaxies and the surrounding large-scale structure, motivating measurements of its environmental dependence.

Protostructures at z≳2z \gtrsim2 allow us to examine this connection during the early stages of cluster assembly (e.g., R. A. Overzier [50]; Y.-K. Chiang et al. [10]; S. Alberts & A. Noble [2]; E. A. Shah et al. [58]). At this epoch, protostructures are still assembling, with their galaxy populations showing enhanced star formation and stellar mass growth (e.g., P. Staab et al. [60]), a greater relative abundance of massive galaxies (e.g., B. Forrest et al. [19]; B. Forrest et al. [21]; D. Sikorski et al. [59]), and elevated rates of galaxy interactions (e.g., F. Giddings et al. [22]) compared to field galaxies at the same epoch. These conditions make protostructures important environments for investigating the connection between supermassive black hole (SMBH) growth and the assembly of large-scale structure.

X-ray observations trace actively accreting SMBHs, including AGN whose emission may be obscured at optical wavelengths (e.g., R. C. Hickox & D. M. Alexander [28]). While X-ray surveys can miss highly obscured AGN, sources with X-ray emission below the detection threshold, and some radiatively inefficient AGN that may be identified through radio emission (e.g., R. C. Hickox et al. [29]), the extensive Chandra coverage of the COSMOS field to moderate depths provides a well-characterized catalog of 4016 X-ray sources, enabling a consistent comparison of X-ray-selected AGN between protoclusters and the coeval field (S. Marchesi et al. [44]; F. Civano et al. [12]).

This relationship between AGN activity and varying environments has been studied in low-redshift (z≲0.30z \lesssim0.30) galaxy clusters, where luminous AGN were found to be less common in dense environments compared to the field (P. Popesso & A. Biviano [52]). At intermediate redshifts (z∼1z \sim1), X-ray studies suggest a more complex relationship, with AGN being found outside the densest cluster environments (D. D. Kocevski et al. [38]) with evidence for a variety of triggering mechanisms (N. Rumbaugh et al. [54]). At even higher redshifts (z∼1z \sim1–33), radio-selected AGN can be reliable signposts of galaxy overdensities, with surveys such as CARLA and COBRA identifying galaxy protoclusters using radio AGN as tracers (N. A. Hatch et al. [27]; D. Wylezalek et al. [70]; R. Paterno-Mahler et al. [51]). At z≥2z \geq2, the relationship between X-ray AGN activity and environment remains less certain. Studies of individual protostructures have found enhanced X-ray AGN fractions relative to the field, including Spiderweb at z=2.16z = 2.16 (P. Tozzi et al. [63]), HS 1700+64 at z=2.30z = 2.30 (J. A. Digby-North et al. [16]), SSA22 at z=3.09z = 3.09 (B. D. Lehmer et al. [41]), and MQN01 at z∼3.25z \sim3.25 (A. Travascio et al. [64]), although other studies find no signs of enhancement (M. Macuga et al. [43]). More recently, statistical studies across samples of high-redshift protostructures have also found evidence for enhanced AGN activity in dense environments (E. A. Shah et al. [58]; M. N. Isla Llave et al. [35]).

In this study, we focus on the Hyperion protosupercluster in the COSMOS field (N. Scoville et al. [55]) at z≈2.5z \approx2.5, corresponding to a lookback time of approximately 11.15 billion years. Although smaller components of the structure were identified in earlier studies (C. Diener et al. [15]; Y.-K. Chiang et al. [11]; C. Casey et al. [9]; K.-G. Lee et al. [40]; T. Wang et al. [67]), O. Cucciati et al. [13] connected these overdensities into a proto-supercluster, estimating a total mass of ∼4.8×1015 M⊙\sim4.8 \times10^{15}\,M_{\odot} and identifying seven overdense peaks. More recent mapping of Hyperion using an expanded spectroscopic and photometric data set has refined the extent and properties of the structure (D. Sikorski et al. [59]). As one of the most massive known protostructures, Hyperion provides an opportunity to study AGN activity across a wide range of local overdensities.

We use the multiwavelength photometric, spectroscopic, HST-grism, and X-ray observations available in the COSMOS field (B. Forrest et al. [20]; B. Forrest et al. [21]) to compare the X-ray AGN activity in Hyperion to a coeval field sample. We account for host-galaxy stellar mass using spectral energy distribution (SED) fitting and a common stellar mass cut, while propagating @page 2/10 starts_mid=1 ends_mid=1 redshift uncertainties through Monte Carlo realizations. To determine whether X-ray AGN activity is associated with the local environment during the early stages of cluster assembly, we compare the X-ray AGN fraction between the field, Hyperion as a whole, its highly overdense peaks, and its lower-density outskirts.

The remainder of this paper is organized as follows. In § 2 we describe the photometric, spectroscopic, and X-ray data used in this study. In § 3 we describe the identification of X-ray AGN, environmental classification, MC analysis, SED fitting, and AGN fraction measurements. We present our results in § 4 and discuss their details in § 5. Throughout this work, we assume a Λ\LambdaCDM cosmology with H0=70H_0 = 70 km s−1^{-1} Mpc−1^{-1}, ΩΛ=0.73\Omega_{\Lambda} = 0.73, and Ωm=0.27\Omega_m = 0.27.

Photometric and Spectroscopic Data

To derive the properties of our AGN host galaxies, we perform robust SED fitting with LePHARE using photometry from the COSMOS2020 Classic Catalog v7.0 [69], which provides far-ultraviolet to mid-infrared photometry for approximately 1.7 million sources across ∼2\sim2 deg2^2 of the COSMOS field. The catalog also provides photometric redshifts with a typical precision of ∼0.01(1+z)\sim0.01(1 + z) for i<22.5i < 22.5 and better than 0.025(1+z)0.025(1 + z) for i<25i < 25. To improve redshift accuracy, we include spectroscopic redshifts compiled by the C3VO collaboration [42, 19], supplemented by spectroscopy from LATIS [49], CLAMATO [31], MOSDEF [39], and additional surveys, as well as HST-grism redshifts from the HST-Hyperion survey [21]. Redshift quality flags are used to determine whether a spectroscopic, grism, or photometric redshift is assigned to a particular source in our MC analysis (§ 3.3). Redshift quality flag and reliability assessment descriptions are available in B. C. Lemaux et al. [42] and B. Forrest et al. [21] for the grism sources.

To define the galaxy population we will use to measure the AGN fraction across different environments, we construct a parent galaxy sample from sources with COSMOS2020 photometry following the spatial and IRAC magnitude selection of D. Sikorski et al. [59]. This selection provides a stellar mass-limited sample to log⁡(M⋆/M⊙)≥10\log(M_{\star}/M_{\odot}) \ge10 across the Hyperion region and surrounding field. Across the MC realizations, the resulting parent catalog contains an average of approximately 8,223 sources, including 6,992 (85.0%) with photometric redshifts, 1,085 (13.2%) with spectroscopic redshifts, and 146 (1.8%) with HST grism redshifts. Hyperion and coeval field galaxies are identified from this parent sample as described in § 3.

Methods

X-ray Matching and AGN Identification

X-ray AGN in this study were identified using the Chandra COSMOS-Legacy Survey catalog [44], which combines new and archival Chandra observations covering 2.2 deg2^2 of the COSMOS field. Effective exposures are ∼160\sim160 ks over the central 1.5 deg2^2 and ∼80\sim80 ks over the remaining area, with a limiting full-band (0.5 – 10 keV) flux of 8.9×10−168.9 \times10^{-16} erg cm−2^{-2} s−1^{-1} [12]. The X-ray luminosities used throughout this study are not corrected for internal absorption. The COSMOS-Legacy catalog provides absorption corrections based on catalog redshifts, which may differ from those adopted in our analysis to recalculate full-band X-ray luminosities (§ 3.3).

We identified X-ray AGN using the Chandra COSMOS-Legacy Survey optical/IR counterparts matched by S. Marchesi et al. [44]. These counterparts were cross-matched to our COSMOS2020-based master spectroscopic and photometric catalog. Of the 4016 X-ray sources in the Chandra COSMOS-Legacy catalog, 3582 (89.2%) had corresponding sources in our master catalog. Unmatched X-ray sources may have ambiguous or undetected optical/IR counterparts, particularly if highly obscured, or may have counterparts identified by S. Marchesi et al. [44] that are absent from our catalog due to its more restrictive source selection criteria. We visually inspected the optical/near-infrared (NIR) counterparts of X-ray sources frequently selected across the MC realizations to assess identification reliability.

After recalculating the full-band X-ray luminosities using the redshift assigned in each MC realization (§ 3.3), we found that the X-ray sources in the Hyperion and coeval field samples spanned log⁡(LX/erg s−1)=43.17−44.85\log(L_X/\mathrm{erg\,s}^{-1}) = 43.17 - 44.85 across the realizations. All sources, therefore, remained within the luminosity regime where X-ray emission is expected to be dominated by SMBH accretion rather than stellar processes [28], and their classification as X-ray AGN did not change between realizations.

Environment Classification

Because high-redshift protostructures have not yet virialized, their density distributions are inhomogeneous and contain highly overdense peaks surrounded by regions of more moderate density. We therefore characterized the galaxy environment using the underlying galaxy overdensity field. We determined the galaxy environments using the Voronoi Monte Carlo (VMC) mapping technique described in detail by B. C. Lemaux et al. [42] and D. Sikorski et al. [59]. Briefly, the VMC maps used available photometric and spectroscopic redshift information to statistically reconstruct the three-dimensional distribution of galaxies. Within overlapping redshift slices, a Voronoi tessellation divided the two-dimensional projected sky plane into cells surrounding each galaxy, with smaller cells corresponding to regions of higher galaxy density. For each redshift slice, we generated 100 MC realizations incorporating the individual galaxy redshift uncertainties, and the resulting maps were combined across redshift slices to construct a three-dimensional density cube. These realizations constructed the published VMC density maps and are distinct from those used to propagate redshift uncertainties in our analysis (§ 3.3). Rather than using galaxy density directly, we characterized environment using the galaxy overdensity relative to the background density at each redshift, expressed as a significance, nσn_{\sigma} ([13, 32, 59]).

Following D. Sikorski et al. [59], we identified protostructures as contiguous regions with nσ≥2.5n_{\sigma} \ge2.5, containing at least one overdense peak (nσ≥4n_{\sigma} \ge4) and an enclosed mass of Mtot≥1013 M⊙M_{\mathrm{tot}} \ge10^{13}\,\mathrm{M}_{\odot}. Within each protostructure, peaks were defined as contiguous sets of voxels with nσ≥4n_{\sigma} \ge4, while the remaining connected voxels (2.5≤nσ<42.5 \le n_{\sigma} < 4) were classified as outskirts.

Galaxies were assigned to Hyperion based on their redshift, sky position, and the corresponding environment classification from the VMC map. The Hyperion sample was further divided into galaxies residing within the overdense peaks and those in the protostructure outskirts. We selected a coeval field sample from galaxies within the same right ascension and declination boundaries as Hyperion but lying in the foreground (2.15≤z≤2.252.15 \le z \le2.25) and background (2.80≤z≤2.902.80 \le z \le2.90) redshift intervals. Restricting both samples to the same angular footprint helps ensure similar spectroscopic sampling and spectroscopic redshift fractions, while the field redshift intervals bracket Hyperion to minimize differences in cosmic time and reduce the likelihood that galaxies are scattered between the Hyperion and field samples by redshift uncertainties (D. Sikorski et al. [59]). Field membership was determined by whether a galaxy was associated with an identified protostructure, rather than by its local nσn_{\sigma} value. Thus, locally overdense galaxies (nσ≥2.5n_{\sigma} \ge2.5) could remain in the field sample if they were not members of an identified structure. However, only (1.8−0.2+0.2%1.8^{+0.2}_{-0.2}\%) of field galaxies were found at nσ≥2.5n_{\sigma} \ge2.5 across the MC realizations. The resulting Hyperion and field samples are shown as a function of redshift and local overdensity in Figure 1, with X-ray sources indicated by larger and darker points.

Distribution of galaxies in redshift and overdensity space

Figure 1. Distribution of galaxies in redshift and overdensity (nσn_{\sigma}) space, using the best available catalog redshift for each galaxy. Galaxies assigned to Hyperion are marked in pink, with corresponding X-ray-identified AGN marked in red. Field galaxies are marked in light blue with corresponding X-ray-identified AGN marked in blue. Unused galaxies are marked in gray (including galaxies outside of the overdensity or redshift range for Hyperion and the field sample).

Monte Carlo Analysis and AGN Fraction Measurements

To propagate redshift uncertainties throughout the analysis, the AGN fraction measurements were performed across 100 MC realizations of the galaxy catalog. In each realization, a redshift was assigned to each galaxy using the available spectroscopic, HSTHST-grism, and/or COSMOS2020 photometric redshift information. Spectroscopic and grism redshift measurements were assigned according to the reliabilities associated with their quality flags, as described in § 2. When neither measurement was selected, or when no spectroscopic or grism redshift was available, a redshift was randomly sampled from the COSMOS2020 photometric redshift probability distribution (J. R. Weaver et al. [69]).

For each MC realization, we fit the 35-band COSMOS2020 photometry at each galaxy’s assigned redshift using LePHARE ([3, 34]) and galaxy templates based on the stellar population models of G. Bruzual & S. Charlot [7], following D. Sikorski et al. [59]. This generated a unique catalog of derived galaxy properties for each realization.

Because the redshift assigned to an X-ray source could differ from its catalog redshift, we recalculated the full-band X-ray luminosity in each realization using the corresponding luminosity-distance ratio and a kk-correction assuming Γ=1.4\Gamma= 1.4, consistent with S. Marchesi et al. [44]. Across the MC realizations, an average of 52.6% of X-ray sources had ∣Δz∣>0.1|\Delta z| > 0.1 relative to the redshift adopted in the Chandra COSMOS-Legacy catalog, with a median ∣Δz∣=0.40|\Delta z| = 0.40 for these sources.

Because absorption corrections from S. Marchesi et al. [44] were modeled using the Chandra COSMOS-Legacy catalog redshifts, we did not apply them to our recalculated luminosities and instead left luminosities as absorption-uncorrected values. We found that 62.2% of Hyperion and 54.9% of field X-ray AGN had catalog-reported absorption correction factors or upper bounds that differed from 1. Among these sources, the mean absorption corrections were comparable between Hyperion and the field, corresponding to luminosity increases of approximately 25% in both environments. Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) tests returned no statistically significant differences in the absorption-correction distributions between the samples across the MC realizations. We also compared their hardness-ratio distributions as a model-independent test of obscuration, adopting limits for 15.4% and 28.1% of Hyperion and field sources, respectively. KS and AD tests again found no statistically significant differences. Although these comparisons cannot rule out differences in highly obscured AGN, the broadly consistent obscuration properties among X-ray-detected sources suggest that our relative comparisons of absorption-uncorrected luminosities would likely remain similar if absorption corrections were applied (i.e., neither sample appears more obscured, at least among the X-ray-detected AGN).

Because AGN activity depends strongly on host-galaxy stellar mass (G. Kauffmann et al. [37]; J. Aird et al. [1]), both the Hyperion and coeval field samples were restricted to galaxies with log⁡(M⋆/M⊙)≥10\log(M_\star/M_\odot) \ge10. We applied this stellar mass cut in each MC realization using the stellar mass derived from that realization’s SED fit. Applying the same stellar mass limit to both samples reduced the possibility that differences in AGN fraction were driven by differences in the underlying stellar mass distributions rather than environment. As shown in § 4, we found no statistically significant difference between the stellar mass distributions of the Hyperion and coeval field samples.

To measure a local nσn_\sigma value, we mapped the assigned MC redshift and sky position of each galaxy into the VMC density cube. Galaxies were assigned to the Hyperion, coeval field, peak, or outskirt samples using the environmental definitions described in § 3.2. After applying the stellar mass selection, we calculated the X-ray AGN fraction independently in each MC realization as the ratio of X-ray AGN to the total number of galaxies in the sample following E. A. Shah et al. [58].

We adopted the median AGN fraction across 100 MC realizations, with MC uncertainties determined from the 16th and 84th percentiles of the AGN fraction distribution. We then combined these uncertainties in quadrature with the uncertainty from the finite number of AGN, calculated assuming binomial statistics following E. Cameron [8]. The AGN enhancement was calculated as the ratio of the median AGN fraction measured in each Hyperion environment to that of the coeval field. We propagated the corresponding upper and lower uncertainties on the Hyperion and field AGN fractions through this ratio to determine the uncertainty on the enhancement.

Results

The X-ray AGN fractions measured for each environment are presented in Figure 2 and summarized in Table 1. The coeval field has an X-ray AGN fraction of 3.2−0.5%+0.6%3.2^{+0.6\%}_{-0.5\%}. In contrast, the full Hyperion sample has an AGN fraction of 7.4−1.8%+3.2%7.4^{+3.2\%}_{-1.8\%}, a factor of 2.3−0.7+1.12.3^{+1.1}_{-0.7} higher than the field.

EnvironmentNgalN_{\mathrm{gal}}X-ray AGNAGN fractionEnhancement
Field1510−51+491510^{+49}_{-51}49−9+949^{+9}_{-9}3.2−0.5+0.6%3.2^{+0.6}_{-0.5}\%—
Hyperion (nσ≥2.5n_{\sigma} \ge2.5)121−12+13121^{+13}_{-12}9−3+49^{+4}_{-3}7.4−1.8+3.2%7.4^{+3.2}_{-1.8}\%2.3−0.7+1.12.3^{+1.1}_{-0.7}
Outskirts (2.5≤nσ<4.02.5 \le n_{\sigma} < 4.0)63−9+1063^{+10}_{-9}1−1+31^{+3}_{-1}2.1−1.1+3.6%2.1^{+3.6}_{-1.1}\%0.6−0.4+1.10.6^{+1.1}_{-0.4}
Peaks (nσ≥4.0n_{\sigma} \ge4.0)58−8+958^{+9}_{-8}7−3+47^{+4}_{-3}12.3−3.4+5.5%12.3^{+5.5}_{-3.4}\%3.8−1.3+1.83.8^{+1.8}_{-1.3}

Table 1. X-ray AGN fractions across the Hyperion environments and coeval field. Galaxy and X-ray AGN counts represent the median values across the Monte Carlo realizations, with uncertainties combining the 16th and 84th percentile MC variations and Poisson counting uncertainties in quadrature. AGN fraction uncertainties combine MC variations and binomial uncertainties following E. Cameron [8] in quadrature.

X-ray AGN fraction as a function of overdensity

Figure 2. X-ray AGN fraction (in %) as a function of local overdensity in the Hyperion proto-supercluster and coeval field. The galaxy population is divided into three environments: the coeval field (blue), the Hyperion outskirts (2.5≤nσ<42.5 \le n_{\sigma} < 4; green), and the dense Hyperion peaks (nσ≥4n_{\sigma} \ge4; red). Error bars represent the total uncertainties on the AGN fractions, combining the Monte Carlo realizations and Bayesian binomial confidence intervals. The highest X-ray AGN fraction is observed in the dense Hyperion peaks, while the outskirts are consistent with the coeval field within the uncertainties.

Notably, the enhanced AGN activity is not distributed uniformly throughout Hyperion. Galaxies in the outskirts have an AGN fraction of 2.1−1.1%+3.6%2.1^{+3.6\%}_{-1.1\%}, consistent with the coeval field. In contrast, galaxies located in the dense Hyperion peaks have an AGN fraction of 12.3−3.4%+5.5%12.3^{+5.5\%}_{-3.4\%}, a factor of 3.8−1.3+1.83.8^{+1.8}_{-1.3} higher than the field and a 2.6σ2.6\sigma difference. These results suggest that the enhanced AGN activity is concentrated in the highest-density regions of the protostructure rather than uniformly throughout Hyperion.

To investigate which Hyperion peaks contribute to the enhanced AGN activity and whether the AGN fraction varies with peak mass, we compare the individual peaks hosting X-ray AGN and the combined more massive and less massive peak populations in Figure 3, using the peak assignments and masses of D. Sikorski et al. [59]. In the individual peak comparison, we find that Peak 2 has the highest AGN fraction at 33.3−11.8%+17.9%33.3^{+17.9\%}_{-11.8\%}, while Peaks 1, 3, and 6 have AGN fractions of 10.5−3.6%+11.1%10.5^{+11.1\%}_{-3.6\%}, 7.7−2.6%+14.0%7.7^{+14.0\%}_{-2.6\%}, and 14.3−6.1%+21.5%14.3^{+21.5\%}_{-6.1\%}, respectively. When grouped by peak mass, the three most massive peaks (Peaks 1–3) have a combined AGN fraction of 15.0−4.4%+6.9%15.0^{+6.9\%}_{-4.4\%}, compared with 6.2−2.1%+11.9%6.2^{+11.9\%}_{-2.1\%} for the less massive peaks (Peaks 4–7). Although the measured AGN fraction is higher in the more massive peaks, the large uncertainties, particularly among the less massive peaks, limit our ability to determine whether the X-ray AGN fraction depends on peak mass.

X-ray AGN fraction as a function of peak mass for the Hyperion density peaks

Figure 3. X-ray AGN fraction as a function of peak mass for the Hyperion density peaks. Individual points are shown for the four peaks that host X-ray AGN (Peaks 1, 2, 3, and 6). We find X-ray AGN fractions of 10.5−3.6+11.1%10.5^{+11.1}_{-3.6}\% (2.0/18.7) in Peak 1, 33.3−11.8+17.9%33.3^{+17.9}_{-11.8}\% (3.0/8.9) in Peak 2, 7.7−2.6+14.0%7.7^{+14.0}_{-2.6}\% (1.1/13.4) in Peak 3, and 14.3−6.1+21.5%14.3^{+21.5}_{-6.1}\% (1.0/7.2) in Peak 6. The values in parentheses give the mean NAGN/NgalN_{\mathrm{AGN}}/N_{\mathrm{gal}} across the Monte Carlo realizations. Peak masses are taken from D. Sikorski et al. [59]. The size of each point scales with the mean number of total galaxies identified in that peak or peak-mass group across the Monte Carlo realizations. The more massive peaks are Peaks 1–3 with Mpeak=0.88M_{\mathrm{peak}} = 0.88–2.02×1014 M⊙2.02 \times10^{14}\ M_{\odot}, while the less massive peaks are Peaks 4–7 with Mpeak=0.13M_{\mathrm{peak}} = 0.13–0.48×1014 M⊙0.48 \times10^{14}\ M_{\odot}. For the two mass groups, the AGN fraction is calculated from the total number of X-ray AGN and galaxies across the included peaks in each Monte Carlo realization, and the point is placed at the mean peak mass of the group. The uncertainty in the mass of each peak is ∼0.25\sim0.25 dex (P. Staab et al. [60]). The dashed line and shaded region show the coeval field AGN fraction and its uncertainty.

Because AGN incidence depends strongly on host-galaxy stellar mass, we examined whether differences in the stellar mass distributions of Hyperion and the coeval field could contribute to the observed environmental trends, despite the imposition of a stellar mass cut. Figure 4 compares the stellar mass distributions of the mass-selected galaxy populations and X-ray AGN hosts, as well as the full-band absorption-uncorrected X-ray luminosity distributions of the AGN. Two-sample KS tests were performed independently for each of the 100 MC realizations. We find median values of D=0.094D = 0.094 and p=0.236p = 0.236 for the full galaxy samples, D=0.222D = 0.222 and p=0.784p = 0.784 for the X-ray AGN host masses, and D=0.319D = 0.319 and p=0.364p = 0.364 for the full-band X-ray luminosities. AD tests on the stellar-mass distributions show consistent results, with median A2=0.907A^2 = 0.907 and p=0.136p = 0.136 for the stellar-mass distribution of the full galaxy samples and A2=−0.683A^2 = -0.683 and p=0.750p = 0.750 for the stellar-mass distribution of the X-ray AGN hosts. Neither test provides evidence that the stellar mass distributions differ significantly between Hyperion and the field, although the small number of X-ray AGN limits the statistical power of these comparisons.

Cumulative distributions of stellar mass, X-ray AGN host stellar mass, and full-band X-ray luminosity

Figure 4. Cumulative distributions of stellar mass (left), X-ray AGN host stellar mass (middle), and full-band (0.5 – 10 keV) X-ray luminosity (right) for the Hyperion and coeval field samples. Two-sample KS tests across the 100 Monte Carlo realizations give median pp-values of 0.236, 0.784, and 0.364, respectively. We find no statistically significant differences between the Hyperion and coeval field distributions in any of the three comparisons.

Discussion and Conclusions

We find an enhanced AGN fraction within Hyperion relative to the coeval field at z≈2.5z \approx2.5. This enhancement is concentrated in the densest regions, with the highly overdense peaks exhibiting an enhanced AGN fraction relative to the field while the lower-density outskirts are consistent with the field. This suggests that local environmental conditions within Hyperion may promote SMBH accretion during the early stages of cluster formation.

As discussed in § 1, studies of low-redshift galaxy clusters find AGN to be less common in dense environments than in the field (e.g., P. Popesso & A. Biviano [52]), while studies at z∼1z \sim1 find AGN outside the densest cluster environments (e.g., D. D. Kocesvski et al. 2009). At higher redshifts (z>2z > 2), studies of AGN activity in protostructures remain limited, but enhanced AGN fractions have been found in individual systems (B. D. Lehmer et al. [41]; J. A. Digby-North et al. [16]; P. Tozzi et al. [63]; A. Travascio et al. [64]) and statistical samples (E. A. Shah et al. [58]; M. N. Isla Llave et al. [35]). Specifically, the 2.3−0.7+1.12.3^{+1.1}_{-0.7} times higher X-ray AGN fraction in Hyperion is consistent with the 2.4−0.7+1.12.4^{+1.1}_{-0.7} enhancement found by M. N. Isla Llave et al. [35] for X-ray AGN in Hyperion among submillimeter-selected galaxies. Despite the different galaxy selections, both studies also find no significant differences in the stellar-mass distributions of galaxies or AGN hosts between protostructures and the field, suggesting that stellar mass alone cannot explain the enhancement.

Although the peaks represent the most overdense regions of Hyperion, they are considerably less dense than evolved galaxy clusters at z∼1z \sim1, and especially those at z≈0z \approx0. These regions are more comparable to grouplike environments, where galaxy mergers and tidal interactions may be more effective at triggering star formation and AGN activity ([38, 54]). Studies of Hyperion support this enhanced galaxy growth. Sikorski et al. in prep. find higher star formation rates (SFRs) for Hype-

rion members relative to the coeval field at fixed stellar mass, suggesting that the elevated star formation cannot be explained solely by the presence of more massive galaxies within the proto-supercluster. This star formation is specifically elevated in the Hyperion peaks, which show an order of magnitude increase in the SFR density relative to lower-density environments. Galaxy interactions may contribute to this enhanced activity. In the z∼0.9z \sim0.9 Cl1604 supercluster, a majority of X-ray AGN hosts show evidence of recent or ongoing interactions, consistent with observations where interactions contribute to triggering AGN activity at z∼1z \sim1 [38, 45]. Interactions may also be important in Hyperion, where galaxies have close kinematic companions ≳2.5\gtrsim2.5 times more frequently than coeval field galaxies [22]. The elevated interaction rate in Hyperion may contribute to its enhanced AGN activity, potentially with a delay between interactions and SMBH fueling that allows accretion to continue after merger signatures fade [30].

However, this rapid growth may begin to change in the highest-density regions. G. Gururajan et al. [25] find increasing SFRs but decreasing molecular gas fractions and depletion timescales toward higher densities at the ∼2σ\sim2\sigma level, suggesting rapid galaxy evolution in the densest regions. At the same time, SFRs are highest in the peaks across much of the stellar-mass range but begin to decline above log⁡(M⋆/M⊙)∼10.5\log(M_\star/M_\odot) \sim10.5 (Sikorski et al. in prep.). The X-ray AGN we identify in Hyperion occupy this higher-mass regime (10.25≲log⁡(M⋆/M⊙)≲11.110.25 \lesssim\log(M_\star/M_\odot) \lesssim11.1) and are preferentially found within the high-density peaks. Enhanced SMBH growth may therefore be occurring alongside declining star formation in the most massive peak galaxies, perhaps suggesting some causal connection. Massive quiescent galaxies have also been identified within dense environments in other high-redshift protostructures [36, 46 b].

Several processes could contribute to this transition during the evolution of Hyperion. The shorter gas depletion timescales observed in the denser environments [25] suggest that rapid gas consumption may limit the fuel available for star formation, and changes in the surrounding environment could decrease the supply of cold gas. Detections of hot intracluster gas in forming clusters at z>2z > 2 also show that a hot gaseous medium can develop early in cluster assembly through a combination of gravitational heating and AGN feedback [23, 67, 68, 14, 72]. AGN feedback could also expel gas, although our observations do not establish that AGN activity causes the suppressed star formation observed among the most massive peak galaxies. Hyperion may therefore be observed during a transitional period where rapid star formation and SMBH growth are still occurring while some of the more massive galaxies begin to show quenching and suppressed star formation.

Future work is needed to connect these processes. Analysis of close kinematic companions and the morphologies of the X-ray AGN hosts will help test whether AGN activity is associated with galaxy interactions. Additionally, extending this analysis to more systems, such as those identified by D. Hung et al. [32], and incorporating broader AGN selection methods, will enable more robust measurements of environmental trends and place stronger constraints on AGN activity across different environments at high redshift.

Acknowledgments

AC acknowledges support from the National Science Foundation Research Experiences for Undergraduates Program at the University of Hawai‘i Institute for Astronomy under NSF grant No. 2548017. AC is grateful to the members of the C3VO collaboration for their guidance, discussions, and support throughout this work.

FV and MIL acknowledge support from 2023 INAF Ricerca Fondamentale Large GO grant.

This work makes use of photometry from the COS-MOS2020 catalog, and we thank the COSMOS collaboration for their efforts in obtaining, reducing, and publicly releasing this data. This work also makes use of spectroscopic data compiled by the C3VO collaboration, including observations from the HST-Hyperion survey. We additionally acknowledge the work of the Chandra COSMOS-Legacy Survey collaboration, which provides the X-ray source catalog used in this analysis.

This work 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. These observations are associated with programs 12328 and 16684. This work is also based in part on observations made with the Chandra X-ray Observatory. This work is based in part on observations collected at the European Southern Observatory under ESO programmes 175.A-0839, 179.A-2005, and 185. A0791, as well as work supported by the National Science Foundation under Grant No. 1908422. Supported by the international Gemini Observatory, a program of NSF NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the U.S. National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. This work is based in part on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/IRFU, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This research is based in part on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. This work is also based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA.

Some of the spectrographic data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation.

The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the Native Hawaiian community. We are most fortunate to have the opportunity to conduct research that benefits from observations obtained from this mountain.

Facilities: CXO, CFHT, GALEX, HST(WFC3), Keck(MOSFIRE), Keck(DEIMOS), Spitzer(IRAC), Subaru(HSC, Suprime-Cam), VISTA(VIRCAM), VLT/ESO

Software: Astropy ([4]), Matplotlib ([33]), NumPy ([26]), SciPy ([66]), pandas ((The pandas development team 2020))

References

  1. [1]Aird, J., Coil, A. L., & Georgakakis, A. 2019, Monthly Notices of the Royal Astronomical Society, 484, 4360–4378, doi: 10.1093/mnras/stz125
  2. [2]Alberts, S., & Noble, A. 2022, Universe, 8, 554, doi: 10.3390/universe8110554
  3. [3]Arnouts, S., Cristiani, S., Moscardini, L., et al. 1999, Monthly Notices of the Royal Astronomical Society, 310, 540
  4. [4]Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, The Astrophysical Journal, 935, 167
  5. [5]Bahé, Y. M., & McCarthy, I. G. 2015, MNRAS, 447, 969, doi: 10.1093/mnras/stu2293
  6. [6]Boselli, A., Fossati, M., & Sun, M. 2022, The Astronomy and Astrophysics Review, 30, doi: 10.1007/s00159-022-00140-3
  7. [7]Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000, doi: 10.1046/j.1365-8711.2003.06897.x
  8. [8]Cameron, E. 2011, PASA, 28, 128, doi: 10.1071/AS10046
  9. [9]Casey, C., Cooray, A., Capak, P., et al. 2015, The Astrophysical Journal Letters, 808, L33
  10. [10]Chiang, Y.-K., Overzier, R. A., Gebhardt, K., & Henriques, B. 2017, ApJL, 844, L23, doi: 10.3847/2041-8213/aa7e7b
  11. [11]Chiang, Y.-K., Overzier, R. A., Gebhardt, K., et al. 2015, The Astrophysical Journal, 808, 37
  12. [12]Civano, F., Marchesi, S., Comastri, A., et al. 2016, The Astrophysical Journal, 819, 62
  13. [13]Cucciati, O., Lemaux, B. C., Zamorani, G., et al. 2018, A&A, 619, A49, doi: 10.1051/0004-6361/201833655
  14. [14]Di Mascolo, L., Saro, A., Mroczkowski, T., et al. 2023, Nature, 615, 809, doi: 10.1038/s41586-023-05761-x
  15. [15]Diener, C., Lilly, S., Ledoux, C., et al. 2015, The Astrophysical Journal, 802, 31
  16. [16]Digby-North, J. A., Nandra, K., Laird, E. S., et al. 2010, MNRAS, 407, 846, doi: 10.1111/j.1365-2966.2010.16977.x
  17. [17]Dressler, A. 1980, ApJ, 236, 351, doi: 10.1086/157753
  18. [18]Fabian, A. 1999, Monthly Notices of the Royal Astronomical Society, 308, L39
  19. [19]Forrest, B., Lemaux, B. C., Shah, E. A., et al. 2024a, The Astrophysical Journal, 971, 169
  20. [20]Forrest, B., Cooper, M. C., Muzzin, A., et al. 2024b, The Astrophysical Journal, 977, 51, doi: 10.3847/1538-4357/ad8b1c
  21. [21]Forrest, B., Shen, L., Lemaux, B. C., et al. 2025, ApJ, 985, 61, doi: 10.3847/1538-4357/adc252
  22. [22]Giddings, F., Lemaux, B. C., Forrest, B., et al. 2026, A&A, 706, A55, doi: 10.1051/0004-6361/202554407
  23. [23]Gobat, R., Daddi, E., Onodera, M., et al. 2011, A&A, 526, A133, doi: 10.1051/0004-6361/201016084
  24. [24]Gunn, J. E., & Gott, III, J. R. 1972, ApJ, 176, 1, doi: 10.1086/151605
  25. [25]Gururajan, G., Cucciati, O., Lemaux, B., et al. 2025, Astronomy & Astrophysics, 698, A312
  26. [26]Harris, C. R., Millman, K. J., Van Der Walt, S. J., et al. 2020, Nature, 585, 357
  27. [27]Hatch, N. A., Röttgering, H. J. A., Miley, G. K., et al. 2013, Monthly Notices of the Royal Astronomical Society, 436, 2244–2253, doi: 10.1093/mnras/stt1734
  28. [28]Hickox, R. C., & Alexander, D. M. 2018, Annual Review of Astronomy and Astrophysics, 56, 625
  29. [29]Hickox, R. C., Jones, C., Forman, W. R., et al. 2009, The Astrophysical Journal, 696, 891–919, doi: 10.1088/0004-637X/696/1/891
  30. [30]Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, The Astrophysical Journal Supplement Series, 175, 356
  31. [31]Horowitz, B., Lee, K.-G., Ata, M., et al. 2022, The Astrophysical Journal Supplement Series, 263, 27, doi: 10.3847/1538-4365/ac982dDOI
  32. [32]Hung, D., Lemaux, B. C., Cucciati, O., et al. 2025, ApJ, 980, 155, doi: 10.3847/1538-4357/ada616DOI
  33. [33]Hunter, J. D. 2007, Computing in science & engineering, 9, 90
  34. [34]Ilbert, O., Arnouts, S., Mccracken, H. J., et al. 2006, Astronomy & Astrophysics, 457, 841
  35. [35]Isla Llave, M. N., Vito, F., Traina, A., et al. 2026, A&A, 714, A12, doi: 10.1051/0004-6361/202659613
  36. [36]Ito, K., Tanaka, M., Valentino, F., et al. 2023, ApJL, 945, L9, doi: 10.3847/2041-8213/acb49bDOI
  37. [37]Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055, doi: 10.1111/j.1365-2966.2003.07154.xDOI
  38. [38]Kocevski, D. D., Lubin, L. M., Lemaux, B. C., et al. 2009, ApJ, 700, 901, doi: 10.1088/0004-637X/700/2/901DOI
  39. [39]Kriek, M., Shapley, A. E., Reddy, N. A., et al. 2015, ApJS, 218, 15, doi: 10.1088/0067-0049/218/2/15DOI
  40. [40]Lee, K.-G., Hennawi, J. F., White, M., et al. 2016, The Astrophysical Journal, 817, 160
  41. [41]Lehmer, B. D., Alexander, D. M., Chapman, S. C., et al. 2009, MNRAS, 400, 299, doi: 10.1111/j.1365-2966.2009.15449.xDOI
  42. [42]Lemaux, B. C., Cucciati, O., Le Fèvre, O., et al. 2022, A&A, 662, A33, doi: 10.1051/0004-6361/202039346DOI
  43. [43]Macuga, M., Martini, P., Miller, E. D., et al. 2019, ApJ, 874, 54, doi: 10.3847/1538-4357/ab0746DOI
  44. [44]Marchesi, S., Civano, F., Elvis, M., et al. 2016, ApJ, 817, 34, doi: 10.3847/0004-637X/817/1/34DOI
  45. [45]Martini, P., Miller, E. D., Brodwin, M., et al. 2013, ApJ, 768, 1, doi: 10.1088/0004-637X/768/1/1DOI
  46. [46]McConachie, I., Wilson, G., Forrest, B., et al. 2025a, ApJ, 978, 17, doi: 10.3847/1538-4357/ad8f36DOI
  47. [48]McNamara, B. R., & Nulsen, P. E. J. 2007, ARA&A, 45, 117, doi: 10.1146/annurev.astro.45.051806.110625DOI
  48. [49]Newman, A. B., Rudie, G. C., Blanc, G. A., et al. 2020, The Astrophysical Journal, 891, 147
  49. [50]Overzier, R. A. 2016, A&A Rv, 24, 14, doi: 10.1007/s00159-016-0100-3DOI
  50. [51]Paterno-Mahler, R., Blanton, E. L., Brodwin, M., et al. 2017, ApJ, 844, 78, doi: 10.3847/1538-4357/aa7b89DOI
  51. [52]Popesso, P., & Biviano, A. 2006, A&A, 460, L23, doi: 10.1051/0004-6361:20066269DOI
  52. [53]Romano, M., Cassata, P., Morselli, L., et al. 2021, A&A, 653, A111, doi: 10.1051/0004-6361/202141306DOI
  53. [54]Rumbaugh, N., Lemaux, B. C., Tomczak, A., et al. 2017, MNRAS, 466, 496, doi: 10.1093/mnras/stw3091DOI
  54. [55]Scoville, N., Aussel, H., Brusa, M., et al. 2007, The Astrophysical Journal Supplement Series, 172, 1
  55. [56]Shah, E. A., Kartaltepe, J. S., Magnagoli, C. T., et al. 2020, ApJ, 904, 107, doi: 10.3847/1538-4357/abbf59DOI
  56. [57]Shah, E. A., Kartaltepe, J. S., Magnagoli, C. T., et al. 2022, ApJ, 940, 4, doi: 10.3847/1538-4357/ac96ebDOI
  57. [58]Shah, E. A., Lemaux, B. C., Forrest, B., et al. 2025, A&A, 704, A101, doi: 10.1051/0004-6361/202452055DOI
  58. [59]Sikorski, D., Forrest, B., Lemaux, B. C., et al. 2026, A&A, 708, A114, doi: 10.1051/0004-6361/202557084DOI
  59. [60]Staab, P., Lemaux, B. C., Forrest, B., et al. 2024, MNRAS, 528, 6934, doi: 10.1093/mnras/stae301DOI
  60. [63]Tozzi, P., Gilli, R., Liu, A., et al. 2022, Astronomy & Astrophysics, 667, A134, doi: 10.1051/0004-6361/202244337DOI
  61. [64]Travascio, A., Cantalupo, S., Tozzi, P., et al. 2025, A&A, 694, A165, doi: 10.1051/0004-6361/202452179DOI
  62. [66]Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261
  63. [67]Wang, T., Elbaz, D., Daddi, E., et al. 2016, The Astrophysical Journal, 828, 56
  64. [68]Wang, T., Elbaz, D., Daddi, E., et al. 2018, The Astrophysical Journal Letters, 867, L29, doi: 10.3847/2041-8213/aaeb2cDOI
  65. [69]Weaver, J. R., Kauffmann, O. B., Ilbert, O., et al. 2022, ApJS, 258, 11, doi: 10.3847/1538-4365/ac3078DOI
  66. [70]Wylezalek, D., Galametz, A., Stern, D., et al. 2013, ApJ, 769, 79, doi: 10.1088/0004-637X/769/1/79DOI
  67. [72]Zhou, D., Chapman, S. C., Aravena, M., et al. 2026, Nature, 649, 1130, doi: 10.1038/s41586-025-09901-3DOI

Paper details

vlnc:W_jhxkcpxw
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 8

downloads · 1

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

Not recorded for this paper · View statement →

Looking for related papers…
Contents