論文要約：21cm 全球信号を用いた第三世代星（Pop III）の性質の制約

1. 背景と問題提起

宇宙の最初の世代の星（第三世代星、Pop III）の性質（典型的な質量、星形成効率など）を直接観測することは、高赤方偏移における暗く微弱な原始銀河の存在により、現在の観測技術では極めて困難です。しかし、Pop III 星は宇宙の再電離（EoR）や金属豊かさの進化に決定的な役割を果たしました。
本研究の目的は、中性水素の 21cm 線による全球信号（Global Signal）の観測が、Pop III 星の**星形成効率（$f_*$）と典型的な恒星質量（$M_s$）**をどの程度制約できるかを理論的に検証することです。特に、前景放射（foreground）の影響を考慮した上で、将来の観測装置（REACH など）による制約の可能性を評価します。

2. 手法とシミュレーションモデル

本研究では、Pop III 星の物理的特性が 21cm 信号に与える影響を半数値シミュレーション（semi-numerical simulation）でモデル化し、フィッシャー行列解析（Fisher analysis）を用いて観測制約を推定しました。

2.1 21cm 信号のモデル化

• 基礎方程式: 21cm 輝度温度の偏移 $\delta T_b$ は、中性水素分率、過密度、スピン温度 $T_S$、CMB 温度 $T_{CMB}$ などの関数として記述されます。
• スピン温度の結合:
  • Wouthuysen-Field (WF) 効果: 最初の星からの赤方偏移した Ly$\alpha$ 光子が $T_S$ をガス運動温度 $T_K$ に結合させます（$T_K < T_{CMB}$ のため吸収信号が生じます）。
  • 加熱: X 線や紫外線（UV）による IGM（銀河間媒質）の加熱が進むと、$T_K$ が上昇し、吸収信号が浅くなり、最終的に放出信号へ移行します。

2.2 物理プロセスの詳細な実装

従来の簡略化されたモデルとは異なり、以下の物理的フィードバックを詳細に組み込みました。

1. 電離光子の脱出率（$f_{esc}$）の質量依存性:
   • 電離光子の脱出率は、ハロー質量（$M_h$）と Pop III 星の質量（$M_s$）に強く依存します。
   • 1 次元放射流体力学（RHD）シミュレーションに基づいたフィッティング式（Tanaka & Hasegawa 2021）を採用し、$f_{esc}(M_h, M_s)$ を計算しました。
   • 一般的に、ハロー質量が大きいほど脱出率は低下し（重力ポテンシャルが深いため）、恒星質量が大きいほど脱出率は上昇します（強い放射圧による電離バブルの急速な膨張のため）。
2. UV による加熱構造:
   • 電離領域の周囲には、UV 光子による加熱を受けた中性ガス層（温かい中性領域）が存在します。
   • 格子サイズより短い平均自由行程を持つため、サブグリッドモデル（電離領域、冷たい中性領域、加熱された中性領域の 3 成分）を導入し、各領域の体積分率を RHD 結果に基づいて算出しました。
3. Lyman-Werner (LW) バックグラウンドによるフィードバック:
   • Pop III 星からの LW 光子は分子水素（H$_2$）を解離し、後の星形成を抑制します。
   • LW 強度（$J_{LW}$）と星形成に必要な最小ハロー質量（$M_{cool}$）の関係を反復計算により自己無撞着に決定し、低質量ハローでの星形成抑制を正確にモデル化しました。

2.3 解析手法：フィッシャー行列

• 将来の観測（REACH など）による制約を評価するため、フィッシャー行列解析を行いました。
• 観測量は全天平均の輝度温度 $T_{sky} = T_{fg} + \delta T_b$ です。
• 誤差源として、前景放射の残差（$\epsilon_0$）と熱雑音（積分時間 $t_{int}$ に依存）を考慮した共分散行列を定義しました。
• 解析対象の赤方偏移範囲は $z=18 \sim 30$ とし、Pop II 星の寄与が支配的になる低赤方偏移領域を除外しています。

3. 主要な結果

3.1 21cm 信号に対するパラメータ依存性

