Search for: All records

Abstract The H i gas content is a key ingredient in galaxy evolution, the study of which has been limited to moderate cosmological distances for individual galaxies due to the weakness of the hyperfine H i 21 cm transition. Here we present a new approach that allows us to infer the H i gas mass M HI of individual galaxies up to z ≈ 6, based on a direct measurement of the [C ii ]-to-H i conversion factor in star-forming galaxies at z ≳ 2 using γ -ray burst afterglows. By compiling recent [C ii ]-158 μ m emission line measurements we quantify the evolution of the H i content in galaxies through cosmic time. We find that M HI starts to exceed the stellar mass M ⋆ at z ≳ 1, and increases as a function of redshift. The H i fraction of the total baryonic mass increases from around 20% at z = 0 to about 60% at z ∼ 6. We further uncover a universal relation between the H i gas fraction M HI / M ⋆ and the gas-phase metallicity, which seems to hold from z ≈ 6 to z = 0. The majority of galaxiesmore »at z > 2 are observed to have H i depletion times, t dep,HI = M HI /SFR, less than ≈2 Gyr, substantially shorter than for z ∼ 0 galaxies. Finally, we use the [C ii ]-to-H i conversion factor to determine the cosmic mass density of H i in galaxies, ρ HI , at three distinct epochs: z ≈ 0, z ≈ 2, and z ∼ 4–6. These measurements are consistent with previous estimates based on 21 cm H i observations in the local universe and with damped Ly α absorbers (DLAs) at z ≳ 2, suggesting an overall decrease by a factor of ≈5 in ρ HI ( z ) from the end of the reionization epoch to the present.« less

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 ∼ 7 remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [C ii]-based redshifts from ALMA. Here, we report MMT/Binospec Ly α spectroscopy of eight UV-bright (MUV ∼ −22) galaxies at z ≃ 7 selected from the ALMA REBELS survey. We detect Ly α in four of eight galaxies and use the [C ii] systemic redshifts to investigate the Ly α velocity profiles. The Ly α lines are significantly redshifted from systemic (average velocity offset = 223 km s–1) and broad (FWHM ≈ 300–650 km s−1), with two sources showing emission extending to ≈750 km s−1. We find that the broadest Ly α profiles aremore »associated with the largest [C ii] 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

Abstract We present results from Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm continuum observations of a sample of 27 star-forming galaxies at z = 2.1–2.5 from the MOSFIRE Deep Evolution Field survey with metallicity and star formation rate measurements from optical emission lines. Using stacks of Spitzer, Herschel, and ALMA photometry (rest frame ∼8–400 μ m), we examine the infrared (IR) spectral energy distributions (SED) of z ∼ 2.3 subsolar-metallicity (∼0.5 Z ⊙ ) luminous infrared galaxies (LIRGs). We find that the data agree well with an average template of higher-luminosity local low-metallicity dwarf galaxies (reduced χ 2 = 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ultraluminous IR galaxies, even in the most favorable case (with reduced χ 2 = 2.8), the templates are rejected at >98% confidence. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar-metallicity galaxies may result from different grain properties or a harder/more intense ionizing radiation field that increases the dust temperature. The obscured star formation rate (SFR) indicated by the far-IR emission of the subsolar-metallicity galaxies is only ∼60% of the total SFR, considerably lower than that of the local LIRGsmore »with ∼96%–97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data overestimate the Rayleigh–Jeans tail measurements by a factor of 2–20. These templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4 dex. At a given stellar mass or metallicity, dust masses at z ∼ 2.3 are an order of magnitude higher than z ∼ 0. Given the predicted molecular gas fractions, the observed z ∼ 2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies with similar metallicities.« less

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

« less

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

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here, we introduce an emissivity framework based on Lyα emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between production and escape of Lyα and LyC photons. Using constraints on the escape fractions of bright LAEs (LLyα > 0.2L*) at z ≈ 2 obtained from resolved Lyα profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at z ≈ 2–6, which is non-trivial given the strong evolution in both the star formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between z ≈ 6−9 under plausible assumptions. Within this framework, the >10 × rise in the UV population-averaged fesc between z ≈ 3–7 naturally arises due to the same phenomena that drive the growing LAE fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm {esc}}=50{{\ \rm per\ cent}}$, ξion = 1025.9 Hz erg−1 formore »half the bright LAEs), a highly ionizing minority of galaxies with MUV < −17 accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a ‘disco’ ionizing background. We conclude proposing tests to further develop our suggested Lyα-anchored formalism.

« less

