Joint Multi-Period Fermi-LAT and LHAASO Constraints on Axion-Like Particles from Mrk 421 Using Profile Likelihood with Gaussian Copula Correlation

技術要約：Mrk 421からのFermi-LATおよびLHAASOを用いたマルチ・ピリオド結合によるアキシオン様粒子の制約

問題提起
アキシオン様粒子（ALP）は、標準模型の拡張として十分に動機付けられた存在であり、外部磁場において光子とALPの振動を引き起こし、遠方の天源からの高エネルギーガンマ線のスペクトルを変化させる可能性がある。これまでの研究ではALPパラメータの制約がなされてきたが、依然として広範なパラメータ空間が未制約のままである。このような探索における大きな課題は、ガンマ線を減衰させる銀河外背景光（EBL）に関連する系統誤差の扱いである。従来の多くのアプローチでは、EBLによる吸収と光子-ALP振動を独立した乗法的な因子として扱い、それらの間の非自明な相互作用を無視している。さらに、マルチ・エポック解析においては、観測期間を統計的に独立したものと仮定することが多いが、これは共有された光源の変動性、計器の系統誤差、およびモデルの仮定から生じる微弱な残留相関を無視している。これらの簡略化は、バイアスのある制約や、人工的に狭められた信頼区間を導く可能性がある。加えて、ブラザルであるMrk 421は、GeV（Fermi-LAT）とTeV（LHAASO）の放射の間に弱い相関を示しており、これが異なる物理的起源を示唆しているため、広帯域データの解釈を複雑にしている。

手法
著者らは、2021年3月から2024年3月までのMrk 421の5つの同時観測エポックを用い、Fermi-LATとLHAもLHAASOのデータを組み合わせることで、ALPパラメータ（$m_a$ および $g_{a\gamma}$）を制約するための包括的なフレームワークを提示している。

1. 伝搬定式化: 因子化近似とは異なり、本研究では、EBL吸収を吸収項（$i/2\lambda_\gamma$）として光子-ALP混合ハミルトニアンに直接組み込んだ、自己整合的な光子伝搬フレームワークを採用している。EBL強度は、不確実性を考慮するために、ガウス事前分布（$\mu=1.0, \sigma=0.2$）を持つ正規化因子（$K_{EBL}$）を伴うナンシング・パラメータとして扱われる。
2. 光源モデリング: Mrk 421における弱いGeV-TeV相関に対処するため、2つの異なるジェット放射モデルをテストしている：
   • コンパクトな内部領域（フレアリングVHE放射）と拡張された外部領域（静穏な放射）を特徴とする二領域ハイブリッドモデル。
   • 弱い相関と異なる変動振幅を考慮すると物理的に妥当な、単一領域ハドロンモデル。
     両モデルは、ヘリカルおよび乱れた磁場成分を利用し、光源、銀河間物質（IGMF）、および銀河磁場（GMF）を通過する光子を伝搬させる。
3. 統計的フレームワーク: 著者らは、結合プロファイル尤度解析を開発した。エポックを独立したものとして扱う限界に対処するため、保守的な相関係数（$\rho = 0.03$）を持つガウス・コピュラ形式を導入している。これにより、変換された $\Delta\chi^2$ 過剰分布を通じて、エポック間の微弱な相関をモデル化している。解析では、ALPのパラメータをプロファイリングする前に、ALPなしの仮説の下でナンシング・パラメータ（固有のスペクトルパラメータおよび $K_{EBL}$）をグローバルに最適化する。
4. データ処理: 光子の生存確率は、gammaALPs コードを用いて、4次元グリッド（エネルギー、ALPパラメータ、EBL正規化）上で事前計算され、尤度評価中に補間される。固有のスペクトルモデルは、有限サンプルサイズを補正した赤池情報量基準（AICc）を用いて、エポックごとに選択される。

主要な貢献

• 統一伝搬モデル: 本論文は、EBL吸収と光子-ALP混合を伝搬方程式内で自己整合的に扱う伝搬モデルを実装しており、乗法的な因子化によるバイアスを回避している。
• 相関マルチ・エポック解析: 観測期間間の微弱な相関を扱うためにガウス・コピュラを導入することで、複数のエポックのデータセットを組み合わせるための統計的に一貫した手法を提供し、人工的な信頼区間の狭まりを防いでいる。
• 広帯域制約: 本研究は、Fermi-LAT（GeV）とLHAASO（マルチTeVからPeV）のユニークな広帯域カバレッジを活用し、広い質量範囲（0.1–500 neV）にわたってALPパラメータを探索している。
• 堅牢なモデリング: 二領域モデルと単一領域ハドロンシナリオの両方をテストすることで、Mrk 421の複雑な放射物理を考慮に入れ、制約が単一の放射仮定に過度に依存しないようにしている。

結果
結合プロファイル尤度フレームワークを用いて、著者らは光子-ALP結合（$g_{a\gamma}$）の95%信頼区間（CL）上限値を導出した：

• 低質量領域 ($m_a \lesssim 1$ neV): 制約は非常に厳しく、下限値は $5.0 \times 10^{-14}$ GeV$^{-1}$、上限値は $8.33 \times 10^{-14}$ GeV$^{-1}$ である。この領域はCASTの制約と感度が重なる。
• 中間から高質量 ($0.1 - 500$ neV):
  • 二領域ハイブリッドモデルの場合、保守的な上限値は $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$ である。
  • 単一領域ハドロンモデルの場合、保守的な上限値は $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$ である。
• 挙動: 除外輪郭は、特徴的な「ジャンプ・プラトー」構造を示す。$m_a \approx 1$ neV 付近で急激な遷移が起こり、上限値は約1桁増加した後、高質量側に向かって緩やかに緩和していく。下限値は全質量範囲にわたって安定している。
• 堅牢性: 代替の相関係数（$\rho = 0, 0.15, 0.3$）を用いたテストでは、制約の変化が1桁未満に収まっており、結果の安定性が確認された。

意義
本論文は、これらの結果がLHAASOの広帯域VHEカバレッジとMrk 421の長期モニタリングから大きな恩恵を受けていると主張している。開発されたフレームワークは、次世代のガンマ線観測装置による将来のマルチメッセンジャーALP探索に適した、堅牢でスケーラブルなパイプラインを提供する。微弱な相関を明示的にモデル化し、EBL吸収を自己整合的に扱うことで、本研究は、従来の因子化アプローチにおける系統的バイアスに対処し、より信頼性が高く保守的なALPパラメータ空間の制約を提供している。本研究は、Mrk 421における弱いGeV-TeV相関が、正確なALP制約のために複雑な放射モデル（二領域またはハドロンモデル）を必要とすることを実証しており、提案された統計的手法がそれに伴う不確実性を効果的に処理できることを示している。