skip to main content

Title: The physics of Lyman-α escape from disc-like galaxies
ABSTRACT

Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$), helium ionization ($f_\text{He} \approx 8.7\,\rm{per\,\,cent}$), and anisotropic escape fractions ($f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$), as these reduce the available budget for hydrogen line emission ($f_\text{H} \approx 55.9\,\rm{per\,\,cent}$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but more » would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.

« less
Authors:
; ; ; ; ; ; ; ; ; ;
Publication Date:
NSF-PAR ID:
10373259
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
517
Issue:
1
Page Range or eLocation-ID:
p. 1-27
ISSN:
0035-8711
Publisher:
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rmmore »per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

    « less
  2. ABSTRACT

    The ionizing photon escape fraction [Lyman continuum (LyC) fesc] of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshifts. Here we circumvent this uncertainty by inferring fesc using resolved (R > 4000) Lyman α (Lyα) profiles from the X-SHOOTER Lyα survey at z = 2 (XLS-z2). With empirically motivated criteria, we use Lyα profiles to select leakers ($f_{\mathrm{ esc}} > 20{{\ \rm per\ cent}}$) and non-leakers ($f_{\mathrm{ esc}} < 5{{\ \rm per\ cent}}$) from a representative sample of >0.2L* Lyman α emitters (LAEs). We use median stacked spectra of these subsets over λrest ≈ 1000–8000 Å to investigate the conditions for LyC fesc. Our stacks show similar mass, metallicity, MUV, and βUV. We find the following differences between leakers versus non-leakers: (i) strong nebular C iv and He ii emission versus non-detections; (ii) [O iii]/[O ii] ≈ 8.5 versus ≈3; (iii) Hα/Hβ indicating no dust versus E(B − V) ≈ 0.3; (iv) Mg ii emission close to the systemic velocity versus redshifted, optically thick Mg ii; and (v) Lyα fesc of ${\approx} 50{{\ \rm per\ cent}}$ versus ${\approx} 10{{\ \rm per\ cent}}$. The extreme equivalent widths (EWs) in leakers ([O iii]+$\mathrm{ H}\beta \approx 1100$ Å rest frame)more »constrain the characteristic time-scale of LyC escape to ≈3–10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers – extremely ionizing stellar populations, low column densities, a dust-free, high-ionization state interstellar medium (ISM) – occur simultaneously in the $f_{\rm esc} > 20{{\ \rm per\ cent}}$ stack, suggesting they are causally connected, and motivating why indicators like [O iii]/[O ii] may suffice to constrain fesc at z > 6 with the James Webb Space Telescope (JWST). The leakers comprise half of our sample, have a median LyC$f_{\rm esc} \approx 50{{\ \rm per\ cent}}$ (conservative range: $20\!-\!55{{\ \rm per\ cent}}$), and an ionizing production efficiency $\log ({\xi _{\rm {ion}}/\rm {Hz\ erg^{-1}}})\approx 25.9$ (conservative range: 25.7–25.9). These results show LAEs – the type of galaxies rare at z ≈ 2, but that become the norm at higher redshift – are highly efficient ionizers, with extreme ξion and prolific fesc occurring in sync.

    « less
  3. ABSTRACT

    We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lyα emitters (LAEs) at 2.36 < z < 2.82 selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as LLyα ≃ 30 × 1042 erg s−1 and MUV ≃ −23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lyα rest-frame equivalent widths, $EW_{0}\,\rm (Ly\alpha)$, ranging from 16 to 50 Å and Lyα escape fractions, $f_{\rm esc}\, \rm (Ly\alpha)$, from 10 per cent to 40 per cent. Large $EW_{0}\, \rm (Ly\alpha)$ and $f_{\rm esc}\, \rm (Ly\alpha)$ are found predominantly in LAEs showing weak low-ionization ISM absorption (EW0 ≲ 1 Å) and narrow Lyα profiles (≲300 km s−1 FWHM) with their peak close (≲80 km s−1) to their systemic redshifts, suggestive of less scatter from low H i column densities that favours the escape of Lyα photons. We infer stellar metallicities of Z/Z⊙ ≃ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines N v1240 Å and C iv1549 Å with those from stellar synthesis models. We also find a trend between MUV and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEsmore »show star formation rates up to ≃180 M⊙ yr−1, yet they appear relatively blue (βUV ≃ −1.8 to −2.0) showing evidence of little dust attenuation [E(B − V) = 0.10–0.14]. These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

    « less
  4. ABSTRACT We measure the Lyman continuum (LyC) escape fraction in 54 faint Lyman-alpha emitters (LAEs) at $z$ ≃ 3.1 in the GOODS-South field. With the average magnitude of R = 26.7 AB (MUV = −18.8 and L ≃ 0.1L*), these galaxies represent a population of compact young dwarf galaxies. Their properties are likely to resemble those in the galaxies responsible for reionizing the Universe at $z$ > 6. We do not detect LyC emission in any individual LAEs in the deep HST F336W images, which covers the rest-frame 820 Å. We do not detect the LyC emission of these LAEs in the stacked F336W images either. The 3σ upper limit of LyC escape fractions is $f_{\rm esc}\lt 14\!-\!32{{\ \rm per\ cent}}$. However, the high Ly α rest-frame equivalent width (EW), low stellar mass, and UV luminosity of these LAEs suggest that they should have $f_{\rm esc}\gt 50{{\ \rm per\ cent}}$. The low LyC escape fraction from this work and other stacking analyses suggests that the LyC-leaking galaxies with $f_{\rm esc}\gt 50{{\ \rm per\ cent}}$ at $z$ = 2–3 do not follow the relation between fesc and UV luminosity and Ly α EW derived from typical galaxies at similar redshifts. Therefore, the UV luminositymore »and Ly α EW are not the best indicators for the LyC escape fraction.« less
  5. 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 <  Δ vmore »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.« less