Extracting the Central Sound Speed of Neutron Stars through Intrinsic Stellar Structure Scaling with Reduced Equation-of-State Model Dependence

技術要約：固有の恒星構造スケーリングを通じた中性子星の中心音速の抽出

問題提起
観測データから中性子星（NS）の核の性質を推論することは、仮定される状態方程式（EOS）に対して極めて敏感である。異なる高密度物質モデルは、同一の核物理学的および天体物理学的制約（例：大質量パルサー質量、GW170817からの潮汐変形能、NICERによる半径測定）を満たすことができる一方で、予測される内部圧力、エネルギー密度、および音速は大きく異なる場合がある。これは、中性子星物理学における逆問題に曖昧さを生じさせる。すなわち、推論された核の性質が、高密度物質の本質的な性質を反映しているのか、それとも解釈に使用された特定のEOSモデルの仮定を反映しているのかが不明確なのである。具体的には、中心部の音速の二乗（$s^2$）およびEOSパラメータ（$\phi = P/\epsilon$）の挙動は、異なる微視的な記述の間で大きく変化するが、マクロな観測量は一致している。

手法
この曖昧さに対処するため、著者らはモデルの多様性と固有の恒星構造解析を組み合わせた二段構えのアプローチを採用している。

1. 多様なEOSアンサンブルの構築：
   著者らは、修正された$\omega$メソンの自己相互作用項（$(g_\omega^2 \omega_\mu \omega^\mu)^q$ の形式）を含む、一般化された非線形ワレカ・モデルを利用している。べき指数 $q$ を変化させ（具体的には $q = 3/2, 2, 3, 4, 5$）、自己相互作用のない参照ケース（$c_\omega = 0$）を含めることで、6つの異なるEOSファミリーを生成した。各ファミリーについて、1,000個のEOSを独立してサンプリングし、以下の広範な制約を満たすようにフィルタリングした：

• 核物質の制約： 飽和密度（$\rho_0$）、結合エネルギー（$E_0$）、非圧縮率（$K_0$）、スキューネス（$J_0$）、対称エネルギー（$E_{sym}$）、およびその傾き（$L$）。
• 天体物理学的制約： 最大質量（$M_{TOV} \gtrsim 2.01 M_\odot$）および潮汐変形能（$\Lambda_{1.4} \lesssim 580$）。
  この設定により、すべてのファミリーが同様のマクロな恒星特性を記述する一方で、異なる高密度挙動と中心音速プロファイルを示すことが保証されている。

2. 無次元TOV方程式の固有および摂動解析（IPAD-TOV）：
   特定のEOSに依存する代わりに、著者らは、無次元のトルマン・オッペンハイマー・ヴォルコフ（TOV）方程式から直接導出された解析的推定値 $\Theta(X, \Psi)$ を用いる。この推定値は、以下の2つの量を用いて中心音速の二乗（$s_c^2$）を探索する：

• $X \equiv P_c/\epsilon_c$：中心部のEOSパラメータ（中心圧力とエネルギー密度の比）。
• $\Psi \equiv 2 d\ln M_{NS} / d\ln \epsilon_c$：中心エネルギー密度に対する恒星質量の対数傾き。
  $\Theta$ は、恒星変数を中心付近で展開し、劣次項の質量項を保持することによって導出され、異なるEOS記述間での中心特性を比較するための統一的な解析的基盤を提供する。

主な貢献と結果

• モデルの曖昧さの証明： 本研究は、6つのEOSファミリーが同じ制約を満たしているにもかかわらず、著しく異なる中心特性を予測することを裏付けている。標準的な$1.4 M_\odot$の中性子星において、中心音速の二乗（$s_c^2$）は、$\approx 0.36$（$q=3/2$の場合）から $\approx 0.50$（$c_\omega=0$の場合）の範囲にある。$c_\omega=0$のファミリーは漸近的に硬い挙動（$s^2 \to 1$）を示すが、$c_\omega$が有限のファミリーは、$q$に応じて異なる漸近限界（例：$1<q<2$では $s^2 \to 1/(2q-1)$、$q \ge 2$では $s^2 \to 1/3$）を示す。
• 内部音速ピークのメカニズム： 著者らは、純粋な核子フレームワーク内において、線形ベクトル反発と非線形$\omega$メソンの自己相互作用の競合が、内部の音速プロファイルのピークを生じさせ得ることを特定した。このピークの位置と大きさは $q$ に依存し、高い $q$ 値はピークを高密度側へシフトさせ、標準的な中性子星の核内でアクセス可能にする。
• 普遍的なスケーリングの発見： 多様な基礎的相互作用にもかかわらず、著者らは解析的推定値 $\Theta(X, \Psi)$ と、数値的に計算された $s_c^2$ との間に、全6ファミリーにわたる緊密な線形相関を見出した。決定係数（$R^2$）は0.94から0.99の範囲にある。
• モデル依存性を低減した抽出： ファミリー固有のスケーリング関係（$s_c^2 \approx \alpha_j \Theta + \beta_j$）の分散を最小化することにより、特定のEOSファミリーに対する感度が低い、共通の中心音速およびEOSパラメータの推定値を抽出した。
  • 全6ファミリーを含む場合： $s_c^2 \approx 0.523^{+0.044}_{-0.044}$ および $X \approx 0.173^{+0.020}_{-0.018}$。
  • 硬い $c_\omega=0$ ファミリーを除外した場合： $s_c^2 \approx 0.496^{+0.092}_{-0.092}$ および $X \approx 0.167^{+0.030}_{-0.029}$。
  • これらの値は、それぞれ約$11.73$ kmおよび$11.92$ kmの標準的な中性子星半径に対応する。
• 整合性チェック： 抽出された値は、全EOSサンプルに基づく結合スケーリング（$R^2 \approx 0.95$）と一致しており、$\Lambda$-$D(X, \Psi)$ 相関に基づく先行研究とも整合している。推論された $s_c^2 \sim 0.5$ は共形限界である $1/3$ を明らかに上回っており、標準的な中性子星の内部核が硬いことを示唆している。

意義と主張
本論文は、現在の観測的制約では異なるEOSモデルやその特定の予測を一意に区別することはできないが、共有された恒星構造の挙動がTOV方程式にコード化されており、モデル依存性を低減した形で中心特性を明らかにできると主張している。

著者らは、本アプローチが単一の「正しい」EOSを特定したり、モデル依存性を完全に排除したりすることを目的としているのではないことを強調している。むしろ、恒星構造方程式に内在する情報（具体的には、大域的な恒星の応答 $\Psi$ と中心の状態変数 $X$ の間のスケーリング関係）を用いることで、中性子星の核に関する頑健な物理情報を分離できることを示している。モデルの多様性をエラーの源ではなく情報の源として扱うことで、IPAD-TOVフレームワークは、採用されたEOS記述のアーティファクトから、固有の核特性を分離するための一般的な経路を提供する。結果は、標準的な中性子星が、EOSの究極的な高密度漸近挙動には比較的影響を受けにくい密度領域を探索しており、共通の値 $s_c^2 \approx 0.5$ に収束していることを示唆している。