skip to main content

Title: The ALMA REBELS survey: the dust-obscured cosmic star formation rate density at redshift 7
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-)millimetre 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 [C ii] and dust emission in UV-selected galaxies at z ∼ 7, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation fobs = SFRIR/SFRUV+IR depends on stellar mass, similar to what is observed at lower redshift, and ranges from fobs ≈ 0.3 − 0.6 for galaxies with log10(M⋆/M⊙) = 9.4–10.4. We further adopt the z ∼ 7 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 ≲ z ≲ 7, with dust-obscured star formation still contributing ${\sim}30{{\ \rm per\ cent}}$ at z ∼ 7. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into more » the epoch of reionization.

« less
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; « less
Publication Date:
NSF-PAR ID:
10385765
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
518
Issue:
4
Page Range or eLocation-ID:
p. 6142-6157
ISSN:
0035-8711
Publisher:
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We present specific star formation rates (sSFRs) for 40 ultraviolet (UV)-bright galaxies at z ∼ 7–8 observed as part of the Reionization Era Bright Emission Line Survey (REBELS) Atacama Large Millimeter/submillimeter Array (ALMA) large programme. The sSFRs are derived using improved star formation rate (SFR) calibrations and spectral energy distribution (SED)-based stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [C ii]-based spectroscopic redshifts. The median sSFR of the sample is $18_{-5}^{+7}$ Gyr−1, significantly larger than literature measurements lacking constraints in the FIR, reflecting the larger obscured SFRs derived 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 in reionization-era galaxies. When large sSFR galaxies (a population that is common at z > 6) are modelled with non-parametric star formation histories, 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 thatmore »is outshined by the recent burst. The [C ii] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses capable of accommodating an old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties among derived parameters, we find that sSFRs increase rapidly toward higher redshifts for massive galaxies (9.6 < log (M*/M⊙) < 9.8), evolving as (1 + z)1.7 ± 0.3, broadly consistent with expectations from the evolving baryon accretion rates.

    « less
  2. ABSTRACT

    We post-process galaxies in the IllustrisTNG simulations with skirt radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at z ≥ 4. The rest-frame K- and z-band galaxy luminosity functions from TNG are overall consistent with observations, despite ${\sim}0.5\, \mathrm{dex}$ underprediction at z = 4 for MK ≲ −25 and Mz ≲ −24. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500\, {\rm arcmin}^{2}$ at z = 6 (z = 8). As opposed to the consistency in the UV, optical, and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshiftmore »are qualitatively consistent with observations, the peak dust temperature of z ≥ 6 galaxies are overestimated by about $20\, {\rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.

    « less
  3. Abstract We make use of sensitive (9.3 μ Jy beam −1 rms) 1.2 mm continuum observations from the Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra-Deep Field (ASPECS) large program to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z  = 1.5–10 (to ∼7–28 M ⊙ yr −1 at 4 σ over the entire range). We find that the fraction of ALMA-detected galaxies in our z  = 1.5–10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 10 9.0 M ⊙ to % at >10 10 M ⊙ . Moreover, on stacking all 1253 low-mass (<10 9.25 M ⊙ ) galaxies over the ASPECS footprint, we find a mean continuum flux of −0.1 ± 0.4 μ Jy beam −1 , implying a hard upper limit on the obscured star formation rate of <0.6 M ⊙ yr −1 (4 σ ) in a typical low-mass galaxy. The correlation between the infrared excess (IRX) of UV-selected galaxies ( L IR / L UV ) and the UV-continuum slope is also seen in our ASPECS data and shows consistency with a Calzetti-like relation at > and an SMC-like relation at lower masses. Using stellar massmore »and β measurements for z  ∼ 2 galaxies over the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, we derive a new empirical relation between β and stellar mass and then use this correlation to show that our IRX– β and IRX–stellar mass relations are consistent with each other. We then use these constraints to express the IRX as a bivariate function of β and stellar mass. Finally, we present updated estimates of star formation rate density determinations at z  > 3, leveraging present improvements in the measured IRX and recent probes of ultraluminous far-IR galaxies at z  > 2.« less
  4. Aims. We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE). These sources, detected in COSMOS and ECDFS, have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z  ≃ 6. Methods. We looked for counterparts of the ALMA sources in all the available multi-wavelength (from HST to VLA) and photometric redshift catalogues. We also made use of deeper UltraVISTA and Spitzer source lists and maps to identify optically dark sources with no matches in the public catalogues. We used the sources with estimated redshifts to derive the 250 μ m rest-frame and total infrared (8–1000 μ m) LFs from z  ≃ 0.5 to 6. Results. Our ALMA blind survey (860 μ m flux density range: ∼0.3–12.5 mJy) allows us to further push the study of the nature and evolution of dusty galaxies at high- z , identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint end of themore »infrared LF, showing little evolution from z  ≃ 2.5 to z ≃ 6, and a “flat” slope up to the highest redshifts (i.e. 4.5 <   z  <  6). The SFRD obtained by integrating the luminosity function remains almost constant between z  ≃ 2 and z  ≃ 6, and significantly higher than the optical or ultra-violet derivations, showing a significant contribution of dusty galaxies and obscured star formation at high- z . About 14% of all the ALPINE serendipitous continuum sources are found to be optically and near-infrared (near-IR) dark (to a depth K s  ∼ 24.9 mag). Six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [C II] emitters at z  ≃ 5. The six HST+near-IR dark galaxies with mid-IR counterparts are found to contribute about 17% of the total SFRD at z  ≃ 5 and to dominate the high-mass end of the stellar mass function at z  >  3.« less
  5. ABSTRACT

    We present 65 extremely dust-obscured galaxies from the UltraVISTA DR3 survey of the COSMOS field at 1 < z < 4. In contrast to other studies of dusty galaxies, we select our sample based on dust attenuation measured by UV–MIR spectral energy distribution (SED) modelling that allows for extreme attenuation levels. We construct our sample by making cuts at 1 ≤ z ≤ 4, AV ≥ 3, and log(M*/M⊙) ≥ 10.5. This method reliably selects galaxies exhibiting independent indicators of significant dust content, including far-infrared detection rates. We perform panchromatic SED modelling with matched Herschel photometry and find stellar and dust properties that differ from typical submillimetre galaxy (SMG) samples as well as Herschel sources matched in redshift and stellar mass. Our sources have lower star formation rates and higher AV than SMGs, but comparable total infrared luminosities. Most of our sample falls on or near the star-forming main sequence for this redshift range. Finally, we perform a morphological analysis with galfit using the KS-band images and Hubble F814W and F160W imaging when available. Typical axial ratios of ∼0.4 suggest disc-like morphology for the majority of our sources, and we note only three apparent merging systems. Our sample generallymore »agrees with the size–mass relation for star-forming galaxies, with a tail extending to smaller sizes. We conclude that the most heavily obscured galaxies in this redshift range share many characteristics with typical star-forming galaxies, forming a population of dusty galaxies that overlaps, but is not encompassed by, those selected through dust emission.

    « less