Analytical and fitting formulae for solutions to Lyman-alpha radiative transfer equations: the effects of geometry, recoil, and velocity gradients
arXiv:2606.27423 · doi:10.1093/mnras/stag1223
Abstract
Lyman-alpha (Ly) radiative transfer (RT) is important in many astrophysical environments and governed by multiple physical processes. In this paper, we provide analytical formulae/procedures for the solutions to Ly RT equations under three simple geometrical symmetries and investigate the effects of atomic recoil and gas bulk motion. We first study Ly spectra by solving Ly RT equations for a static, uniform gas cloud under cylindrical geometry. The solution is verified through Ly Monte Carlo RT simulations, and compared to those under slab and spherical geometries in literature. Second, to characterise the recoil effect, we empirically modify recoil-free Ly spectra. The method is motivated by Ly RT equations with recoil and justified by simulations. Finally, we account for constant velocity gradients in Ly RT equations and obtain series solutions for Ly spectra. The solutions demonstrate good agreement to Ly spectra from simulations for small velocity gradients (i.e. edge velocity of a cloud being comparable to the thermal velocity ) but become less accurate for large ones. To characterise Ly spectra under large velocity gradients, we empirically extend the functional form of solutions and constrain them from fitting simulated Ly spectra. The resulting fitting formulae show significant improvement for large velocity gradients () under large optical depths. The analytical study of Ly spectra in this work completes the set of solutions under simple geometries, provides physical insights for Ly RT under recoil and velocity gradient, and develops analytical tools for theoretical studies that require inputs from Ly RT.
20 pages, 9 figures, 2 tables, 2 github repositories: https://github.com/zhengzheng-astro/RandomGenerator and https://github.com/PengfeiLiAstro/LyaRTAnalytical, accepted by MNRAS