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

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.