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

技術要約：第一通過転送作用素によるヘニェイ・グリーンシュタイン散乱の有限スラブ反射率および透過率

問題提起
本論文は、ヘニェイ・グリーンシュタイン（HG）散乱が生じる光学厚さ $\tau$ の面平行スラブにおける反射率 ($R$) および透過率 ($T$) の計算について扱う。系は非対称パラメータ $g$、単一散乱アルベド $a$、および入射角 $\theta_0$ によって定義される。標準的な手法はモンテカルロ（MC）サンプリングや離散オルディネート法、あるいはH関数に依存するが、本研究が目指すのは、定式化において（輸送モデルを超えて）厳密であり、ランダムサンプリングを必要とせず、単一の作用素評価から全アルベド軸にわたる結果を提供できる手法である。

手法：第一通過転送作用素
核心となる手法は、深さ方向の状態 $(z, \mu)$（ここで $z$ は光学深さ、$\mu$ は方向余弦）に作用する**第一通過転送作用素（first-passage transfer operator）**である。このアプローチは、方位角の相関を破棄したり拡散近似を導入したりすることなく、完全な3次元ランダムウォークを $(z, \mu)$ における閉じたマルコフ連鎖へと還元する。

1. 状態の簡約化： 著者らは、脱出確率が深さステップ $\{z_k\}$ と方向余弦 $\{\mu_k\}$ の配列のみに依存することを証明している。方位角は進化するものの、面平行な脱出に関連する横方向の座標には影響を与えないため、この過程は $(z, \mu)$ のみのマルコフ連鎖としてモデル化できる。
2. 作用素の構築：
   • 自由飛行（Free Flight）： 光子は固定された $\mu$ で深さ方向への自由飛行を行い、その経路長は指数分布に従う。
   • 角度再分配： 衝突時、方向はHG位相関数に従って更新される。一様な散乱方位角について平均化することにより、著者らは入射方向余弦と放出方向余弦のみに依存する再分配カーネル $\Psi(\mu \to \mu')$ を導出している。
   • 吸収境界： 作用素はスラブ $(0, \tau)$ 内に限定され、吸収境界を持つ。境界に到達した質量は除去され、これにより、ちょうど $n$ 回の衝突を経て脱出した光子の次数分解された反射則 $P_R(n)$ および透過則 $P_T(n)$ が定義される。
   • アルベドの重み付け： 任意のアルベド $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$ の計算が可能となる。
3. 数値実装： 本手法は、一様な深さグリッドとガウス・レジェンドル求積法を用いている。作用素は、保存的なスラブ（$a=1$）において質量を十分に捉えることができるカットオフ次数 $n_{max}$ まで反復される。

主要な結果と検証

• 精度： 本手法は、広範なパラメータ空間（$g \in [0, 0.95], \tau \in [0.5, 16], a \in [0.5, 1]$）において、フル3Dモンテカルロ結果を絶対誤差 $\le 1.6 \times 10^{-3}$ 以内で再現する。エネルギー保存（$R+T=1$）は、$a=1$ のとき $10^{-4}$ 未満の精度で満たされる。
• 出現分布： 本手法は、追加の計算コストなしに、出現角分布 $R(\mu)$ および $T(\mu)$ を算出する。これはMCデータと、反射については誤差 $<1\%$、透過については（鋭い前方ピークのMCビニングに起因する制限はあるが）$\sim 2\%$ の誤差で一致する。
• 厚さと異方性： 前方散乱（$g > 0$）が反射を抑制し、透過を促進することを確認している。厚さが増加するにつれ、$a=1$ の保存的媒体において $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 の裾の挙動を回復する。

第II部：生存因子に関する作業ノート
本論文には、$S(n, \tau)$ の解析的形態に関する予備的な経験的ノート（第II部）が含まれている。

• スケーリング： 経験的なデータは、$S(n, \tau)$ が輸送平均自由行程 $\ell^* = 1/(1-g)$ と加法的な深さオフセット $c_g$ に基づくスケーリング崩壊（scaling collapse）に従うことを示唆している。
• ブラウン運動との類似： 生存因子は、$g$ に依存するオフセットと緩やかな形状補正を受けた、**ブラウン運動の極大値（Brownian excursion maximum）**の累積分布関数に従っているように見える。
• 現状： これらの知見は、明示的に「作業ノート」および「経験的」であるとラベル付けられており、第一原理からの導出を待つ段階にある。これらは、第I部のフォワードモデルとして検証された結果とは主張されていない。

限界と範囲
著者らは、本手法が単一の放射輝度を計算するための最適化された離散オルディネート・コード（DISORTなど）を凌駕することを目的としていないことを述べている。むしろ、本手法は、他のソルバーが直接的に露出させることのない、次数分解された法則と半無限統計との明示的な接続を提供する補完的なものである。現在の実装は方位角平均された角度分布を提供するが、斜め入射に対する完全なBRDF/BTDF依存性への拡張は、単純な拡張として示唆されている。論文は、スラブ光学と半無限帰還統計の第一通過作用素による統一が、主要な概念的進歩であると結論づけている。