Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics

技術要約：詳細な原子物理学を用いたキロノバ噴出物におけるベータ粒子輸送と熱化

問題提起
中性子星合体による放出物中の$r$-プロセス核の放射性崩壊によって駆動されるキロノバ（KN）放射は、核合成、連星進化、および高密度物質の状態方程式を理解するための重要なプローブである。キロノバの光度曲線を正確にモデル化するには、崩壊生成物（$\beta$粒子、$\alpha$粒子、$\gamma$線、および核分裂破片）の熱化によって提供される加熱率の精密な計算が必要である。$\gamma$線輸送はフル放射輸送で扱われることが多いが、荷電粒子（$\beta$および$\alpha$）の熱化は、歴史的に簡略化された記述に依存してきた。これらの簡略化は、エネルギー堆積の効率、加熱の空間分布、および放出物の電離状態に関して、重大な不確実性を導入する。さらに、荷電粒子の輸送は磁場幾何学に敏感であるが、これは依然として未解明である。本論文は、詳細な原子物理学の欠如と、キロノバの膨張する放出物における$\beta$粒子熱化をより良くモデル化するための相対論的輸送フレームワークの必要性に対処するものである。

手法
著者らは、詳細な原子物理学と相対論的モンテカルロ輸送コードを組み合わせた多角的なアプローチを採用している。

1. 微視的データと断面積:

   • 核生成量とスペクトル: 強、中、弱の$r$-プロセス軌道に対する核生成量は、Prism反応ネットワークを用いて計算される。$\beta$崩壊率とスペクトルは、IAEA Livechart (ENSDF) から取得され、安定に近い核については BetaShape を用いて QRPA+HF モデルにより計算される。
   • 電子－原子相互作用: 電子衝突電離（EII）、電子衝突励起（EIE）、制動放射、およびモット散乱の断面積を決定するために、EPICS2023 データベースの Evaluated Electron Data Library (EEDL) を使用する。
   • 電離補正: 中性原子の断面積は、経験的な Lotz スケーリングと Hartree-Fock 結合エネルギーを用いてイオン化標的に対して調整される。
   • 二次電子: 著者らは、一次$\beta$粒子からの二次および三次電子の生成を追跡し、総電離率を決定するために「電離カスケード」を計算する。
   • その他の相互作用: モラー散乱（電子－電子）は解析的に計算される。シンクロトロン放射およびチェレンコフ損失も評価されるが、衝突損失と比較して無視できることが示されている。

2. グローバル輸送モデル:

   • 放出物の幾何学: 放出物は、べき乗則の密度プロファイル（$\rho \propto r^{-\alpha}$）を持つ、球対称で斉次的に膨張する流れとしてモデル化される。
   • 輸送シナリオ: 2つの限定的な磁場幾何学が検討される：
     • 捕捉限界 (Trapped Limit): 高度に絡み合った、あるいはトロイダル状の磁場に対応し、粒子は特定の速度座標に閉じ込められる。
     • 半径方向磁場限界 (Radial Field Limit): コヒーレントな半径方向磁場に対応し、粒子は磁力線に沿って自由に流れる。著者らは、コヒーレントな半径方向磁場の場合、磁気ミラー効果と幾何学的効果が相殺されるため、1次元輸送方程式が非磁化自由ストリーミングの場合と数学的に等価であることを示している。
   • モンテカルロ・シミュレーション: 完全な相対論的モンテカルロコードが、離散的な$\beta$粒子パケットを追跡する。このコードは、膨張による断熱冷却、衝突エネルギー損失（阻止能）、および散乱（輸送オパシティ）を考慮する。粒子は、熱化（$E < 1$ eV）、脱出（$u > v_{ej}$）、または最大時間（1年）への到達時に終了する。

3. 電離構造:

   • 輸送と電離カスケードが放出物の温度と電離状態に与える影響を示すために、ネオジウム（Nd）を用いた簡略化された速度方程式モデルが適用される。ここでは、異なる元素間の光電離および結合は無視される。

主な貢献と結果

• 原子微物理学: 著者らは、オンラインで利用可能な、種ごとの阻止および輸送断面積、ならびに EII 断面積を提供している。MeV エネルギー領域では EII がエネルギー損失を支配するが、低エネルギー（熱化フェーズ）では、自由熱電子とのクーロン相互作用（モラー散乱）が支配的になることを示している。
• 電離カスケード: 著者らは、総電離率に対する二次電子の寄与を定量化している。カスケードは電離断面積を大幅に増強し、1 MeV においては一次電子と同等の寄与をする。
• 熱化効率:
  • 非局所的効果: 非局所的なエネルギー堆積と粒子の脱出は、特に内側および外側の放出物において、熱化効率を低下させる。これにより、局所的な堆積モデルと比較して、効率は放出物全体でより均一になる。
  • 輸送限界: 半径方向磁場限界では、粒子の脱出が不効率の主な要因となる。捕捉限界では、粒子が放出物中に長く留まることにより、断熱膨張損失が不効率の要因となる。
  • パラメータ依存性: 熱化効率は、放出物質量、速度、$r$-プロセス存在量パターン、および密度プロファイルの勾配に敏感である。電離状態への感度は低く、これは、ある文脈においてはエネルギー堆積を自己整合的な電離計算なしにモデル化できることを示唆している。
• 解析的処方: 著者らは、不効率時間（$t_{ineff}$）と漸近的勾程（$n$）によってパラメータ化された、空間依存的な熱化効率 $f_{thm}(t) = \ln(1+\tau)/\tau$ の解析式を導出した。これらの値の適合パラメータは、放出物質量、速度、および位置の関数として提供される。モンテカルロ結果とこの処方との中央値の残差は、半径方向で $\sim 0.04$ dex、捕捉限界で $\sim 0.02$ dex である。
• 電離構造: Nd に関する例示的な計算は、再結合がより効率的でないために、ネブラ期において「逆転した」電離構造（密度の低い外側領域の方が電離度が高い）が生じることを明らかにしている。半径方向磁場は、粒子の脱出により電離状態を低下させるが、電離カスケードを含めることは、放出物全体で電離度を上昇させる。

意義と主張
本論文は、次世代のキロノバモデルは、光度曲線の予測における不確実性を減らすために、評価された原子データと明示的な荷電粒子輸送を取り入れる必要があると主張している。著者らの主張は以下の通りである：

1. 輸送物理学の重要性: 非局所的な輸送効果を無視すると、内側および外側の放出物における加熱が過大評価され、電離構造が誤ったものになる。
2. 原子データの重要性: 評価された原子データ (EEDL) の使用と電離カスケードの包含は、正確な電離率を得るために必要であり、これらはオパシティや冷却に影響を与える。
3. 解析的有用性: 提供された解析的処方は、半解析的および放射輸送モデルに適した、計算コストの低い空間依存的な熱化効率の近似を提供する。ただし、高精度なシミュレーションには明示的な輸送が推奨される。
4. 限界: 著者らは、自らの結果が理想化された球対称および限定的な磁場幾何学に基づいていることを謙虚に述べている。彼らは、現実的な磁場形態、非球形放出物、および集団プラズマ効果（例：ストリーミング不安定性）が将来の研究領域であることを示唆している。また、彼らの結果は $\beta$ 粒子に焦点を当てており、完全な熱化の全体像を描くためには、$\alpha$ 粒子および核分裂破片に関する将来の研究と統合される必要があることも強調している。