改进极化辐射传输用于天体物理应用(I):允许非局域Magnus解的形式化方法
Reformulating polarized radiative transfer for astrophysical applications (I). A formalism allowing non-local Magnus solutions
摘要 Abstract
太阳大气通过求解处于非局域热动平衡(NLTE)等离子体的极化辐射传输问题进行诊断。在多维NLTE诊断中的一个关键挑战是高效地积分辐射传输方程(RTE),但目前的方法是局部的,即限于恒定传播矩阵。本文介绍了利用Magnus展开对RTE进行非局域积分的形式化方法。我们首先从洛伦兹/庞加莱群内的旋转角度出发(Stokes形式化),从而引出Magnus展开的应用。结合后者与传播矩阵的高度详细的代数特征,我们推导出一个紧凑的解析演化算符,该算符支持传播矩阵的任意变化,并允许逐步考虑Magnus展开的任意阶次。此外,我们还重新表述了RTE的非齐次部分,再次利用Magnus展开,并引入新的非齐次演化算符概念。这提供了RTE的第一个一致、通用且非局域的形式化解,而且效率高且分离了积分与代数形式解。我们的框架通过解析和计算验证,产生了一组新的数值辐射传输方法及其潜在应用,例如加速NLTE计算。经过微小调整,我们的结果适用于其他在狭义相对论和电磁学中共享洛伦兹/庞加莱代数的普适物理问题。
The solar atmosphere is diagnosed by solving the polarized radiative transfer problem for plasmas in Non-Local Thermodynamic Equilibrium (NLTE). A key challenge in multidimensional NLTE diagnosis is to integrate efficiently the radiative transfer equation (RTE), but current methods are local, i.e. limited to constant propagation matrices. This paper introduces a formalism for non-local integration of the RTE using the Magnus expansion. We begin by framing the problem in terms of rotations within the Lorentz / Poincare group (Stokes formalism), motivating the use of the Magnus expansion. By combining the latter with a highly detailed algebraic characterization of the propagation matrix, we derive a compact analytical evolution operator that supports arbitrary variations of the propagation matrix and allows to increasingly consider any order in the Magnus expansion. Additionally, we also reformulate the inhomogeneous part of the RTE, again using the Magnus expansion, and introducing the new concept of inhomogeneous evolution operator. This provides the first consistent, general, and non-local formal solution to the RTE that is furthermore efficient and separates integration from the algebraic formal solution. Our framework is verified analytically and computationally, leading to a new family of numerical radiative transfer methods and potential applications such as accelerating NLTE calculations. With minor adjustments, our results apply to other universal physical problems sharing the Lorentz / Poincare algebra in special relativity and electromagnetism.