論文要約：ポスト・グリッチ回復と中性子星構造：ヴェラパルサー

論文タイトル: Post-glitch Recovery and the Neutron Star Structure: The Vela Pulsar
掲載誌: Publications of the Astronomical Society of Australia (2025)
著者: Himanshu Grover ら

1. 研究の背景と問題提起

パルサー（特に中性子星）の回転速度の急激な増加（グリッチ）と、その後の緩やかな回復過程は、中性子星内部の超流体ダイナミクスを理解する上で極めて重要である。特に、ヴェラパルサー（PSR J0835–4510）は頻繁に大きなグリッチを起こすことで知られ、その回復過程は「渦糸クリープモデル（vortex creep model）」を用いて説明されることが多い。

しかし、従来のモデルでは説明しきれない「渦糸残差（vortex residuals：観測値とモデル予測値の差分）」や、グリッチの大きさ（マグニチュード）と次のグリッチ発生までの時間との相関関係など、未解明の現象が存在する。本研究では、2016 年 9 月から 2025 年 1 月までの約 100 ヶ月にわたる観測データ（4 回のグリッチを含む）を用いて、ヴェラパルサーの回転挙動を詳細に分析し、以下の点を解明することを目的とした：

1. 渦糸クリープモデルによるポスト・グリッチ回復の精密な記述。
2. 渦糸残差を「渦糸曲げモデル（vortex bending model）」の文脈で解釈し、振動現象の解明。
3. グリッチの大きさと次のグリッチ発生までの時間の相関の検証。
4. ブレーキング指数（braking index）の推定。

2. 研究方法

観測データ

• 対象: ヴェラパルサー（PSR B0833–45）。
• 期間: 2016 年 9 月〜2025 年 1 月（MJD 57651〜60700）。
• 施設:
  • Ooty 電波望遠鏡 (ORT): 326.5 MHz で観測。主にグリッチ検出とタイミングノイズのモデル化に使用。
  • アップグレード・ジャイアント・メータ波電波望遠鏡 (uGMRT): 250-500 MHz および 550-950 MHz で観測。分散測定（DM）の補正や、ORT の観測空白期間の埋め合わせに使用。
  • パルサーアーカイブデータ (ATNF/Parkes): 特定の期間（特に 2 番目のグリッチ周辺）のデータ欠損を補完するために使用。
• データ処理: DSPSR, PINTA, PSRCHIVE, TEMPO2 を用いて到達時刻（ToAs）を生成し、精密タイミング解析を実施。

解析手法

1. ベイズ統計解析:
   • 異なる回復モデル（指数関数的回復、線形緩和、およびその組み合わせ）を比較するために、DYNESTY（ネストド・サンプリング）と EMCEE（MCMC）を使用。
   • ベイズ証拠（Bayesian evidence）に基づき、最も適切なモデルを選択。
2. 渦糸クリープモデルの適用:
   • 観測されたスピンダウン率の変化を、指数関数項（複数の成分）と線形項（非線形クリープ領域）の和としてモデル化。
   • 各層の慣性モーメントの分数（FMI）と再結合タイムスケールを推定。
3. 渦糸残差の分析（渦糸曲げモデル）:
   • クリープモデルで説明しきれない残差（残差）を解析。
   • 一般化ロムブ・スカルグル周期図（Generalized Lomb-Scargle periodogram）を用いて、残差中の周期的振動を検出。
   • 減衰正弦波モデル（Hypothesis II）などを仮定し、振動の振幅、位相、減衰定数を推定。
4. 相関分析とブレーキング指数:
   • 過去の 26 回のグリッチデータを用いて、グリッチマグニチュードと次のグリッチまでの時間との相関を調査（ピアソン、スピアマン、ケンダル相関係数）。
   • グリッチ前のスピンダウン率の 1 階微分値の線形フィットから、2 階微分値を推定し、ブレーキング指数 $n$ を算出。

3. 主要な結果

4 回のグリッチの回復特性

