Search

Our knowledge of relations between supermassive black holes and their host galaxies at z ≳ 1 is still limited, even though being actively sought out to z ∼ 6. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 107 X-ray-selected type-I active galactic nuclei (AGNs) at 0.68

We present the properties of two candidate massive ( M _⋆ ∼ 10 ^11 M _⊙ ) and dusty ( A _v > 2.5 mag) galaxies at z = 5–7 in the first 0.28 deg ^2 of the COSMOS-Web survey. One object is spectroscopically confirmed at z _spec = 5.051, while the other has a robust z _phot = 6.7 ± 0.3. Thanks to their extremely red colors (F277W–F444W ∼ 1.7 mag), these galaxies satisfy the nominal color selection for the widely studied “little red dot” (LRD) population with the exception of their spatially resolved morphologies. The morphology of our targets allows us to conclude that their red continuum is dominated by highly obscured stellar emission and not by reddened nuclear activity. Using a variety of spectral energy distribution fitting tools and star formation histories, we estimate the stellar masses to be $\mathrm{log}({M}_{\star })={11.32}_{-0.15}^{+0.07}\,{M}_{\odot }$ and $\mathrm{log}({M}_{\star })={11.2}_{-0.2}^{+0.1}\,{M}_{\odot }$ , respectively, with a red continuum emission dominated by a recent episode of star formation. We then compare their number density to the halo mass function to infer stellar baryon fractions of ϵ _⋆ ∼ 0.25 and ϵ _⋆ ∼ 0.5. Both are significantly higher than what is commonly observed in lower- z galaxies or more dust-obscured galaxies at similar redshifts. With very bright ultra-high- z Lyman-Break Galaxies and some non-AGN-dominated LRDs, such “extended” LRDs represent another population that may require very efficient star formation at early times.

The role of internal and environmental factors in the star formation activity of galaxies is still a matter of debate, in particular at higher redshifts. Leveraging the most recent release of the COSMOS catalog, COSMOS2020, as well as density measurements from our previous study, we disentangle the impact of environment and stellar mass on the star formation rate (SFR) and specific SFR (sSFR) of a sample of ∼210,000 galaxies within a redshift range of 0.4 < z < 4, and present our findings in three cosmic epochs: (1) out to z ∼ 1, the average SFR and sSFR decline in extremely dense environments and at the high-mass end of the distribution, which is mostly due to the presence of the massive quiescent population; (2) at 1 < z < 2, the environmental dependence diminishes, while mass is still the dominant factor in star formation activity; and (3) beyond z ∼ 2, our sample is dominated by star-forming galaxies and we observe a reversal of the trends seen in the local Universe—the average SFR increases with increasing environmental density. Our analysis shows that both environmental and mass quenching efficiencies increase with stellar mass at all redshifts, with mass being the dominant quenching factor in massive galaxies compared to environmental effects. At 2 < z < 4, negative values of environmental quenching efficiency suggest that the fraction of star-forming galaxies in dense environments exceeds that in less-dense regions, likely due to the greater availability of cold gas, higher merger rates, and tidal effects that trigger star formation activity.

We present a statistical study of 180 dust continuum sources identified in 33 massive cluster fields by the Atacama Large Millimeter/submillimeter Array Lensing Cluster Survey (ALCS) over a total of 133 arcmin ^2 area, homogeneously observed at 1.2 mm. ALCS enables us to detect extremely faint millimeter sources by lensing magnification, including near-infrared (NIR) dark objects showing no counterparts in existing Hubble Space Telescope and Spitzer images. The dust continuum sources belong to a blind sample ( N = 141) with signal-to-noise ratio (S/N) ≳ 5.0 (a purity of >0.99) or a secondary sample ( N = 39) with S/N = 4.0–5.0 screened by priors. With the blind sample, we securely derive 1.2 mm number counts down to ∼7 μ Jy, and find that the total integrated 1.2 mm flux is ${20.7}_{-6.5}^{+8.5}$ Jy deg ^−2 , resolving ≃80% of the cosmic infrared background light. The resolved fraction varies by a factor of 0.6–1.1 due to the completeness correction depending on the spatial size of the millimeter emission. We also derive infrared (IR) luminosity functions (LFs) at z = 0.6–7.5 with the $1/{V}_{\max }$ method, finding the redshift evolution of IR LFs characterized by positive luminosity and negative density evolution. The total (= UV + IR) cosmic star formation rate density (SFRD) at z > 4 is estimated to be ${161}_{-21}^{+25}$ % of the Madau and Dickinson measurements mostly based on rest-frame UV surveys. Although our general understanding of the cosmic SFRD is unlikely to change beyond a factor of 2, these results add to the weight of evidence for an additional (≈60%) SFRD component contributed by the faint millimeter population, including NIR-dark objects.

