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 polarization signals produced by the scattering of anistropic radiation in strong resonance lines encode important information about the elusive magnetic fields in the outer layers of the solar atmosphere. An accurate modeling of these signals is a very challenging problem from the computational point of view, in particular when partial frequency redistribution (PRD) effects in scattering processes are accounted for with a general angle-dependent treatment. Aims. We aim at solving the radiative transfer problem for polarized radiation in nonlocal thermodynamic equilibrium conditions, taking angle-dependent PRD effects into account. The problem is formulated for a two-level atomic model in the presence of arbitrary magnetic and bulk velocity fields. The polarization produced by scattering processes and the Zeeman effect is considered. Methods. The proposed solution strategy is based on an algebraic formulation of the problem and relies on a convenient physical assumption, which allows its linearization. We applied a nested matrix-free GMRES iterative method. Effective preconditioning is obtained in a multifidelity framework by considering the light-weight description of scattering processes in the limit of complete frequency redistribution (CRD). Results. Numerical experiments for a one-dimensional (1D) atmospheric model show near optimal strong and weak scaling of the proposed CRD-preconditioned GMRES method, which converges in few iterations, independently of the discretization parameters. A suitable parallelization strategy and high-performance computing tools lead to competitive run times, providing accurate solutions in a few minutes. Conclusions. The proposed solution strategy allows the fast systematic modeling of the scattering polarization signals of strong resonance lines, taking angle-dependent PRD effects into account together with the impact of arbitrary magnetic and bulk velocity fields. Almost optimal strong and weak scaling results suggest that this strategy is applicable to realistic 3D models. Moreover, the proposed strategy is general, and applications to more complex atomic models are possible.

No abstract