[論文レビュー] Cyclotron resonant scattering feature simulations

Fritz-Walter Schwarm, Ralf Ballhausen+13 more|arXiv (Cornell University)|Jan 26, 2017

Astrophysical Phenomena and Observations参考文献 46被引用数 13

ひとこと要約

本論文は、X線パulsarにおける合成サイクロトロン共鳴散乱特徴（CRSFs）を生成するための新しいモンテカルロシミュレーションコードを提示する。これにより、複雑な円対称な降着幾何構造の正確なモデリングが可能になる。この手法は、事前に計算された平均自由行程（MFP）テーブルとグリーン関数を用い、ラインプロファイルを効率的にシミュレートする。その後、XSPEC互換のモデル（cyclofs）を用いてNuSTARデータにフィットさせ、Cep X-4の磁場強度は3.6×10¹² Gであると導出される。これはガウス線フィットとは顕著に異なる。物理的に現実的なCRSF形状が得られるためである。

ABSTRACT

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green’s functions, which can be used later to generate spectra for arbitrary continua. Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green’s function tables are available online (see link in footnote, page 1).

研究の動機と目的

• サイクロトロン共鳴散乱特徴（CRSFs）のモンテカルロシミュレーションにおける従来の制限、特に単純で事前に定義された幾何構造に限定されていた点を克服すること。
• 任意に複雑な円対称な降着柱配置をモデル化できる柔軟性と並列処理が可能なシミュレーションコードの開発。
• モノエネルギーフォトン注入に対するグリーン関数の生成により、任意の連続スペクトルモデルと効率的なスペクトル適合が可能になる。
• 合成ラインプロファイルと観測データを直接比較できるXSPEC互換モデル（cyclofs）の提供。
• 古典的幾何構造および先行研究との妥当性検証により、今後の物理的解釈に向けた一貫性と信頼性を確保すること。

提案手法

• 本シミュレーションは、同伴論文（Schwarm et al. 2017, Paper I）で得られた平均自由行程（MFP）テーブルを用い、熱平均散乱断面積の高速補間を実現する。
• 並列化されたモンテカルロアルゴリズムにより、ユーザー定義の磁場、温度、密度勾配を有する3次元円柱幾何構造内でのフォトンの伝播、散乱、エネルギー損失をシミュレートする。
• モノエネルギーフォトンを注入し、外部スペクトル応答を記録することでグリーン関数を生成し、時間のかかるシミュレーションと連続スペクトルの選択を分離する。
• XSPECモデル「cyclofs」は、事前に計算されたグリーン関数テーブルを用いて実装され、合成CRSFプロファイルを観測スペクトルに直接適合可能となる。
• 相対論的効果（重力赤方偏移など）を考慮し、視線角依存性を組み込み、ライン形状と幅に顕著な影響を与える。
• 複数の視線角の平均化を近似的な手法で実装し、パルス位相平均データに対応可能となる。

実験結果

リサーチクエスチョン

• RQ1 複雑で非一様な降着幾何構造は、X線パulsarにおけるサイクロトロン共鳴散乱特徴の形状と深さにどのように影響を与えるか？
• RQ2 視線角と温度勾配は、観測されたライン幅と中心エネルギーにどの程度の影響を及ぼすか？
• RQ3 物理的根拠に基づく非ガウス型ラインプロファイルモデルは、経験的ガウス吸収モデルと比較して、磁場推定の精度を向上させ得るか？
• RQ4 現実的な散乱断面積と伝播効果を組み込むことで、観測フィットから導かれる物理的パラメータはどのように変化するか？
• RQ5 物理的に整合性のあるライン形状モデルを用いる場合、連続スペクトルモデルの選択が、導出されるCRSFパラメータにどのような影響を及えるか？

主な発見

• 新規シミュレーションコードは、従来の幾何制限に縛られることなく、複雑な円対称幾何構造における合成サイクロトロンラインプロファイルの生成に成功した。
• コードは視線角依存性が強く、磁場軸からの角度に応じてライン幅と深さが顕著に変化することが示された。
• Cep X-4のNuSTARデータにフィットした結果、物理的XSPECモデル「cyclofs」により磁場強度が3.6×10¹² Gと導出された。これはガウス吸収ラインモデルで得られる値とは顕著に異なる。
• 最良のフィットでは、自由度875のとき赤色カイ二乗（χ²）が1.17となり、ガウスモデル（χ² = 1.57）を上回り、固定中心エネルギーを仮定した場合に明確な残差が観測された。
• 基本的なサイクロトロンライン形状は本質的に複雑で非対称的であり、非ガウス型プロファイルを用いることで、滑らかで対称的に見えるラインであっても磁場推定に系統的差が生じることが明らかになった。
• XSPECモデル「cyclofs」と関連するグリーン関数テーブルは公開されており、広範な利用と標準X線スペクトル適合ワークフローへの統合が可能である。