We report on the spectroscopic confirmation of a massive quiescent galaxy at z _spec = 4.53 in the COSMOS field. The object was first identified as a galaxy with suppressed star formation at z _phot ∼ 4.65 from the COSMOS2020 catalog. The follow-up spectroscopy with Keck/MOSFIRE in the K band reveals faint [O ii ] emission and the Balmer break, indicative of evolved stellar populations. We fit the spectral energy distribution using photometry and a spectrum to infer physical properties. The obtained stellar mass is high ( M _* ∼ 10 ^10.8 M _⊙ ) and the current star formation rate is more than 1 dex below that of main-sequence galaxies at z = 4.5. Its star formation history suggests that this galaxy experienced rapid quenching from z ∼ 5. The galaxy is among the youngest quiescent galaxies confirmed so far at z _spec > 3 with z _form ∼ 5.2 (200 Myr ago), which is the epoch when 50% of the total stellar mass was formed. A unique aspect of the galaxy is that it is in an extremely dense region; there are four massive star-forming galaxies at 4.4 < z _phot < 4.7 located within 150 physical kpc from the galaxy. Interestingly, three of them have virial radii that strongly overlap with that of the central quiescent galaxy (∼70 kpc), suggesting that the overdensity region is likely the highest-redshift candidate of a dense group with a spectroscopically confirmed quiescent galaxy at the center. The group provides us with a unique opportunity to gain insights into the role of the group environment in quenching at z ∼ 5, which corresponds to the formation epoch of massive elliptical galaxies in the local Universe.

We present Atacama Large Millimeter/submillimeter Array (ALMA) deep spectroscopy for a lensed galaxy at z _spec = 8.496 with $\mathrm{log}({M}_{\mathrm{star}}/{M}_{\odot })\sim 7.8$ whose optical nebular lines and stellar continuum are detected by JWST/NIRSpec and NIRCam Early Release Observations in the field of SMACS J0723.3–7327. Our ALMA spectrum shows [O iii ] 88 μ m and [C ii ] 158 μ m line detections at 4.0 σ and 4.5 σ , respectively. The redshift and position of the [O iii ] line coincide with those of the JWST source, while the [C ii ] line is blueshifted by 90 km s ^−1 with a spatial offset of 0.″5 (≈0.5 kpc in the source plane) from the centroid of the JWST source. The NIRCam F444W image, including [O iii ] λ 5007 and H β line emission, spatially extends beyond the stellar components by a factor of >8. This indicates that the z = 8.5 galaxy has already experienced strong outflows as traced by extended [O iii ] λ 5007 and offset [C ii ] emission, which would promote ionizing photon escape and facilitate reionization. With careful slit-loss corrections and the removal of emission spatially outside the galaxy, we evaluate the [O iii ] 88 μ m/ λ 5007 line ratio, and derive the electron density n _e by photoionization modeling to be ${220}_{-130}^{+230}$ cm ^−3 , which is comparable with those of z ∼ 2–3 galaxies. We estimate an [O iii ] 88 μ m/[C ii ] 158 μ m line ratio in the galaxy of >4, as high as those of known z ∼ 6–9 galaxies. This high [O iii ] 88 μ m/[C ii ] 158 μ m line ratio is generally explained by the high n _e as well as the low metallicity ( ${Z}_{\mathrm{gas}}/{Z}_{\odot }={0.04}_{-0.02}^{+0.02}$ ), high ionization parameter ( $\mathrm{log}U\gt -2.27$ ), and low carbon-to-oxygen abundance ratio (log(C/O) = [−0.52: −0.24]) obtained from the JWST/NIRSpec data; further [C ii ] follow-up observations will constrain the covering fraction of photodissociation regions.

