Quenching and the UVJ Diagram in the SIMBA Cosmological Simulation

Over the past decade, rest-frame color–color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color–color selection tool—the rest-frame U − V versus V − J diagram—by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties in UVJ space and the color–color evolution of massive galaxies as they quench at redshifts z ∼ 1–2. We find that our models broadly reproduce the observed UVJ diagram at z = 1–2, including (for the first time in a cosmological simulation) reproducing the population of extremely dust-reddened galaxies in the top right of the UVJ diagram. However, our models primarily populate this region with low-mass galaxies and do not produce as clear a bimodality between star-forming and quiescent galaxies as is seen in observations. The former issue is due to an excess of dust in low-mass galaxies and relatively gray attenuation curves in high-mass galaxies, while the latter is due to the overpopulation of the green valley in simba. When investigating the time evolution of galaxies on the UVJ diagram, we find that the quenching pathway on the UVJ diagram is independent of the quenching timescale, and instead dependent primarily on the average specific star formation rate in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies

Accurate dust temperature determination in a z = 7.13 galaxy

We report ALMA Band 9 continuum observations of the normal, dusty star-forming galaxy A1689-zD1 at z = 7.13, resulting in a ∼4.6 σ detection at 702 GHz. For the first time, these observations probe the far-infrared spectrum shortward of the emission peak of a galaxy in the Epoch of Reionization (EoR). Together with ancillary data from earlier works, we derive the dust temperature, Td, and mass, Md, of A1689-zD1 using both traditional modified blackbody spectral energy density fitting, and a new method that relies only on the [C ii] 158 μm line and underlying continuum data. The two methods give Td = (42+13-7, 40+13-) K, and Md} = (1.7+1.3-0.7, 2.0+1.8-1.0), 7, 107, M⊙. Band 9 observations improve the accuracy of the dust temperature (mass) estimate by ∼50 per cent (6 times). The derived temperatures confirm the reported increasing Td-redshift trend between z = 0 and 8; the dust mass is consistent with a supernova origin. Although A1689-zD1 is a normal UV-selected galaxy, our results, implying that ∼85 per cent of its star-formation rate is obscured, underline the non-negligible effects of dust in EoR galaxies

Broad emission lines in optical spectra of hot dust-obscured galaxies can contribute significantly to JWST/NIRCam photometry

Selecting the first galaxies at z>7-10 from JWST surveys is complicated by z<6 contaminants with degenerate photometry. For example, strong optical nebular emission lines at z7-10 Lyman Break Galaxies (LBGs). Dust-obscured 3<z<6 galaxies in particular are potentially important contaminants, and their faint rest-optical spectra have been historically difficult to observe. A lack of optical emission line and continuum measures for 3<z<6 dusty galaxies now makes it difficult to test their expected JWST/NIRCam photometry for degenerate solutions with NIRCam dropouts. Towards this end, we quantify the contribution by strong emission lines to NIRCam photometry in a physically motivated manner by stacking 21 Keck II/NIRES spectra of hot, dust-obscured, massive ($\log\mathrm{M_*/M_\odot}\gtrsim10-11$) and infrared (IR) luminous galaxies at z~1-4. We derive an average spectrum and measure strong narrow (broad) [OIII]5007 and H$\alpha$ features with equivalent widths of $130\pm20$ A ($150\pm50$ A) and $220\pm30$ A ($540\pm80$ A) respectively. These features can increase broadband NIRCam fluxes by factors of 1.2-1.7 (0.2-0.6 mag). Due to significant dust-attenuation ($A_V\sim6$), we find H$\alpha$+[NII] to be significantly brighter than [OIII]+H$\beta$, and therefore find that emission-line dominated contaminants of high-z galaxy searches can only reproduce moderately blue perceived UV continua of $S_\lambda\propto\lambda^\beta$ with $\beta>-1.5$ and z>4. While there are some redshifts (z~3.75) where our stack is more degenerate with the photometry of z>10 LBGs between $\lambda_{rest}\sim0.3-0.8\,\mu$m, redder filter coverage beyond $\lambda_{obs}>3.5\,\mu$m and far-IR/sub-mm follow-up may be useful for breaking the degeneracy and making a crucial separation between two fairly unconstrained populations, dust-obscured galaxies at z~3-6 and LBGs at z>10.Comment: 8 pages, 3 figures, 1 table, submitted to ApJ

