COSMOS Field Research Articles

Probing magnetic fields in the circumgalactic medium using polarization data from MIGHTEE

Context.The properties and evolution of magnetic fields surrounding galaxies are observationally largely unconstrained. The detection and study of these magnetic fields is important to understand galaxy evolution since magnetic fields are tracers for dynamical processes in the circumgalactic medium (CGM) and can have a significant impact on the evolution of the CGM.Aims.The Faraday rotation measure (RM) of the polarized light of background radio sources passing through the magnetized CGM of intervening galaxies can be used as a tracer for the strength and extent of magnetic fields around galaxies.Methods.We used rotation measures observed by the MIGHTEE-POL (MeerKAT International GHz Tiered Extragalactic Exploration polarization) survey by MeerKAT in the XMM-LSS and COSMOS fields to investigate the RM around foreground star-forming galaxies. We used spectroscopic catalogs of star-forming and blue cloud galaxies to measure the RM of MIGHTEE-POL sources as a function of the impact parameter from the intervening galaxy. In addition, we examined the dependence of the RM on redshift. We then repeated this procedure using a deeper galaxy catalog with photometric redshifts.Results.For the spectroscopic star-forming sample, we find a redshift-corrected |RM| excess of 5.6 ± 2.3 rad m−2which corresponds to a 2.5σsignificance around galaxies with a median redshift ofz = 0.46 for impact parameters below 130 kpc only selecting the intervenor with the smallest impact parameter. Making use of a photometric galaxy catalog and taking into account all intervenors withMg < −13.6 mag, the signal disappears. We find no indication for a correlation between redshift and RM, nor do we find a connection between the total number of intervenors to the total |RM|.Conclusions.We have presented tentative evidence that the CGM of star-forming galaxies is permeated by coherent magnetic fields within the virial radius. We conclude that mostly bright, star-forming galaxies with impact parameters less than 130 kpc significantly contribute to the RM of the background radio source.

Identifying Lyα emitter candidates with Random Forest: Learning from galaxies in the CANDELS survey

The physical processes that make a galaxy a Lyman alpha emitter have been extensively studied over the past 25 yr. However, the correlations between physical and morphological properties of galaxies and the strength of the Lyα emission line are still highly debated. Here, we investigate the correlations between the rest-frame Lyα equivalent width and stellar mass, star formation rate, dust reddening, metallicity, age, half-light semi-major axis, Sérsic index, and projected axis ratio in a sample of 1578 galaxies in the redshift range of 2 ≤ z ≤ 7.9 from the GOODS-S, UDS, and COSMOS fields. From the large sample of Lyα emitters (LAEs) in the dataset, we find that LAEs are typically common main sequence (MS) star-forming galaxies that show a stellar mass ≤109 M⊙, star formation rate ≤ 100.5 M⊙ yr−1, E(B − V)≤0.2, and half-light semi-major axis ≤1 kpc. Building on these findings, we have developed a new method based on a random forest (RF) machine learning (ML) classifier to select galaxies with the highest probability of being Lyα emitters. When applied to a population in the redshift range z ∈ [2.5, 4.5], our classifier holds a (80 ± 2)% accuracy and (73 ± 4)% precision. At higher redshifts (z ∈ [4.5, 6]), we obtained an accuracy of 73% and precision of 80%. These results highlight the possibility of overcoming the current limitations in assembling large samples of LAEs by making informed predictions that can be used for planning future large-scale spectroscopic surveys.

A Value-added COSMOS2020 Catalog of Physical Properties: Constraining Temperature-dependent Initial Mass Function

