Search

The dearth of high quality spectroscopy of dusty star-forming galaxies (DSFGs) -- the main drivers of the assembly of dust and stellar mass at the peak of activity in the Universe -- greatly hinders our ability to interpret their physical processes and evolutionary pathways. We present JWST/NIRSpec observations from RUBIES of four submillimeter-selected, ALMA-detected DSFGs at cosmic noon, $z\sim2.3-2.7$. While photometry uniformly suggests vigorous ongoing star formation for the entire sample in line with canonical DSFGs, the spectra differ: one source has spectroscopic evidence of an evolved stellar population, indicating a recent transition to a post-starburst phase, while the remainder show strong spectroscopic signatures of ongoing starbursts. All four galaxies are infrared-luminous (log$_{10}$$L_{\rm{IR}}$/L$_{\rm \odot}$ $>12.4$), massive (log$_{10}\,M_\star$/M$_{\rm \odot}$ $>11$), and very dust-obscured ($A_V\sim3-4$ ABmag). Leveraging detections of multiple Balmer and Paschen lines, we derive an optical attenuation curve consistent with Calzetti overall, yet an optical extinction ratio $R_V\sim2.5$, potentially indicating smaller dust grains or differences in star-dust geometry. This case study provides some of the first detailed spectroscopic evidence that the DSFGs encompass a heterogeneous sample spanning a range of star formation properties and evolutionary stages, and illustrates the advantages of synergistic JWST and ALMA analysis of DSFGs., Comment: 20 pages, 4 figures; submitted to ApJ

JWST has now revealed a population of broad-line AGN at $z>4$ characterized by a distinctive SED shape, with very red rest-frame optical and very blue rest-frame UV continuum. While the optical continuum is thought to originate from the accretion disk, the origin of the UV continuum has been largely unclear. We report the detection of the strong rest-frame UV emission lines of CIII]$\lambda\lambda$1907,1909 and CIV$\lambda\lambda$1549,1551 in a "little red dot" AGN, COS-66964. Spectroscopically confirmed at $z=7.0371$, COS-66964 exhibits broad H$\alpha$ emission (FWHM $\sim 2000$ km s$^{-1}$), and weak broad H$\beta$, implying significant dust attenuation to the BLR ($A_V = 3.9^{+1.7}_{-0.9}$). The H$\alpha$ line width implies a central SMBH mass of $M_{\rm BH} = \left(1.9^{+1.6}_{-0.7}\right)\times10^{7}$ M$_\odot$, and an Eddington ratio $\lambda\sim0.3$-$0.5$. While marginal HeII$\lambda4687$ and [FeX]$\lambda6376$ detections further indicate that the AGN dominates in the rest-frame optical, the non-detection of HeII$\lambda1640$ in the UV despite high EW CIII] and CIV ($\sim 35$ {\AA}) is more consistent with photoionization by massive stars. The non-detection of MgII$\lambda\lambda$2800 is similarly inconsistent with an AGN scattered light interpretation. Assuming the rest-frame UV is dominated by stellar light, we derive a stellar mass of $\log M_\star/M_\odot\sim8.5$, implying an elevated $M_{\rm BH}/M_\star$ ratio $\sim2$ orders of magnitude above the local relation, but consistent with other high-$z$ AGN discovered by JWST. The source is unresolved in all bands, implying a very compact size $\lesssim200$ pc in the UV. This suggests that the simultaneous buildup of compact stellar populations (i.e., galaxy bulges) and the central SMBH is ongoing even at $z>7$., Comment: 14 pages, 6 figures; submitted to ApJL

The current paradigm for the co-evolution of galaxies and their supermassive black holes postulates that dust-obscured active galactic nuclei (AGNs) represent a transitional phase towards a more luminous and unobscured state. However, our understanding of dusty AGNs and their host galaxies at early cosmic times is inadequate due to observational limitations. Here, we present JWST observations of CID-931, an X-ray-detected AGN at a spectroscopic redshift of $z_{\rm spec}=4.91$. Multiband NIRCam imaging from the COSMOS-Web program reveals an unresolved red core, similar to JWST-discovered dusty AGNs. Strikingly, the red core is surrounded by at least eight massive star-forming clumps spread over $1.\!\!^{\prime\prime}6 \approx 10~{\rm kpc}$, each of which has a stellar mass of $10^9-10^{10}M_\odot$ and $\sim0.1-1~{\rm kpc}$ in radius. The whole system amounts to $10^{11}M_\odot$ in stellar mass, higher than typical star-forming galaxies at the same epoch. In this system, gas inflows and/or complex merger events may trigger clump formation and AGN activity thus leading to the rapid formation of a massive galaxy hosting a supermassive black hole. Future follow-up observations will provide new insights into the evolution of the galaxy-black hole relationship during such transitional phases in the early universe., Comment: 13 pages, 10 figures, submitted to PASJ

