Full non-LTE multi-level radiative transfer I. An atom with three bound infinitely sharp levels
Tristan Lagache, Frédéric Paletou, Malali Sampoorna
TL;DR
The paper tackles the full non-LTE (FNLTE) radiative-transfer problem for a multi-level atom while allowing deviations of the velocity distribution of massive particles from Maxwellian. It develops a numerical scheme that iterates the radiative-transfer equation, kinetic-equilibrium equations, and Boltzmann equations for the velocity distributions, starting from a MALI-CRD initialization and using a three-level, infinitely sharp- level atom. It validates the method against standard two-level PRD, multi-level CRD, and cross-redistribution XRD benchmarks, and then presents new FNLTE results showing non-Maxwellian features in the excited-level velocity distributions and the impact on emergent line profiles. The work highlights the conceptual and practical importance of fully coupled FNLTE treatment and discusses limitations and directions for future improvements, including extending broadening, additional levels, and computational efficiency.
Abstract
The standard nonlocal thermodynamic equilibrium (non-LTE) multi-level radiative transfer problem only takes into account the deviation of the radiation field and atomic populations from their equilibrium distribution. We aim to show how to solve for the full non-LTE (FNLTE) multi-level radiative transfer problem, also accounting for deviation of the velocity distribution of the massive particles from Maxwellian. We considered, as a first step, a three-level atom with zero natural broadening. In this work, we present a new numerical scheme. Its initialisation relies on the classic, multi-level approximate Lambda-iteration (MALI) method for the standard non-LTE problem. The radiative transfer equations, the kinetic equilibrium equations for atomic populations, and the Boltzmann equations for the velocity distribution functions were simultaneously iterated in order to obtain self-consistent particle distributions. During the process, the observer's frame absorption and emission profiles were re-computed at every iterative step by convolving the atomic frame quantities with the relevant velocity distribution function. We validate our numerical strategy by comparing our results with the standard non-LTE solutions in the limit of a two-level atom with Hummer's partial redistribution in frequency, and with a three-level atom with complete redistribution. In this work, we considered the so-called cross-redistribution problem. We then show new FNLTE results for a simple three-level atom while evaluating the assumptions made for the emission and absorption profiles of the standard non-LTE problem with partial and cross-redistribution.
