Abstract:The calculation of the emerging radiation from a model atmosphere requires knowledge of the emissivity and absorption coefficients, which are proportional to the atomic level population densities of the levels involved in each transition. Due to the intricate interdependency of the radiation field and the physical state of the atoms, iterative methods are required in order to calculate the atomic level population densities. A variety of different methods have been proposed to solve this problem, which is known as the Non-Local Thermodynamical Equilibrium (NLTE) problem. In this study we have developed a Jacobian-Free Newton-Krylov method (JFNK) to solve multi-level NLTE radiative transfer problems. Using the Rybicki & Hummer (1992) method as a reference (Rybicki, G. B. & Hummer, D. G. 1992, A&A, 262, 209), our results show that our JFNK solver can achieve up to a factor two speed up when using local approximate operators / preconditioners, while also achieving a lower residual error in the statistical equilibrium equations. Another advantage of this method is that the addition of charge conservation and partial redistribution effects should be straight forward. Our method can help accelerating the calculation of the emerging spectra from numerical models and also the reconstruction of chromospheric datasets through NLTE inversions.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to calculate the atomic level population density more efficiently and rapidly when solving the radiative transfer problem in model atmospheres under non - local thermodynamic equilibrium (NLTE) conditions. Specifically: 1. **Background**: - When calculating the radiation emitted from the model atmosphere, the emissivity and absorption coefficient need to be known, and these quantities are proportional to the atomic level population density involved in the transition. - Due to the complex interdependence between the radiation field and the atomic physical state, an iterative method must be used to calculate the atomic level population density. 2. **Objectives**: - The objective of the paper is to develop an efficient and rapidly convergent method under the assumption of statistical equilibrium to solve the NLTE problem. - In particular, to explore the applicability of the Jacobian - Free Newton - Krylov (JFNK) method, which estimates new correction terms by the Krylov subspace method without the need to explicitly construct the Jacobian matrix. 3. **Limitations of existing methods**: - Early attempts used the Lambda iteration method, but its convergence was very poor and it was not feasible in practical applications. - Other methods such as linearization and operator splitting techniques have certain effects, but the convergence speed is slow in the case of high scattering, and it is difficult to handle complex overlapping transitions and other problems. 4. **Advantages of the new method**: - The JFNK method does not depend on the operator splitting of the local approximation operator, so it can improve the convergence of the NLTE problem in highly scattering cases. - This method allows for the direct implementation of more complex problems, such as overlapping transitions of different active atoms, charge conservation, or more effective partial redistribution processing, without the need to explicitly linearize equations. In summary, this paper aims to provide a new solution to more efficiently solve the radiative transfer problem under NLTE conditions by introducing the JFNK method.

Jacobian-Free Newton-Krylov method for multilevel NLTE radiative transfer problems