We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at z ≥ 3 using the Near-Infrared Camera (NIRCam) on board the James Webb Space Telescope (JWST). Based on a sample selected by UVJ color or NUVUVJ color, we focus on 26 quiescent galaxies with $9.8\lt \mathrm{log}({M}_{\star }/{M}_{\odot })\lt 11.4$ at 2.8 < z _phot < 4.6 with publicly available JWST data. Their sizes are constrained by fitting the Sérsic profile to all available NIRCam images. We see a negative correlation between the observed wavelength and the size and derive their size at the rest frame 0.5 μ m using size measurements in multiple bands. Our quiescent galaxies show a significant correlation between the rest-frame 0.5 μ m size and the stellar mass at z ≥ 3. The analytical fit for them at $\mathrm{log}({M}_{\star }/{M}_{\odot })\gt 10.3$ implies that our size–stellar mass relations are below those at lower redshifts, with the amplitude of ∼0.6 kpc at M _⋆ = 5 × 10 ^10 M _⊙ . This value agrees with the extrapolation of the size evolution of quiescent galaxies at z < 3 in the literature, implying that the size of quiescent galaxies increases monotonically from z ∼ 3–5. Our sample mainly comprises galaxies with bulge-like structures according to their median Sérsic index and axis ratio of n ∼ 3–4 and q ∼ 0.6–0.8, respectively. On the other hand, there is a trend of increasing fraction of galaxies with low Sérsic index at higher redshift, suggesting 3 < z < 5 might be the epoch of onset of morphological transformation with a fraction of very notable disky quenched galaxies.

We report the discovery of 15 exceptionally luminous 10 ≲ z ≲ 14 candidate galaxies discovered in the first 0.28 deg ^2 of JWST/NIRCam imaging from the COSMOS-Web survey. These sources span rest-frame UV magnitudes of −20.5 > M _UV > −22, and thus constitute the most intrinsically luminous z ≳ 10 candidates identified by JWST to date. Selected via NIRCam imaging, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine z ≳ 10 sources and 3/15 (20%) likely low-redshift contaminants. Three of our z ∼ 12 candidates push the limits of early stellar mass assembly: they have estimated stellar masses ∼ 5 × 10 ^9 M _⊙ , implying an effective stellar baryon fraction of ϵ _⋆ ∼ 0.2−0.5, where ϵ _⋆ ≡ M _⋆ /( f _b M _halo ). The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales < 100 Myr where the star formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for M _⋆ ∼ 10 ^10 M _⊙ galaxies relative to M _⋆ ∼ 10 ^9 M _⊙ —both about 10 ^−6 Mpc ^−3 —implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UV luminosity function from a double power law to a Schechter function at z ≈ 8. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understand how, and if, such early massive galaxies push the limits of galaxy formation in the Lambda cold dark matter paradigm.

To study the role of environment in galaxy evolution, we reconstruct the underlying density field of galaxies based on COSMOS2020 (The Farmer catalog) and provide the density catalog for a magnitude-limited ( K _s < 24.5) sample of ∼210,000 galaxies at 0.4 < z < 5 within the COSMOS field. The environmental densities are calculated using a weighted kernel density estimation approach with the choice of a von Mises–Fisher kernel, an analog of the Gaussian kernel for periodic data. Additionally, we make corrections for the edge effect and masked regions in the field. We utilize physical properties extracted by LePhare to investigate the connection between star formation activity and the environmental density of galaxies in six mass-complete subsamples at different cosmic epochs within 0.4 < z < 4. Our findings confirm a strong anticorrelation between star formation rate (SFR)/specific SFR (sSFR) and environmental density out to z ∼ 1.1. At 1.1 < z < 2, there is no significant correlation between SFR/sSFR and density. At 2 < z < 4, we observe a reversal of the SFR/sSFR–density relation such that both SFR and sSFR increase by a factor of ∼10 with increasing density contrast, δ , from −0.4 to 5. This observed reversal at higher redshifts supports the scenario where an increased availability of gas supply, along with tidal interactions and a generally higher star formation efficiency in dense environments, could potentially enhance star formation activity in galaxies located in rich environments at z > 2.

