Simulating the Convective Urca Process with Multiple Urca Pairs in a Simmering White Dwarf

技術サマリー：燃焼中の白色矮星における複数のウルカ対を用いた対流ウルカ過程のシミュレーション

問題提起
Ia 型超新星（SNe Ia）は白色矮星の熱核爆発であるが、その前駆星システムと爆発メカニズムは未だ完全には解明されていない。チャンドラセカール質量に近い質量モデルでは、爆発に先行して「燃焼相（simmering phase）」が存在し、これは対流を駆動する緩やかな中心部での炭素燃焼を特徴とする。この相の重要なかつ未解明な構成要素が、弱い核反応（ベータ崩壊と電子捕獲）と対流混合を結びつける「対流ウルカ過程」である。従来の理論研究および 2 次元シミュレーションでは、この過程が対流を著しく妨げ、対流領域のサイズを制限し、白色矮星の進化を変化させる可能性が示唆されていた。しかし、最近の 3 次元シミュレーション（Boyd et al. 2025a; Poca-Amorós et al. 2026）は、その効果がより制限的ではない可能性を示したが、これらの研究は中性子過剰に関するエネルギー生成率や組成変化を正確にモデル化できなかった簡略化された核反応ネットワークを利用していた。

手法
著者らは、燃焼中の白色矮星の対流領域の 3 次元シミュレーションを実行するために、低マッハ数流体力学コード MAESTROeX を用いた。主要な手法上の進歩は以下の通りである。

• 包括的な反応ネットワーク: 本研究では、以前の 3 次元研究と比較して、大幅に拡張された核反応ネットワーク（33 の反応率、21 の原子核）を実装した。このネットワークは炭素燃焼を正確に追跡し、3 つのウルカ対、すなわち A=21（$^{21}$Ne–$^{21}$F）、A=23（$^{23}$Na–$^{23}$Ne）、および A=25（$^{25}$Mg–$^{25}$Na）を含んでいる。
• シミュレーション設計: 対流ウルカ過程の効果を分離するために、2 つのシミュレーションが実行された。
  • フルネットワーク（FN）: ベータ崩壊を含むすべての反応率を含む。
  • ベータなし（NB）: FN と同一であるが、ベータ崩壊の反応率を除去しており、これにより対流ウルカサイクルが実質的に破断されている。
• 初期条件: モデルは、中心密度 $4.5 \times 10^9$ g cm$^{-3}$、温度 $5.5 \times 10^8$ K を用いた。初期組成には、小さな $^{12}$C-$^{16}$O 勾配と、迅速な混合限界近似に基づいたウルカ対の分布が含まれていた。数値的な励起を最小化するために、断熱温度勾配と等温温度勾配の間に「滑らかな」遷移が導入された。
• 数値設定: シミュレーションは、対流領域で有効分解能 5 km を達成する 3 層のメッシュ細化を用いた「全星」幾何学で行われた。シミュレーションは定常状態に達するまで 2500 秒以上（対流ターンオーバーの数十回分）進化させた。

主要な結果

1. 対流領域のサイズ: いくつかの以前の 2 次元および 1 次元の予測とは異なり、対流ウルカ過程は対流領域の全体的なサイズを制限しなかった。FN および NB の両方のシミュレーションにおいて、混合は初期の等エントロピー領域の端（約 700 km）まで広がった。
2. 混合の効率: 領域のサイズは同様のままだったが、FN シミュレーションでは対流境界付近の混合効率が低下した。550 km から 700 km の間の明確な領域は、組成勾配によって安定化された（ルドルフ安定だがシュワルツシルト不安定な）超断熱温度勾配を特徴とする非効率な混合を示した。
3. A=23 対の優位性: A=23 ウルカ対が対流ウルカ過程の最も重要な駆動因子として特定された。これは、中心部での炭素燃焼の直接生成物であるその高い存在量と、ウルカ殻に対するその位置（電子捕獲とベータ崩壊の両方の実質的な領域を可能にする）によるものである。A=21 および A=25 の対は、存在量が低く、殻の位置が不利であるため、無視できる影響しか及ぼさなかった。
4. 以前の研究との比較（PA2026）: 結果は、著しく小さな対流領域と潜在的な「分裂」構造を見出した Poca-Amorós et al. (2026) とは異なる。著者らは、この不一致を PA2026 モデルにおけるA=23 対のより高い存在量とより大きな初期等エントロピー領域に起因すると帰属させた。本研究は、ウルカ過程の制限的な性質がこれらの初期前駆星条件に敏感であることを示唆している。
5. エネルギー収支: 外部領域におけるウルカ反応からの正味のエネルギー生成は modest（中心部での炭素燃焼よりも桁違いに小さい）であったが、この過程は依然として、組成勾配と浮力の変化を通じて対流構造を変化させる上で重要な役割を果たした。

意義と主張
本論文は、炭素燃焼と中性子化の物理をよりよく反映する包括的な反応ネットワークを利用することで、対流ウルカ過程に関するより正確な理解を提供すると主張している。主要な発見は、対流ウルカ過程が境界付近の対流混合の効率を低下させるものの、モデル化された特定の初期条件下では、対流領域の空間的範囲を必ずしも制限しないということである。

著者らは、結果が初期前駆星モデル、特に A=23 ウルカ対の存在量と等エントロピー領域のサイズに極めて敏感であることを強調している。彼らは、A=23 対がこれらのシミュレーションにおける支配的要因であると結論づけている。この研究は対流境界の複雑さを浮き彫りにし、将来の 1 次元恒星進化モデルが、Ia 型超新星の最終的な構造と爆発をよりよく予測するために、これらの非局所平衡効果と半対流領域の可能性を考慮しなければならないことを示唆している。本研究は燃焼の全時間スケールの解決を主張するものではないが、この重要な相における流体力学的および核的な相互作用の洗練されたスナップショットを提供する。