This work presents and releases a catalog of new photometrically derived physical properties for the ∼105 most well-measured galaxies in the COSMOS field on the sky. Using a recently developed technique, spectral energy distributions are modeled assuming a stellar initial mass function (IMF) that depends on the temperature of gas in star-forming regions. The method is applied to the largest current sample of high-quality panchromatic photometry, the COSMOS2020 catalog, that allows for testing this assumption. It is found that the galaxies exhibit a continuum of IMF and gas temperatures, most of which are bottom-lighter than measured in the Milky Way. As a consequence, the stellar masses and star formation rates of most galaxies here are found to be lower than those measured by traditional techniques in the COSMOS2020 catalog by factors of ∼1.6–3.5 and 2.5–70.0, respectively, with the change being the strongest for the most active galaxies. The resulting physical properties provide new insights into variation of the IMF-derived gas temperature along the star-forming main sequence and at quiescence, produce a sharp and coherent picture of downsizing, as seen from the stellar mass functions, and hint at a possible high-temperature and high-density stage of early galactic evolution.

COSMOS-Web: An Overview of the JWST Cosmic Origins Survey

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hr treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg2 NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5σ point-source depths ranging ∼27.5–28.2 mag. In parallel, we will obtain 0.19 deg2 of MIRI imaging in one filter (F770W) reaching 5σ point-source depths of ∼25.3–26.0 mag. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization (6 ≲ z ≲ 11) and map reionization’s spatial distribution, environments, and drivers on scales sufficiently large to mitigate cosmic variance, (2) to identify hundreds of rare quiescent galaxies at z > 4 and place constraints on the formation of the universe’s most-massive galaxies (M ⋆ > 1010 M ⊙), and (3) directly measure the evolution of the stellar-mass-to-halo-mass relation using weak gravitational lensing out to z ∼ 2.5 and measure its variance with galaxies’ star formation histories and morphologies. In addition, we anticipate COSMOS-Web’s legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool subdwarf stars in the Galactic halo, and possibly the identification of z > 10 pair-instability supernovae. In this paper we provide an overview of the survey’s key measurements, specifications, goals, and prospects for new discovery.

Elentári:a massive proto-supercluster at z ∼ 3.3 in the cosmos field

ABSTRACT Motivated by spectroscopic confirmation of three overdense regions in the COSMOS field at z ∼ 3.35, we analyse the uniquely deep multiwavelength photometry and extensive spectroscopy available in the field to identify any further related structure. We construct a three-dimensional density map using the Voronoi tesselation Monte Carlo method and find additional regions of significant overdensity. Here, we present and examine a set of six overdense structures at 3.20 < z < 3.45 in the COSMOS field, the most well-characterized of which, PCl J0959 + 0235, has 80 spectroscopically confirmed members and an estimated mass of 1.35 × 1015 M⊙, and is modelled to virialize at z ∼ 1.5−2.0. These structures contain 10 overdense peaks with >5σ overdensity separated by up to 70 cMpc, suggestive of a proto-supercluster similar to the Hyperion system at z ∼ 2.45. Upcoming photometric surveys with JWST such as COSMOS-Web, and further spectroscopic follow-up will enable more extensive analysis of the evolutionary effects that such an environment may have on its component galaxies at these early times.

CONCERTO: Extracting the power spectrum of the [CII] emission line

