Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON: Comparison to the Galactic Systems

POSYDON連星進化計算コードを用いた本研究は、銀河系の二重中性子星が、特定のコア定義または高い放出効率と低い超新星キックを必要とする2つの異なる共通外層サブチャネル（ケースBおよびケースC）を介して形成されることを明らかにしており、これにより観測された軌道周期の分裂を説明し、これらの系の詳細な進化物理学を制約している。

要約

全体像：宇宙の探偵物語

宇宙を、巨大で混沌としたダンスフロアだと想像してみてください。ほとんどの場合、星たちは一人で踊っています。しかし時折、二つの巨大な星が共に生まれ、手を取り合うことがあります。数百万年をかけて、彼らは激しくドラマチックなダンスを繰り広げますが、その結末は二通りです。永遠に離れ離れになるか、あるいは衝突して合体するかです。

この論文は、**連星中性子星（DNS）**に関する探偵物語です。これらは宇宙のダンスフロアにおける「ヘビー級」の存在です。非常に高密度で死んだ星（中性子星）が二つ、互いの周りを回っています。私たちの銀河系（天の川銀河）には、これらが約25組存在することが分かっています。著者たちは、POSYDONと呼ばれるスーパーコンピュータによるシミュレーションを用いて、これらのペアがどのように形成され、なぜ今のような姿をしているのかを解明しました。

主な発見：同じ目的地への二つの異なる道

この論文における最大の驚きは、連星中性子星を作る方法は一つではないということです。著者たちは、これらが**二つの異なる「サブチャネル（経路）」**を通じて形成されることを発見しました。それは、同じ街へと続く二つの異なる高速道路のようなものです。

その違いは、ドラマが始まる時の二番目の星の「年齢」や「成熟度」によって決まります。

1. 「若い」経路（ヘリウム核チャネル）:

   • シナリオ: まだ十分に成長していない若い星を想像してください。その核はまだヘリウムでできています（まだ高校を卒業していないティーンエイジャーのようなものです）。
   • ドラマ: この若い星が膨張しようとすると、パートナー（最初の中性子星）に近づきすぎます。そして、彼らは**共通外層（コモン・エンベロープ）**と呼ばれる、巨大で乱雑なガスの雲の中に絡め取られます。
   • 結末: このもつれを生き延びるためには、ガス雲を吹き飛ばすための多くのエネルギーが必要です。もし成功すれば、彼らは非常にタイトで速い軌道に入ります。これらが、最終的に衝突して重力波を生み出す（「合体する」）ペアとなります。

2. 「年老いた」経路（炭素・酸素核チャネル）:

   • シナリオ: 今度は、より成熟した年老いた星を想像してください。燃料を使い果たし、核は炭素と酸素でできています（キャリアを終えた退職者のようなものです）。
   • ドラマ: この年老いた星もまた、パートナーと共にガスの雲の中に絡め取られます。
   • 結末: 年老いた星のガス雲は緩やかで押し除けやすいため、システムはこのもつれを非常に容易に生き延びることができます。しかし、彼らははるかに広く、緩やかな軌道に入ります。これらのペアは、宇宙の寿命内に衝突するほど近くはありません。

比喩: これは、混雑した部屋の中で二組のカップルが喧嘩を解こうとしている様子に似ています。

• 若いカップル（ヘリウム核）は、狭く窮屈な場所で喧嘩をしています。群衆を押し退けるには多大な努力（エネルギー）が必要ですが、もし成功すれば、彼らは非常に固く抱き合います。
• 年老いたカップル（炭素・酸素核）は、広々とした部屋で喧嘩をしています。群衆を押し退けるのは簡単ですが、彼らは互いに離れた場所に立つことになります。

「キック」の問題：なぜダンスフロアは滑りやすいのか

星が死ぬとき、それは超新星爆発を起こします。この爆発はしばしば、新しい中性子星に対して、ビリヤードの球を突いた時のように「キック（衝撃）」を与えます。

• 旧来の理論: かつて科学者たちは、これらのキックは巨大でランダムなもの（スレッジハンマーでビリヤードの球を叩くようなもの）だと考えていました。
• 新しい発見: 著者たちは、もしキックがそれほど大きければ、ほとんどのペアはバラバラに飛び散ってしまい、連星中性子星にはなれないことを発見しました。
• 解決策: 私たちが実際に空で見ている姿と一致させるためには、「キック」はもっと穏やか（ソフトなタップのようなもの）でなければなりません。論文は、二番目の星は爆発する前にその質量の大部分を剥ぎ取られているため、非常に穏やかなキックを受けるのだと示唆しています。これにより、ペアが離れずに結びつき続けることが可能になります。

最初の一星の「リサイクル」

二番目の星が爆発する前に、最初の中性子星（古い方）は「メイクアップ（改造）」を受けます。二番目の星がガスを放出する際、最初の一星はそのガスを食べてしまうのです。これは、星の回転を速め、より明るくする「宇宙的なダイエット計画」のようなもので、星を「リサイクル・パルサー」へと変貌させます。

