Bayesian Inference of Neutron Star Properties in $f(Q)$ Gravity Using NICER Observations

技术摘要：利用 NICER 观测数据进行 $f(Q)$ 引力中中子星属性的贝叶斯推断

问题陈述
虽然广义相对论（GR）成功描述了各种尺度下的引力现象，但在概念上仍面临奇异点、暗物质和暗能量方面的挑战，这促使了对修正引力（MG）理论的探索。在这些理论中，对称平行引力（STEGR）及其扩展形式 $f(Q)$ 引力提供了一种独特的几何表述，其中引力是通过时空的非度量性（non-metricity）而非曲率或挠率来编码的。尽管 $f(Q)$ 模型已在宇宙学领域得到了研究，但其在中子星（NSs）强场机制中的应用在很大程度上仍处于空白。具体而言，目前缺乏系统性的贝叶斯推断研究来利用观测数据约束 $f(Q)$ 模型参数，因为以往的分析往往依赖于现象学参数的选择，而缺乏与数据进行严谨统计对抗的过程。

方法论
作者通过考虑三种代表性的函数形式（线性、对数和指数模型），研究了 $f(Q)$ 引力框架下的中子星内部结构。为了将修正引力的效应与致密物质物理的不确定性分离，本研究采用了现实的、密度依赖的相对论平均场 DDME2 物态方程（EoS）。

研究方法如下：

1. 理论框架： 作者推导了 $f(Q)$ 引力中静态、球对称时空的修正 Tolman–Oppenheimer–Volkoff (TOV) 方程。这些方程控制恒星的平衡，并使用四阶龙格-库塔（RK4）方案进行数值求解。
2. 贝叶斯推断： 采用贝叶斯框架来约束三种 $f(Q)$ 模型的自由参数（$\alpha, \beta, Q_0$）。该分析利用了来自 NICER 对四个脉冲星（PSR J0030+0451, PSR J0740+6620, PSR J0437+4715, 和 PSR J0614+3329）的质量-半径后验分布。
3. 统计分析： 研究计算了模型参数及衍生中子星属性（质量、半径、潮汐形变度和紧凑度）的后验分布。通过使用贝叶斯因子（BF）进行模型比较，以确定哪种函数形式在统计上更受数据青睐。

核心贡献

• 首次实现 $f(Q)$ 的贝叶斯约束： 本工作代表了首次利用中子星观测量对 $f(Q)$ 引力模型参数进行的系统性贝叶斯推断分析，将研究范围从宇宙学测试扩展到了强场机制。
• 对比模型分析： 本文对线性、对数及指数 $f(Q)$ 模型进行了严谨的统计比较，评估了它们在面对高精度 NICER 数据时的可行性。
• 参数简并与相关性： 研究绘制了模型参数与中子星属性之间的相关性图谱，揭示了不同的函数形式如何影响引力修正与恒星结构之间的耦合关系。

结果

• 模型偏好： 贝叶斯分析（由贝叶斯因子比较支持）在统计上更倾向于指数 $f(Q)$ 模型，而非线性或对数模型。指数模型表现出参数约束良好且接近高斯、单峰的后验分布。相比之下，线性模型显示出明显的双峰性和宽泛的参数分布，表明存在强烈的简并性；而对数模型则表现出中等程度的不确定性。
• 中子星属性： 对于首选的指数模型，研究得出的 $1.4 M_\odot$ 处的半径和潮汐形变为 $R_{1.4} = 11.27^{+0.53}_{-0.36}$ km 和 $\Lambda_{1.4} = 156.95^{+84.02}_{-41.73}$。这些数值与当前的观测约束是一致的。
• 最大质量与质量间隙： 一个显著的发现是，所有三个受约束的模型都预测最大中子星质量可达 $M_{\text{max}} \simeq 2.98 M_\odot$。质量-半径关系的 95% 置信区域延伸到了较低的质量间隙区域（$\sim 2.5\text{--}5 M_\odot$）。这一预测是在没有预设构造的情况下，从贝叶斯约束的参数空间中自然产生的。
• 相关性： 指数模型展示了最大质量与其他中子星属性之间强烈的结构相互依赖性，而线性模型则显示出最大质量与其他观测量之间的统计独立性。

意义
本文声称，这些结果强调了中子星作为对称平行引力在强场机制下有力探测工具的潜力。在低质量间隙区域出现的稳定构型被视为 $f(Q)$ 框架的一个明确且可测试的预测，这将其与通常产生较低最大质量的标准广义相对论场景区分开来。研究结论认为，指数 $f(Q)$ 模型是三种测试模型中对中子星结构描述在统计上最受支持的一种，这表明修正引力效应可以显著增强高密度引力支撑。作者指出，未来在质量间隙区域的精确质量测量可以作为区分广义相对论与扩展引力场景的有力判据。他们还建议，未来的工作可以涉及与灵活物态方程参数化的联合贝叶斯推断，以进一步解决强场引力与微观物理之间的简并问题。