One of the greatest challenges in galaxy evolution over the last decade has been constraining the prevalence of heavily dust-obscured galaxies in the early Universe. At $z>3$, these galaxies are increasingly rare, and difficult to identify as they are interspersed among the more numerous dust-obscured galaxy population at $z=1-3$, making efforts to secure confident spectroscopic redshifts expensive, and sometimes unsuccessful. In this work, we present the Extended Mapping Obscuration to Reionization with ALMA (Ex-MORA) Survey -- a 2mm blank-field survey in the COSMOS-Web field, and the largest ever ALMA blank-field survey to-date covering 577 arcmin$^2$. Ex-MORA is an expansion of the MORA survey designed to identify primarily $z>3$ dusty, star-forming galaxies while simultaneously filtering out the more numerous $z<3$ population by leveraging the very negative $K$-correction at observed-frame 2mm. We identify 37 significant ($>$5$\sigma$) sources, 33 of which are robust thermal dust emitters. We measure a median redshift of $\langle z \rangle = 3.6^{+0.1}_{-0.2}$, with two-thirds of the sample at $z>3$, and just under half at $z>4$, demonstrating the overall success of the 2mm-selection technique. The integrated $z>3$ volume density of Ex-MORA sources is $\sim1-3\times10^{-5}$ Mpc$^{-3}$, consistent with other surveys of infrared luminous galaxies at similar epochs. We also find that techniques using rest-frame optical emission (or lack thereof) to identify $z>3$ heavily dust-obscured galaxies miss at least half of Ex-MORA galaxies. This supports the idea that the dusty galaxy population is heterogeneous, and that synergies across observatories spanning multiple energy regimes are critical to understanding their formation and evolution at $z>3$., Comment: 13 pages, 5 figures, submitted to ApJ; fully reduced mosaic will be shared upon publication