We report on discovery of a concentration of massive quiescent galaxies located at z = 4. The concentration is first identified using high-quality photometric redshifts based on deep, multiband data in Subaru/XMM-Newton Deep Field. Follow-up near-infrared spectroscopic observations with MOSFIRE on Keck confirm a massive (∼10 ^11 M _⊙ ) quiescent galaxy at z = 3.99. Our spectral energy distribution analyses reveal that the galaxy experienced an episode of starburst about 500 Myr prior to the observed epoch, followed by rapid quenching. Since its spectrum is sufficiently good to measure the stellar velocity dispersion, we infer its dynamical mass and find that it is consistent with its stellar mass. The galaxy is surrounded by four massive (>10 ^10 M _⊙ ) quiescent galaxies on a ∼1 physical Mpc scale, all of which are consistent with being located at the same redshift based on high-accuracy spectrophotometric redshifts. This is likely a (proto)cluster dominated by quiescent galaxies, the first of the kind reported at such a high redshift as z = 4. Interestingly, it is in a large-scale structure revealed by spectroscopic redshifts from VANDELS. Furthermore, it also exhibits a red sequence, adding further support to the physical concentration of the galaxies. We find no such concentration in the Illustris-TNG300 simulation; it may be that the cluster is such a rare system that the simulation box is not sufficiently large to reproduce it. The total halo mass of the quiescent galaxies is ∼10 ^13 M _⊙ , suggesting that they form a group-sized halo once they collapse together. We discuss the implications of our findings for the quenching physics and conclude with future prospects.

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 pMpc ^2 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 ( $\mathrm{log}({M}_{\star }/{M}_{\odot })\gt 11.0$ ) and located in a narrow redshift range (2.76 < z < 2.79). Moreover, three of them are in the 1 × 1 pMpc ^2 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 $\mathrm{log}({M}_{\mathrm{halo}}/{M}_{\odot })\gt 13.2$ from its stellar mass. Similar structures found in the IllustrisTNG simulation are expected to evolve into massive galaxy clusters with $\mathrm{log}({M}_{\mathrm{halo}}/{M}_{\odot })\geqslant 14.8$ 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.

We present the results of a systematic search for candidate quiescent galaxies in the distant universe in 11 JWST fields with publicly available observations collected during the first 3 months of operations and covering an effective sky area of ∼145 arcmin ^2 . We homogeneously reduce the new JWST data and combine them with existing observations from the Hubble Space Telescope. We select a robust sample of ∼80 candidate quiescent and quenching galaxies at 3 < z < 5 using two methods: (1) based on their rest-frame UVJ colors, and (2) a novel quantitative approach based on Gaussian mixture modeling of the near-UV − U , U − V , and V − J rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ( M _⋆ ≥ 10 ^10.6 M _⊙ ) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of 2–3, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at z > 3, as could be expected for highly clustered massive systems. Importantly, JWST enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available.

Atacama Large Millimeter/submillimeter Array (ALMA) observations of the long-wavelength dust continuum are used to estimate the gas masses in a sample of 708 star-forming galaxies at z = 0.3−4.5. We determine the dependence of gas masses and star formation efficiencies (SFEs; SFR per unit gas mass) on redshift (z), M _* , and star formation rate (SFR) relative to the main sequence (MS). We find that 70% of the increase in SFRs of the MS is due to the increased gas masses at earlier epochs, while 30% is due to increased efficiency of star formation (SF). For galaxies above the MS this is reversed—with 70% of the increased SFR relative to the MS being due to elevated SFEs. Thus, the major evolution of star formation activity at early epochs is driven by increased gas masses, while the starburst activity taking galaxies above the MS is due to enhanced triggering of star formation (likely due to galactic merging). The interstellar gas peaks at z = 2 and dominates the stellar mass down to z = 1.2. Accretion rates needed to maintain continuity of the MS evolution reach >100 M _⊙ yr ^−1 at z > 2. The galactic gas contents are likely the driving determinant for both the rise in SF and AGN activity from z = 5 to their peak at z = 2 and subsequent fall at lower z . We suggest that for self-gravitating clouds with supersonic turbulence, cloud collisions and the filamentary structure of the clouds regulate the star formation activity.

