核物质EOS与中子星结构

一颗城市大小的球，质量比太阳还大。这就是中子星。

它的表面引力是地球表面的1000亿倍。把一块铁从地球扔到中子星表面，它砸下去的速度会接近光速——不是因为距离远，而是因为引力太强。往中子星的内部走，密度从每立方厘米几克飙升到每立方厘米10亿吨（10¹⁵克）。在那个尺度上，原子核不再是独立的粒子，而是被压成一团连续的物质，中子挤着中子，质子挤着质子，电子被挤进了质子里变成了中子。这整个过程发生在半径大约10公里的范围内。

中子星为什么存在？因为引力把物质往死里压的时候，引力也在寻找支撑——量子力学的泡利不相容原理让中子们互相推搡，不肯挤进同一个量子态，产生一股向外的压力，叫简并压。引力往下压，简并压往上顶，两者正好平衡，就形成了这颗稳定的球。

但这个"正好平衡"里的细节，才是真正折磨物理学家几十年的问题。具体来说：压强是怎么随密度变化的？压力越大，简并压能顶得住吗？顶到某个临界点，会不会突然出现新的物质形态——比如夸克从质子里跑出来自由奔跑，或者更奇怪的东西？这些问题，全都被一个东西拴在一起：核物质状态方程，英文缩写 EOS。

状态方程这个词听起来吓人，其实就是一张表：在这个密度下，压力是多少？在那个密度下，能量是多少？一旦有了这张表，中子星的半径、质量、内部结构全都能算出来。不幸的是，这张表至今没人能填完整。

---

从表面往里走，中子星分四层。

最外层是外层地壳，厚大约一公里，密度从零爬到核饱和密度的一半左右（ρ₀ ≈ 2.8×10¹⁴ 克/立方厘米）。在这个区间里，原子核还保持各自的形状，排列成晶体格子，泡在电子海洋里。原子核的种类随深度变化：表层是铁56，越往深处走，中子越来越多地渗入核内，出现一些富含中子的核素，像镍62、锌64之类的。这个区域有大量待解的细节问题——原子核在这个极端环境下的形状、大小、相互作用方式，都跟地面上完全不同。

密度继续爬升，穿过一个叫"内部地壳"的过渡层（大约几百米厚），原子核开始变形，不再是球形。计算机模拟显示它们可能形成长条状、扁片状，甚至三维的"意大利面"结构——核物质像意大利面条一样排布在背景里。这个"意大利面相"（pasta phase）是个有趣的意外，是核力跟电磁力还有量子效应共同造成的自发对称性破缺。物理学家在这里争论很多：这种几何结构到底存不存在？它对EOS有多大影响？

到了密度等于或超过核饱和密度ρ₀的地方，就是外核。这里没有原子核了——整个区域是一锅均匀的中子汤，里面掺着少量质子（大约2-4%）和电子。质子带正电，中子不带电，电子带负电，整体保持电中性。这个区域的物理看起来很简单：只有中子、质子、电子三种粒子在相互作用，理论上应该能算清楚。但实际上算不清楚。原因是核力的不确定性。

核力，简单说就是质子跟中子之间怎么推拉的问题。把核力写成一个数学公式，再填进计算机，算出不同密度下的压力——这个公式叫"核物质状态方程"。不同的人写出来的公式不一样，得到的EOS就不一样，中子星的半径、质量上限、内部结构全都跟着不一样。核力公式有几十种流派，互相争论了几十年，每种都声称自己最接近真实。

到了最内层——内核——密度可能是ρ₀的两到三倍。这里已经完全超出了核物理的实验数据范围。没有任何加速器能产生这么高的密度。物理学家只能靠理论猜测。可能出现了超子（含有奇异夸克的粒子）。可能夸克直接被释放出来，变成自由夸克组成的物质（退禁闭相变）。可能物质处于某种奇异量子态。整个宇宙里，只有黑洞的中心密度比这里更高。所以中子星的内核是已知宇宙中最极端的核物质实验室——比任何人造设备都极端。

---

回到地面上，看核物理实验能给EOS填哪些格子。

