We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the Atacama Large Millimeter/submillimeter Array (ALMA) public archive in the ...COSMOS deep field (A3COSMOS). Our A3COSMOS data set contains ∼700 galaxies (0.3 z 6) with high-confidence ALMA detections in the (sub)millimeter continuum and multiwavelength spectral energy distributions. Multiple gas mass calibration methods are compared, and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans tail continuum) have been tested. Combining our A3COSMOS sample with ∼1000 CO-observed galaxies at 0 z 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time ( ) and molecular gas to stellar mass ratio ( ) each as a function of the stellar mass ( ), offset from the star-forming main sequence ( ) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature) and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semianalytic modeling. These together provide a coherent picture of cold molecular gas, star formation rate, and stellar mass evolution in galaxies across cosmic time.
The Atacama Large Millimeter Array (ALMA) Large Program to INvestigate CII at Early times (ALPINE) targets the CII 158
μ
m line and the far-infrared continuum in 118 spectroscopically confirmed ...star-forming galaxies between
z
= 4.4 and
z
= 5.9. It represents the first large CII statistical sample built in this redshift range. We present details regarding the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates (SFRs), we measured the conversion factor from the ALMA 158
μ
m rest-frame dust continuum luminosity to the total infrared luminosity (
L
IR
) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median
L
IR
of 4.4 × 10
11
L
⊙
. We also detected 57 additional continuum sources in our ALMA pointings. They are at a lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5 ± 0.2. We measured the 850
μ
m number counts between 0.35 and 3.5 mJy, thus improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40% of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the CII line in 75 of our targets. Their median CII luminosity is 4.8 × 10
8
L
⊙
and their median full width at half maximum is 252 km s
−1
. After measuring the mean obscured SFR in various CII luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR–
L
CII
relations.
Abstract
The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous ...star-forming galaxies known in the
z
> 6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7 deg
2
area for bright C
ii
158
μ
m
and O
iii
88
μ
m
lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at
z
> 6.5 by ∼4–5× to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the first-year observations, 18 highly significant ≥ 7
σ
C
ii
158
μ
m
lines have already been discovered, the bulk of which (13/18) also show ≥3.3
σ
dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the
z
> 6.5 universe, such that the number of ALMA-derived redshifts at
z
> 6.5 rival Ly
α
discoveries. An analysis of the completeness of our search results versus star formation rate (SFR) suggests an ∼79% efficiency in scanning for C
ii
158
μ
m
when the SFR
UV+IR
is >28
M
⊙
yr
−1
. These new LP results further demonstrate ALMA’s efficiency as a “redshift machine,” particularly in the Epoch of Reionization.
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 ≃ 7. We find that the galaxy dust mass, Md, is regulated by the combination of Type II supernova dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard time-scale τ0 = 30 Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07\!-\!0.1{{\ \rm per\ cent}}$ and a UV-to-total star formation rate relation such that log(ψUV) = −0.05 log(ψ)2 + 0.86 log(ψ) − 0.05 (implying that 55–80 per cent of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9}\!-\!10^{10} \rm 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{{\ \rm per\ cent}}$) 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 of ∼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.
ABSTRACT
ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature Td 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 C $\scriptstyle \rm II$ 158-μm line and underlying dust continuum measurements, we derive Td in the continuum and C $\scriptstyle \rm II$ detected $z$ ≈ 7 galaxies in the ALMA Large Project REBELS sample. We find 39 < Td < 58 K, and dust masses in the narrow range Md = (0.9−3.6) × 107 M⊙. These results allow us to extend for the first time the reported Td($z$) relation into the Epoch of Reionization. We produce a new physical model that explains the increasing Td($z$) trend with the decrease of gas depletion time, tdep = Mg/SFR, induced by the higher cosmological accretion rate at early times; this hypothesis yields Td ∝ (1 + $z$)0.4. The model also explains the observed Td 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, Td only depends on the gas column density (metallicity), $T_{\rm d} \propto N_{\rm H}^{1/6}$ (Td ∝ Z−1/6). REBELS galaxies are on average relatively transparent, with effective gas column densities around NH ≃ (0.03−1) × 1021 cm−2. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated Td ≫ 60 K, are significantly obscured, low-metallicity systems. In fact, Td is higher in metal-poor systems due to their smaller dust content, which for fixed LIR results in warmer temperatures.
We study the efficiency of galactic feedback in the early Universe by stacking the C II 158
μ
m emission in a large sample of normal star-forming galaxies at 4 <
z
< 6 from the ALMA Large Program ...to INvestigate C II at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked C II profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked C II spectrum, with velocities of |
v
|≲500 km s
−1
. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFR
med
= 25
M
⊙
yr
−1
), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at
z
> 4. From the stacking analysis of the datacubes, we find that the combined C II core emission (|
v
|< 200 km s
−1
) of the higher-SFR galaxies is extended on physical sizes of ∼30 kpc (diameter scale), well beyond the analogous C II core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.
ABSTRACT
We present a method that derives the dust temperatures and infrared (IR) luminosities of high-redshift galaxies assuming radiation equilibrium in a simple dust and stellar distribution ...geometry. Using public data from the Atacama Large Millimeter/submillimeter Array (ALMA) archive, we studied dust temperatures assuming a clumpy interstellar medium (ISM) model for high-redshift galaxies, then tested the consistency of our results with those obtained using other methods. We find that a dust distribution model assuming a clumpiness of ${\rm log}\, \xi _{\rm clp}=-1.02\pm 0.41$ may accurately represent the ISM of high-redshift star-forming galaxies. By assuming a value of ξclp, our method enables the derivation of dust temperatures and IR luminosities of high-redshift galaxies from dust continuum fluxes and emission sizes obtained from single-band ALMA observations. To demonstrate the method proposed herein, we determined the dust temperature ($T_{\rm d}=95^{+13}_{-17}\, \rm {K}$) of a z ∼ 8.3 star-forming galaxy, MACS0416-Y1. Because the method only requires a single-band dust observation to derive a dust temperature, it is more easily accessible than multiband observations or high-redshift emission line searches and can be applied to large samples of galaxies in future studies using high-resolution interferometers such as ALMA.
ABSTRACT
We analyse FIR dust continuum measurements for 14 galaxies (redshift z ≈ 7) in the ALMA Reionization Era Bright Emission Line Survey (REBELS) Large Program to derive their physical ...properties. Our model uses three input data, i.e. (a) the UV spectral slope, β, (b) the observed UV continuum flux at 1500 Å, F1500, (c) the observed continuum flux at $\approx 158\, \mu$m, F158, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have 28−90.5 per cent of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 \lt {\rm SFR}/({\rm M}_\odot \, {\rm yr}^{-1}) \lt 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 \, \mathrm{M}_\odot$ and 52 ± 11 K, respectively. However, in some galaxies dust is particularly abundant (REBELS-14, $M^{\prime }_{\rm d} \approx 3.4 \times 10^7 \, \mathrm{M}_\odot$), or hot (REBELS-18, $T^{\prime }_{\rm d} \approx 67$ K). The dust distribution is compact (<0.3 kpc for 70 per cent of the galaxies). The inferred dust yield per supernova is $0.1 \le y_{\rm d}/\, \mathrm{M}_\odot \le 3.3$, with 70 per cent of the galaxies requiring $y_{\rm d} \lt 0.25 \, \mathrm{M}_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_{\rm d} \gt 1 \, \mathrm{M}_\odot$, which is likely inconsistent with pure SN production, and might require dust growth via accretion of heavy elements from the interstellar medium. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in C ii nicely follow the local LCII−SFR relation, and are approximately located on the Kennicutt–Schmidt relation. The sample-averaged gas depletion time is $0.11\, y_{\rm P}^{-2}$ Gyr, where yP is the ratio of the gas-to-stellar distribution radius. For some systems, a solution simultaneously matching the observed (β, F1500, F158) values cannot be found. This occurs when the index Im = (F158/F1500)/(β − βint), where βint is the intrinsic UV slope, exceeds $I_m^{*}\approx 1120$ for an 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.
ABSTRACT
While the kinematics of galaxies up to z ∼ 3 have been characterized in detail, only a handful of galaxies at high redshift (z > 4) have been examined in such a way. The Atacama Large ...Millimeter/submillimeter Array (ALMA) Large Program to INvestigate C ii at Early times (ALPINE) survey observed a statistically significant sample of 118 star-forming main-sequence galaxies at z = 4.4–5.9 in C ii158 $\mu$m emission, increasing the number of such observations by nearly 10×. A preliminary qualitative classification of these sources revealed a diversity of kinematic types (i.e. rotators, mergers, and dispersion-dominated systems). In this work, we supplement the initial classification by applying quantitative analyses to the ALPINE data: a tilted ring model (TRM) fitting code (3Dbarolo), a morphological classification (Gini-M20), and a set of disc identification criteria. Of the 75 C ii-detected ALPINE galaxies, 29 are detected at sufficient significance and spatial resolution to allow for TRM fitting and the derivation of morphological and kinematic parameters. These 29 sources constitute a high-mass subset of the ALPINE sample ($M_*\gt 10^{9.5}\, \mathrm{M}_{\odot }$). We robustly classify 14 of these sources (six rotators, five mergers, and three dispersion-dominated systems); the remaining sources showing complex behaviour. By exploring the G-M20 of z > 4 rest-frame far-infrared and C ii data for the first time, we find that our 1 arcsec ∼ 6 kpc resolution data alone are insufficient to separate galaxy types. We compare the rotation curves and dynamical mass profiles of the six ALPINE rotators to the two previously detected z ∼ 4–6 unlensed main-sequence rotators, finding high rotational velocities (∼50–250 km s−1) and a diversity of rotation curve shapes.
Abstract
We present results on the dust attenuation of galaxies at redshift ∼3–6 by studying the relationship between the UV spectral slope (βUV) and the infrared excess (IRX; L
IR/L
UV) using ...Atacama Large Millimeter/submillimeter Array (ALMA) far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of M
* ∼ 1010.7 M
⊙ spanning a redshift range z = 2.6–3.7 (median z = 3.2) that were observed with ALMA at
$\lambda _{\text{rest}}=300\,{\rm \mu m}$
. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX–βUV relationship at z ∼ 3 is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to ±0.5 dex. Nevertheless, the locally calibrated dust correction factors based on the IRX–βUV relation are on average applicable to main-sequence z ∼ 3 galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of z ∼ 4–6 galaxies we find evidence for a significant evolution in the IRX–βUV and the IRX–M
* relations beyond z ∼ 3 towards lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger z > 3 galaxy sample spanning a wide range of physical parameters (e.g. lower stellar mass) will be important to investigate this intriguing redshift evolution further.