Evolution and formation of ultramassive white dwarf stars: The case for a 9Msun progenitor

この論文は、MESA コードを用いて 9 太陽質量の進化星から誕生した 1.313 太陽質量の酸素・ネオン白色矮星の進化を予主系列から冷却曲線まで詳細に追跡し、その組成や熱パルス、冷却遅延などを報告した世界初の研究である。

要約

巨大な星の最期：9 倍の太陽から生まれた「超巨大な白い亡霊」の物語

この論文は、天文学者が**「9 倍の太陽質量を持つ星」が、どのように進化して「超巨大な白色矮星（ホワイト・ドワーフ）」**と呼ばれる死んだ星になるのかを、コンピューターシミュレーションで詳しく描いた研究です。

まるで、星の一生を「映画」のように最初から最後まで再生し、その最期の姿を詳細に分析したような内容です。以下に、専門用語を排して、身近な例え話で解説します。

---

1. 物語の舞台：9 倍の太陽

まず、主人公は**「9 倍の太陽質量」**を持つ星です。
私たちの太陽は、一生を終えると「白色矮星」という小さく重い星になります。しかし、この主人公は太陽の 9 倍も重いため、通常の星とは違う過酷な運命をたどります。

• 通常の星（太陽など）： 燃え尽きて、炭素と酸素の芯（コア）を残して静かに消えます。
• この主人公： 燃え尽きる過程で、さらに高温・高圧になり、**「ネオン（Ne）」や「マグネシウム（Mg）」**といった、もっと重い元素まで作ってしまいます。

2. 過酷な修行：熱いパルスと「鉄の壁」

星の進化の過程で、この星は**「熱パルス（TP-AGB）」**と呼ばれる激しい爆発的な現象を約 140 回も繰り返します。
これは、星の内部で燃料が突然燃え上がり、外側が膨らんで収縮を繰り返す「呼吸」のようなものです。

しかし、ここで大きな壁にぶつかります。

• 鉄の壁（Iron Peak Instability）： 星の表面近くで「鉄」の元素が特定の役割を果たし、計算が破綻してしまう現象が起きます。まるで、コンピューターが「エラー！これ以上進めません！」と叫んで止まってしまうような状態です。
• 研究者の工夫： 過去の研究では、この壁にぶつかって計算が止まってしまい、星がどうなるか最後までわかっていませんでした。しかし、この論文の著者たちは、**「計算が止まる前に、あえて風を強く吹かせて（質量放出を加速させ）、強制的にこの過酷な修行を終了させる」**という大胆な作戦をとりました。
  • これにより、星の進化を「ポスト-AGB（白色矮星になる直前）」までスムーズに導くことができました。

3. 結末：超巨大な「白い亡霊」

最終的に、この星は**「1.313 倍の太陽質量」を持つ白色矮星になりました。
これは、白色矮星が理論的に持てる最大質量（チャンドラセカール限界：約 1.44 倍）に非常に近い、「超巨大な白色矮星（UMWD）」**です。

• 中身（芯）： 炭素や酸素だけでなく、**「酸素（47.7%）」と「ネオン（39.7%）」**が主成分です。まるで、炭酸飲料の瓶の中に、ネオンガスと酸素が混ざり合ったような状態です。
• 外装： 表面には非常に薄いヘリウムの層（太陽の質量の 0.000015 倍程度）が覆っています。
• 特徴： 水素がほとんどない、非常に重く、密度の高い「死んだ星」です。

4. 重要な発見：2 つの驚き

この研究でわかった、2 つの重要なポイントは以下の通りです。

① 「修行（熱パルス）」を省略しても大丈夫？

「鉄の壁」にぶつかる前に、あえて熱パルスを途中で終わらせる（または飛ばす）計算もしました。

• 結果： 最終的な星の重さや中身の成分は、ほとんど変わりませんでした。
• 意味： 計算が難しい「熱パルス」の全過程をシミュレーションしなくても、ある程度正確なモデルを作れることがわかりました。これは、将来の研究者にとって「時短・高効率」な方法として役立ちます。

② 結晶化による「遅延」はわずか

白色矮星は冷えていくと、内部が「結晶化」します（液体から固体になる）。このとき、熱が放出されて冷却が遅れます。