We present the properties of two candidate massive ($M_\star\sim10^{11}M_\odot$) and dusty ($A_{\rm 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_{\rm spec}=5.051$, while the other has a robust $z_{\rm phot}=6.7\pm0.3$. Thanks to their extremely red colors ($F277W-F444W\sim1.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 SED-fitting tools and star formation histories, we estimate the stellar masses to be $\log(M_\star)=11.32^{+0.07}_{-0.15}$ $M_\odot$ and $\log(M_\star)=11.2^{+0.1}_{-0.2}$ $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 $\epsilon_\star\sim0.25$ and $\epsilon_\star\sim0.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., Comment: 11 pages, 4 figures, 1 table. Accepted for publication in ApJL

JWST has enabled detecting and spatially resolving the heavily dust-attenuated stellar populations of sub-millimeter galaxies, revealing detail that was previously inaccessible. In this work we construct a sample of 289 sub-millimeter galaxies with detailed joint ALMA and JWST constraints in the COSMOS field. Sources are originally selected using the SCUBA-2 instrument and have archival ALMA observations from various programs. Their JWST NIRCam imaging is from COSMOS-Web and PRIMER. We extract multi-wavelength photometry in a manner that leverages the unprecedented near-infrared spatial resolution of JWST, and fit the data with spectral energy distribution models to derive photometric redshifts, stellar masses, star-formation rates and optical attenuation. The sample has an average z=2.6, A_V=2.5, SFR=270 and log(M*)=11.1. There are 81 (30%) galaxies that have no previous optical/near-infrared detections, including 75% of the z>4 sub-sample (n=28). The faintest observed near-infrared sources have the highest redshifts and largest A_V=4. In a preliminary morphology analysis we find that ~10% of our sample exhibit spiral arms and 5% host stellar bars, with one candidate bar found at z>3. Finally, we find that the clustering of JWST galaxies within 10 arcseconds of a sub-mm galaxy is a factor of 2 greater than what is expected based on either random clustering or the distribution of sources around any red galaxy irrespective of a sub-mm detection., Comment: 37 pages (15 for RGBs + references), 14 figures, submitted to ApJ

The James Webb Space Telescope (JWST) has uncovered a ubiquitous population of dust-obscured compact sources at $z\gtrsim 4$. Many of these objects exhibit signs of active galactic nucleus (AGN) activity, making their study crucial for understanding the formation of supermassive black holes (SMBHs) and their growth with host galaxies. In this work, we examine low and medium resolution JWST/NIRSpec spectra from the JADES GTO public data release in the GOODS-N field of a red, luminous ($M_{\rm B}\sim-22.2$ mag) and compact ($<500$ pc) source at $z=4.13$. The rest-optical ($\lambda_{\rm rest} > 4000$ A) continuum of this source is strongly dominated by a massive (log$_{10}$[$M_*/M_\odot] \sim 10.6$), quenched (log$_{10}$[sSFR/yr$^{-1}$] $< -11$) galaxy, as indicated by the clear presence of a Balmer break and stellar absorption lines. Star-formation history modeling reveals a starburst episode followed by rapid quenching about 200 Myr ago. The spectrum shows extremely broad (FWHM $\sim 2500$ km/s) H$\alpha$ emission and elevated optical line ratios, indicating an actively accreting SMBH. Moreover, our work has potentially revealed clear AGN signatures in the rest-UV in LRDs for the first time, via a detection of a strong Ly$\alpha$ emission and a broad MgII, doublet. The derived black hole mass of log$_{10}$($M_{\rm BH}/M_\odot) \sim 7.3$ results in $M_{\rm BH}/M_*\sim 0.04$ %, consistent with the local relations, unlike the elevated ratios in other high-$z$ reddened AGN. Finally, we use JWST data from AGN at $z=4-10$ to explore an evolutionary link between high-$z$ reddened AGN, early quiescent galaxies, and local ellipticals., Comment: 20 pages, 7 figures, 2 tables. Submitted to ApJ

JWST has revealed a ubiquitous population of ``little red dots'' (LRDs) at $z\gtrsim4$, selected via their red rest-frame optical emission and compact morphologies. They are thought to be reddened by dust, whether in tori of active galactic nuclei or the interstellar medium (ISM), though none have direct dust detections to date. Informed by the average characteristics of 675 LRDs drawn from the literature, we provide ballpark constraints on the dust characteristics of the LRD population and estimate they have average dust masses of $\langle M_{\rm dust}\rangle=(1.6^{+4.8}_{-0.9})\times10^{4} M_\odot$, luminosities of $\langle L_{\rm IR}\rangle=(8^{+3}_{-5})\times10^{10} L_\odot$ and temperatures of $\langle T_{\rm dust}\rangle=110^{+21}_{-36}$ K. Notably, the spectral energy distributions are thought to peak at $\sim$100 K (rest-frame 20-30 $\mu$m) regardless of heating mechanism, whether AGN or star formation. LRDs' compact sizes $R_{\rm eff}\sim100$ pc are the dominant factor contributing to their low estimated dust masses. Our predictions likely mean LRDs have, on average, submillimeter emission a factor of $\sim$100$\times$ fainter than current ALMA limits provide. The star-to-dust ratio is a factor $\sim$100$\times$ larger than expected from dust formation models if one assumes the rest-optical light is dominated by stars; this suggests stars do not dominate. Despite their high apparent volume density, LRDs contribute negligibly (0.1%) to the cosmic dust budget at $z\gtrsim4$ due to their low dust masses., Comment: 11 pages, 3 figures; accepted to ApJL

JWST has revealed a population of compact and extremely red galaxies at $z>4$, which likely host active galactic nuclei (AGN). We present a sample of 434 ``little red dots'' (LRDs), selected from the 0.54 deg$^2$ COSMOS-Web survey. We fit galaxy and AGN SED models to derive redshifts and physical properties; the sample spans $z\sim5$-$9$ after removing brown dwarf contaminants. We consider two extreme physical scenarios: either LRDs are all AGN, and their continuum emission is dominated by the accretion disk, or they are all compact star-forming galaxies, and their continuum is dominated by stars. If LRDs are AGN-dominated, our sample exhibits bolometric luminosities $\sim10^{45-47}$ erg\,s$^{-1}$, spanning the gap between JWST AGN in the literature and bright, rare quasars. We derive a bolometric luminosity function (LF) $\sim100$ times the (UV-selected) quasar LF, implying a non-evolving black hole accretion density of $\sim10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ from $z\sim2$-$9$. By contrast, if LRDs are dominated by star formation, we derive stellar masses $\sim10^{8.5-10}\,M_\odot$. MIRI/F770W is key to deriving accurate stellar masses; without it, we derive a mass function inconsistent with $\Lambda$CDM. The median stellar mass profile is broadly consistent with the maximal stellar mass surface densities seen in the nearby universe, though the most massive $\sim50$\% of objects exceed this limit, requiring substantial AGN contribution to the continuum. Nevertheless, stacking all available X-ray, mid-IR, far-IR/sub-mm, and radio data yields non-detections. Whether dominated by dusty AGN, compact star-formation, or both, the high masses/luminosities and remarkable abundance of LRDs implies a dominant mode of early galaxy/SMBH growth., Comment: 27 pages, 13 figures. Submitted to ApJ, comments welcome! Data access at https://github.com/hollisakins/akins24_cw

We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN., Comment: 23 pages, 17 figures, submitted to ApJ

GHZ2/GLASS-z12 has been recently observed by JWST with both NIRSpec and MIRI spectrographs, making it the most distant galaxy ($z_{spec}=12.34$) with complete spectroscopic coverage from rest-frame UV to optical. It is identified as a strong CIV$_{1549}$ emitter with many detected emission lines (NIV], HeII, OIII], NIII], CIII], [OII], [NeIII], [OIII], and H$\alpha$), including a remarkable OIII$_{1333}$ Bowen fluorescence line. We analyze in this paper the joint NIRSpec+MIRI spectral data set. Combining six optical diagnostics (R2, R3, R23, O32, Ne3O2, and Ne3O2Hd), we find extreme ionization conditions, with O32 $=1.39 \pm 0.19$ and Ne3O2 $=0.37 \pm 0.18$ in stark excess compared to typical values in the ISM at lower redshifts. These line properties are compatible either with an AGN or with a compact, dense star-forming environment ($\Sigma_{\rm SFR}$ $\sim 10^2$-$10^3$ Msun/yr/kpc$^2$), with a high ionization parameter ($\log_{10}$(U) $=-1.75 \pm 0.16$), a high ionizing photon production efficiency $\log(\xi_{\rm ion}) = 25.7_{-0.1}^{+0.2}$, and a low, although not pristine, metal content between $5\%$ and $11\%$ Z$_\odot$ (confirmed by the T$_e$ method), indicating a rapid metal enrichment in the last few Myrs. These properties also suggest that a substantial amount of ionizing photons ($\sim 10\%$) are leaking outside. The general lessons learned from GHZ2 are the following: (i) the UV to optical combined nebular indicators are broadly in agreement with UV-only or optical-only indicators. (ii) UV+optical diagnostics fail to discriminate between an AGN and star-formation in a low metallicity, high density, and extreme ionization environment. (iii) comparing the nebular line ratios with local analogs may be approaching its limits at $z \gtrsim 10$, as this approach is potentially challenged by the unique conditions of star formation experienced by galaxies at these extreme redshifts., Comment: 14 pages, 3 figures, and 1 table. Resubmitted to ApJ. A section (3.4) has been added to the previous version

