skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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 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
  3. 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
  4. 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
  5. ABSTRACT The observability of Lyα emitting galaxies (LAEs) during the Epoch of Reionization can provide a sensitive probe of the evolving neutral hydrogen gas distribution, thus setting valuable constraints to distinguish different reionization models. In this study, we utilize the new thesan suite of large-volume ($$L_\text{box} = 95.5\, \text{cMpc}$$) cosmological radiation-hydrodynamic simulations to directly model the Lyα emission from individual galaxies and the subsequent transmission through the intergalactic medium. thesan combines the arepo-rt radiation-hydrodynamic solver with the IllustrisTNG galaxy formation model and includes high- and medium-resolution simulations designed to investigate the impacts of halo-mass-dependent escape fractions, alternative dark matter models, and numerical convergence. We find important differences in the Lyα transmission based on reionization history, bubble morphology, frequency offset from line centre, and galaxy brightness. For a given global neutral fraction, Lyα transmission reduces when low-mass haloes dominate reionization over high-mass haloes. Furthermore, the variation across sightlines for a single galaxy is greater than the variation across all galaxies. This collectively affects the visibility of LAEs, directly impacting observed Lyα luminosity functions (LFs). We employ Gaussian Process Regression using SWIFTEmulator to rapidly constrain an empirical model for dust escape fractions and emergent spectral-line profiles to match observed LFs. We find that dust strongly impacts the Lyα transmission and covering fractions of MUV ≲ −19 galaxies in $$M_\text{vir} \gtrsim 10^{11}\, \text{M}_{\bigodot }$$ haloes, such that the dominant mode of removing Lyα photons in non-LAEs changes from low-IGM transmission to high dust absorption around z ∼ 7. 
    more » « less