論文は、「若い経路」（ヘリウム核）が最初の一星に対してより優れた食事を提供し、回転をより速めることを示しています。これが、空に見える合体するペア（サブポピュレーション i）が非常に速く回転し、非合体ペア（サブポピュレーション ii）がより遅く回転している理由を説明しています。

変わり者たちについては？

論文では、これら二つの経路のどちらにも当てはまらない、三つの連星中性子星の小さなグループについても言及しています。彼らは説明の難しい奇妙な軌道を持っています。著者たちは、自分たちのモデルではこれらを簡単に説明できないことを認めており、これらは全く別の方法（例えば、ペアとして生まれたのではなく、混雑した星団の中で衝突して形成された可能性など）で形成された可能性を示唆しています。

結果の要約

• 分岐: 連星中性子星は主に二つの経路を経て形成されます。一つは「若い」ドナー（合体するペアにつながる）のため、もう一つは「年老いた」ドナー（合体しないペアにつながる）のためです。
• エネルギー: 「若い経路」を成立させるためには、宇宙がガス雲を放出することにおいて非常に効率的である必要があります。
• キック: これらの星を生み出す爆発は、ペアを繋ぎ止めておくために、暴力的ではなく穏やかでなければなりません。
• 合体: 最終的に衝突して重力波を生み出す（有名なGW170817のような）ペアは、ほぼ間違いなく「若い経路」から来ます。

この論文が述べていないこと

• 次の合体がいつ起こるかを予測するものではありません。
• （上述の「変な星」のように）宇宙のすべての奇妙な星を説明できると主張しているわけではありません。
• 星がどのように「回転が減衰（スピンダウン）」したり、次第に衰退していくかについてのモデルを含んでいないため、現在の現実世界の観測結果と完全に一致させることには限界があります。

要するに、この論文は詳細なコンピュータシミュレーションを用いて、宇宙には連星中性子星を作るための二つの明確なレシピがあり、そのうちの一つのレシピだけが、最終的に衝突して時空の構造を揺るがすペアを生み出すのだということを示しています。

技術概要

技術要約：POSYDONを用いた詳細な連星進化モデルによる二重中性子星の形成

問題提起
二重中性子星（DNS）は、大質量連星進化、超新星メカニズム、および重力波源を理解するための極めて重要な実験場である。銀河系内では2ダース以上のDNSが特定されているが、その形成チャネルは依然として複雑であり、議論の対象となっている。従来のポピュレーション・シンセシス研究の多くは、フィッティング公式や単独星のトラック（軌跡）のライブラリを利用した高速な連星ポピュレーション・シンセシス（BPS）コードに依存してきた。これらの手法は、特に共通外層（CE）フェーズやコア崩壊超新星（SN）における、連星相互作用の自己整合的な取り扱いを欠いていることが多い。さらに、観測される銀河系DNS集団は、軌道周期と離心率において明確な分岐を示しているが、標準的なモデルでは微調整なしにこれを再現することが困難である。本論文は、より厳密で自己整合的な枠組みを用いて、DNS形成に至る詳細な進化段階を制約する必要性に対処するものである。

手法
著者らは、MESAに基づく1次元コードを用いた詳細な単独星および連星の進化モデルを統合した、最先端の連星ポピュレーション・シンセシス・コードであるPOSYDON（POpulation SYnthesis with Detailed binary-evolution simulatiONs）を利用している。高速なBPSコードとは異なり、POSYDONは事前計算された恒星進化グリッドとオンザフライの計算を採用しており、磁気ブレーキ、重力波放射、CE進化、および超新星キックを含む複雑な物理学を扱う。

本研究では、$2 \times 10^6$ 個の連星集団を25の異なるモデルバリエーション（表2に要約）にわたってシミュレートしている。操作された主な変数は以下の通りである：

1. 共通外層（CE）の取り扱い： CE放出効率（$\alpha_{CE}$）のバリエーション、コア・エンベロープ境界の定義（水素質量分率：1%、10%、30%）、およびポストCE質量放出の処理（安定質量移動 vs 風による質量放出 vs 単相処理）。
2. 超新星（SN）モデル： コア崩壊の結果（Sukhbold et al. 2016; Patton & Sukhbold 2020; Fryer et al. 2012）の異なる処方に基づき、鉄コア崩壊（CCSN）と電子捕獲超新星（ECSN）を区別している。
3. 誕生キック（Natal Kick）の処方： 標準的なHobbs-Maxwell分布（$\sigma \approx 265$ km/s）と、「物理に基づいた」モデルとの比較：
   • AsymEj： 非対称な噴出メカニズム（Janka 2017）に由来するキックであり、噴出質量に応じてスケールする。
   • Linear： 噴出質量に線形にスケールするキック（Richards et al. 2023）。
   • Low Maxwellian： 低減された分散モデル（$\sigma \approx 20$ km/s）。

モデルは太陽金属量を仮定し、一定の星形成史を想定している。得られた合成集団は、軌道周期（$P_{orb}$）および離心率（$e$）の分布、ならびに合体率に焦点を当てて、観測された25個の銀河系DNSサンプル（表1）と定性的に比較される。

