[論文レビュー] On the Sound Speed in Neutron Stars

Sinan Altiparmak, Christian Ecker+1 more|arXiv (Cornell University)|Mar 28, 2022

Experimental and Theoretical Physics Studies被引用数 10

ひとこと要約

本研究では、核理論、摂動的QCD、および天体的制約を組み合わせることで、中性子星における1000万個以上の連続的音速状態方程式（EOS）を構築した。その結果、音速が1/3未塔（c²ₛ < 1/3）である状態は統計的に稀（0.03%）であり、星内部でc²ₛ > 1/3であるモデルが強く支持される。1.4 M☉の中性子星の95%信用区間推定値はR₁.₄ = 12.42⁺⁰.⁵²₋₀.⁹⁹ km、2.0 M☉の中性子星ではR₂.₀ = 12.12⁺¹.¹¹₋₁.²³ kmであり、GW170817のための潮汐歪み係数の95%信用区間推定値はŨ₁.₁₈₆ = 485⁺²²⁵₋₂₁₁であった。

ABSTRACT

Determining the sound speed $c_s$ in compact stars is an important open question with numerous implications on the behaviour of matter at large densities and hence on gravitational-wave emission from neutron stars. To this scope, we construct more than $10^7$ equations of state (EOSs) with continuous sound speed and build more than $10^8$ nonrotating stellar models consistent not only with nuclear theory and perturbative QCD, but also with astronomical observations. In this way, we find that EOSs with sub-conformal sound speeds, i.e. with $c^2_s &lt; 1/3$ within the stars, are possible in principle but very unlikely in practice, being only $0.03\%$ of our sample. Hence, it is natural to expect that $c^2_s &gt; 1/3$ somewhere in the stellar interior. Using our large sample, we obtain estimates at $95\%$ credibility of neutron-star radii for representative stars with $1.4$ and $2.0$ solar masses, $R_{1.4}=12.42^{+0.52}_{-0.99}\,{ m km}$, $R_{2.0}=12.12^{+1.11}_{-1.23}\,{ m km}$, and for the binary tidal deformability of the GW170817 event, $ ildeΛ_{1.186}=485^{+225}_{-211}$. Interestingly, our lower-bounds on the radii are in very good agreement with the prediction derived from very different arguments, namely, the threshold mass. Finally, we provide simple analytic expressions to determine the minimum and maximum values of $ ildeΛ$ as a function of the chirp mass.

研究の動機と目的

• 実際の物理的および観測的制約下で、中性子星における非共形的音速（c²ₛ < 1/3）が生じる可能性を特定すること。
• 質量が2 M☉を超える中性子星が存在することを踏まえ、星内部で音速が1/3を超えることが必須であるかどうかを評価すること。
• 大規模で物理的に整合性のある状態方程式のサンプルを用いて、中性子星の半径および潮汐歪み係数の堅牢で観測的に制約された推定値を提供すること。
• 区分的多項式モデルに見られる不連続性を回避するため、音速パラメータ化に基づく、状態方程式の性質を推定する統計的に堅牢なフレームワークを構築すること。

提案手法

• Annalaら（2020）のパラメータ化法を用いて、熱力学的および因果的安定性を保証する連続的音速を持つ10⁷個以上の状態方程式（EOS）を構築した。
• 核理論（n < 0.5n₀でBPSモデル）、有効場理論（n ≲ n₀）、摂動的QCD（n ≫ n₀）の制約を統合し、物理的整合性を確保した。
• 最大7段階の密度セグメントを用いたパッチワーク的手法を採用し、物理的に不適切な解を除外するフィルタリングを施して収束性を向上させた。
• 全EOSサンプルを対象に非回転中性子星の星体モデルを計算し、質量-半径関係および潮汐歪み係数関係を導出した。
• ベイズ的信用区間を用いて、統計的ゆらぎを最小限に抑える大規模サンプルを用いて、半径および潮汐歪み係数の95%信頼区間を推定した。
• 音速の確率密度関数（PDF）を分析し、単調増加・非単調・非共形的・超共形的など、さまざまな状況の統計的尤もらしさを評価した。

実験結果

リサーチクエスチョン

• RQ1 実際の物理的および観測的制約下で、中性子星における音速が共形限界（c²ₛ < 1/3）未満である確率はどの程度か？
• RQ2 質量が2 M☉を超える中性子星が存在することは、星内部のある密度で音速が1/3を超える必要があることを示唆するか？
• RQ3 GW170817および中性子星質量からの観測的制約が、中性子星の半径および潮汐歪み係数の信用区間の境界にどのように影響するか？
• RQ4 高密度領域におけるQCD領域の境界条件に、音速プロファイルの性質（例：単調性、ピーク位置）がどの程度依存するか？

主な発見

• 非共形的音速（c²ₛ < 1/3）は、検討された10⁷個以上のEOSのうちたった0.03%にしか現れず、自然に実現されにくいことが示された。
• 1.4 M☉の中性子星の95%信用区間推定値はR₁.₄ = 12.42⁺⁰.⁵²₋₀.⁹⁹ kmであり、2.0 M☉の中性子星ではR₂.₀ = 12.12⁺¹.¹¹₋₁.²³ kmであった。
• GW170817二重星イベントの潮汐歪み係数の95%信用区間推定値はŨ₁.₁₈₆ = 485⁺²²⁵₋₂₁₁であり、現在の観測的制約と整合的であった。
• 本分析から得られた中性子星半径の下限値は、閾値質量議論から導かれる予測と非常に良好に一致しており、独立した理論的手法間での一貫性が示された。
• 音速が密度とともに増加する単調なプロファイル（c²ₛが密度とともに増加）では、1.99 M☉を超える星を支持できないため、このようなモデルは観測された質量の大きな中性子星と整合しない。
• 中性子星密度付近に音速PDFのピークが現れる主な要因は、共形値自体ではなく、因果的制約に従って摂動的QCD極限に近づく必要があることによる。