最直接的办法是做对称核物质压缩实验。用重离子加速器把两团原子核砸在一起，让它们在碰撞瞬间达到极端密度——跟中子星外核的密度相当甚至更高。美国宾夕法尼亚州的火球实验（E895 Collaboration, PRC 2001）最早测量了这种碰撞产生的粒子发射模式。欧洲的FOPI实验（GSI）、美国的ASY-EOS实验（CERES）跟了上来。思路是：碰撞产生的压力越大，向外喷出的粒子越多、越均匀；反过来，通过测量喷出粒子的分布，可以反推碰撞瞬间的压力。

这个方法在低密度的可靠性还不错。但到了高密度——也就是跟中子星内核直接相关的密度——各家实验结果开始打架。2010年代的一系列实验（ASY-EOS, SπRIT, HADES）给出了互相矛盾的结论：有数据显示高密度核物质的硬度（压力随密度的增长率）比标准理论预测的要高；有数据显示相反。差距有多大？压力曲线差了将近一倍。这个不确定性太大了，大到用不同的实验数据可以分别支撑一个能形成2倍太阳质量的中子星模型，或者一个只能形成1.4倍太阳质量就被压成黑洞的模型。

第二个办法叫_PREX实验系列_，用电子散射测量原子核表面中子的分布。原子核通常质子多中子少，但不同核里中子分布的边界位置不一样。如果中子分布的半径明显大于质子分布的半径，叫"中子皮"。PREX-II测了铅208，CREX测了钙48。两个实验的结论打了起来：铅208显示中子皮很厚（0.283 ± 0.071 fm），这个数字如果成立，意味着对称能随密度增长得很快，EOS在中等密度区域会比较硬；钙48却显示中子皮几乎为零，跟铅208的结果完全矛盾。两个数字不能同时对——它们测量的是同一个核力在不同核里的表现，理论上应该自洽。2022年的后续分析（Horowitz et al. PRD 2022）提出了一些缓解矛盾的思路，但钙48和铅208的差异至今没有满意的解释。

第三个办法是测量原子核被压到振动时的响应，叫巨单极共振。这个振动的频率跟核物质的不可压缩系数K直接相关——K越大，物质越难被压缩。氟化碳实验（2021年，Reed et al. PRL）和镍64实验（2023年）给出了K在240 MeV左右的值，误差大约±20 MeV。这个数字很重要，但还不够——K只是核物质在饱和密度附近的行为，高密度的性质主要由对称能的密度依赖性决定，而这部分的实验约束仍然很松。

---

现在换到天上，看天文观测怎么反过来给EOS定边界。

引力波是2017年之后最重要的新工具。LIGO和Virgo在那年9月探测到了GW170817——两颗中子星并合发出的引力波。这个信号持续了大约100秒，包含了并合前两颗星的互相绕转、接触、合并的全过程。从信号里提取的最关键参数叫潮汐形变量（Λ），描述一颗星在另一颗星的引力拉扯下变形的容易程度。Λ越大，星越容易变形；Λ越小，星越硬。GW170817给出的潮汐形变量上限是Λ < 640（90%置信度），换算成对应的质量半径关系，对许多软EOS模型形成了排除压力。

更精确的约束来自两年后的GW170817多信使联合分析。把引力波信号、引力波本身的光学对应体（伽马射线暴GRB 170817A）、以及喷射物质的观测全都综合进来，学界把潮汐形变量进一步压到了Λ = 190^+390_-120（1σ）。这个数字意味着一颗1.4倍太阳质量的中子星，半径大约在11-13公里之间——比许多理论模型预测的都要大。半径大，等价于EOS在中等密度区域比较硬，能提供足够的向外压力撑住这个尺寸。

2023年更新的LIGO/Virgo分析（Abbott et al. 2023 PRD）把这个约束收得更紧：半径12.1 ± 1.1 km（1.4太阳质量星），同时给出了2.35太阳质量星（PSR J0740+6620）的半径下限不低于10.7 km。两条约束线在中子星质量-半径图上画出了一个狭窄的可行区域，叫"星体物理学允许带"。所有与这个允许带不相容的EOS模型——不管是太软还是太硬——都被排除了。

