Showing 1–50 of 139 results for author: Stefanon, M

Carbon and Iron Deficiencies in Quiescent Galaxies at z=1-3 from JWST-SUSPENSE: Implications for the Formation Histories of Massive Galaxies

Authors: Aliza G. Beverage, Martje Slob, Mariska Kriek, Charlie Conroy, Guillermo Barro, Rachel Bezanson, Gabriel Brammer, Chloe M. Cheng, Anna de Graaff, Natascha M. Förster Schreiber, Marijn Franx, Brian Lorenz, Pavel E. Mancera Piña, Danilo Marchesini, Adam Muzzin, Andrew B. Newman, Sedona H. Price, Alice E. Shapley, Mauro Stefanon, Katherine A. Suess, Pieter van Dokkum, David Weinberg, Daniel R. Weisz

Abstract: We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log M/M$_\odot$=10.2-11.2) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 dex in [Mg/H], implying rapid formation and quenching befor… ▽ More We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log M/M$_\odot$=10.2-11.2) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 dex in [Mg/H], implying rapid formation and quenching before significant enrichment from asymptotic giant branch stars and Type Ia supernovae. Additionally, we find that galaxies that form at higher redshift have higher [Mg/Fe] and lower [Fe/H] and [Mg/H], irrespective of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily explained by lower redshift samples naturally including galaxies with longer star-formation timescales. On the other hand, the lower [Mg/H] can be explained by galaxies forming at earlier epochs expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is also lower at z=1-3 compared to the lower redshift relation, though the slopes are similar. Finally, we compare our results to standard stellar population modeling approaches employing solar abundance patterns and non-parametric star-formation histories (using Prospector). Our SSP-equivalent ages agree with the mass-weighted ages from Prospector, while the metallicities disagree significantly. Nonetheless, the metallicities better reflect [Fe/H] than total [Z/H]. We also find that star-formation timescales inferred from elemental abundances are significantly shorter than those from Prospector, and we discuss the resulting implications for the early formation of massive galaxies. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ; 18 pages, 6 figures, 1 table

Bulge+disc decomposition of HFF and CANDELS galaxies: UVJ diagrams and stellar mass-size relations of galaxy components at $0.2 \leq z \leq 1.5$

Authors: Kalina V. Nedkova, Boris Häußler, Danilo Marchesini, Gabriel B. Brammer, Adina D. Feinstein, Evelyn J. Johnston, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Nicholas S. Martis, Adam Muzzin, Marc Rafelski, Heath V. Shipley, Rosalind E. Skelton, Mauro Stefanon, Arjen van der Wel, Katherine E. Whitaker

Abstract: Using deep imaging from the CANDELS and HFF surveys, we present bulge+disc decompositions with GalfitM for $\sim$17,000 galaxies over $0.2 \leq z\leq 1.5$. We use various model parameters to select reliable samples of discs and bulges, and derive their stellar masses using an empirically calibrated relation between mass-to-light ratio and colour. Across our entire redshift range, we show that disc… ▽ More Using deep imaging from the CANDELS and HFF surveys, we present bulge+disc decompositions with GalfitM for $\sim$17,000 galaxies over $0.2 \leq z\leq 1.5$. We use various model parameters to select reliable samples of discs and bulges, and derive their stellar masses using an empirically calibrated relation between mass-to-light ratio and colour. Across our entire redshift range, we show that discs follow stellar mass-size relations that are consistent with those of star-forming galaxies, suggesting that discs primarily evolve via star formation. In contrast, the stellar mass-size relations of bulges are mass-independent. Our novel dataset further enables us to separate components into star-forming and quiescent based on their specific star formation rates. We find that both star-forming discs and star-forming bulges lie on stellar mass-size relations that are similar to those of star-forming galaxies, while quiescent discs are typically smaller than star-forming discs and lie on steeper relations, implying distinct evolutionary mechanisms. Similar to quiescent galaxies, quiescent bulges show a flattening in the stellar mass-size relation at $\sim$10$^{10}$M$_\odot$, below which they show little mass dependence. However, their best-fitting relations have lower normalisations, indicating that at a given mass, bulges are smaller than quiescent galaxies. Finally, we obtain rest-frame colours for individual components, showing that bulges typically have redder colours than discs, as expected. We visually derive UVJ criteria to separate star-forming and quiescent components and show that this separation agrees well with component colour. HFF bulge+disc decomposition catalogues used for these analyses are publicly released with this paper. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 30 pages, 18 figures, and 6 tables. Resubmitted to MNRAS after addressing a thorough and constructive referee report

REBELS-25: Discovery of a dynamically cold disc galaxy at z = 7.31

Authors: Lucie E. Rowland, Jacqueline Hodge, Rychard Bouwens, Pavel Mancera Piña, Alexander Hygate, Hiddo Algera, Manuel Aravena, Rebecca Bowler, Elisabete da Cunha, Pratika Dayal, Andrea Ferrara, Thomas Herard-Demanche, Hanae Inami, Ivana van Leeuwen, Ilse de Looze, Pascal Oesch, Andrea Pallottini, Siân Phillips, Matus Rybak, Sander Schouws, Renske Smit, Laura Sommovigo, Mauro Stefanon, Paul van der Werf

