[论文解读] Nontrivial features in the speed of sound inside neutron stars

Débora Mroczek, M. Coleman Miller+2 more|arXiv (Cornell University)|Sep 5, 2023

Pulsars and Gravitational Waves Research被引用 5

一句话总结

本文提出一种基于高斯过程的方法，用于在中子星核心中生成具有非平凡声速特征（如平台、凸起和相变）的方程态（EoS）模型。通过结合X射线、引力波和微扰QCD数据进行贝叶斯推断，发现全局声速最大值与当前观测一致但非必需，暗示可能存在奇异相，但目前数据尚无法对其加以约束。

ABSTRACT

Measurements of neutron star masses, radii, and tidal deformability have direct connections to nuclear physics via the equation of state (EoS), which for the cold, catalyzed matter in neutron star cores is commonly represented as the pressure as a function of energy density. Microscopic models with exotic degrees of freedom display nontrivial structure in the speed of sound ($c_s$) in the form of first-order phase transitions and bumps, oscillations, and plateaus in the case of crossovers and higher-order phase transitions. We present a procedure based on Gaussian processes to generate an ensemble of EoSs that include nontrivial features. Using a Bayesian analysis incorporating measurements from X-ray sources, gravitational wave observations, and perturbative QCD results, we show that these features are compatible with current constraints. We investigate the possibility of a global maximum in $c_s$ that occurs within the densities realized in neutron stars -- implying a softening of the EoS and possibly an exotic phase in the core of massive stars -- and find that such a global maximum is consistent with, but not required by, current constraints.

研究动机与目标

• 开发一种灵活的非参数化方法，用于生成具有复杂非单调声速（cₛ）行为的方程态（EoS）模型。
• 检验当前天体物理和理论约束是否能排除或支持cₛ中非平凡特征（如相变或平台）的存在。
• 评估在中子星密度范围内cₛ中全局最大值的可行性，该最大值可能标志EoS软化及大质量恒星中奇异物质的存在。
• 整合来自X射线观测、引力波事件（如GW170817）以及高密度下微扰QCD的约束条件。
• 评估结果对高密度QCD行为不确定性的敏感性，并判断当前数据是否能够区分具有与不具有奇异相变的EoS。

提出的方法

• 采用多尺度高斯过程（mGP）框架，非参数化地建模压力随能量密度变化的方程态（EoS），以灵活表示cₛ² = dp/dε。
• 使用具有多条相关长度的分层先验结构，以捕捉不同密度区域中cₛ²的尖锐转变、振荡和平台等特征。
• 通过贝叶斯模型比较，整合来自X射线定时（如NICER）、引力波事件（如GW170817）和潮汐可变形性测量的观测约束。
• 在高密度区域（n_B ≳ 40n_sat）引入微扰QCD（pQCD）约束，使cₛ² → 1/3，以确保与高密度QCD预期的一致性。
• 通过马尔可夫链蒙特卡洛（MCMC）采样探索EoS实现后验分布，评估cₛ²全局最大值等特征的可能性。
• 通过评估违反pQCD极限的采样EoS比例，量化EoS与pQCD约束之间的张力，发现仅3.5%存在部分张力，0.0083%完全不一致。

实验结果

研究问题

• RQ1 当前天体物理和理论约束能否排除中子星核心内声速的非单调行为？
• RQ2 在中子星密度范围（2–6 n_sat）内，声速（cₛ²）的全局最大值是否与现有观测和理论数据一致？
• RQ3 高密度下微扰QCD约束在多大程度上影响具有奇异特征的EoS后验分布？
• RQ4 当前数据在多大程度上能区分cₛ²中的一阶相变、交叉转变或平滑行为的EoS？
• RQ5 在当前观测不确定性与模型灵活性下，潜在相变在EoS中的统计显著性如何？

主要发现

• 在中子星密度范围内，cₛ²的全局最大值与当前观测和理论约束一致，但并非必需。
• 仅0.0083%的采样EoS在高密度下违反微扰QCD约束，3.5%存在部分张力，表明pQCD对后验分布影响极小。
• 贝叶斯分析表明，cₛ²中的非单调特征（如平台、凸起和相变）未被当前数据（包括X射线和引力波观测约束）所排除。
• 该方法成功生成了具有复杂cₛ²行为的EoS，包括一阶相变和类似交叉转变的特征，且无需假设特定函数形式。
• 即使拥有约100个未来中子星观测的目录，对强一阶相变的支持的贝叶斯因子也仅能达到约10:1，表明当前数据灵敏度有限。
• 结果支持需要来自NICER、LIGO/Virgo/KAGRA以及未来核物理实验的更精确数据，以明确探测中子星核心中的奇异相。