Context. CONCERTO is the first experiment to perform a [CII] line intensity mapping (LIM) survey on the COSMOS field to target z > 5.2. Measuring the [CII] angular power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the epochs of Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds: the dust continuum emission and atomic and molecular lines from foreground galaxies at z ≲ 3. Aims. We evaluate our ability to retrieve the [CII] signal in mock observations of the sky using the Simulated Infrared Dusty Extragalactic Sky (SIDES), which covers the mid-infrared to millimetre range. We also measure the impact of field-to-field variance on the residual foreground contamination. Methods. We compared two methods for dealing with the dust continuum emission from galaxies (i.e. the cosmic infrared background fluctuations): the standard principal component analysis (PCA) and the asymmetric re-weighted penalized least-squares (arPLS) method. For line interlopers, the strategy relies on masking low-redshift galaxies using the instrumental beam profile and external catalogues. As we do not have observations of CO or deep-enough classical CO proxies (such as LIR), we relied on the COSMOS stellar mass catalogue, which we demonstrate to be a reliable CO proxy for masking. To measure the angular power spectrum of masked data, we adapted the P of K EstimatoR (POKER) from cosmic infrared background studies and discuss its use on LIM data. Results. The arPLS method achieves a reduction in the cosmic infrared background fluctuations to a sub-dominant level of the [CII] power at z ∼ 7, a factor of > 70 below our fiducial [CII] model. When using the standard PCA, this factor is only 0.7 at this redshift. The masking lowers the power amplitude of line contamination down to 2 × 10−2 Jy2 sr−1. This residual level is dominated by faint undetected sources that are not clustered around the detected (and masked) sources. For our [CII] model, this results in a detection at z = 5.2 with a power ratio [CII]/(residual interlopers) = 62 ± 32 for a 22% area survey loss. However, at z = 7, [CII]/(residual interlopers) = 2.0 ± 1.4, due to the weak contrast between [CII] and the residual line contamination. Thanks to the large area covered by SIDES-Uchuu, we show that the power amplitude of line residuals varies by 12–15% for z = 5.2 − 7, which is less than the field-to-field variance affecting [CII] power spectra. Conclusions. We present an end-to-end simulation of the extragalactic foreground removal that we ran to detect the [CII] at high redshift via its angular power spectrum. We show that cosmic infrared background fluctuations are not a limiting foreground for [CII] LIM. On the contrary, the CO and [CI] line contamination severely limits our ability to accurately measure the [CII] angular power spectrum at z ≳ 7.

Lyman α emitters in ionized bubbles: constraining the environment and ionized fraction

ABSTRACT Lyman α emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to overdense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assumed to mark overdense, ionized bubbles. Here, we develop the first framework to quantify the ionization and density fields of high-z galaxy associations. We explore the interplay between (i) the large-scale density of a survey field, (ii) Poisson noise due to the small number density of bright sources at high redshifts (z ∼ 7), and (iii) the effects of the ionized fraction on the observation of LAEs. We use Bayesian statistics, a simple model of reionization, and a Monte Carlo simulation to construct a more comprehensive method for calculating the large-scale density of LAE regions than previous works. We find that Poisson noise has a strong effect on the inferred density of a region and show how the ionized fraction can be inferred. We then apply our framework to the strongest association yet identified: Hu et al. found 14 LAEs in a volume of ∼50 000 cMpc3 inside the COSMOS field at z ∼ 7. We show that this is most likely a 2.4σ overdensity inside of an ionized or nearly ionized bubble. We also show that this LAE association implies that the global ionized fraction is $\bar{Q} = 0.59^{+0.10}_{-0.11}$, within the context of a simple reionization model.

The PAU survey: classifying low-z SEDs using Machine Learning clustering

ABSTRACT We present an application of unsupervised Machine Learning clustering to the PAU survey of galaxy spectral energy distribution (SED) within the COSMOS field. The clustering algorithm is implemented and optimized to get the relevant groups in the data SEDs. We find 12 groups from a total number of 5234 targets in the survey at 0.01 < z < 0.28. Among the groups, 3545 galaxies (68 per cent) show emission lines in the SEDs. These groups also include 1689 old galaxies with no active star formation. We have fitted the SED to every single galaxy in each group with CIGALE. The mass, age, and specific star formation rates (sSFR) of the galaxies range from 0.15 < age/Gyr <11; 6 < log (M⋆/M⊙) <11.26, and −14.67 < log (sSFR/yr−1) <−8. The groups are well-defined in their properties with galaxies having clear emission lines also having lower mass, are younger and have higher sSFR than those with elliptical like patterns. The characteristic values of galaxies showing clear emission lines are in agreement with the literature for starburst galaxies in COSMOS and GOODS-N fields at low redshift. The star-forming main sequence, sSFR versus stellar mass and UVJ diagram show clearly that different groups fall into different regions with some overlap among groups. Our main result is that the joint of low- resolution (R ∼ 50) photometric spectra provided by the PAU survey together with the unsupervised classification provides an excellent way to classify galaxies. Moreover, it helps to find and extend the analysis of extreme ELGs to lower masses and lower SFRs in the local Universe.

