Modeling the scattering polarization in the solar Ca I 4227Å line with angle-dependent PRD effects and bulk velocities

Abstract

Modeling the scattering polarization signals of strong chromospheric lines requires solving the radiative transfer (RT) problem for polarized radiation, out of local thermodynamic equilibrium, taking partial frequency redistribution (PRD) effects into account. This problem is extremely challenging from a computational standpoint and, so far, most studies have been carried out by either modeling PRD effects under the angle-average approximation or by considering academic models of the solar atmosphere. Thanks to a new solution strategy, applicable to atomic systems that allow for a linearization of the problem, accurate solutions can now be routinely obtained in realistic 1D models, taking angle-dependent (AD) PRD effects into account. This work is aimed at assessing the suitability and performance of this new approach to handling dynamic scenarios. At the same time, it aims to explore the joint impact of magnetic fields and bulk velocities on the scattering polarization profiles of strong resonance lines, accounting for AD PRD effects and considering more realistic atmospheric models than in previous investigations. By using a two-level atomic model for neutral calcium, we synthesized the intensity and polarization profiles of the Ca i AA line. Our calculations were performed in 1D atmospheric models, both semi-empirical and extracted from 3D MHD simulations, including vertical bulk velocities and magnetic fields of arbitrary strength and orientation, both constant and varying with height. We obtained accurate solutions after only a few iterations across all considered scenarios. Even when formulating the problem in the observer's reference frame, the frequency and angular grids required for accurate results were easily manageable. The calculated profiles showed the expected signatures of bulk velocities: wavelength shifts, enhancement of the line-core polarization amplitude, and prominent asymmetries in the wing signals. The results obtained in atmospheric models with complex thermal, dynamic, and magnetic structures unveiled the broad diversity of features in the emergent radiation that can be expected from realistic scenarios. The presented results assess the suitability of the proposed solution strategy and its parallel implementation, thus supporting its generalization to the 3D case. Our applications in increasingly realistic atmospheric models showed the difficulty related to precisely establishing the individual weight of bulk velocities and magnetic fields in the shape of the emergent profiles. This highlights the need to account for both these physical ingredients to perform reliable inversions of observed scattering polarization profiles.