[論文レビュー] Magnetic power spectra from Faraday rotation maps - REALMAF and its use on Hydra A

P. Kuchar, T. A. Enßlin|arXiv (Cornell University)|Dec 19, 2009

Geomagnetism and Paleomagnetism Studies参考文献 46被引用数 45

ひとこと要約

本稿では、球対称な磁場自己相関モデルと磁場の散発的制約のフーリエ空間的取り扱いを組み込んだ、Faraday回折マップから磁場パワースペクトルを再構成する新規の最大後確信度手法であるREALMAFを提案する。ハイドラAに適用した結果、0.3–8 kpcのスケールでコルモゴロフに類似したパワースペクトル（スペクトル指数 ≈ 1.8）が観測され、折りたたみがないこと、および45°のジェット投影角における中心磁場強度が36 μGに一致することが判明した。これは電子密度スケーリング（α = 1）と整合的である。

ABSTRACT

We develop a novel maximum-a-posteriori method to measure magnetic power spectra from Faraday rotation data and implement it in the REALMAF code. A sophisticated model for the magnetic autocorrelation in real space permits us to alleviate previously required simplifying assumptions in the processing. We also introduce a way to treat the divergence relation of the magnetic field with a multiplicative factor in Fourier space, with which we can model the magnetic autocorrelation as a spherically symmetric function. Applied to the dataset of Hydra A north, we find a power law power spectrum on spatial scales between 0.3 kpc and 8 kpc, with no visible turnover at large scales within this range and a spectral index consistent with a Kolmogorov-like power law regime. The magnetic field strength profile seems to follow the electron density profile with an index alpha=1. A variation of alpha from 0.5 to 1.5 would lead to a spectral index between 1.55 and 2.05. The extrapolated magnetic field strength in the cluster centre highly depends on the assumed projection angle of the jet. For an angle of 45 degree we derive extrapolated 36 muG in the centre and directly probed 16 muG at 50 kpc radius.

研究の動機と目的

• 従来の簡略化仮定に依存しない、Faraday回折データから磁場パワースペクトルを推定するための堅牢な手法の開発。
• 物理的制約をよりよく反映するため、磁場の自己相関を実空間で高次元の球対称関数でモデル化すること。
• フーリエ空間における乗法的要因を用いて磁場の散発的性質（∇·B = 0）を組み込むことで、スペクトル再構成の正確性を向上させること。
• 本手法をハイドラA銀河団に適用し、実宇宙環境における磁場構造とスペクトル的性質を調査すること。
• 仮定されたジェット投影角に応じて推定される中心磁場強度の感度をテストし、クラスターコアにおける推定値の不確実性を定量化すること。

提案手法

• パワースペクトルにスケール不変の事前分布を適用することで、情報のない正則化を保証するベイズ的最尤推定フレームワークを採用する。
• 磁場の実空間自己相関関数を、修正ベッセル関数に基づく形式でモデル化し、フーリエ空間における球対称表現を可能にする。
• 散発的制約（∇·B = 0）はフーリエ空間における乗法的要因を介して強制され、発散なしの磁場の一貫したモデル化が可能になる。
• 測定誤差を考慮した尤度モデルを用いて、回折測定マップからパワースペクトルを再構成し、磁場およびRMデータの両方にガウス事前分布を仮定する。
• スペクトル指数は、再構成されたスペクトルにパワーロー ε(k) = qk⁻ʳ をフィットさせることで決定され、不確実性は後騒のヘッセ行列から導出される。
• 全逆ヘッセ行列を投影して、スペクトルパラメータの共分散行列を抽出し、勾配と振幅の統計的推論を可能にする。

実験結果

リサーチクエスチョン

• RQ1 ベイズ的最尤推定法を用いることで、制限的な簡略化仮定に依存せずに、Faraday回折データから磁場パワースペクトルを再構成できるか？
• RQ2 球対称な磁場自己相関モデルを組み込むことで、従来の手法と比較してパワースペクトル再構成の正確性がどのように向上するか？
• RQ3 ハイドラA銀河団における磁場パワースペクトルのスペクトル指数は何か？ 0.3–8 kpcのスケールでコルモゴロフに類似した振る舞い（指数 ≈ 1.8）を示すか？
• RQ4 推定される中心磁場強度は、ハイドラAにおける仮定されたジェット投影角にどれほど敏感か？
• RQ5 ハイドラAにおける磁場強度プロファイルは、一部の磁場生成モデルが予測するように、指数α = 1のべき乗則に従うか？

主な発見

• ハイドラAにおける磁場パワースペクトルは、0.3 kpcから8 kpcのスケールでパワーロー的挙動を示し、スペクトル指数は約1.8であり、コルモゴロフに類似した状態と整合的である。
• 調査範囲内の大きなスケールでスペクトルの折り返し（turnover）は観測されず、最小分解スケールまで連続的なエネルギー分布が存在することを示唆している。
• 磁場強度プロファイルは、電子密度プロファイルにべき乗則指数α = 1に従い、αを0.5から1.5の範囲で変化させた場合、スペクトル指数は1.55から2.05の間で変動する。
• 45°のジェット投影角において、中心磁場強度の外挿値は36 μGであり、50 kpcの半径で直接測定された値は16 μGであった。
• 本手法は、フーリエ空間における乗法的要因を用いて磁場の散発的制約を効果的に組み込み、物理的整合性を向上させた。
• ベイズフレームワークによりヘッセ行列を用いた堅牢な誤差推定が可能となり、スペクトルパラメータおよびその不確実性の信頼性のある推論が可能になった。