• 結果： この超巨大な星でも、結晶化による冷却の遅れは**「たった 1,600 万年」**でした。
• 意味： 宇宙の時間尺度（数十億年）からすれば、これは「一瞬」です。つまり、この星の年齢を測る時計として、結晶化の影響はあまり気にしなくてよいことがわかりました。

5. 全体のまとめ：なぜこの研究が重要？

この論文は、**「9 倍の太陽質量を持つ星が、どうやって超巨大な白色矮星になるか」**という、これまで完全には解明されていなかった「最後のページ」を埋めました。

• 従来の問題： 計算が途中で止まっていたため、このタイプの星の「中身」や「年齢」が不明だった。
• 今回の解決： 工夫を凝らして計算を完走させ、**「酸素とネオンが混ざった、1.3 倍の太陽質量を持つ超巨大な白色矮星」**の姿を初めて詳細に描き出しました。

比喩で言うと：
これまで、星の進化という「長編小説」の、最も盛り上がるクライマックス（超巨大な星の最期）のページが破れていて、読めませんでした。この論文は、その破れたページを補完し、**「主人公がどんな結晶のような心（コア）を持って、静かに宇宙を漂うことになるのか」**を、初めて完全な形で読み解いたのです。

この研究成果は、将来、宇宙の年齢を測る「時計」として使われる白色矮星の理解を深め、さらには超新星爆発のメカニズム解明にも役立つ重要な一歩となります。

技術概要

以下は、提供された論文「Evolution and formation of ultramassive white dwarf stars: The case for a 9M⊙progenitor（超巨大白色矮星の進化と形成：9太陽質量の progenitor のケース）」に関する詳細な技術的サマリーです。

1. 研究の背景と課題 (Problem)

白色矮星（WD）は銀河系内の恒星の 95% 以上を占める最終進化段階ですが、特に 1.05 太陽質量（$M_\odot$）を超える「超巨大白色矮星（UMWD）」の形成と進化には未解決の課題が多く残されています。

• 最大質量の不明確さ: 白色矮星の progenitor（元となる恒星）の最大質量は、金属量や対流・質量損失の物理的仮定によって 7〜12 $M_\odot$ の間で大きく変動しており、定説がありません。
• 数値計算の困難さ: 超巨大白色矮星は、(S)AGB 段階で炭素燃焼（炭素フラッシュ）を起こし、熱パルス（TP）を数十回から数百回繰り返します。しかし、TP-AGB 段階の終盤、特に「鉄ピーク不安定性（FeI）」や「水素再結合不安定性（HRI）」により、1 次元の静水圧コード（MESA など）は計算が収束しなくなり、ポスト-AGB 段階や白色矮星冷却軌道への移行を計算できないという問題があります。
• 既存モデルの限界: 過去の研究では、TP-AGB 段階を完全にスキップするか、熱パルスを無視して極端な質量損失を仮定するか、あるいは冷却軌道の途中からポリトロープモデルを開始するなどの近似手法が用いられてきました。これにより、TP-AGB とポスト-AGB の両段階を完全な進化計算で扱った UMWD モデルが存在しませんでした。

2. 研究方法 (Methodology)

本研究では、MESA (Modules for Experiments in Stellar Astrophysics) コード（バージョン r24.08.01）を用いて、初期質量 9 $M_\odot$、金属量 $Z=0.02$ の恒星の進化を、主系列前（pre-ZAMS）から白色矮星冷却軌道まで完全計算しました。

• 物理設定:
  • 核反応ネットワーク：mesa45.net（45 同位体、367 反応）。
  • 対流処理：混合長さ理論（MLT）、半対流、オーバーシュート（指数関数的）、熱塩対流。
  • 質量損失：AGB 段階では Bloecker 定式化を採用。
  • 不透明度：OPLIB テーブルと Ferguson テーブル。
• TP-AGB 段階の突破戦略:
  • 鉄ピーク不安定性（FeI）により計算が停止するのを防ぐため、139 回の熱パルスを計算した後、意図的に質量損失率（$\dot{m}$）を劇的に増加させ（$\eta_B=10$、最大 $1 M_\odot$/yr）、TP-AGB 段階からの脱出を強制しました。
  • これにより、HRI や FeI による不安定性を回避しつつ、ポスト-AGB 段階を経て白色矮星冷却軌道に至る計算を完了させました。
