Finite-Slab Reflectance and Transmittance for Henyey-Greenstein Scattering via a First-Passage Transfer Operator

技术摘要：通过第一类通过传输算子计算有限厚度平板的 Henyey–Greenstein 散射反射率与透射率

问题陈述
本文研究了计算具有光学厚度 $\tau$ 的平面平行平板在经历 Henyey–Greenstein (HG) 散射时的反射率 ($R$) 和透射率 ($T$) 的问题。该系统由不对称参数 $a$、单次散射反照率 $a$ 以及入射角 $\theta_0$ 定义。虽然标准方法依赖于蒙特卡洛 (MC) 采样或离散序数/H-函数机制，但本研究寻求一种在公式化上是精确的（超越传输模型本身）、无需随机采样、且能够通过单次算子求值获得整个反照率轴结果的方法。

方法论：第一类通过传输算子 (First-Passage Transfer Operator)
核心方法是一个作用于深度方向状态 $(z, \mu)$ 的第一类通过传输算子，其中 $z$ 是光学深度，$\mu$ 是方向余弦。该方法将完整的 3D 随机游走简化为 $(z, \mu)$ 空间内的闭合马尔可夫链，而无需丢弃方位角相关性或调用扩散近似。

1. 状态简化： 作者证明，逃逸概率仅取决于深度步长序列 $\{z_k\}$ 和方向余弦序列 $\{\mu_k\}$。尽管方位角在演化，但它并不影响与平面平行逃逸相关的横向坐标，因此该过程可以仅作为 $(z, \mu)$ 上的马尔可夫链进行建模。
2. 算子构建：
   • 自由飞行： 光子在固定的 $\mu$ 下进行深度的自由飞行，其路径长度遵循指数分布。
   • 角度重新分布： 在发生碰撞时，方向通过 HG 相函数进行更新。通过对均匀散射方位角进行平均，作者推导出了一个仅取决于入射和出射方向余弦的重新分布核 $\Psi(\mu \to \mu')$。
   • 吸收边界： 算子被限制在具有吸收边界的平板 $(0, \tau)$ 内。到达边界的质量被移除，从而定义了在恰好经过 $n$ 次碰撞后逃逸的光子的阶数解析反射 ($P_R(n)$) 和透射 ($P_T(n)$) 定律。
3. 反照率加权： 对于任何反照率 $a$，总反射率和透射率通过对阶数解析定律进行加权获得：
   $$R(\tau, g, a, \theta_0) = \sum_{n \ge 1} P_R(n) a^n$$
   $$T(\tau, g, a, \theta_0) = e^{-\tau/\mu_0} + \sum_{n \ge 1} P_T(n) a^n$$
   这使得通过单次离散化算子求值即可计算所有 $a \in [0, 1]$ 的 $R$ 和 $T$。
4. 数值实现： 该方法采用均匀深度网格和 Gauss–Legendre 积分进行角度处理。算子被迭代至截止阶数 $n_{max}$，以捕捉保守平板（$a=1$）中的质量。

关键结果与验证

• 准确性： 在广泛的参数空间内（$g \in [0, 0.95]$, $\tau \in [0.5, 16]$, $a \in [0.5, 1]$），该算子重现全 3D 蒙特卡洛结果，绝对误差 $\le 1.6 \times 10^{-3}$。对于 $a=1$，能量守恒（$R+T=1$）满足度达 $< 10^{-4}$。
• 涌现分布： 该方法可以获得涌现角分布 $R(\mu)$ 和 $T(\mu)$，且不增加额外计算成本，其与 MC 数据匹配度在反射方面误差 $< 1\%$，在透射方面误差 $\sim 2\%$（受限于 MC 对尖锐前向峰值的分箱处理）。
• 厚度与各向异性： 结果证实，前向散射 ($g > 0$) 会抑制反射并增强透射。随着厚度增加，$R$ 对于保守介质趋于 1，恢复了扩散性的 $T \sim 1/\tau$ 尾部行为。
• 吸收介质： 在 $\tau \to \infty$ 的极限下，该算子恢复了吸收介质的精确半无限空间反照率，与精确的 H-函数解一致，误差 $\le 6 \times 10^{-4}$。

结构性贡献与意义
论文强调其主要贡献是结构性的，而非仅仅是计算性的。核心结果是有限平板反射率的精确分解：
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
其中：

• $P_\infty(n)$ 是半无限空间第一类返回定律（即在第 $n$ 次碰撞时逃逸半无限介质的概率）。
• $S(n, \tau)$ 是深度生存因子，代表一个阶数为 $n$ 的返回路径保持在 $\tau$ 之内的概率。

这种分解确立了有限平板问题与半无限空间返回统计是“同一个算子的两面”。当 $\tau \to \infty$ 时，平板反射定律 $P_R(n, \tau)$ 项项收敛于半无限空间定律 $P_\infty(n)$，从而恢复了 Sparre–Andersen/Cauchy-kernel 尾部行为。

第二部分：关于生存因子的工作笔记
论文包含关于 $S(n, \tau)$ 解析形式的初步、经验性笔记（第二部分）。

• 定标： 经验数据表明，$S(n, \tau)$ 遵循基于传输平均自由程 $\ell^* = 1/(1-g)$ 和附加深度偏移 $c_g$ 的定标塌缩。
• 布朗运动类比： 生存因子似乎遵循布朗运动偏离最大值的累积分布函数，并经过一个与 $g$ 相关的偏移量和轻微形状修正的修改。
• 状态： 这些发现被明确标记为“工作笔记”和“经验性”内容，尚需等待从基本原理出发的推导。它们并不属于第一部分前向模型中经过验证的结果。

局限性与范围
作者指出，该方法并非旨在超越优化的离散序数代码（如 DISORT）来计算单个辐射亮度。相反，它是互补性的，提供了阶数解析定律以及与其他求解器不同、能直接揭示的与半无限空间统计之间的显式联系。目前的实现提供了方位角平均的角度分布；对于斜入射，通过引入全 BRDF/BTDF 依赖性，可以将其作为一个简单的扩展。论文结论认为，通过第一类通过算子将平板光学与半无限空间返回统计统一起来，才是主要的概念性进展。