Multi-Resonant-Line Radiative Transfer: Lyman-Alpha Fine Structure and Deuterium Coupling

Resonance lines encode important information about astrophysical sources and their environments, and the Lyman-alpha (Lya) transition of neutral hydrogen is especially powerful in this regard. We present exact, closed-form analytic solutions for steady-state resonant-line radiative transfer in "V-shaped" networks, i.e. media where a single ground state couples to multiple transitions, such as Lya fine-structure and deuterium scattering. Starting from the full angle-dependent transfer equation, we combine (i) moment closure in the optically thick limit and (ii) a modified Fokker–Planck expansion of the frequency-redistribution integral to obtain a single partial differential equation governing the mean intensity. A judicious change of variables reduces this equation to a free-space Helmholtz problem with point-like sources in frequency space, admitting analytic solutions for arbitrary sets of lines with fixed frequency offsets and strengths. We implement coupled V-shaped line networks in the COLT Monte Carlo radiative transfer code and find excellent agreement with our analytic predictions across line separations, optical depths, and damping parameters. Modern Lya codes increasingly include fine-structure and deuterium physics, which can shape the emergent spectra and internal radiation field. Our framework provides novel analytic and numerical insights, rigorous validation benchmarks, and facilitates progress in general modeling of multi-lined radiative transfer in astrophysics.