skip to main content


Title: Deciphering the Lyman- α emission line: towards the understanding of galactic properties extracted from Ly α spectra via radiative transfer modelling
ABSTRACT Existing ubiquitously in the Universe with the highest luminosity, the Lyman-α (Lyα) emission line encodes abundant physical information about the gaseous medium it interacts with. Nevertheless, the resonant nature of the Lyα line complicates the radiative transfer (RT) modelling of the line profile. We revisit the problem of deciphering the Lyα emission line with RT modelling. We reveal intrinsic parameter degeneracies in the widely used shell model in the optically thick regime for both static and outflowing cases, which suggest the limitations of the model. We also explore the connection between the more physically realistic multiphase, clumpy model, and the shell model. We find that the parameters of a ‘very clumpy’ slab model and the shell model have the following correspondences: (1) the total column density, the effective temperature, and the average radial clump outflow velocity of the clumpy slab model are equal to the H i column density, effective temperature, and expansion velocity of the shell model, respectively; (2) large intrinsic linewidths are required in the shell model to reproduce the wings of the clumpy slab models; (3) adding another phase of hot interclump medium increases peak separation, and the fitted shell expansion velocity lies between the outflow velocities of two phases of gas. Our results provide a viable solution to the major discrepancies associated with Lyα fitting reported in previous literature, and emphasize the importance of utilizing information from additional observations to break the intrinsic degeneracies and interpreting the model parameters in a more physically realistic context.  more » « less
Award ID(s):
2009278
NSF-PAR ID:
10340936
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
513
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5034 to 5051
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We present new observations of Lyman-α (Lyα) Blob 1 (LAB1) in the SSA22 protocluster region (z = 3.09) using the Keck Cosmic Web Imager and Keck Multi-object Spectrometer for Infrared Exploration. We have created a narrow-band Lyα image and identified several prominent features. By comparing the spatial distributions and intensities of Lyα and Hβ, we find that recombination of photo-ionized H i gas followed by resonant scattering is sufficient to explain all the observed Lyα/Hβ ratios. We further decode the spatially resolved Lyα profiles using both moment maps and radiative transfer modelling. By fitting a set of multiphase, ‘clumpy’ models to the observed Lyα profiles, we manage to reasonably constrain many parameters, namely the H i number density in the interclump medium (ICM), the cloud volume filling factor, the random velocity and outflow velocity of the clumps, the H i outflow velocity of the ICM, and the local systemic redshift. Our model has successfully reproduced the diverse Lyα morphologies, and the main results are: (1) the observed Lyα spectra require relatively few clumps per line of sight as they have significant fluxes at the line centre; (2) the velocity dispersion of the clumps yields a significant broadening of the spectra as observed; (3) the clump bulk outflow can also cause additional broadening if the H i in the ICM is optically thick; (4) and the H i in the ICM is responsible for the absorption feature close to the Lyα line centre. 
    more » « less
  2. ABSTRACT

    We present a new semi-analytical formalism for modelling metal absorption lines that emerge from a clumpy galactic environment, ALPACA. We predict the “down-the-barrel” (DTB) metal absorption line profiles and the equivalent width (EW) of absorption at different impact parameters (b) as a function of the clump properties, including clump kinematics, clump volume filling factor, clump number density profile, and clump ion column densities. With ALPACA, we jointly model the stacked DTB C ii λ1334 spectrum of a sample of z ∼ 3 Lyman break galaxies and the EW versus b profile of a sample of z ∼ 2 star-forming galaxy–galaxy pairs. ALPACA successfully reproduced two data sets simultaneously, and the best fit prefers a low clump volume filling factor (∼3 × 10−3). The radial velocities of the clumps are a superposition of a rapidly accelerated outflow with a maximum velocity of $\sim 400 \, {\mathrm{km}\, \mathrm{s}^{-1}}$ and a velocity dispersion of $\sigma \sim 120 \, {\mathrm{km}\, \mathrm{s}^{-1}}$. The joint modelling reveals a physical scenario where the absorption observed at a particular velocity is contributed by the clumps distributed over a fairly broad range of radii. We also find that the commonly adopted Sobolev approximation is at best only applicable within a narrow range of radii where the clumps are undergoing rapid acceleration in a non-volume-filling clumpy medium. Lastly, we find that the clump radial velocity profile may not be fully constrained by the joint modelling and spatially resolved Ly α emission modelling may help break the degeneracy.

     
    more » « less
  3. 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 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.

     
    more » « less
  4. 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
  5. Abstract

    The resonantly scattered Lyαline illuminates the extended halos of neutral hydrogen in the circumgalactic medium of galaxies. We present integral field Keck Cosmic Web Imager observations of double-peaked, spatially extended Lyαemission in 12 relatively low-mass (M∼ 109M)z∼ 2 galaxies characterized by extreme nebular emission lines. Using individual spaxels and small bins as well as radially binned profiles of larger regions, we find that for most objects in the sample the Lyαblue-to-red peak ratio increases, the peak separation decreases, and the fraction of flux emerging at line center increases with radius. We use new radiative transfer simulations to model each galaxy with a clumpy, multiphase outflow with radially varying outflow velocity, and self-consistently apply the same velocity model to the low-ionization interstellar absorption lines. These models reproduce the trends of peak ratio, peak separation, and trough depth with radius, and broadly reconcile outflow velocities inferred from Lyαand absorption lines. The galaxies in our sample are well-described by a model in which neutral, outflowing clumps are embedded in a hotter, more highly ionized inter-clump medium (ICM), whose residual neutral content produces absorption at the systemic redshift. The peak ratio, peak separation, and trough flux fraction are primarily governed by the line-of-sight component of the outflow velocity, the Hicolumn density, and the residual neutral density in the ICM respectively. The azimuthal asymmetries in the line profile further suggest nonradial gas motions at large radii and variations in the Hicolumn density in the outer halos.

     
    more » « less