Cyclotron lines in subcritical X-ray pulsars: Monte Carlo simulations reveal the origin of the observed variability

技術的サマリー：臨界未満の X 線パルサーにおけるサイクロトロン線

問題提起
X 線パルサー（XRP）は、そのスペクトルにサイクロトロン共鳴散乱特徴（CRSF）を示し、これは中性子星（NS）の磁場強度の直接的な診断手段となる。観測により、CRSF の重心エネルギー（$E_{CRSF}$）および線幅（$\sigma_{CRSF}$）が X 線光度（$L$）およびパルス位相の関数として有意に変動することが明らかになっている。臨界未満の光度領域（$L \lesssim 10^{37}$ erg s$^{-1}$）では、$E_{CRSF}$ と $L$ の間に正の相関がしばしば観測されるが、これは臨界超過源で見られる負の相関と対照的である。定性的なモデルは、この現象が降着漏斗内の落下プラズマから散乱される光子のドップラーシフトに起因すると示唆しているが、この変動に対する定量的かつ厳密な放射輸送処理は欠如していた。具体的には、共鳴散乱、放射圧による流れ速度への影響、および観測幾何学との相互作用は、観測された傾向を説明するために詳細なシミュレーションを必要とする。

手法
著者らは、NS 高温スポット上の降着漏斗内で詳細な放射輸送計算を実行するための相対論的モンテカルロ（MC）コードを開発した。このモデルは、臨界未満の源に対する Mushtukov ら（2015b）が提案した物理的枠組みを採用している：

• 幾何学： 磁気軸（$B = B\hat{z}$）に沿った半無限の円筒状降着チャネル。
• プラズマプロファイル： 電子密度（$n_e$）およびバルク速度（$\beta$）の解析的式は、質量および運動量保存則から導出され、放射圧による落下プラズマの減速を組み込んでいる。速度プロファイルは、臨界光度 $L^*$ に対する光度に依存する。
• 放射輸送： 黒体スペクトルを持つシード光子が高温スポットから等方的に注入される。光子パケットはプラズマ中を伝播し、コンプトン散乱を受ける。
• 散乱物理： コードは、$b \ll 1$ および $E \ll m_ec^2$ に対して有効な、近似された偏光平均共鳴磁気散乱断面積（Loudas ら 2021）を採用する。散乱は電子静止系で処理され、プラズマのバルク運動を考慮してローレンツ変換が適用される。
• 境界条件： 2 つの方式がテストされた。「反射」（基底での幾何学的反射、エネルギー保存）と「熱化」（吸収と再放射）である。主要な結果は反射方式を用いている。
• パラメータ： シミュレーションは、代表的な NS（$M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G）に対して実行され、特定の XRP GX 304−1 のケースに適用された。降着光度は $L/L^* \in (0, 1)$ の範囲で変えられた。

主要な貢献と結果
本研究は、Mushtukov ら（2015b）のモデルに対する最初の定量的 MC 検証を提供し、以下の結果をもたらした：

1. スペクトル形態： 出現スペクトルは、広範な青翼を伴う顕著で非対称な CRSF を示す。CRSF は局所的な古典的サイクロトロンエネルギー（$E_c$）に対して系統的に赤方偏移している。この赤方偏移は、実験室系と落下電子系との間のドップラー効果に起因する。
2. 光度依存性： 全ての観測角度において、$E_{CRSF}$ および $\sigma_{CRSF}$ は光度と正の相関を示す。$L$ が増加すると、放射圧が落下流れをより効果的に減速させ、バルク速度（$\beta$）を減少させる。この速度の減少はドップラー赤方偏移を減じ、$E_{CRSF}$ を $E_c$ に近づける。同時に、より高い光度における光学的深さの増加は、線を深くし、青翼を広げる。
3. 観測角度依存性： CRSF の特性は、観測角度 $\theta$（磁気軸に対する）によって著しく変化する。
   • $E_{CRSF}$ は $\cos \theta$ が増加するにつれて減少する（つまり、視界がより「正面」になるにつれて）。これは、光子運動量と流れ速度の間の角度に依存するドップラーシフトによるものである。
   • $\sigma_{CRSF}$ は $\cos \theta$ とともに増加する。正面からの視界は視線方向に沿ったより広範な速度勾配をサンプリングし、より大きなドップラー広がりをもたらす。
   • したがって、このモデルはパルスサイクル中（位相分解変動）に $E_{CRSF}$ と $\sigma_{CRSF}$ の間の負の相関を予測するのに対し、位相平均スペクトルは光度と正の相関を示す。
4. GX 304−1 への適用： GX 304−1 の観測データに適用した際、このモデルはほぼ 1 桁の範囲にわたる $E_{CRSF}$ および $\sigma_{CRSF}$ の光度に対する観測された正の相関を成功裡に再現した。最良の適合は、降着漏斗が主に側面から観測される角度（$\cos \theta \approx 0.0-0.1$）で得られ、これは特定の幾何学的配置（例えば、高傾斜角におけるほぼ整列したローテーター、または低傾斜角におけるほぼ直交するローテーター）と一致する。

意義と主張
本論文は、磁気漏斗内の降着プラズマを伝播する放射の共鳴散乱が、臨界未満の XRP における観測された CRSF 変動に対する自然かつ定量的な説明を提供すると主張している。著者らは、その結果が以下の点を強調していると述べている：

• 降着チャネルが臨界未満の源における線形成領域であるという仮説を検証すること。
• バルク運動およびそれに伴うドップラーシフトが、衝突のない衝撃波（CS）を主要な線形成サイトとして必要とすることなく、光度依存変動の主要な駆動因子であることを示すこと。
• 位相分解データを解釈する際の観測幾何学の決定的な役割を強調し、位相平均傾向とは異なる線エネルギーと線幅の間の特定の負の相関を予測すること。
• 彼らのシナリオを CS モデル（線形成が衝撃波と表面の間で起こる）と区別すること。CS シナリオは、彼らのシミュレーションで観測された強い位相依存変動を予測しない点に注目すること。

著者らは、単純化された円筒幾何学の使用、1 次元解析プロファイルの使用、および一般相対論的効果や偏光の無視といった限界を指摘している。しかし、彼らはこれらの単純化が変動の起源に関する根本的な結論を変更しないとの見解を示している。この研究は、降着漏斗における共鳴散乱を用いて臨界未満の XRP スペクトルをモデル化するための概念実証として機能する。