Bayesian Inference of Stellar r-Process Abundances

技術要約：ベイズ推論による恒星のr過程存在量

問題提起
r過程（迅速中性子捕獲過程）は、宇宙における重元素の合成を担っている。金属欠乏のr過程強化星や太陽のr過程残留物における重元素（A > 130）の存在量パターンは驚くべき堅牢性を示している一方で、より軽い重元素（Sr, Y, Zr、および第2ピークに至るまで）については、顕著な星間変動が存在する。これらのパターンを説明するためのこれまでの試みでは、一次的プロセス、核分裂、弱いr過程、あるいは複数の天体サイトからの寄与など、多様なメカニズムが示唆されてきた。しかし、観測された恒星の存在量を再現するために必要な具体的な天体物理学的条件（電子分率、エントロピー、膨張タイムスケール）を特定することは依然として困難である。既存の研究は、多くの場合、固定された数の成分や特定の軌跡に依存しており、不確実性を定量化したり、寄与する条件の必要な数を決定したりするための厳密な統計的枠組みを欠いている。

手法
著者らは、r過程の天体物理学的条件を推論するために、マルコフ連鎖モンテカルロ（MCMC）サンプリングを利用したベイズ推論フレームワークを提示している。この手法は以下のコンポーネントを統合している：

• 核合成グリッド： 解析は、120,000個のパラメータ化された核合成計算を含む、サイトに依存しない大規模なサーベイ（KAR2_25）に基づいている。これらの計算は、WinNet反応ネットワークとFRDM2012質量モデルを使用し、初期電子分率（$Y_e$）、初期エントロピー（$s$）、および膨張タイムスケール（$\tau$）によってパラメータ化されている。
• 統計モデル： 観測された恒星の存在量パターンは、$N_{comp}$ 個の核合成成分の線形重ね合わせとしてモデル化される。各成分は、4つのパラメータ（$Y_{e,i}$, $s_i$, $\tau_i$, および重み $w_i$）によって定義される。全存在量は $Y_{sample}(Z) = \sum w_i \cdot Y_{model,i}$ で与えられる。
• 推論の設定： 著者らは、事後確率分布をサンプリングするために emcee Pythonパッケージを採用している。尤度関数は、観測された存在量（ストロンチウムに対して正規化）を比較する簡約カイ二乗指標を介して定義される。
• 事前分布： 2つの異なる事前分布がテストされている：
  1. アグノスティック（無作為）一様事前分布： 物理的な範囲内（$Y_e \in [0.005, 0.5]$、および $s$ と $\tau$ の特定の範囲）のすべてのパラメータに等しい確率を割り当てる。
  2. シミュレーション情報に基づく事前分布： 22種類の異なるシミュレーション（中性子星合体、円盤噴出物、および磁気回転超新星をカバー）から得られた29,000以上の流体解析軌跡を用いたガウスカーネル密度推定（KDE）を用いて構築される。この事前分布は、現在のシミュレーションにおける頻度に基づいて条件を重み付けする。
• 識別可能性： 同一の成分の重複カウントを防ぐため、MCMCサンプラーは、成分が増加する電子分量を持つ必要があるという順序制約（$Y_{e,1} \le \dots \le Y_{e,N_{comp}}$）を課している。

主な結果
本フレームワークは、ほぼ完全なr過程テンプレート星であるHD222925、限定的r過程星（HD128279, HD122563）、および太陽のr過程残留物に適用された。

1. 2成分の支配性： HD222925に対しては、2つの成分があれば観測された存在量パターンを高精度に再現するのに十分である。

   • H成分： 第2ピークから第3ピークまでの元素を生成する、より低い $Y_e$ 条件（$Y_e \approx 0.16-0.18$）（重いr過程）。
   • L成分： 第1ピークから第2ピークまでの元素を生成する、より高い $Y_e$ 条件（$Y_e \approx 0.32-0.33$）（軽いr過程）。
   • 成分の数（$N_{comp} = 3, 4, 5$）を増やしても、追加の成分が主に2成分解と同じパラメータ空間領域をサンプリングするため、観測との中央値の合致度は有意には改善されない。

2. 条件の普遍性： HD222925、HD128279、HD122563、および太陽のr過程残留物のr過程パターンを支配するのは、同じ2組の天体物理学的条件（H成分およびL成分）である。これらの星の違いは、根本的に異なる核合成サイトの存在ではなく、主にこれら2つの成分の相対的な重み付け（$w_2/w_1$）に起因している。

   • HD222925: $w_2/w_1 \approx 1.7$
   • HD128279: $w_2/w_1 \approx 2.7$
   • HD122563: $w_2/w_1 \approx 7.4$
   • 太陽r過程: $w_2/w_1 \approx 4.9$

3. 太陽r過程の不一致： 2成分モデルはHD222925にはよく適合するが、太陽のr過程残留物におけるより軽い中性子捕獲元素（$Z \approx 30-35$）の再現には失敗する。これらはHD222925と比較して著しく強化されている。太陽の第1ピークをフィットするには、より高い $Y_e$（$\approx 0.42$）と長い膨張タイムスケールを持つ第3の成分が必要である。

4. 残留する不一致： モデルと観測の間の残りの不一致は、特定の元素に集中している：

   • 観測の系統誤差： AgおよびCd（$Z=47, 48$）のような元素は、3D効果および非局所熱平衡（NLTE）効果に起因すると思われる不一致を示している。
   • 核物理学： 第3ピーク付近（Re, Os, Ir）およびThに関する不一致は、核質量モデルや核分裂障壁への感度を示唆している。
   • Ga, Ge, As： これらの元素の太陽存在量はHD222925よりも高く、太陽における追加の生成チャネル、あるいは太陽r過程の抽出における問題を示唆している。

5. 事前分布の比較： アグノスティックな一様事前分布を用いて得られた結果と、シミュレーション情報に基づく事前分布を用いて得られた結果はほぼ同一である。これは、HD222925を再現するために必要な条件は、現在の流体解析シミュレーションにおいて既に一般的であり、シミュレートされたパラメータ空間の外にあるエキゾチックな条件を必要としないことを示している。

意義と主張
本論文は、ベイズ推論を介して、大規模な核合成計算と恒星の存在量観測を組み合わせることで、r過程の条件を制約できることを示す力を証明したと主張している。主な意義は以下の通りである：

• 2つの明確な天体物理学的条件（軽い成分と重い成分）の重ね合わせが、金属欠乏星および太陽系で見られる堅牢なr過程パターンを説明できることを確立し、r過程の「普遍性」を支持すること。
• 恣意的な仮定を置くことなく、不確実性を定量化し、必要な成分の数を決定する統計的枠組みを提供すること。
• 残された不一致は、根本的な天体物理学的サイトの理解不足ではなく、観測の系統誤差（NLTE、3D効果）や核物理学の入力（質量、核分裂）に起因する可能性が高いことを強調すること。
• 制約をさらに精緻化するために、より均質な解析が行われたr過程強化星のサンプルと、改良された核・原子データの必要性を動機付けること。

著者らは、特定の天体サイト（例：中性子星合体 vs 超新星）を特定したと主張しているのではなく、むしろ必要な熱力学的条件を特定しており、これらの条件は現在の様々なシミュレーションシナリオで見出されるものであると述べている。