Temporal evolution of the circumstellar disk orientation in the transient X-ray pulsar GRO J1008-57

技術的サマリー：遷移性 X 線パルサー GRO J1008−57 における星周円盤の方位の時間的進化

問題提起
Be/X 線連星（BeXRB）は、通常ペリオン付近で発生する周期的なタイプ I（通常）アウトバーストと、 sporadic なタイプ II（巨大）アウトバーストという、2 種類の明確な X 線アウトバーストを示す。タイプ I アウトバーストは一般的に、中性子星が Be 星の星周脱出円盤に最も接近することによって引き起こされると考えられているが、これらのアウトバーストの軌道位相は常に安定しているわけではない。EXO 2030+375 などの一部のシステムでは、タイプ I アウトバーストの軌道位相が、位相の遅れと進みの周期的な交互を伴う複雑な長期的変動を示し、これはしばしばタイプ II アウトバーストによって引き起こされる。これらの変動は、Be 円盤の構造的・物理的進化を反映していると考えられている。しかし、これらの位相シフトを駆動するメカニズムが、円盤の歳差運動、離心率振動（Kozai-Lidov）、あるいは密度擾乱のいずれに関与するかについては、依然として議論の余地がある。長い軌道周期（249.48 日）を持ち、複数のタイプ II アウトバーストの歴史を有する BeXRB である GRO J1008−57 は、より短い周期のシステムで観測された位相進化のメカニズムが、異なる力学的時間スケールを持つシステムにも適用されるかどうかを検証するための重要なケーススタディとして機能する。

手法
著者らは、2005 年から 2026 年までの期間をカバーする、Swift/BAT（15–150 keV）および MAXI/GSC（2–30 keV）装置から得られた GRO J1008−57 の長期的 X 線監視データを分析した。

1. 光変曲線解析: 日次光変曲線を用いて、31 個のタイプ I アウトバースト（Swift 観測）と 24 個（MAXI 観測）、ならびに 2012 年、2015 年、2017 年、2020 年に発生した 4 つの主要なタイプ II アウトバーストを同定した。
2. タイミングモデリング: タイプ I アウトバーストの時間プロファイルを、簡略化された急上昇指数減衰（FRED）関数を用いてモデル化し、正確な開始時刻（$t_{start}$）および極大時刻（$t_{peak}$）を決定した。
3. 位相計算: 軌道位相は、Kühnel ら（2013）の暦から導出されたペリオン通過時刻を基準として計算した。
4. 進化モデリング: アウトバースト開始の軌道位相の長期的進化を、線形ドリフト、二次ドリフト、およびステップ状のモデルの 3 つのモデルに対してフィッティングした。最適なフィッティングモデルを統計的に選択するために、赤池情報量基準（AIC）を用いた。
5. エネルギー収支: 各タイプ I アウトバーストの全エネルギー放出（フラエンス）は、極大を中心とした$\pm0.15$軌道周期のウィンドウ内で、フィットされたモデル曲線を積分することで推定した。

主要な結果

1. ステップ状の位相進化: 解析により、GRO J1008−57 におけるタイプ I アウトバーストの軌道位相は、連続的な線形または二次ドリフトに従わないことが明らかになった。代わりに、ステップ状の進化を示す。位相は多くの軌道周期にわたって安定しているが、観測された 4 つのタイプ II アウトバーストのそれぞれと時間的に一致する、急激で小振幅のジャンプ（位相で$\sim -0.003$から$-0.016$の範囲）を経験する。
2. モデル選択: ステップ状モデルは、線形モデルおよび二次モデルと比較して、Swiftおよび MAXI の両方のデータに対して統計的に優れたフィッティングを提供し、$\Delta$AIC 値は$5\sigma$を超える有意性を示した。
3. エネルギーサイクル: タイプ I アウトバーストの全エネルギー放出は、体系的なパターンを示す：タイプ II アウトバーストに先立って緩やかに増加し、それに続いて急激な減少、そしてその後の緩やかな回復が見られる。これは、タイプ II 事象によって駆動される円盤質量の蓄積、枯渇、および再構築のサイクルを示唆している。
4. ジャンプ間の安定性: 大きな振幅の位相シフトと長期的な回復期間（数百日）を示すシステム EXO 2030+375 とは異なり、GRO J1008−57 では、タイプ II 事象間の長期にわたる位相回復や連続的なドリフトの証拠は見られない。このシステムは、タイプ II アウトバーストが終了するとほぼ直ちに、新しい安定した位相状態に戻る。

意義と提案されるメカニズム
本論文は、GRO J1008−57 において観測されたステップ状の進化が、Kozai-Lidov（KL）振動メカニズム、全球的な片腕円盤振動、あるいは長期的な円盤歳差運動など、他の BeXRB に対して提案されている既存のモデルに挑戦するものであると主張する。

• KL 振動の棄却: GRO J1008−57 に対する推定 KL 振動の時間スケールは約 79 年であり、観測ベースラインを大幅に上回る。さらに、KL メカニズムは通常、大きな位相シフトと離心率の増大を予測するが、ここでは観測されていない。
• 歳差運動/振動の棄却: 長期間の位相回復期間（EXO 2030+375 で見られるような）の欠如は、長期的な円盤歳差運動または全球的な密度擾乱が主要な駆動因子であるという説を否定する。

提案される説明:
著者らは、GRO J1008−57 の独特な振る舞いは、その**長い軌道周期（249.48 日）**によって駆動されており、これにより中性子星と星周円盤との間の相互作用が稀になることが原因であると提案する。

• 円盤再構築サイクル: タイプ II アウトバーストは、Be 円盤を枯渇させる事象として解釈される。長い軌道周期により、システムは相互作用の間に十分な時間を有し、円盤が密度を完全に再構築し、アウトバースト前の状態と幾何学的に類似した構造を回復する。
• 小さなジャンプ: したがって、各タイプ II 事象の後、円盤はわずかに異なるが安定した構成に落ち着き、これに続くタイプ I アウトバーストの軌道位相に小さく急激なシフトを引き起こす。システムは、次の枯渇事象が発生するまで、この新しい安定状態に留まる。
• EXO 2030+375 との対比: 対照的に、EXO 2030+375 は短い軌道周期（$\sim46$日）を持ち、これにより頻繁な相互作用が生じ、円盤が完全に緩和するのを妨げる。その結果、擾乱が蓄積し、大きく連続的な位相シフトを駆動する。

結論
本研究は、タイプ I アウトバースト位相の時間的進化が普遍的なものではなく、システムの軌道力学に依存することを示すことで、BeXRB のアウトバーストメカニズムに関する理解を拡張する。GRO J1008−57 におけるステップ状の進化は、タイプ II アウトバーストによって駆動される円盤の枯渇と再構築サイクルの証拠を提供する。このサイクルにおいて、長い軌道周期が事象間で円盤が安定した幾何学を再確立することを可能にし、連続的なドリフトではなく離散的な位相ジャンプをもたらす。