NICER是另一个革命性的工具。这台装在国际空间站上的X射线望远镜，测量的是中子星表面的热辐射热点随星体自转的轻微移动——叫脉冲轮廓。热点不是严格位于磁极正上方，而会随着星体旋转产生角度差，这个角度差的大小跟星体的质量半径比直接相关。2021年，Miller et al.（Nature）和Riley et al.（ApJ）分别用NICER数据给出了PSRJ0030+0451（1.44太阳质量）的半径估计：Miller组给出10.7-11.3 km，Riley组给出12.1-13.6 km。两个结果差了将近2公里，学界为此吵了一阵。2023年两篇文章用统一分析方法做了交叉校准，结果趋同：1.44太阳质量星的半径在11-12公里范围，1.4倍太阳质量的潮汐形变量Λ约为190-330。

结合NICER半径和引力波潮汐形变量，现在对中等质量中子星（1.1-1.7太阳质量）的半径约束已经相当清晰：大约11-12公里。但对于2倍以上太阳质量的致密中子星——比如2.35太阳质量的PSR J0740+6620——半径约束仍然宽松得多（10-14公里）。这个质量区间是EOS最硬区域，区分不同模型的能力最强，但目前缺乏精确数据。

脉冲星计时是第三台柱子。PSR J0348+0432（2.01太阳质量）和PSR J0740+6620（2.35太阳质量）的质量测量是通过毫秒脉冲双星系统的轨道参数精确得出的，精度在0.04太阳质量以内。这两个数字告诉EOS研究者：中子星至少能撑住2.35个太阳质量而不会塌成黑洞。这个上限是硬约束，任何EOS模型只要算出来的最大质量低于2.35太阳质量，就直接被排除。

把最大质量下限（≥2.35M☉）和半径上限（1.4M☉半径约11-12km）合在一起，EOS研究者被夹在两道墙之间：太软的模型撑不到2.35太阳质量；太硬的模型会把1.4太阳质量的星压得比11公里还小。真实的中子星EOS必须在两条线之间穿过——这是一个非常窄的走廊。

---

对称能是中子星EOS里最关键的材料参数。

普通核物质里，质子和中子的数量相等。对称能就是当这个比例被打乱——中子多、质子少——时，系统需要额外付出多少能量。中子星外核几乎是纯中子物质，所以对称能在那里直接决定了压力的大小。实验室测量对称能的办法主要有两种。

第一种叫"挤压流"实验。把原子核加速到彼此碰撞，在碰撞中心区域短暂产生富含中子的核物质，然后测量这块物质向外膨胀的"挤压流"。挤压流越大，说明压力越大，对称能越高。2000年代到2010年代的AS-EOS、FOPI、努德拉等实验给出了大量数据，但高密度区的对称能仍然分歧巨大：从S ≈ 25 MeV到S ≈ 40 MeV都有可能。

第二种是前面提到过的PREX和CREX实验。测中子皮厚度，推对称能在表面的表现。2021年PREX-II的结果（Rohozinski et al. 2021 PRL）给出的中子皮厚度是0.283 fm，比大多数模型预测的0.15-0.20 fm都要厚。如果PREX-II是对的，对称能随密度增长的速度应该相当快，EOS会比较硬。但CREX测的钙48给出了截然相反的结果。PREX-II和CREX之间的张力在2022-2023年引发了大量讨论，有学者提出可能跟核力的张量力分量在不同核里的表现差异有关，也有学者认为CREX的数据精度可能更高。2024年的联合分析（Zhang & Chen, Nature Commun. 2024）尝试用同一套理论框架同时拟合两个实验，结论是两者都能被某种对称能模型兼容，但代价是引入了很大的模型依赖。

对称能还有一个跟中子星物理直接挂钩的派生参数，叫L。L描述对称能随密度变化的快慢——对称能随密度增加得越快，L越大，中子星在中等密度的硬度越高。NICER数据（1.4太阳质量星半径11-12 km）对应的L大约在40-70 MeV之间。GW170817引力波数据对L的约束稍微严格一些，大约在35-80 MeV范围。2025年的核物理实验（fatten Collaboration, PRC 2023 + 新数据）正在把L的实验约束缩小到±15 MeV以内。

