Hydrodynamic simulations of the recurrent nova T Coronae Borealis: Nucleosynthesis predictions

技術要約：再発新星 T Coronae Borealis の流体解析シミュレーション：核合成の予測

問題と背景
再発新星（RN）は、白色矮星（WD）の表面で熱核暴走（TNR）が発生する連星系であり、その再発周期は通常1年から100年と極めて短い。この短い間隔は、非常に高い質量降着率、チャンドラセカール限界に近い白色矮星質量、および高い初期白色矮星光度を必要とする。T Coronae Borealis（T CrB）は、1866年と1946年に噴火した、銀河系内で著名な再発新星である。次回の噴火が間近に迫っていることから、これらの噴火を駆動する条件を理解するために、系の物理パラメータ（具体的には白色矮星質量、初期光度、降着率、および金属量）を制約することが喫緊の課題となっている。さらに、再発新星はIa型超新星のプロジェニター（前駆体）となる可能性があり、その核合成と質量蓄積の特性評価は極めて重要である。

手法
著者らは、T CrB の爆発をシミュレートするために、ラグランジュ型の差分法を用いた流体計算コード SHIVA を用いた。シミュレーションでは、縮退電子ガス、多成分イオンプラズマ、および放射を含む一般状態方程式を使用し、OPAL ロスランド平均不透明度を採用した。核エネルギー生成は、$^1$H から $^{48}$Ti までの120種（630のプロセスに関連）からなる反応ネットワークを用いて計算され、STARLIB データベースから更新された反応率が組み込まれた。

本研究では、様々な条件下で80年の再発周期（$\tau_{rec}$）を再現するように設計された11の新しい流体モデルを提示している：

• 白色矮星質量 ($M_{WD}$): 1.30, 1.35, および 1.38 $M_{\odot}$。
• 初期白色矮星光度 ($L_{WD}$): 0.01, 0.1, および 1.0 $L_{\odot}$。
• 降着物質の金属量 ($Z_{acc}$): 太陽組成 ($1.0 Z_{\odot}$)、低金属量 ($0.1 Z_{\odot}$)、および高金属量 ($10 Z_{\odot}$)。
• 混合: 降着包層と下層の白色矮星（炭素・酸素または酸素・ネオン組成）との間の混合（10%および25%）の影響を評価するため、1.35 および 1.38 $M_{\odot}$ の白色矮星に対して8つの追加モデルを計算した。
• 低質量白色矮星: 低質量の推定値を持つ再発新星を考慮するため、1.20 $M_{\odot}$ の白色矮星について3つのモデルを計算した。

各構成において、8年間の降着後に正確に爆発が起こるように質量降着率（$\dot{M}_{acc}$）を調整した。

主要な結果

• 降着率と物理的条件: 太陽組成の物質を用いて80年の再発周期を実現するには、$M_{WD}$ が 1.30 – 1.38 $M_{\odot}$ かつ $L_{WD}$ が 0.01 – 1.0 $L_{\odot}$ の場合、$\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$ の質量降着率が必要となる。初期白色矮星光度が低い、白色矮星質量が低い、あるいは降着物質の金属量が減少しているモデルでは、8年間の窓の中で臨界点火圧に達するために、より高い降着率が必要となる。
• 爆発ダイナミクス: $L_{WD}$ または $Z_{acc}$ が低いモデルは、より多くの蓄積された包層質量と高い点火圧をもたらす。その結果、これらのモデルは、より高いピーク温度（$T_{max}$）、より大きな放出質量（$\Delta M_{eje}$）、およびより大きな運動エネルギーを伴う、より激しい噴火を生じさせる。例えば、1.38 $M_{\odot}$ の白色矮星で $L_{WD}=0.01 L_{\odot}$ かつ太陽金属量の場合、$T_{max} \approx 2.92 \times 10^8$ K に達するが、1.30 $M_{\odot}$ の白色矮星で $L_{WD}=1.0 L_{\odot}$ の場合は $2.01 \times 10^8$ K にしか達しない。
• 核合成と元素存在量:
  • 同位体比: 放出物（エジェクタ）は CNO サイクル同位体（$^{12,13}$C, $^{14,15}$N）によって支配されている。異なるモデル間で、同位体比（例：$^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si）に顕著な差異が見られる。
  • 白色矮星質量の依存性: 同一の初期光度を持つ異なる白色矮星質量モデル間では、Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, および Sc の存在量に大きな違いが見られる。具体的には、1.38 $M_{\odot}$ モデルは、太陽組成と比較して S, Cl, K, Ca, および Sc が顕著に濃縮されている。一方、1.20 $M_{\odot}$ および 1.30 $M_{\odot}$ モデルは、一般にこれらの元素に対して太陽組成またはやや低めの組成を示す（モデル130Cを除く）。シリコン（Si）は 1.20 および 1.30 $M_{\odot}$ モデルで増強されるが、1.38 $M_{\odot}$ モデルでは枯渇または太陽組成に近い。
  • 金属量の影響: 1.38 $M_{\odot}$ モデルにおいて、金属量の変化は組成を劇的に変える。低金属量（$0.1 Z_{\odot}$）は顕著な $^{40}$Ca 生成をもたらすが、高金属量（$10 Z_{\odot}$）は N, C, Ne, および Si への優位性をシフトさせる。
• 混合シナリオ: 包層と下層の白色矮星基質（CO または ONe）との間の混合を取り入れたモデルは、非混合モデルと比較して、降着フェーズが短縮され、総放出質量が減少する。決定的なことに、混合はガンマ線放出体（$^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al）の生成を大幅に強化する。非混合モデルは無視できるほどのガンマ線放出を予測するが、混合モデルは検出可能なシグネチャーを生じる可能性があり、これにより CO 白色矮星と ONe 白色矮星を区別できる可能性がある。
• ガンマ線放射: 標準的な非混合モデルについては、放射性同位体の収率が低いため、ガンマ線線（511 keV および 1275 keV）および関連する連続成分は、現在の観測機器の検出限界を数桁下回ると著者らは予測している。

意義と主張
本論文は、放出物の組成の違い、特に Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, および Sc といった元素の存在量が、潜在的な診断ツールとして機能すると主張している。これらの変動は、パラメータ空間を制約し、本研究で提案された様々な T CrB モデルを識別することができる。具体的には、著者らは、間近に迫った T CrB の噴火における初期の分光観測によって、白色矮星質量と混合の程度を決定でき、それによってシステムの噴火を支配する物理的条件を明らかにできると主張している。本研究は、ガンマ線の検出は（大幅な混合がない限り）困難であるが、元素存在量パターンは、基礎となる白色矮星の性質と作用している混合プロセスを特定するための堅牢な手法を提供すると結論付けている。