A numerical approach for modelling the polarisation signals of strong resonance lines with partial frequency redistribution. Numerical applications to two-term atoms and plane-parallel atmospheres
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866913942027108352 |
|---|---|
| author | Riva, Fabio Janett, Gioele Belluzzi, Luca Alemán, Tanausú del Pino Ballester, Ernest Alsina Bueno, Javier Trujillo Benedusi, Pietro Riva, Simone Krause, Rolf |
| author_facet | Riva, Fabio Janett, Gioele Belluzzi, Luca Alemán, Tanausú del Pino Ballester, Ernest Alsina Bueno, Javier Trujillo Benedusi, Pietro Riva, Simone Krause, Rolf |
| contents | Aims. The main goal of this paper is to present an accurate and efficient numerical strategy for solving the radiative transfer problem for polarised radiation in strong resonance lines forming out of local thermodynamic equilibrium, taking angle-dependent (AD) partial frequency redistribution (PRD) effects and J-state interference into account. We consider the polarisation produced both by the Zeeman effect and by the scattering of anisotropic radiation, along with its sensitivity to the Hanle and magneto-optical effects. Methods. We introduce a formalism that allows treating both a two-level and a two-term atom in the presence of arbitrary magnetic and bulk velocity fields. The problem is formulated by treating the population of the lower level/term as a fixed input parameter. This approach makes the problem linear with respect to the radiation field, enabling the application of efficient matrix-free preconditioned iterative methods for its solution. Additionally, the computation of the scattering emissivity in the comoving frame, together with a careful choice of the angular and spectral quadrature nodes, allow us to speed up the calculations by reducing the number of evaluations of the redistribution functions. Results. The proposed solution strategy is applied to synthesise the Stokes profiles of the Mg ii h&k doublet and the H i Ly-α line in 1D semi-empirical models. The results demonstrate that the method is both fast and accurate. A comparison with calculations from HanleRT-TIC displays an overall good agreement, thereby validating our solution strategy. Moreover, for the wavelength-integrated polarisation profiles of the H i Ly-α line, we find an excellent agreement between the results obtained including PRD effects in their general AD description and those obtained considering the angle-averaged simplifying approximation. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2505_20968 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A numerical approach for modelling the polarisation signals of strong resonance lines with partial frequency redistribution. Numerical applications to two-term atoms and plane-parallel atmospheres Riva, Fabio Janett, Gioele Belluzzi, Luca Alemán, Tanausú del Pino Ballester, Ernest Alsina Bueno, Javier Trujillo Benedusi, Pietro Riva, Simone Krause, Rolf Solar and Stellar Astrophysics Instrumentation and Methods for Astrophysics Aims. The main goal of this paper is to present an accurate and efficient numerical strategy for solving the radiative transfer problem for polarised radiation in strong resonance lines forming out of local thermodynamic equilibrium, taking angle-dependent (AD) partial frequency redistribution (PRD) effects and J-state interference into account. We consider the polarisation produced both by the Zeeman effect and by the scattering of anisotropic radiation, along with its sensitivity to the Hanle and magneto-optical effects. Methods. We introduce a formalism that allows treating both a two-level and a two-term atom in the presence of arbitrary magnetic and bulk velocity fields. The problem is formulated by treating the population of the lower level/term as a fixed input parameter. This approach makes the problem linear with respect to the radiation field, enabling the application of efficient matrix-free preconditioned iterative methods for its solution. Additionally, the computation of the scattering emissivity in the comoving frame, together with a careful choice of the angular and spectral quadrature nodes, allow us to speed up the calculations by reducing the number of evaluations of the redistribution functions. Results. The proposed solution strategy is applied to synthesise the Stokes profiles of the Mg ii h&k doublet and the H i Ly-α line in 1D semi-empirical models. The results demonstrate that the method is both fast and accurate. A comparison with calculations from HanleRT-TIC displays an overall good agreement, thereby validating our solution strategy. Moreover, for the wavelength-integrated polarisation profiles of the H i Ly-α line, we find an excellent agreement between the results obtained including PRD effects in their general AD description and those obtained considering the angle-averaged simplifying approximation. |
| title | A numerical approach for modelling the polarisation signals of strong resonance lines with partial frequency redistribution. Numerical applications to two-term atoms and plane-parallel atmospheres |
| topic | Solar and Stellar Astrophysics Instrumentation and Methods for Astrophysics |
| url | https://arxiv.org/abs/2505.20968 |