REQUIEM-2D (Resolving Quiescent Magnified Galaxies with 2D Grism Spectroscopy) comprises a sample of eight massive ( ${\rm{log}}({M}_{\ast }/{M}_{\odot })\gt 10.6$ ) strongly lensed quiescent galaxies at z ∼ 2. REQUIEM-2D combines the natural magnification from strong gravitational lensing with the high-spatial-resolution grism spectroscopy of the Hubble Space Telescope through a spectrophotometric fit to study spatially resolved stellar populations. We show that quiescent galaxies in the REQUIEM-2D survey have diverse formation histories with age gradients at the 1 σ –3 σ level, including examples of (1) a younger central region supporting outside-in formation, (2) flat age gradients that show evidence for both spatially uniform early formation and inside-out quenching, and (3) regions at a fixed radial distance having different ages (such asymmetries cannot be recovered when averaging stellar population measurements azimuthally). The typical dust attenuation curve for the REQUIEM-2D galaxies is constrained to be steeper than Calzetti’s law in the UV and generally consistent with A _V < 1. Combined together and accounting for the different physical radial distances and formation timescales, we find that the REQUIEM-2D galaxies that formed earlier in the universe exhibit slow and uniform growth in their inner core, whereas the galaxies that formed later have rapid inner growth in their inner core with younger ages than the outskirts. These results challenge the currently accepted paradigm of how massive quiescent galaxies form, where the earliest galaxies are thought to form most rapidly. Significantly larger samples close to the epoch of formation with similar data quality and higher spectral resolution are required to validate this finding.

We select and characterize a sample of massive (log( M _* / M _⊙ ) > 10.6) quiescent galaxies (QGs) at 3 < z < 5 in the latest Cosmological Evolution Survey catalog (COSMOS2020). QGs are selected using a new rest-frame color-selection method, based on their probability of belonging to the quiescent group defined by a Gaussian mixture model (GMM) trained on rest-frame colors ( NUV − U , U − V , V − J ) of similarly massive galaxies at 2 < z < 3. We calculate the quiescent probability threshold above which a galaxy is classified as quiescent using simulated galaxies from the shark semi-analytical model. We find that, at z ≥ 3 in shark , the GMM/ NUVU − VJ method outperforms classical rest-frame UVJ selection and is a viable alternative. We select galaxies as quiescent based on their probability in COSMOS2020 at 3 < z < 5, and compare the selected sample to both UVJ - and NUVrJ -selected samples. We find that, although the new selection matches UVJ and NUVrJ in number, the overlap between color selections is only ∼50%–80%, implying that rest-frame color commonly used at lower-redshift selections cannot be equivalently used at z > 3. We compute median rest-frame spectral energy distributions for our sample and find the median QG at 3 < z < 5 has a strong Balmer/4000 Å break, and residual NUV flux indicating recent quenching. We find the number densities of the entire quiescent population (including post-starbursts) more than doubles from 3.5 ± 2.2 × 10 ^−6 Mpc ^−3 at 4 < z < 5 to 1.4 ± 0.4 × 10 ^−5 Mpc ^−3 at 3 < z < 4, confirming that the onset of massive galaxy quenching occurs as early as 3 < z < 5.

We present the Texas Euclid Survey for Ly α (TESLA), a spectroscopic survey in the 10 deg ^2 of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Ly α emitters (LAEs) correlate with Ly α emission to understand the escape of Ly α emission from galaxies at redshifts of 2–3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 deg ^2 (H20) survey, and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby–Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust to the Ly α rest-frame equivalent width ( W _Ly _α ). We uncover marginal (1 σ significance) correlations between stellar mass and W _Ly _α , and SFR and W _Ly _α , with a Spearman correlation coefficient of −0. ${34}_{-.14}^{+.17}$ and −0. ${37}_{-.14}^{+.16}$ , respectively. We show that the W _Ly _α distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of W _0 = 150 Å. Once complete the TESLA survey will enable the study of ≳50,000 LAEs to explore more correlations between galaxy properties and W _Ly _α . The large sample size will allow the construction of a predictive model for W _Ly _α as a function of SED-derived galaxy properties, which could be used to improve Ly α -based constraints on reionization.