Abstract: We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$μ$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive (… ▽ More We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$μ$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{M_{\odot}}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{M_{\odot}} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low velocity dispersion ($\barσ = 33 \pm 9$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\barσ = 11 ^{+8}_{-4}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a Sérsic index, $n$, of $1.3 \pm 0.2$, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

JWST FRESCO: a comprehensive census of H$β$+[OIII] emitters at 6.8<z<9.0 in the GOODS fields

Authors: R. A. Meyer, P. A. Oesch, E. Giovinazzo, A. Weibel, G. Brammer, J. Matthee, R. P. Naidu, R. J. Bouwens, J. Chisholm, A. Covelo-Paz, Y. Fudamoto, M. Maseda, E. Nelson, I. Shivaei, M. Xiao, T. Herard-Demanche, G. D. Illingworth, J. Kerutt, I. Kramarenko, I. Labbe, E. Leonova, D. Magee, J. Matharu, G. Prieto Lyon, N. Reddy , et al. (5 additional authors not shown)

Abstract: We present the census of H$β$+[O III] $4960,5008$ Åemitters at $6.8<z<9.0$ from the JWST FRESCO survey over 124 arcmin$^2$ in the GOODS-North and GOODS-South fields. Our unbiased spectroscopic search results in 137 spectroscopically-confirmed galaxies at $6.8<z<9.0$ with observed [O III] fluxes $f_{[O III]}\gtrsim 1\times 10^{-18}\ \rm{erg}\ \rm{s}^{-1} \ \rm{cm}^{-2}$. The rest-frame optical line… ▽ More We present the census of H$β$+[O III] $4960,5008$ Åemitters at $6.8<z<9.0$ from the JWST FRESCO survey over 124 arcmin$^2$ in the GOODS-North and GOODS-South fields. Our unbiased spectroscopic search results in 137 spectroscopically-confirmed galaxies at $6.8<z<9.0$ with observed [O III] fluxes $f_{[O III]}\gtrsim 1\times 10^{-18}\ \rm{erg}\ \rm{s}^{-1} \ \rm{cm}^{-2}$. The rest-frame optical line ratios of the median stacked spectrum indicate negligible dust attenuation, low metallicity ($12+\log(\rm{O/H})= 7.2-7.7$) and a high ionisation parameter $\log_{10}U \simeq -2.5$ at a median UV magnitude $M_{\rm{UV}}=-19.65^{+0.59}_{-1.05}$. We find a factor $\times\ 1.3$ difference in the number density of $6.8<z<9.0$ galaxies between GOODS-South and GOODS-North, which is caused by single overdensity at $7.0<z<7.2$ in GOODS-North. The bright end of the UV luminosity function of spectroscopically-confirmed [O III] emitters is in good agreement with that from pre-JWST dropout-selected samples. Discrepancies between the observed [O III] LF, [O III] /UV ratio and [O III] equivalent widths distribution and that predicted by theoretical models suggest burstier star-formation histories and/or more heterogeneous metallicity and ionising conditions in $z>7$ galaxies. We report a rapid decline of the [O III] luminosity density at $z\gtrsim 6-7$ which cannot be explained solely by the evolution of the cosmic star-formation rate density. Finally, we find that FRESCO, in only $2$h, captures star-forming galaxies likely accounting for $\sim 10-20\%$ of the ionising budget at $z=7$ and $z=8$, raising the prospect of detecting directly all the sources of reionisation with JWST. △ Less

Submitted 16 May, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: 18 pages + appendices. Submitted to MNRAS. Public catalogue release at https://github.com/rameyer/fresco . V2: minor edits to Figs. 14, 15, 19

MAGAZ3NE: Massive, Extremely Dusty Galaxies at $z\sim2$ Lead to Photometric Overestimation of Number Densities of the Most Massive Galaxies at $3<z<4$

Authors: Ben Forrest, M. C. Cooper, Adam Muzzin, Gillian Wilson, Danilo Marchesini, Ian McConachie, Percy Gomez, Marianna Annunziatella, Z. Cemile Marsan, Joey Braspenning, Wenjun Chang, Gabriella de Lucia, Fabio Fontanot, Michaela Hirschmann, Dylan Nelson, Annalisa Pillepich, Joop Schaye, Stephanie M. Urbano Stawinski, Mauro Stefanon, Lizhi Xie

Abstract: We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distribution… ▽ More We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distributions (SEDs) from the UltraVISTA and VIDEO surveys. In contrast to previous spectroscopic observations of blue star-forming and post-starburst ultramassive galaxies, candidates in this sample have very red SEDs implying significant dust attenuation, old stellar ages, and/or active galactic nuclei (AGN). Of these galaxies, eight are revealed to be heavily dust-obscured $2.0<z<2.7$ galaxies with strong emission lines, some showing broad features indicative of AGN, three are Type I AGN hosts at $z>3$, one is a $z\sim1.2$ dusty galaxy, and four galaxies do not have a confirmed spectroscopic redshift. In fact, none of the sample has |$z_{\rm spec}-z_{\rm phot}$|$<0.5$, suggesting difficulties for photometric redshift programs in fitting similarly red SEDs. The prevalence of these red interloper galaxies suggests that the number densities of high-mass galaxies are overestimated at $z\gtrsim3$ in large photometric surveys, helping to resolve the `impossibly early galaxy problem' and leading to much better agreement with cosmological galaxy simulations. A more complete spectroscopic survey of ultramassive galaxies is required to pin down the uncertainties on their number densities in the early universe. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ. Comments welcome!

The JWST-SUSPENSE Ultradeep Spectroscopic Program: Survey Overview and Star-Formation Histories of Quiescent Galaxies at 1 < z < 3

Authors: Martje Slob, Mariska Kriek, Aliza G. Beverage, Katherine A. Suess, Guillermo Barro, Rachel Bezanson, Gabriel Brammer, Chloe M. Cheng, Charlie Conroy, Anna de Graaff, Natascha M. Förster Schreiber, Marijn Franx, Brian Lorenz, Pavel E. Mancera Piña, Danilo Marchesini, Adam Muzzin, Andrew B. Newman, Sedona H. Price, Alice E. Shapley, Mauro Stefanon, Pieter van Dokkum, Daniel R. Weisz

Abstract: We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quie… ▽ More We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quiescent galaxy candidates ($z=1-3$, $H_{AB}\le23$), as well as 53 star-forming galaxies at $z=1-4$. With 16~hr of integration and the G140M-F100LP dispersion-filter combination, we observe numerous Balmer and metal absorption lines for all quiescent candidates. We derive stellar masses (log$M_*/M_{\odot}\sim10.2-11.5$) and detailed star-formation histories (SFHs) and show that all 20 candidate quiescent galaxies indeed have quenched stellar populations. These galaxies show a variety of mass-weighted ages ($0.8-3.3$~Gyr) and star formation timescales ($\sim0.5-4$~Gyr), and four out of 20 galaxies were already quiescent by $z=3$. On average, the $z>1.75$ $[z<1.75]$ galaxies formed 50\% of their stellar mass before $z=4$ $[z=3]$. Furthermore, the typical SFHs of galaxies in these two redshift bins ($z_{\text{mean}}=2.2~[1.3]$) indicate that galaxies at higher redshift formed earlier and over shorter star-formation timescales compared to lower redshifts. Although this evolution is naturally explained by the growth of the quiescent galaxy population over cosmic time, number density calculations imply that mergers and/or late-time star formation also contribute to the evolution. In future work, we will further unravel the early formation, quenching, and late-time evolution of these galaxies by extending this work with studies on their chemical abundances, resolved stellar populations and kinematics. △ Less

Submitted 18 July, 2024; v1 submitted 18 April, 2024; originally announced April 2024.

Comments: Accepted in ApJ; 25 pages, 14 figures, 2 tables (excluding appendices)

Quiescent or dusty? Unveiling the nature of extremely red galaxies at $z>3$

Authors: L. Barrufet, P. Oesch, R. Marques-Chaves, K. Arellano-Cordova, J. F. W. Baggen, A. C. Carnall, F. Cullen, J. S. Dunlop, R. Gottumukkala, Y. Fudamoto, G. D. Illingworth, D. Magee, R. J. McLure, D. J. McLeod, M. J. Michałowski, M. Stefanon, P. G. van Dokkum, A. Weibel

Abstract: The advent of the JWST has revolutionised our understanding of high-redshift galaxies. In particular, the NIRCam instrument on-board JWST has revealed a population of Hubble Space Telescope (HST)-dark galaxies that had previously evaded optical detection, potentially due to significant dust obscuration, quiescence, or simply extreme redshift. Here, we present the first NIRSpec spectra of 23 HST-da… ▽ More The advent of the JWST has revolutionised our understanding of high-redshift galaxies. In particular, the NIRCam instrument on-board JWST has revealed a population of Hubble Space Telescope (HST)-dark galaxies that had previously evaded optical detection, potentially due to significant dust obscuration, quiescence, or simply extreme redshift. Here, we present the first NIRSpec spectra of 23 HST-dark galaxies ($\mathrm{H-F444W>1.75}$), unveiling their nature and physical properties. This sample includes both dusty and quiescent galaxies with spectroscopic data from NIRSpec/PRISM, providing accurate spectroscopic redshifts with $\mathrm{\overline{z}_{spec} = 4.1 \pm 0.7}$. The spectral features demonstrate that, while the majority of HST-dark galaxies are dusty, a substantial fraction, $\mathrm{13^{+9}_{-6} \%}$, are quiescent. For the dusty galaxies, we have quantified the dust attenuation using the Balmer decrement ($\mathrm{Hα/ Hβ}$), finding attenuations $\mathrm{A_{V} > 2\ mag}$. We find that HST-dark dusty galaxies are $\mathrm{Hα}$ emitters with equivalent widths spanning the range $\mathrm{ 68 A < EW_{Hα} < 550 A }$, indicative of a wide range of recent star-formation activity. Whether dusty or quiescent, we find that HST-dark galaxies are predominantly massive, with 85\% of the galaxies in the sample having masses $\mathrm{log(M_{*}/M_{\odot}) > 9.8}$. This pilot NIRSpec program reveals the diverse nature of HST-dark galaxies and highlights the effectiveness of NIRSpec/PRISM spectroscopic follow-up in distinguishing between dusty and quiescent galaxies and properly quantifying their physical properties. Upcoming research utilising higher-resolution NIRSpec data and combining JWST with ALMA observations will enhance our understanding of these enigmatic and challenging sources. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 13 pages, 11 figures

Galaxy Build-up in the first 1.5 Gyr of Cosmic History: Insights from the Stellar Mass Function at $z\sim4-9$ from JWST NIRCam Observations

Authors: Andrea Weibel, Pascal A. Oesch, Laia Barrufet, Rashmi Gottumukkala, Richard S. Ellis, Paola Santini, John R. Weaver, Natalie Allen, Rychard Bouwens, Rebecca A. A. Bowler, Gabe Brammer, Adam C. Carnall, Fergus Cullen, Pratika Dayal, Callum T. Donnan, James S. Dunlop, Mauro Giavalisco, Norman A. Grogin, Garth D. Illingworth, Anton M. Koekemoer, Ivo Labbe, Danilo Marchesini, Derek J. McLeod, Ross J. McLure, Rohan P. Naidu , et al. (4 additional authors not shown)

Abstract: Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim$500 arcmin$^2$, we obtain a sample of $>$30,000 galaxies at $z\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $β$ distributions between JWST- and HST-selected samples, we generally… ▽ More Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim$500 arcmin$^2$, we obtain a sample of $>$30,000 galaxies at $z\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $β$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a small population of UV-red galaxies that was missing from previous HST-based LBG samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ and show that they are broadly consistent with existing literature results. However, UV-red galaxies dominate the high-mass end of the SMF at least out to $z\sim6$. In particular the most massive galaxies typically show very red colors between $λ_{obs}\sim1.5μ$m and $\sim4.5μ$m, and thus JWST's unprecedented resolution and sensitivity at these wavelengths yields more accurate constraints on their abundance and masses. The implied redshift evolution of the high-mass end of the SMF suggests a rapid build-up of massive dust-obscured as well as quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). We find the SMFs to be steep over the entire redshift range, and slightly steepening with redshift from $z\sim 4-6$, reaching values of $\approx-2$ at $z\gtrsim6$. Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ in the $\sim1$ Gyr of cosmic time from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 29 pages, 17 figures, submitted to MNRAS

Metal and dust evolution in ALMA REBELS galaxies: insights for future JWST observations

Authors: Marco Palla, Ilse De Looze, Monica Relaño, Stefan van der Giessen, Pratika Dayal, Andrea Ferrara, Raffaella Schneider, Luca Graziani, Hiddo S. B. Algera, Manuel Aravena, Rebecca A. A. Bowler, Alexander P. S. Hygate, Hanae Inami, Ivana van Leeuwen, Rychard Bouwens, Jacqueline Hodge, Renske Smit, Mauro Stefanon, Paul van der Werf

Abstract: ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed c… ▽ More ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA REBELS survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the buildup of stellar mass by using non-parametric SFHs for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [CII]158micron line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other ISM dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust buildup in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 16 pages + appendices, 9 Figures, 1 Table. Resubmitted to MNRAS after moderate revision

Unveiling the hidden universe with JWST: The contribution of dust-obscured galaxies to the stellar mass function at $z\sim3-8$

Authors: R. Gottumukkala, L. Barrufet, P. A. Oesch, A. Weibel, N. Allen, B. Alcalde Pampliega, E. J. Nelson, C. C. Williams, G. Brammer, Y. Fudamoto, V. González, K. E. Heintz, G. Illingworth, D. Magee, R. P. Naidu, M. Shuntov, M. Stefanon, S. Toft, F. Valentino, M. Xiao

Abstract: With the advent of JWST, we can probe the rest-frame optical emission of galaxies at $z>3$ with high sensitivity and spatial resolution, making it possible to accurately characterise red, optically-faint galaxies and thus move towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/CEERS survey, colour… ▽ More With the advent of JWST, we can probe the rest-frame optical emission of galaxies at $z>3$ with high sensitivity and spatial resolution, making it possible to accurately characterise red, optically-faint galaxies and thus move towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/CEERS survey, colour-selected using solely JWST bands. With deep JWST/NIRCam data from 1.15$μ$m to 4.44$μ$m and ancillary HST/ACS and WFC3 data, we determine the physical properties of our sample using spectral energy distribution fitting with BAGPIPES. We demonstrate that our selection method efficiently identifies massive ($\mathrm{\langle \log M_\star/M_\odot \rangle \sim 10}$) and dusty ($\mathrm{\langle A_V\rangle \sim 2.7\ mag}$) sources, with a majority at $z>3$ and predominantly lying on the galaxy main-sequence. The main results of this work are the stellar mass functions (SMF) of red, optically-faint galaxies from redshifts between $3<z<8$: these galaxies make up a significant relative fraction of the pre-JWST total SMF at $3<z<4$ and $4<z<6$, and dominate the high-mass end of the pre-JWST SMF at $6<z<8$, suggesting that our census of the galaxy population needs amendment at these epochs. While larger areas need to be surveyed in the future, our results suggest already that the integrated stellar mass density at $\mathrm{\log M_\star/M_\odot\geq9.25}$ may have been underestimated in pre-JWST studies by up to $\sim$15-20\% at $z\sim3-6$, and up to $\sim$45\% at $z\sim6-8$, indicating the rapid onset of obscured stellar mass assembly in the early universe. △ Less

Submitted 13 June, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures, published in MNRAS; view published article at https://academic.oup.com/mnras/article/530/1/966/7630203

The UNCOVER Survey: A First-look HST+JWST Catalog of Galaxy Redshifts and Stellar Population Properties Spanning $0.2 \lesssim z \lesssim 15$

Authors: Bingjie Wang, Joel Leja, Ivo Labbé, Rachel Bezanson, Katherine E. Whitaker, Gabriel Brammer, Lukas J. Furtak, John R. Weaver, Sedona H. Price, Adi Zitrin, Hakim Atek, Dan Coe, Sam E. Cutler, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Danilo Marchesini, Marijn Franx, Natascha Förster Schreiber, Seiji Fujimoto, Marla Geha, Karl Glazebrook, Anna de Graaff, Jenny E. Greene, Stéphanie Juneau , et al. (19 additional authors not shown)

Abstract: The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 to create the deepest view of our universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to $z \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame c… ▽ More The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 to create the deepest view of our universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to $z \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of $0.2 \lesssim z \lesssim 15$. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector-$β$ priors on masses and star formation histories to produce joint redshift and stellar population posteriors, and additionally lensing magnification is performed on-the-fly to ensure consistency with the scale-dependent priors. We show that this approach produces excellent photometric redshifts with $σ_{\rm NMAD} \sim 0.03$, of a similar quality to the established photometric redshift code EAzY. In line with the open-source scientific objective of the Treasury survey, we publicly release the stellar population catalog with this paper, derived from the photometric catalog adapting aperture sizes based on source profiles. This release includes posterior moments, maximum-likelihood spectra, star-formation histories, and full posterior distributions, offering a rich data set to explore the processes governing galaxy formation and evolution over a parameter space now accessible by JWST. △ Less

Submitted 16 April, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: Corrected typos: Eq.1 should've been (1-kappa)^2, and the lens maps are normalized to D_ds/D_s=1. These errors were only in the writing; no data products or results were affected. The SPS catalogs are accessible via the UNCOVER survey webpage: https://jwst-uncover.github.io/DR2.html#SPSCatalogs, with a copy deposited to Zenodo: https://doi.org/10.5281/zenodo.8401181

Journal ref: The Astrophysical Journal Supplement Series, 270, 12 (2024)

The ALMA REBELS survey: obscured star formation in massive Lyman-break galaxies at z = 4-8 revealed by the IRX-$β$ and $M_{\star}$ relations

Authors: R. A. A. Bowler, H. Inami, L. Sommovigo, R. Smit, H. S. B. Algera, M. Aravena, L. Barrufet, R. Bouwens, E. da Cunha, F. Cullen, P. Dayal, I. de Looze, J. S. Dunlop, Y. Fudamoto, V. Mauerhofer, R. J. McLure, M. Stefanon, R. Schneider, A. Ferrara, L. Graziani, J. A. Hodge, T. Nanayakkara, M. Palla, S. Schouws, D. P. Stark , et al. (1 additional authors not shown)

Abstract: We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-… ▽ More We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($β\simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $β_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$β$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$β$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$β$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs. △ Less

Submitted 28 November, 2023; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Updated to match MNRAS accepted version after minor corrections

Mapping dusty galaxy growth at $z>5$ with FRESCO: Detection of H$α$ in submm galaxy HDF850.1 and the surrounding overdense structures

Authors: Thomas Herard-Demanche, Rychard J. Bouwens, Pascal A. Oesch, Rohan P. Naidu, Roberto Decarli, Erica J. Nelson, Gabriel Brammer, Andrea Weibel, Mengyuan Xiao, Mauro Stefanon, Fabian Walter, Jorryt Matthee, Romain A. Meyer, Stijn Wuyts, Naveen Reddy, Pablo Arrabal Haro, Helmut Dannerbauer, Alice E. Shapley, John Chisholm, Pieter van Dokkum, Ivo Labbe, Garth Illingworth, Daniel Schaerer, Irene Shivaei

Abstract: We report the detection of a 13$σ$ H$α$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$α$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and s… ▽ More We report the detection of a 13$σ$ H$α$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$α$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and southern component, mirroring that seen in [CII] from the Plateau de Bure Interferometer. Modeling the SED of each component separately, we find that the northern component has a higher mass, star formation rate (SFR), and dust extinction than the southern component. The observed H$α$ emission appears to arise entirely from the less-obscured southern component and shows a similar $Δ$v$\sim$+130 km/s velocity offset to that seen for [CII] relative to the source systemic redshift. Leveraging H$α$-derived redshifts from FRESCO observations, we find that HDF850.1 is forming in one of the richest environments identified to date at $z>5$, with 100 $z=5.17-5.20$ galaxies distributed across 10 structures and a $\sim$(15 cMpc)$^3$ volume. Based on the evolution of analogous structures in cosmological simulations, the $z=5.17-5.20$ structures seem likely to collapse into a single $>$10$^{14}$ $M_{\odot}$ cluster by $z\sim0$. Comparing galaxy properties forming within this overdensity with those outside, we find the masses, SFRs, and $UV$ luminosities inside the overdensity to be clearly higher. The prominence of H$α$ line emission from HDF850.1 and other known highly-obscured $z>5$ galaxies illustrates the potential of NIRCam-grism programs to map both the early build-up of IR-luminous galaxies and overdense structures. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS. 20 pages, 10 figures and 8 tables (including appendices)

Massive Optically Dark Galaxies Unveiled by JWST Challenge Galaxy Formation Models

Authors: Mengyuan Xiao, Pascal Oesch, David Elbaz, Longji Bing, Erica Nelson, Andrea Weibel, Rohan Naidu, Emanuele Daddi, Rychard Bouwens, Jorryt Matthee, Stijn Wuyts, John Chisholm, Gabriel Brammer, Mark Dickinson, Benjamin Magnelli, Lucas Leroy, Pieter van Dokkum, Daniel Schaerer, Thomas Herard-Demanche, Laia Barrufet, Ryan Endsley, Yoshinobu Fudamoto, Carlos Gómez-Guijarro, Rashmi Gottumukkala, Garth Illingworth , et al. (12 additional authors not shown)

Abstract: Over the past decade, the existence of a substantial population of optically invisible, massive galaxies at $z\gtrsim3$ has been implied from mid-infrared to millimeter observations. With the unprecedented sensitivity of the JWST, such extremely massive galaxy candidates have immediately been identified even at $z>7$, in much larger numbers than expected. These discoveries raised a hot debate. If… ▽ More Over the past decade, the existence of a substantial population of optically invisible, massive galaxies at $z\gtrsim3$ has been implied from mid-infrared to millimeter observations. With the unprecedented sensitivity of the JWST, such extremely massive galaxy candidates have immediately been identified even at $z>7$, in much larger numbers than expected. These discoveries raised a hot debate. If confirmed, early, high-mass galaxies challenge the current models of galaxy formation. However, the lack of spectroscopic confirmations leads to uncertain stellar mass ($M_{\star}$) estimates, and the possible presence of active galactic nuclei (AGN) adds further uncertainty. Here, we present the first sample of 36 dust-obscured galaxies with robust spectroscopic redshifts at $z_{\rm spec}=5-9$ from the JWST FRESCO survey. The three most extreme sources at $z\sim5-6$ ($\sim$1 billion years after the Big Bang) are so massive (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) that they would require, on average, about 50% of the baryons in their halos to be converted into stars -- two to three times higher than even the most efficient galaxies at later times. The extended emission of these galaxies suggests limited contribution by AGN. This population of ultra-massive galaxies accounts for 20% of the total cosmic star formation rate density at $z\sim5-6$, suggesting a substantial proportion of extremely efficient star formation in the early Universe. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: Submitted to Nature. 22 pages, 4 main figures, 7 supplementary figures, 3 supplementary tables. Comments are welcome

UNCOVER: A NIRSpec Identification of a Broad Line AGN at z = 8.50

Authors: Vasily Kokorev, Seiji Fujimoto, Ivo Labbe, Jenny E. Greene, Rachel Bezanson, Pratika Dayal, Erica J. Nelson, Hakim Atek, Gabriel Brammer, Karina I. Caputi, Iryna Chemerynska, Sam E. Cutler, Robert Feldmann, Yoshinobu Fudamoto, Lukas J. Furtak, Andy D. Goulding, Anna de Graaff, Joel Leja, Danilo Marchesini, Tim B. Miller, Themiya Nanayakkara, Pascal Oesch, Richard Pan, Sedona H. Price, David J. Setton , et al. (7 additional authors not shown)

Abstract: Deep observations with JWST have revealed an emerging population of red point-like sources that could provide a link between the postulated supermassive black hole seeds and observed quasars. In this work we present a JWST/NIRSpec spectrum from the JWST Cycle 1 UNCOVER Treasury survey, of a massive accreting black hole at $z=8.50$, displaying a clear broad-line component as inferred from the H$β$… ▽ More Deep observations with JWST have revealed an emerging population of red point-like sources that could provide a link between the postulated supermassive black hole seeds and observed quasars. In this work we present a JWST/NIRSpec spectrum from the JWST Cycle 1 UNCOVER Treasury survey, of a massive accreting black hole at $z=8.50$, displaying a clear broad-line component as inferred from the H$β$ line with FWHM = $3439\pm413$ km s$^{-1}$, typical of the broad line region of an active galactic nucleus (AGN). The AGN nature of this object is further supported by high ionization, as inferred from emission lines, and a point-source morphology. We compute the black hole mass of log$_{10}(M_{\rm BH}/M_\odot)=8.17\pm0.42$, and a bolometric luminosity of $L_{\rm bol}\sim6.6\times10^{45}$ erg s$^{-1}$. These values imply that our object is accreting at $\sim 40\%$ of the Eddington limit. Detailed modeling of the spectral energy distribution in the optical and near-infrared, together with constraints from ALMA, indicate an upper limit on the stellar mass of log$_{10}(M_{\rm *}/M_\odot)<8.7$, which would lead to an unprecedented ratio of black hole to host mass of at least $\sim 30 \%$. This is orders of magnitude higher compared to the local QSOs, but is consistent with recent AGN studies at high redshift with JWST. This finding suggests that a non-negligible fraction of supermassive black holes either started out from massive seeds and/or grew at a super-Eddington rate at high redshift. Given the predicted number densities of high-$z$ faint AGN, future NIRSpec observations of larger samples will allow us to further investigate the galaxy-black hole co-evolution in the early Universe. △ Less

Submitted 15 October, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: 14 pages, 6 figures, 2 tables. Accepted for a publication in ApJL

UNCOVER: A NIRSpec Census of Lensed Galaxies at z=8.50-13.08 Probing a High AGN Fraction and Ionized Bubbles in the Shadow

Authors: Seiji Fujimoto, Bingjie Wang, John Weaver, Vasily Kokorev, Hakim Atek, Rachel Bezanson, Ivo Labbe, Gabriel Brammer, Jenny E. Greene, Iryna Chemerynska, Pratika Dayal, Anna de Graaff, Lukas J. Furtak, Pascal A. Oesch, David J. Setton, Sedona H. Price, Tim B. Miller, Christina C. Williams, Katherine E. Whitaker, Adi Zitrin, Sam E. Cutler, Joel Leja, Richard Pan, Dan Coe, Pieter van Dokkum , et al. (11 additional authors not shown)

Abstract: We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$α$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reporte… ▽ More We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$α$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reported in Atek et al. (2023). Using six sources with multiple emission line detections, we find that the offset of the redshift estimates between the lines and the Ly$α$ break alone with prism can be as large as $\pm0.2$, raising caution in designing future follow-up spectroscopy for the break-only sources. With spec-$z$ confirmed sources in UNCOVER and the literature, we derive lower limits on the rest-frame ultraviolet (UV) luminosity function (LF) at $z\simeq9$-12 and find these lower limits to be consistent with recent photometric measurements. We identify at least two unambiguous and several possible active galactic nucleus (AGN) systems based on X-ray emission, broad line (BL) H$β$, high ionization line (e.g., NIV]1487, CIV1549) detections, and excess in UVLF. This requires the AGN LFs at $z\simeq$ 9-10 to be comparable or even higher than the X-ray AGN LF estimated at $z\sim6$ and indicates a plausible cause of the high abundance of $z>9$ galaxies claimed in recent photometric studies may be AGNs. One UV-luminous source is confirmed at the same redshift as a dusty BL AGN at $z=8.50$ with a physical separation of 380 kpc in the source plane. These two sources show blueward Ly$α$ line or continuum emission, suggesting that they reside in the same ionized bubble with a radius of $7.69\pm0.18$ pMpc. Our results imply that AGNs have a non-negligible contribution to cosmic reionization. △ Less

Submitted 25 August, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: 27 pages, 11 figures, 4 tables, submitted to ApJ (See also arXiv:2308.11610)

Most of the photons that reionized the Universe came from dwarf galaxies

Authors: Hakim Atek, Ivo Labbé, Lukas J. Furtak, Iryna Chemerynska, Seiji Fujimoto, David J. Setton, Tim B. Miller, Pascal Oesch, Rachel Bezanson, Sedona H. Price, Pratika Dayal, Adi Zitrin, Vasily Kokorev, John R. Weaver, Gabriel Brammer, Pieter van Dokkum, Christina C. Williams, Sam E. Cutler, Robert Feldmann, Yoshinobu Fudamoto, Jenny E. Greene, Joel Leja, Michael V. Maseda, Adam Muzzin, Richard Pan , et al. (8 additional authors not shown)

Abstract: The identification of sources driving cosmic reionization, a major phase transition from neutral Hydrogen to ionized plasma around 600-800 Myr after the Big Bang (Dayal et al. 2018, Mason et al. 2019, Robertson et al. 2022), has been a matter of intense debate (Robertson et al. 2022). Some models suggest that high ionizing emissivity and escape fractions ($f_{\rm esc}$) from quasars support their… ▽ More The identification of sources driving cosmic reionization, a major phase transition from neutral Hydrogen to ionized plasma around 600-800 Myr after the Big Bang (Dayal et al. 2018, Mason et al. 2019, Robertson et al. 2022), has been a matter of intense debate (Robertson et al. 2022). Some models suggest that high ionizing emissivity and escape fractions ($f_{\rm esc}$) from quasars support their role in driving cosmic reionization (Madau & Haardt 2015, Mitra et al. 2018). Others propose that the high $f_{\rm esc}$ values from bright galaxies generates sufficient ionizing radiation to drive this process (Naidu et al. 2020). Finally, a few studies suggest that the number density of faint galaxies, when combined with a stellar-mass-dependent model of ionizing efficiency and $f_{\rm esc}$, can effectively dominate cosmic reionization (Finkelstein et al. 2019, Dayal et al. 2020). However, so far, low-mass galaxies have eluded comprehensive spectroscopic studies owing to their extreme faintness. Here we report an analysis of eight ultra-faint galaxies (in a very small field) during the epoch of reionization with absolute magnitudes between $M_{\rm UV}$ $\sim -17$ to $-15$ mag (down to 0.005 $L^{\star}$. We find that faint galaxies during the Universe's first billion years produce ionizing photons with log($ξ_{\mathrm{ion}}$/ Hz erg$^{-1}$) =$25.80\pm 0.14$, a factor of 4 higher than commonly assumed values (Robertson et al. 2015). If this field is representative of the large scale distribution of faint galaxies, the rate of ionizing photons exceeds that needed for reionization, even for escape fractions of order five per cent. △ Less

Submitted 30 April, 2024; v1 submitted 16 August, 2023; originally announced August 2023.

Comments: 29 pages, 7 figures, 2 tables. Published in Nature

Journal ref: Nature Volume 626, 2024, 975-978

The Spitzer Coverage of HSC-Deep with IRAC for Z studies (SHIRAZ) I: IRAC mosaics

Authors: Marianna Annunziatella, Anna Sajina, Mauro Stefanon, Danilo Marchesini, Mark Lacy, Ivo Labbe, Lilianna Houston, Rachel Bezanson, Eiichi Egami, Xiaohui Fan, Duncan Farrah, Jenny Greene, Andy Goulding, Yen-Ting Lin, Xin Liu, Thibaud Moutard, Yoshiaki Ono, Masami Ouchi, Marcin Sawicki, Jason Surace, Katherine Whitaker

Abstract: We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μm mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improve… ▽ More We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μm mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improved photometric redshifts and stellar population parameters at cosmic noon and earlier epochs. The total area mapped by this work is {\sim} 17 deg2 with a mean integration time of {\sim}1200s, providing a median 5σ depth of 23.7(23.3) at 3.6(4.5)μm in AB. We perform SExtractor photometry both on the combined mosaics as well as the single-epoch mosaics taken {\sim}6 months apart. The resultant IRAC number counts show good agreement with previous studies. In combination with the wealth of existing and upcoming spectro-photometric data in these fields, our IRAC mosaics will enable a wide range of galactic evolution and AGN studies. With that goal in mind, we make the combined IRAC mosaics and coverage maps of these three fields publicly available. counts show good agreement with previous studies. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: accepted for publication in AJ

The ALMA REBELS Survey: Discovery of a massive, highly star-forming and morphologically complex ULIRG at $z =7.31$

Authors: A. P. S. Hygate, J. A. Hodge, E. da Cunha, M. Rybak, S. Schouws, H. Inami, M. Stefanon, L. Graziani, R. Schneider, P. Dayal, R. J. Bouwens, R. Smit, R. A. A. Bowler, R. Endsley, V. Gonzalez, P. A. Oesch, D. P. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, A. Ferrara, Y. Fudamoto, J. H. A, I. De Looze, T. Nanayakkara , et al. (5 additional authors not shown)

Abstract: We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμm$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a sig… ▽ More We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμm$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS. 21 pages, 8 figures

The Hα Luminosity Function of Galaxies at z {\sim} 4.5

Authors: Victoria Bollo, Valentino González, Mauro Stefanon, Pascal A. Oesch, Rychard J. Bouwens, Renske Smit, Garth D. Illingworth, Ivo Labbé

Abstract: We present the Hα luminosity function (LF) derived from a large sample of Lyman break galaxies at z {\sim} 4.5 over the GOODS-South and North fields. This study makes use of the new, full-depth Spitzer/IRAC [3.6] and [4.5] imaging from the GOODS Re-ionization Era wide-Area Treasury from the Spitzer program. The Hα flux is derived from the offset between the continuum flux estimated from the best-f… ▽ More We present the Hα luminosity function (LF) derived from a large sample of Lyman break galaxies at z {\sim} 4.5 over the GOODS-South and North fields. This study makes use of the new, full-depth Spitzer/IRAC [3.6] and [4.5] imaging from the GOODS Re-ionization Era wide-Area Treasury from the Spitzer program. The Hα flux is derived from the offset between the continuum flux estimated from the best-fit spectral energy distribution, and the observed photometry in IRAC [3.6]. From these measurements, we build the Hα LF and study its evolution providing the best constraints of this property at high redshift, where spectroscopy of Hα is not yet available. Schechter parameterizations of the Hα LF show a decreasing evolution of {Φ^\star} with redshift, increasing evolution in L{^\star}, and no significant evolution in the faint-end slope at high z. We find that star formation rates (SFRs) derived from Hα are higher than those derived from the rest-frame UV for low SFR galaxies but the opposite happens for the highest SFRs. This can be explained by lower mass galaxies (also lower SFR) having, on average, rising star formation histories (SFHs), while at the highest masses the SFHs may be declining. The SFR function is steeper, and because of the excess SFR(Hα) compared to SFR(UV) at low SFRs, the SFR density estimated from Hα is higher than the previous estimates based on UV luminosities. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: Astrophysical Journal, Accepted, 17 pages, 16 figures

Journal ref: The Astrophysical Journal, 946:117 (17pp), 2023 April 1

The JWST FRESCO Survey: Legacy NIRCam/Grism Spectroscopy and Imaging in the two GOODS Fields

Authors: P. A. Oesch, G. Brammer, R. P. Naidu, R. J. Bouwens, J. Chisholm, G. D. Illingworth, J. Matthee, E. Nelson, Y. Qin, N. Reddy, A. Shapley, I. Shivaei, P. van Dokkum, A. Weibel, K. Whitaker, S. Wuyts, A. Covelo-Paz, R. Endsley, Y. Fudamoto, E. Giovinazzo, T. Herard-Demanche, J. Kerutt, I. Kramarenko, I. Labbe, E. Leonova , et al. (16 additional authors not shown)

Abstract: We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W fi… ▽ More We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R~1600 covering 3.8 to 5.0 $μ$m for most galaxies in the NIRCam field-of-view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z~0.2-0.5, to Pa$α$ and Pa$β$ at z~1-3, HeI and [SIII] at z~2.5-4.5, H$α$ and [NII] at z~5-6.5, up to [OIII] and H$β$ for z~7-9 galaxies, and possibly even [OII] at z~10-12. FRESCO's grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ~28.2 mag (5 $σ$ in 0.32" diameter apertures). Together with this publication, the v1 imaging mosaics are released as high-level science products via MAST. Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 years ago, the discovery of broad-line active galactic nuclei at z>4, and resolved Pa$α$ maps of galaxies at z~1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations. △ Less

Submitted 16 August, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 12 pages, 11 figures; MNRAS in press; for more information on the survey and data releases, see http://jwst-fresco.astro.unige.ch/ and https://archive.stsci.edu/hlsp/fresco

The ALMA REBELS Survey: The First Infrared Luminosity Function Measurement at $\mathbf{z \sim 7}

Authors: L. Barrufet, P. A. Oesch, R. Bouwens, H. Inami, L. Sommovigo, H. Algera, E. da Cunha, M. Aravena, P. Dayal, A. Ferrara, Y. Fudamoto, V. Gonzalez, L. Graziani, A. Hygate, I. de Looze, T. Nanayakkara, A. Pallottini, R. Schneider, M. Stefanon, M. Topping, P. van Der Werf

Abstract: We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit… ▽ More We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 9 pages, 5 figures

Cold Dust and Low [OIII]/[CII] Ratios: an Evolved Star-forming Population at Redshift 7

Authors: Hiddo Algera, Hanae Inami, Laura Sommovigo, Yoshinobu Fudamoto, Raffaella Schneider, Luca Graziani, Pratika Dayal, Rychard Bouwens, Manuel Aravena, Elisabete da Cunha, Andrea Ferrara, Alexander Hygate, Ivana van Leeuwen, Ilse De Looze, Marco Palla, Andrea Pallottini, Renske Smit, Mauro Stefanon, Michael Topping, Paul van der Werf

Abstract: We present new ALMA Band 8 (rest-frame $90\,μ$m) observations of three massive ($M_\star \approx 10^{10}\,M_\odot$) galaxies at $z\approx7$ previously detected in [CII]$158\,μ$m and underlying dust continuum emission in the Reionization Era Bright Emission Line Survey (REBELS). We detect the dust continuum emission of two of our targets in Band 8 (REBELS-25 and REBELS-38), while REBELS-12 remains… ▽ More We present new ALMA Band 8 (rest-frame $90\,μ$m) observations of three massive ($M_\star \approx 10^{10}\,M_\odot$) galaxies at $z\approx7$ previously detected in [CII]$158\,μ$m and underlying dust continuum emission in the Reionization Era Bright Emission Line Survey (REBELS). We detect the dust continuum emission of two of our targets in Band 8 (REBELS-25 and REBELS-38), while REBELS-12 remains undetected. Through modified blackbody fitting we determine cold dust temperatures ($T_\mathrm{dust} \approx 30 - 35\,$K) in both of the dual-band detected targets, given a fiducial model of optically thin emission with $β= 2.0$. Their dust temperatures are lower than most $z\sim7$ galaxies in the literature, and consequently their dust masses are higher ($M_\mathrm{dust} \approx 10^{8}\,M_\odot$). Nevertheless, these large dust masses are still consistent with predictions from models of dust production in the early Universe. In addition, we target and detect [OIII]$88\,μ$m emission in both REBELS-12 and REBELS-25, and find $L_\mathrm{[OIII]} / L_\mathrm{[CII]}$ ratios of approximately unity, low compared to the $L_\mathrm{[OIII]} / L_\mathrm{[CII]} \gtrsim 2 - 10$ observed in the known $z\gtrsim6$ population thus far. We argue the lower line ratios are due to a comparatively weaker ionizing radiation field resulting from the less starbursty nature of our targets. This low burstiness supports the cold dust temperatures and below average $\mathrm{[OIII]}λ\lambda4959,5007 + \mathrm{H}β$ equivalent widths of REBELS-25 and REBELS-38, compared to the known high-redshift population. Overall, this provides evidence for the existence of a massive, dust-rich galaxy population at $z\approx7$ which has previously experienced vigorous star formation, but is currently forming stars in a steady, as opposed to bursty, manner. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 18 pages + appendices; 8 figures in main text; Submitted to MNRAS on 27 Dec 2022

The UNCOVER Survey: A first-look HST+JWST catalog of 60,000 galaxies near Abell 2744 and beyond

Authors: John R. Weaver, Sam E. Cutler, Richard Pan, Katherine E. Whitaker, Ivo Labbe, Sedona H. Price, Rachel Bezanson, Gabriel Brammer, Danilo Marchesini, Joel Leja, Bingjie Wang, Lukas J. Furtak, Adi Zitrin, Hakim Atek, Dan Coe, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Natascha Forster Schreiber, Marijn Franx, Seiji Fujimoto, Yoshinobu Fudamoto, Karl Glazebrook, Anna de Graaff, Jenny E. Greene , et al. (19 additional authors not shown)

Abstract: In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolut… ▽ More In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 60,000 objects across the 49 arcmin$^2$ JWST footprint down to a $5\,σ$ limiting magnitude of $\sim$30 mag in 0.32" apertures. Photometry is performed using circular apertures on images matched to the point spread function of the reddest NIRCam broad band, F444W, and cleaned of bright cluster galaxies and the related intra-cluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a $σ_{\rm NMAD}\approx0.03$ based on known, but relatively small, spectroscopic samples. With this paper, we publicly release our HST and JWST PSF-matched photometric catalog with optimally assigned aperture sizes for easy use, along with single aperture catalogs, photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow-up of some of the first JWST-selected samples in Summer 2023. △ Less

Submitted 2 October, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 28 pages, 19 figures, resubmitted to ApJS following significant data product improvements. Comments welcome. Catalogs can be accessed at https://jwst-uncover.github.io/DR2.html#PhotometricCatalogs

UV Luminosity Density Results at z>8 from the First JWST/NIRCam Fields: Limitations of Early Data Sets and the Need for Spectroscopy

Authors: Rychard Bouwens, Garth Illingworth, Pascal Oesch, Mauro Stefanon, Rohan Naidu, Ivana van Leeuwen, Dan Magee

Abstract: We have derived luminosity functions, and set constraints on the UV luminosity and SFR density from z~17 to z~8, using the three most-studied JWST/NIRCam data sets, the SMACS0723, GLASS Parallel, and CEERS fields. We first used our own selections on two independent reductions of these datasets using the latest calibrations. 18 z~8, 12 z~10, 5 z~13, and 1 z~17 candidate galaxies are identified over… ▽ More We have derived luminosity functions, and set constraints on the UV luminosity and SFR density from z~17 to z~8, using the three most-studied JWST/NIRCam data sets, the SMACS0723, GLASS Parallel, and CEERS fields. We first used our own selections on two independent reductions of these datasets using the latest calibrations. 18 z~8, 12 z~10, 5 z~13, and 1 z~17 candidate galaxies are identified over these fields in our primary reductions, with a similar number of candidates in our secondary reductions. We then use these two reductions, applying a quantitative discriminator, to segregate the full set of z>~8 candidates reported over these fields from the literature, into three different samples, ``robust,'' ``solid,'' and ``possible''. Using all of these samples we then derive UV LF and luminosity density results at $z\geq8$, finding substantial differences. For example, including the full set of ``solid'' and ``possible'' z>~12 candidates from the literature, we find UV luminosity densities which are ~7x and ~20x higher than relying on the ``robust'' candidates alone. These results indicate the evolution of the UV LF and luminosity densities at z>~8 is still extremely uncertain, emphasizing the need for spectroscopy and deeper NIRCam+optical imaging to obtain reliable results. Nonetheless, even with the very conservative ``robust'' approach to selections, both from our own and those of other studies, we find the luminosity density from luminous (M(UV)<-19) galaxies to be ~2x larger than is easily achievable using constant star-formation efficiency models, similar to what other early JWST results have suggested. △ Less

Submitted 1 April, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: 26 pages, 10 figures, 15 tables, accepted for publication in MNRAS, figure 7 shows the key results

The JWST UNCOVER Treasury survey: Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization

Authors: Rachel Bezanson, Ivo Labbe, Katherine E. Whitaker, Joel Leja, Sedona H. Price, Marijn Franx, Gabe Brammer, Danilo Marchesini, Adi Zitrin, Bingjie Wang, John R. Weaver, Lukas J. Furtak, Hakim Atek, Dan Coe, Sam E. Cutler, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Natascha Forster Schreiber, Seiji Fujimoto, Marla Geha, Karl Glazebrook, Anna de Graaff, Jenny E. Greene, Stephanie Juneau , et al. (18 additional authors not shown)

Abstract: In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in November 2022. The UNCOVER survey includes ultradeep ($\sim29-30\mathrm{AB}$) imaging of $\sim$45 arcmin$^2$ on and around the well-studied Abell 2744 galaxy cluster at $z=0.308$ and wil… ▽ More In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in November 2022. The UNCOVER survey includes ultradeep ($\sim29-30\mathrm{AB}$) imaging of $\sim$45 arcmin$^2$ on and around the well-studied Abell 2744 galaxy cluster at $z=0.308$ and will follow-up ${\sim}500$ galaxies with extremely deep low-resolution spectroscopy with the NIRSpec/PRISM during the summer of 2023. We describe the science goals, survey design, target selection, and planned data releases. We also present and characterize the depths of the first NIRCam imaging mosaic, highlighting previously unparalleled resolved and ultradeep 2-4 micron imaging of known objects in the field. The UNCOVER primary NIRCam mosaic spans 28.8 arcmin$^2$ in seven filters (F115W, F150W, F200W, F277W, F356W, F410M, F444W) and 16.8 arcmin$^2$ in our NIRISS parallel (F115W, F150W, F200W, F356W, and F444W). To maximize early community use of the Treasury data set, we publicly release full reduced mosaics of public JWST imaging including 45 arcmin$^2$ NIRCam and 17 arcmin$^2$ NIRISS mosaics on and around the Abell 2744 cluster, including the Hubble Frontier Field primary and parallel footprints. △ Less

Submitted 9 December, 2022; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 17 pages, 8 figures, 4 tables, submitted to ApJ, comments welcome (v2 with full author list in metadata)

Evolution of the UV LF from z~15 to z~8 Using New JWST NIRCam Medium-Band Observations over the HUDF/XDF

Authors: Rychard J. Bouwens, Mauro Stefanon, Gabriel Brammer, Pascal A. Oesch, Thomas Herard-Demanche, Garth D. Illingworth, Jorryt Matthee, Rohan P. Naidu, Pieter G. van Dokkum, Ivana F. van Leeuwen

Abstract: We present the first constraints on the prevalence of z>10 galaxies in the Hubble Ultra Deep Field (HUDF) leveraging new NIRCam observations from JEMS (JWST Extragalactic Medium-band Survey). These NIRCam observations probe redward of 1.6$μ$m, beyond the wavelength limit of HST, allowing us to search for galaxies to z>10. These observations indicate that the highest redshift candidate identified i… ▽ More We present the first constraints on the prevalence of z>10 galaxies in the Hubble Ultra Deep Field (HUDF) leveraging new NIRCam observations from JEMS (JWST Extragalactic Medium-band Survey). These NIRCam observations probe redward of 1.6$μ$m, beyond the wavelength limit of HST, allowing us to search for galaxies to z>10. These observations indicate that the highest redshift candidate identified in the HUDF09 data with HST, UDFj-39546284, has a redshift of z>11.5, as had been suggested in analyses of the HUDF12/XDF data. This has now been confirmed with JWST NIRSpec. This source is thus the most distant galaxy discovered by HST in its >30 years of operation. Additionally, we identify nine other z~8-13 candidate galaxies over the HUDF, two of which are new discoveries that appear to lie at z~11-12. We use these results to characterize the evolution of the UV luminosity function (LF) from z~15 to z~8.7. While our LF results at z~8.7 and z~10.5 are consistent with previous findings over the HUDF, our new LF estimates at z~12.6 are higher than other results in the literature, potentially pointing to a milder evolution in the UV luminosity density from z~12.6. We emphasize that our LF results are uncertain given the small number of z~12.6 sources and limited volume probed. The new NIRCam data also indicate that the faint z~8-13 galaxies in the HUDF/XDF show blue UV-continuum slopes beta~-2.7, high specific star formation rates ~24.5 Gyr**-1, and high EW (~1300A) [OIII]+H$β$ emission, with two z~8.5 sources showing [OIII]+H$β$ EWs of ~2300 A. △ Less

Submitted 8 March, 2023; v1 submitted 4 November, 2022; originally announced November 2022.

Comments: 21 pages, 14 figures, 7 tables, resubmitted to MNRAS, after replying to the referee report

The ALMA REBELS Survey: The Dust-obscured Cosmic Star Formation Rate Density at Redshift 7

Authors: Hiddo Algera, Hanae Inami, Pascal Oesch, Laura Sommovigo, Rychard Bouwens, Michael Topping, Sander Schouws, Mauro Stefanon, Daniel Stark, Manuel Aravena, Laia Barrufet, Elisabete da Cunha, Pratika Dayal, Ryan Endsley, Andrea Ferrara, Yoshinobu Fudamoto, Valentino Gonzalez, Luca Graziani, Jacqueline Hodge, Alexander Hygate, Ilse de Looze, Themiya Nanayakkara, Rafaella Schneider, Paul van der Werf

Abstract: Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and… ▽ More Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization. △ Less

Submitted 9 November, 2022; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: MNRAS Accepted (01 Nov 2022)

MAGAZ3NE: High Stellar Velocity Dispersions for Ultra-Massive Quiescent Galaxies at $z\gtrsim3$

Authors: Ben Forrest, Gillian Wilson, Adam Muzzin, Danilo MArchesini, M. C. Cooper, Z. Cemile Marsan, Marianna Annunziatella, Ian McConachie, Kumail Zaidi, Percy Gomez, Stephanie M. Urbano Stawinski, Wenjun Chang, Gabriella de Lucia, Francesco La Barbera, Lori Lubin, Julie Nantais, Theodore Peña, Paolo Saracco, Jason Surace, Mauro Stefanon

Abstract: In this work we publish stellar velocity dispersions, sizes, and dynamical masses for 8 ultra-massive galaxies (UMGs; log($M$/M$_\odot>11$, $z\gtrsim3$) from the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) Survey, more than doubling the number of such galaxies with velocity dispersion measurements at this epoch. Using the deep Keck/MOSFIRE and Keck/NIRES spectroscopy of these object… ▽ More In this work we publish stellar velocity dispersions, sizes, and dynamical masses for 8 ultra-massive galaxies (UMGs; log($M$/M$_\odot>11$, $z\gtrsim3$) from the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) Survey, more than doubling the number of such galaxies with velocity dispersion measurements at this epoch. Using the deep Keck/MOSFIRE and Keck/NIRES spectroscopy of these objects in the $H$- and $K$-bandpasses, we obtain large velocity dispersions of $\sim400$ km s$^{-1}$ for most of the objects, which are some of the highest stellar velocity dispersions measured, and $\sim40$\% larger than those measured for galaxies of similar mass at $z\sim1.7$. The sizes of these objects are also smaller by a factor of 1.5-3 compared to this same $z\sim1.7$ sample. We combine these large velocity dispersions and small sizes to obtain dynamical masses. The dynamical masses are similar to the stellar masses of these galaxies, consistent with a Chabrier initial mass function (IMF). Considered alongside previous studies of massive quiescent galaxies across $0.2<z<4.0$, there is evidence for an evolution in the relation between the dynamical mass - stellar mass ratio and velocity dispersion as a function of redshift. This implies an IMF with fewer low mass stars (e.g., Chabrier IMF) for massive quiescent galaxies at higher redshifts in conflict with the bottom-heavy IMF (e.g., Salpeter IMF) found in their likely $z\sim0$ descendants, though a number of alternative explanations such as a different dynamical structure or significant rotation are not ruled out. Similar to data at lower redshifts, we see evidence for an increase of IMF normalization with velocity dispersion, though the $z\gtrsim3$ trend is steeper than that for $z\sim0.2$ early-type galaxies and offset to lower dynamical-to-stellar mass ratios. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 17 pages and 8 figures in the main text, plus references and Appendices for a total of 25 pages and 14 figures. Accepted for publication in The Astrophysical Journal

Schrodinger's Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

Authors: Rohan P. Naidu, Pascal A. Oesch, David J. Setton, Jorryt Matthee, Charlie Conroy, Benjamin D. Johnson, John R. Weaver, Rychard J. Bouwens, Gabriel B. Brammer, Pratika Dayal, Garth D. Illingworth, Laia Barrufet, Sirio Belli, Rachel Bezanson, Sownak Bose, Kasper E. Heintz, Joel Leja, Ecaterina Leonova, Rui Marques-Chaves, Mauro Stefanon, Sune Toft, Arjen van der Wel, Pieter van Dokkum, Andrea Weibel, Katherine E. Whitaker

Abstract: $JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}… ▽ More $JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $Λ$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2μ$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($Δz\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

Unveiling the Nature of Infrared Bright, Optically Dark Galaxies with Early JWST Data

Authors: L. Barrufet, P. A. Oesch, A. Weibel, G. Brammer, R. Bezanson, R. Bouwens, Y. Fudamoto, V. Gonzalez, R. Gottumukkala, G. Illingworth, K. E. Heintz, B. Holden, I. Labbe, D. Magee, R. P. Naidu, E. Nelson, M. Stefanon, R. Smit, P. van Dokkum, J. Weaver, C. Williams

Abstract: Over the last few years, both ALMA and Spitzer/IRAC observations have revealed a population of likely massive galaxies at $z>3$ that was too faint to be detected in HST rest-frame ultraviolet imaging. However, due to the very limited photometry for individual galaxies, the true nature of these so-called HST-dark galaxies has remained elusive. Here, we present the first sample of such galaxies obse… ▽ More Over the last few years, both ALMA and Spitzer/IRAC observations have revealed a population of likely massive galaxies at $z>3$ that was too faint to be detected in HST rest-frame ultraviolet imaging. However, due to the very limited photometry for individual galaxies, the true nature of these so-called HST-dark galaxies has remained elusive. Here, we present the first sample of such galaxies observed with very deep, high-resolution NIRCam imaging from the Early Release Science Program CEERS. 30 HST-dark sources are selected based on their red colours across 1.6 $μ$m to 4.4 $μ$m. Their physical properties are derived from 12-band multi-wavelength photometry, including ancillary HST imaging. We find that these galaxies are generally heavily dust-obscured ($A_{V}\sim2$ mag), massive ($\log (M/M_{\odot}) \sim10$), star-forming sources at $z\sim2-8$ with an observed surface density of $\sim0.8$ arcmin$^{-2}$. This suggests that an important fraction of massive galaxies may have been missing from our cosmic census at $z>3$ all the way into the Reionization epoch. The HST-dark sources lie on the main sequence of galaxies and add an obscured star formation rate density (SFRD) of $\mathrm{3.2^{+1.8}_{-1.3} \times 10^{-3} M_{\odot}/yr/Mpc^{3}}$ at $z\sim7$ showing likely presence of dust in the Epoch of Reionization. Our analysis shows the unique power of JWST to reveal this previously missing galaxy population and to provide a more complete census of galaxies at $z=2-8$ based on rest-frame optical imaging. △ Less

Submitted 15 March, 2023; v1 submitted 29 July, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures, accepted to MNRAS

A population of red candidate massive galaxies ~600 Myr after the Big Bang

Authors: Ivo Labbe, Pieter van Dokkum, Erica Nelson, Rachel Bezanson, Katherine Suess, Joel Leja, Gabriel Brammer, Katherine Whitaker, Elijah Mathews, Mauro Stefanon, Bingjie Wang

Abstract: Galaxies with stellar masses as high as $\sim 10^{11}$ solar masses have been identified out to redshifts $z \sim 6$, approximately one billion years after the Big Bang. It has been difficult to find massive galaxies at even earlier times, as the Balmer break region, which is needed for accurate mass estimates, is redshifted to wavelengths beyond $2.5\mum$. Here we make use of the $1-5\mum$ covera… ▽ More Galaxies with stellar masses as high as $\sim 10^{11}$ solar masses have been identified out to redshifts $z \sim 6$, approximately one billion years after the Big Bang. It has been difficult to find massive galaxies at even earlier times, as the Balmer break region, which is needed for accurate mass estimates, is redshifted to wavelengths beyond $2.5\mum$. Here we make use of the $1-5\mum$ coverage of the JWST early release observations to search for intrinsically red galaxies in the first ~750 million years of cosmic history. In the survey area, we find six candidate massive galaxies (stellar mass $>10^{10}$ solar masses) at $7.4 < z < 9.1$, 500 - 700 Myr after the Big Bang, including one galaxy with a possible stellar mass of $\sim 10^{11}$ solar masses. If verified with spectroscopy, the stellar mass density in massive galaxies would be much higher than anticipated from previous studies based on rest-frame ultraviolet-selected samples. △ Less

Submitted 31 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication in Nature. Updated to match the accepted version. Significant revision with updated calibration and stellar masses. Coordinates and fluxes of the sources can be found here: https://github.com/ivolabbe/red-massive-candidates

Two Remarkably Luminous Galaxy Candidates at $z\approx10-12$ Revealed by JWST

Authors: Rohan P. Naidu, Pascal A. Oesch, Pieter van Dokkum, Erica J. Nelson, Katherine A. Suess, Gabriel Brammer, Katherine E. Whitaker, Garth Illingworth, Rychard Bouwens, Sandro Tacchella, Jorryt Matthee, Natalie Allen, Rachel Bezanson, Charlie Conroy, Ivo Labbe, Joel Leja, Ecaterina Leonova, Dan Magee, Sedona H. Price, David J. Setton, Victoria Strait, Mauro Stefanon, Sune Toft, John R. Weaver, Andrea Weibel

Abstract: The first few hundred Myrs at $z>10$ mark the last major uncharted epoch in the history of the Universe, where only a single galaxy (GNz11 at $z\approx11$) is currently spectroscopically confirmed. Here we present a search for luminous $z>10$ galaxies with $JWST$/NIRCam photometry spanning $\approx1-5μ$m and covering 49 arcmin$^{2}$ from the public Early Release Science programs (CEERS and GLASS).… ▽ More The first few hundred Myrs at $z>10$ mark the last major uncharted epoch in the history of the Universe, where only a single galaxy (GNz11 at $z\approx11$) is currently spectroscopically confirmed. Here we present a search for luminous $z>10$ galaxies with $JWST$/NIRCam photometry spanning $\approx1-5μ$m and covering 49 arcmin$^{2}$ from the public Early Release Science programs (CEERS and GLASS). Our most secure candidates are two $M_{\rm{UV}}\approx-21$ systems: GLASS-z12 and GLASS-z10. These galaxies display abrupt $\gtrsim1.8$ mag breaks in their spectral energy distributions, consistent with complete absorption of flux bluewards of Lyman-$α$ that is redshifted to $z=12.4^{+0.1}_{-0.3}$ and $z=10.4^{+0.4}_{-0.5}$. Lower redshift interlopers such as quiescent galaxies with strong Balmer breaks would be comfortably detected at $>5σ$ in multiple bands where instead we find no flux. From SED modeling we infer that these galaxies have already built up $\sim 10^9$ solar masses in stars over the $\lesssim300-400$ Myrs after the Big Bang. The brightness of these sources enable morphological constraints. Tantalizingly, GLASS-z10 shows a clearly extended exponential light profile, potentially consistent with a disk galaxy of $r_{\rm{50}}\approx0.7$ kpc. These sources, if confirmed, join GNz11 in defying number density forecasts for luminous galaxies based on Schechter UV luminosity functions, which require a survey area $>10\times$ larger than we have studied here to find such luminous sources at such high redshifts. They extend evidence from lower redshifts for little or no evolution in the bright end of the UV luminosity function into the cosmic dawn epoch, with implications for just how early these galaxies began forming. This, in turn, suggests that future deep $JWST$ observations may identify relatively bright galaxies to much earlier epochs than might have been anticipated. △ Less

Submitted 25 October, 2022; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: Accepted for publication in ApJL. Analysis updated with in-flight NIRCam chip-wise zero-points -- Figure 6 shows a comparison with v1

Dual constraints with ALMA: new [O III] 88 ${\rm μ}$m and dust-continuum observations reveal the ISM conditions of luminous LBGs at $z \sim 7$

Authors: Joris Witstok, Renske Smit, Roberto Maiolino, Nimisha Kumari, Manuel Aravena, Leindert Boogaard, Rychard Bouwens, Stefano Carniani, Jacqueline A. Hodge, Gareth C. Jones, Mauro Stefanon, Paul van der Werf, Sander Schouws

Abstract: We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison wi… ▽ More We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [O III] tracing ionised gas. By contrast, we find more spatially extended [C II] emission likely produced in neutral gas, as indicated by a [N II] 205 ${\rm μ}$m non-detection in one source. We find a positive correlation between the equivalent width of the optical [O III] and H${\rm β}$ lines and the [O III]/[C II] ratio, as seen in local metal-poor dwarf galaxies. Cloudy models of a nebula of typical density harbouring a young stellar population with a high ionisation parameter appear to adequately reproduce the far-infrared lines. Surprisingly, however, our models fail to reproduce the strength of [O III] 88 ${\rm μ}$m, unless we assume an ${\rm α}$/Fe enhancement and a near-solar nebular oxygen abundance. On spatially resolved scales, we find [O III]/[C II] shows a tentative anti-correlation with infrared excess, $L_{\rm IR}/L_{\rm UV}$, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code MERCURIUS to show that dust-continuum measurements of one source appear to favour a low dust temperature coupled with a high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae. △ Less

Submitted 5 July, 2022; v1 submitted 30 June, 2022; originally announced July 2022.

Comments: 23 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Deep Spitzer/IRAC Data for z~10 galaxies Reveal Blue Balmer Break Colors: Young Stellar Populations at ~500 Myr of Cosmic Time

Authors: Mauro Stefanon, Rychard J. Bouwens, Ivo Labbé, Garth D. Illingworth, Valentino Gonzalez, Pascal A. Oesch

Abstract: We present the deepest constraints yet on the median rest-UV+optical SED of $z\sim10$ galaxies, prior to JWST science operations. We constructed stacks based on four robust $J_{125}$-dropouts, previously identified across the GOODS fields. We used archival HST/WFC3 data and the full depth Spitzer/IRAC mosaics from the GREATS program, the deepest coverage at $\sim3-5μ$m to date. The most remarkable… ▽ More We present the deepest constraints yet on the median rest-UV+optical SED of $z\sim10$ galaxies, prior to JWST science operations. We constructed stacks based on four robust $J_{125}$-dropouts, previously identified across the GOODS fields. We used archival HST/WFC3 data and the full depth Spitzer/IRAC mosaics from the GREATS program, the deepest coverage at $\sim3-5μ$m to date. The most remarkable feature of the SED is a blue IRAC $[3.6]-[4.5]=-0.18\pm0.25$ mag color. We also find a nearly flat $H_{160}-[3.6]=0.07\pm0.22$ mag color, corresponding to a UV slope $β= -1.92\pm0.25$. This is consistent with previous studies, and indicative of minimal dust absorption. The observed blue IRAC color and SED fitting suggest that $z\sim10$ galaxies have very young (few $\times10$ Myr) stellar populations, with $80\%$ of stars being formed in the last $\lesssim 160$ Myr ($2σ$). While an exciting result, the uncertainties on the SED are too large to allow us to place strong constraints on the presence of a nebular continuum in $z\sim10$ galaxies (as might be suggested by the blue $[3.6]-[4.5] < 0$ mag color). The resulting sSFR is consistent with the specific accretion rate of dark matter halos, indicative of a star-formation efficiency showing quite limited evolution at such early epochs. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Submitted to ApJ. 9 pages, 3 figures. Comments welcome

The ALMA REBELS Survey: The Cosmic HI Gas Mass Density in Galaxies at $z\approx 7$

Authors: K. E. Heintz, P. A. Oesch, M. Aravena, R. J. Bouwens, P. Dayal, A. Ferrara, Y. Fudamoto, L. Graziani, H. Inami, L. Sommovigo, R. Smit, M. Stefanon, M. Topping, A. Pallottini, P. van der Werf

Abstract: The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine HI 21-cm transition. Here we estimate the HI gas mass of a sample of main-sequence star-forming galaxies at $z\sim 6.5 - 7.8$ surveyed for [CII]$-158μ$m emission as part of the Reionization Era Bright Emission Line Survey (REBELS),… ▽ More The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine HI 21-cm transition. Here we estimate the HI gas mass of a sample of main-sequence star-forming galaxies at $z\sim 6.5 - 7.8$ surveyed for [CII]$-158μ$m emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [CII]-to-HI conversion factor. We find that the HI gas mass excess in galaxies increases as a function of redshift, with an average of $M_{\rm HI} / M_\star \approx 10$, corresponding to HI gas mass fractions of $f_{\rm HI} = M_{\rm HI} / (M_\star + M_{\rm HI}) = 90\%$, at $z\approx 7$. Based on the [CII]-$158μ$m luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic HI gas mass density in galaxies ($ρ_{\rm HI}$) at this redshift, which we measure to be $ρ_{\rm HI} = 7.1^{+6.4}_{-3.0} \times 10^{6}\,M_{\odot}\,{\rm Mpc^{-3}}$. This estimate is substantially lower by a factor of $\approx 10$ than that inferred from an extrapolation of damped Lyman-$α$ absorber (DLA) measurements, and largely depend on the exact [CII] LF adopted. However, we find this decrease in $ρ_{\rm HI}$ to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sightlines predominantly probing the substantial fraction of HI gas in high-$z$ galactic halos, whereas [CII] traces the HI in the ISM associated with star formation. We make predictions for this build-up of neutral gas in galaxies as a function of redshift, showing that at $z\gtrsim 5$ only $\approx 10\%$ of the cosmic HI gas content is confined in galaxies and associated with the star-forming ISM. △ Less

Submitted 27 July, 2022; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: The accepted version of the manuscript in press for ApJL

The ALMA REBELS Survey: Average [CII] $158\,\rm{μm}$ sizes of Star-Forming Galaxies from $z\sim 7$ to $z\sim 4$

Authors: Y. Fudamoto, R. Smit, R. A. A. Bowler, P. A. Oesch, R. Bouwens, M. Stefanon, H. Inami, R. Endsley, V. Gonzalez, S. Schouws, D. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, E. da Cunha, P. Dayal, A. Ferrara, L. Graziani, J. A. Hodge, A. P. S. Hygate, A. K. Inoue, T. Nanayakkara, A. Pallottini, E. Pizzati, R. Schneider , et al. (14 additional authors not shown)

Abstract: We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective ra… ▽ More We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm μm}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range. △ Less

Submitted 15 June, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 9pages, 4 figures, submitted to ApJ; minor corrections were made for typos in the author list and the title

z~2-9 Galaxies Magnified by the Hubble Frontier Field Clusters II: Luminosity Functions and Constraints on a Faint-End Turnover

Authors: R. J. Bouwens, G. D. Illingworth, R. S. Ellis, P. A. Oesch, M. Stefanon

Abstract: We present new determinations of the rest-UV luminosity functions (LFs) at z=2-9 to extremely low luminosities (>-14 mag) from a sample of >2500 lensed galaxies found behind the HFF clusters. For the first time, we present faint-end slope results from lensed samples that are fully consistent with blank-field results over the redshift range z=2-9, while reaching to much lower luminosities than poss… ▽ More We present new determinations of the rest-UV luminosity functions (LFs) at z=2-9 to extremely low luminosities (>-14 mag) from a sample of >2500 lensed galaxies found behind the HFF clusters. For the first time, we present faint-end slope results from lensed samples that are fully consistent with blank-field results over the redshift range z=2-9, while reaching to much lower luminosities than possible from the blank-field studies. Combining the deep lensed sample with the large blank-field samples allows us to set the tight constraints on the faint-end slope alpha of the z=2-9 UV LFs and its evolution. We find a smooth flattening in alpha from -2.28+/-0.10 (z=9) to -1.53+/-0.03 (z=2) with cosmic time (d(alpha)/dz=-0.11+/-0.01), fully consistent with dark matter halo buildup. We utilize these new results to present new measurements of the evolution in the UV luminosity density rho(UV) brightward of -13 mag from z~9 to z~2. Accounting for the SFR densities to faint luminosities implied by our LF results, we find that unobscured star formation dominates the SFR density at z>~4, with obscured star formation dominant thereafter. Having shown we can quantify the faint-end slope alpha of the LF accurately with our lensed HFF samples, we also quantify the apparent curvature in the shape of the UV LF through a curvature parameter delta. The constraints on the curvature delta strongly rule out the presence of a turn-over brightward of -13.1 mag at z~3, -14.3 mag at z~6, and -15.5 mag at all other redshifts between z~9 to z~2. △ Less

Submitted 1 April, 2023; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: 32 pages, 16 figures, 7 tables, ApJ, in press

High Equivalent Width of Hα+[N II] Emission in z~8 Lyman-break Galaxies from IRAC 5.8μm Observations: Evidence for Efficient Lyman-continuum Photon production in the Epoch of Re-ionization

Authors: Mauro Stefanon, Rychard J. Bouwens, Garth D. Illingworth, Ivo Labbé, Pascal A. Oesch, Valentino Gonzalez

Abstract: We measure, for the first time, the median equivalent width (EW) of H$α$+[N II] in star-forming galaxies at $z\sim8$. Our estimate leverages the unique photometric depth of the Spitzer/IRAC $5.8μ$m-band mosaics (probing $\approx 5500 - 7100$ A at $z\sim8$) of the GOODS Reionization Era Wide Area Treasury from Spitzer (GREATS) program. We median stacked the stamps of $102$ Lyman-break galaxies in t… ▽ More We measure, for the first time, the median equivalent width (EW) of H$α$+[N II] in star-forming galaxies at $z\sim8$. Our estimate leverages the unique photometric depth of the Spitzer/IRAC $5.8μ$m-band mosaics (probing $\approx 5500 - 7100$ A at $z\sim8$) of the GOODS Reionization Era Wide Area Treasury from Spitzer (GREATS) program. We median stacked the stamps of $102$ Lyman-break galaxies in the $3.6, 4.5, 5.8$ and $8.0μ$m bands, after carefully removing potential contamination from neighbouring sources. We infer an extreme rest-frame EW$_0$(H$α$+[N II])$=2328^{+1326}_{-1127}$ A from the measured red $[3.6]-[5.8]=0.82\pm0.27$ mag, consistent with young ($\lesssim10^7$ yr) average stellar population ages at $z\sim8$. This implies an ionizing photon production efficiency of $\log(ξ_{\mathrm{ion},0}/\mathrm{erg\ Hz}^{-1})=25.97^{+0.18}_{-0.28}$. Such a high value for photo production, similar to the highest values found at $z\lesssim4$, indicates that only modest escape fractions $f_\mathrm{esc}\lesssim0.3$ (at $2σ$) are sufficient for galaxies brighter than $M_\mathrm{UV}<-18$ mag to re-ionize the neutral Hydrogen at $z\sim8$. This requirement is relaxed even more to $f_\mathrm{esc}\le 0.1$ when considering galaxies brighter than $M_\mathrm{UV}\approx -13$ mag, consistent with recent luminosity functions and as typically assumed in studies addressing re-ionization. These exceptional results clearly indicate that galaxies can be the dominant source of reionizing photons, and provide us with an exciting glimpse into what we might soon learn about the early universe, and particularly about the Reionization Epoch, from forthcoming JWST/MIRI and NIRCam programs. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ. 12 pages, 5 figures (main) + appendix. Comments welcome

The ALMA REBELS Survey: Dust Continuum Detections at z > 6.5

Authors: Hanae Inami, Hiddo S. B. Algera, Sander Schouws, Laura Sommovigo, Rychard Bouwens, Renske Smit, Mauro Stefanon, Rebecca A. A. Bowler, Ryan Endsley, Andrea Ferrara, Pascal Oesch, Daniel Stark, Manuel Aravena, Laia Barrufet, Elisabete da Cunha, Pratika Dayal, Ilse De Looze, Yoshinobu Fudamoto, Valentino Gonzalez, Luca Graziani, Jacqueline A. Hodge, Alexander P. S. Hygate, Themiya Nanayakkara, Andrea Pallottini, Dominik A. Riechers , et al. (3 additional authors not shown)

Abstract: We report 18 dust continuum detections ($\geq 3.3σ$) at $\sim88{\rm μm}$ and $158{\rm μm}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as pa… ▽ More We report 18 dust continuum detections ($\geq 3.3σ$) at $\sim88{\rm μm}$ and $158{\rm μm}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as part of REBELS. In addition, 15 of the dust continuum detected galaxies also show a [CII]$_{\rm 158{\rm μm}}$ emission line, providing us with accurate redshifts. We anticipate more line emission detections from six targets (including three continuum detected targets) where observations are still ongoing. The dust continuum detected sources in our sample tend to have a redder UV spectral slope than the ones without a dust continuum detection. We estimate that all of the sources have an infrared (IR) luminosity ($L_{\rm IR}$) in a range of $3-8 \times 10^{11} L_\odot$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\,L_{\odot}$. Their fraction of obscured star formation is significant at $\gtrsim 50\%$. Some of the dust continuum detected galaxies show spatial offsets ($\sim 0.5-1.5''$) between the rest-UV and far-IR emission peaks. These separations appear to have an increasing trend against an indicator that suggests spatially decoupled phases of obscured and unobscured star formation. REBELS offers the best available statistical constraints on obscured star formation in UV-bright, massive galaxies at $z > 6.5$. △ Less

Submitted 26 September, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 17 pages, 9 figures, 4 tables, Accepted for publication in MNRAS

Journal ref: MNRAS, 515, 3126 (2022)

z~2-9 Galaxies magnified by the Hubble Frontier Field Clusters I: Source Selection and Surface Density-Magnification Constraints from >2500 galaxies

Authors: R. J. Bouwens, G. Illingworth, R. S. Ellis, P. Oesch, A. Paulino-Afonso, B. Ribeiro, M. Stefanon

Abstract: We assemble a large comprehensive sample of 2534 z~2, 3, 4, 5, 6, 7, 8, and 9 galaxies lensed by the six clusters from the Hubble Frontier Fields (HFF) program. Making use of the availability of multiple independent magnification models for each of the HFF clusters and alternatively treating one of the models as the "truth," we show that the median magnification factors from the v4 parametric mode… ▽ More We assemble a large comprehensive sample of 2534 z~2, 3, 4, 5, 6, 7, 8, and 9 galaxies lensed by the six clusters from the Hubble Frontier Fields (HFF) program. Making use of the availability of multiple independent magnification models for each of the HFF clusters and alternatively treating one of the models as the "truth," we show that the median magnification factors from the v4 parametric models are typically reliable to values of 30 to 50, and in one case to 100. Using the median magnification factor from the latest v4 models, we estimate the UV luminosities of the 2534 lensed z~2-9 galaxies, finding sources as faint as -12.4 mag at z~3 and -12.9 mag at z~7. We explicitly demonstrate the power of the surface density-magnification relations Sigma(z) vs. mu in the HFF clusters to constrain both distant galaxy properties and cluster lensing properties. Based on the Sigma(z) vs. mu relations, we show that the median magnification estimates from existing public models must be reliable predictors of the true magnification mu to mu<15 (95% confidence). We also use the observed Sigma(z) vs. mu relations to derive constraints on the evolution of the luminosity function faint-end slope from z~7 to z~2, showing that faint-end slope results can be consistent with blank-field studies if, and only if, the selection efficiency shows no strong dependence on the magnification factor mu. This can only be the case if very low luminosity galaxies are very small, being unresolved in deep lensing probes. △ Less

Submitted 1 May, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 22 pages, 13 figures, 6 tables, accepted for publication in the Astrophysical Journal, updated to be consistent with the proofs

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

Authors: Michael W. Topping, Daniel P. Stark, Ryan Endsley, Rychard J. Bouwens, Sander Schouws, Renske Smit, Mauro Stefanon, Hanae Inami, Rebecca A. A. Bowler, Pascal Oesch, Valentino Gonzalez, Pratika Dayal, Elisabete da Cunha, Hiddo Algera, Paul van der Werf, Andrea Pallottini, Laia Barrufet De Soto, Raffaella Schneider, Ilse De Looze, Laura Sommovigo, Lily Whitler, Luca Graziani, Yoshinobu Fudamoto, Andrea Ferrara

Abstract: We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we de… ▽ More We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 18 pages, 10 figures, submitted to MNRAS

The ALMA REBELS Survey: the dust content of $z \sim 7$ Lyman Break Galaxies

Authors: P. Dayal, A. Ferrara, L. Sommovigo, R. Bouwens, P. A. Oesch, R. Smit, V. Gonzalez, S. Schouws, M. Stefanon, C. Kobayashi, J. Bremer, H. S. B. Algera, M. Aravena, R. A. A. Bowler, E. da Cunha, Y. Fudamoto, L. Graziani, J. Hodge, H. Inami, I. De Looze, A. Pallottini, D. Riechers, R. Schneider, D. Stark, R. Endsley

Abstract: We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destructi… ▽ More We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $τ_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (ψ_{\rm UV}) = -0.05 ~[log (ψ)]^{2} + 0.86 ~log(ψ) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $τ_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted to MNRAS

The ALMA REBELS Survey. Epoch of Reionization giants: properties of dusty galaxies at $z \approx 7$

Authors: A. Ferrara, L. Sommovigo, P. Dayal, A. Pallottini, R. J. Bouwens, V. Gonzalez, H. Inami, R. Smit, R. A. A. Bowler, R. Endsley, P. Oesch, S. Schouws, D. Stark, M. Stefanon, M. Aravena, E. da Cunha, I. De Looze, Y. Fudamoto, L. Graziani, J. Hodge, D. Riechers, R. Schneider, H. S. B. Algera, L. Barrufet, A. P. S. Hygate , et al. (5 additional authors not shown)

Abstract: We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different d… ▽ More We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($β, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(β-β_{\rm int})$, where $β_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$β$ relation. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

ALMA as a Redshift Machine: Using [CII] to Efficiently Confirm Galaxies in the Epoch of Reionization

Authors: Sander Schouws, Rychard Bouwens, Renske Smit, Jacqueline Hodge, Mauro Stefanon, Joris Witstok, Juliëtte Hilhorst, Ivo Labbe, Hiddo Algera, Leindert Boogaard, Michael Maseda, Pascal Oesch, Huub Röttgering, Paul van der Werf

Abstract: The [CII]$_{158μm}$ line has long been proposed as a promising line to spectroscopically confirm galaxies in the epoch of reionization. In this paper we present the results of new ALMA observations spectral scanning for [CII] in six particularly luminous Lyman Break Galaxies at $z\sim7$. The six sources were drawn from a sample of bright $z\sim7$ galaxies identified using the wide-area optical, ne… ▽ More The [CII]$_{158μm}$ line has long been proposed as a promising line to spectroscopically confirm galaxies in the epoch of reionization. In this paper we present the results of new ALMA observations spectral scanning for [CII] in six particularly luminous Lyman Break Galaxies at $z\sim7$. The six sources were drawn from a sample of bright $z\sim7$ galaxies identified using the wide-area optical, near-IR, and Spitzer/IRAC data over the COSMOS/UltraVISTA field and were targeted on the basis of tight constraints on their redshifts from their IRAC [3.6]-[4.5] colors. We detect significant ($>9σ$) [CII] lines in three of our six targets ($50\%$) co-spatial with the rest-$UV$ emission from the ground/space-based near-IR imaging. The luminosities of the [CII] lines lie in the range $5.6$ to $8.8\times10^{8}L_{\odot}$, consistent with the local [CII]-SFR relation. Meanwhile, their [CII]/$L_{IR}\sim1-3\times10^{-3}$ ratios are slightly elevated compared to local (U)LIRGS. This could be due to lower dust-to-gas or dust-to-metal ratios. We also find that our sources display a large kinematic diversity, with one source showing signs of rotation, one source a likely major merger and one dispersion dominated source that might contain a bright star-forming clump. Our results highlight the effectiveness of spectral scans with ALMA in spectroscopically confirming luminous galaxies in the epoch of reionization, something that is being be applied on a significantly larger sample in the on-going REBELS large program. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 16 pages, 7 figures, 3 tables. Submitted to ApJ

The REBELS ALMA Survey: cosmic dust temperature evolution out to z $\sim$ 7

Authors: L. Sommovigo, A. Ferrara, A. Pallottini, P. Dayal, R. J. Bouwens, R. Smit, E. da Cunha, I. De Looze, R. A. A. Bowler, J. Hodge, H. Inami, P. Oesch, R. Endsley, V. Gonzalez, S. Schouws, D. Stark, M. Stefanon, M. Aravena, L. Graziani, D. Riechers, R. Schneider, P. van der Werf, H. Algera, L. Barrufet, Y. Fudamoto , et al. (5 additional authors not shown)

Abstract: ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star format… ▽ More ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$μ$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

The ALMA REBELS Survey: Efficient Ly$α$ Transmission of UV-Bright z$\simeq$7 Galaxies from Large Velocity Offsets and Broad Line Widths

Authors: Ryan Endsley, Daniel P. Stark, Rychard J. Bouwens, Sander Schouws, Renske Smit, Mauro Stefanon, Hanae Inami, Rebecca A. A. Bowler, Pascal Oesch, Valentino Gonzalez, Manuel Aravena, Elisabete da Cunha, Pratika Dayal, Andrea Ferrara, Luca Graziani, Themiya Nanayakkara, Andrea Pallottini, Raffaella Schneider, Laura Sommovigo, Michael Topping, Paul van der Werf, Anne Hutter

Abstract: Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$α$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galax… ▽ More Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$α$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly$α$ profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at $z\sim7$ remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [CII]-based redshifts from ALMA. Here we report MMT/Binospec Ly$α$ spectroscopy of eight UV-bright ($\mathrm{M_{UV}}^{}\sim-22$) galaxies at $z\simeq7$ selected from the ALMA REBELS survey. We detect Ly$α$ in 4 of 8 galaxies and use the [CII] systemic redshifts to investigate the Ly$α$ velocity profiles. The Ly$α$ lines are significantly redshifted from systemic (average velocity offset=223 km/s) and broad (FWHM$\approx$300$-$650 km/s), with two sources showing emission extending to $\approx$750 km/s. We find that the broadest Ly$α$ profiles are associated with the largest [CII] line widths, suggesting a potential link between the Ly$α$ FWHM and the dynamical mass. Since Ly$α$ photons at high velocities transmit efficiently through the $z=7$ IGM, our data suggest that velocity profiles play a significant role in boosting the Ly$α$ visibility of the most UV-luminous reionization-era galaxies. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 17 pages, 7 figures, submitted to MNRAS

ALMA measures molecular gas reservoirs comparable to field galaxies in a low-mass galaxy cluster at z=1.3

Authors: Christina C. Williams, Stacey Alberts, Justin S. Spilker, Allison G. Noble, Mauro Stefanon, Christopher N. A. Willmer, Rachel Bezanson, Desika Narayanan, Katherine E. Whitaker

Abstract: We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas… ▽ More We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas measurements, bridging the observational gap between galaxies in the more extreme, well-studied clusters (Log$_{10}$~M$_{\rm halo}/M_{\odot}\gtrsim14$) and those in group or field environments at cosmic noon. The CO sources are concentrated on the sky (within ~1-arcmin diameter) and phase space analysis indicates the gas resides in galaxies already within the cluster environment. We find that CO sources sit in similar phase space as CO-rich galaxies in more massive clusters at similar redshifts (have similar accretion histories) while maintaining field-like molecular gas reservoirs, compared to scaling relations. This work presents the deepest CO survey to date in a galaxy cluster at z>1, uncovering gas reservoirs down to M$_{\rm H_{2}}>1.6\times10^{10}$M$_{\odot}$ (5$σ$ at 50% primary beam). Our deep limits rule out the presence of gas content in excess of the field scaling relations; however, combined with literature CO detections, cluster gas fractions in general appear systematically high, on the upper envelope or above the field. This study is the first demonstration that low mass clusters at z~1-2 can host overdensities of CO emitters with surviving gas reservoirs, in line with the prediction that quenching is delayed after first infall while galaxies consume the gas bound to the disk. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: Submitted to AAS Journals on December 15, 2021. Comments welcome

The Prevalence of Galaxy Overdensities Around UV-Luminous Lyman $\mathbfα$ Emitters in the Epoch of Reionization

Authors: E. Leonova, P. A. Oesch, Y. Qin, R. P. Naidu, J. S. B. Wyithe, S. de Barros, R. J. Bouwens, R. S. Ellis, R. M. Endsley, A. Hutter, G. D. Illingworth, J. Kerutt, I. Labbe, N. Laporte, D. Magee, S. J. Mutch, G. W. Roberts-Borsani, R. Smit, D. P. Stark, M. Stefanon, S. Tacchella, A. Zitrin

Abstract: Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$α$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$α$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the… ▽ More Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$α$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$α$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the Ly$α$ transmission of such sources. We present an analysis of dedicated Hubble Space Telescope (HST) imaging in the CANDELS/EGS field to search for fainter neighbours around three of the most UV luminous and most distant spectroscopically confirmed Ly$α$ emitters: EGS-zs8-1, EGS-zs8-2 and EGSY-z8p7 at $z_\mathrm{spec}=7.73$, 7.48, and 8.68, respectively. We combine the multi-wavelength HST imaging with Spitzer data to reliably select $z\sim7-9$ galaxies around the central, UV-luminous sources. In all cases, we find a clear enhancement of neighbouring galaxies compared to the expected number in a blank field (by a factor $\sim 3-9\times$). Our analysis thus reveals ubiquitous overdensities around luminous Ly$α$ emitting sources in the heart of the cosmic reionization epoch. We show that our results are in excellent agreement with expectations from the Dragons simulation, confirming the theoretical prediction that the first ionized bubbles preferentially formed in overdense regions. JWST follow-up observations of the neighbouring galaxies identified here will be needed to confirm their physical association and to map out the ionized regions produced by these sources. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 13 pages, 14 figures, Submitted to MNRAS

Sizes of Lensed Lower-luminosity z=4-8 Galaxies from the Hubble Frontier Field Program

Authors: R. J. Bouwens, G. D. Illingworth, P. G. van Dokkum, P. A. Oesch, M. Stefanon, B. Ribeiro

Abstract: We constrain the rest-UV size-luminosity relation for star-forming galaxies at z~4 and z~6, 7, and 8 identified behind clusters from the Hubble Frontier Fields (HFF) program. The size-luminosity relation is key to deriving accurate luminosity functions (LF) for faint galaxies. Making use of the latest lensing models and full data set for these clusters, lensing-corrected sizes and luminosities are… ▽ More We constrain the rest-UV size-luminosity relation for star-forming galaxies at z~4 and z~6, 7, and 8 identified behind clusters from the Hubble Frontier Fields (HFF) program. The size-luminosity relation is key to deriving accurate luminosity functions (LF) for faint galaxies. Making use of the latest lensing models and full data set for these clusters, lensing-corrected sizes and luminosities are derived for 68 z~4, 184 z~6, 93 z~7, and 53 z~8 galaxies. We show that size measurements can be reliably measured up to linear magnifications of 30x, where the lensing models are well calibrated. The sizes we measure span a >1-dex range, from <50 pc to >~500 pc. Uncertainties are based on both the formal fit errors and systematic differences between the public lensing models. These uncertainties range from ~20 pc for the smallest sources to 50 pc for the largest. Using a forward-modeling procedure to model the impact of incompleteness and magnification uncertainties, we characterize the size-luminosity relation at both z~4 and z~6-8. We find that the source sizes of star-forming galaxies at z~4 and z~6-8 scale with luminosity L as L^{0.54\pm0.08} and L^{0.40+/-0.04}, respectively, such that lower luminosity (>~-18 mag) galaxies are smaller than expected from extrapolating the size-luminosity relation at high luminosities (<~-18 mag). The new evidence for a steeper size-luminosity relation (3 sigma) adds to earlier evidence for small sizes based on the prevalence of highly magnified galaxies in high shear regions, theoretical arguments against upturns in the LFs, and other independent determinations of the size-luminosity relation from the HFF clusters. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Comments: 22 pages, 18 figures, 2 tables, accepted for publication into ApJ. arXiv admin note: substantial text overlap with arXiv:1711.02090