The slippery slope of dust attenuation curves

Aims. We investigate the dust attenuation of 122 heavily dust-obscured galaxies detected with the Atacama Large Millimeter Array (ALMA) and Herschel in the COSMOS field. We search for correlations between dust attenuation recipes and the variation of physical parameters, namely, the effective radii of galaxies, their star formation rates, and stellar masses. We aim to understand which of the commonly used laws best describes dust attenuation in dusty star-forming galaxies (DSFGs) at high redshift. Methods. We made use of the extensive photometric coverage of the COSMOS data combined with highly resolved dust continuum maps from ALMA. We used CIGALE to estimate various physical properties of these dusty objects, namely: their star formation rates (SFR), their stellar masses, and their attenuation at short wavelengths. We inferred the effective radii (Re) of galaxies using GALFIT in the Y band of HSC and ALMA continuum maps. We used these radii to investigate the relative compactness of the dust continuum and the extension of the rest-frame UV/optical Re(y)/Re(ALMA). Results. We find that the physical parameters calculated from our models strongly depend on the assumption of the dust attenuation curve. As expected, the most impacted parameter is the stellar mass, which leads to a change in the “starburstiness” of the objects. We find that taking into account the relative compactness of star-to-dust emission prior to SED fitting is crucial, especially when studying dust attenuation of dusty star-forming galaxies. Shallower attenuation curves do not show a clear preference of compactness with attenuation, while the Calzetti attenuation curve is shown to prefer a comparable spatial extent of unattenuated stellar light and dust emission. The evolution of the Re(UV)/Re(ALMA) ratio with redshift peaks around the cosmic noon in our sample of DSFGs, showing that this compactness is correlated with the cosmic SFR density of these dusty sources.

Searching for Intragroup Light in Deep U-band Imaging of the COSMOS Field

We present the results of deep, ground based U-band imaging with the Large Binocular Telescope of the Cosmic Evolution Survey field as part of the near-UV imaging program, UVCANDELS. We utilize a seeing sorted stacking method along with night-to-night relative transparency corrections to create optimal depth and optimal resolution mosaics in the U-band, which are capable of reaching point source magnitudes of AB∼26.5 mag at 3σ. These ground-based mosaics bridge the wavelength gap between the Hubble Space Telescope WFC3 F275W and ACS F435W images and are necessary to understand galaxy assembly in the last 9–10 Gyr. We use the depth of these mosaics to search for the presence of U-band intragroup light (IGrL) beyond the local universe. Regardless of how groups are scaled and stacked, we do not detect any U-band IGrL to unprecedented U-band depths of ∼29.1–29.6 mag arcsec−2, which corresponds to an IGrL fraction of ≲1% of the total group light. This stringent upper limit suggests that IGrL does not contribute significantly to the Extragalactic Background Light at short wavelengths. Furthermore, the lack of UV IGrL observed in these stacks suggests that the atomic gas observed in the intragroup medium is likely not dense enough to trigger star formation on large scales. Future studies may detect IGrL by creating similar stacks at longer wavelengths or by pre-selecting groups which are older and/or more dynamically evolved similar to past IGrL observations of compact groups and loose groups with signs of gravitational interactions.

Photometric IGM tomography with Subaru/HSC: the large-scale structure of Lyαemitters and IGM transmission in the COSMOS field atz˜ 5

