On the triple nature of the PSR J0435+3233 system

技术摘要：关于 PSR J0435+3233 三星系统的研究

问题陈述
本文针对 PSR J0435+3233 脉冲星计时历元中的显著异常进行了研究。该毫秒脉冲星（MSP）由 Commensal Radio Astronomy FAST 巡天（Wu et al. 2026，以下简称 W26）发现。W26 报告观测到的自转自降率（$\dot{P}$）为 $4.88 \times 10^{-17}$ s s$^{-1}$，这比典型的银河系毫秒脉冲星高出两个数量级。如果该数值是内在的，则意味着其特征年龄约为 $10^6$ 年，表面磁场为 $1.26 \times 10^{10}$ G，且自转能耗（$\dot{E}$）约为 $\sim 5.9 \times 10^{37}$ erg s$^{-1}$，这将挑战当前的毫秒脉冲星形成与演化模型。此外，尽管 Fermi 大面积望远镜（LAT）探测到了 $\gamma$ 射线脉冲，但其通量与 W26 预测的高 $\dot{E}$ 值不符，表明其 $\gamma$ 射线效率异常低。此外，W26 中观测到的自转频率和轨道频率导数表现出显著的变率（例如，在五年内 $|\dot{P}|$ 增加了约 60%）以及相关的演化关系，这难以用简单的双星系统或标准的“蜘蛛”脉冲星模型来解释。

方法论
作者提出，这种异常的计时行为并非内在属性，而是由于内层双星系统（脉冲星及其低质量伴星）在远距离第三体（tertiary companion）的引力场中发生加速所致。

1. 频率导数分析： 作者分析了 W26 中报道的自转频率（$f$）与轨道频率（$f_B$）导数的关联演化。他们证明了 $\dot{f}$ 与 $\dot{f}_B$ 的近乎一致的变化遵循关系式 $\dot{f}/f \simeq \dot{f}_B/f_B \simeq -a_l/c$，其中 $a_l$ 是由于第三体引起的视线方向加速度。这排除了圆轨道（圆轨道要求高阶导数出现符号反转），但支持了一个偏心外轨道。
2. 开普勒参数估计： 利用 W26 中的五个自转频率导数和三个轨道频率导数，作者使用马尔可夫链蒙特卡洛（MCMC）采样器估计了外轨道的五个开普勒参数（投影半长轴 $x_O$、周期 $P_{B,O}$、离心率 $e_O$、近星点幅角 $\omega_O$ 和真近点角 $\nu_O$）。
3. 分级三体计时模型： 作者在 PINT 中实现了一个分级三体系统计时模型（使用 JAX 进行基于梯度的优化）。该模型在计算内轨道相位之前，先考虑了外轨道的罗默延迟（Rømer delay）。作者利用该模型对来自 W26 的重心到达时间（ToAs）进行了拟合。
4. $\gamma$ 射线整合： 为了精细化解并延长计时基线，作者利用射电计时解对 Fermi-LAT 的 $\gamma$ 射线光子进行了折叠。作者根据 $\gamma$ 射线 H-统计量对来自射电计时的后验样本进行重新加权，以最大化联合对数似然，从而将计时基线从 2008 年有效延长至 2024 年。
5. 光学对应体识别： 作者将脉冲星位置与 Gaia DR3 和 2MASS 星表进行交叉匹配，以识别潜在的第三体光学对应体，并分析了其自行、视差及光度特性。

主要结果

• 三星系统确认： 数据最符合一个分级三体系统的特征。内层双星由该毫秒脉冲星及其低质量伴星（可能为白矮星）组成，轨道周期约为 8 天。该双星绕共同质心运行，并围绕一个具有宽阔偏心轨道的第三体运行（$P_{B,O} \approx 72$ 年，$e_O \approx 0.60$）。
• 解决 $\dot{P}$ 异常： 观测到的巨大 $\dot{P}$ 主要是由外轨道变化的视线方向加速度造成的伪影。其内在的自转降率被限制在比观测值低至少两个数量级的水平，这与典型的银河系毫秒脉冲星一致。
• 扩展的 $\gamma$ 射线探测： 三星系统计时模型成功预测了回到 2008 年 8 月 Fermi-LAT 数据开始时的脉冲相位，这相比于以往在射电数据跨度之外会失去相干性的模型有了显著改进。联合射电/$\gamma$ 射线解得到了极高的 H-统计量（$H=171.9$）。
• 第三体识别： 位于脉冲星 11 mas 处的恒星 2MASS J04353375+3233080 被确定为第三体的光学对应体。根据 Gaia 视差和光度，该恒星是一颗主序星（光谱类型 G0V–F0V），质量约为 $1.2 M_\odot$，距离约为 $1.5$–$2.2$ kpc。
• 轨道几何结构： 外轨道的投影半长轴约为 $2477$s，倾角$i_O \approx 31^\circ$或$149^\circ$。模型预测外轨道的近星点经过将发生在 2036 年初。
• 计时残差： 新模型将 FAST ToAs 的残差均方根（RMS）从 W26 中的 $1.79 \mu$s 降低至 $1.49 \mu$s，约化 $\chi^2$ 为 1.002。然而，作者指出，残差中可能仍存在高阶效应（例如，潮汐相互作用、相对论性时间膨胀）。

意义与主张
论文得出结论，PSR J0435+3233 是一个分级三体系统，解决了其自转降率和 $\gamma$ 射线效率之间看似矛盾的问题。作者主张：

1. 该系统为此前异常的计时导数提供了一个自然的解释，无需诉诸奇异的形成通道或极端的磁场。
2. 第三体伴星的识别以及外轨道参数的确定，使得该系统能够建立起一个一致的物理模型。
3. 该系统对于测试强等效原理（SEP）和自由落体普遍性（UFF）具有重要潜力。由于其内轨道的投影半长轴大于著名的三星系统 PSR 0337+1715，PSR J0435+3233 有潜力为违反 UFF（特别是诺德特效应，Nordtvedt effect）提供更严格的约束。
4. 未来的观测，特别是通过高分辨率光学光谱测量第三体恒星的径向速度变化，以及持续接近 2036 年近星点经过时的脉冲计时，将能够精确测定系统的质量比、倾角和距离，从而进一步完善这些引力测试。