[论文解读] 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），发现亚共形声速（c²ₛ < 1/3）在统计上极为罕见（仅0.03%），强烈支持星体内部声速满足c²ₛ > 1/3的模型，1.4 M☉中子星的95%可信估计半径为R₁.₄ = 12.42⁺⁰.⁵²₋₀.⁹⁹ km，2.0 M☉中子星为R₂.₀ = 12.12⁺¹.¹¹₋₁.²³ km，对GW170817事件的潮汐去形因子估计为Ũ₁.₁₈₆ = 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）出现的可能性。
• 在存在大质量中子星（M ≳ 2 M☉）的前提下，评估星体内部声速是否必须超过1/3。
• 利用大规模、物理解释一致的方程样本，提供稳健且受观测约束的中子星半径与潮汐去形因子估计。
• 开发一种统计上稳健的框架，通过声速参数化方法估计EOS特性，避免分段多相模型中常见的不连续性。

提出的方法

• 采用Annala等人（2020）的参数化方法，构建超过10⁷条具有连续声速的方程（EOS），确保热力学与因果稳定性。
• 结合核理论（n < 0.5n₀时的BPS模型）、有效场论（n ≲ n₀）和微扰QCD（n ≫ n₀）的约束，以保证物理一致性。
• 采用分段拼接方法，最多使用七个密度段建模EOS，通过筛选去除非物理解，并提升收敛性。
• 对整个EOS样本中的非旋转中子星计算星体模型，推导质量-半径关系与潮汐去形因子关系。
• 应用贝叶斯可信区间，利用大规模样本最小化统计波动，估计半径与潮汐去形因子的95%置信区间。
• 分析声速概率密度函数（PDF），以确定不同情形（单调、非单调、亚共形或超共形）的统计可能性。

实验结果

研究问题

• RQ1 在真实的物理与观测约束下，中子星中声速低于共形极限（c²ₛ < 1/3）的统计可能性是多少？
• RQ2 大质量中子星（M ≳ 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极限所驱动，而非共形值本身。