On the braking index evolution of PSR B0540-69: wind braking of pulsars

技術要約：PSR B0540-69の制動指数進化について：パルサーのウィンド制動

問題提起
パルサーPSR B0540-69（本文中ではPSR B0540-69とされているが、アブストラクトにはPSR B0450-69という誤植がある）は、実測された制動指数と、記録されたスピンダウン状態の変化の両方を備えているため、パルサー物理学における特異なケーススタディを提示している。スピンダウン率が36%増加した状態変化に続き、このパルサーは制動指数（$n$）の時間依存的な進化を示した。先行研究（Paper I, Kou et al. 2016）では、もしこの状態変化が粒子ウィンドの増強によって引き起こされたのであれば、新しい高スピンダウン状態における制動指数は、真空ギャップモデルに基づくと、$n \approx 1.79$（永続状態の値である$n=2.13$よりも低い値）に落ち着くと予測されていた。しかし、観測データは、制動指数が非常に低い値（$n \approx 0.03$）から$n \approx 1.6$まで進化し、理論的予測に近づいていることを示している。本研究で取り組む中心的な問題は、この遷移過程における制動指数の過渡的な進化を説明すること、および制動メカニズムを制約するために最終的な定常状態の値を決定することである。

手法
本論文は、中性子星の回転エネルギー損失を磁気双極放射と粒子ウィンドの流出の両方に帰属させる現象論的な枠組みであるウィンド制動モデルを採用している。このモデルは、傾斜角（$\alpha$）、ゴールドリッヒ・ジュリアン密度に対する粒子数密度の比（$\kappa$）、および粒子加速電位（$\Delta\phi$）を組み合わせた無次元パラメータ$\eta$を定義する。

著者らは、観測される制動指数の進化は、平均スピンダウン率の変化（高状態の状態は維持されている）によるものではなく、むしろ粒子数密度の微小かつ過渡的な変動によるものであると提案している。具体的には、新しい状態における粒子密度$\kappa$には、以下の減衰摂動が含まれると仮定している：
$$\kappa(t) = \kappa_1(1 + \delta e^{-t/\tau})$$
ここで、$\delta$は微小な振幅であり、$\tau$は減衰時定数である。スピンダウン方程式の時間微分を取ることで、著者らは、この過渡的な密度変動が角速度の二階微分（$\ddot{\Omega}$）にどのように寄与するかを導出した。これは、実測される制動指数に直接影響を与える。全制動指数は、定常状態の値（$n_{new}$）と過渡的な指数項の和としてモデル化される：
$$n(t) = n_{new} + \Delta n e^{-t/\tau}$$

主要な貢献と結果

1. 過渡的進化のモデリング： 本論文は、粒子密度の指数関数的に減衰する微小な変動（$\delta \ll 1$）が、パルサーの特性年齢と減衰時定数の比（$\tau_c/\tau$）が大きい場合、制動指数の顕著な過渡的変動を引き起こし得ることを導き出した。これにより、スピンダウン率が安定している一方で、なぜ制動指数が進化するのかが説明される。
2. 定量的適合： 観測データ（Ho et al. 2026）へのモデルのフィッティングにより、制動指数の進化は指数関数的減衰の形式に従うことが判明した。最良適合パラメータは、初期の過渡的振幅 $\Delta n \approx -2.1$、および減衰時定数 $\tau \approx 1708$ 日である。これは、粒子密度の変動振幅 $\delta \sim 10^{-3}$ に相当する。
3. 理論的予測への収束： モデルは、制動指数が近い将来、一定の定常状態の値に落ち着くことを予測している。この最終値は、真空ギャップモデルに基づくと $n_{new} \approx 1.79$ と予測され、これは永続状態の値（$n=2.13$）よりも低く、ウィンド制動の仮説と一致している。
4. 代替モデルとの差別化： 本論文は、制動指数の進化に関する代替的な説明（進化する磁場 $\dot{B}$、慣性モーメント $\dot{I}$、または傾斜角 $\dot{\alpha}$ など）は、最終的な定常状態の値に関する予測能力を欠いていると主張している。対照的に、ウィンド制動モデルは、粒子加速電位に基づく、具体的かつ検証可能な最終制動指数の予測を提供する。

意義
本論文は、PSR B0540-69の制動指数の進化が、パルサーの制動メカニズムの重要な診断ツールとなることを主張している。制動指数の進化が、ウィンド制動の枠組みにおける粒子密度の過渡的な変動と一致することを実証することで、本研究は、スピンダウン状態の変化が間欠的パルサーで見られるメカニズムと同様に、強化された粒子ウィンドによって引き起こされたという仮説を支持している。

主な意義は、このモデルの将来的な検証可能性にある。本論文は、PSR B0540-69の継続的なモニタリングによって、最終的に定常状態の制動指数が明らかになるであろうと断じている。もし、測定された値が予測された $n \approx 1.79$（あるいは特定の加速電位によって決定される値）に収束すれば、それはウィンド制動メカニズムを明白に裏付け、パルサ・マグネットスフィア内の粒子加速電位を制約することになる。これにより、最終的な制動指数に対する具体的な制約を与えない、磁場や幾何学的進化に関わる競合理論から、ウィンド制動モデルを区別することが可能となる。