ABSTRACTWe present a novel technique called ‘photometric IGM tomography’ to map the intergalactic medium (IGM) at z ≃ 4.9 in the COSMOS field. It utilizes deep narrow-band (NB) imaging to photometrically detect faint Ly α forest transmission in background galaxies across the Subaru/Hyper-Suprime Cam (HSC)’s $1.8\rm \, sq.\, deg$ field of view and locate Ly α emitters (LAEs) in the same cosmic volume. Using ultra-deep HSC images and Bayesian spectral energy distribution fitting, we measure the Ly α forest transmission at z ≃ 4.9 along a large number (140) of background galaxies selected from the DEIMOS10k spectroscopic catalogue at 4.98 < z < 5.89 and the SILVERRUSH LAEs at z ≃ 5.7. We photometrically measure the mean Ly α forest transmission and achieve a result consistent with previous measurements based on quasar spectra. We also measure the angular LAE-Ly α forest cross-correlation and Ly α forest autocorrelation functions and place an observational constraint on the large-scale fluctuations of the IGM around LAEs at z ≃ 4.9. Finally, we present the reconstructed 2D tomographic map of the IGM, co-spatial with the large-scale structure of LAEs, at a transverse resolution of $11 \, h^{-1}\rm cMpc$ across $140\, h^{-1}\rm cMpc$ in the COSMOS field at z ≃ 4.9. We discuss the observational requirements and the potential applications of this new technique for understanding the sources of reionization, quasar radiative history, and galaxy–IGM correlations across z ∼ 3–6. Our results represent the first proof-of-concept of photometric IGM tomography, offering a new route to examining early galaxy evolution in the context of the large-scale cosmic web from the epoch of reionization to cosmic noon.

Comparing the host galaxy ages of X-ray selected AGN in COSMOS

We explore the properties of the host galaxies of X-ray selected AGN in the COSMOS field using the Chandra Legacy sample and the LEGA-C survey VLT optical spectra. Our main goal is to compare the relative ages of the host galaxies of the obscured and unobscured AGN by means of the calcium break Dn(4000) and the Hδ Balmer line. The host galaxy ages are examined in conjunction with other properties such as the galaxy stellar mass, and star-formation rate as well as the AGN Eddington ratio. Our sample consists of 50 unobscured or mildly obscured (NH < 1023 cm−2) and 23 heavily obscured AGN (NH > 1023 cm−2) in the redshift range z = 0.6 − 1. We take specific caution to create control samples in order to match the exact luminosity and redshift distributions for the obscured and unobscured AGN. The majority of unobscured AGN appear to live in young galaxies in contrast to the obscured AGN which appear to live in galaxies located between the young and old galaxy populations. This finding may be in contrast to those evolutionary AGN unification models which postulate that the AGN begin their life in a heavy obscuration phase. The host galaxies of the obscured AGN have significantly lower levels of specific star-formation. At the same time the obscured AGN have lower Eddington ratios indicating a link between the star-formation and the black hole accretion. We find that the distribution of the stellar masses of the host galaxies of obscured AGN is skewed towards higher stellar masses in agreement with previous findings. Our results on the relative age of obscured AGN are valid when we match our obscured and unobscured AGN samples according to the stellar mass of their host galaxies. All the above results become less conspicuous when a lower column density (log NH(cm−2) = 21.5 or 22) is used to separate the obscured and unobscured AGN populations.

Probabilistic mass-mapping with neural score estimation

