Convective stability analysis of massive neutron stars formed in binary mergers

技術要約：連星合体によって形成された大質量中性子星の対流安定性解析

問題提起
連星中性子星（BNS）合体後の進化は、その結果として生じる大質量中性子星（NS）の力学によって支配される。この進化における極めて重要な側面は、流体力学的不安定性、特に対流不安定性が発達する可能性であり、これが大規模な流体運動を駆動し、熱構造を変化させ、重力波（GW）モードを励起する可能性がある。先行研究では相反する結論が得られている。シュワルツシルトの判定基準を用いたシミュレーションの中には、大質量NSが大規模な対流不安定性を発達させ、その結果として観測可能な慣性モードを励起することを示唆したものがある [74, 75]。一方、ニュートリノ放射を組み込んだ研究では、安定な成層構造が見出されている [55]。これまでの解析における主要な限界は、回転による復元力が頻繁に無視されている点にある。合体残骸は急速かつ差動的な回転を示すため、回転の寄与（しばしばソルベルグ・ホーランドの判定基準として記述される）を考慮することが不可欠である。本論文は、完全相対論的流体解析シミュレーションを実行し、熱的な浮力と回転の効果の両方を考慮した包括的な安定性基準を適用することで、これらの相違を解決することを目的とする。

手法
著者らは、SACRA-MPI数値相対論コードを用いて、等質量BNS系の合体を合体後約100 msの期間にわたってシミュレートする。

• 状態方程式（EOS）： 衝撃加熱をモデル化するために$\Gamma$則の熱成分（$\Gamma_{\text{th}} = 1.8$）で補強された3つの区分的多項式EOS（APR4, SLy4, MPA1）を使用する。
• 数値設定： シミュレーションでは、軌道面に対して鏡像対称性と2:1の細分化を持つデカルト座標の入れ子状格子を利用する。収束性を確認するため、APR4モデルに対して3つの解像度をテストする。
• 安定性解析： 著者らは、高温の差動回転星に対する相対論的ソルベルグ・ホーランド安定性基準（基準Iおよび基準II）を導出し、適用する。これらの基準は、熱的浮力（シュワルツシルト・判別式）と回転による復元力（比角運動量勾配）の相互作用を取り入れている。解析は、比エントロピーと比角運動量の角度平均プロファイルに対して行われる。
• モード解析： 流体モードの励起とGW放射を分析し、特に$m=1$の一腕モードと、四重極$f$モードの減衰後に発生する可能性のある$m=2$慣性モードの励起に焦点を当てる。

主要な貢献

1. 相対論的安定性基準の導出： 本論文は、シュワルツシルトの判定基準のみに依存していた研究で無視されていた回転による復元力を明示的に考慮し、相対論的ソルベルグ・ホーランド基準（式46および47）を大質量合体後NSに対して厳密に適用している。
2. 独立したシミュレーション・フレームワーク： 著者らは、直接的な一対一の比較を避けて物理的な調査を自己整合的に行うために、先行する相反する研究とは異なる数値設定（EOSの選択、解像度、流体スキーム）を用いた独立した解析を実施している。
3. 対称性と運動量保存のテスト： 本研究は、数値コードが人工的に$\pi$対称性を破っていないことを検証するための対称保持テストを導入している。さらに、$m=1$モードの振幅の成長と線形運動量保存の破れの間の相関関係を調査しており、これは物理的な不安定性と数値的なアーティファクトを区別するための重要なチェックである。

結果

• 対流安定性： シミュレーションは、大規模な対流不安定性の証拠を示さない。大質量NS内では、角度平均された比エントロピーと比角運動量の両方が外側に向かって増加している。局所的な領域はシュワルツシルトの判定基準のみでは不安定に見える場合があるが、回転の効果（特に比角運動量の外向きへの増加）を考慮すると、これらの領域は安定化される。全体的な安定性基準（基準IおよびII）は、星の大部分において満たされている。
• 慣性モード： 先行研究 [74, 75] の知見とは対照的に、四重極$f$モードの減衰後に観測可能な慣性モードの励起は検出されなかった。GW信号はシミュレーション全体を通じて$f$モードによって支配されており、キロヘルツ周波数帯における離散的な慣性モードの出現は見られない。
• $m=1$ 不安定性： $m=1$の一腕モードが増幅し飽和することが観察され、これは先行文献と一致している。しかし、著者らはこのモードの振幅の成長と線型運動量保存の破れとの間に強い相関があることを見出した。これは、これらの特定のシミュレーションにおける一腕モードの励起が、純粋な物理的不安定性ではなく、数値的な起源を持つ可能性を示唆している。
• 熱および回転プロファイル： 大質量NSは20〜30 ms以内に準平衡状態に落ち着く。コアは低温で低速回転しており、エンベロープはスパイラルアームの相互作用によって加熱され、差動回転している。比角運動量は外側に向かって増加しており、強力な安定化遠心力を提供している。

意義と主張
本論文は、回転による復元力の導入が、合体後の大質量NSの安定性評価を根本的に変えることを主張している。著者らは、安定な構成を示す自らの結果が、大規模な対流とそれに続く慣性モードの励起を予測した先行研究とは相容れないものであると断言している。彼らは、シミュレーションにおける対流不安定性の欠如が、励起された慣性モードの不在を説明すると論じている。

$m=1$モードに関しては、著者らは「一腕モードの励起が数値的な起源を持つ可能性を否定できない」とし、線形運動量非保存との相関に基づき慎重な立場をとっている。彼らは、コードが対称保持テストに合格している一方で、メッシュ細分化シミュレーションにおける線形運動量の保存が依然として課題であり、一腕不安定性の物理的性質を確定させるためにはさらなる調査が必要であることを強調している。

本研究は、微物理学（組成勾配）や磁場を取り入れた将来の研究が、より決定的な結論を提供するために必要であると結論付けているが、現在の流体力学シミュレーションの範囲内では、連星合体によって形成される大質量NSは対流的に安定していると考えられる。