⚛️ nuclear theory

Bayesian approach study of hybrid neutron stars

本研究通过使用贝叶斯分析将用于强子物质的相对论平均场模型与用于夸克物质的矢量-MIT袋模型相结合，证明了天文质量与半径观测结果支持在密度低于 $2.0\,n_{0}$ 时发生向解禁闭夸克物质的相变，这意味着存在具有大夸克核的混合中子星。

原作者： Fábio Köpp, César H. Lenzi, César V. Flores, and Débora P. Menezes

发布于 2026-09-15

📖 1 分钟阅读🧠 深度阅读

原作者： Fábio Köpp, César H. Lenzi, César V. Flores, and Débora P. Menezes

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在遥远的宇宙深处，在我们的望远镜无法触及的地方，存在着一些宇宙中最极端的物体：中子星。这些是质量巨大的恒星坍缩后的核心，被压缩成直径仅约二十公里的球体，却包含了比我们整个太阳还要大的质量。在它们内部，引力将物质挤压得极其紧密，以至于原子本身都被压碎，留下了由质子和中子组成的致密“汤”。几十年来，物理学家一直在思考这种压力是否会进一步推进，迫使质子和中中子分解为其自身的更基本构建模块：夸克。如果这种情况发生，恒星的核心将转化为一种新的物态，创造出科学家所称的“混合星”——一种拥有由游离夸克构成的核心，并被普通核物质外壳包围的天体。确定这些奇异核心是否存在是现代天体物理学面临的重大挑战之一，因为这将揭示物质在地球上任何实验室都无法重现的条件下是如何表现的。

为了解决这个谜题，研究小组转向了一种被称为贝叶斯分析（Bayesian analysis）的强大统计方法。他们并没有猜测哪种物质模型可能是正确的，而是使用这种方法，根据先进仪器收集到的现实世界数据来测试几种不同的理论。他们专注于中子星质量和大小的测量，这些数据来源于监听恒星碰撞引起的时空涟漪的引力波探测器，以及通过计时旋转恒星脉冲光信号的X射线望远镜。研究人员构建了这些恒星的计算机模型，为每一颗恒星创建了两个截然不同的层次：一个由普通核物质组成的外壳，以及一个由去禁闭夸克物质构成的中心核心。他们测试了三种描述外壳的不同方式，每种方式都基于关于粒子相互作用的略微不同的假设，并将它们与一个包含特定参数（用于描述夸克如何相互推挤以及如何被禁闭）的夸克核心模型相匹配。

随后，研究小组对他们的模型进行了一个严格的排除过程。他们要求计算机模拟寻找特定的夸克核心设置，使得生成的混合星能够与观测到的真实中子星的质量和大小相匹配。至关重要的是，他们还检查了这些恒星是否能保持稳定，以及从普通物质到夸克物质的转变是否发生在符合物理常理的密度水平。结果指向了一个非常特定的结果。最成功的模型——它结合了某种特定的外壳描述与一套特定的夸克核心规则——表明，向夸克物质的相变可能发生在小于正常原子核密度的两倍的密度下。这一发现意味着，拥有大型夸克核心的混合星是一种与数据一致的可能构型，尽管支持这一特定模型的统计证据仍不确定。

然而，这项研究也强调了其他可能性所面临的显著挑战。当研究人员测试了在外壳中包含被称为“超子”（hyperons）的奇异粒子的模型，或者使用不同的核物质数学描述时，生成的恒星往往无法匹配观测数据，或违反已知的物理极限。特别是，许多替代模型（包括含有超子的模型）预测的恒星要么过于“软”，要么过于“硬”，或者它们暗示向夸克物质的转变会在密度过高的情况下发生，导致恒星在达到必要质量之前就会发生坍缩。分析表明，虽然超子模型作为一种理论可能性并未被完全排除，但它们违反了源自摄动量子色动力学（pQCD）的具体约束。只有极窄范围的夸克核心参数才能产生既稳定又符合近期天文观测精确测量值的恒星，尽管最好的拟合模型在某些方面仍与理论极限存在张力。