Context.Weak lensing mass-mapping is a useful tool for accessing the full distribution of dark matter on the sky, but because of intrinsic galaxy ellipticies, finite fields, and missing data, the recovery of dark matter maps constitutes a challenging, ill-posed inverse problemAims.We introduce a novel methodology that enables the efficient sampling of the high-dimensional Bayesian posterior of the weak lensing mass-mapping problem, relying on simulations to define a fully non-Gaussian prior. We aim to demonstrate the accuracy of the method to simulated fields, and then proceed to apply it to the mass reconstruction of the HST/ACS COSMOS field.Methods.The proposed methodology combines elements of Bayesian statistics, analytic theory, and a recent class of deep generative models based on neural score matching. This approach allows us to make full use of analytic cosmological theory to constrain the 2pt statistics of the solution, to understand any differences between this analytic prior and full simulations from cosmological simulations, and to obtain samples from the full Bayesian posterior of the problem for robust uncertainty quantification.Results.We demonstrate the method in theκTNG simulations and find that the posterior mean significantly outperfoms previous methods (Kaiser–Squires, Wiener filter, Sparsity priors) both for the root-mean-square error and in terms of the Pearson correlation. We further illustrate the interpretability of the recovered posterior by establishing a close correlation between posterior convergence values and the S/N of the clusters artificially introduced into a field. Finally, we apply the method to the reconstruction of the HST/ACS COSMOS field, which yields the highest-quality convergence map of this field to date.Conclusions.We find the proposed approach to be superior to previous algorithms, scalable, providing uncertainties, and using a fully non-Gaussian prior.

The hierarchical clustering method: abundance and properties of local satellite populations

ABSTRACT The faint satellites of the local Universe provide an important benchmark for our understanding of structure formation and galaxy formation, but satellite populations are hard to identify beyond the Local Group. We recently developed an iterative method to quantify satellite abundance using galaxy clustering and tested it on a local sample in the COSMOS field, where accurate photometric redshifts are available for a large number of faint objects. In this paper, we consider the properties of these satellite populations in more detail, studying the satellite stellar mass function (SSMF), the satellite-central connection, and quenching as a function of satellite and central mass and colour. Despite the limited sample size, our results show good consistency with those from much larger surveys and constrain the SSMF down to some of the lowest primary masses considered to date. We reproduce several known trends in satellite abundance and quenching, and find evidence for one new one, a dependence of the quiescent fraction on the primary-to-secondary halo mass ratio. We discuss the prospects for the clustering method in current and forthcoming surveys.

The PAU Survey and Euclid: Improving broadband photometric redshifts with multi-task learning

Current and future imaging surveys require photometric redshifts (photo-zs) to be estimated for millions of galaxies. Improving the photo-z quality is a major challenge but is needed to advance our understanding of cosmology. In this paper we explore how the synergies between narrow-band photometric data and large imaging surveys can be exploited to improve broadband photometric redshifts. We used a multi-task learning (MTL) network to improve broadband photo-z estimates by simultaneously predicting the broadband photo-z and the narrow-band photometry from the broadband photometry. The narrow-band photometry is only required in the training field, which also enables better photo-z predictions for the galaxies without narrow-band photometry in the wide field. This technique was tested with data from the Physics of the Accelerating Universe Survey (PAUS) in the COSMOS field. We find that the method predicts photo-zs that are 13% more precise down to magnitude iAB < 23; the outlier rate is also 40% lower when compared to the baseline network. Furthermore, MTL reduces the photo-z bias for high-redshift galaxies, improving the redshift distributions for tomographic bins with z > 1. Applying this technique to deeper samples is crucial for future surveys such as Euclid or LSST. For simulated data, training on a sample with iAB < 23, the method reduces the photo-z scatter by 16% for all galaxies with iAB < 25. We also studied the effects of extending the training sample with photometric galaxies using PAUS high-precision photo-zs, which reduces the photo-z scatter by 20% in the COSMOS field.

COSMOS2020: Discovery of a Protocluster of Massive Quiescent Galaxies at z = 2.77

