A New Scaling of Neutron Star Tidal Deformability for Directly Probing the Core Equation of State

技術要約：中性子星の核となる状態方程式を直接探るための、潮汐変形能の新たなスケーリング

問題提起
中性子星（NS）の無次元潮汐変形能（$\Lambda$）は、インスパイラル（接近）過程における重力波観測から推測されるものであり、従来は核飽和密度の約2倍付近の物質の圧力を制約するために用いられてきた。その結果、中性子星コアの深部における状態方程式（EOS）は、従来の恒星構造および潮汐応答方程式に基づく解析では、インスパイラル期の重力波観測からはほとんどアクセスできない領域となっていた。将来のポスト・マージャー（合体後）の観測は、高密度EOSをより直接的に探ることができる可能性があるが、既存のインスパイラル期のデータに、これまで認識されていなかったコアEOSに関する情報が含まれているかどうかは未解決の問題である。高密度EOSの推論に関する現在の取り組みは、パラメータ化されたモデルに大きく依存しており、その制約はモデル依存性を引き継いでしまう。本論文は、$\Lambda$ 自体が、コアEOSに関する堅牢でモデルに依存しない情報を含んでいるかどうかを検証するものである。

手法
著者らは、無次元の恒星構造および相対論的潮速応答方程式の摂動解析を行うために、IPAD-TOV（Intrinsic and Perturbative Analyses of the Dimensionless TOV equations）フレームワークを採用している。

1. 理論的導出： 本研究では、特定のEOSモデルを導入することなく、無次元の潮汐変形能 $\Lambda = 2k_2/3\xi^5$（ここで $k_2$ は四重極潮汐ラブ数、$\xi$ はコンパクトネス）の基礎となる方程式から直接、スケーリング関係を導出している。解析は、潮汐応答関数 $y(r)$ の中心近傍の挙動に焦点を当て、簡約変数を用いた展開を利用している。
2. 主要変数： 導出過程において、スケーリングは主に以下の2つの無次元量によって支配されることが特定された。
   • 中心部のEOSパラメータ $X \equiv P_c/\epsilon_c$ （中心圧力を中心エネルギー密度で除したもの。これは中心部までのエネルギー密度平均音速の二乗を表す）。
   • 対数安定性勾配 $\Psi = 2 d\ln M_{NS} / d\ln \epsilon_c$。
3. スケーリング関数： 著者らは、$\Lambda$ の本質的なEOS依存性を包含する、固有の潮汐応答関数 $D(X, \Psi)$ を定義している。
4. 検証： 導出されたスケーリングの普遍性を確立するため、核理論、重イオン衝突データ、および天体観測に基づき校正された、非線形ワレカ・モデルおよびメタモデル記述に基づく広範なEOSアンサンブルを構築した。さらに、相転移、ハイペロン、または非局在化したクォークを含む文献中の33種類の代表的な微視的EOSに対して、このスケーリングを検証している。

主な貢献と結果

1. 新たなスケーリング関係の発見： 本論文は、潮汐変形能 $\Lambda$ と中心部EOSパラメータ $X$ の間の、驚くべき、かつEOSにほとんど依存しないスケーリング関係を明らかにした。この関係は $\Lambda \approx a D(X, \Psi) + b$ と表される。
   • 標準的な $1.4 M_\odot$ の中性子星については、適合の結果 $\Lambda_{1.4} \approx 6.41 D(X, \Psi) + 23.50$ となり、決定係数は $R^2 \approx 0.961$ である。
   • 最大質量安定構成（TOV中性星、$\Psi=0$ の場合）については、関係式は $\Lambda_{TOV} \approx 2.36 D_{TOV} + 3.26$ であり、$R^2 \approx 0.958$ である。
2. コアパラメータの直接抽出： このスケーリングをGW170817から推測された潮汐変形能（$190 \lesssim \Lambda_{1.4} \lesssim 580$）に適用することで、中心部EOSパラメータ $X$ を直接決定することが可能となった。標準的な中性子星に対する解析では、$0.133 \lesssim X \lesssim 0.185$ という結果が得られた。これにより、特定のEOSパラメータ化に依存することなく、潮汐変形能の測定値からコアEOSの情報を直接抽出できることが示された。
3. 中心圧力およびエネルギー密度の制約： この新しい $\Lambda$-$D$ スケーリングを、以前に確立された質量スケーリング関係と組み合わせることで、標準的な中性子星の中心エネルギー密度（$\epsilon_c$）および中心圧力（$P_c$）を制約している。結果は、標準的な中性子星において $2.48 \lesssim \epsilon_c/\epsilon_0 \lesssim 4.55$ および $49.4 \lesssim P_c \lesssim 125.9$ MeV/fm$^3$ であることを示している。
4. 潮汐変形能の下限値： 最大質量安定中性星について、このスケーリングは、ロバストな下限値 $\Lambda_{TOV} \gtrsim 9.2$（因果律の限界 $X \lesssim 0.381$ を仮定した場合）を予測している。これは、最もコンパクトな安定中性星であっても、依然として測定可能な変形能を持ち、ブラックホール（$\Lambda_{BH} = 0$）とは明確に区別されることを定量的に示している。
5. 音速プロファイル： 推測された値は、標準的な中性子星における平均的な超核密度での硬さ（$s^2 \sim 0.3\text{--}0.4$）が、共形限界（$s^2 = 1/3$）に接近しているか、あるいは一時的に超えている可能性を示唆しており、中性子星コアにおける物質の著しい硬化を示している。

意義
本論文は、インスパイラル期の重力波観測量と、強重力領域における超高密度物質の微物理学との間の直接的なつながりを確立したと主張している。主な意義は、無次元潮汐変形能が、核飽和密度の約2倍付近の物質のプローブであるだけでなく、コアEOSに関する直接的な情報を持っていることを示した点にある。相転移やエキゾチックな物質を含む多様なEOSモデルにわたるスケーリングの堅牢性は、相対論的な恒星構造および潮汐応答方程式が、本質的な普遍性の構造を持っていることを示唆している。この構造により、微視的なEOSへの依存性が、少数のマクロな無次元量（$X$ および $\Psi$）へと効果的に圧縮されている。したがって、導出されたスケーリングは、現在および将来の重力波データを用いて、中性子星コアの熱力学的性質を探索するための、新しい、かつモデルにほとんど依存しない枠組みを提供するものである。