Maximal mass of neutron stars constrained by neutron star observations

想象宇宙中充满了一种名为“中子星物质”的神秘超致密材料。它重得惊人，仅一茶匙的重量就相当于一座山。长期以来，物理学家一直在试图弄清这种物质的“游戏规则”——具体而言，就是它有多硬或多易被压缩。这套规则被称为状态方程（EOS）。

本文试图回答的核心问题是：中子星在坍缩成黑洞之前，绝对能达到的最大质量是多少？

以下是作者如何解开这一谜题的故事，以通俗易懂的方式呈现：

1. 两个起点（两种“配方”）

为了找出规则，科学家们从两种不同的“配方”出发，描述这种致密物质在较低密度下的行为。你可以将其理解为关于“原料”如何混合的两种不同理论：

配方 A（SFHo）： 一种“较软”的配方，意味着物质相对更容易被压缩。
配方 B（DD2）： 一种“较硬”的配方，意味着物质更抗拒被压缩。

他们知道这些配方在密度尺度的“起点”处表现良好，但不知道在中子星中心那种极端超高密度下会发生什么。为了填补这一空白，他们利用数学上的“桥梁”，将这些配方与我们在最高可能能量下所知的粒子物理知识连接起来。

2. 侦探工作（利用真实线索）

作者没有仅仅依靠猜测，而是像侦探一样行事。他们拿这两种配方，用望远镜和引力波探测器收集到的真实世界线索进行检验。他们使用了一种特殊的统计方法（称为贝叶斯加权），来观察哪些版本的配方通过了测试。

以下是他们使用的线索：

“大撞击”（GW170817）： 当两颗中子星相互撞击时，它们在时空中激起涟漪。这些涟漪的行为方式告诉科学家这些星体有多“软”。
“探照灯”（NICER）： 一台空间望远镜拍摄了旋转中子星上热点的照片。通过测量这些星体看起来的大小及其重量，他们获得了直接的尺寸与重量之比。
“轻量级”候选者（HESS J1731–347）： 一个非常小且轻的天体，可能是一颗中子星。
“重量级”候选者（GW190814）： 一个神秘的天体，比大多数中子星重，但比大多数黑洞轻。科学家们问道：这实际上会不会是一颗超重的中子星？

3. 结果：线索告诉了他们什么

科学家们将这两种配方代入这些线索中，并观察结果。

质量上限（最大质量）：

令人惊讶的是： 无论他们使用哪种起始配方（软或硬），结果差别不大。真实世界的线索如此有力，迫使两种配方都得出了相同的结论。
裁决： 当他们使用最可靠的线索（“大撞击”和“探照灯”）时，中子星所能承载的最大质量约为太阳质量的 2.2 到 2.3 倍。
“重量级”的转折： 如果他们假设那个神秘的重天体（GW190814）是一颗中子星，那么上限会跃升至约太阳质量的 2.6 到 2.7 倍。然而，这与来自“大撞击”的“软度”线索相冲突，使得情况变得棘手。

尺寸上限（半径）：

差异： 与质量不同，恒星的尺寸确实取决于他们使用了哪种起始配方。
裁决： “软”配方预测的半径约为11.8 公里，而“硬”配方预测的半径约为12.4 公里。
最佳区间： 当结合所有最佳线索时，这些恒星最可能的尺寸约为12 公里（上下浮动 1 公里）。

4. 大局观

该论文得出结论，通过观察“终点”（最重和最大的可能恒星），并利用混合的真实天文数据，我们可以缩小宇宙中最致密物质的规则范围。

质量： 宇宙似乎对中子星能有多重设有一个“速度限制”，稳定在2.2 到 2.3 倍太阳质量左右。这与我们迄今为止实际观测到的最重中子星相符。
尺寸： 它们的大小大致相当于一个小城市，直径约12 公里。
要点： 真实世界的观测（线索）远比理论上的起始假设强大得多。无论你从哪种理论出发，来自恒星本身的数据都会迫使答案收敛到相同的数值。

简而言之，宇宙给了我们一个非常明确的答案：中子星可以变得极其沉重，但存在一个硬性上限，而且相对于它们的重量而言，它们小得令人惊讶。

以下是 Gábor Kasza 和 György Wolf 所著论文《中子星观测对中子星最大质量的约束》的详细技术总结。

1. 问题陈述

核物理与重离子物理领域的核心挑战在于确定中等密度和化学势下强相互作用物质的状态方程（EOS）。虽然低密度物质受核实验和格点 QCD 的约束，高密度物质受微扰 QCD（pQCD）的约束，但与中子星核心相关（数倍于核饱和密度 $\rho_0$ ）的区域仍知之甚少。

本文解决的具体问题是确定中子星的最大质量（ $M_{TOV}$ ）及其对应的半径（ $R_{TOV}$ ）。质量 - 半径（ $M(R)$ ）序列的这些端点对理解超核密度下 EOS 的刚度至关重要。作者旨在量化当前的多信使观测如何约束这些端点，以及结果在多大程度上依赖于所选的基础强子基准模型。

2. 方法论

A. 混合 EOS 的构建

作者通过连接低密度强子物质和高密度夸克物质，构建了因果混合 EOS 族：