We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wave bands for measuring the physical properties of galaxies in a Hawaii Two-0- (H20) and UVISTA-like survey for a sample of i < 25 AB mag galaxies. We find that with the available i -band fluxes, r , u , IRAC/ ch 2, and z bands provide most of the information regarding the redshift with importance decreasing from r band to z band. We also find that for the same sample, IRAC/ ch 2, Y , r , and u bands are the most relevant bands in stellar-mass measurements with decreasing order of importance. Investigating the intercorrelation between the bands, we train a model to predict UVISTA observations in near-IR from H20-like observations. We find that magnitudes in the YJH bands can be simulated/predicted with an accuracy of 1 σ mag scatter ≲0.2 for galaxies brighter than 24 AB mag in near-IR bands. One should note that these conclusions depend on the selection criteria of the sample. For any new sample of galaxies with a different selection, these results should be remeasured. Our results suggest that in the presence of a limited number of bands, a machine-learning model trained over the population of observed galaxies with extensive spectral coverage outperforms template fitting. Such a machine-learning model maximally comprises the information acquired over available extensive surveys and breaks degeneracies in the parameter space of template fitting inevitable in the presence of a few bands.

We conduct a systematic search for protocluster candidates at z ≥ 6 in the Cosmic Evolution Survey (COSMOS) field using the recently released COSMOS2020 source catalog. We select galaxies using a number of selection criteria to obtain a sample of galaxies that have a high probability of being inside a given redshift bin. We then apply overdensity analysis to the bins using two density estimators, a Weighted Adaptive Kernel estimator and a Weighted Voronoi Tessellation estimator. We have found 15 significant (>4 σ ) candidate galaxy overdensities across the redshift range 6 ≤ z ≤ 7.7. The majority of the galaxies appear to be on the galaxy main sequence at their respective epochs. We use multiple stellar-mass-to-halo-mass conversion methods to obtain a range of dark matter halo mass estimates for the overdensities in the range of ∼10 ^11 –10 ^13 M _⊙ , at the respective redshifts of the overdensities. The number and the masses of the halos associated with our protocluster candidates are consistent with what is expected from the area of a COSMOS-like survey in a standard Λ cold dark matter cosmology. Through comparison with simulation, we expect that all of the overdensities at z ≃ 6 will evolve into Virgo-/Coma-like clusters at present (i.e., with masses ∼10 ^14 –10 ^15 M _⊙ ). Compared to other overdensities identified at z ≥ 6 via narrowband selection techniques, the overdensities presented appear to have ∼10× higher stellar masses and star formation rates (SFRs). We compare the evolution in the total SFR and stellar mass content of the protocluster candidates across the redshift range 6 ≤ z ≤ 7.7 and find agreement with the total average SFR from simulations.

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 deg ^2 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 deg ^2 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 _⋆ > 10 ^10 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.

