$^1$S$_0$ pairing gaps, chemical potential and entrainment matrix in superfluid neutron-star cores for the Brussels-Montreal functionals

想象一下宇宙中最极端的游乐场：中子星。这些是死掉的恒星坍缩后形成的、密度极高的核心，其质量之大，以至于在地球上，仅仅一茶匙这种物质的重量就高达十亿吨。在深处，压力如此巨大，以至于原子被压扁，只留下由亚原子粒子组成的“汤”。但这不仅仅是一个静态、冰冷的肿块；它是一个动态的、翻腾的流体。事实上，其核心被预期是一种“超流体”，这是一种粒子流动时摩擦力为零的状态，就像一种永不停歇的液体。

为了理解这些恒星是如何旋转、摇摆甚至发生“突跳”（突然加速）的，科学家需要了解内部不同粒子是如何相互作用的。把核心想象成一个拥挤的舞池，有两个主要的群体：中子和质子。在超流体中，这些粒子会配对并保持完美的同步起舞。然而，它们并非独自起舞；它们会互相影响。如果中子开始旋转，它们就会带动质子一起旋转，反之亦然。这种“拖拽”被称为纠缠（entrainment）。这就像两个手拉手的滑冰者；如果其中一个试图转身，即使另一个并不想转，也必须跟着转。你即将阅读的论文将深入探讨描述这种抓握力究竟有多强的数学细节，包括温度如何改变这场舞蹈，以及粒子在配对破碎前能移动多快。

---

宇宙舞池：中子星的新地图

中子星是自然界终极的实验室，但它们也非常难以建模。中子星的核心是由中子和质子（连同一些电子和缪子）组成的致密混合物，其行为类似于超流体。在这种状态下，粒子形成配对并无阻力地流动。但问题在于：中子和质子并不只是独立流动。由于它们之间强大的核力作用，它们是“纠缠”在一起的。这意味着如果你试图推动中子，质子就会被拖动；如果你推动质子，中子也会随之跟随。这是宇宙版的“一团乱麻的湿意面”，移动其中一根，就会拉动整个乱团。

长期以来，科学家一直试图绘制出这种拖拽作用的具体图谱。他们需要了解三个主要要素：打破粒子配对需要多少能量（配对能隙/pairing gap）、维持它们运动需要多少“推力”（化学势/chemical potential），以及中子和质子结合得究竟有多紧密（纠缠矩阵/entrainment matrix）。问题在于，这些数值会随着恒星的热度和粒子的流动速度而变化。以往的模型通常不得不进行猜测，或者使用无法同时兼顾所有这些变化条件的简化规则。

这篇由 Valentin Allard 和 Nicolas Chamel 撰写的论文通过不折不扣地求解完整的复杂数学方程（即随时间变化的 Hartree–Fock–Bogoliubov 方程），实现了巨大的跨越。他们使用了一套关于核物质行为的特定规则，称为 Brussels–Montreal functional BSk24，该规则以高度精确著称。他们模拟了中子星外壳核心中超流性可以存在的整个温度和速度范围。

他们的发现：舞蹈的规则

该团队的主要发现是关于中子和质子在这种极端环境下行为的完整且自洽的图谱。他们并没有靠猜测，而是计算出了配对能隙、化学势和纠缠矩阵在整个范围内的精确数值。

他们最有趣的发现之一涉及配对能隙。这是将粒子对连接在一起的能量“胶水”。他们发现，对于质子来说，这种胶水出奇地弱——比许多其他研究假设的要弱得多。事实上，他们计算出的质子配对能隙明显小于中子能隙。这一点很重要，因为如果胶水太弱，质子可能会比之前认为的在更低的温度或速度下就停止超流行为。他们的结果与近期一些考虑了粒子在介质中如何相互作用的复杂计算相一致，表明质子的“弱胶水”特性很可能是真实的物理情况。

他们还发现，这些粒子的行为具有惊人的普适性。当你使用正确的“尺子”（通过将温度和速度归一化为它们的临界极限）来观察数据时，中子和质子的行为几乎是完全一样的。无论是中子还是质子，配对能隙随着恒星变热或流动变快而缩小的过程都遵循相同的模式。这表明底层的物理机制是相同的，受相同的弱耦合规则支配。

另一个关键结果涉及纠缠矩阵，它告诉我们中子在多大程度上拖拽质子。作者发现，只要流动速度低于某个特定的临界极限（朗道临界速度），这种拖拽效应就不会发生太大变化，即使粒子移动得很快。然而，一旦流动过快或温度过高，超流性就会开始崩溃。在那时，中子和质子不再互相拖拽，而是变得“动力学解耦”。中子可能会完全失去超流性，而质子则继续独自起舞，不再受中子的束缚。

他们排除了什么以及仍有哪些未知

论文非常明确地说明了其未涵盖的内容。作者明确指出，他们仅研究中子星的外壳（outer core），那里的密度很高，但还不足以产生像超子或夸克这样的奇异粒子。他们也仅关注了 1S0 配对相（一种中子和质子配对的特定方式）。他们排除了 3PF2 中子超流性，指出在两者可能共存的区域，除非存在极其强大的磁场，否则 1S0 相会完全将 3PF2 相挤出局外。因此，对于外壳而言，3PF2 相实际上是不存在的。

他们也并未仅仅依赖假设；他们针对各种近似模型进行了测试。他们发现，虽然一些较旧、较简单的模型（如朗道的理论）在已知临界温度的情况下对于纠缠矩阵是有效的，但不足以准确计算化学势。要获得真实的化学势，必须求解他们所使用的完整、高强度的方程。

为什么这很重要

这项工作为构建更好的中子星模型提供了“原料”。就像厨师需要精确测量面粉和糖的用量才能烤出完美的蛋糕一样，天体物理学家需要精确的配对能隙和纠缠数值，才能理解中子星如何旋转、如何冷却以及如何对引力波做出反应。通过提供基于 BSk24 泛函的、一致的微观输入，Allard 和 Chamel 为模拟这些神秘、致密恒星心脏的科学界提供了一个更可靠的“食谱”。他们并没有解决中子星的所有奥秘，但他们澄清了关于核心舞池中两位主要舞者如何共同起舞的迷雾。