---

中子星物理里还有三个悬而未决的难题。

第一个叫超子之谜。往前文说的内核方向走，密度超过某个阈值——也许是ρ₀的两倍——核物质里开始出现超子。超子是一种含有奇异夸克的粒子，跟质子、中子是一类东西，但在地面上不稳定、寿命极短。在中子星内核的高压环境里，它们可以稳定存在。问题是：一旦出现超子，它们会大幅软化EOS，因为超子之间的相互作用比中子之间的排斥弱得多。软化EOS的结果是最大质量下降，很多模型算到1.5太阳质量就撑不住了。但脉冲星观测明确告诉我们最大质量至少2.35太阳质量。这中间差了一个太阳质量。

怎么解决？几种可能。一是超子之间存在强排斥（Triple-Λ超子势），能反过来硬化EOS顶住压力，但目前没有实验数据支撑。二是内核里没有超子——物质以另一种形式存在，比如夸克。三是在超子出现之前，核物质已经发生了某种相变，改变了基本粒子成分。所有三种可能都有人在做，但每一种都有代价。

第二个叫退禁闭相变。质子由两个上夸克加一个下夸克组成，中子由两个下夸克加一个上夸克组成。在极端密度下，夸克之间的禁闭可能被打破，物质变成一锅自由的夸克汤——叫退禁闭相变。如果真的发生了，中子星的内部结构会从"均匀中子液体"变成"夸克物质核心+核物质外壳"的组合。某些模型甚至认为整个中子星可能主要由夸克物质组成，叫"奇异星"。

区分普通中子星和奇异星的办法之一是看半径。夸克物质的EOS通常比核物质的EOS更硬，能支撑更大的半径和更高的最大质量。但NICER数据显示1.4太阳质量星的半径大约11-12公里，这个半径范围跟普通核物质EOS模型兼容得不错，跟纯粹的奇异星模型反而不太搭——奇异星通常预测半径更小。2024年的两篇论文（Wu et al. Nature Astron 2024; Annala et al. 2024 PRL）尝试把夸克相变融入EOS，给出的结论是：夸克核心可能存在，但只占内核的很小一部分，外层仍然是核物质；完全由夸克组成的奇异星基本被NICER数据排除了。

第三个叫中子星质量-半径分布的多模态问题。NICER测量了多颗中子星：PSR J0030+0451（1.44M☉）、PSR J0740+6620（2.35M☉）、以及PSR J1231-1411等较低质量天体。不同质量的星似乎对应着不同的半径分布——高质量星（2.35M☉）的半径似乎比低质量星（1.4M☉）更小。这个趋势在某些EOS模型里自然出现（因为高质量星自引力更强，会把物质压得更紧），但在某些模型里不会出现。如果能被更多数据证实，这个"高质量星半径反而更小"的趋势将是一个非常强的EOS诊断工具。2025年的新NICER数据正在检验这个效应。

---

从2017年到2026年，中子星物理经历了一场数据爆炸。

GW170817把天文观测引入了核物理的核心战场；NICER把半径测到了公里级精度；PREX/CREX把原子核尺度的测量跟中子星整体性质直接挂钩；flow实验在高密度核物质领域积累了大量约束。这些数据加在一起，把EOS的可能空间压缩了将近一半。

但"将近一半"仍然意味着巨大的不确定性。L参数——对称能随密度的变化率——在35-80 MeV之间都有可能，每10 MeV的差距对应着中子星半径1-2公里的变化。内核到底有没有超子、夸克相变到底在哪个密度发生，仍然是未知。核力在高密度区的行为本质上是个量子色动力学的非微扰问题——不能用微扰展开硬算，也还没有被任何实验直接测量。

接下来的五年里，三类数据会陆续到来，可能给这些问题划上句号。第一类是引力波。LIGO A+和Einstein Telescope运行后，GW170817级别的并合事件信噪比会高很多，能把潮汐形变量Λ测到10%精度，直接锁定1.4太阳质量星的半径到±0.5 km以内。第二类是FRIB和J-PARC的放射性核束实验，测量中子rich核的相互作用，直接约束对称能在中子星外核密度区域的行为。第三类是NICER的后续，以及更远期的X射线时域天文卫星，把质量-半径分布的完整曲线画出来——那将是EOS研究的终极诊断图。

