Bayesian Smooth-Fit Extrapolation of the $^{12}\mathrm{C}+{}^{12}\mathrm{C}$ Astrophysical $S$ Factor

技術要約：$^{12}\text{C} + ^{12}\text{C}$ 天体物理学的Sファクターのベイズ平滑適合外挿

問題提起
$^{12}\text{C} + ^{12}\text{C}$ 核融合反応は、恒星の進化の後期段階や、タイプIa超新星およびスーパーバーストのような爆発的な天体現象において極めて重要である。しかし、低エネルギー領域（$E \lesssim 2$ MeV）における天体物理学的Sファクターの信頼できる外挿は、断面積の急速な減少と顕著な実験的不確定性のために、依然として困難な課題となっている。既存のデータセットは不均一であり、異なるエネルギー範囲をカバーし、それぞれ異なる系統的な不確定性と正規化を持っている。さらに、低エネルギー挙動に関する理論的記述には大きな相違があり、低エネルギーでの急激な増大または抑制を示唆するモデルもあれば、そうでないモデルもある。中心となる課題は、これらの異質なデータセットを統計的に一貫した方法で組み合わせ、個々の共鳴構造を解明しようとすることなく、Sファクターの滑らかなグローバル成分を制約することであり、それによって熱核反応率の計算のための堅牢な入力を提供することである。

手法
著者らは、修正天体物理学的因子 $S^*(E) = S(E)\exp(0.46E)$（$E$ はMeV単位）のグローバルな平滑適合を行うためのベイズフレームワークを採用している。解析は、データのダイナミックレンジを縮小し、ガウス型の不確定性に近似するために、対数変数 $y(E) = \log_{10} S^*(E)$ に焦点を当てている。

• モデル: 関数 $y(E)$ は、ピボットエネルギー $E_0 = 1.5$ MeV を中心とした二次多項式によって表される：
  $$y(E; \mathbf{a}) = a_0 + a_1(E - E_0) + a_2(E - E_0)^2$$
  ここで、$\mathbf{a} = (a_0, a_1, a_2)$ は関心対象となるパラメータである。滑らかな関数形式という要件を除き、これらの係数に対して理論的な事前制約（抑制モデルなど）は課されていない。

• データの統合: リカルト関数（尤度関数）には、ノートルダム・グループ、STELLAコラボレーション、および初期の研究からの直接測定値が含まれる。クーロン・核再正規化トロイの木馬法（THM）のデータは、逆キネマティクス・データによって既に制約されている共鳴情報の二重計上を避けるため、外部比較としてのみ含まれている。

• 不確定性の取り扱い: 保守的な不確定性モデルが採用されている。ほとんどのデータセットにおいて、有効な不確定性は、引用された点ごとの統計誤差に、正規化の違いや残留する系統的効果を考慮した固定の分数系統成分（$\delta_{\text{sys}} = 0.25$）を加えたものである。逆キネマティクス・データについては、公表されているエネルギー依存のバンドを直接使用し、最低5%の不確定性を課している。解析は対数空間で行われ、不確定性はガウス分布として扱われる。

• 推論: 平坦な事前分布を用いて、係数の事後分布が決定される。モデルは係数に対して線形であり、不定性が対数空間でガウス型であるため、事後分布は最大事後確率（MAP）推定値と一致する多変量ガウス分布によって近似される。これは最小$\chi^2$解と一致する。係数の共分散行列のモンテカルロ・サンプリングを用いて、不確定性を $S^*(E)$ へと伝播させ、中央値の曲線と中心68%信頼区間を生成する。

• 反応率の計算: 熱核反応率 $N_A\langle\sigma v\rangle$ はハイブリッド入力を用いて構築される。低エネルギー領域（$E \lesssim 2.2$ MeV）では、マックスウェル積分において測定された共鳴構造を明示的に保持するために、直接の逆キネマティクス・データが使用される。高エネルギー領域では、この入力はベイズ平滑事後分布と結合される。

主要な貢献と結果

• 滑らかなグローバル傾向: 解析により、$0.5 \le E \le 3.5$ MeV のエネルギー区間における $S^*(E)$ の滑らかな事後分布が得られた。得られた中央値の曲線は、従来のFowler–Caughlan–Zimmerman (FCZ75) の正規化値を系統的に下回っている。
• 定量的比較: 診断エネルギー $E = 1.5$ MeV において、中央値の事後値は $S^*(1.5 \text{ MeV}) = 1.20 \times 10^{16}$ keV b である。これは、FCZ75 の参照値である $3.0 \times 10^{16}$ keV b よりも大幅に低い。68%信頼区間は $(1.13, 1.28) \times 10^{16}$ keV b である。
• 熱核反応率: 更新された中央値の反応率は、温度範囲 $0.3 \le T_9 \le 1.2$ において、CF88の解析的反応率よりも系統的に低い。比率 $R_{\text{present}}/R_{\text{CF88}}$ は、$T_9 = 0.3$ で約0.33から、$T_9 \approx 0.95–1.0$ 付近の最大値 0.46 まで変化し、$T_9 = 2.0$ では $\approx 0.16$ まで減少する。
• 共鳴の取り扱い: 本論文は、二次適合が個々の共鳴をモデル化していないことを明確にしている。代わりに、共鳴の影響を受けたデータ点は、滑らかな傾向のグローバルな正規化、傾き、および曲率に影響を与える。逆キネマティクス・データで観察される低エネルギー共鳴は、ハイブリッド入力を通じて反応率計算において保持される一方、滑らかな事後分布は高エネルギー側への継続を提供する。

意義と主張
本論文は、特定の抑制モデルの事前分布に依存しない、$^{12}\text{C} + ^{12}\text{C}$ の滑らかなグローバル成分の透明かつ再現可能で統計的に一貫した決定を提供すると主張している。不均一なデータセットを保守的な不確定性モデルと共に組み合わせることで、著者らは更新された事後分布が、従来の参照値（FCZ75）と比較して、より低い滑らかな低エネルギー傾向を支持していることを示している。

著者らは、その結果が低エネルギーにおける $S^*(E)$ の急激な増大を支持するものでも、元の平面波THMによる増大を依拠するものでもないことを強調している。主な物理的出力は、改訂された熱核反応率であり、これは従来の標準的な評価よりも低い。本研究は、ガンモ窓が積分限界を超える温度（$T_9 \gtrsim 1.2–1.3$）において、計算された反応率は完全な恒星反応率ではなく、切り詰められた反応率として見なされるべきであると述べている。このフレームワークはモジュール式であり、将来の直接測定や微視的理論の制約を取り入れて不確定性をさらに減少させることが可能である。