On the triple nature of the PSR J0435+3233 system

技術要約：PSR J0435+3233系の三重連星性について

問題提起
本論文は、Commensal Radio Astronomy FAST Surveyによって発見されたミリ秒パルサー（MSP）であるPSR J0435+3233のタイミング・エフェメリスにおける重大な異常を取り上げている（Wu et al. 2026、以下W26）。W26は、$4.88 \times 10^{-17}$ s s$^{-1}$という観測されたスピンダウン率（$\dot{P}$）を報告したが、これは典型的な銀河系内MSPよりも2桁高い値である。もしこれが固有のものであるならば、特性年齢は約$10^6$年、表面磁場は$1.26 \times 10^{10}$ G、スピンダウン出力（$\dot{E}$）は$\sim 5.9 \times 10^{37}$ erg s$^{-1}$となり、現在のMSPの形成および進化モデルに疑問を投げかけることになる。さらに、Fermi Large Area Telescope (LAT) はガンマ線パルスを検出しているものの、そのフラックスはW26が予測する高い$\dot{E}$とは一致しておらず、ガンマ線効率が異常に低いことを示唆している。加えて、W26で観測されたスピン周波数および軌道周波数の微分は、顕著な変動（例：5年間で$|\dot{P}|$が約60%増加）と相関した進化を示しており、これは単純な連星系や標準的な「スパイダー」パルサーモデルでは説明が困難である。

手法
著者らは、この異常なタイミング挙動は固有のものではなく、遠方の第三天体の重力場内における内側連星系（パルサーとその低質量伴星）の加速によって引き起こされていると提案している。

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から得られた5つのスピン周波数微分と3つの軌道周波数微分を用いて、著者らはMCMCサンプラーを用い、外側軌道の5つのケプラー要素（投影軌道長半径$x_O$、公転周期$P_{B,O}$、離心率$e_O$、近日点引数$\omega_O$、真近点角$\nu_O$）を推定した。
3. 階層的三重連星タイミングモデル： 著者らは、勾配ベースの最適化にJAXを用いたPINTにおける階層的三重連星タイミングモデルを実装した。このモデルは、内側軌道の位相を計算する前に、外側軌道のロマー遅延を考慮するものである。彼らはこのモデルをW26の重心時刻到着（ToA）に適合させた。
4. ガンマ線統合： 解を精緻化し、タイミングの基底期間を延長するために、著者らはFermi-LATのガンマ線光子をラジオ・タイミング解を用いてフォールディングした。彼らは、ラジオ・タイミングからの事後分布サンプルをガンマ線のH統計量に基づいて再重み付けし、結合対数尤度を最大化することで、タイミングの基底期間を2008年から2024年まで実質的に延長した。
5. 光学的伴星の特定： 著者らは、パルサーの位置をGaia DR3および2MASSカタログと照合して、第三天体の潜在的な光学的伴星を特定し、その固有運動、年周視差、および光度特性を分析した。

主な結果

• 三重連星系の確認： データは階層的三重連星系として最もよく記述される。内側連星は、MSPとその低質量伴星（おそらく白色矮星）からなり、約8日の軌道を持つ。この連星は、離心率$e_O \approx 0.60$、公転周期$P_{B,O} \approx 72$年の広大な離心軌道を持つ第三天体との共通重心の周りを回っている。
• $\dot{P}$の異常の解決： 大きな観測された$\dot{P}$は、主に外側軌道による視線方向の加速度の変化によるアーティファクトである。固有のスピンダウン率は、観測値よりも少なくとも2桁低いことが制約されており、典型的な銀河系内MSPと一致している。
• 拡張されたガンマ線検出： 三重連星タイミングモデルは、Fermi-LATデータの開始時期である2008年8月に遡ってパルサーの位相を予測することに成功した。これは、ラジオデータの範囲外でコヒーレンスを失っていた従来のモデルと比較して大幅な改善である。結合されたラジオ/ガンマ線解は、高いH統計量（$H=171.9$）を示す。
• 第三天体の特定： パルサーから11 masの位置にある恒星2MASS J04353375+3233080が、第三天体の光学的伴星として特定された。Gaiaの年周視差と光度に基づくと、これは距離約1.5〜2.2 kpcに位置する質量約$1.2 M_\odot$の主系列星（スペクトル型G0V–F0V）である。
• 軌道幾何学： 外側軌道の投影軌道長半径は約2477秒であり、軌道傾斜角は$i_O \approx 31^\circ$または$149^\circ$である。モデルは、外側軌道の近日点通過が2036年初頭に起こることを予測している。
• タイミング残差： 新しいモデルは、FASTのToAの残差RMSを、W26の$1.79 \mu$sから$1.49 \mu$sへと減少させ、簡約$\chi^2$は1.002となった。しかし、著者らは高次の効果（例：潮汐相互作用、相対論的時間遅延）が依然として残差に存在している可能性があると指摘している。

意義と主張
本論文は、PSR J0435+3233が階層的三重連星系であることを結論付けており、これによりスピンダウン率とガンマ線効率に関する一見矛盾する問題が解決される。著者らは以下のことを主張している：

1. 本システムは、特異な形成経路や極端な磁場を想定することなく、以前の異常なタイミング微分を自然に説明する。
2. 第三天体の特定と外側軌道のパラメータ決定により、システムの整合性のある物理モデルが可能となる。
3. 本システムは、強重力等価原理（SEP）および自重天体に対する自由落下の普遍性（UFF）をテストするための重要な潜在能力を有している。内側軌道の投影軌道長半径が、よく知られた三重連星系PSR J0337+1715よりも大きいため、PSR J0435+3233は、ノルトヴェト効果に関するUFFの破れに対して、より厳格な制約を与える可能性がある。
4. 今後の観測、特に第三天体の視線速度変動を測定するための高分解能分光観測、および2036年の近日点通過に向けた継続的なパルサー・タイミングは、システムの質量比、傾斜角、および距離を精密に決定し、これらの重力テストをさらに精緻化することを可能にする。