Discovery of a New Spectral Transition in Swift J0243.6+6124 in the Sub-Eddington Regime

技術要約：Swift J0243.6+6124における新たなスペクトル遷移の発見

問題提起
降着X線パルサー（AXRP）は、質量降着率の変化によって駆動される複雑なスペクトルおよびタイミングの進化を示す。理論的枠組みは、クーロン光度（$L_{\text{Coul}}$）や臨界光度（$L_{\text{crit}}$）といった特徴的な遷移光度を予測しているが、ほとんどの天体ではそのような遷移は1つまたは2つしか観測されない。Swift J0243.6+6124は、最初の銀河系超高光度X線パルサーであり、6桁（$10^{34}$–$10^{39}$ erg s$^{-1}$）に及ぶ降着光度を持つユニークな実験場を提供している。これまでの研究で本天体におけるいくつかの遷移が特定されているが、狭帯域（NICER）と広帯域（Insight-HXMT）のスペクトル解析の間には不一致が存在し、サブエディントン領域における複数の遷移の物理的起源は依然として不明である。本論文は、サブエディントン領域における新たなスペクトル遷移の特定に取り組み、この天体の複雑な挙動を説明するために必要な基礎となる磁場構成について調査するものである。

方法論
著者らは、3つのアウトバースト（2018年、2019年、2023年）の減衰フェーズにおけるInsight-HXMTおよびNICER望遠鏡のデータを用いて、Swift J0243.6+6124の詳細なスペクトル解析を行った。

• データ選択: サブエディントン領域（$L_X \sim 7.9 \times 10^{36}$ to $1.2 \times 10^{38}$ erg s$^{-1}$）をカバーするように観測データを選択した。解析には、最新のGaia DR3による距離 $5.2 \pm 0.3$ kpc を採用した。
• スペクトルモデリング: スペクトルはXSPECを用い、constant*tbabs*(bbodyrad+cutoffpl) モデルを使用してモデリングを行った。本研究では、3つのデータ構成（Insight-HXMTのみ（2–100 keV）、NICERのみ（0.7–10 keV）、および結合フィット（0.7–100 keV））を比較した。
• 統計解析: パラメータの不確実性はマルコフ連鎖モンテカルロ（MCMC）サンプリングを介して推定した。スペクトルの進化を特徴付けるため、ブラックボディ半径と光度の関係を単一線形モデルおよび折れ線モデルの両方を用いてフィッティングした。モデル比較は、$\Delta$AICc、$\Delta$BIC、および10,000回のブートストラップ反復を用いた尤度比検定（LRT）を用いて行った。
• 物理的解釈: 導出された遷移光度は、クーロン減速に関する理論的予測と比較された。著者らは、複数の遷移の存在と、弱磁場および強磁場の推定値間の不一致を解決するために、多重極磁場モデル（双極子＋四重極子）を提案した。

主な結果

1. 新たな遷移（$L_t$）の発見: 解析の結果、遷移光度 $L_t \approx 4.46^{+0.53}_{-0.53} \times 10^{37}$ erg s$^{-1}$ において、スペクトルパラメータの進化における明確な破れが明らかになった。これは、Swift J0243.6+6124において特定された5番目の遷移光度である。
2. スペクトル進化のシグネチャ: $L_t$ において、ブラックボディ半径と光度の相関が正から負へと逆転する。具体的には、ブラックボディの規格化（半径）は顕著なターンオーバーを示し、$L_t$ 以下では光度とともに増加し、$L_t$ 以上では減少する。同時に、光子指数はほぼ一定の状態から光度に対して正の相関を持つよう進化し、カットオフエネルギー（$E_{\text{cut}}$）の進化は鈍化する。
3. 装置間の不一致: 本研究は、NICERのみのフィット（狭帯域）が非物理的な結果（低温度 $kT \sim 0.4$ keV、大半径 $R \sim 20$ km の解へのドリフト）をもたらすことを強調している。対照的に、Insight-HXMTおよび結合フィットは、極キャップのホットスポットと一致する高温（$kT \sim 3$ keV）かつ小半径（$R \sim 1$ km）のブラックボディ成分を制約している。著者らは、狭帯域のフィッティングのみでは、連続成分の推論に固有のバイアスが生じることを結論付けている。
4. 磁場構成: 観測された $L_t$ は、典型的な中性子星のパラメータに対する標準的なクーロン光度の予測と矛盾しており、これを満たすには極めて低い磁場（$B \sim 10^{10}$ G）が必要となる。代わりに、著者らは多重極磁場構成を提案している。
   • 弱磁極: $L_t$ は、表面磁場強度 $B_S \sim 2.8 \times 10^{12}$ G を持つ「弱い」磁極の臨界光度と解釈される。
   • 強磁極: 以前に検出された 146 keV のサイクロトロン共鳴散乱特徴（CRSF）は、$B_N \sim 1.6 \times 10^{13}$ G の「強い」磁極を示唆している。
   • 幾何学: 双極子（$6.6 \times 10^{12}$ G）と四重極子（$9.4 \times 10^{12}$ G）の軸対称な重ね合わせによる磁場は、この非対称性を再現できる。この構成は、なぜ異なる極が異なる光度で臨界領域に入るのか、そしてなぜ複数のスペクトル遷移が発生するのかを説明する。

意義と主張
本論文は、$L_t$ の発見が Swift J0243.6+6124 における降着物理の複雑さを大幅に増大させ、同天体を最も多くの遷移が特定された天体にしたと主張している。著者らは、5つの異なる遷移の存在は、通常 $L_{\text{Coul}}$ と $L_{\text{crit}}$ のみを予測する標準的な双極子磁場モデルでは説明できないと断言している。

本研究の主要な意義は、多重極磁場がこの天体の放射を支配する基礎的なメカニズムであるという提案にある。本研究は、双極子磁場と四重極子磁場の組み合わせが非対称な極キャップを作り出し、各極において異なる臨界光度を許容することを実証している。このモデルは、低磁場遷移（$L_t$）と高磁場CRSF検出、およびより高光度の遷移（$L_4, L_5$）との整合性をとることに成功している。さらに、本論文は、AXRPのスペクトルパラメータを正しく制約するためには広帯域の観測が不可ло欠であることを強調しており、狭帯域データは誤った物理的解釈を招く可能性があることを示している。これらの知見は、磁気圏スケールでの有効磁場がより弱く見える場合でも、中性子星の局所的な磁場は $10^{13}$ G を超え得ることを示唆しており、高度に磁化された系における降着柱の理解に新たな枠組みを提供している。