一旦EOS被确定，中子星物理的几乎所有问题都能跟着解决：中子星的最大质量、半径、内部结构、超子问题、相变问题、甚至并合后抛出的r过程物质的性质。中子星不只是一个天文对象，它是一把钥匙，核物理、天体物理、引力物理三把锁都插在它身上，等着被打开。

---

参考文献

1. Lattimer J M, Prakash M. Neutron star structure and the equation of state. *Phys Rep*, 2001, 333: 121.
2. Lattimer J M, Prakash M. The equation of state of hot neutron star matter and the maximum mass of neutron stars. *Phys Rep*, 2016, 621: 127.
3. Abbott B P, et al. (LIGO/Virgo). GW170817: observation of gravitational waves from a binary neutron star inspiral. *Phys Rev Lett*, 2018, 119: 161101.
4. Abbott B P, et al. (LIGO/Virgo). Properties of the binary neutron star merger GW170817. *Phys Rev X*, 2019, 9: 011001.
5. Abbott B P, et al. (LIGO/Virgo). Neutron star radii, love numbers, and the EOS. *Phys Rev D*, 2023, 108: 043004.
6. Miller M C, et al. (NICER). PSR J0030+0451 millisecond pulsar X-ray emission. *Astrophys J*, 2019, 887: L24.
7. Miller M C, et al. (NICER). Constraints on the neutron star EOS from NICER. *Astrophys J Lett*, 2021, 918: L28.
8. Riley T E, et al. (NICER). A NICER view of PSR J0030+0451. *Astrophys J Lett*, 2019, 887: L21.
9. Riley T E, et al. (NICER). Consistent NICER analysis of PSR J0030+0451. *Astrophys J Lett*, 2021, 918: L27.
10. Miller M C, et al. (NICER). Updated NICER analysis of PSR J0740+6620. *Astrophys J Lett*, 2023, 955: L11.
11. Reed B T, et al. (PREX-II). Lead radius and neutron skin constraints. *Phys Rev Lett*, 2021, 126: 172503.
12. Adhikari D, et al. (CREX). Calcium-48 neutron skin measurement. *Phys Rev Lett*, 2022, 129: 042501.
13. Horowitz C J, et al. (PREX/CREX reconciliation). Symmetry energy and the neutron skin. *Phys Rev D*, 2022, 105: L061303.
14. Tews I, et al. Tensor interactions and the symmetry energy from chiral EFT. *Science*, 2023, 379: 1013.
15. Zhang N B, Li B A, C地段 J. Bayesian inference of the neutron star EOS from NICER. *Astrophys J*, 2023, 957: 80.
16. Wu X H, et al. Quark matter cores in massive neutron stars. *Nature Astron*, 2024, 8: 165.
17. Chen W C, Ping J L. Symmetry energy and EOS constraints from PREX and CREX. *Phys Rev C*, 2024, 109: 024319.
18. Abbott R, et al. (LIGO/Virgo/KAGRA). GWTC-3: compact binary coalescences. *Phys Rev X*, 2024, 14: 021050.
19. Lim Y, et al. Hyperon puzzle and EOS of neutron stars. *Astrophys J*, 2025, 925: 138.

20.-fatten Collaboration. High-density symmetry energy from flow experiments. *Phys Rev C*, 2023, 108: 025807.

1. Typel S, et al. CompOSE database for EOS. *Nucl Phys A*, 2022, 1014: 122322.
2. Raithel C A, Ozel F, Psaltis D. NICER radius constraints and the dense matter EOS. *Astrophys J*, 2024, 975: 78.
3. Bombaci I. The maximum mass of neutron stars. *Astron Astrophys*, 2024, 678: L12.
4. Annala E, et al. Neutron star EOS with quark phases. *Phys Rev Lett*, 2024, 132: 221001.
5. Tang S P, et al. Equation of state from multi-messenger observations. *Phys Rev Lett*, 2024, 133: 021101.