Extremely Red Galaxies at z = 5-9 with MIRI and NIRSpec:Dusty Galaxies or Obscured Active Galactic Nuclei?

We study a new population of extremely red objects (EROs) recently discovered by the James Webb Space Telescope (JWST) based on their NIRCam colors F277W − F444W &gt; 1.5 mag. We find 37 EROs in the Cosmic Evolution Early Release Science Survey (CEERS) field with F444W &lt; 28 mag and photometric redshifts between 5 &lt; z &lt; 7, with median z = 6.9 − 1.6 + 1.0 . Surprisingly, despite their red long-wavelength colors, these EROs have blue short-wavelength colors (F150W − F200W ∼ 0 mag) indicative of bimodal spectral energy distributions (SEDs) with a red, steep slope in the rest-frame optical, and a blue, flat slope in the rest-frame UV. Moreover, all these EROs are unresolved, point-like sources in all NIRCam bands. We analyze the SEDs of eight of them with MIRI and NIRSpec observations using stellar population models and active galactic nucleus (AGN) templates. We find that dusty galaxies or obscured AGNs provide similarly good SED fits but different stellar properties: massive and dusty, log M ⋆ / M ⊙ ∼ 10 and A V ≳ 3 mag, or low mass and obscured, log M ⋆ / M ⊙ ∼ 7.5 and A V ∼ 0 mag, hosting an obscured quasi-stellar object (QSO). SED modeling does not favor either scenario, but their unresolved sizes are more suggestive of AGNs. If any EROs are confirmed to have log M ⋆ / M ⊙ ≳ 10.5, it would increase the pre-JWST number density at z &gt; 7 by up to a factor ∼60. Similarly, if they are QSOs with luminosities in the L bol &gt; 1045-46 erg s−1 range, their number would exceed that of bright blue QSOs by more than three orders of magnitude. Additional photometry at mid-infrared wavelengths will reveal the true nature of the red continuum emission in these EROs and will place this puzzling population in the right context of galaxy evolution.</p

The Web Epoch of Reionization Lyman-α\alpha Survey (WERLS) I. MOSFIRE Spectroscopy of z∼7−8\mathbf{z \sim 7-8} Lyman-α\alpha Emitters

We present the first results from the Web Epoch of Reionization Lyman-$\alpha$ Survey (WERLS), a spectroscopic survey of Lyman-$\alpha$ emission using Keck I/MOSFIRE and LRIS. WERLS targets bright ($J<26$) galaxy candidates with photometric redshifts of $5.5\lesssim z \lesssim 8$ selected from pre-JWST imaging embedded in the Epoch of Reionization (EoR) within three JWST deep fields: CEERS, PRIMER, and COSMOS-Web. Here, we report 11 $z\sim7-8$ Lyman-$\alpha$ emitters (LAEs; 3 secure and 8 tentative candidates) detected in the first five nights of WERLS MOSFIRE data. We estimate our observed LAE yield is $\sim13$%, broadly consistent with expectations assuming some loss from redshift uncertainty, contamination from sky OH lines, and that the Universe is approximately half-ionized at this epoch, whereby observable Lyman-$\alpha$ emission is unlikely for galaxies embedded in a neutral intergalactic medium. Our targets are selected to be UV-bright, and span a range of absolute UV magnitudes with $-23.1 < M_{\text{UV}} < -19.8$. With two LAEs detected at $z=7.68$, we also consider the possibility of an ionized bubble at this redshift. Future synergistic Keck+JWST efforts will provide a powerful tool for pinpointing beacons of reionization and mapping the large scale distribution of mass relative to the ionization state of the Universe.Comment: 27 pages, 8 figures; ApJ submitte

COSMOS-Web: Intrinsically Luminous z≳\gtrsim10 Galaxy Candidates Test Early Stellar Mass Assembly

