[Paper Review] 3D Lyman-alpha radiation transfer. I. Understanding Lyman-alpha line profile morphologies

Anne Verhamme, D. Schaerer+1 more|ArXiv.org|Aug 3, 2006

Photocathodes and Microchannel Plates266 citations

TL;DR

This paper presents a 3D Monte Carlo radiation transfer code for Lyman-alpha (Ly𝛼) photons in galaxies, accounting for arbitrary density, ionization, temperature, dust, and velocity structures. It reveals that in expanding shells with high column densities (NH ≳ 10²⁰ cm⁻²), the Ly𝛼 emission peak is systematically redshifted by approximately twice the expansion velocity—matching observations of z ≈ 3 Lyman break galaxies and confirming symmetric outflows as a key physical model.

ABSTRACT

Using a Monte Carlo technique, we have developed a 3D lyman-alpha radiation transfer code allowing for prescribed arbitrary hydrogen density, ionisation, temperature structures, and dust distribution, and arbitrary velocity fields and UV photon sources. We have examined the lyman-alpha line profiles predicted for several simple geometrical configurations and their dependence on the main input parameters. Overall, we find line profiles reaching from doubly peaked symmetric emission to symmetric Voigt (absorption) in static configurations with increasing dust content, and asymmetric red-(blue-) shifted emission lines with a blue (red) counterpart ranging from absorption to emission (with increasing line/continuum strength) in expanding (infalling) media. The following results are of interest for the interpretation of lya profiles from galaxies. 1/ Standard lya absorption line fitting of global spectra of galaxies may lead to an underestimate of the true hydrogen column density in certain geometrical conditions. 2/ Normal (inverted) P-Cygni like lya profiles can be obtained in expanding (infalling) media from objects without any intrinsic lya emission, as a natural consequence of radiation transfer effects. 3/ The formation and the detailed shape of lya profiles resulting from expanding shells has been thoroughly revised: for sufficiently large column densities, the position of the main lya emission peak is redshifted by twice the expansion velocity.This is in excellent agreement with the observations of z~3 LBGs. This finding indicates also that large scale, fairly symmetric shell structures must be a good description for the outflows in LBGs.(shortened abstract)

Motivation & Objective

• To develop a general 3D Ly𝛼 radiation transfer code capable of modeling complex astrophysical geometries and physical conditions in star-forming galaxies.
• To resolve discrepancies between observed Ly𝛼 line profiles and simplistic Voigt-profile fitting by incorporating full radiation transfer effects.
• To clarify the physical origin of asymmetric and P-Cygni-like Ly𝛼 profiles in outflowing or infalling media without intrinsic emission.
• To test how dust, column density, and kinematics shape emergent Ly𝛼 line morphologies in realistic 3D configurations.
• To provide a physically grounded tool for interpreting Ly𝛼 profiles in high-redshift galaxies and Lyman break galaxies (LBGs).

Proposed method

• A 3D Monte Carlo radiation transfer code is developed to simulate Ly𝛼 photon propagation through arbitrary 3D distributions of H I, dust, and velocity fields.
• The code handles arbitrary UV source spectra, including pure line emission, pure continuum, and intermediate cases with variable line-to-continuum strength.
• Radiation transfer includes resonant scattering, absorption, and dust attenuation, with full treatment of frequency redistribution and Doppler shifts.
• The model incorporates prescribed hydrogen density, ionization, temperature, and dust distributions, along with arbitrary velocity fields (e.g., expanding or infalling shells).
• Simulations are performed for static slabs, disks, externally illuminated slabs, and expanding/infalling shells to explore diverse morphologies.
• The code is validated against known limits and used to re-express earlier results, particularly on expanding shell profiles.

Experimental results

Research questions

• RQ1 How do 3D radiation transfer effects shape Ly𝛼 line profiles in expanding or infalling gas structures without intrinsic emission?
• RQ2 To what extent do standard Voigt-profile fits underestimate the true hydrogen column density in static or dusty systems?
• RQ3 What physical conditions produce the observed redshifted Ly𝛼 emission peaks in z ≈ 3 Lyman break galaxies?
• RQ4 How does the interplay between dust, column density, and kinematics transform Ly𝛼 line profiles from emission to absorption?
• RQ5 Can P-Cygni-like profiles arise purely from radiation transfer in outflowing media, even without intrinsic Ly𝛼 emission?

Key findings

• In static configurations with increasing dust content, Ly𝛼 line profiles evolve from symmetric double-peaked emission to symmetric Voigt-like absorption, with standard Voigt fitting underestimating the true hydrogen column density.
• Normal P-Cygni-like profiles can emerge in expanding media and inverted P-Cygni profiles in infalling media purely due to radiation transfer, even without intrinsic Ly𝛼 emission.
• For expanding shells with NH ≳ 10²⁰ cm⁻², the main Ly𝛼 emission peak is consistently redshifted by approximately twice the expansion velocity (2V_exp), matching observations of z ≈ 3 LBGs.
• This redshifted peak at 2V_exp strongly supports the presence of large-scale, symmetric outflowing shells in high-redshift Lyman break galaxies.
• The radiation transfer in expanding shells is fundamentally revised: the peak emission position is robustly shifted redward by 2V_exp for high column densities, independent of detailed velocity structure.
• The model demonstrates that degeneracies between column density, Doppler parameter, and intrinsic line width are real and must be resolved with full 3D radiation transfer, not simple profile fitting.