Numerical simulations of black hole-neutron star mergers with equal and near-equal mass ratios

技術要約：等質量および近等質量比におけるブラックホール-中性子星合体の数値シミュレーション

問題と動機
重力波イベントGW230529の検出は、質量比 $q \approx 1/2.6$ のブラックホール-中性子星（BHNS）合体と解釈されており、これまで予想されていたよりも対称的なBHNS系の存在を示唆している。標準的な恒星進化モデルは $q \lesssim 1/3.5$ の質量比を予測するが、高密度星団における動力学的形成チャネルや第二世代の合体は、$q \sim 1/2$ や $q \approx 1$ にさえ達する系を生み出す可能性がある。これらの対称的な系は、理論的に大規模な降着円盤を生成し、キロノバのような明るい電磁波過渡現象を駆動することが期待される。しかし、等質量および近等質量比のパラメータ空間は、数値相対論（NR）シミュレーションによって十分に探索されていない。このデータの不足は、重力波波形モデルや残留物の特性に関する解析的フィッティング公式の精度を制限し、実際のイベントの誤分類や、電磁波対応天体の特定機会の喪失を招く可能性がある。

手法
著者らは、BAMコードを用いて、非回転ブラックホール-中性子星合体の12個のNRシミュレーションを行う。シミュレーションは、質量比 $q \in \{1, 1/2, 1/3\}$ および中性子星（NS）質量 $m_{NS} \in \{0.8, 1.2, 1.4, 2.2\} M_\odot$ をカバーしている。2つの状態方程式（EOS）として、軟らかいSLyと硬いDD2が用いられている。シミュレーションでは、EOSに区分的多項式（piecewise-polytropic）表現を用い、時空に対しては共形Z4定式化、流体力学に対してはバレンシア定式化を用いて、約9軌道の進化を行う。初期データはFUKAコードを用いて構築され、離心率は $10^{-3}$ 未満に低減されている。

本研究は以下の項目を含む：

1. 重力波波形解析： 複数の半径で波形を抽出し、二次外挿を行い、重心ドリフトを補正する。これらをIMRPhenomNSBHおよびSEOBNv4_ROM_NRTidalv2_NSBH波形モデルと比較する。
2. 残留物特性の検証： 動的放出物（dynamical ejecta）の質量と降着円盤の質量、および合体後10 msにおける残留ブラックホール（BH）の質量とスピンを算出する。これらを文献[28, 31, 40, 41]の解析的フィッティング公式と比較する。
3. 円盤力学： 有意な円盤質量を持つ構成について、形態、回転則、および残留円盤の全域的な密度振動を分析する。
4. 降着モデリング： 降着挙動を分類し、円盤振動による変調を特定する。
5. キロノバ・モデリング： 動的放出物のプロファイルと、追加された球状風成分（円盤質量の30%）に基づき、POSSIS放射伝達コードを用いてキロノバの光度曲線を描く。

主要な結果

• 重力波波形： 現在のBHNS波形モデル（IMRPhenomNSBHおよびSEOBNv4）は、対称的な質量比の領域において、NR波形との一致が不十分である。ほとんどのケースで、合体時に約1ラジアンの位相差が生じている。これらのモデルは、NRデータと比較して、合体時刻を早く予測し、リンダウン段階を長く予測する傾向がある。特に、亜太陽質量NSのケース（$m_{NS}=0.8 M_\odot$）では、潮汐変形能が既存モデルの有効なパラメータ空間を超えており、外挿した場合、モデルは顕著な位相および振幅の不一致を示す。
• 残留物特性：
  • 放出物および円盤質量： 動的放出物質量および円盤質量に関する解析的フィッティング公式は、シミュレーションデータと良好に一致している。
  • 残留BH質量： 解析的公式は、残留BH質量を極めて正確に予測している（誤差は一般に1%未満）。
  • 残留BHスピン： 更新されたフィッティング公式は、残留スピンの予測においては精度が低く、特に潮汐変形能が大きい系（DD2 EOS）で顕著である。亜太陽NSのケースでは、スピン予測が約50%逸脱している。著者らは、これが潮汐変形能がより低く、より非対称な質量比を持つ系への過学習に起因している可能性を示唆している。
• 円盤構造と力学：
  • 形態： 等質量（$q=1$）の円盤は、合体直後（$\sim 3$ ms）にほぼ円環化するが、$q=1/2$ の円盤は螺旋密度波や落下物質による摂動をより長く受け、$\sim 8$ ms経過して初めて軸対称へと落ち着く。
  • 回転則： 比角運動量は冪乗則 $j \propto \Omega^{-n}$ に従う。指数 $n$ は半径方向に変化し、ISCO内部では一定角運動量（$n=0$）に従うが、外側円盤では劣ケプラー回転へと遷移する。これは、長期間の軸対称シミュレーションにおける固定された回転則の一般的な仮定と矛盾する。
  • 振動： 円盤質量が $10^{-2} M_\odot$ を超える場合、全域的な捕捉された $g$ モード密度振動が特定される。その周波数は、残留BHの最大順行エピサイクリック周波数と一致する。
  • 降着： 3つの降着タイプが区別される。タイプB（$q=1/2$）は、全域的な $g$ モード振動によって直接変調される降着率の変動を示す。タイプA（$q=1/2$）およびタイプC（直接落下）は異なる挙動を示す。一部の系における $g$ モード周波数は、ガンマ線バーストで観測される高周波準周期振動（QPO）と互換性がある。
• キロノバ放射： モデル化されたキロノバは、主に円盤風による放出物が支配的である。硬いDD2 EOSを持つ系は、軟らかいSLy EOSを持つ系よりも、より大規模な円盤を持つため、より明るいキロノバを生成する。距離200 Mpcにおいて、ほとんどの系はVera C. Rubin天文台によって4日以内に、またDECamによって $r$ および $i$ バンドで2日以内に検出可能である。

意義と主張
本論文は、対称的なBHNS合体におけるNRシミュレーションのパラメータ空間における決定的な空白を埋めるものである。主な意義は、現在の重力波波形モデルが等質量および近等質量比において不十分な精度であることを示し、GW230529のような将来の検出におけるパラメータ推定にバイアスをもたらす可能性があることを実証した点にある。本研究は、残留質量に関する解析的フィッティング公式を検証する一方で、高潮汐変形能を持つ系における残留スピン予測の限界を浮き彫りにした。さらに、円盤力学の詳細な分析は、複雑な回転則と振動駆動の降着変調を明らかにしており、$g$ モードがガンマ線バーストの光度曲線に観測可能な痕跡を残す可能性を示唆している。最後に、キロノバ・モデリングは、対称的なBHNS合体が現在の広視野サーベイにとって有望な電磁波フォローアップの対象であり、従来の観測で支配的であった非対称な系と比較して、検出確率が高いことを示している。