The James Webb Space Telescope (JWST) has discovered a surprising population of bright galaxies in the very early universe (< 500 Myrs after the Big Bang) that is hard to explain with conventional galaxy formation models and whose physical properties remain to be fully understood. Insight into the internal physics of galaxies is captured best via observations of excited-state atomic transitions of ionized gas, but beyond z~7-9, the brightest spectral signatures are redshifted into the mid-infrared regime, where observations are increasingly more difficult. Here, we present the first detection of a hydrogen recombination line (H{\alpha}) and doubly-ionized oxygen ([OIII]4959,5007{\AA}) at z>10 using the JWST Mid-Infrared Instrument, MIRI. These detections place the bright galaxy GHZ2/GLASS-z12 at z=12.33+/-0.02, making it the most distant astronomical object with direct spectroscopic detection of these lines and the brightest confirmed object at this epoch. These observations provide key insights into the conditions of this primeval galaxy, which shows hard ionizing conditions rarely seen in the local Universe and likely driven by compact, young (<30 Myr) star formation. Its oxygen-to-hydrogen abundance is close to a tenth of the solar value, indicating a rapid metal enrichment during the earliest phases of galaxy formation. This study confirms the unique conditions of the brightest and most distant galaxies recently discovered by JWST and the huge potential of mid-IR observations to characterize these systems, opening a range of new possibilities in the study of the very early Universe., Comment: 14 pages, 6 figures. Submitted for publication. Comments are most welcome

We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalactic surveys, focusing on COSMOS-Web, as well as JADES, UNCOVER, CEERS, and PRIMER. Our search strategy includes two stages: applying stringent photometric cuts to catalog-level data and detailed spectral energy distribution fitting. These techniques effectively isolate the quasar candidates while mitigating contamination from low-redshift interlopers, such as brown dwarfs and nearby galaxies. The selected candidates indicate physical traits compatible with low-luminosity active galactic nuclei, likely hosting $\approx10^5$--$10^7~M_\odot$ supermassive black holes (SMBHs) living in galaxies with stellar masses of $\approx10^8$--$10^{10}~M_\odot$. The SMBHs selected in this study, on average, exhibit an elevated mass compared to their hosts, with the mass ratio distribution slightly higher than those of galaxies in the local Universe. As with other high-$z$ studies, this is at least in part due to the selection method for these quasars. An extensive Monte Carlo analysis provides compelling evidence that heavy black hole seeds from the direct collapse scenario appear to be the preferred pathway to mature this specific subset of SMBHs by $z\approx7$. This work underscores the significance of further spectroscopic observations, as the quasar candidates presented here offer exceptional opportunities to delve into the nature of the earliest galaxies and SMBHs that formed during cosmic infancy., Comment: Accepted for publication in the Astronomy & Astrophysics journal. 19 pages, 10 figures, and 4 tables. We welcome comments from the reader