A growing number of far-infrared (FIR) bright sources completely invisible in deep extragalactic optical surveys hint at an elusive population of z > 4 dusty, star-forming galaxies. Cycle 1 JWST surveys are now detecting their rest-frame optical light, which provides key insight into their stellar properties and statistical constraints on the population as a whole. This work presents the JWST Near Infrared Camera (NIRCam) counterpart from the COSMOS-Web survey to an FIR SCUBA-2 and Atacama Large Millimeter/submillimeter Array (ALMA) source, AzTECC71, which was previously undetected at wavelengths shorter than 850 μ m. AzTECC71, among the reddest galaxies in COSMOS-Web with F277W − F444W ∼ 0.9, is undetected in NIRCam/F150W and F115W and fainter in F444W than other submillimeter galaxies identified in COSMOS-Web by 2–4 magnitudes. This is consistent with the system having both a lower stellar mass and higher redshift than the median dusty, star-forming galaxy. With deep ground- and space-based upper limits combined with detections in F277W, F444W, and the FIR including ALMA Band 6, we find a high probability (99%) that AzTECC71 is at z > 4 with ${z}_{\mathrm{phot}}={5.7}_{-0.7}^{+0.8}$ . This galaxy is massive ( $\mathrm{log}\,{M}_{* }/{M}_{\odot }\sim 10.7$ ) and infrared-luminous ( $\mathrm{log}\,{L}_{\mathrm{IR}}/{L}_{\odot }\sim 12.7$ ), comparable to other optically undetected but FIR-bright dusty, star-forming galaxies at z > 4. This population of luminous, infrared galaxies at z > 4 is largely unconstrained but comprises an important bridge between the most extreme dust-obscured galaxies and more typical high-redshift star-forming galaxies. If further FIR-selected galaxies that drop out of the F150W filter in COSMOS-Web have redshifts z > 4 like AzTECC71, then the volume density of such sources may be ∼3–10 × greater than previously estimated.

We present a search for extremely red, dust-obscured, z > 7 galaxies with JWST/NIRCam+MIRI imaging over the first 20 arcmin ^2 of publicly available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W−F444W (∼2.5 mag) and detection in MIRI/F770W (∼25 mag), we identify two galaxies, COS-z8M1 and CEERS-z7M1, that have best-fit photometric redshifts of $z={8.4}_{-0.4}^{+0.3}$ and ${7.6}_{-0.1}^{+0.1}$ , respectively. We perform spectral energy distribution fitting with a variety of codes (including bagpipes , prospector , beagle , and cigale ) and find a >95% probability that these indeed lie at z > 7. Both sources are compact ( R _eff ≲ 200 pc) and highly obscured ( A _V ∼ 1.5–2.5) and, at our best-fit redshift estimates, likely have strong [O iii ]+H β emission contributing to their 4.4 μ m photometry. We estimate stellar masses of ∼10 ^10 M _⊙ for both sources; by virtue of detection in MIRI at 7.7 μ m, these measurements are robust to the inclusion of bright emission lines, for example, from an active galactic nucleus. We identify a marginal (2.9 σ ) Atacama Large Millimeter/submillimeter Array detection at 2 mm within 0.″5 of COS-z8M1, which, if real, would suggest a remarkably high IR luminosity of ∼10 ^12 L _⊙ . These two galaxies, if confirmed at z ∼ 8, would be extreme in their stellar and dust masses and may be representative of a substantial population of highly dust-obscured galaxies at cosmic dawn.

We present the characteristics of 2 mm selected sources from the largest Atacama Large Millimeter/submillimeter Array (ALMA) blank-field contiguous survey conducted to date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184 arcmin ^2 at 2 mm. Twelve of 13 detections above 5 σ are attributed to emission from galaxies, 11 of which are dominated by cold dust emission. These sources have a median redshift of $\langle {z}_{2\,\mathrm{mm}}\rangle ={3.6}_{-0.3}^{+0.4}$ primarily based on optical/near-infrared photometric redshifts with some spectroscopic redshifts, with 77% ± 11% of sources at z > 3 and 38% ± 12% of sources at z > 4. This implies that 2 mm selection is an efficient method for identifying the highest-redshift dusty star-forming galaxies (DSFGs). Lower-redshift DSFGs ( z < 3) are far more numerous than those at z > 3 yet are likely to drop out at 2 mm. MORA shows that DSFGs with star formation rates in excess of 300 M _⊙ yr ^−1 and a relative rarity of ∼10 ^−5 Mpc ^−3 contribute ∼30% to the integrated star formation rate density at 3 < z < 6. The volume density of 2 mm selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at z > 2. Analysis of MORA sources’ spectral energy distributions hint at steeper empirically measured dust emissivity indices than reported in typical literature studies, with $\langle \beta \rangle ={2.2}_{-0.4}^{+0.5}$ . The MORA survey represents an important step in taking census of obscured star formation in the universe’s first few billion years, but larger area 2 mm surveys are needed to more fully characterize this rare population and push to the detection of the universe’s first dusty galaxies.