• 白色矮星冷却軌道:
  • 相分離（Phase Separation）と結晶化を Skye モジュールを用いて自己無撞着に処理。
  • 元素の沈降（sedimentation）と拡散を考慮。
  • 一般相対性理論（GR）や Urca 過程によるニュートリノ冷却は、今回の質量範囲では無視可能と判断し含めませんでした。

3. 主要な貢献 (Key Contributions)

• 初の完全進化計算: 既知の文献において、TP-AGB 段階とポスト-AGB 段階の両方を計算した最初の超巨大白色矮星系列の提供。
• 高質量 progenitor の詳細なモデル: 9 $M_\odot$ progenitor から生成される 1.313 $M_\odot$ の白色矮星について、冷却時間と詳細な元素存在量プロファイルを初めて公表。
• TP-AGB スキップの妥当性検証: 熱パルスを完全にスキップする手法（早期の極端な質量損失）が、最終質量や白色矮星の化学組成に与える影響が微小であることを示し、UMWD モデル作成の効率的な手法を確立。

4. 主要な結果 (Results)

進化の特性

• 最終質量: 1.313 $M_\odot$ の白色矮星が形成されました。
• コア構造: 炭素燃焼（EAGB 段階）を経て、酸素（O）、ネオン（Ne）、マグネシウム（Mg）を主成分とするコアが形成されました。
  • 質量割合：$^{16}\text{O}$ (47.7%), $^{20}\text{Ne}$ (39.7%), $^{24}\text{Mg}$ (4.2%), $^{23}\text{Na}$ (3.3%), $^{12}\text{C}$ (0.386%)。
  • 炭素の総質量は $5 \times 10^{-3} M_\odot$ です。
• 外層: ヘリウム層は $1.5 \times 10^{-5} M_\odot$、水素はほぼ完全に燃焼・消失しており、水素欠乏（H-deficient）の白色矮星となりました。
• 熱パルスの影響: 139 回の熱パルスを計算しましたが、最終質量への寄与はわずかでした。熱パルスをスキップした場合でも、最終質量と組成への影響は最小限（質量差は約 0.05 $M_\odot$ 以内）であることが確認されました。

冷却軌道と相分離

• 冷却時間: 有効温度 10,000 K 到達までの冷却時間は約 1.667 Gyr、恒星全体の年齢は約 1.698 Gyr でした。
• 相分離の影響: O/Ne 混合における相分離を考慮した場合、冷却遅延はわずか 16 Myr でした。これは、この質量帯では相分離が冷却時間に与える影響が小さいことを示唆しています。
• 結晶化: 結晶化は約 48,530 K で開始され、50% 結晶化までに約 160 Myr、90% 結晶化までにさらに約 270 Myr 要しました。

既存研究との比較

• Lauffer et al. (2018) や Camisassa et al. (2019) のモデルと比較し、冷却時間や結晶化開始温度において概ね良好な一致（冷却年齢で 20% 以内の差異など）を確認しましたが、組成の違い（特に Ne/O 比）が冷却曲線に影響を与えることが示されました。
• 観測対象である磁気白色矮星 ZTF J1901+1458 に対して、本研究のモデルは半径 2406 km、冷却年齢 506 Myr を予測し、既存の観測値と整合性があることを示しました。

5. 意義と結論 (Significance)

本研究は、チャンドラセカール限界に近い質量を持つ白色矮星の形成メカニズムを、TP-AGB 段階の不安定性を回避する新しいアプローチで解明した点で画期的です。

• 理論的基盤の強化: 9 $M_\odot$ 恒星が超巨大白色矮星（ONeMg コア）として進化しうることを示し、IFMR（初期質量 - 最終質量関係）の上限に関する議論に重要なデータを提供しました。
• 計算手法の確立: 鉄ピーク不安定性に直面する高質量 AGB 星の計算において、意図的な質量損失増加による「脱出」手法が有効であることを実証しました。
• 将来の応用: 得られた詳細な化学組成プロファイルと冷却曲線は、白色矮星の観測データ（質量、年齢、組成）の解釈や、Ia 型超新星の progenitor 候補としての評価、および銀河系化学進化モデルの精度向上に不可欠なリソースとなります。

要約すると、この論文は数値計算の困難さを克服し、超巨大白色矮星の完全な進化シナリオを初めて提示することで、高質量白色矮星の物理学における重要なブレイクスルーを達成したものです。