Bayesian inferences on covariant density functionals from multimessenger astrophysical data: Influences of parametrizations of density dependent couplings

本研究采用结合多信使天体物理数据的贝叶斯框架，证明尽管协变密度泛函中密度依赖耦合的不同参数化形式会得出大体相似的推断，但具体的函数形式会显著影响超饱和密度下的状态方程与对称能，因此需要在同位旋矢量通道中引入直至曲率系数$K_{sym}$的扩展灵活性以实现精确建模。

要点

想象宇宙中充满了存在于中子星内部的宇宙“超级物质”。这种物质密度极高，仅一茶匙的重量就相当于一座山。物理学家称这种物质为致密核物质。为了理解这种物质的行为，他们使用称为**协变密度泛函（CDFs）**的数学“配方”。可以将这些配方视为构建恒星内部模型的蓝图。

然而，这些蓝图并不完美。它们依赖于科学家必须调节的“旋钮”和“表盘”（参数）。本文提出的核心问题是：我们编写这些表盘指令的具体方式，究竟是否重要？

以下是研究人员所做工作及发现结果的简明概述：

1. 问题：编写配方的方式过多

过去，科学家主要使用一种特定的指令来说明物质密度如何变化。他们假设“旋钮”仅对粒子数量（即像数房间里有多少人那样）做出反应。

但测量密度的另一种方式是观察粒子之间的相互作用（即像看它们抱得有多紧那样）。研究人员希望探究：改变密度测量的类型（计数与拥抱）或改变指令的数学形状（使用直线与曲线），是否会极大地改变我们对中子星的认知图景。

2. 实验：贝叶斯“味觉测试”

团队使用了一种称为贝叶斯推断的强大统计方法。想象你是一位试图完善汤食谱的厨师。你有一份约束清单：

• 汤必须足够咸（如同大质量脉冲星的质量）。
• 汤必须足够浓稠（如同通过 X 射线望远镜测量的中子星大小）。
• 汤在搅拌时必须表现出特定的行为（如同引力波数据）。

他们尝试了六种不同版本的食谱（不同的密度依赖性数学公式）。他们将所有最新的天文数据（来自引力波、X 射线望远镜和粒子实验）输入计算机，以查看哪些食谱能做出满足所有约束条件的“汤”。

3. 结果：什么变了，什么没变？

“大局”变化不大
令人惊讶的是，无论是计数粒子还是测量相互作用，中子星的最终图景几乎相同。

• 类比：想象你试图猜测一个神秘盒子的重量。无论你使用电子秤还是弹簧秤，得到的结果都是一样的。
• 发现：所有不同食谱预测的中子星最大重量（质量）和尺寸（半径）几乎完全相同。恒星基本结构的“旋钮”足够灵活，无论使用何种具体数学方法，都能根据数据进行调节。

“隐藏成分”确实发生了变化
虽然恒星的外部看起来相同，但“汤”内部发生的情况却不同。

• 类比：两块蛋糕外观可能完全一样，但一块是用黄油做的，另一块是用油做的。仅凭肉眼无法分辨，但它们的质地和冷却方式却不同。
• 发现：不同的食谱对对称能（一种决定混合物中质子与中子比例的性质）的行为做出了不同的预测。
  • 有些食谱暗示恒星核心会有大量质子（如同高糖蛋糕）。
  • 另一些则暗示质子极少（如同低糖蛋糕）。
  • 这至关重要，因为质子的数量决定了恒星冷却的速度。如果有足够的质子，恒星就能非常迅速地“尖叫”着释放能量（这一过程称为直接 Urca 过程）。

4. 结论：我们需要更好的工具

本文结论如下：

1. 当前数据足以告诉我们中子星的一般尺寸和重量，无论我们使用哪种具体的数学食谱。
2. 当前数据不足以告诉我们深藏内部的“隐藏成分”（对称能）究竟在做什么。不同的食谱都能拟合当前的观测结果，但它们关于恒星内部成分的叙述却各不相同。

核心启示：
要真正理解中子星内部致密物质的“风味”，仅靠尺寸和重量测量是不够的。我们需要观察恒星的新方法，例如观测它们随时间冷却的方式。在此之前，恒星内部的“食谱”仍是一个谜，有几个不同的版本看起来都同样合理。

技术摘要

技术摘要：基于多信使天体物理数据的协变密度泛函贝叶斯推断

问题陈述
致密强子物质的行为，特别是致密星体（CSs）内部的行为，仍然是核物理与天体物理领域的核心挑战。核物质的状态方程（EOS）构成了微观核子动力学与宏观恒星性质之间的桥梁。尽管协变密度泛函（CDF）理论为描述从有限原子核到致密星体的核数据提供了统一框架，但在超饱和密度下的状态方程仍存在显著的不确定性。这种不确定性的一个关键来源在于密度依赖的介子 - 核子耦合的参数化。现有模型通常依赖于特定的函数形式（如有理函数或指数函数），并取决于矢量密度或标量密度。在当前的多信使约束下，这些特定选择（函数形式及密度依赖类型：矢量与标量）在多大程度上影响了对致密物质性质和核对称能的推断，需要进行系统性的调查。

方法论
作者采用贝叶斯推断框架，系统地比较了不同类型的“密度依赖”（DD）CDF 模型。该研究利用了一致的计算框架，并基于以下推导出的相同多物理约束：

