研究背景：脉冲星“时钟”的突然加速
脉冲星被誉为宇宙中最精准的时钟。然而，这些高度稳定的“时钟”偶尔会突然加速，表现为自转频率在短时间内骤然增加，这一现象被称为周期跃变（glitch）。目前，约6%的已知脉冲星曾表现出周期跃变，其中Vela脉冲星跃变频繁、幅度较大，是研究这一现象最重要的观测对象之一。
目前普遍认为，周期跃变源于中子星内部超流成分与可观测“壳层”之间的角动量突然转移。旋转超流体以量子化涡旋（vortex）的形式承载角动量。部分涡旋可以钉扎在内壳层的原子核晶格上，使超流体与持续减速的壳层暂时脱耦。当作用于涡旋的力超过钉扎所能承受的极限时，大量涡旋可能发生集体解钉扎（unpinning），将长期积累的角动量在极短时间内传递给壳层，从而导致可观测自转频率突然增加。
2016年，Vela脉冲星发生了一次大幅度跃变。脉冲到脉冲的高时间分辨率观测为研究跃变上升阶段的瞬态动力学提供了前所未有的机会。相关观测将跃变上升时间限制在12.6秒以内，并捕捉到跃变后的短暂过冲特征，对中子星内部不同组分之间的角动量交换施加了严格约束。
为什么跃变上升阶段如此难以研究？
尽管跃变已被观测了数十年，其上升阶段仍然是中子星研究中最难解析的过程之一。
首先，跃变上升通常发生在数秒甚至更短的时间尺度内，而脉冲星计时观测通常需要累积多个脉冲以获得足够的信噪比。因此，绝大多数跃变只能观测到发生前后的频率变化，真正的上升过程往往无法被直接分辨。只有像Vela这样明亮、跃变幅度较大且被持续高时间分辨率监测的脉冲星，才有可能提供逐脉冲尺度的动力学信息。
其次，跃变涉及跨越巨大尺度范围的多层次物理过程：微观尺度上的涡旋钉扎、解钉扎、拖曳和Kelvin波激发，会共同决定宏观尺度上的角动量转移以及整颗中子星的自转响应。如何将核物理尺度上的相互作用与天文观测中的计时信号定量联系起来，是一项极具挑战性的任务。
此外，中子星内部包含壳层超流、核心超流、原子核晶格以及带电组分等多个相互耦合的成分。完整描述这些成分需要求解复杂的三维多流体动力学问题，其中还可能涉及涡旋弯曲、缠结、蠕动及跨越壳层—核心界面的运动。现阶段的观测数据也往往只能约束若干微观参数的组合，而难以对每一种机制进行完全独立的测量。
因此，跃变上升阶段既是最难观测和建模的过程之一，也是连接中子星微观物理与宏观观测最直接的窗口之一。本研究聚焦于以下问题：钉扎能、夹带效应和互摩擦机制等微观输入如何影响跃变上升过程？现有观测能够在多大程度上约束这些参数及其组合？这些动力学信息能否进一步用于检验致密物质状态方程？
核心创新：从微观物理出发建立跃变动力学模型
本研究的核心创新在于建立了一套以微观物理参数为输入、自洽计算跃变动力学并与高时间分辨率观测进行统计比较的方法。
研究团队采用基于相对论平均场理论构建的统一状态方程，并选取DD-ME2和PKDD两种有效相互作用，自洽计算中子星从壳层到核心的内部结构、物质成分、转动惯量以及内壳层涡旋钉扎性质。在此基础上，研究通过柱坐标投影，将复杂的三维星体近似简化为一维动力学问题，构建了由壳层超流、核心超流和“壳层”组分组成的三组分跃变模型。其中，“壳层”组分包括固态壳层以及与其紧密耦合、共同参与可观测自转的其他成分。
这一处理采用了轴对称以及涡旋保持笔直、刚性的近似。研究强调，刚性涡旋是假设在圆柱模型中常用的一种几何简化，其作用是使复杂的多流体动力学问题在计算上可处理，并不意味着它能够完全真实地描述中子星内部的涡旋构型。因此，在解释模型给出的定量约束时，需要充分考虑这一近似的适用范围。
在该框架中，涡旋运动产生的互摩擦（mutual friction）是控制不同组分之间角动量转移速率的关键物理机制（见文章表3）。壳层超流的互摩擦系数 B_cr 主要由涡旋与原子核相互作用所激发的Kelvin波决定；核心超流的互摩擦系数 B_co 则主要源于电子对磁化涡旋的散射。
研究还将夹带效应（entrainment）纳入动力学模型。夹带效应是指部分自由中子的运动受到原子核晶格影响，从而降低能够独立传递角动量的有效超流密度。这一效应直接影响角动量库的大小，因而是解释跃变幅度和上升过程的重要微观因素。
论文表3 中子星内部摩擦机制的比较
关键发现：跃变上升形态携带微观物理信息
通过系统的动力学模拟和贝叶斯参数推断，研究获得了以下主要结果：
壳层超流耦合时标约100秒。 跃变触发后，壳层超流与可观测壳层在约100秒的时间尺度上逐步达到耦合。推断得到的壳层超流互摩擦系数 B_cr 约为10⁻³—10⁻²。模型显示，其中，赤道区域的耦合作用最强，而弱耦合区则位于赤道面纯壳层区域的中部。
核心超流呈现“过冲”行为。 核心中央区域较强的互摩擦使该区域的超流体在跃变后迅速与壳层耦合，其角速度可以短暂超过最终平衡值，形成动力学过冲（overshoot）。随后，角动量向核心外层重新分布，系统最终趋于整体平衡。这一结果表明，跃变后的过冲不仅是计时曲线中的瞬态特征，也可能反映中子星内部不同区域具有不均匀的耦合强度和响应时标。
观测到的过冲要求壳层摩擦强而核心摩擦弱。 模型成功再现了2016年Vela跃变中观测到的壳层自转频率过冲。研究发现，在该模型框架下，产生这一特征要求壳层互摩擦整体上强于核心互摩擦。贝叶斯分析推断，核心互摩擦系数 B_co 的数量级约为10⁻⁴，与电子散射机制给出的理论预期相符。
内壳层夹带效应较弱。 分析结果表明，内壳层夹带效应可能相对较弱，约70%的自由中子仍可保持有效超流状态。推断得到的夹带参数 m_n*/m_n 约为1.76—6.13，与近期相关微观理论计算的结果总体一致。
关键参数间的相关性。 贝叶斯分析揭示了跃变前减速持续时间 Δt 与减速幅度 Δν_slow 之间的内在反相关关系，以及Vela脉冲星质量 M 与夹带强度 m_n*/m_n 之间的反相关关系。
研究还发现，部分跃变动力学参数会随状态方程的核物质参数呈现系统性变化。这些结果表明，即使在经过简化的动力学框架下，跃变上升曲线及其过冲形态仍然对特定的微观物理参数组合具有敏感性，并可能为约束中子星状态方程提供新的观测途径。
论文图7｜中子夹带效应（entrainment effects）对周期跃变动力学的系统性抑制。
(a) 不同中子有效质量比 mₙ*/mₙ 下，脉冲星壳层频率变化 Δν 随时间的演化。随着夹带效应增强，即 mₙ*/mₙ 增大，频率过冲（overshoot）的峰值显著降低，系统最终达到的平衡频率增量也随之减小。(b) 相应的计时残差（timing residual）演化。不同颜色曲线代表不同的 mₙ*/mₙ，结果表明计时残差对夹带强度具有较高敏感性。因此，对自转突变后早期演化过程进行高时间分辨率观测，有望为约束中子星内部的夹带效应提供重要信息。其中，mₙ* 表示夹带效应修正后的中子有效质量（neutron effective mass），mₙ 表示裸中子质量（bare neutron mass）；mₙ*/mₙ 越大，代表夹带效应越强。
科学意义与未来展望
本文从理论方法、状态方程框架与未来观测三个层面，阐述其科学意义与展望：
以往的跃变研究通常采用唯象耦合时标描述不同组分之间的角动量交换。本研究进一步将这些耦合时标与涡旋钉扎、互摩擦、夹带效应以及状态方程联系起来，尝试建立从核物理尺度到天文观测尺度的定量桥梁。研究结果显示，在特定的动力学假设下，观测到的跃变上升和过冲特征确实对若干微观物理参数组合具有敏感性。这为未来发展更真实的三维多流体动力学模拟、改进涡旋钉扎与蠕动（creep）的处理，以及自洽纳入壳层和核心夹带效应奠定了基础。
（二）建立统一的状态方程框架
研究采用同一种有效相互作用，对中子星从壳层到核心的物质性质进行统一描述，使钉扎能、互摩擦系数、夹带效应和转动惯量等关键物理量能够在相互一致的理论框架下计算。这有助于减少不同微观输入彼此不相容所造成的系统误差，并为比较不同状态方程的动力学效应提供统一基础。
（三）为未来高精度观测做好理论准备随着中国天眼FAST、平方公里阵列SKA和增强型X射线时变与偏振空间天文台eXTP等观测设施的发展，更高时间分辨率和更高精度的脉冲星计时观测将逐步成为可能。未来观测有望以秒级乃至更高时间分辨率解析跃变上升和瞬态过冲，为检验更真实的涡旋动力学及三维多流体模型提供关键数据。如果能够同时获得更多脉冲星的跃变上升时间、过冲幅度和恢复过程，跃变动力学有望从解释单次事件的理论工具，进一步发展为探测中子星内部超流性质和致密物质状态方程的新型“微观物理探针”。
该工作由厦门大学天文学系涂中豪博士后（第一作者）与李昂教授（通讯作者）合作完成。研究获得了国家自然科学基金等项目的支持。
本研究是团队在致密物态与脉冲星跃变交叉方向上的系列工作之一。从统一状态方程建模和微观钉扎力计算，到跃变动力学的自洽模拟与贝叶斯推断，研究团队持续致力于将微观核物理与宏观天体物理观测相结合，为理解中子星内部结构及致密物质性质提供理论支持。欢迎关注相关研究工作（公众号：Nucl_Astrophys_xmu）。
论文信息：
Zhonghao Tu, Ang Li. A Microphysical Probe of Neutron Star Interiors: Constraining the Equation of State with Glitch Dynamics, 2026 ApJ 1007, 118
论文链接：https://iopscience.iop.org/article/10.3847/1538-4357/ae8768
预印本链接：https://arxiv.org/abs/2510.20791
