論文要約：GRAPES-3 ミューオン望遠鏡による雷暴現象の観測とミューオン電荷非対称性の役割

論文タイトル: The muon charge asymmetry and the directional distribution of thunderstorm events observed by the GRAPES-3 muon telescope
投稿誌: JCAP (Journal of Cosmology and Astroparticle Physics) への提出準備中
著者: B. Hariharan ら (GRAPES-3 共同研究チーム)

1. 研究の背景と課題 (Problem)

GRAPES-3 実験（インド、ウーティ）のミューオン望遠鏡（G3MT）は、2011 年 4 月から 2020 年 12 月までの 10 年間にわたり、雷暴（スコール）に伴う二次宇宙線ミューオン強度の変動（$\Delta I_\mu$）を継続的に観測してきました。雷雲内のギガボルト（GV）規模の電場は、通過するミューオンを加速または減速させ、地上でのミューオン強度に変化をもたらします。

しかし、観測データには極めて顕著な方向性非対称性が存在していました。具体的には、東方向からの雷暴イベントの検出数が、西方向からのそれよりも約 6 倍多いという結果でした。雷暴自体が特定の方向に偏って発生する物理的な理由はないため、この偏りの原因は、雷雲の電場そのものではなく、ミューオンの検出プロセスや大気中での伝播に起因する何らかの要因であると考えられました。本論文の目的は、この方向性非対称性の物理的メカニズムを解明し、ミューオン電荷比（$R_\mu$）の役割を定量的に検証することです。

2. 研究方法 (Methodology)

本研究では、以下の手法を組み合わせて分析を行いました。

• 観測データの解析:

  • G3MT によって記録された 487 件の雷暴イベント（2011 年 4 月〜2020 年 12 月）を分析。
  • ミューオン強度変動の閾値 $|\Delta I_\mu| > 0.3\%$ を用いてイベントを選別。
  • 観測領域を 9 方向（北西、北、北東、西、垂直、東、南西、南、南東）に分割し、各方向のイベント分布を評価。

• モンテカルロシミュレーション (CORSIKA):

  • 大気中の宇宙線シャワー発展をシミュレートするコード「CORSIKA」を使用。
  • 高エネルギー相互作用生成器（SIBYLL, EPOS-LHC, QGSJETIII, DPMJETIII）と低エネルギー生成器（FLUKA, UrQMD）の 7 種類の組み合わせを用いて、ミューオン生成を詳細にモデル化。
  • 雷雲電位（-3 GV 〜 +3 GV）を仮定し、ミューオン強度変化と電荷比 $R_\mu$（正ミューオン数 $N_{\mu+}$ / 負ミューオン数 $N_{\mu-}$）の関係を計算。
  • 地磁気の影響を考慮し、西から東へ向かう方向ごとのカットオフ剛性（15 GV 〜 24 GV）の違いをシミュレーションに反映。

• 対照実験:

  • $R_\mu$ を一定値（西方向の値）に固定した場合のシミュレーションを行い、方向性非対称性が消滅するかどうかを確認。

3. 主要な成果と結果 (Key Contributions & Results)

3.1 観測された強い東西非対称性

• 観測された 487 件のイベントのうち、東方向（東、北東、南東）の合計は全体の**82%を占め、西方向（西、北西、南西）は14%**に過ぎませんでした。
• 東西方向のイベント数の比率は 5.9 ± 0.8 であり、東方向が圧倒的に多いことが確認されました。

3.2 ミューオン電荷比 ($R_\mu$) の方向依存性

• CORSIKA シミュレーションにより、G3MT の視野内における $R_\mu$ が方向によって系統的に変化することが明らかになりました。
  • 西方向 (WG): $R_\mu \approx 1.14$
  • 東方向 (EG): $R_\mu \approx 1.37$
• この差は、地球磁場による一次宇宙線（主に陽子）の偏向効果（東西効果）に起因します。東から来る陽子は西から来るものよりも強い磁気的遮断（高いカットオフ剛性）を受けるため、高エネルギーの陽子成分が相対的に増加し、結果として正ミューオンの生成割合が高まります。

3.3 非対称性のメカニズム解明

• 雷雲は主に正の電位（正電荷）を持つことが多く、これは正ミューオン（$\mu^+$）を減速させ、負ミューオン（$\mu^-$）を加速させます。
• $R_\mu > 1$（正ミューオンの方が多い）であるため、正ミューオンの減速による強度減少が、負ミューオンの加速による増加を上回り、結果として地上でのミューオン強度は**減少（欠乏）**します。
• 東方向は $R_\mu$ が大きいため、同じ雷雲電位に対してミューオン強度の減少幅（$\Delta I_\mu$）が西方向よりも大きくなります。
• 検出閾値（$|\Delta I_\mu| > 0.3\%$）を越えるイベント数が、東方向では西方向よりもはるかに多くなるため、検出バイアスとして東方向のイベントが過剰に観測されることになります。

3.4 シミュレーションによる再現

• 実際の $R_\mu$ の方向依存性を用いたシミュレーションでは、東西イベント数の比が 5.09 ± 0.02 となり、観測値（5.9 ± 0.8）とよく一致しました。
• 逆に、$R_\mu$ を一定値に固定したシミュレーションでは、東西のイベント数比は 1.00 となり、非対称性が消失しました。
• これにより、観測された非対称性が機器の偏りではなく、地磁気によるミューオン電荷非対称性の物理的効果であることが証明されました。

3.5 ハドロン相互作用モデルの影響評価

• 7 種類のハドロン相互作用生成器の組み合わせを用いた比較により、$R_\mu$ の絶対値はモデル間で最大 7% 程度のばらつきがあることが示されました。
• これに伴い、導出される雷雲電位の推定値には最大 14% 程度の不確実性が生じますが、「東方向で $R_\mu$ が大きく、イベント数が多い」という傾向はすべてのモデルで共通しており、結論の頑健性が確認されました。

4. 意義と結論 (Significance & Conclusions)

本研究は、以下の点で重要な科学的貢献を果たしました。

1. 雷暴観測の方向性偏りの解明: 雷暴イベントの東西非対称性が、雷雲の物理的な偏りではなく、地磁気によるミューオン電荷比（$R_\mu$）の方向依存性によって引き起こされる「検出バイアス」であることを初めて定量的に証明しました。
2. ミューオン観測の信頼性向上: 雷雲電位の推定において、$R_\mu$ の変化を適切に考慮することが不可欠であることを示しました。特に、地磁気の影響を無視すると、電位の推定に系統的な誤差が生じる可能性があります。
3. 大気電気構造のプローブとしての G3MT: ミューオン観測が、大気中のギガボルト規模の電場構造を研究するための強力なプローブであることを再確認しました。地磁気と局所的な気象条件の相互作用を理解することで、より精密な大気電気学の研究が可能になります。

結論として、GRAPES-3 による雷暴観測データの方向性偏りは、地磁気場が一次宇宙線および二次ミューオンに及ぼす選択的フィルタリング効果（東西効果）に起因するミューオン電荷非対称性によって説明可能です。この発見は、宇宙線ミューオンを用いた大気電気現象の研究における基礎物理の理解を深め、将来の観測データの解釈精度を高める上で極めて重要です。