Protoclusters of galaxies have been found in the last quarter-century. However, most of them have been found through the overdensity of star-forming galaxies, and there have been no known structures identified by more than two spectroscopically confirmed quiescent galaxies at z > 2.5. In this letter, we report the discovery of an overdense structure of massive quiescent galaxies with the spectroscopic redshift z = 2.77 in the COSMOS field, QO-1000. We first photometrically identify this structure as a 4.2σ overdensity with 14 quiescent galaxies in 7 × 4 pMpc2 from the COSMOS2020 catalog. We then securely confirm the spectroscopic redshifts of four quiescent galaxies by detecting multiple Balmer absorption lines with Keck/MOSFIRE. All the spectroscopically confirmed members are massive () and located in a narrow redshift range (2.76 < z < 2.79). Moreover, three of them are in the 1 × 1 pMpc2 in the transverse direction at the same redshift (z = 2.760–2.763). Such a concentration of four spectroscopically confirmed quiescent galaxies implies that QO-1000 is >68 times denser than the general field. In addition, we confirm that they form a red sequence in the J − K s color. This structure’s halo mass is estimated as from its stellar mass. Similar structures found in the IllustrisTNG simulation are expected to evolve into massive galaxy clusters with at z = 0. These results suggest that QO-1000 is a more mature protocluster than the other known protoclusters. It is likely in a transition phase between star-forming protoclusters and quenched galaxy clusters.

Observational Evidence for Large-scale Gas Heating in a Galaxy Protocluster at z = 2.30

We report a z = 2.30 galaxy protocluster (COSTCO-I) in the COSMOS field, where the Lyα forest as seen in the CLAMATO IGM tomography survey does not show significant absorption. This departs from the transmission–density relationship (often dubbed the fluctuating Gunn–Peterson approximation; FGPA) usually expected to hold at this epoch, which would lead one to predict strong Lyα absorption at the overdensity. For comparison, we generate mock Lyα forest maps by applying the FGPA to constrained simulations of the COSMOS density field and create mocks that incorporate the effects of finite sight-line sampling, pixel noise, and Wiener filtering. Averaged over r = 15 h −1 Mpc around the protocluster, the observed Lyα forest is consistently more transparent in the real data than in the mocks, indicating a rejection of the null hypothesis that the gas in COSTCO-I follows the FGPA (p = 0.0026, or 2.79σ significance). It suggests that the large-scale gas associated with COSTCO-I is being heated above the expectations of the FGPA, which might be due to either large-scale AGN jet feedback or early gravitational shock heating. COSTCO-I is the first known large-scale region of the IGM that is observed to be transitioning from the optically thin photoionized regime at cosmic noon to eventually coalesce into an intracluster medium (ICM) by z = 0. Future observations of similar structures will shed light on the growth of the ICM and allow constraints on AGN feedback mechanisms.

ALMA confirmation of an obscured hyperluminous radio-loud AGN at z = 6.853 associated with a dusty starburst in the 1.5 deg2 COSMOS field

ABSTRACT We present band 6 ALMA observations of a heavily obscured radio-loud (L1.4 GHz = 1025.4 W Hz−1) active galactic nucleus (AGN) candidate at zphot = 6.83 ± 0.06 found in the 1.5 deg2 COSMOS field. The ALMA data reveal detections of exceptionally strong [C ii]158 $\mu$m (z[C ii] = 6.8532) and underlying dust continuum emission from this object (COS-87259), where the [C ii] line luminosity, line width, and 158 $\mu$m continuum luminosity are comparable to those seen from z ∼ 7 sub-mm galaxies and quasar hosts. The 158 $\mu$m continuum detection suggests a total infrared luminosity of $9\times 10^{12}\, \mathrm{ L}_\odot$ with corresponding very large obscured star formation rate (1300 M⊙ yr−1) and dust mass ($2\times 10^9\, \mathrm{ M}_\odot$). The strong break seen between the VIRCam and IRAC photometry perhaps suggests that COS-87259 is an extremely massive reionization-era galaxy with $M_\ast \approx 1.7\times 10^{11}\, \mathrm{ M}_\odot$. Moreover, the MIPS, PACS, and SPIRE detections imply that this object harbours an AGN that is heavily obscured ($\tau _{_{\mathrm{9.7\,\mu m}}}=2.3$) with a bolometric luminosity of approximately $5\times 10^{13}\, \mathrm{ L}_\odot$. Such a very high AGN luminosity suggests that this object is powered by an ≈1.6 × 10$^9\, \mathrm{ M}_\odot$ black hole if accreting near the Eddington limit, and is effectively a highly obscured version of an extremely ultraviolet (UV)-luminous (M1450 ≈ −27.3) z ∼ 7 quasar. Notably, these z ∼ 7 quasars are an exceedingly rare population (∼0.001 deg−2), while COS-87259 was identified over a relatively small field. Future very wide area surveys with e.g. Roman and Euclid have the potential to identify many more extremely red yet UV-bright z ≳ 7 objects similar to COS-87259, providing richer insight into the occurrence of intense obscured star formation and supermassive black hole growth among this population.

