Title: The Effects of Stellar Population and Gas Covering Fraction on the Emergent Lyα Emission of High-redshift Galaxies*
Abstract We perform joint modeling of the composite rest-frame far-UV and optical spectra of redshift 1.85 ≤ z ≤ 3.49 star-forming galaxies to deduce key properties of the massive stars, ionized interstellar medium (ISM), and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Ly α photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Ly β through C iii ] λλ 1907, 1909 and [O ii ], [Ne iii ], H β , [O iii ], H α , [N ii ], and [S ii ]. Spectral and photoionization modeling indicates that the galaxies are uniformly consistent with stellar population synthesis models that include the effects of stellar binarity. Over the dynamic range of our sample, there is little variation in stellar and nebular abundance with Ly α equivalent width, W λ (Ly α ), and only a marginal anticorrelation between age and W λ (Ly α ). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W λ (Ly α ); rather, the covering fraction of optically thick H i appears to be the principal factor modulating the escape of Ly α , with most of the Ly α photons in down-the-barrel observations of galaxies escaping through low column density or ionized channels in the ISM. Our analysis shows that a high star-formation-rate surface density, Σ SFR , particularly when coupled with a low galaxy potential (i.e., low stellar mass), can aid in reducing the covering fraction and ease the escape of Ly α photons. We conclude with a discussion of the implications of our results for the escape of ionizing radiation at high redshift.  more » « less
Award ID(s):
2009313 2009278 2009085
PAR ID:
10336789
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
926
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
31
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We present an updated model of the cosmic ionizing background from the UV to the X-rays. Relative to our previous model, the new model provides a better match to a large number of up-to-date empirical constraints, including: (1) new galaxy and AGN luminosity functions; (2) stellar spectra including binary stars; (3) obscured and unobscured AGN; (4) a measurement of the non-ionizing UV background; (5) measurements of the intergalactic H i and He ii photoionization rates at z ∼ 0−6; (6) the local X-ray background; and (7) improved measurements of the intergalactic opacity. In this model, AGN dominate the H i ionizing background at z ≲ 3 and star-forming galaxies dominate it at higher redshifts. Combined with the steeply declining AGN luminosity function beyond z ∼ 2, the slow evolution of the H i ionization rate inferred from the high-redshift H i Ly α forest requires an escape fraction from star-forming galaxies that increases with redshift (a population-averaged escape fraction of $$\approx 1{{\ \rm per\ cent}}$$ suffices to ionize the intergalactic medium at z = 3 when including the contribution from AGN). We provide effective photoionization and photoheating rates calibrated to match the Planck 2018 reionization optical depth and recent constraints from the He ii Ly α forest in hydrodynamic simulations. 
    more » « less
  2. Abstract We analyze several key rest-optical emission-line ratios in a sample of 828 galaxies as well as composite spectra from JADES DR3 in the range 1.4 < z < 7. These emission-line ratios include: [Oiii]λ5008/Hβ, [Nii]λ6585/Hα, [Sii]λλ6718, 6733/Hα, [Oi]λ6302/Hα, O32, R23, Ne3O2, and RO2Ne3. We find evidence for a harder ionizing spectrum atz ∼ 3.5 compared toz ∼ 2 at fixed gas-phase metallicity, resulting in a pronounced shift in the star-forming galaxy locus on the [Nii]/HαBaldwin-Phillips-Terlevich (BPT) diagram and the O32 versus R23 diagram. Atz ≳ 3.5, star-forming galaxies occupy a common locus, indicating that interstellar medium (ISM) ionizing conditions at fixed gas-phase metallicity do not evolve strongly at these early cosmic times. There is a connection between ISM ionizing conditions and the chemical abundance patterns (i.e.,α/Fe) in massive stars, providing the ionizing radiation field. Therefore, the lack of evolution in ISM ionizing conditions atz ≳ 3.5, followed by evolution towards a softer ionizing spectrum at fixed nebular metallicity as cosmic time proceeds toz ∼ 2 and lower redshift mirrors the chemical abundance patterns in Milky Way stars as a function of iron abundance. Our results highlight the diagnostic power of emission-line diagrams in the era of JWST to further our understanding of the ISM conditions into the Epoch of Reionization. 
    more » « less
  3. Context. The Lyman- α line in the ultraviolet (UV) and the [CII] line in the far-infrared (FIR) are widely used tools to identify galaxies in the early Universe and to obtain insights into interstellar medium (ISM) properties in high-redshift galaxies. By combining data obtained with ALMA in band 7 at ∼320 GHz as part of the ALMA Large Program to INvestigate [CII] at Early Times (ALPINE) with spectroscopic data from DEIMOS at the Keck Observatory, VIMOS and FORS2 at the Very Large Telescope, we assembled a unique sample of 53 main-sequence star-forming galaxies at 4.4 <   z  <  6 in which we detect both the Lyman- α line in the UV and the [CII] line in the FIR. Aims. The goal of this paper is to constrain the properties of the Ly α emission in these galaxies in relation to other properties of the ISM. Methods. We used [CII], observed with ALMA, as a tracer of the systemic velocity of the galaxies, and we exploited the available optical spectroscopy to obtain the Ly α -[CII] and ISM-[CII] velocity offsets. Results. We find that 90% of the selected objects have Ly α -[CII] velocity offsets in the range 0 <  Δ v Ly α  − [CII]  <  400 km s −1 , in line with the few measurements available so far in the early Universe, and significantly smaller than those observed at lower redshifts. At the same time, we observe ISM-[CII] offsets in the range −500 <  Δ v ISM−[CII]  <  0 km s −1 , in line with values at all redshifts, which we interpret as evidence for outflows in these galaxies. We find significant anticorrelations between Δ v Ly α −[CII] and the Ly α rest-frame equivalent width EW 0 (Ly α ) (or equivalently, the Ly α escape fraction f esc (Ly α )): galaxies that show smaller Δ v Ly α −[CII] have larger EW 0 (Ly α ) and f esc (Ly α ). Conclusions. We interpret these results in the framework of available models for the radiative transfer of Ly α photons. According to the models, the escape of Ly α photons would be favored in galaxies with high outflow velocities, producing large EW 0 (Ly α ) and small Δ v Ly α -[CII] , in agreement with our observations. The uniform shell model would also predict that the Ly α escape in galaxies with slow outflows (0 <   v out  <  300 km s −1 ) is mainly determined by the neutral hydrogen column density (NHI) along the line of sight, while the alternative model by Steidel et al. (2010, ApJ, 717, 289) would more highly favor a combination of NHI at the systemic velocity and covering fraction as driver of the Ly α escape. We suggest that the increase in Ly α escape that is observed in the literature between z  ∼ 2 and z  ∼ 6 is not due to a higher incidence of fast outflows at high redshift, but rather to a decrease in average NHI along the line of sight, or alternatively, a decrease in HI covering fraction. 
    more » « less
  4. Abstract We present spatially resolved Hubble Space Telescope grism spectroscopy of 15 galaxies at z ∼ 0.8 drawn from the DEEP2 survey. We analyze H α +[N ii ], [S ii ], and [S iii ] emission on kiloparsec scales to explore which mechanisms are powering emission lines at high redshifts, testing which processes may be responsible for the well-known offset of high-redshift galaxies from the z ∼ 0 locus in the [O iii ]/H β versus [N ii ]/H α Baldwin—Phillips—Terlevich (BPT) excitation diagram. We study spatially resolved emission-line maps to examine evidence for active galactic nuclei (AGN), shocks, diffuse ionized gas (DIG), or escaping ionizing radiation, all of which may contribute to the BPT offsets observed in our sample. We do not find significant evidence of AGN in our sample and quantify that, on average, AGN would need to contribute ∼25% of the H α flux in the central resolution element in order to cause the observed BPT offsets. We find weak (2 σ ) evidence of DIG emission at low surface brightnesses, yielding an implied total DIG emission fraction of ∼20%, which is not significant enough to be the dominant emission line driver in our sample. In general we find that the observed emission is dominated by star-forming H ii regions. We discuss trends with demographic properties and the possible role of α -enhanced abundance patterns in the emission spectra of high-redshift galaxies. Our results indicate that photoionization modeling with stellar population synthesis inputs is a valid tool to explore the specific star formation properties which may cause BPT offsets, to be explored in future work. 
    more » « less
  5. Abstract We use Hubble Space Telescope Wide Field Camera 3 G102 and G141 grism spectroscopy to measure rest-frame optical emission-line ratios of 533 galaxies at z ∼ 1.5 in the CANDELS Ly α Emission at Reionization survey. We compare [O iii ]/H β versus [S ii ]/(H α + [N ii ]) as an “unVO87” diagram for 461 galaxies and [O iii ]/H β versus [Ne iii ]/[O ii ] as an “OHNO” diagram for 91 galaxies. The unVO87 diagram does not effectively separate active galactic nuclei (AGN) and [Ne v ] sources from star-forming galaxies, indicating that the unVO87 properties of star-forming galaxies evolve with redshift and overlap with AGN emission-line signatures at z > 1. The OHNO diagram does effectively separate X-ray AGN and [Ne v ]-emitting galaxies from the rest of the population. We find that the [O iii ]/H β line ratios are significantly anticorrelated with stellar mass and significantly correlated with log ( L H β ) , while [S ii ]/(H α + [N ii ]) is significantly anticorrelated with log ( L H β ) . Comparison with MAPPINGS V photoionization models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-star formation rate (SFR) galaxies. We do not find evidence for redshift evolution of the emission-line ratios outside of the correlations with mass and SFR. Our results suggest that the OHNO diagram of [O iii ]/H β versus [Ne iii ]/[O ii ] will be a useful indicator of AGN content and gas conditions in very high-redshift galaxies to be observed by the James Webb Space Telescope. 
    more » « less