主要な貢献および結果

• 形成チャネルの分岐： 本研究は、第一誕生中性子星とのロッシュ・ローブ・オーバーフロー（RLO）開始時におけるドナー星の進化状態によって決定される、CE形成チャネルにおける根本的な分岐を特定した。

  • Heコア・チャネル（Case B RLO）： ドナーはヘリウムコアを保持した状態で質量移動を開始する。これらの系は通常、短い軌道周期と高いエンベロープ結合エネルギーを持つ。
  • C/Oコア・チャネル（Case C RLO）： ドナーはRLOの前に、炭素・酸素コアを発達させるまでさらに進化する。これらの系は緩やかに結合したエンベロープを持ち、より広い軌道周期をもたらす。
  • 意義： この分岐は、短周期の合体システム（サブ集団 i）と広軌道の非合体システム（サブ集団 ii）の間にある、銀河系DNS集団の観察された分裂を自然に説明する。

• 共通外層（CE）物理への制約： （合体するDNSを生成する）Heコア・チャネルの形成は、CE効率に強く依存する。モデルは、以下の条件を満たす場合にのみ、このチャネルが生存することを明らかにしている：

  • 寛大なコア定義が使用されている場合（H分率 $\simeq 30\%$）、または
  • 高いCE放出効率が採用されている場合（$\alpha_{CE} \gtrsim 1.2$）。
    これらの条件がない場合、Heコア系はCE内で合体してしまう。一方、C/Oコア・チャネルはより堅牢であり、より低い $\alpha_{CE}$ 値でも生存する。

• 誕生キックの制約： 標準的なHobbs-Maxwell分布（$\sigma \approx 265$ km/s）は、第一超新星の際に大多数の連星を破壊するため、観測される合体率よりも数桁低い値となることから、否定された。観測される銀河系DNS集団は、低速度キック処方（$\lesssim 50$ km/s）によって最もよく再現される。「物理に基づいた」AsymEjおよびLinearモデルは、剥ぎ取られた前駆体に自然に低いキックをもたらし、軌道分布と局所的な合体率密度（$O(10) - O(100)$ Gpc$^{-3}$ yr$^{-1}$）の両方を正常に再現している。

• パルサーのリサイクリングと質量降着： 本研究では、パルサーのリサイクリング（第一誕生NSの回転加速）は、CEフェーズではなく、主にヘリウム星の伴星との安定したCase BB質量移動中に起こることを発見した。Heコア・チャネルを通じて形成された系は、より効率的なリサイクリングが行われ、これは観測される短いスピン周期（サブ集団 (i)）と一致している。

• 超新星の種類と前駆体：

  • 第一超新星は、主に鉄コア崩壊（CCSN）である。
  • 第二超新星は、しばしば高度に剥ぎ取られた（ultra-stripped）前駆体を伴う。合体するDNSの前駆体の約60%は、第二SNの時点でロッシュ・ローブを満たしている。
  • 合体するDNSは、$\sim 0.2 M_\odot$ の中間的なSN前駆体エンベロープ質量から誕生している。これは、合体するDNSのかなりの部分が、Type IbまたはType Ic超新星として観測され得ることを示唆している。

• 合体率と重力波源： モデルは、Heコア・チャネルを通じて形成されたDNSのみがハッブル時間内に合体すると予測している。したがって、GW170817やGW190425のような重力波イベントは、このチャネルに由来すると推論される。物理に基づいたキックを用いたモデルによる予測される局所合体率は、LIGO-Virgo-Kaleraの制約および銀河系DNSサンプルと一致している。

意義および主張
本論文は、詳細な進化モデルがDNSの形成履歴を解明するために不可欠であると主張している。フィッティング公式を超えて自己整合的なグリッドを用いることで、著者らは、観測される銀河系DNSの多様性が、CEフェーズ開始時におけるドナー星の進化状態の自然な帰結であることを示した。

著者らは、現在のPOSYDON内にパルサー進化モデル（選択効果やスピンダウンを考慮するもの）が存在しないため、自らの結論は定性的な比較に基づいていることを控えめに述べている。また、結果が太陽金属量に特有であることも認めている。しかし、本研究は以下の事項を確立している：

1. 形成チャネルの分岐（Heコア vs C/Oコア）は、連星進化の堅牢な特徴であり、観測された軌道周期の分裂を説明する。
2. 合体する集団を形成するには、効率的なCE放出（寛大なコア定義または高い $\alpha_{CE}$ による）が必要である。
3. 低速の誕生キックは、DNS連星の生存のための必要条件であり、高速度キック分布からの最近の理論的転換を支持している。
4. 不安定なCase BB質量移動は、これらの詳細なモデルにおいては合体するDNSを生成せず、超高速合体系の形成に関する従来の仮説に疑問を投げかけている。

本研究は、POSYDONフレームワーク内にパルサー進化モデルが統合された後の、将来的な定量的ベイズ解析のための基礎を提供するものである。