Effect of Thermal Emission in Isotropic Scattering Atmospheres: An Invariant-Embedding Extension of Chandrasekhar's $H(\mu)$-Function

論文技術要約：等方散乱大気における熱放射の影響：Chandrasekharの$H(\mu)$関数の不変埋め込み拡張

1. 背景と問題設定 (Problem)

従来の放射伝達理論、特にChandrasekhar (1960) による半無限大の大気における散乱問題の基礎となる「$H(\mu)$関数」は、外部からの照射（星光など）による拡散反射を記述するものでした。しかし、この古典的な定式化には、大気内部の熱源からの**熱放射（Thermal Emission）**が含まれていません。

ホット・ジュピターや褐色矮星、強照射を受ける系外惑星のような天体では、内部熱や再放射されたエネルギーが放射場に大きな影響を与えます。従来の理論では、散乱と熱放射が同程度に重要な環境を正確に記述できないという課題がありました。

2. 研究手法 (Methodology)

本研究では、Bellmanら(1967)による**不変埋め込み法（Invariant-embedding formalism）**を用い、Chandrasekharの拡散反射フレームワークを拡張しました。

• M関数の導入: 熱放射をソース関数内の埋め込まれた不変な寄与として組み込むことで、従来の$H(\mu)$関数を一般化した新しい角度再分配関数 $M(\mu, U, \tilde{\omega}_0)$ を定義しました。
• パラメータの設定:
  • $\tilde{\omega}_0$: 単一散乱アルベド（散乱の強さを制御）。
  • $U(T) = B(T)/F$: 無次元の熱放射係数（内部熱放射と入射フラックスの比）。
• 数学的導出: $M(\mu)$ に関する非線形積分方程式を導出し、そのモーメント（$A_n$）に関する積分定理（Theorem 1, 2）を数学的に証明しました。
• 数値計算: Gauss-Legendre求積法を用いた安定した反復スキームにより、物理的に意味のあるパラメータ範囲における$M$関数の高精度な数値解を算出しました。

3. 主な貢献 (Key Contributions)

1. 理論的拡張: 散乱と熱放射が同時に作用する半無限大大気における、新しい角度分布関数 $M(\mu)$ の数学的定式化を完了させたこと。
2. 積分定理の確立: $M(\mu)$ のゼロ次および高次モーメントに関する厳密な制約条件（定理）を導出し、数値解の妥当性を検証する手段を提供したこと。
3. 収束条件の特定: 熱放射が物理的に意味を持つ（解が実数となる）ための条件（Corollary 1.1: $U < 0.721$ 付近の制約）を明らかにしたこと。
4. 古典理論との整合性: 熱放射が消失する極限（$U \to 0$）において、本モデルがChandrasekharの古典的な$H(\mu)$関数に正確に一致することを証明したこと。

4. 結果 (Results)

• $M(\mu)$ の振る舞い: $M(\mu)$ は方向余弦 $\mu$、アルベド $\tilde{\omega}_0$、および熱放射係数 $U$ のすべてに対して単調増加する傾向を示します。特に、法線方向（$\mu=1$）において熱放射と散乱の影響が最も顕著になります。
• 放射レジームの分類: 系外惑星K2-137bをモデルケースとして、波長に基づいた3つの放射領域を特定しました。
  • 散乱支配領域 ($\lambda < 0.85\,\mu\text{m}$): $U \ll 1$。古典的な$H(\mu)$関数が適用可能。
  • 同時放射・散乱領域 ($0.85 < \lambda < 2.5\,\mu\text{m}$): $0.1 \lesssim U \lesssim 0.7$。本研究の$M(\mu)$関数が最も有効な領域。JWSTやARIELなどの観測機器の能力に合致。
  • 熱放射支配領域 ($\lambda > 2.5\,\mu\text{m}$): $U > 0.7$。熱放射が支配的となり、拡散反射理論の適用限界となる。

5. 意義 (Significance)

本研究は、系外惑星の分光観測（反射・透過・二次食スペクトル）を解析するための強力な理論的ツールを提供します。特に、ホット・ジュピターのような高温の惑星において、散乱と熱放射の両方を考慮した正確な放射強度計算が可能になるため、惑星の大気組成、エネルギー収支、および熱再分配の理解を深める上で極めて重要な進展といえます。