On the impact of the carbon fusion rate over the properties of superbursts -- Numerical simulations of superbursts with MESA

技術要約：炭素核融合率がスーパーバーストの特性に与える影響について

問題提起
スーパーバーストは、低質量X線連星（LMXB）内の中性子星の地殻で発生するエネルギーの大きな熱核爆発であり、$T \le 10^9$ K の温度における不安定な炭素燃焼によって引き起こされる。従来のマルチゾーン・シミュレーション（例：Keek & Heger 2011, 2012）は、スーパーバーストの再発時間や光度曲線をモデル化することに成功してきたが、これらは古典的な炭素核融合率である Caughlan & Fowler (1988, CF88) に依存していた。近年の核物理学の研究は、この反応率に異議を唱えており、関連する温度域（$T \le 10^9$ K）において、実際の反応率が大幅に異なる可能性を示唆している。具体的には、LUNAコラボレーションによる実験データ（Tumino et al. 2018）は、反応率が $\sim 10^3$ 倍に増強される可能性を示唆している一方で、阻害モデル（Monpribat et al. 2022）は、同様の倍率で減少する可能性を示唆している。これらの大幅な変動が、再発時間、減衰時間、点火深さ、および核合成収量といったスーパーバーストの観測量に与える影響は、現代の恒星進化コードの文脈において未解明のままである。

手法
著者らは、質量降着を受ける中性子星の外層の進化をシミュレートするために、公開されている恒星進化コード MESA v.24.08.1 を用いた。シミュレーションには、縮小された核反応ネットワークである Approx_A56（$A=56$ までの核種をカバー）を使用し、一部の実行では切断効果の検証のために Approx_A56 を用いた。

主なパラメータと設定は以下の通りである：

• 降着率： 2つの標準的なレートをテストした。「中程度の」レート（$\dot{M}_{moderate} \approx 5.88 \times 10^{-9} M_\odot \text{yr}^{-1}$）と「高い」レート（$\dot{M}_{high} \approx 2 \times 10^{-8} M_\odot \text{yr}^{-1}$）。
• 組成： 降着物質は、スーパーバーストのトリガーに必要な文献上の要件と一致するように、20% の $^{12}\text{C}$ と 80% の $^{56}\text{Fe}$ としてモデル化した。初期の外層は純粋な $^{56}\text{Fe}$ である。
• 基底加熱： 以前の KEPLER シミュレーションに一致するように校正された、外層の下部に固定された基底ルミノシティ（$L_b$）を適用し、不安定な炭素点火を確実にした。
• 炭素核融合率： 4つの異なる反応率を比較した：
  1. CF88： 古典的な反応率（標準モデル）。
  2. LUNA： LUNA の測定に基づくもので、$T \sim 2\text{--}3 \times 10^8$ K において CF88 より $\sim 10^3$ 倍高い反応率を予測する。
  3. HIN： 阻害モデルに基づくもので、CF88 より $\sim 10^3$ 倍低い反応率を予測する。
  4. HIN-RES： 低エネルギー共鳴を含む HIN モデルであり、HIN と CF88 の中間的な結果をもたらす。
• 物理的取り扱い： シミュレーションには、放射不透明度（Wolf 処方）、対流のための混合長理論、および固定された表面圧力境界条件が含まれる。一般相対論的補正は省略されている。

主な結果
本研究は、炭素核融合率の変動がスーパーバーストの特性をどのように変化させるかを定量化している：

1. タイムスケールと点火深さ：

   • 増強された反応率 (LUNA)： 再発時間と減衰時間を短縮させ、より浅い柱密度（$y_{ignition}$）での点火を引き起こす。$\dot{M}_{moderate}$ の場合、再発時間は $\sim 241$ 日（CF88）から $\sim 100$ 日へと減少する。
   • 減少した反応率 (HIN)： 再発時間と減衰時間を長くし、より深い柱密度での点火をもたらす。$\dot{M}_{moderate}$ の場合、再発時間は $\sim 364$ 日へと増加する。
   • 中間的な反応率 (HIN-RES)： CF88 の標準モデルとほぼ区別がつかない結果を生む。