1. 地面约束： 饱和态下的核物质特性（有效质量、饱和密度、结合能、不可压缩性、对称能斜率）以及手征有效场论（$\chi$EFT）对纯中子物质的ab initio预测，范围高达约 1.5-2 倍饱和密度。
2. 天体物理约束： 大质量脉冲星的质量测量（如 PSR J0348+0432）、双中子星并合事件（GW170817、GW190425）的潮汐形变，以及 NICER X 射线天文台对四颗脉冲星（PSR J0740+6620、J0030+0451、J0437-4715、J1231-1411）的同步质量 - 半径估算。

该研究改变了 CDF 参数化的三个关键方面：

• 密度类型： 耦合是依赖于矢量密度（$n_v$）、标量密度（$n_s$），还是混合依赖。
• 函数形式： 密度依赖采用有理函数（R）还是指数函数（E）。
• 参数自由度： 同位旋标量（$\sigma, \omega$）和同位旋矢量（$\rho$）通道中允许的独立参数数量。

研究探讨了六种具体模型（论文中的表 I），包括标准 VRE 模型（矢量密度、同位旋标量有理函数、同位旋矢量指数函数）及其扩展模型，如 MRE（混合密度）、MRE2（增加同位旋标量自由度）、VRR（同位旋矢量有理函数）和 VRR2（增加同位旋矢量自由度）。贝叶斯分析对耦合强度和函数参数赋予均匀先验，生成了核特征系数（如 $K_{sat}, Q_{sat}, L_{sym}, K_{sym}$）和恒星观测量的后验分布。

主要贡献

• 参数化的系统比较： 该工作全面比较了密度依赖类型和函数形式不同的 CDF 模型，隔离了这些选择对高密度物质性质的影响。
• 同位旋通道的有理函数参数化： 作者引入并测试了同位旋矢量 $\rho$-介子耦合的有理函数参数化（VRR 和 VRR2 模型），与标准的指数形式相比，这允许在高密度下更灵活地描述对称能。
• 模型依赖性的量化： 该研究量化了推断出的致密物质性质在多大程度上依赖于密度依赖的具体函数形式，而非观测数据提供的约束。

结果

• 全局恒星性质： 分析表明，具有不同密度依赖类型（矢量与标量）或同位旋标量通道不同函数形式的 CDF 模型，在 95.4% 置信区间内，对全局致密星体性质（最大质量、半径、潮汐形变）产生的后验分布几乎相同。同位旋标量与同位旋矢量通道之间的相互作用对热力学的影响，弱于单一通道内的变化。
• 同位旋标量通道： 允许同位旋标量通道中的核饱和特性自由调整，特别是调整至偏度系数 $Q_{sat}$，提供了足以重现当前数据的灵活性。增加同位旋标量自由度的模型（如 MRE2）预测大质量恒星（$M \gtrsim 1.4 M_\odot$）的半径略大，且高密度下的 EOS 更硬，这主要由 $Q_{sat}$ 和峰度系数 $Z_{sat}$ 的更宽分布所驱动。
• 同位旋矢量通道： 同位旋矢量耦合参数化的变化对全局恒星性质影响适中，但导致对称能的密度依赖和粒子组分出现显著差异。
  • 与指数形式（VRE）相比，有理函数参数化（VRR、VRR2）允许恒星核心具有更大的质子分数。
  • VRR2 模型为 $\rho$-介子耦合引入了第三个独立参数，允许曲率系数 $K_{sym}$ 变化。该模型允许在低至 $1.4 M_\odot$ 质量的恒星中开启直接 Urca（DU）冷却过程，而 VRE 模型则倾向于在 $M \lesssim 2.0 M_\odot$ 时抑制该过程。
• 对系数的约束： 当前的天体物理数据对低阶饱和参数的约束较弱。高阶系数（$K_{sym}, Q_{sym}$）是同位旋矢量部门中区分模型的主要依据。在 VRR2 模型中，发现对称能斜率 $L_{sym}$ 与曲率 $K_{sym}$ 之间的相关性可忽略不计。

意义与主张
该论文声称，尽管不同的参数化对全局恒星性质的推断产生了大致可比的结论，但在超饱和密度下，状态方程和对称能的差异仍然显著，反映了其对所选函数形式的敏感性。作者得出结论：

1. 同位旋标量建模： 当允许同位旋标量饱和性质（直至 $Q_{sat}$）自由调整时，当前对核物质和中子星物质的建模具有足够的灵活性。
2. 同位旋矢量细化： 同位旋矢量通道需要进一步细化。将自由度扩展至曲率系数 $K_{sym}$ 对于捕捉高密度下的对称能和粒子组分变化是必要的。
3. 未来方向： 同位旋矢量介子耦合的有理函数参数化为未来旨在探测同位旋不对称性的致密星体和重离子碰撞模拟提供了合适的框架。然而，约束高密度对称能需要除积分恒星参数之外的互补策略，例如研究致密星体冷却和直接 Urca 过程的开启。

该工作通过实施由多信使数据启发的有理函数参数化，推进了先前的研究，证明了尽管低阶参数受到约束，但高密度行为仍然对密度依赖的具体函数假设保持敏感。