Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

技術的サマリー：太陽圏磁場が宇宙線ミューオンの影に及ぼす影響のモデル化

問題定義
宇宙線（CR）は、主に陽子とアルファ粒子から構成され、太陽系を通過する際に電磁場によって偏向される。地球の地磁気（GMF）は詳細にマッピングされているが、太陽圏磁場（HMF）は複雑でモデル化が困難である。太陽と月は入射する宇宙線を吸収し、地球から観測可能な宇宙線フラックスにおける「影」を形成する。しかし、これらの影は静的ではなく、HMF と GMF によってシフトし歪められる。MINOS+ コラボレーションは、MINOS 遠方検出器からのミューオンデータを用いて、太陽活動極小期、太陽活動極大期、および 13 年間の全期間にわたるこれらの太陽影の位置を特定した。中心的な課題は、観測された太陽ミューオン影のシフトと歪みが、HMF の標準的なパーカー・スパイラルモデルと GMF の双極子モデルの組み合わせと整合するか、そして整合しない場合、理論と観測を調和させるためにモデルやパラメータにどのような修正が必要かという点である。

手法
本研究は、Soudan 鉱山地下 2070 メートル水相当の深さに位置する MINOS 遠方検出器からの 13 年間のミューオンデータ（2003 年 7 月～2016 年 7 月）を利用した。解析は、入射天頂角が 80 度未満で、軌跡を再構成するのに十分な軌道長を持つ単一ミューオン事象に焦点を当てた。

1. 影の局所化: チームは、影の位置を特定するために $\Lambda$ 確率検定（対数尤度比）を採用した。これには、太陽近傍の観測粒子数（$N_{obs}$）を、理論的な平坦な背景（$N_{back}$）および理論的な影テンプレート（$N_{th}$）と比較することが含まれる。テンプレートは、太陽の角直径（0.52 度）および岩石被覆層における散乱効果を考慮した。
2. 統計的検証: 観測された影の位置の有意性を決定するために、モンテカルロシミュレーションを用いてランダムな背景ヒストグラムを生成した。影を模倣するランダムな変動の確率を計算し、高い信頼水準（全データセットで最大 6$\sigma$）を得た。
3. 計算モデル化: 影のシフトを解釈するために、著者は検出器から太陽へ向かう宇宙線の時間逆転軌跡をシミュレーションした。
   • 磁場モデル: シミュレーションでは、HMF に対するパーカー・スパイラルモデルと GMF に対する双極子モデルを重ね合わせた。パーカー・スパイラルのパラメータには、一定の太陽風速度（$u_{sowi} = 400$ km/s）および太陽自転周波数が含まれた。GMF は双極子モーメントを用いてモデル化され、太陽中心の座標系に変換された。
   • 粒子力学: ローレンツ力方程式は、4 次ルンゲ・クッタ法を用いて数値積分された。粒子は相対論的として扱われ、エネルギーはべき乗則分布（$E^{-2.75}$）に従う。
   • シナリオ: シミュレーションでは、様々な磁場極性比（50/50、33/67、67/33）と、2 種類の影テンプレート（単純な円形影と、実際の $\Lambda$ 確率ヒストグラムから導出された「歪んだ」影）をテストした。
   • パラメータ変異: 一次宇宙線のエネルギーカットオフ（$E_0$）は自由パラメータとして扱われ、異なるエネルギー分布がモデルと観測された影の位置を調和させるかどうかを決定した。

主要な貢献

• 影の特性評価: 本論文は、3 つの異なる時間枠（13 年間の全データセット、太陽活動極小期（2007–2009）、太陽活動極大期（2012–2014））において、高い統計的有意性で太陽ミューオン影を特定することに成功した。
• テンプレート開発: 著者は、観測された $\Lambda$ 確率分布に基づいて「歪んだ」影サンプリングテンプレートを作成する方法を導入し、完全な円形影という仮定を超えた。
• モデル検証: 本研究は、パーカー・スパイラルモデルを観測データに対して厳密にテストし、シミュレーションされた影のシフトと観測された影のシフトの一致度を、赤経（$\Delta RA$）および赤緯（$\Delta DEC$）において定量化した。
• エネルギー分布解析: 本論文は、影の位置が一次宇宙線のエネルギー分布に対してどの程度敏感か、特に低エネルギーカットオフ（$E_0$）の影響をテストして調査した。

結果

• 影のシフト: 観測された影は、太陽活動サイクルの段階に応じて明確なシフトを示した。太陽活動極小期の影は、太陽活動極大期の影とは異なるシフトを示し、全データセットはこれらのシフトの組み合わせを表した。
• モデルの整合性:
  • パーカー・スパイラルモデルは、太陽活動極小期のデータと最もよく一致した。
  • 同モデルは、特に $\Delta DEC$ シフトに関して、太陽活動極大期のデータと最も一致しなかった。
  • 歪んだ影テンプレートは、太陽活動極大期および全データ期間においてより良好な性能を示したが、単純な円形テンプレートは太陽活動極小期のデータにおいてより効果的であった。
• 極性効果: シミュレーションは、磁場極性比が影の $\Delta DEC$ 位置に直接影響することを確認した。データは、特に太陽活動極大期において、正極性が支配的な HMF を支持するように見えた。
• エネルギーカットオフの不一致: パーカー・スパイラルモデルを用いてシミュレーションされた影の平均を理論的中心（0,0）に整合させるためには、一次宇宙線のエネルギーカットオフ（$E_0$）を22.8 $\pm$ 1.2 TeVに調整する必要があった。これは、実験的に決定された約 10 TeV の値よりも著しく高い（2～3 倍）ものである。
• 天頂角補正: エネルギーカットオフの天頂角依存性に対する補正が適用された。これによりスペクトルは観測されたミューオンエネルギー分布に一致するよう調整されたが、観測された $\Delta RA$ シフトを完全に説明するには不十分であった。

意義と結論
本論文は、宇宙線太陽ミューオン影が観測可能であり、HMF の構造を探るために利用できると結論づけている。結果は、太陽活動極小期の影が、標準的なパーカー・スパイラルモデルと双極子 GMF の組み合わせと最も整合することを示している。しかし、非現実的に高いエネルギーカットオフを仮定しない限り、このモデルは太陽活動極大期の影および全データセットを完全に再現できない。

著者は、この不一致は、太陽活動極大期の条件下ではパーカー・スパイラルモデルが過度の単純化であることを示唆していると主張している。彼らは、モデル化されていない要因、すなわち中性電流シート傾斜（「バレリーナモデル」）、衝撃波前面、および複雑な極性変動が、観測されたシフトに寄与している可能性があると提案している。本研究は、偏向の一般的な方向がパーカー・スパイラルと一致している一方で、太陽活動極大期における HMF の詳細な構造は、単純なスパイラルモデルが提供するものよりも複雑なモデル化を必要とすることを確立している。この作業は、HMF モデルを制約する手段としてのミューオン影観測の有効性を示しているが、太陽活動サイクルの変動を説明するためには、より洗練されたシミュレーションが必要であることを浮き彫りにしている。