Slowly rotating white dwarfs in f(R,Lm, T) = R + α LmT gravity: moment of inertia, rotational mass shift, and quadrupole

技术摘要：$f(R, L_m, T) = R + \alpha L_m T$ 引力理论中的慢转白矮星

问题陈述
白矮星是 Ia 型超新星的前身，然而观测到的超亮事件表明其前身质量可能超过了标准的钱德拉塞卡极限（$M_{Ch} \approx 1.44 M_\odot$）。虽然已有研究提出通过修改引力部门来解释超钱德拉塞卡质量，但以往在 $f(R, L_m, T)$ 引力中的研究仅限于静态构型。本文旨在填补理解慢转如何影响 $f(R, L_m, T) = R + \alpha L_m T$ 特定模型下白矮星的结构、稳定性及观测性质方面的空白，其中 $\alpha$ 是物质-几何耦合常数，且物质拉格朗日量选取为 $L_m = p$。

方法论
作者采用 Hartle-Thorne 形式，将度规和场方程展开至角速度 $\Omega$ 的二阶。

1. 场方程： 该理论由一个其引力拉格朗日量依赖于里奇标量 $R$、物质拉格朗日量 $L_m$ 和迹 $T$ 的作用量定义。对于特定的模型 $f = R + \alpha L_m T$，场方程简化为 $G_{\mu\nu} = S_{\mu\nu}$，其中有效源 $S_{\mu\nu}$ 在代数上依赖于物质变量（$\rho, p$）和耦合项 $\alpha$。至关重要的是，由于 $f_R = 1$，该理论不会传播额外的标量自由度，且外部真空解与广义相对论（GR）完全一致。
2. 结构方程： 作者推导出了完整的微分方程组：
   • 零阶 ($O(\Omega^0)$)： 一个描述静态背景的修正 Tolman-Oppenheimer-Volkoff (TOV) 系统。
   • 一阶 ($O(\Omega)$)： 关于局部惯性系角速度的框架拖拽方程。
   • 二阶 ($O(\Omega^2)$)： 包括质量偏移和形状函数的单极（$\ell=0$）与四极（$\ell=2$）形变方程。
3. 验证： 所有推导出的闭合关系均通过基于精确有理算术的独立符号残差测试进行了验证，确认场方程和守恒律（$\nabla_\mu S^\mu_\nu = 0$）在字面意义上完全满足。
4. 数值积分： 使用针对全电离 $^{12}\text{C}$ 的现实零温状态方程（EoS）对系统进行积分，该 EoS 包含了相对论简并电子气体和体心立方（bcc）晶格库仑修正。计算是在中心密度和七个耦合参数 $\alpha/u_1 \in \{-1.6, \dots, +1.6\}$ 的网格上进行的。
5. 稳定性分析： 通过求解有效 బ罗托罗皮克流体的基本径向脉动模（$\omega_0^2$）的 Sturm-Liouville 特征值问题来评估径向稳定性，确认转折点判据（$dM/d\rho_c = 0$）在该理论中依然严格成立。

关键结果

• 质量-半径关系： 白矮星的最大质量随正 $\alpha$ 而增加。对于 $\alpha = -1.6 u_1$， $M_{max} \approx 1.34 M_\odot$；对于 $\alpha \approx +0.8 u_1$， $M_{max} \approx 1.51 M_\odot$。
• 质量最大值的本质： 对于 $\alpha \leq 0$，质量序列终止于一个真实的转折点（$dM/d\rho_c = 0$）。然而，对于所有 $\alpha > 0$，序列终止于一个内部正则性边界，即在质量仍在上升时，修正后的静力学方程的分母变为零。因此，所报告的“最大”质量是正则性极限而非动力学不稳定性点。
• 转动惯量： 无量纲转动惯量（$\bar{I} = I/M^3$）随质量和耦合的变化平滑变化。对于一个代表性恒星（$\rho_c = 2 \times 10^9 \text{ g cm}^{-3}$），$\bar{I} \approx 1.20 \times 10^5$。
• 旋转质量偏移： 报告了两个常数：
  • $\delta M_N$（固定重子数）：代表旋转动能（$\delta M_N = \frac{1}{2}\bar{I}$）。该值很小（$\sim 6 \times 10^4$），表明仅靠旋转无法显著提升白矮星质量以接近钱德拉塞卡极限。
  • $\delta M_c$（固定中心密度）：包括在固定密度下支持的额外重子静止质量，其随 $C^{-3}$ 缩放，达到 $\sim 1.5 \times 10^8$。
• 四极矩： 由于在白矮星紧凑度（$C \sim 10^{-4}$）下，渐近模匹配无法恢复物理四极矩——这是因为衰减的物理模相对于增长模被抑制了约 $C^4 \sim 10^{-12}$ 倍，导致数值噪声占主导。相反，作者通过 Radau 方程的良定体积积分计算了无量纲四极矩 $\bar{Q}$，得到代表性恒星的 $\bar{Q} \approx 2.0 \times 10^3$，且沿序列保持 $\bar{Q}C \approx 2$ 的关系。
• 稳定性： 所有报告的构型均处于稳定分支（$\omega_0^2 > 0$）。对于正 $\alpha$，模型在正则性边界前保持稳定。

意义与主张
本文声称提供了首个包含转动惯量、旋转质量偏移和四极矩的 $f(R, L_m, T)$ 引力中的慢转白矮星模型。

• 理论严谨性： 通过证明 $f_R=1$ 将动力学简化为由守恒有效流体驱动的 GR，作者确立了标准的径向脉动稳定性判据适用，避免了通常在非最小耦合理论中出现的歧义。
• 观测一致性： 模型重现了观测到的白矮星质量-半径序列（包括超大质量的 ZTF J1 1901），并与来自 Gaia 数据的关于 $\alpha$ 的静态约束保持一致。
• 方法论贡献： 本工作强调了通过外部匹配提取低紧凑度恒星四极矩时的数值限制，并提出了将 Radau 方程体积积分作为一种鲁棒的替代方案。
• 局限性： 作者明确指出，该理论中的超钱德拉塞卡质量是由耦合项 $\alpha$ 支持的，而非旋转，因为旋转质量偏移极小（$\delta M/M \sim 10^{-4}$）。本研究受限于零温 EoS，且忽略了磁场，而磁场被认为是未来为了完全分离物质-曲率耦合与磁流体静力学支撑所必需的研究方向。