シミュレーション結果から、以下の傾向が確認されました（図 3, 4, 5 参照）：

• 星形成効率（$f_*$）の影響: $f_*$ が大きいほど、より多くの星が形成され、Ly$\alpha$ 光子が増加します。これにより WF 結合が強化され、スピン温度がガス温度に近づき、吸収信号が深く（振幅が大きく）なります。
• 恒星質量（$M_s$）の影響: $M_s$ が大きいほど、電離光子の脱出率が高まり、LW 強度が増加します。その結果、$M_{cool}$ が上昇して低質量ハローでの星形成が抑制され、星形成率密度（SFRD）が低下します。また、強い電離放射による IGM 加熱も促進されます。これらは吸収信号を浅くする方向に働きます。
• パラメータ間の相関: $f_*$ と $M_s$ は、吸収の深さに対して逆の符号で影響を与えるため、両方が同時に増加すると信号の変化が打ち消し合い、正の相関（デジェネラシー）が生じます。

3.2 将来観測による制約能力（フィッシャー解析）

REACH などの将来の観測を想定した解析結果（図 9, 11 参照）：

• 前提条件: 前景残差 $\epsilon_0 = 10^{-4}$、積分時間 $t_{int} = 100 \sim 1000$ 時間、$M_s = 200 M_\odot$、$f_* = 0.01$ のケース。
• 制約精度:
  • 前景残差を $10^{-4}$ に抑え、積分時間を 100 時間とすれば、$f_*$ と $M_s$ をそれぞれ約 24%、23% の精度で制約可能です。
  • 前景残差を $10^{-4}$ に抑え、積分時間を 1000 時間に延長すれば、精度は 11%、12% まで向上します。
  • 積分時間の増加による精度向上は限定的であり、前景放射の除去精度（$\epsilon_0$）が制約の鍵となります。
• 高効率ケース（$f_* = 0.1$）:
  • $f_*$ が大きい場合、UV 加熱が早期に効き始め、吸収の深さに対する $f_*$ と $M_s$ の影響が同じ方向（浅くする方向）に働く赤方偏移領域が現れます。
  • これにより、両パラメータのデジェネラシーが弱まり、より独立した制約が可能になります。

4. 考察と意義

4.1 物理的解釈

• デジェネラシーのメカニズム: 低 $f_*$ の場合、吸収の深さに対する $f_*$ と $M_s$ の影響が逆であるため、両者の相関が強く現れます。しかし、UV 加熱が支配的になる高赤方偏移側では、両者が同様の効果（信号を浅くする）を持つため、異なる赤方偏移領域の観測を組み合わせることでデジェネラシーを解くことが可能になります。
• 初期質量関数（IMF）への拡張: 本研究では単一質量を仮定しましたが、IMF の分布を考慮した平均化を行うことで、IMF の傾きなどのパラメータ制約も可能であると考えられます。特に、X 線加熱（大規模）と UV 加熱（小規模）の両方を考慮することは、パワースペクトル解析において重要です。

4.2 科学的意義

• Pop III 星の直接観測の代替: 直接観測が困難な Pop III 星の性質を、21cm 全球信号という間接的なプローブを通じて定量的に制約できる可能性を示しました。
• 将来実験の指針: REACH や DARE（Dark Ages Radio Explorer）などの次世代実験が、Pop III 星の星形成効率や質量分布を明らかにする上で極めて重要であることを示唆しています。
• JWST 観測との相補性: James Webb 宇宙望遠鏡（JWST）による高赤方偏移銀河の観測結果と組み合わせることで、宇宙初期の星形成史や再電離プロセスの理解が飛躍的に進むことが期待されます。

結論

本研究は、詳細な放射流体力学的フィードバック（脱出率の質量依存性、UV 加熱構造、LW フィードバック）を組み込んだシミュレーションとフィッシャー解析を通じて、将来の 21cm 全球信号観測が Pop III 星の**星形成効率（$f_*$）と典型的な恒星質量（$M_s$）**を 10%〜20% 程度の精度で制約できることを示しました。特に、前景放射の高精度な除去が実現されれば、宇宙の最初の星の性質を解明する有力な手段となり得ます。