Abstract:Making the conversion from the geometrical spatial scale to the optical depth spatial scale is useful in obtaining numerical solutions for the radiative transfer equation. This is because it allows for the use of exponential integrators, while enforcing numerical stability. Such a conversion involves the integration of the total opacity of the medium along the considered ray path. This is usually approximated by applying a piecewise quadrature in each spatial cell of the discretized medium. However, a rigorous analysis of this numerical step is lacking. This work is aimed at clearly assessing the performance of different optical depth conversion schemes with respect to the solution of the radiative transfer problem for polarized radiation, out of the local thermodynamic equilibrium. We analyzed different optical depth conversion schemes and their combinations with common formal solvers, both in terms of the rate of convergence as a function of the number of spatial points and the accuracy of the emergent Stokes profiles. The analysis was performed in a 1D semi-empirical model of the solar atmosphere, both in the absence and in the presence of a magnetic field. We solved the transfer problem of polarized radiation in different settings: the continuum, the photospheric Sr i AA modeled under the assumption of complete frequency redistribution, and the chromospheric Ca i AA taking the partial frequency redistribution effects into account during the modeling. High-order conversion schemes generally outperform low-order methods when a sufficiently high number of spatial grid points is considered. In the synthesis of the emergent Stokes profiles, the convergence rate, as a function of the number of spatial points, is impacted by both the conversion scheme and formal solver. The use of low-order conversion schemes significantly reduces the accuracy of high-order formal solvers. In practical applications, the use of low-order optical depth conversion schemes introduces large numerical errors in the formal solution. To fully exploit high-order formal solvers and obtain accurate synthetic emergent Stokes profiles, it is necessary to use high-order optical depth conversion schemes.

Numerical Simulation of Polarized Light and Temperature with a Refractive Interface

Polarized Radiative Transfer in an Arbitrary Multilayer Semitransparent Medium

Theoretical Extension of Universal Forward and Backward Monte Carlo Radiative Transfer Modeling for Passive and Active Polarization Observation Simulations

Monte Carlo method for polarized radiative transfer in gradient-index media

Polarized Radiative Transfer Considering Thermal Emission in Semitransparent Media

Approximate Solution to Vector Radiative Transfer in Gradient-Index Medium

Analysis of Radiative Heat Transfer in Two-Dimensional Semitransparent Medium with Piece-Wise Constant Refractive Index

Boundary conditions for NLTE polarized radiative transfer with incident radiation

Transient/Time-Dependent Radiative Transfer in A Two-Dimensional Scattering Medium Considering the Polarization Effect

Vector Radiative Transfer in a Vertically Inhomogeneous Scattering and Emitting Atmosphere. Part I: A New Discrete Ordinate Method

Transient Radiative Transfer in A Scattering Slab Considering Polarization

Impact of the numerical conversion to optical depth on the transfer of polarized radiation

Vector Radiative Transfer in A Multilayer Medium by Natural Element Method

Ray-Tracing/Nodal-Analyzing Model for Transient Thermal Behavior of A Scattering Medium with A Variable Refractive Index

Transient Radiative Heat Transfer in an Inhomogeneous Participating Medium with Fresnel’s Surfaces

Polarized Radiative Transfer Simulations: A Tutorial Review and Upgrades of the Vector Discrete Ordinate Radiative Transfer Computational Tool

Transient Coupled Heat Transfer Inside a Scattering Mmedium with Graded Refractive Index

On the Derivation of Vector Radiative Transfer Equation for Polarized Radiative Transport in Graded Index Media

Analytical ray-tracing in planetary atmospheres

Analysis of the Rigorous Coupled Wave Approach for p-polarized light in gratings

Polarized radiative transfer in two-dimensional scattering medium with complex geometries by natural element method