低密度区域： 使用两个代表性的强子模型作为基准：
- SFHo： 相对论平均场模型（较软的 EOS，不可压缩性 $K=245$ MeV）。
- DD2： 较硬的强子 EOS。
高密度区域： 基于 $U(3) \times U(3)$ 手征对称性的扩展线性西格玛模型（eLSM），包含组分夸克和介子九重态。该模型被约束为在低化学势下匹配格点 QCD 热力学结果，在高化学势下匹配 pQCD 结果。
插值： 为了连接这两个区域，对能量密度随重子密度的函数应用五阶多项式插值。这确保了热力学一致性，以及压力和声速（ $c_s$ ）在过渡区域的连续性，避免了大的一阶相变。
渐近约束： 夸克 - 介子 EOS 被要求在高渐近密度（ $\mu \approx 2.6$ GeV）下平滑连接到 pQCD 结果。

B. 贝叶斯框架

本研究采用贝叶斯加权框架来分析 $M(R)$ 序列端点的概率分布。

先验： 生成一个广泛的因果混合 EOS 集合，涵盖超核密度下的刚度范围。
似然加权： 根据观测数据对集合进行加权，生成 $M_{TOV}$ 和 $R_{TOV}$ 的后验概率分布。

C. 观测约束

以下约束被依次并组合应用：

pQCD 约束： 确保 EOS 在保持因果性（ $c_s < 1$ ）的同时，趋近于 N3LO pQCD 参考点（ $\mu=2.6$ GeV, $\rho=6.47$ fm $^{-3}$ , $p=3823$ MeV/fm $^3$ ）。
GW170817： 潮汐形变约束（ $\tilde{\Lambda} < 720$ ）以及更严格的组合约束（" $\Lambda_{190}$ "： $70 < \Lambda_{1.4} < 580$ 且 $9.1 < R_{1.4} < 12.8$ km）。
NICER： 五颗脉冲星（J0030+0451, J0740+6620, J0614–3329, J1231–1411, J0437–4715）的质量 - 半径测量，精度约为 10%。
HESS J1731–347： 一个极低质量中子星候选体（ $M \approx 0.77 M_\odot$ ），由于解释存在不确定性，处理时需谨慎。
GW190814： “质量间隙”天体（ $M \approx 2.59 M_\odot$ ）。作者通过应用以 $2.59 M_\odot$ 为中心的 Gaussian 约束，检验该天体为中子星的假设。

3. 主要贡献

统一诊断： 证明 $M(R)$ 序列端点的概率分布是多信使框架下约束高密度 EOS 行为的敏感且互补的诊断工具。
基准独立性与依赖性： 系统地将观测数据的影响与强子基准模型（SFHo 与 DD2）的理论选择区分开来。
张力识别： 识别出硬 DD2 EOS、“质量间隙”中子星假设与严格潮汐形变约束（ $\Lambda_{190}$ ）之间的特定张力。

4. 关键结果

最大质量（ $M_{TOV}$ ）分布

观测的主导性： $M_{TOV}$ 分布的峰值主要由观测约束驱动，而非基准强子 EOS（SFHo 或 DD2）的选择。
约束的演进：
- 仅 pQCD： 产生略高于 $1.5 M_\odot$ 的峰值。
- GW170817/HESS： 峰值位置无显著偏移。
- NICER： 将峰值移至 2.1–2.2 $M_\odot$ 并锐化分布。
- 质量间隙（GW190814）： 将峰值显著移至 2.6–2.7 $M_\odot$ 。
稳健结论： 当应用最稳健的约束（pQCD + GW170817 + NICER）时，两种 EOS 模型均收敛于 2.2–2.3 $M_\odot$ 的优选最大质量。这与观测到的 PSR J0952–0607 的质量（ $2.35 \pm 0.17 M_\odot$ ）一致。

半径（ $R_{TOV}$ ）分布

基准依赖性： 与质量不同，半径分布显示出对基础强子 EOS 的强烈依赖性。
- SFHo+eLSM： 在稳健约束下，峰值位于 11.8 km 附近。
- DD2+eLSM： 在稳健约束下，峰值位于 12.4 km 附近。
总体偏好： 组合约束倾向于 $12 \pm 1$ km 的半径范围。
潮汐形变的影响： 应用更严格的 $\Lambda_{190}$ 约束会收窄半径分布并降低优选值，因为较低的潮汐形变与较小的半径相关。

“质量间隙”张力

如果将 GW190814 天体解释为中子星，则 基于 DD2 的 EOS 在与 $\Lambda_{190}$ 潮汐约束结合时无法提供一致的描述。DD2 模型中支撑 $2.6 M_\odot$ 恒星所需的刚度会导致潮汐形变过大，无法满足 GW170817 的界限。
基于 SFHo 的 EOS（较软）处理这种组合的情况稍好，但仍面临挑战。

5. 意义

这项工作对中子星最大质量进行了严格的统计评估，超越了单点估计，转向概率分布。

当前模型的验证： 结果支持存在高达 $\sim 2.2-2.3 M_\odot$ 的中子星，验证了当前核物理模型与天体物理数据的一致性。
对相变的约束： 平滑匹配程序及由此产生的分布表明，在当前观测所探测的密度范围内，不太可能发生大的一阶相变，或者至少受到严格约束。
未来方向： 研究强调，未来半径测量的精度（例如来自 NICER 或未来的 X 射线任务）以及 GW190814 天体的确切性质（中子星还是黑洞）对于解决高密度 EOS 中剩余的不确定性至关重要。质量间隙假设与潮汐形变约束之间的张力表明，如果 GW190814 确实是中子星，那么 EOS 必须比 DD2 模型预测的更软，或者潮汐约束需要重新评估。