Black Hole-Neutron Star Binaries near Neutron Star Disruption Limit in the Mass Regime of Event GW230529

技術要約：GW230529の質量領域における中性子星潮汐破壊限界近傍のブラックホール–中性子星連星

問題提起
LIGO–Virgo–KAGRA (LVK) コラボレーションによる重力波イベント GW230529 の検出は、電磁波（EM）対応天体を生じる可能性が高い、最初のブラックホール–中性子星（BHNS）合体である可能性を示した。このイベントは、主成分となるブラックホール（BH）が 2.5–4.5 $M_\odot$、副成分となるコンパクト天体が 1.2–2.0 $M_\odot$ であり、中性子星（NS）がブラックホールへ直接落下するのではなく、潮汐破壊される可能性がある質量領域に位置している。潮汐破壊は、降着円盤の形成と中性子に富む物質の放出を必要とするキロノバやコンパクト連星ガンマ線バースト（cbGRB）などの観測可能な電磁波信号を生成するための前提条件である。

しかし、特定の BHNS 合体が潮汐破壊をもたらすかどうかを予測することは、合体後の残骸円盤の質量を推定する FNH18 モデル（Foucart et al. 2018）のような解析的なフィッティング公式に大きく依存している。これらのモデルは、主に質量比 $q \gtrsim 3$ かつ、潮汐破壊に必要な臨界閾値から遠いスピンを持つシミュレーションに基づいて校正されている。その結果、GW230529 と一致する質量およびスピンパラメータを持つ、潮汐破壊限界付近の系に関する高精度な数値相対論（NR）データの不足が生じている。このギャップは、EM信号の発生率を正確に予測し、EM観測から中性子星の状態方程式（EOS）を制約する能力を制限している。

手法
著者らは、Spectral Einstein Code (SpEC) を用いて、BHNS 合体の高解像度数値シミュレーションを実施した。研究では、GW230529 イベントと一致するパラメータである質量比 $q=3$ および中性子星質量 $M_{NS} = 1.35 M_\odot$ に焦点を当てている。中性子星のコンパクト性のレベルをテストするために、2 つの異なる状態方程式を採用している：

1. SFHo: 半径 $R_{NS} \approx 11.9$ km の軟らかい EOS。
2. DD2: 半径 $R_{NS} \approx 13.2$ km の硬い EOS。

シミュレーションでは、ブラックホールスピンの大きさ（$\chi_{BH}$）を変化させ、FNH18 モデルが予測する臨界スピン（$\chi_{crit}$）の直近の領域を調査している。SFHo の場合、$\chi_{crit} \approx 0.037$ であり、DD2 の場合は $\chi_{crit} \approx -0.23$ である。研究には、これらの臨界閾値をわずかに下回るスピンと、上回るスピンの両方を含むシミュレーションが含まれている。

主な手法の特徴は以下の通りである：

• モンテカルロ・ニュートリノ輸送: これは、BHNS 合体シミュレーションにおけるモンテカルロ・ニュートリノ放射輸送の初の適用である。リークモデルとは異なり、この手法は吸収、放出、および散乱を正確にモデル化するために、位置と運動量空間におけるニュートリノの軌跡を追跡する。これは、合体後の残骸および流出物の組成（電子分率、$Y_e$）を決定するために極めて重要であり、それは r プロセス核合成の可能性を左右する。
• 解像度研究: 数値的な不確実性を定量化するために、SFHo のケースについて 3 つの解像度で収束テストが行われた。
• 初期条件: 初期データは Spells ソルバーを用いて生成され、軌道離心率は $e \lesssim 0.009$ に低減された。

主な結果

1. 円盤質量 vs. 解析的予測:

   • SFHo EOS: 残骸の総質量に関する数値結果は、一般に FNH18 公式の予想されるモデリング誤差範囲内に収まっている。しかし、シミュレーションは一貫して、モデルが予測するよりもわずかに高い残骸質量を示している。
   • DD2 EOS: 重大な不一致が観察された。潮汐破壊限界付近のスピンにおいて、FNH18 モデルは、モデルに記載された誤差範囲を超える係数で残骸円盤質量を過小評価している。例えば、$\chi_{BH} = -0.030$ のとき、モデルは $\sim 0.012 M_\odot$ の残骸質量を予測するが、シミュレーションは $\sim 0.049 M_\odot$ を算出する。
   • 臨界スピン閾値: 質量の不一致にもかかわらず、潮汐破壊を誘発するために必要な「臨界スピン」の値は、解析モデルとシミュレーションの間で一致していることが判明した。主な違いは、閾値付近における残骸質量の曲線の傾き（$dM_{rem}/d\chi_{BH}$）にあり、シミュレーションはモデルよりも急激な上昇を示している。

2. 低質量円盤の出現:

   • シミュレーションにより、予測された潮汐破壊限界をわずかに下回るスピン（例：$\chi_{BH} = -0.0130$ の SFHo ケース）であっても、ゼロではない降着円盤が形成されることが明らかになった。これらの円盤は低質量（$\sim 10^{-3} M_\odot$）であるが、先行研究（Gottlieb et al. など）が示唆するように、短時間の cbGRB を駆動するのに十分である可能性がある。

3. 質量放出と組成:

   • 動力学的放出物: 合成中に動力学的に放出される質量の量は非常に少なく、通常 $\lesssim 0.004 M_\odot$ である。最も低いスピン構成では、非束縛質量は数値解像度の限界に近く、不確実性は $\sim 100\%$ に達する可能性がある。
   • 円盤風: 本研究は、これらの系において、後期の円盤風（粘性および磁気駆動）を通じて放出される総質量が、初期の動力学的放出量を上回る可能性が高いことを示唆している。
   • 電子分率 ($Y_e$): 動力学的放出物は高度に中性子に富んでいる（$Y_e < 0.1$）。しかし、流出物および円盤物質は、弱い相互作用（ニュートリノ加熱）により、より広い $Y_e$ の範囲（最大 $\sim 0.3$）を示す。DD2 EOS は SFHo よりも中性子に富む放出物を生成する。

4. 観測的意義:

   • これらの特定の GW230529 のような合成からのキローバ信号は、動力学的放出物による赤外線放射よりも、円盤風の光学的放出によって支配されるため、暗いと予測される。
   • r プロセス核合成に利用可能な総放出質量は $\lesssim 0.02 M_\odot$ である。

意義と主張
本論文は、FNH18 解析モデルは潮汐破壊の存在（臨界スピン）を予測する上では概して信頼できるものの、特に硬い状態方程式（DD2）において、結果として生じる残骸円盤の質量を系統的に過小評価していると主張している。

著者らは、自らの結果が、潮汐破壊閾値に近い BHNS 合への現在の校正の限界を浮き彫りにしていることを強調している。具体的には：

• 予測された臨界スピン以下において、残骸質量を「ゼロ」とするモデルの仮定は厳密には正確ではない。低質量の円盤も形成され得る。
• これらの低質量円盤は、明るいキロノバは生成しないとしても、短時間のガンマ線バーストの潜在的なプロジェニターとして物理的に重要である。
• これらの知見は、低質量比かつ臨界スピン近傍の領域におけるシミュレーションがさらに蓄積されるにつれ、フィッティング公式の緩やかな再校正が必要になる可能性があることを示唆している。

本研究は、LVK による GW230529 が潮汐破壊型であるかどうかの確率評価を決定的に変えるものではないとし、観察された差異が LVK の解析結果を質的に変える可能性は低いと述べている。しかし、この特定の、観測的に極めて重要なパラメータ空間における BHNS 連星の挙動を理解するための、より厳密な数値的基礎を提供している。