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 the epoch of reionization.

 
more » « less
NSF-PAR ID:
10385765
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; « less
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
518
Issue:
4
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 6142-6157
Size(s):
["p. 6142-6157"]
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We investigate the degree of dust obscured star formation in 49 massive (log10(M⋆/M⊙) > 9) Lyman-break galaxies (LBGs) at z = 6.5–8 observed as part of the Atacama Large Millimeter/submillimeter Array (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 ultraviolet (UV), optical, and far-infrared (FIR) properties which reveal a significant fraction (fobs = 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 (β ≃ −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 βd = 2.0), we find that the REBELS sources are in agreement with the local ‘Calzetti-like’ starburst Infrared-excess (IRX)–β relation. By re-analysing the data available for 108 galaxies at z ≃ 4–6 from the ALMA Large Program to Investigate C+ at Early Times (ALPINE) using a consistent methodology and assumed FIR SED, we show that from z ≃ 4–8, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX–β relation. When comparing the IRX–M⋆ 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 ≳ 3 at a given a M⋆. Our IRX–β results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for z ≃ 7 galaxies with similar stellar masses and star formation rates.

     
    more » « less
  2. 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 that 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.

     
    more » « less
  3. 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 redshift 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.

     
    more » « 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 the 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. 
    more » « less
  5. 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 mass 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. 
    more » « less