研究人员发现，最符合的模型预测向夸克物质的转变发生在约为正常核密度的1.9倍的密度处。在此处，恒星的核心将开始改变其性质，在内部创造出一个巨大的夸克物质区域。研究表明，这些混合星的半径约为12公里，并能支撑高达约2.1倍太阳质量的质量。这些数值与我们观测到的最重的中子星相一致，例如已知质量超过两个太阳质量的PSR J0740+6620，尽管研究指出，最好的拟合模型在某些图表中仍违反了特定的pQCD约束。研究结果表明，这些大质量恒星的内部很可能是一个复杂的相位混合体，极端的压力迫使物质发生了根本性的性质变化。

虽然这项研究为揭示宇宙中最致密物体的隐藏内部提供了更清晰的图景，但也凸显了我们当前知识的局限性。研究人员指出，转变发生的准确时刻在很大程度上取决于用于恒星外壳的具体模型。有些模型表明转变发生得更早，而另一些则将其推向更高的密度，但只有一种组合满足了来自天文数据的绝大多数约束，尽管它仍然违反了某些图表中的特定pQCD约束。团队还检查了恒星在受到伴星引力牵引时形状的变化，这一属性被称为潮汐可变形性（tidal deformability）。他们的结果显示，最成功的模型所产生的恒星在形变方式上，与2017年探测到的两颗中子星碰撞产生的引力波是一致的。

最终，这项工作为揭示宇宙中最致密物体的隐藏内部提供了更清晰的图景。通过将精确的天文测量与复杂的统计工具相结合，研究人员缩小了关于中子星核心可能存在的范围。他们证明了拥有大型夸克核心是这些宇宙巨兽的一种合理的特征，前提是控制普通物质与夸克物质之间转变的规律遵循特定的路径，即便对于这一特定情景的统计证据尚未定论。这并不证明每一颗中子星都有夸克核心，但它证明了这种构型与我们目前对宇宙的认知是相容的。随着来自更灵敏仪器的更多数据到来，这些模型将被进一步测试，可能从而揭示存在于存在之边缘的物质的真实本质。

技术摘要：混合中子星的贝叶斯方法研究

问题陈述
中子星的内部组成，特别是是否存在从强子物质到解禁夸克物质的相变，仍然是天体物理学中的一个开放性问题。虽然量子色动力学（QCD）是强相互作用的基础理论，但目前尚无法实现对高密度物质的从头算（ab initio）描述。因此，需要使用有效模型来描述状态方程（EOS）。一个关键的挑战是确定混合星（含有夸克核）是否能够存在，以及在何种条件下存在。本研究探讨了如何利用天文观测——特别是来自 NICER 的质量与半径测量以及来自 GW170817 的引力波数据——来约束夸克物质模型的参数以及相变的发生条件。

方法论
作者采用贝叶斯推断框架来约束混合星模型的参数。该方法包括以下组成部分：

微观物理模型：
- 强子物质： 使用相对论平均场（RMF）近似进行建模。使用了三组参数集：
  - NL3*ωρ：仅含核子（硬 EOS）。
  - EL3ωρ：仅含核子（中等硬度）。
  - EL3ωρY：包含超子（奇异重子）。