本研究では、以下の 4 つのグリッチ（MJD 57734.4, 58517, 59417.6, 60429.9）の回復過程を詳細に解析した。

• モデル選択: ベイズ解析により、すべてのグリッチにおいて「指数関数的回復（3 成分）＋線形緩和」を組み合わせたモデル（Model 3c）が最も優れていることが示された。
• パラメータ推定:
  • 各指数成分の慣性モーメント分数（FMI）と減衰タイムスケール（数日から数百日）を高精度で推定。
  • 非線形クリープ領域の FMI と再結合タイムスケールも推定され、これらが次のグリッチ発生までの時間予測に寄与することが確認された。
• 次のグリッチ予測: 渦糸クリープモデルに基づく理論予測とベイズ推定は、観測されたグリッチ間隔とよく一致した。最新のグリッチ（MJD 60429.9）に対する次のグリッチ発生予測は、MJD 61377.7（2026 年 12 月 3 日頃）であり、68% 信頼区間は MJD 61249〜61506 である。

渦糸残差と減衰振動

• 渦糸クリープモデルからの残差を解析した結果、明確な正弦波状の振動（渦糸曲げに起因すると考えられる）が検出された。
• G1 (MJD 57734.4): 周期約 315 日の減衰振動（減衰定数 $\tau \approx 266$ 日）が検出され、現在のグリッチに起因することが確認された。
• G3 (MJD 59417.6): 周期約 344 日と 153.5 日の 2 つの有意な振動成分が検出された。
• これらの振動は、渦糸の曲げ（vortex bending）や非線形クリープ領域の励起によるものである可能性が高い。

グリッチマグニチュードと間隔の相関

• 大きなグリッチ（マグニチュード $> 10^{-7}$）に限定した場合、グリッチの大きさと次のグリッチまでの時間間に正の相関（ピアソン相関係数 0.647）が確認された。
• これは、大きなグリッチほど次のグリッチまでの時間が長くなる傾向があることを示唆する。
• ただし、小さなグリッチを含めるとこの相関は消失するため、小さなグリッチがタイミングノイズと区別できない場合があることへの注意が必要である。

ブレーキング指数

• グリッチ前の安定したスピンダウン率のデータから、ヴェラパルサーのブレーキング指数を $n = 2.94 \pm 0.55$ と推定した。
• この値は、プラズマ充填された磁気圏における磁気双極子放射モデルの予測に近い。以前の研究（1.4〜2.8 の範囲）と比較して、より高い精度で 3 に近い値が得られた。

4. 結論と意義

本研究は、ヴェラパルサーの 4 回のグリッチに対する包括的な分析を行い、以下の重要な知見をもたらした：

1. モデルの精緻化: ポスト・グリッチ回復は、単一の指数関数ではなく、複数の指数成分と線形緩和を組み合わせたハイブリッドモデルで最もよく記述されることを実証した。
2. 内部構造への洞察: 渦糸残差における減衰振動の検出は、渦糸曲げモデルの妥当性を支持し、中性子星内部の超流体と結晶核の結合メカニズムに関する新たな制約を提供した。
3. 予測可能性: 渦糸クリープモデルとベイズ推定を用いることで、次のグリッチ発生時期を高い精度で予測できる可能性を示した。
4. 相関の発見: 大きなグリッチのマグニチュードと次のグリッチまでの時間の正の相関を発見し、グリッチ発生メカニズムの予測モデル（Antonelli et al. 2022 などのモデル）への重要な制約条件を提供した。
5. ブレーキング指数: 従来の研究よりも制約の厳しい値（$2.94 \pm 0.55$）を報告し、パルサーの回転進化と磁場構造に関する理解を深めた。

これらの結果は、中性子星の超流体ダイナミクス、内部構造、およびグリッチの発生・回復メカニズムの理論的枠組みを強化するものであり、将来の他のパルサー（蟹座パルサーなど）への応用や、自動化されたタイミング解析ツールの開発への道を開くものである。