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

TL;DR

This work introduces the first analytic framework for steady-state radiative transfer across multiple resonance lines in a V-shaped atomic network, deriving closed-form solutions that reduce to a weighted superposition of single-line results in a transformed frequency coordinate. By generalizing the RTE, the authors define a multi-line Voigt profile and a frequency mapping that render the problem diffusion-dominated and tractable for arbitrary line counts, focusing on $N=2$. They validate the analytic solutions against a revised COLT Monte Carlo code and apply the method to Lyα fine structure and deuterium injection, finding that multi-line coupling largely preserves single-line behavior in typical conditions while enabling precise predictions of subtle spectral features. The results provide analytic benchmarks and a versatile numerical framework for multi-line resonant transfer, with implications for interpreting Lyα observations and modeling line-rich astrophysical environments, and they outline clear avenues for extending to more complex networks and additional physics such as recoil, velocity fields, and branching transitions.

Abstract

Resonance lines encode rich information about astrophysical sources and their environments, yet fully analytic treatments of multi-line radiative transfer remain almost entirely unexplored. We present exact, closed-form solutions for steady-state resonant-line radiative transfer in "V-shaped" atomic networks, where a single ground state couples to multiple transitions. Starting from the full angle-dependent transfer equation, we generalize absorption and emission coefficients to an arbitrary number of lines, derive a modified Fokker-Planck expansion of the frequency-redistribution COLT Monte Carlo radiative transfer code and find excellent agreement with the analytic predictions across a wide range of line separations, optical depths, and damping parameters, establishing our solutions as stringent validation benchmarks. For concrete applications related to the Lyman-alpha transition of neutral hydrogen, we examine how fine-structure splitting and deuterium injection modify the emergent spectra, internal radiation field, and radiative force multiplier. We show that these effects leave previous conclusions about Lyman-alpha feedback in the early universe essentially unchanged. Even when direct observational diagnostics are subtle, our framework provides novel analytic and numerical insights into coupled resonance-line transport and facilitates progress in general modeling of multi-line radiative transfer in diverse astrophysical settings.

Figures (10)

• Figure 1: Multi-line Voigt profile $\bar{H}(x)$ from Eq. (\ref{['eq:H']}) for a two-line system with asymmetric line strengths. The stronger line dominates except in the immediate vicinity of the weaker line. The normalization also affects the weak-line heights.
• Figure 2: A two-line V-shaped network with various optical depths and line separations. The asymmetry in chosen line parameters ($\omega_1 = 1/3$ and $\omega_2 = 2/3$) causes a difference between the peaks, with the second line being much stronger than the first. Upper Panel: Internal spectrum for a plane-parallel slab at $\tilde{r} = 0.5$, changing the optical depth over $\tau_0 \in \{10^4, 10^5, 10^6, 10^7\}$ at a temperature of $T = 10^4$ K, based on Eq. (\ref{['eq:PS-interior-solution-slab']}). Higher optical depth enables more diffusion across frequency space, since photons must typically reach a critical frequency of order $x_\text{esc} \sim (a \tau_0)^{1/3}$ before escape, increasing the likelihood of coupling between nearby lines. Lower Panel: Internal spectrum for a spherical cloud geometry with different line separations from the center as parametrized by $\Delta \in \pm \{0.5, 1.25, 2, 2.75\}\,(a\tau_0)^{1/3}$ at $\tilde{r} = 0.5$, based on Eq. (\ref{['eq:PS-interior-solution-sphere']}). As the spacing increases the solution separates into two distinct profiles, decoupling from each other. While there is still significant overlap the relative heights are similar due to mixing between lines.
• Figure 3: Internal spectrum $(\tilde{r} = 0.5)$ contributions of individual lines, i.e. $J_1$ and $J_2$ from Eq. (\ref{['eq:PS-interior-solution-sphere']}), in comparison to standard single-line solutions. We emphasize that the terms exhibit qualitatively different behavior as only the sum $J$ is physically meaningful across the full spectrum, rather than $J_1$ and $J_2$ in isolation.
• Figure 4: A toy network with two arbitrary lines of equal oscillator strength, illustrating the transition from well-separated lines with distinct peaks to overlapping lines with merged profiles as the frequency separation is reduced. We compare emergent spectra for slab geometry with different separations between the lines (top: widely separated; middle: intermediate; bottom: nearly overlapping). The analytic solution from Eq. (\ref{['eq:PS-emergent-solution-slab']}) and numerical colt simulations are in excellent agreement even as the lines become highly coupled to produce nontrivial spectra.
• Figure 5: A physical Ly$\alpha$ fine-structure setup where we vary the temperature of the system. Even for very low gas temperatures the two lines do not cleanly separate in the emergent spectrum. We compare emergent Ly$\alpha$ spectra for spherical geometry at $T = \{1, 10, 100\}\,\text{K}$. Decreasing the temperature reduces the Doppler width and causes the line separation to increase. Even at $1$ K, however, the two components overlap strongly and the profiles remain double-peaked rather than fully splitting. Atomic recoil becomes important at extremely low temperatures.