[Paper Review] Generalized sqrt(epsilon)-law. The role of unphysical source terms in resonance line polarization transfer and its importance as an additional test of NLTE radiative transfer codes

Jiří Štěpán, V. Bommier|ArXiv.org|Apr 12, 2007

Ionosphere and magnetosphere dynamics9 references3 citations

TL;DR

This paper generalizes the $√\epsilon$-law for non-LTE radiative transfer by introducing unphysical source terms into the statistical equilibrium equations for polarization transfer in a two-level atom under non-thermal collisional excitation. By artificially suppressing population excitation while enhancing alignment, the method creates a test case where polarization dominates over intensity, enabling sensitive validation of NLTE radiative transfer codes through an analytical $√\epsilon$-law for higher-rank source function components.

ABSTRACT

Context. A derivation of a generalized sqrt(epsilon)-law for nonthermal collisional rates of excitation by charged perturbers is presented. Aims. Aim of this paper is to find a more general analytical expression for a surface value of the source function which can be used as an addtional tool for verification of the non-LTE radiative transfer codes. Methods. Under the impact approximation hypothesis, static, one-dimensional, plane-parallel atmosphere, constant magnetic field of arbitrary strength and direction, two-level atom model with unpolarized lower level and stimulated emission neglected, we introduce the unphysical terms into the equations of statistical equilibrium and solve the appropriate non-LTE integral equations. Results. We derive a new analytical condition for the surface values of the source function components expressed in the basis of irreducible spherical tensors.

Motivation & Objective

• To develop a more general analytical condition for surface source function components in non-LTE radiative transfer, particularly for polarized radiation.
• To introduce unphysical source terms into the statistical equilibrium equations to create extreme test cases for NLTE radiative transfer codes.
• To enable sensitive verification of polarization transfer accuracy in codes by favoring higher-rank tensor components over intensity.
• To extend the $√\epsilon$-law to non-thermal collisional rates and arbitrary magnetic fields, enhancing its utility as a benchmark.

Proposed method

• Adopts the irreducible spherical tensor formalism for the density matrix, using $\rho^{K}_{Q}$ components to describe atomic level populations and polarizations.
• Introduces unphysical collisional rates by swapping excitation contributions between $K=0$ (population) and $K=2, Q=0$ (alignment), setting $(C_{jj'}^{0})_{00} \to 0$ and $(C_{jj'}^{2})_{20} = C_{jj'}/\sqrt{2j'+1}$.
• Solves the non-LTE integral equations under the impact approximation, Wien limit, and constant magnetic field, neglecting stimulated emission.
• Derives a generalized $√\epsilon$-law for the surface source function components $S^{0}_{0}$ and $S^{2}_{0}$, with $\epsilon^{K}_{Q}$ as the photon destruction probability tensor.
• Uses the kernel integrals $\int \widetilde{K}_{KQ,K'Q'} d\tau'$ to reduce the system to a solvable algebraic form, yielding Eq. (33).
• Applies the resulting formula (38) to test codes under conditions where polarization dominates, making the test highly sensitive to polarization accuracy.

Experimental results

Research questions

• RQ1 How can the $√\epsilon$-law be generalized to include non-thermal collisional excitation and arbitrary magnetic fields in polarized radiative transfer?
• RQ2 What analytical condition emerges for the surface source function components when unphysical collisional rates are introduced to enhance polarization transfer?
• RQ3 Can the resulting test case with dominant alignment ($S^{2}_{0}$) over intensity ($S^{0}_{0}$) serve as a sensitive benchmark for NLTE radiative transfer codes?
• RQ4 How does the inclusion of unphysical source terms in the statistical equilibrium equations affect the validity and utility of the generalized $√\epsilon$-law?

Key findings

• The generalized $√\epsilon$-law is derived in the form $\sqrt{[S^{0}_{0}(0)]^{2} + [S^{2}_{0}(0)]^{2}} = \frac{\epsilon'}{\sqrt{1 - \frac{7}{10}W_{2}(1 - \epsilon')}} B_{\rm P}$, valid for unphysical collisional rates favoring alignment.
• By setting $(C_{jj'}^{0})_{00} = 0$ and $(C_{jj'}^{2})_{20} = C_{jj'}/\sqrt{2j'+1}$, the source function component $S^{2}_{0}(0)$ becomes dominant, unlike in physical cases where $S^{0}_{0}(0)$ dominates.
• The derived formula provides a sensitive test for NLTE radiative transfer codes, particularly for polarization transfer accuracy, due to the enhanced role of $S^{2}_{0}(0)$.
• The method successfully validated a multilevel non-LTE radiative transfer code, demonstrating its utility as a benchmark for polarization modeling.
• The generalized $√\epsilon$-law reduces to known cases (e.g. Mihalas, 1970; Ivanov, 1990) under symmetric physical conditions, confirming consistency.
• The approach enables analytical solutions for highly anisotropic velocity distributions and unphysical collisional rates, expanding the range of testable models.