The ionizing photon escape fraction [Lyman continuum (LyC) fesc] of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshifts. Here we circumvent this uncertainty by inferring fesc using resolved (R > 4000) Lyman α (Lyα) profiles from the X-SHOOTER Lyα survey at z = 2 (XLS-z2). With empirically motivated criteria, we use Lyα profiles to select leakers ($f_{\mathrm{ esc}} > 20{{\ \rm per\ cent}}$) and non-leakers ($f_{\mathrm{ esc}} < 5{{\ \rm per\ cent}}$) from a representative sample of >0.2L* Lyman α emitters (LAEs). We use median stacked spectra of these subsets over λrest ≈ 1000–8000 Å to investigate the conditions for LyC fesc. Our stacks show similar mass, metallicity, MUV, and βUV. We find the following differences between leakers versus non-leakers: (i) strong nebular C iv and He ii emission versus non-detections; (ii) [O iii]/[O ii] ≈ 8.5 versus ≈3; (iii) Hα/Hβ indicating no dust versus E(B − V) ≈ 0.3; (iv) Mg ii emission close to the systemic velocity versus redshifted, optically thick Mg ii; and (v) Lyα fesc of ${\approx} 50{{\ \rm per\ cent}}$ versus ${\approx} 10{{\ \rm per\ cent}}$. The extreme equivalent widths (EWs) in leakers ([O iii]+$\mathrm{ H}\beta \approx 1100$ Å rest frame)more »constrain the characteristic time-scale of LyC escape to ≈3–10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers – extremely ionizing stellar populations, low column densities, a dust-free, high-ionization state interstellar medium (ISM) – occur simultaneously in the $f_{\rm esc} > 20{{\ \rm per\ cent}}$ stack, suggesting they are causally connected, and motivating why indicators like [O iii]/[O ii] may suffice to constrain fesc at z > 6 with the James Webb Space Telescope (JWST). The leakers comprise half of our sample, have a median LyC$f_{\rm esc} \approx 50{{\ \rm per\ cent}}$ (conservative range: $20\!-\!55{{\ \rm per\ cent}}$), and an ionizing production efficiency $\log ({\xi _{\rm {ion}}/\rm {Hz\ erg^{-1}}})\approx 25.9$ (conservative range: 25.7–25.9). These results show LAEs – the type of galaxies rare at z ≈ 2, but that become the norm at higher redshift – are highly efficient ionizers, with extreme ξion and prolific fesc occurring in sync.

« less

ABSTRACT The infrared (IR) spectral energy distributions (SEDs) of main-sequence galaxies in the early Universe (z > 4) is currently unconstrained as IR continuum observations are time-consuming and not feasible for large samples. We present Atacama Large Millimetre Array Band 8 observations of four main-sequence galaxies at z ∼ 5.5 to study their IR SED shape in detail. Our continuum data (rest-frame 110 $\rm \mu m$, close to the peak of IR emission) allows us to constrain luminosity-weighted dust temperatures and total IR luminosities. With data at longer wavelengths, we measure for the first time the emissivity index at these redshifts to provide more robust estimates of molecular gas masses based on dust continuum. The Band 8 observations of three out of four galaxies can only be reconciled with optically thin emission redward of rest-frame $100\, {\rm \mu m}$. The derived dust peak temperatures at z ∼ 5.5 ($30\!-\!43\, {\rm K}$) are elevated compared to average local galaxies, however, $\sim 10\, {\rm K}$ below what would be predicted from an extrapolation of the trend at z < 4. This behaviour can be explained by decreasing dust abundance (or density) towards high redshifts, which would cause the IR SED at the peakmore »to be more optically thin, making hot dust more visible to the external observer. From the $850{\hbox{-}}{\rm \mu m}$ dust continuum, we derive molecular gas masses between 1010 and $10^{11}\, {\rm M_{\odot }}$ and gas fractions (gas over total mass) of $30\!-\!80{{\ \rm per\ cent}}$ (gas depletion times of $100\!-\!220\, {\rm Myr}$). All in all, our results provide a first measured benchmark SED to interpret future millimetre observations of normal, main-sequence galaxies in the early Universe.« less

Abstract The Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Band 6 scan (212–272 GHz) covers potential [C ii ] emission in galaxies at 6 ≤ z ≤ 8 throughout a 2.9 arcmin 2 area. By selecting on known Ly α emitters (LAEs) and photometric dropout galaxies in the field, we perform targeted searches down to a 5 σ [C ii ] luminosity depth L [C II ] ∼ 2.0 × 10 8 L ⊙ , corresponding roughly to star formation rates (SFRs) of 10–20 M ⊙ yr −1 when applying a locally calibrated conversion for star-forming galaxies, yielding zero detections. While the majority of galaxies in this sample are characterized by lower SFRs, the resulting upper limits on [C ii ] luminosity in these sources are consistent with the current literature sample of targeted ALMA observations of z = 6–7 LAEs and Lyman-break galaxies (LBGs), as well as the locally calibrated relations between L [C ii ] and SFR—with the exception of a single [C ii ]-deficient, UV-luminous LBG. We also perform a blind search for [C ii ]-bright galaxies that may have been missed by optical selections, resulting in an upper limitmore »on the cumulative number density of [C ii ] sources with L [C II ] > 2.0 × 10 8 L ⊙ (5 σ ) to be less than 1.8 × 10 −4 Mpc −3 (90% confidence level). At this luminosity depth and volume coverage, we present an observed evolution of the [C ii ] luminosity function from z = 6–8 to z ∼ 0 by comparing the ASPECS measurement to literature results at lower redshift.« less