We report Atacama Large Millimeter Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-infrared (FIR)/submillimeter observations in the COSMOS ...field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z = 3.62-5.85 using ALMA 3 mm line scans, detecting multiple CO and/or C i transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z = 5.847. These redshifts are lower than we had expected, as these galaxies have substantially colder dust temperatures (i.e., their spectral energy distributions peak at longer rest-frame wavelengths) than most literature sources at z > 4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, cosmic microwave background (CMB) plays a significant role in biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High-resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were underrepresented in current surveys and are even colder than typical main-sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift.
Abstract
We investigate the fraction of close pairs and morphologically identified mergers on and above the star-forming main sequence (MS) at 0.2 ≤ z ≤2.0. The novelty of our work lies in the use of ...a non-parametric morphological classification performed on resolved stellar mass maps, reducing the contamination by non-interacting, high-redshift clumpy galaxies. We find that the merger fraction rapidly rises to ≥70 per cent above the MS, implying that – already at z ≳ 1 – starburst (SB) events (ΔMS ≥ 0.6) are almost always associated with a major merger (1:1 to 1:6 mass ratio). The majority of interacting galaxies in the SB region are morphologically disturbed, late-stage mergers. Pair fractions show little dependence on MS offset and pairs are more prevalent than late-stage mergers only in the lower half of the MS. In our sample, major mergers on the MS occur with a roughly equal frequency of ∼5–10 per cent at all masses ≳ 1010 M⊙. The MS major merger fraction roughly doubles between z = 0.2 and 2, with morphological mergers driving the overall increase at z ≳ 1. The differential redshift evolution of interacting pairs and morphologically classified mergers on the MS can be reconciled by evolving observability time-scales for both pairs and morphological disturbances. The observed variation of the late-stage merger fraction with ΔMS follows the perturbative 2-Star Formation Mode model, where any MS galaxy can experience a continuum of different star formation rate enhancements. This points to an SB–merger connection not only for extreme events, but also more moderate bursts which merely scatter galaxies upward within the MS, rather than fully elevating them above it.
The emergence of passive galaxies in the early Universe results from the delicate interplay among the different physical processes responsible for their rapid assembly and the abrupt shut-down of ...their star formation activity. Investigating the individual properties and demographics of early passive galaxies improves our understanding of these mechanisms. In this work we present a follow-up analysis of the
z
> 3 passive galaxy candidates selected by Merlin et al. (2019, MNRAS, 490, 3309) in the CANDELS fields. We begin by first confirming the accuracy of their passive classification by exploiting their sub-millimetre emission to demonstrate the lack of ongoing star formation. Using archival ALMA observations we are able to confirm at least 61% of the observed candidates as passive. While the remainder lack sufficiently deep data for confirmation, we are able to validate the entire sample in a statistical sense. We then estimate the stellar mass function (SMF) of all 101 passive candidates in three redshift bins from
z
= 5 to
z
= 3. We adopt a stepwise approach that has the advantage of taking into account photometric errors, mass and selection completeness issues, as well as the Eddington bias, without any a posteriori correction. We observe a pronounced evolution in the SMF around
z
∼ 4, indicating that we are witnessing the emergence of the passive population at this epoch. Massive (
M
> 10
11
M
⊙
) passive galaxies, only accounting for a small (< 10%) fraction of galaxies at
z
> 4, become dominant at later epochs. Thanks to a combination of photometric quality, sample selection, and methodology, we overall find a higher density of passive galaxies than in previous works. The comparison with theoretical predictions, despite a qualitative agreement (at least for some of the models considered), denotes a still incomplete understanding of the physical processes responsible for the formation of these galaxies. Finally, we extrapolate our results to predict the number of early passive galaxies expected in surveys carried out with future facilities.
Abstract
We present the first rest-frame optical size–luminosity relation of galaxies at
z
> 7, using the NIRCam imaging data obtained by the GLASS James Webb Space Telescope Early Release Science ...(GLASS-JWST-ERS) program, providing the deepest extragalactic data of the ERS campaign. Our sample consists of 19 photometrically selected bright galaxies with
m
F444W
≤ 27.8 at 7 <
z
< 9 and
m
F444W
< 28.2 at
z
∼ 9−15. We measure the size of the galaxies in five bands, from rest-frame optical (∼4800 Å) to the UV (∼1600 Å) based on the Sérsic model, and analyse the size–luminosity relation as a function of wavelength. Remarkably, the data quality of the NIRCam imaging is sufficient to probe the half-light radius
r
e
down to ∼100 pc at
z
> 7. Given the limited sample size and magnitude range, we first fix the slope to that observed for larger samples in rest-frame UV using Hubble Space Telescope samples. The median size
r
0
at the reference luminosity
M
= −21 decreases slightly from rest-frame optical (600 ± 80 pc) to UV (450 ± 130 pc). We then refit the size–luminosity relation allowing the slope to vary. The slope is consistent with
β
∼ 0.2 for all bands except F150W, where we find a marginally steeper slope of
β
= 0.53 ± 0.15. The steep UV slope is mainly driven by the smallest and faintest galaxies. If confirmed by larger samples, it implies that the UV size–luminosity relation breaks toward the faint end, as suggested by lensing studies.
Abstract
We present the first James Webb Space Telescope/NIRCam-led determination of 7 <
z
< 9 galaxy properties based on broadband imaging from 0.8 to 5
μ
m as part of the GLASS-JWST Early Release ...Science program. This is the deepest data set acquired at these wavelengths to date, with an angular resolution ≲0.″14. We robustly identify 13 galaxies with signal-to-noise ratio ≳ 8 in F444W from 8 arcmin
2
of data at
m
AB
≤ 28 from a combination of dropout and photometric redshift selection. From simulated data modeling, we estimate the dropout sample purity to be ≳90%. We find that the number density of these F444W-selected sources is broadly consistent with expectations from the UV luminosity function determined from Hubble Space Telescope data. We characterize galaxy physical properties using a Bayesian spectral energy distribution fitting method, finding a median stellar mass of 10
8.5
M
⊙
and age 140 Myr, indicating they started ionizing their surroundings at redshift
z
> 9.5. Their star formation main sequence is consistent with predictions from simulations. Lastly, we introduce an analytical framework to constrain main-sequence evolution at
z
> 7 based on galaxy ages and basic assumptions, through which we find results consistent with expectations from cosmological simulations. While this work only gives a glimpse of the properties of typical galaxies that are thought to drive the reionization of the universe, it clearly shows the potential of JWST to unveil unprecedented details of galaxy formation in the first billion years.
ABSTRACT
We present a study designed to measure the average Lyman-continuum escape fraction (〈fesc〉) of star-forming galaxies at z ≃ 3.5. We assemble a sample of 148 galaxies from the VANDELS ...spectroscopic survey at 3.35 ≤ zspec ≤ 3.95, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, U-band imaging and Hubble Space Telescope V-band imaging to robustly measure the distribution of $\mathcal {R_{\rm obs}}\, =(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the $\mathcal {R_{\rm obs}}$ distribution as a function of 〈fesc〉, carefully accounting for attenuation by dust, the intergalactic medium and the circumgalactic medium. A maximum likelihood fit to the $\mathcal {R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07^{+0.02}_{-0.02}$, a result confirmed using an alternative Bayesian inference technique (both techniques exclude 〈fesc〉 = 0.0 at >3σ). By splitting our sample in two, we find evidence that 〈fesc〉 is positively correlated with Ly α equivalent width (Wλ(Ly α)), with high and low Wλ(Lyα) subsamples returning values of $\langle f_{\rm esc}\rangle =0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle =0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that 〈fesc〉 is anticorrelated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation subsamples both returning best fits in the range 0.10 ≤ 〈fesc〉 ≤ 0.22. We do not find a clear correlation between fesc and galaxy stellar mass, suggesting stellar mass is not a primary indicator of fesc. Although larger samples are needed to further explore these trends, our results suggest that it is entirely plausible that the low dust, low-metallicity galaxies found at z ≥ 6 will display the 〈fesc〉 ≥ 0.1 required to drive reionization.
Abstract
Star-forming galaxies can exhibit strong morphological differences between the rest-frame far-UV and optical, reflecting inhomogeneities in star formation and dust attenuation. We exploit ...deep, high-resolution, NIRCAM seven-band observations to take a first look at the morphology of galaxies in the epoch of reionization (
z
> 7), and its variation in the rest-frame wavelength range between Ly
α
and 6000–4000 Å, at
z
= 7–12. We find no dramatic variations in morphology with wavelength—of the kind that would have overturned anything we have learned from the Hubble Space Telescope. No significant trends between morphology and wavelengths are detected using standard quantitative morphology statistics. We detect signatures of mergers/interactions in 4/19 galaxies. Our results are consistent with a scenario in which Lyman-break galaxies—observed when the universe is only 400–800 Myr old—are growing via a combination of rapid, galaxy-scale star formation supplemented by the accretion of star-forming clumps and interactions.
Abstract
We exploit James Webb Space Telescope (JWST) NIRCam observations from the GLASS-JWST-Early Release Science program to investigate galaxy stellar masses at
z
> 7. We first show that JWST ...observations reduce the uncertainties on the stellar mass by a factor of at least 5–10, when compared with the highest-quality data sets available to date. We then study the UV mass-to-light ratio, finding that galaxies exhibit a a two orders of magnitude range of
M
/
L
UV
values for a given luminosity, indicative of a broad variety of physical conditions and star formation histories. As a consequence, previous estimates of the cosmic stellar-mass density—based on an average correlation between UV luminosity and stellar mass—can be biased by as much as a factor of ∼6. Our first exploration demonstrates that JWST represents a new era in our understanding of stellar masses at
z
> 7 and, therefore, of the growth of galaxies prior to cosmic reionization.
The estimate of stellar metallicities (
Z
*
) of high-
z
galaxies are of paramount importance in order to understand the complexity of dust effects and the reciprocal interrelations among stellar ...mass, dust attenuation, stellar age, and metallicity. Benefiting from uniquely deep far-UV spectra of > 500 star-forming galaxies at redshifts 2 <
z
< 5 extracted from the VANDELS survey and stacked in bins of stellar mass (
M
*
) and UV continuum slope (
β
), we estimate their stellar metallicities
Z
*
from stellar photospheric absorption features at 1501 and 1719 Å, which are calibrated with Starburst99 models and are largely unaffected by stellar age, dust, IMF, nebular continuum, or interstellar absorption. Comparing them to photometric-based spectral slopes in the 1250–1750 Å range, we find that the stellar metallicity increases by ∼0.5 dex from
β
∼ −2 to
β
∼ −1 (1 ≲
A
1600
≲ 3.2), and a dependence with
β
holds at fixed UV absolute luminosity
M
UV
and stellar mass up to ∼10
9.65
M
⊙
. As a result, metallicity is a fundamental ingredient for properly rescaling dust corrections based on
M
UV
and
M
*
. Using the same absorption features, we analyzed the mass-metallicity relation (MZR), and find it to be consistent with the previous VANDELS estimation based on a global fit of the FUV spectra. Similarly, we do not find a significant evolution between
z
∼ 2 and
z
∼ 3.5. Finally, the slopes of our MZR and
Z
*
−
β
relation are in agreement with the predictions of well-studied semi-analytic models (SAM) of galaxy formation, while some tensions remain concerning the absolute metallicity normalization. The relation between the UV slope and stellar metallicity is fundamental to the exploitation of large volume surveys with next-generation telescopes and for the physical characterization of galaxies in the first billion years of our Universe.
ABSTRACT
The physical properties of Epoch of Reionization (EoR) galaxies are still poorly constrained by observations. To better understand the ionizing properties of galaxies in the EoR, we ...investigate deep, rest-frame ultraviolet (UV) spectra of ≃500 star-forming galaxies at 3 ≤ z ≤ 5 selected from the public ESO-VANDELS spectroscopic survey. The absolute ionizing photon escape fraction ($f_{\rm esc}^{\rm abs}$, i.e. the ratio of leaking against produced ionizing photons) is derived by combining absorption line measurements with estimates of the UV attenuation. The ionizing production efficiency (ξion, i.e. the number of ionizing photons produced per non-ionizing UV luminosity) is calculated by fitting the far-UV (FUV) stellar continuum of the VANDELS galaxies. We find that the $f_{\rm esc}^{\rm abs}$ and ξion parameters increase towards low-mass, blue UV-continuum slopes and strong Ly α emitting galaxies, and both are slightly higher-than-average for the UV-faintest galaxies in the sample. Potential Lyman Continuum Emitters (LCEs, $f_{\rm esc}^{\rm abs} \ge 5{{\ \rm \, per\ cent}}$) and selected Lyman Alpha Emitters (LAEs, WLyα ≤ −20 Å) show systematically higher ξion (log ξion(Hz erg−1) ≈ 25.38, 25.41) than non-LCEs and non-LAEs (log ξion(Hz erg−1) ≈ 25.18, 25.14) at similar UV magnitudes. This indicates very young underlying stellar populations (≈10 Myr) at relatively low metallicities (≈0.2 Z⊙). The FUV non-ionizing spectra of potential LCEs is characterized by blue UV slopes (≤−2), enhanced Ly α emission (≤−25 Å), strong UV nebular lines (e.g. high ${\rm C\, \small {IV}}$1550/${\rm C\, \small {III}}$1908 ≥0.75 ratios), and weak absorption lines (≤1 Å). The latter suggests the existence of low gas-column-density channels in the interstellar medium, which enables the escape of ionizing photons. By comparing our VANDELS results against other surveys in the literature, our findings imply that the ionizing budget in the EoR was likely dominated by UV-faint, low-mass, and dustless galaxies.