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.
Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio active galactic nuclei (AGNs) to be used to ...their fullest potential. We study the composition of the Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of ∼350 arcmin2. This ultra-deep survey consists of a single pointing in the well-studied COSMOS field at both 3 and 10 GHz and reaches rms sensitivities of 0.53 and 0.41 Jy beam−1, respectively. We find multiwavelength counterparts for 97% of radio sources, based on a combination of near-UV/optical to sub-millimeter data, and through a stacking analysis at optical/near-IR wavelengths we further show that the sources lacking such counterparts are likely to be high-redshift in nature (typical z ∼ 4−5). Utilizing the multiwavelength data over COSMOS, we identify AGNs through a variety of diagnostics and find these to make up 23.2 1.3% of our sample, with the remainder constituting uncontaminated star-forming galaxies. However, more than half of the AGNs exhibit radio emission consistent with originating from star formation, with only 8.8 0.8% of radio sources showing a clear excess in radio luminosity. At flux densities of ∼30 Jy at 3 GHz, the fraction of star formation-powered sources reaches ∼90%, and this fraction is consistent with unity at even lower flux densities. Overall, our findings imply that ultra-deep radio surveys such as COSMOS-XS constitute a highly effective means of obtaining clean samples of star formation-powered radio sources.
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.
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
We present ultradeep, matched-resolution Karl G. Jansky Very Large Array observations at 10 and 3 GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and 3 GHz images cover ∼16 and
...and reach median rms values at the phase center of 0.41 and 0.53
μ
Jy beam
−1
, respectively. Both images have an angular resolution of ∼20. To account for the spectral shape and resolution variations across the broad bands, we image all data with a multiscale, multifrequency synthesis algorithm. We present source catalogs for the 10 and 3 GHz image with 91 and 1498 sources, respectively, above a peak brightness threshold of 5
σ
. We present source counts with completeness corrections included that are computed via Monte Carlo simulations. Our corrected counts at 3 GHz are consistent within the uncertainties with other results at 3 and 1.4 GHz but extend to fainter flux densities than previous direct detections. The 3 GHz number counts exceed the counts predicted by the semiempirical simulations developed in the framework of the SKA Simulated Skies project, consistent with previous
P
(
D
) analyses. Our source counts suggest a steeper luminosity function evolution for faint star-forming sources. The semiempirical Tiered Radio Extragalactic Continuum Simulation predicts this steeper evolution and is in better agreement with our results at 10 and 3 GHz within the expected variations from cosmic variance. In summary, the multiband, matched-resolution COSMOS-XS survey in the COSMOS field provides a high-resolution view of the ultrafaint radio sky that can help guide next-generation radio facilities.
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
Attempts to trace star formation with rest-frame UV/optical observations at redshifts
z
> 2 are affected by the presence of potentially substantial, yet uncertain, dust attenuation. Recent ...studies have demonstrated the existence of a population of galaxies that are virtually invisible in the observed optical/near-infrared (NIR) due to dust obscuration, but which could contribute substantially to the star formation history at 2 <
z
< 6. Here, we make use of ultradeep 3 GHz Karl G. Jansky Very Large Array observations from the COSMOS-XS survey to investigate the contribution 20of radio-selected “optically dark” galaxies (undetected to a depth of
K
S
∼ 25.9 mag) to the cosmic star formation rate density (SFRD). We identify 19 such “optically dark” sources and utilize recent deblended far-infrared photometry to determine photometric redshifts based on IR and radio information for 11 of them. Through stacking, we infer that the remaining eight sources reside predominantly at high redshift (
z
> 4). Therefore, we conservatively assume these sources lie between
z
= 2 and
z
= 5. We derive the radio luminosity function (LF) for the sample with and without “optically dark” sources by fixing the faint and bright end shape of the radio LF to the local values and allowing for luminosity evolution. By integrating both LFs, we estimate the contribution of the “optically dark” galaxies to the radio SFRD to be
∼
15
−
7
+
7
%
at
z
∼ 5. This is consistent with constraints from NIR-dark and UV-dark sources, while being in disagreement with some estimates using
H
-dropouts. This result implies that “optically dark” sources play a nonnegligible role at high redshift.
Abstract
We present Atacama Large Millimeter/submillimeter Array (ALMA) 2 mm continuum observations of a complete and unbiased sample of 99 870
μ
m selected submillimeter galaxies (SMGs) in the ...Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53
μ
Jy beam
−1
. We measure the flux densities for 70 sources, for which we obtain a typical 870
μ
m-to-2 mm flux ratio of 14 ± 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2 mm observations, the dust emissivity index is well constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is
β
≃ 1.9 ± 0.4, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust-grain model predictions. We also find a negative correlation between the dust temperature and
β
, similar to low-redshift observational and theoretical studies. Our results indicate that
β
≃ 2 in high-redshift dusty star-forming galaxies, implying little evolution in dust-grain properties between our SMGs and local dusty galaxy samples, and suggesting that these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their interstellar medium.
Abstract
The high-frequency radio sky has historically remained largely unexplored due to the typical faintness of sources in this regime, and the modest survey speed compared to observations at ...lower frequencies. However, high-frequency radio surveys offer an invaluable tracer of high-redshift star formation, as they directly target the faint radio free–free emission. We present deep continuum observations at 34 GHz in the COSMOS and GOODS-North fields from the Karl G. Jansky Very Large Array (VLA), as part of the COLD
z
survey. The deep COSMOS mosaic spans
down to
σ
= 1.3
μ
Jy beam
−1
, while the wider GOODS-N observations cover
to
σ
= 5.3
μ
Jy beam
−1
. We detect a total of 18 galaxies at 34 GHz, of which nine show radio emission consistent with being powered by star formation; although for two sources, this is likely due to thermal emission from dust. Utilizing deep ancillary radio data at 1.4, 3, 5, and 10 GHz, we decompose the spectra of the remaining seven star-forming galaxies into their synchrotron and thermal free–free components, and find typical thermal fractions and synchrotron spectral indices comparable to those observed in local star-forming galaxies. We further determine free–free star formation rates (SFRs), and show that these are in agreement with SFRs from spectral energy distribution-fitting and the far-infrared/radio correlation. Our observations place strong constraints on the high-frequency radio emission in typical galaxies at high redshift, and provide some of the first insights into what is set to become a key area of study with future radio facilities, such as the Square Kilometer Array Phase 1 and next-generation VLA.
Abstract
We present the average C
ii
158
μ
m emission line sizes of UV-bright star-forming galaxies at
z
∼ 7. Our results are derived from a stacking analysis of C
ii
158
μ
m emission lines and ...dust continua observed by the Atacama Large Millimeter/submillimeter Array (ALMA), taking advantage of the large program Reionization Era Bright Emission Line Survey. We find that the average C
ii
emission at
z
∼ 7 has an effective radius
r
e
of 2.2 ± 0.2 kpc. It is ≳2× larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for
z
≲ 6 galaxies. Additionally, we compared the average C
ii
size with 4 <
z
< 6 galaxies observed by the ALMA Large Program to INvestigate C
ii
at Early times (ALPINE). By analyzing C
ii
sizes of 4 <
z
< 6 galaxies in two redshift bins, we find an average C
ii
size of
r
e
= 2.2 ± 0.2 kpc and
r
e
= 2.5 ± 0.2 kpc for
z
∼ 5.5 and
z
∼ 4.5 galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the C
ii
158
μ
m emitting regions at redshift between
z
∼ 7 and
z
∼ 4. This finding suggests that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.