3D radiative transfer modeling of scattering polarization with partial frequency redistribution I. Verification and disk-center results for the solar Ca I 4227 {\AA} line

Abstract

Several strong solar resonance lines show observable linear scattering polarization signals, holding a great potential for investigating the magnetism of the outer solar atmosphere. Accurately modeling these signals requires solving the radiative transfer (RT) problem for polarized radiation in comprehensive 3D models of the solar atmosphere, in non-local thermodynamic equilibrium, accounting for partial frequency redistribution (PRD) effects. This problem has so far been computationally inaccessible. We present the first scientific application of TRIP, a novel software for the massively parallel solution of the 3D non-LTE RT problem for polarized radiation, including scattering polarization and PRD. We aim to verify the code and explore the combined action of PRD and the 3D structure of the solar atmosphere on scattering polarization. We run TRIP to synthesize the Stokes profiles of the Ca I line at 4227 {\AA} in a 3D model of the solar atmosphere extracted from a radiation magneto-hydrodynamic simulation. We efficiently solve the resulting large-scale problem, with up to $4 \times 10^{10}$ degrees of freedom, with a state-of-the-art preconditioned Krylov method, using up to 20 thousand parallel CPUs. After including verification tests, we find that the joint impact of PRD effects and the detailed 3D structure of the atmospheric model produce disk-center scattering polarization signals in the line wings. These signals are sensitive to the magnetic field, via magneto-optical effects, and to bulk velocity gradients. We also show that the CRD approximation underestimates the amplitude of disk-center line-core signals. This achievement represents a crucial step forward for diagnosing the magnetism of the solar chromosphere and transition region through the quantitative comparisons of synthetic and observational data.