With their high angular resolutions of 30--100 mas, large fields of view, and complex optical systems, imagers on next-generation optical/near-infrared space observatories, such as the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST), present both new opportunities for science and also new challenges for empirical point spread function (PSF) characterization. In this context, we introduce ShOpt, a new PSF fitting tool developed in Julia and designed to bridge the advanced features of PIFF (PSFs in the Full Field of View) with the computational efficiency of PSFEx (PSF Extractor). Along with ShOpt, we propose a suite of non-parametric statistics suitable for evaluating PSF fit quality in space-based imaging. Our study benchmarks ShOpt against the established PSF fitters PSFEx and PIFF using real and simulated COSMOS-Web Survey imaging. We assess their respective PSF model fidelity with our proposed diagnostic statistics and investigate their computational efficiencies, focusing on their processing speed relative to the complexity and size of the PSF models. We find that ShOpt can already achieve PSF model fidelity comparable to PSFEx and PIFF while maintaining competitive processing speeds, constructing PSF models for large NIRCam mosaics within minutes., Comment: 53 pages, 27 figures, submitted to Astronomical Journal. Same as v2 with a corrected typo, which contained revisions from v1

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 $\sim 210 \, k$ galaxies at $0.42$., Comment: 22 pages, 11 figures, 2 tables, Submitted to ApJ

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the first set of observations (Epoch 1, 50$\%$ of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3$\sigma$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star-formation rates, and metallicity ([OIII]/H$\beta$), and dust content, of galaxies at $1

We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO J=5--4 emission with ALMA, PACS-819 is rotation-dominated thus has a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components, not discernable in the previous HST images. The CO J = 5 -- 4 and FIR dust continuum emission are co-spatial with a heavily-obscured starbursting core (<1 kpc) which is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally-distorted dwarf galaxy, and a clump, either a sign of an ongoing violent disk instability or a recently accreted low-mass satellite. With spatially-resolved maps, we find a high molecular gas fraction in the central area reaching $\sim3$ ($M_{\text{gas}}$/$M_*$) and short depletion times ($M_{\text{gas}}/SFR\sim$ 120 Myrs) across the entire system. These observations provide insights into the complex nature of starbursts in the distant universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST., Comment: 21 pages(including appendix), 13 figures, accepted by ApJ

We present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z~9, 11 and 14, and show that the space density of bright (M_UV=-20) galaxies changes only modestly from z~14 to z~9, compared to a steeper increase from z~8 to z~4. While our candidates are photometrically selected, spectroscopic followup has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z>10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models which naturally have enhanced star-formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic followup of a large sample of early galaxies can distinguish between these competing scenarios., Comment: Submitted to ApJL. Main paper is 33 pages, 13 figures, 5 tables. Two appendices with additional figures and tables

We present the results of our Keck/DEIMOS spectroscopic follow-up of candidate galaxies of i-band-dropout protocluster candidate galaxies at $z\sim6$ in the COSMOS field. We securely detect Lyman-$\alpha$ emission lines in 14 of the 30 objects targeted, 10 of them being at $z=6$ with a signal-to-noise ratio of $5-20$, the remaining galaxies are either non-detections or interlopers with redshift too different from $z=6$ to be part of the protocluster. The 10 galaxies at $z\approx6$ make the protocluster one of the riches at $z>5$. The emission lines exhibit asymmetric profiles with high skewness values ranging from 2.87 to 31.75, with a median of 7.37. This asymmetry is consistent with them being Ly$\alpha$, resulting in a redshift range of $z=5.85-6.08$. Using the spectroscopic redshifts, we re-calculate the overdensity map for the COSMOS field and find the galaxies to be in a significant overdensity at the $4\sigma$ level, with a peak overdensity of $\delta=11.8$ (compared to the previous value of $\delta=9.2$). The protocluster galaxies have stellar masses derived from Bagpipes SED fits of $10^{8.29}-10^{10.28} \rm \,M_{\rm \odot}$ and star formation rates of $2-39\,\rm M_{\rm \odot}\rm\,yr^{-1}$, placing them on the main sequence at this epoch. Using a stellar-to-halo-mass relationship, we estimate the dark matter halo mass of the most massive halo in the protocluster to be $\sim 10^{12}\rm M_{\rm \odot}$. By comparison with halo mass evolution tracks from simulations, the protocluster is expected to evolve into a Virgo- or Coma-like cluster in the present day., Comment: 26 pages, 14 figues, 5 tables, main text is 16 pages, appendix is 10 pages, published in MNRAS

This paper presents the results of 475 hours of interferometric observations with the Australia Telescope Compact Array towards the Spiderweb protocluster at \(z=2.16\). We search for large, extended molecular gas reservoirs among 46 previously detected CO(1-0) emitters, employing a customised method we developed. Based on the CO emission images and position-velocity diagrams, as well as the ranking of sources using a binary weighting of six different criteria, we have identified 14 robust and 7 tentative candidates that exhibit large extended molecular gas reservoirs. These extended reservoirs are defined as having sizes greater than 40 kpc or super-galactic scale. This result suggests a high frequency of extended gas reservoirs, comprising at least \(30 \%\) of our CO-selected sample. An environmental study of the candidates is carried out based on N-th nearest neighbour and we find that the large molecular gas reservoirs tend to exist in denser regions. The spatial distribution of our candidates is mainly centred on the core region of the Spiderweb protocluster. The performance and adaptability of our method are discussed. We found 13 (potentially) extended gas reservoirs located in nine galaxy (proto)clusters from the literature. We noticed that large extended molecular gas reservoirs surrounding (normal) star-forming galaxies in protoclusters are rare. This may be attributable to the lack of observations low-J CO transitions and the lack of quantitative analyses of molecular gas morphologies. The large gas reservoirs in the Spiderweb protocluster are a potential source of the intracluster medium seen in low redshift Virgo- or Coma-like galaxy clusters., Comment: 25 pages, 16 figures. Accepted for publication in MNRAS

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 submitted

Low mass (sub)stellar objects represent the low end of the initial mass function, the transition to free-floating planets and a prominent interloper population in the search for high-redshift galaxies. Without proper motions or spectroscopy, can one identify these objects photometrically? JWST/NIRCam has several advantages over HST/WFC3 NIR: more filters, a greater wavelength range, and greater spatial resolution. Here, we present a catalogue of (sub)stellar dwarfs identified in the Cosmic Evolution Early Release Science Survey (CEERS). We identify 518 stellar objects down to $m_F200W \sim 28$ using half-light radius, a full three magnitudes deeper than typical HST/WFC3 images. A kNN nearest neighbour algorithm identifies and types these sources, using four HST/WFC3 and four NIRCam filters, trained on SpeX spectra of nearby brown dwarfs. The kNN with four neighbors classifies well within two subtypes: e.g M2$\pm$2 or T4$\pm$2, achieving $\sim$95% precision and recall. More granular typing results in worse metrics. In CEERS, we find 9 M8$\pm$2, 2 L6$\pm$2, 1 T4$\pm$2, and 15 T8$\pm$2. We compare the observed long wavelength NIRCam colours -- not used in the kNN -- to those expected for brown dwarf atmospheric models. The NIRCam F356W-F444W and F410M-F444W colours are redder by a magnitude for the type assigned by the kNN, hinting at a wider variety of atmospheres for these objects. We find a 300-350pc scale-height for M6$\pm$2 dwarfs plus a second structural component and a 150-200pc scale-height for T6$\pm$2 type dwarfs, consistent with literature values., Comment: 18 pages, 24 figures, 3 tables, submitted to MNRAS

In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud AGN 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.28 \mathrm{GHz}} \sim 2$ mJy, $q_{24\mu m} = -1.1$, $\alpha_{1.28-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 AGN 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.7$^{+0.4}_{-0.3}$ and estimate an extremely massive host-galaxy ($\log M_{\star} = 11.4 -12\,\mathrm{M}_{\odot}$) hosting a powerful, growing SMBH (L$_{\mathrm{Bol}} = 4-12 \times 10^{46}$ erg s$^{-1}$). This source represents the furthest known obscured RL AGN candidate, and its level of obscuration aligns with the most representative but observationally scarce population of AGN at these epochs., Comment: Accepted to ApJL

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 submitted

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 ApJ

We present a spectroscopic survey of Ly$\alpha$ emitters in the Extended Groth Strip (EGS) field, targeting the regime near the Epoch of Reionization. Using Keck/DEIMOS, we observed 947 high-$z$ candidates with photometric redshifts from 3 $< z_\text{phot} <$ 7 and down to an $H$-band (HST/WFC3 F160W) magnitude limit of < 27.5. Observations were taken over the course of 8 nights, with integration times ranging from 4 to 7.8 hours. Our survey secured 137 unique redshifts, 126 of which are Ly$\alpha$ emitters at 2.8 $< z <$ 6.5 with a mean redshift of $\overline{z} = 4.3$. We provide a comprehensive redshift catalog for our targets, as well as the reduced one- and two- dimensional spectra for each object. These observations will provide an important auxiliary dataset for the JWST Directors Discretionary Early Release Science (DD-ERS) program the Cosmic Evolution Early Release Science Survey (CEERS), which recently completed near- and mid-IR imaging and spectroscopy of galaxies in the EGS field., Comment: 8 pages, 7 figures, 2 tables, Accepted in MNRAS

We present a sample of 30 massive (log$(M_{\ast}/M_\odot) >11$) $z=3-5$ quiescent galaxies selected from the \textit{Spitzer-}HETDEX Exploratory Large Area (SHELA) Survey and observed at 1.1mm with Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations. These ALMA observations would detect even modest levels of dust-obscured star-formation, on order of $\sim 20 \ M_\odot \textrm{yr}^{-1}$ at $z\sim4$ at a $1\sigma$ level, allowing us to quantify the amount of contamination from dusty star-forming sources in our quiescent sample. Starting with a parent sample of candidate massive quiescent galaxies from the Stevans et al. 2021 v1 SHELA catalog, we use the Bayesian \textsc{Bagpipes} spectral energy distribution fitting code to derive robust stellar masses ($M_*$) and star-formation rates (SFRs) for these sources, and select a conservative sample of 36 candidate massive ($M_* > 10^{11}M_\odot$) quiescent galaxies, with specific SFRs at $>2\sigma$ below the star-forming main sequence at $z\sim4$. Based on ALMA imaging, six of these candidate quiescent galaxies have the presence of significant dust-obscured star-formation, thus were removed from our final sample. This implies a $\sim 17\%$ contamination rate from dusty star-forming galaxies with our selection criteria using the v1 SHELA catalog. This conservatively-selected quiescent galaxy sample at $z=3-5$ will provide excellent targets for future observations to better constrain how massive galaxies can both grow and shut-down their star-formation in a relatively short time period., Comment: 16 pages, 13 figures

Substantial populations of massive quiescent galaxies at $z\ge3$ challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at $z>3$. As a result, studies report that as much as 70\% of quiescent galaxies at $z>3$ may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at $3\lesssim z \lesssim 6$ using JWST NIRCam imaging data. We use empirically-constrained galaxy SED templates to define a region in the $F277W-F444W$ vs. $F150W-F277W$ color plane that captures quiescent galaxies at $z>3$. We apply this color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and identify 44 candidate $z\gtrsim3$ quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of post-starburst galaxies. We derive volume density estimates of $n\sim1-4\times10^{-5}$\,Mpc$^{-3}$ at $3< z <5$, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam's wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies., Comment: 12 pages, 4 figures, 2 tables, submitted to ApJ

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{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 ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED 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_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

The most luminous quasars at $z > 6$ are suspected to be both highly clustered and reside in the most massive dark matter halos in the early Universe, making them prime targets to search for galaxy overdensities and/or protoclusters. We search for Lyman-break dropout-selected galaxies using HST WFC3/ACS broadband imaging in the fields of three $6 < z < 7$ quasars, as well as their simultaneously observed coordinated-parallel fields, and constrain their photometric redshifts using EAZY. One field, J0305-3150, shows a volume density 10$\times$ higher than the blank-field UV luminosity function (UVLF) at M$_{UV} < -20$, with tentative evidence of a 3$\sigma$ overdensity in its parallel field located 15 cMpc away. Another field, J2054-0005, shows an angular overdensity within 500 ckpc from the quasar but still consistent with UVLF predictions within 3$\sigma$, while the last field, J2348-3054, shows no enhancement. We discuss methods for reducing uncertainty in overdensity measurements when using photometric selection and show that we can robustly select LBGs consistent with being physically associated with the quasar, corroborated by existing JWST/NIRCam WFSS data in the J0305 field. Even accounting for incompleteness, the overdensities in J0305 and J2054 are higher for brighter galaxies at short angular separations, suggesting preferential enhancement of more massive galaxies in the immediate vicinity of the quasar. Finally, we compare the LBG population with previously-identified [CII] and mm-continuum companions; the LBG overdensities are not accompanied by an enhanced number of dusty galaxies, suggesting that the overdense quasar fields are not in the bursty star-forming phase sometimes seen in high-redshift protoclusters., Comment: 22 pages (main text), 12 figures, 10 tables, 2 appendices. Final version after addressing referee report, accepted to ApJ May 2023

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

A growing number of far-infrared 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/NIRCam counterpart from the COSMOS-Web survey to a far-infrared SCUBA-2 and ALMA source, AzTECC71, which was previously undetected at wavelengths shorter than 850 microns. AzTECC71, amongst the reddest galaxies in COSMOS-Web with F277W - F444W~0.9, is undetected in NIRCam/F150W and F115W and fainter in F444W than other sub-millimeter 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 far-IR including ALMA Band 6, we find a high probability (99%) that AzTECC71 is at z>4 with z_phot=5.7(+0.8,-0.7). This galaxy is massive (logM*/Msun~10.7) and IR-luminous (logLIR/Lsun~12.7), comparable to other optically-undetected but far-IR 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 far-IR-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-10x greater than previously estimated., Comment: 13 pages, 7 figures, submitted to ApJ

We analyze rest-frame ultraviolet to optical spectra of three $z\simeq7.47$ - $7.75$ galaxies whose Ly$\alpha$-emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) survey. From NIRSpec data, we confirm the systemic redshifts of these Ly$\alpha$ emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal poor. We investigate Ly$\alpha$ line properties, including the line flux, velocity offset, and spatial extension. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements ($\sim5\times$ greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly$\alpha$ emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly$\alpha$ in one galaxy is significantly more extended than the non-resonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly$\alpha$ sources, discussing viable scenarios for the creation of sizable ionized bubbles ($>$1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, requiring additional sources of ionizing photons. Therefore, the fact that we detect Ly$\alpha$ from these galaxies suggests diverse scenarios on escape of Ly$\alpha$ during the epoch of reionization. High spectral resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization., Comment: Submitted to ApJ (18 pages, 7 figures, 2 tables)

We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 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 square degree (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 Lyman-alpha rest-frame equivalent widths (EW). We uncover marginal (1 sigma significance) correlations between stellar mass and Lyman-alpha EW, and star formation rate (SFR) and Lyman-alpha EW, with a Spearman correlation coefficient of -0.$34_{-.14}^{+.17}$ and -0.$37_{-.14}^{+.16}$ respectively. We show that the Lyman-alpha distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of 150 Angstrom. Once complete the TESLA survey will enable the study of ~ thousands of LAEs to explore correlations between galaxy properties and Lyman-alpha EW. The large sample size will allow the construction of a predictive model for the Lyman-alpha EW as a function of SED-derived galaxy properties, which could be used to improve Lyman-alpha based constraints on reionization.

During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far beyond pre-JWST limits. While generally robust, such photometric redshifts can suffer from degeneracies and occasionally catastrophic errors. Spectroscopic measurement is required to validate these sources and to reliably quantify physical properties that can constrain galaxy formation models and cosmology. Here we present JWST spectroscopy that confirms redshifts for two very luminous galaxies with $z > 11$, but also demonstrates that another candidate with suggested $z\approx 16$ instead has $z = 4.9$, with an unusual combination of nebular line emission and dust reddening that mimics the colors expected for much more distant objects. These results reinforce evidence for the early, rapid formation of remarkably luminous galaxies, while also highlighting the necessity of spectroscopic verification. The large abundance of bright, early galaxies may indicate shortcomings in current galaxy formation models, or deviation from physical properties (such as the stellar initial mass function) that are generally believed to hold at later times., Comment: Accelerated Article Preview (AAP) version published in Nature

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $\mu$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch., Comment: 26 pages, 14 figures, 6 tables, published in ApJL

We present CEERS JWST/NIRCam imaging of a massive galaxy group at z=1.85, to explore the early JWST view on massive group formation in the distant Universe. The group contains >16 members (including 6 spectros. confirmations) down to log10(Mstar/Msun)=8.5, including the brightest group galaxy (BGG) in the process of actively assembling at this redshift. The BGG is comprised of multiple merging components extending ~3.6" (30kpc) across the sky. The BGG contributes 69% of the group's total galactic stellar mass, with one of the merging components containing 76% of the total mass of the BGG and a SFR>1810Msun/yr. Most importantly, we detect intra-halo light (IHL) in several HST and JWST/NIRCam bands, allowing us to construct a state-of-the-art rest-frame UV-NIR Spectral Energy Distribution of the IHL for the first time at this high redshift. This allows stellar population characterisation of both the IHL and member galaxies, as well as the morphology distribution of group galaxies vs. their star-formation activity when coupled with Herschel data. We create a stacked image of the IHL, giving us a sensitivity to extended emission of 28.5 mag/arcsec2 at rest-frame 1um. We find that the IHL is extremely dust poor (Av~0), containing an evolved stellar population of log10(t50/yr)=8.8, corresponding to a formation epoch for 50% of the stellar material 0.63Gyr before z=1.85. There is no evidence of ongoing star-formation in the IHL. The IHL in this group at z=1.85 contributes ~10% of the total stellar mass, comparable with what is observed in local clusters. This suggests that the evolution of the IHL fraction is more self-similar with redshift than predicted by some models, challenging our understanding of IHL formation during the assembly of high-redshift clusters. JWST is unveiling a new side of group formation at this redshift, which will evolve into Virgo-like structures in the local Universe., Comment: 14 pages + appendix, 11 figures, 4 tables. Accepted to A&A on 15th May 2023

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.512, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products., Comment: 21 pages, 9 figures, submitted to AAS Journals

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5\mu$m either with the NIRSpec prism or its three medium resolution gratings. For $z\simeq9$ photometric candidates, we achieve a high confirmation rate of $\simeq$90\%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at $z>10$, where the strong [OIII] and H$\beta$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$\alpha$ continuum break is not detected with the current sensitivity. Compared with HST-selected bright galaxies ($M_{\rm\,UV}\simeq-22$) that are similarly spectroscopically confirmed at $z\gtrsim8$, these NIRCam-selected galaxies are characterized by lower star formation rates (SFR$\simeq4\,M_{\odot}$~yr$^{-1}$) and lower stellar masses ($\simeq10^{8}\,M_{\odot}$), but with higher [OIII]+H$\beta$ equivalent widths ($\simeq$1100$\r{A}$), and elevated production efficiency of ionizing photons ($\log(\xi_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($<10$~Myrs) accounting for $\simeq20\%$ of the galaxy mass, highlighting the key contribution of faint galaxies to cosmic reionization. Taking advantage of the homogeneous selection and sensitivity, we also investigate metallicity and ISM conditions with empirical calibrations using the [OIII]/H$\beta$ ratio. We find that galaxies at $z\sim8-9$ have higher SFRs and lower metallicities than galaxies at similar stellar masses at $z\sim2-6$, which is generally consistent with the current galaxy formation and evolution models., Comment: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour 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$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. 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 ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.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 sub-dwarf 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., Comment: 46 pages, 16 figures, ApJ accepted

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

We present ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$\mu$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr$^{-1}$. We detect a 5.1$\sigma$ line feature at $338.726\pm0.007$~GHz exactly coinciding with the JWST source position, with a 2\% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$\mu$m at $z=16.01$ or [CII]158$\mu$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11--13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11--17 is favored in every candidate as a very blue (UV continuum slope of $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several candidates that reasonable SED fits ($\Delta$ $\chi^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$--5, where such populations predicted from their luminosity functions and EW([OIII]+H$\beta$) distributions are abundant in survey volumes used for the identification of the $z\sim$11--17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search., Comment: 25 pages, 8 figures, 7 tables. ApJ in press

Due to their extremely dust-obscured nature, much uncertainty still exists surrounding the stellar mass growth and content in dusty, star-forming galaxies (DSFGs) at $z>1$. In this work, we present a numerical model built using empirical data on DSFGs to estimate their stellar mass contributions across the first $\sim$10 Gyr of cosmic time. We generate a dust-obscured stellar mass function that extends beyond the mass limit of star-forming stellar mass functions in the literature, and predict that massive DSFGs constitute as much as $50-100\%$ of all star-forming galaxies with M $\ge10^{11}$M$_\odot$ at $z>1$. We predict the number density of massive DSFGs and find general agreement with observations, although more data is needed to narrow wide observational uncertainties. We forward model mock massive DSFGs to their quiescent descendants and find remarkable agreement with observations from the literature demonstrating that, to first order, massive DSFGs are a sufficient ancestral population to describe the prevalence of massive quiescent galaxies at $z>1$. We predict that massive DSFGs and their descendants contribute as much as $25-60\%$ to the cosmic stellar mass density during the peak of cosmic star formation, and predict an intense epoch of population growth during the $\sim1$ Gyr from $z=6$ to 3 during which the majority of the most massive galaxies at high-$z$ grow and then quench. Future studies seeking to understand massive galaxy growth and evolution in the early Universe should strategize synergies with data from the latest observatories (e.g. JWST and ALMA) to better include the heavily dust-obscured galaxy population., Comment: 22 pages, 9 figures, submitted to ApJ

We conduct a systematic search for protocluster candidates at $z \geq 6$ in the 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\sigma$) candidate galaxy overdensities across the redshift range $6\le z\le7.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 $\sim10^{11-13}\,M_{\rm \odot}$, 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 $\Lambda$CDM cosmology. Through comparison with simulation, we expect that all the overdensities at $z\simeq6$ will evolve into a Virgo-/Coma-like clusters at present (i.e., with masses $\sim 10^{14}-10^{15}\,M_{\rm \odot}$). Compared to other overdensities identified at $z \geq 6$ via narrow-band selection techniques, the overdensities presented appear to have $\sim10\times$ higher stellar masses and star-formation rates. We compare the evolution in the total star-formation rate and stellar mass content of the protocluster candidates across the redshift range $6\le z\le7.7$ and find agreement with the total average star-formation rate from simulations., Comment: 52 pages, 32 figues, 18 tables, main text is 30 pages, appendix is 22 pages, to be published in ApJ