2. 温度と光度曲線：

   • 爆発中に到達する最大温度は、反応率に対して敏感である。HIN 反応率は最も高いピーク温度（$\sim 5 \times 10^9$ K）をもたらし、LUNA はより低いピーク（$\sim 3 \times 10^9$ K）をもたらす。
   • ルミノシティの減衰プロファイルもこれに応じて変化する。LUNA シミュレーションはより速い減衰を示し、HIN シミュレーションはより遅い減衰を示す。
   • $\dot{M}_{high}$ において、ピークルミノシティは反応率の選択にはほとんど影響されないが、減衰時間と再発間隔は依然として大きく変化する。

3. 核合成と灰の組成：

   • 灰の分布は、すべてのモデルにおいて特徴的な二峰性の構造（$^{28}\text{Si}$ と $^{56}\text{Fe}$ 付近）を示す。
   • 反応率への感度： $A \le 44$ について、$\alpha$ 核種の存在量は反応率に非常に敏感である。LUNA 反応率（高温／高効率）は、CF88/HIN-RES と比較して有意に多くの $\alpha$ 核種を生成する一方、HIN 反応率はより少ない。
   • 高降着率の影響： $\dot{M}_{high}$ では、$28 \le A \le 50$ の核種の合成における顕著な増強と、$\sim 0.1\%$ の $^{16}\text{O}$ の特異な蓄積が見られる。これは中程度のレートでは目立たない特徴である。
   • 光分解： 不活性な $^{56}\text{Fe}$ を $^{72}\text{Ge}$ に置き換えたシミュレーションにより、鉄族の灰の分布は炭素核融合単独ではなく、主に $^{56}\text{Fe}$ の光分解によって駆動されることが確認された。光分解チャネルを除去すると、$^{40}\text{Ca}$ を超える核合成が停止した。

4. 基底加熱との縮退：

   • 本研究は、炭素核融合率の変化が、外層の基底ルミノシティ（$L_b$）を変化させるのと同等の影響を、スーパーバーストのタイムスケールと点火深さに与えることを示している。例えば、反応率を減少させる（HIN）ことと $L_b$ を増加させることは、同様の再発時間をもたらす。これは、光度曲線のフィッティングのみでは、核融合率と基底加熱の両方を一意に制約できないという縮退が存在することを示唆している。

意義と結論
本論文は、炭素核融合率が $T \le 10^9$ K においてスーパーバーストの物理を支配する極めて重要なパラメータであると結論付けている。著者らは以下の通り主張している：

• 炭素核融合率の不確定性は、質量降着率や基底加熱の変化を必ずしも想定することなく、観測された再発時間や減衰時間の変動を説明できる可能性がある。
• 従来の KEPLER シミュレーション（CF88 使用）から得られた「標準的な」灰の組成は、二峰性の構造については堅牢であるが、$\alpha$ 核種の存在量については特定の反応率に敏感である。
• 炭素核融合率と基底加熱の間の縮退は、単一のスーパーバーストの光度曲線をフィッティングするだけでは、基礎となる核物理学や中性子星の熱状態を明確に決定できないことを意味する。
• 本研究の結果は、スーパーバーストモデリングにおける MESA の使用を正当化するものであり、核物理学の不確定性を探索するための公開コードの選択肢を提供しつつ、従来の KEPLER の結果との一貫性を示している。

著者らは、彼らの知見が、もし LUNA の反応率が正しいならば、高い基底加熱を必要とするモデルの必要性を緩和できる可能性を示唆しているものの、特定の天体（4U 1636-53 など）における質量降着率の不確定性が、どの反応率が物理的に実現されているかについての決定的な結論を妨げていることを控えめに述べている。さらに、彼らは、炭素燃焼を純粋に熱核的なものとして扱うことは、到達する密度と温度において正当化される（ピクノ核的領域とは異なる）ことを明確にしている。