The discussion regarding the numerical integration of the polarized radiative transfer equation is still open and the comparison between the different numerical schemes proposed by different authors in the past is not fully clear. Aiming at facilitating the comprehension of the advantages and drawbacks of the different formal solvers, this work presents a reference paradigm for their characterization based on the concepts of order of accuracy, stability, and computational cost. Special attention is paid to understand the numerical methods belonging to the Diagonal Element Lambda Operator family, in an attempt to highlight their specificities.

Efficient numerical approximation of the polarized radiative transfer equation is challenging because this system of ordinary differential equations exhibits stiff behavior, which potentially results in numerical instability. This negatively impacts the accuracy of formal solvers, and small step-sizes are often necessary to retrieve physical solutions. This work presents stability analyses of formal solvers for the radiative transfer equation of polarized light, identifies instability issues, and suggests practical remedies. In particular, the assumptions and the limitations of the stability analysis of Runge-Kutta methods play a crucial role. On this basis, a suitable and pragmatic formal solver is outlined and tested. An insightful comparison to the scalar radiative transfer equation is also presented.

Context. The correct modeling of the scattering polarization signals observed in several strong resonance lines requires taking partial frequency redistribution (PRD) phenomena into account. Modeling scattering polarization with PRD effects is very computationally demanding and the simplifying angle-averaged (AA) approximation is therefore commonly applied. Aims. This work aims to assess the impact and the range of validity of the AA approximation with respect to the general angle-dependent (AD) treatment of PRD effects in the modeling of scattering polarization in strong resonance lines, with a focus on the solar Ca I 4227 Å line. Methods. Spectral line polarization was modeled by solving the radiative transfer problem for polarized radiation, under nonlocal thermodynamic equilibrium conditions, taking PRD effects into account in static one-dimensional semi-empirical atmospheric models presenting arbitrary magnetic fields. The problem was solved through a two-step approach. In step 1, the problem was solved for the intensity only, considering a multilevel atom. In step 2, the problem was solved including polarization, considering a two-level atom with an unpolarized and infinitely sharp lower level, and fixing the lower level population calculated at step 1. Results. The results for the Ca I 4227 Å line show a good agreement between the AA and AD calculations for the Q/I and U/I wings’ signals. However, AA calculations reveal an artificial trough in the line-core peak of the linear polarization profiles, whereas AD calculations show a sharper peak in agreement with the observations. Conclusions. An AD treatment of PRD effects is essential to correctly model the line-core peak of the scattering polarization signal of the Ca I 4227 Å line. By contrast, in the considered static case, the AA approximation seems to be suitable to model the wing scattering polarization lobes and their magnetic sensitivity through magneto-optical effects.

Context. The numerical modeling of the generation and transfer of polarized radiation is a key task in solar and stellar physics research and has led to a relevant class of discrete problems that can be reframed as linear systems. In order to solve such problems, it is common to rely on efficient stationary iterative methods. However, the convergence properties of these methods are problemdependent, and a rigorous investigation of their convergence conditions, when applied to transfer problems of polarized radiation, is still lacking. Aims. After summarizing the most widely employed iterative methods used in the numerical transfer of polarized radiation, this article aims to clarify how the convergence of these methods depends on different design elements, such as the choice of the formal solver, the discretization of the problem, or the use of damping factors. The main goal is to highlight advantages and disadvantages of the different iterative methods in terms of stability and rate of convergence. Methods. We first introduce an algebraic formulation of the radiative transfer problem. This formulation allows us to explicitly assemble the iteration matrices arising from different stationary iterative methods, compute their spectral radii and derive their convergence rates, and test the impact of different discretization settings, problem parameters, and damping factors. Results. Numerical analysis shows that the choice of the formal solver significantly affects, and can even prevent, the convergence of an iterative method. Moreover, the use of a suitable damping factor can both enforce stability and increase the convergence rate. Conclusions. The general methodology used in this article, based on a fully algebraic formulation of linear transfer problems of polarized radiation, provides useful estimates of the convergence rates of various iterative schemes. Additionally, it can lead to novel solution approaches as well as analyses for a wider range of settings, including the unpolarized case.