Simulating the LOcal Web (SLOW): VII. Intergalactic magnetic field models for multi-messenger applications

技術要約：Simulating the LOcal Web (SLOW) VII. マルチメッセンジャー・アプリケーションのための銀河間磁場モデル

問題提起
極高エネルギー宇宙線（UHECR）および極高エネルギーガンマ線の伝播は、銀河間磁場（IGMF）の影響を強く受ける。マルチメッセンジャー観測の解釈および磁場生成シナリオの制約には、宇宙の大規模構造にわたるIGMFを正確に特徴付けることが不可欠である。しかし、マルチメッセンジャー研究のための既存の磁場モデルには、主に以下の2つの制限がある：

1. 制約のない初期条件： 多くのシミュレーションは観測された局所的な大規模構造を再現しておらず、特定の天体ソースに対する視線方向のプロファイルが不正確になっている。
2. スケールと忠実度の問題： 制約付きのシミュレーションであっても、ミラーリングやボリュームの連結といった人工的な手法を用いなければ、関連する天体（ブレーザーなど）までの距離をカバーするには規模が小さすぎる場合があり、それらの手法は磁場プロファイルの忠実度を低下させる。さらに、多くの研究は物理的なシミュレーションから抽出されたモデルではなく、簡略化された表現（例：ガウス型ランダム場）に依存している。加えて、個々の銀レット質ハロー（潜在的なUHECR/TeVソース）のシミュレーション上の対応物を特定することは困難である。なぜなら、初期条件は線形再構成されており、約3〜4 Mpc以上のスケールのみを制約しているためである。

手法
著者らは、SLOW (Simulating the LOcal Web) シリーズの制約付き宇宙論的シミュレーションから導出された磁場モデルを提示しており、具体的には以下の2つのコンポーネントを利用している：

• SLOW-AGN15363： 活動銀河核（AGN）および超巨大ブラックホール（SMBH）を含む高解像度シミュレーションであり、銀河や銀河団の位置を特定できるが、磁場は含まれていない。
• SLOW-CR30723： 非放射MHD（磁気流体力学）シミュレーションであり、より高い解像度を持ち、磁場と宇宙線進化を含んでいるが、AGNフィードバックは欠いている。
  これらのシミュレーションは、Cosmicflows-2 カタログ（特異速度）に基づく初期条件を共有しており、ハローの相互照合が可能である。

本研究では、以下の6つの磁場モデルを分析している：

1. $B_{sim}$：MHDシミュレーションからの生の磁場（$10^{-14}$ Gの一様な原始磁場によってシードされたもの）。
2. 再スケーリング・モデル： フィラメントにおける解像度の限界に対処するために、物理的仮定に基づいて再スケーリングされた5つのモデル：
   • $B_{\beta}$：一定のプラズマ-$\beta$（熱圧と磁気圧の比）。
   • $B_{F}$：乱流圧との平衡。
   • $B_{ff}$：理想的なMHDの磁束凍結（電子密度）によるスケーリング。
   • $B_{dyn, \downarrow}$ および $B_{dyn, \uparrow}$：乱流ダイナモ増幅シナリオに基づくスケーリング（下限および上限）。

新規アルゴリズム：「理想的な位置（Ideal Position）」の決定
初期条件が個々の銀レット質ハローを制約できないという問題に対処するため、著者らは、特定の観測ソースに対する「理想的な位置」($\mathbf{x}_{ideal}$) を決定するための新しいアルゴリズムを導入した：

1. 三角形分割： アルゴリズムは、ソースに最も近い3つの観測された銀河団を特定する。
2. ファジーな三角形分割： 観測された相対位置ベクトルをシミュレーション上の銀河団の位置に加えることで、「シフトされた重心」を計算する。
3. 最適化： 重心の周囲のグリッド上で、$R$（観測された距離とシミュレーション上の距離のガウス重み付き差分の和）を最小化する。この「ファジー」なアプローチにより、シミュレーションと観測された構造間の厳密な一致を強制するのではなく、大規模なドリフトを考慮したグローバルな最小値を求めることができる。
4. 不確実性の推定： 著者らは、理想的な位置からオフセットさせたランダムな位置への視線方向に沿って磁場プロファイルを計算することで、不確実性を推定している（2D天球面シフト、3D空間シフト、およびSMBHを宿す現実的なサブハローの選択を使用）。

主な結果

• 磁場特性： モデルは幅広い充填率（filling factor）を示す。すべてのモデルは銀河団のコアにおいて収束する（Coma銀河団の観測と一致する）が、フィラメントやボイドにおいては大きく乖離する。生のシミュレーション（$B_{sim}$）は、再スケーリングされたモデルと比較してフィラメントにおける磁場強度が著しく低く、しばしば0.1 nGを下回る。
• ガンマ線カスケード解析： モンテカルロ・コード ELMAG3.03 を用い、Mrk 501の「理想的な位置」における電磁ガンマ線カスケードをシミュレートした。
  • ボイドにおける充填率が高いモデル（$B_{dyn, \uparrow}$, $B_F$, $B_{ff}$, $B_{\beta}$）は、カスケード粒子を完全に偏向させ、GeV放射を抑制する。
  • 生のシミュレーション（$B_{sim}$）および下限ダイナモモデル（$B_{dyn, \downarrow}$）は、有意な二次放射を許容し、これはガンマ線カスケード観測から得られた下限値（例：$B > 1.3 \times 10^{-17}$ G）と一致している。
  • 結果は、$B_{sim}$ が電磁ガンマ線カスケードの制約から導かれたIGMFを最もよく再現していることを示唆しており、これは、粒子を偏向させるために必要な充填率が、いくつかの再スケーリングモデルが示唆するほど高くはないことを意味している。
• 不確実性の低減： 「理想的な位置」法は、ランダムな視線方向と比較して、より低い分散（カスケンドスペクトルの平均四分位距離で測定）を持つ磁場プロファイルをもたらす。具体的には、現実的な銀河ハロー（「サブハロー」法）で終わる視線方向を使用することが、最も正確で不確実性の低いカスケンドスペクトルを提供した。

意義と主張
本論文は、制約付き宇宙論的シミュレーションから導出された、物理的に動機付けられた正確なIGMFモデルをアストロパーティクル物理学コミュニティに提供することを主張している。

• 手法的な進展： 著者らは、初期条件の制約スケール以下の銀河に対する「理想的な位置」を特定するアルゴリズムを初めて提案し、これにより特定のソースに対する正確な視線方向の磁場抽出を可能にしたと主張している。
• 堅牢性： 本研究は、制約付きの「理想的な位置」モデルを使用することで、ランダムまたは制約のない磁場プロレットを使用する場合と比較して、ガンマ線カスケードシミュレーションにおける不確実性が減少することを実証している。
• 将来の有用性： アルゴリズムが利用している、シミュレーション上のクラスターと実際のクラスター間の大規模構造のドリフトは、将来の作業においてシミュレーションの制約付き初期条件を精緻化するための潜在的なルートとして示唆されている。
• 制限事項： 著者らは、現在のモデルが一様な原始磁場シードに基づいていることを謙虚に認めている。IGMFのトポロジー（ガンマ線ハローなどの観測量に極めて重要）を完全に捉えるためには、将来の研究において、天体物理学的なシード生成シナリオや、ヘリシティを伴う（あるいは伴わない）確率的な原始シードを探索する必要があると述べている。現在のモデルは、局所宇宙（ボックスサイズ 500 $h^{-1}$ Mpc）に限定されており、すべての磁場生成シナリオをまだカバーしていない。