We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ 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_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble ACS/F814W, 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\gtrsim10$ sources and 3/15 (20%) likely low-redshift contaminants. Three of our $z\sim12$ candidates push the limits of early stellar mass assembly: they have estimated stellar masses $\sim5\times10^{9}\,M_\odot$, implying an effective stellar baryon fraction of $\epsilon_{\star}\sim0.2-0.5$, where $\epsilon_{\star}\equiv M_{\star}/(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_\star\sim10^{10}\,M_\odot$ galaxies relative to $M_\star\sim10^{9}\,M_\odot$ -- 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 UVLF from a double powerlaw to Schechter at $z\approx8$. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understanding how, and if, such early massive galaxies push the limits of galaxy formation in $\Lambda$CDM.Comment: 30 pages, 9 figures; ApJ submitte

Uncovering a Massive z~7.65 Galaxy Hosting a Heavily Obscured Radio-Loud QSO Candidate in COSMOS-Web

In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud QSO candidate selected using JWST NIRCam/MIRI, mid-IR, sub-mm, and radio imaging in the COSMOS-Web field. Using multi-frequency radio observations and mid-IR photometry, we identify a powerful, radio-loud (RL), growing supermassive black hole (SMBH) with significant spectral steepening of the radio SED ($f_{1.32 \mathrm{GHz}} \sim 2$ mJy, $q_{24\mu m} = -1.1$, $\alpha_{1.32-3\mathrm{GHz}}=-1.2$, $\Delta \alpha = -0.4$). In conjunction with ALMA, deep ground-based observations, ancillary space-based data, and the unprecedented resolution and sensitivity of JWST, we find no evidence of QSO contribution to the UV/optical/NIR data and thus infer heavy amounts of obscuration (N$_{\mathrm{H}} > 10^{23}$ cm$^{-2}$). Using the wealth of deep UV to sub-mm photometric data, we report a singular solution photo-z of $z_\mathrm{phot}$ = 7.65$^{+0.4}_{-0.3}$ and estimate an extremely massive host-galaxy ($\log M_{\star} = 11.92 \pm 0.06\,\mathrm{M}_{\odot}$). This source represents the furthest known obscured RL QSO candidate, and its level of obscuration aligns with the most representative but observationally scarce population of QSOs at these epochs.Comment: Submitted to ApJL, Comments welcom

Unveiling the distant Universe: Characterizing z≥9z\ge9 Galaxies in the first epoch of COSMOS-Web

We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ ($\langle z\rangle=10.0$), UV-magnitudes between M$_{\rm UV}$ = $-$21.2 and $-$19.5 (with $\langle$M$_{\rm UV}\rangle=-20.2$) and rest-frame UV slopes ($\langle \beta\rangle=-2.4$). These galaxies are, on average, more luminous than most $z\ge9$ candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED-fitting are blue ($\beta\sim$[$-$2.0, $-$2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower redshift Universe. The derived stellar masses with $\langle \log_{\rm 10} ($M$_\star/$M$_\odot)\rangle\approx8-9$ are not in tension with the standard $\Lambda$CDM model and our measurement of the volume density of such UV luminous galaxies aligns well with previously measured values presented in the literature at $z\sim9-10$. Our sample of galaxies, although compact, are significantly resolved.Comment: Submitted to Ap

Unveiling the Distant Universe: Characterizing z ≥ 9 Galaxies in the First Epoch of COSMOS-Web

We report the identification of 15 galaxy candidates at z ≥ 9 using the initial COSMOS-Web JWST observations over 77 arcmin2 through four Near Infrared Camera filters (F115W, F150W, F277W, and F444W) with an overlap with the Mid-Infrared Imager (F770W) of 8.7 arcmin2. We fit the sample using several publicly available spectral energy distribution (SED) fitting and photometric redshift codes and determine their redshifts between z = 9.3 and z = 10.9 (〈z〉 = 10.0), UV magnitudes between M UV = −21.2 and −19.5 (with 〈M UV〉 = −20.2), and rest-frame UV slopes (〈β〉 = −2.4). These galaxies are, on average, more luminous than most z ≥ 9 candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED fitting are blue (β ∼ [−2.0, −2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower-redshift Universe. The derived stellar masses with 〈log10( M ⋆/M ⊙)〉 ≈ 8–9 are not in tension with the standard Lambda cold dark matter (ΛCDM) model, and our measurement of the volume density of such UV-luminous galaxies aligns well with previously measured values presented in the literature at z ∼ 9–10. Our sample of galaxies, although compact, is significantly resolved