Monte Carlo Simulation of Lyman Alpha Scattering and Application to Damped Lyman Alpha Systems

Zheng Zheng,Jordi Miralda-Escude
DOI: https://doi.org/10.1086/342400
2002-06-26
Abstract:A Monte Carlo code to solve the transfer of Lyman alpha (Lya) photons is developed, which can predict the Lya image and two-dimensional Lya spectra of a hydrogen cloud with any given geometry, Lya emissivity, neutral hydrogen density distribution, and bulk velocity field. We apply the code to several simple cases of a uniform cloud to show how the Lya image and emitted line spectrum are affected by the column density, internal velocity gradients, and emissivity distribution. We then apply the code to two models for damped Lya absorption systems: a spherical, static, isothermal cloud, and a flattened, axially symmetric, rotating cloud. If the emission is due to fluorescence of the external background radiation, the Lya image should have a core corresponding to the region where hydrogen is self-shielded. The emission line profile has the characteristic double peak with a deep central trough. We show how rotation of the cloud causes the two peaks to shift in wavelength as the slit is perpendicular to the rotation axis, and how the relative amplitude of the two peaks is changed. In reality, damped Lya systems are likely to have a clumpy gas distribution with turbulent velocity fields, which should smooth the line emission profile, but should still leave the rotation signature of the wavelength shift across the system.
Astrophysics
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to study the Lyman - α (Lyα for short) scattering of hydrogen atoms by the Monte Carlo simulation method and apply it to Damped Lyman - α Systems (DLAs). Specifically, the author has developed a Monte Carlo code to calculate the transmission of Lyα photons under any geometric shape, emissivity distribution, neutral hydrogen density distribution and overall velocity field. This code can predict the two - dimensional Lyα image and spectral line in a given direction. #### Main research problems include: 1. **Understanding the properties of high - column - density absorption systems**: - High - column - density absorption systems (such as DLAs and Lyman limit systems) are a crucial part of understanding galaxy formation. These systems may represent the gas halos in galaxies in the process of formation, rather than fully - formed spiral disks. - Studying the internal structure of these systems can help reveal the specific processes of galaxy formation. 2. **Studying the characteristics of Lyα emission**: - Lyα emission is produced by the recombination of hydrogen atoms and then scattered through the gas. This emission can provide information about the internal structure and dynamic state of gas clouds. - By simulating Lyα emission under different conditions, the author has explored the effects of column density, internal velocity gradient and emissivity distribution on Lyα images and emission - line spectra. 3. **Verifying model assumptions**: - The paper has verified the accuracy of the code by comparing the Monte Carlo simulation results with analytical solutions. - Two different models of DLAs have been simulated: static isothermal spherical clouds and flat axisymmetric rotating clouds, to study the effect of rotation on Lyα emission. 4. **Exploring applications in actual observations**: - By simulating different physical conditions (such as temperature, velocity dispersion, etc.), the author has shown how to detect the rotation and other dynamic features of DLAs through Lyα emission. - Research has shown that even if there may be inhomogeneous gas distributions and turbulent velocity fields in actual DLAs, the rotation features can still be manifested through wavelength shifts. #### Summary: The main objective of this paper is to study the transmission behavior of Lyα photons under different conditions through Monte Carlo simulation, thereby providing a new perspective for understanding the properties of high - column - density absorption systems and their roles in the galaxy formation process.