Abstract
Star-forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy ...population with significant ionizing emission have proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star-forming galaxies in the 3 < z < 4 redshift range with deep rest-frame observations of the Lyman continuum (LyC). In our sample, we do not find any galaxy with significant emission of LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images, we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample, we measure meaningful 1σ upper limits of f
esc, rel < 0.07, 0.2 and 0.6 at
$L \sim L_{z=3}^{*}, 0.5L_{z=3}^{*}$
and
$0.1L_{z=3}^{*}$
, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of f
esc(M
UV). We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to M
UV ∼ −18 mag in the future.
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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. One of the main challenges in observational cosmology is related to the redshift evolution of the average hydrogen (HI) ionization in the Universe, as evidenced by the changing in ionization ...level of the intergalactic medium (IGM) through cosmic time. Starting from the first cosmic reionization, the rapid evolution of the IGM physical properties in particular poses severe constraints for the identification of the sources responsible for maintaining its high level of ionization up to lower redshifts. Aims. In order to probe the ionization level of the IGM and the ionization capabilities of bright quasi-stellar objects (QSOs) at z = 4, we selected a sample of 2508 QSOs drawn from the Sloan Digital Sky Survey (SDSS, DR14) in the redshift interval 3.6 ≤ z ≤ 4.6 and absolute magnitude range −29.0 ≲ M1450 ≲ −26.0. Particularly, we focus on the estimate of the escape fraction of HI-ionizing photons and their mean free path (MFP), which are fundamental for characterizing the surrounding IGM. Methods. Starting from UV/optical rest-frame spectra of the whole QSO sample from the SDSS survey, we estimated the escape fraction and free path individually for each of the QSOs. We calculated the Lyman continuum (LyC) escape fraction as the flux ratio blueward (∼900 Å rest frame) and redward (∼930 Å rest frame) of the Lyman limit. We then obtained the probability distribution function (PDF) of the individual free paths of the QSOs in the sample and studied its evolution in luminosity and redshift, comparing our results with those in literature. Results. We find a lower limit to the mean LyC escape fraction of 0.49, in agreement with the values obtained for both brighter and fainter sources at the same redshift. We show that the free paths of ionizing photons are characterized by a skewed distribution function that peaks at low values, with an average of ∼49 − 59 proper Mpc at z ∼ 4, after possible associated absorbers (AAs) were excluded. This value is higher than the one obtained at the same redshift by many authors in the literature using different techniques. Moreover, the PDF of free path gives information that is complementary to the MFP derived through the stacking technique. Finally, we also find that the redshift evolution of this parameter might be milder than previously thought. Conclusions. Our new determination of the MFP at z ∼ 4 implies that previous estimates of the HI photoionization rate ΓHI available in the literature should be corrected by a factor of 1.2−1.7. These results have important implications when they are extrapolated at the epoch of reionization.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. Finding the sources responsible for the hydrogen reionization is one of the most pressing issues in observational cosmology. Bright quasi-stellar objects (QSOs) are known to ionize their ...surrounding neighborhood, but they are too few to ensure the required HI ionizing background. A significant contribution by faint active galactic nuclei (AGNs), however, could solve the problem, as recently advocated on the basis of a relatively large space density of faint active nuclei at z > 4. Aims. This work is part of a long-term project aimed at measuring the Lyman Continuum escape fraction for a large sample of AGNs at z ~ 4 down to an absolute magnitude of M1450 ~ −23. We have carried out an exploratory spectroscopic program to measure the HI ionizing emission of 16 faint AGNs spanning a broad U − I color interval, with I ~ 21–23, and 3.6 < z < 4.2. These AGNs are three magnitudes fainter than the typical SDSS QSOs (M1450 ≲−26) which are known to ionize their surrounding IGM at z ≳ 4. Methods. We acquired deep spectra of these faint AGNs with spectrographs available at the VLT, LBT, and Magellan telescopes, that is, FORS2, MODS1-2, and LDSS3, respectively. The emission in the Lyman Continuum region, close to 900 Å rest frame, has been detected with a signal to noise ratio of ~10–120 for all 16 AGNs. The flux ratio between the 900 Å rest-frame region and 930 Å provides a robust estimate of the escape fraction of HI ionizing photons. Results. We have found that the Lyman Continuum escape fraction is between 44 and 100% for all the observed faint AGNs, with a mean value of 74% at 3.6 < z < 4.2 and − 25.1 ≲ M1450 ≲−23.3, in agreement with the value found in the literature for much brighter QSOs (M1450 ≲−26) at the same redshifts. The Lyman Continuum escape fraction of our faint AGNs does not show any dependence on the absolute luminosities or on the observed U − I colors of the objects. Assuming that the Lyman Continuum escape fraction remains close to ~75% down to M1450 ~ − 18, we find that the AGN population can provide between 16 and 73% (depending on the adopted luminosity function) of the whole ionizing UV background at z ~ 4, measured through the Lyman forest. This contribution increases to 25–100% if other determinations of the ionizing UV background are adopted from the recent literature. Conclusions. Extrapolating these results to z ~ 5–7, there are possible indications that bright QSOs and faint AGNs can provide a significant contribution to the reionization of the Universe, if their space density is high at M1450 ~ −23.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims. The goal of this work is to infer the star formation properties and the mass assembly process of high redshift ($0.3 \leq z < 2.5$) galaxies from their IR emission using the 24 μm band of ...MIPS-Spitzer. Methods. We used an updated version of the GOODS-MUSIC catalog, which has multiwavelength coverage from 0.3 to 24 μm and either spectroscopic or accurate photometric redshifts. We describe how the catalog has been extended by the addition of mid-IR fluxes derived from the MIPS 24 μm image. We compared two different estimators of the star formation rate (SFR hereafter). One is the total infrared emission derived from 24 μm, estimated using both synthetic and empirical IR templates. The other one is a multiwavelength fit to the full galaxy SED, which automatically accounts for dust reddening and age-star formation activity degeneracies. For both estimates, we computed the SFR density and the specific SFR. Results. We show that the two SFR indicators are roughly consistent, once the uncertainties involved are taken into account. However, they show a systematic trend, IR-based estimates exceeding the fit-based ones as the star formation rate increases. With this new catalog, we show that: a) at $z>0.3$, the star formation rate is correlated well with stellar mass, and this relationship seems to steepen with redshift if one relies on IR-based estimates of the SFR; b) the contribution to the global SFRD by massive galaxies increases with redshift up to $\simeq $2.5, more rapidly than for galaxies of lower mass, but appears to flatten at higher z; c) despite this increase, the most important contributors to the SFRD at any z are galaxies of about, or immediately lower than, the characteristic stellar mass; d) at $z\simeq 2$, massive galaxies are actively star-forming, with a median ${\it SFR} \simeq 300~ M_\odot$ yr-1. During this epoch, our targeted galaxies assemble a substantial part of their final stellar mass; e) the specific SFR (SSFR) shows a clear bimodal distribution. Conclusions. The analysis of the SFR density and the SSFR seems to support the downsizing scenario, according to which high mass galaxies have formed their stars earlier and more rapidly than their low mass counterparts. A comparison with renditions of theoretical simulations of galaxy formation and evolution indicates that these models follow the global increase in the SSFR with redshift and predict the existence of quiescent galaxies even at $z>1.5$. However, the average SSFR is systematically underpredicted by all models considered.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
ABSTRACT
We present results from the NIRVANDELS survey on the gas-phase metallicity (Zg, tracing O/H) and stellar metallicity (Z⋆, tracing Fe/H) of 33 star-forming galaxies at redshifts 2.95 < z < ...3.80. Based on a combined analysis of deep optical and near-IR spectra, tracing the rest-frame far-ultraviolet (FUV; 1200–2000 Å) and rest-frame optical (3400–5500 Å), respectively, we present the first simultaneous determination of the stellar and gas-phase mass–metallicity relationships (MZRs) at z ≃ 3.4. In both cases, we find that metallicity increases with increasing stellar mass (M⋆) and that the power-law slope at M⋆ ≲ 1010M⊙ of both MZRs scales as $Z \propto M_{\star }^{0.3}$. Comparing the stellar and gas-phase MZRs, we present direct evidence for super-solar O/Fe ratios (i.e. α-enhancement) at z > 3, finding (O/Fe) = 2.54 ± 0.38 × (O/Fe)⊙, with no clear dependence on M⋆.
ABSTRACT
We present Lyman continuum (LyC) radiation escape fraction (fesc) measurements for 183 spectroscopically confirmed star-forming galaxies in the redshift range 3.11 < z < 3.53 in the Chandra ...Deep Field South. We use ground-based imaging to measure fesc, and use ground- and space-based photometry to derive galaxy physical properties using spectral energy distribution (SED) fitting. We additionally derive O iii + H β equivalent widths (that fall in the observed K band) by including nebular emission in SED fitting. After removing foreground contaminants, we report the discovery of 11 new candidate LyC leakers at ≳ 2σ level, with fesc in the range 0.14−0.85. From non-detections, we place 1σ upper limits of fesc <0.12, where the Lyman-break selected galaxies have fesc <0.11 and ‘blindly’ discovered galaxies with no prior photometric selection have fesc <0.13. We find a slightly higher 1σ limit of fesc <0.20 from extreme emission line galaxies with rest-frame O iii + H β equivalent widths >300 Å. For candidate LyC leakers, we find a weak negative correlation between fesc and galaxy stellar masses, no correlation between fesc and specific star-formation rates (sSFRs) and a positive correlation between fesc and EW0(O iii + H β). The weak/no correlations between stellar mass and sSFRs may be explained by misaligned viewing angles and/or non-coincident time-scales of starburst activity and periods of high fesc. Alternatively, escaping radiation may predominantly occur in highly localized star-forming regions, or fesc measurements may be impacted by stochasticity of the intervening neutral medium, obscuring any global trends with galaxy properties. These hypotheses have important consequences for models of reionization.
We present a direct comparison between the observed star formation rate functions (SFRFs) and the state-of-the-art predictions of semi-analytic models (SAMs) of galaxy formation and evolution. We use ...the PACS Evolutionary Probe Survey and Herschel Multi-tiered Extragalactic Survey data sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z ∼ 4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z ∼ 2, when the observational errors on the SFR are taken into account. However, all the models seem to underpredict the bright end of the SFRF at z ≳ 2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fallback of gas, caused by weak feedback and outflows at earlier epochs.
Aims. We wish to investigate the physical properties of a sample of Lyα emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in ...the escape and spatial distribution of Lyα photons. Methods. From all the Lyα emitting galaxies in the VANDELS Data Release 2 at 3.5 ≲ z ≲ 4.5, we selected a sample of 52 galaxies that also have a precise systemic redshift determination from at least one nebular emission line (HeII or CIII). For these galaxies, we derived different physical properties (stellar mass, age, dust extinction, and star formation rate) from spectral energy distribution (SED) fitting of the exquisite multiwavelength photometry available in the VANDELS fields, using the dedicated spectral modeling tool BEAGLE and the UV β slope from the observed photometry. We characterized the Lyα emission in terms of kinematics, equivalent width (EW), full width at half-maximum, and spatial extension and then estimated the velocity of the neutral outflowing gas. The ultra-deep VANDELS spectra (up to 80 h on-source integration) enable this for individual galaxies without the need to rely on stacks. We then investigated the correlations between the Lyα properties and the other measured properties to study how they affect the shape and intensity of Lyα emission. Results. We reproduce some of the well-known correlations between Lyα EW and stellar mass, dust extinction, and UV β slope, in the sense that the emission line appears brighter in galaxies with lower mass that are less dusty and bluer. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Lyα tend to be more compact in both UV and in Lyα. Our data reveal an interesting correlation between the Lyα velocity offset and the shift of the interstellar absorption lines with respect to the systemic redshift, observed for the first time at high redshifts: galaxies with higher interstellar medium (ISM) outflow velocities show smaller Lyα velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Lyα kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Lyα spatial profiles; this is a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrowband imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f esc) of ionizing ...radiation for L >= L* z=3 Lyman-break galaxies (LBGs) at redshift 3.4-4.5. Only one LBG, at redshift z = 3.795, is detected in its Lyman continuum (LyC; S/N 5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R eff = 0.8 kpc physical). Three out of seven active galactic nuclei are also detected in their LyC, including one at redshift z = 3.951 and z 850 = 26.1. From stacked data (LBGs), we set an upper limit to the average f esc in the range 5%-20%, depending on how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction, and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal, and Gaussian) are explored. We find that the median f esc is lower than 6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f esc, and estimates of f esc from a sample of faint LBGs from HUDF (i 775 <=28.5) are in broad quantitative agreement with such a scenario.