Searching Far and Long I: Pilot ALMA 2mm Follow-up of Bright Dusty Galaxies as a Redshift Filter

A complete census of dusty star-forming galaxies (DSFGs) at early epochs is necessary to constrain the obscured contribution to the cosmic star formation rate density (CSFRD), however DSFGs beyond $z \sim 4$ are both rare and hard to identify from photometric data alone due to degeneracies in submillimeter photometry with redshift. Here, we present a pilot study obtaining follow-up Atacama Large Millimeter Array (ALMA) $2\,$mm observations of a complete sample of 39 $850\,\rm\mu m$-bright dusty galaxies in the SSA22 field. Empirical modeling suggests $2\,$mm imaging of existing samples of DSFGs selected at $850\,\rm\mu m - 1\,$mm can quickly and easily isolate the "needle in a haystack" DSFGs that sit at $z>4$ or beyond. Combining archival submillimeter imaging with our measured ALMA $2\,$mm photometry ($1\sigma \sim 0.08\,$mJy$\,$beam$^{-1}$ rms), we characterize the galaxies' IR SEDs and use them to constrain redshifts. With available redshift constraints fit via the combination of six submillimeter bands, we identify 6/39 high-$z$ candidates each with $>50\%$ likelihood to sit at $z > 4$, and find a positive correlation between redshift and $2\,$mm flux density. Specifically, our models suggest the addition of $2\,$mm to a moderately constrained IR SED will improve the accuracy of a millimeter-derived redshift from $\Delta z/(1+z) = 0.3$ to $\Delta z/(1+z) = 0.2$. Our IR SED characterizations provide evidence for relatively high emissivity spectral indices ($\langle \beta \rangle = 2.4\pm0.3$) in the sample. We measure that especially bright ($S_{850\rm\mu m}>5.55\,$mJy) DSFGs contribute $\sim10$% to the cosmic-averaged CSFRD from $2<z<5$, confirming findings from previous work with similar samples.Comment: 22 pages, 7 figures, accepted for publication in Ap

A Mixture of LBG Overdensities in the Fields of Three 6<z<76 < z < 7 Quasars: Implications for the Robustness of Photometric Selection

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 202

Missing Giants: Predictions on Dust-Obscured Galaxy Stellar Mass Assembly Throughout Cosmic Time

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 Ap

Emergence of an Ultra-Red Ultra-Massive Galaxy Cluster Core at z=4

Recent simulations and observations of massive galaxy cluster evolution predict that the majority of stellar mass buildup happens within cluster members by z = 2, before cluster virialization. Protoclusters rich with dusty, star-forming galaxies (DSFGs) at z > 3 are the favored candidate progenitors for these massive galaxy clusters at z ~ 0. We present here the first study analyzing stellar emission along with cold dust and gas continuum emission in a spectroscopically confirmed z = 4.002 protocluster core rich with DSFGs, the Distant Red Core (DRC). We combine new Hubble Space Telescope and Spitzer data with existing Gemini, Herschel, and Atacama Large Millimeter/submillimeter Array observations to derive individual galaxy-level properties and compare them to coeval field and other protocluster galaxies. All of the protocluster members are massive (>1010 M ⊙), but not significantly more so than their coeval field counterparts. Within uncertainty, all are nearly indistinguishable from galaxies on the star-forming versus stellar mass main-sequence relationship and the star formation efficiency plane. Assuming no future major influx of fresh gas, we estimate that these gaseous DSFGs will deplete their gas reservoirs in ~300 Myr, becoming the massive quiescent ellipticals dominating cluster cores by z ~ 3. Using various methodologies, we derive a total z = 4 halo mass of ~1014 M ⊙ and estimate that the DRC will evolve to become an ultramassive cluster core of mass 1015 M ⊙ by z = 0

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

Evidence for a Shallow Evolution in the Volume Densities of Massive Galaxies at z=4z=4 to 88 from CEERS

We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam wide-field slitless spectroscopy. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei (AGN). We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for $z \sim$ 1--4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.Comment: 20 pages, 10 figure

JWST CEERS probes the role of stellar mass and morphology in obscuring galaxies

A population of massive galaxies invisible or very faint in deep optical/near-infrared surveys, but brighter at longer wavelengths has been uncovered in the past years. However, the nature of these optically dark/faint galaxies (OFGs, among other nomenclatures) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. Particularly, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on understanding why galaxies like OFGs are so faint at optical/near-infrared wavelengths. We find that stellar mass is the primary proxy of dust attenuation among those studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation 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$), exhibiting 30% smaller effective radius than the average SFGs at the same stellar mass and redshift. On the contrary, they do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show stellar mass as a primary proxy of dust attenuation and compact stellar light profiles behind thick dust columns obscuring typical massive SFGs.Comment: Submitted to A&A. 13 pages, 9 figure

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

CEERS Key Paper. V. Galaxies at 4 &lt; z &lt; 9 Are Bluer than They Appear-Characterizing Galaxy Stellar Populations from Rest-frame ∼1 μm Imaging

We present results from the Cosmic Evolution Early Release Survey on the stellar population parameters for 28 galaxies with redshifts 4 &lt; z &lt; 9 using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 μm data extend the coverage of the rest-frame spectral energy distribution to nearly 1 μm for galaxies in this redshift range. By modeling the galaxies’ SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å)—I-band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower ratios of stellar mass to light. The MIRI data reduce the stellar masses by 〈 Δ log M * 〉 = 0.25 dex at 4 &lt; z &lt; 6 and 0.37 dex at 6 &lt; z &lt; 9. This also reduces the star formation rates (SFRs) by 〈ΔlogSFR〉 = 0.14 dex at 4 &lt; z &lt; 6 and 0.27 dex at 6 &lt; z &lt; 9. The MIRI data also improve constraints on the allowable stellar mass formed in early star formation. We model this using a star formation history that includes both a “burst” at z f = 100 and a slowly varying (“delayed-τ”) model. The MIRI data reduce the allowable stellar mass by 0.6 dex at 4 &lt; z &lt; 6 and by ≈1 dex at 6 &lt; z &lt; 9. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early Universe (z ≈ 9). Therefore, observations of rest-frame ≳1 μm are paramount for constraining the stellar-mass buildup in galaxies at very high redshifts.</p