- 夸克物质： 使用矢量 MIT 包模型（vMIT）进行建模。拉格朗日量包含了排斥性矢量相互作用、包常数（ $B$ ）以及狄拉克海贡献项（ $b_4$ ）。
- 相变： 通过麦克斯韦构造（Maxwell construction）实现，假设存在一个一阶相变，其中压力和化学势在界面处连续，而重子密度发生跳变。
宏观物理计算：
- 通过求解托尔曼-奥本海默-沃尔科夫（TOV）方程来计算恒星结构（质量与半径）。
- 计算潮汐形变度（ $\Lambda$ ），以便与引力波约束进行比较。
- 使用近似普适关系（例如 $C-\Lambda$ 、 $\Upsilon-C$ 、 $\Upsilon-\Lambda$ ）来对照已建立的核物质行为进行结果交叉验证。
贝叶斯推断：
- 参数 ( $\theta$ )： 包常数 ( $B^{1/4}$ )、矢量耦合常数 ( $G_v$ ) 以及狄拉克海贡献 ( $b_4$ )。
- 先验分布： 假设 $B^{1/4} \in [150, 180]$ MeV、 $G_v \in [0.2, 0.3]$ fm $^{-2}$ 和 $b_4 \in [0, 2]$ 为均匀分布。
- 似然函数： 由以下部分构建：
  - 来自 PSR J0030+0451、PSR J0740+6620 和 GW170817 的质量与半径约束。
  - 一个相变密度约束（中心位于 $n_b = 0.275$ fm $^{-3}$ ），以倾向于产生大的夸克核。
  - 一个关于最大质量的严格下界（ $M \ge 2.006 M_\odot$ ）。
- 证据计算： 使用 dynesty 采样器执行，以计算贝叶斯证据并利用贝叶斯因子比较模型。

主要贡献与结果

参数约束： 研究确定了满足观测约束的具体 vMIT 参数集。最一致的结果是在 EL3ωρ 参数化（仅含核子）下通过两个参数（ $G_v, B^{1/4}$ ）得到的，其数值约为 $(0.22, 157.47)$ 。
- 该配置产生的相变密度为 $n_b \approx 1.9 n_0$ ，相变质量约为 $\approx 0.8 M_\odot$ 。
- 它产生的无量纲潮汐形变度 $\Lambda_{1.4} \approx 397$ ，与 pQCD 约束（ $\Lambda_{1.4} < 400$ ）一致。
- 该模型的自然对数证据为 $-4.54$。
超子与第三个参数的影响：
- 超子 (EL3ωρY)： 包含超子的模型导致的相变密度（ $\approx 2.13 n_0$ ）违反了 pQCD 约束，且产生的 $\Lambda_{1.4} \approx 488$ 与 pQCD 极限不符。
- 三参数模型： 包含狄拉克海贡献 ( $b_4$ ) 通常会导致结果违反 pQCD 约束（例如 $\Lambda_{1.4} > 400$ ），或未能满足无量纲转动惯量的普适关系。
- 硬 EOS (NL3*ωρ)： 使用较硬的 NL3*ωρ 参数化模型也无法满足 pQCD 约束，产生的相变质量在 $1.0-1.2 M_\odot$ 之间，且 $\Lambda_{1.4}$ 超过 400。
贝叶斯因子分析：
- 将两参数 EL3ωρ 模型与三参数版本进行比较，得到的对数贝叶斯因子为 $0.92$，表明对于更简单的模型而言，证据是“不确定的”。
- 然而，当将贝叶斯因子与物理约束（pQCD 和普适关系）结合起来时，两参数 EL3ωρ 模型是唯一完全符合所有标准的模型。
相变特征：
- 研究发现，对于 EL3ωρ 和 NL3*ωρ 参数化，相变可以在低于 $2.0 n_0$ 的密度下发生。
- 推断出的参数暗示了具有大夸克核的混合星的存在。

意义与主张
本文声称，通过将贝叶斯推断与高精度天体物理数据相结合，可以显著约束解禁夸克物质的性质。主要结论是，只要强子物质由 EL3ωρ 参数化描述，且夸克物质由两参数矢量 MIT 包模型描述，那么具有大夸克核的混合星是符合当前质量、半径和潮汐形变度数据的可行方案。

作者强调，虽然引入超子或额外的矢量通道参数（狄拉克海）在理论上是有动机的，但在特定的似然函数背景下，这些扩展目前会导致违反 pQCD 约束或普适关系的模型。该研究表明，将矢量 MIT 包模型的“两参数”方法与 EL3ωρ 强子 EOS 相结合，为当前的观测图景提供了最稳健的混合星描述。这项工作旨在完善可行的 EOS 模型，并量化控制强子-夸克相变的参数不确定性。

技术摘要：混合中子星的贝叶斯方法研究

类似论文