GALFIT-ing AGN Host Galaxies in COSMOS: HST versus Subaru

The COSMOS field has been extensively observed by most major telescopes, including Chandra, HST, and Subaru. HST imaging boasts very high spatial resolution and is used extensively in morphological studies of distant galaxies. Subaru provides lower spatial resolution imaging than HST but a substantially wider field of view with greater sensitivity. Both telescopes provide near-infrared imaging of COSMOS. Successful morphological fitting of Subaru data would allow us to measure morphologies of over 104 known active galactic nucleus (AGN) hosts, accessible through Subaru wide-field surveys, currently not covered by HST. The morphological parameters indicate the types of galaxies that host AGNs. For 4016 AGNs between 0.03 < z < 6.5, we study the morphology of their galaxy hosts using GALFIT, fitting components representing the AGN and host galaxy simultaneously using the i-band imaging from both HST and Subaru. Comparing the fits for the differing telescope spatial resolutions and image signal-to-noise ratios, we identify parameter regimes for which there is strong disagreement between distributions of fitted parameters for HST and Subaru. In particular, the Sérsic index values strongly disagree between the two sets of data, including sources at lower redshifts. In contrast, the measured magnitude and radius parameters show reasonable agreement. Additionally, large variations in the Sérsic index have little effect on the of each fit, whereas variations in other parameters have a more significant effect. These results indicate that the Sérsic index distributions of high-redshift galaxies that host AGNs imaged at ground-based spatial resolution are not reliable indicators of galaxy type and should be interpreted with caution.

Relaxed blue ellipticals: accretion-driven stellar growth is a key evolutionary channel for low mass elliptical galaxies

ABSTRACT How elliptical galaxies form is a key question in observational cosmology. While the formation of massive ellipticals is strongly linked to mergers, the low mass (M⋆/M⊙ &amp;lt; 109.5) regime remains less well explored. In particular, studying elliptical populations when they are blue, and therefore rapidly building stellar mass, offers strong constraints on their formation. Here, we study 108 blue low-mass ellipticals (which have a median stellar mass of 108.7 M⊙) at z &amp;lt; 0.3 in the COSMOS field. Visual inspection of extremely deep optical HSC images indicates that less than 3 per cent of these systems have visible tidal features, a factor of 2 less than the incidence of tidal features in a control sample of galaxies with the same distribution of stellar mass and redshift. This suggests that the star formation activity in these objects is not driven by mergers or interactions but by secular gas accretion. We combine accurate physical parameters from the COSMOS2020 catalogue, with measurements of local density and the locations of galaxies in the cosmic web, to show that our blue ellipticals reside in low-density environments, further away from nodes and large-scale filaments than other galaxies. At similar stellar masses and environments, blue ellipticals outnumber their normal (red) counterparts by a factor of 2. Thus, these systems are likely progenitors of not only normal ellipticals at similar stellar mass but, given their high star formation rates, also of ellipticals at higher stellar masses. Secular gas accretion therefore likely plays a significant (and possibly dominant) role in the stellar assembly of elliptical galaxies in the low-mass regime.