Line-of-sight magnetic-field propagation effects on axion-like particle constraints from GRB 221009A

技術要約：GRB 221009Aからのアキシオン様粒子制約における視線方向の磁場伝搬効果

問題提起
LHAASOによるGRB 221009Aからの、エネルギーが最大18 TeVに達する5,000個以上の超高エネルギー（VHE）光子の検出は、標準的な天体物理学的期待との間に緊張関係を生じさせている。従来のモデルでは、極限外背景光（EBL）による対生成（$\gamma\gamma \to e^+e^-$）を通じて、これらの光子が強く減衰すると予測されている。この緊張は、EBLモデリングの系統的な不確実性や限界に起因する場合もあるが、同時にアキシオン様粒子（ALP）を探索する動機ともなっている。ALPは、ホスト銀河内で光子へと変換され、宇宙空間をALPとして通過することで（EBLによる減衰を回避し）、その後、天の川銀河内で再変換されることで、宇宙の透明性を高めることができる。しかし、この光子-ALP振動の効率は、視線方向の磁場環境に極めて敏感である。GRB 221009Aを用いてALPパラメータ（$m_a$, $g_{a\gamma\gamma}$）を制約した先行研究では、ホスト銀河磁場（HGMF）、銀河間磁場（IGMF）、および天の川銀河磁場（MWMF）の不確実性が、導出される排除限界にどのように影響するかを体系的に検討していない。

手法
著者らは、GRB 221009AのLHAASOスペクトル（具体的には、0.23–11.61 TeVをカバーする$T_0 + 230$ sから$T_0 + 300$ sまでの時間間隔）を用いて、ALP制約を再検討している。解析には、厳密な伝搬フレームワークを用いている：

1. 定式化： ビームの進化は、ALP質量、プラズマ振動数、QED真空偏極、CMB分散、および光子吸収を考慮した混合行列を含む、シュレディンガー型の状態ベクトルに対するシュレディンガー型方程式によって記述される。複数の磁化ドメインを通じた進化は、転送行列の順序付き積を介して計算される。
2. 磁場モデル：
   • ホスト銀河（HGMF）： 簡略化された均一磁場、天の川銀河を模した現象論的なプロキシ、およびHSTの観測とNGC 891に似たシミュレーションに基づく現実的なモデル（モデル3）の3つのモデルをテストしている。
   • 天の川銀河（MWMF）： Jansson–Farrar（JF12, JF12b, JF12c）モデルと、より新しいUnger–Farrar（UF）モデルを比較し、UF「ベース」モデルをベンチマークとして選択している。
   • 銀河間空間（IGMF）： IGMFは、空隙に点在するランダムな方位を持つ磁化ドメインの連鎖としてモデル化されている。DLSHE（Domain-Like Sharp-Edges）モデルとDLSME（Domain-Like Smooth-Edges）モデルの2つの実現形態を検討している。研究では、ドメイン長（$L_{dom}$）、充填率（$f$）、および磁場強度（$B_{IG}$）を変化させている。
3. 統計解析： 観測されたスペクトルに対して95%信頼区間での最小$\chi^2$適合を用いて制約を導出している。固有のGRBスペクトルは、対数放物線関数（パラメータ $J_0, a, b$）としてモデル化されている。光子の生存確率（$P_{ALP}$）は、様々なALPパラメータと磁場構成に対して計算され、$m_a$の関数としての$g_{a\gamma\gamma}$の上限値を決定する。

主要な貢献と結果
本論文は、ALP排除等高線における磁場不確実性の影響を、$m_a$–$g_{a\gamma\gamma}$平面において体系的に定量化している：

• ホスト銀河および銀河磁場への依存性： HGMFおよびMWMFモデルの選択は、排除限界に中程度の変動しか及ぼさない。低質量領域（$m_a \lesssim 3 \times 10^{-8}$ eV）では、排除限界は質量に依存せず、異なるモデル間で約1.5倍程度の変動にとどまる。モデル3（HGGF）およびUFベースモデル（MWMF）が、堅牢なベンチマークとして特定された。
• IGMF不確実性の支配性： IGMFが支配的な天体物理学的不確実性である。IGMFパラメータを変化させると、排除等高線に顕著な変化が生じ、明確な振動構造が現れる。
• 3つの質量領域： IGMFの導入により、排除限界に3つの明確な領域が現れることが明らかになった：
  1. 大質量領域 ($m_a \gtrsim 10^{-8}$ eV)： 振動長が短いため（$l_{osc} \ll L_{dom}$）、IGMFにおける光子-ALP変換は非効率的である。制約は、$B_{IG}=0$の場合の結果と密接に一致する。
  2. 中間質量領域 ($10^{-10} \text{ eV} \lesssim m_a \lesssim 10^{-8}$ eV)： 振動長がドメインサイズと同程度になる（$l_{osc} \sim L_{dom}$）。マルチドメイン伝搬は、位相 $\omega = \pi L_{dom}/l_{osc}$ の正弦関数と余弦関数の積を含む項を生成する。これにより、排除等高線には、$m_a$と$L_{dom}$に依存した特徴的な振動構造が生じる。
  3. 小質量領域 ($m_a \lesssim 10^{-10}$ eV)： 振動長が長くなり、質量にほとんど依存しなくなる（$l_{osc} \gg L_{dom}$）。振動的な特徴は消失するが、IGMFは依然として、ALPsがEBL減衰を回避することを可能にすることで、高エネルギー光子の生存確率を大幅に高める。その結果、排除限界は$B_{IG}=0$の場合よりも強くなり（より低い$g_{a\gamma\gamma}$）、質量に対する滑らかな依存性を示す。
• ドメイン長への感度： 振動領域の発生は、IGMFのコヒーレンススケール（$L_{dom}$）に直接結びついている。$L_{dom}$を大きくすると、遷移がより高い質量へとシフトする。

意義
本論文は、GRB 221009Aからの堅牢なALP制約を導出するためには、視線方向の磁場、特に不透明な部分が多いIGMFの現実的なモデリングが不可欠であると結論付けている。ホスト銀河および銀河磁場の不確実性は副次的であるが、IGMFの確率的な性質（磁場強度、コヒーレンススケール、および充填率）は、導出される限界を劇的に変え得る。排除等高線に現れる振動的な特徴は、IGMFの確率的な性質の可視化された痕跡である。著者らは、将来のガンマ線によるALP探索において、天体物理学的な不確実性を新物理学と誤認したり、あるいはその逆が生じたりすることを避けるために、これらの伝搬効果を考慮する必要があると強調している。