Search for a solar-bound axion halo using the Global Network of Optical Magnetometers for Exotic physics searches

技術要約：GNOMEを用いた太陽束縛アクシオンハローの探索

問題と動機
ダークマターの性質は、天体物理学および粒子物理学における中心的な未解決問題である。超軽量ボソン、特に質量 $m_a \ll 1$ eV のアクシオンまたはアクシオン様粒子（ALP）は、十分に動機付けられた候補である。近年の理論的研究は、太陽のような高密度の天体が、重力的にアクシオンを捕捉できる可能性を示唆している。もしアクシオンが四次自己相互作用を持つならば、これらの相互作用は、アクシオンが太陽の重力ポテンシャルの基底状態へと熱化することを可能にするために必要なエネルギー散逸メカニズムを提供し得る。このプロセスにより、「太陽アクシオンハロー」が形成される。これは、太陽を中心とする基底状態の構成（ground-state configuration）である。

標準的な銀河ダークマターハロー（熱平衡状態にあり、確率論的である）とは異なり、基底状態（$n=1, \ell=0$）にある太陽アクシオンハローは、明確な密度プロファイル（$\rho \propto e^{-2r/R_\star}$）と、アクシオンのコンプトン周波数 $\omega_a = m_a c^2 / \hbar$ におけるコヒーレントで単色的な振動を示す。このようなハローの存在は、平均的な銀河ダークマター密度と比較して、地球近傍におけるアクシオンの過剰密度を示唆しており、検出感度を高める可能性がある。このハロー内のアクシオン場は、勾配相互作用（線形または二次）を介して原子スピンに結合し、地球の自転および太陽に対する公転運動によって生じる特定の日常的な変調パターンを持つ擬似磁場（$B_p$）を生成する。

手法
本研究では、GNOME（Global Network of Optical Magnetometers for Exotic physics searches）の第5回サイエンスラン（SR5）のデータを使用している。SR5は、地理的に分散した12箇所の磁気シールドされた原子磁力計による69日間のデータで構成されている。

1. 信号モデル: 解析は、アクシオンが太陽の重力ポテンシャルの基底状態にあると仮定している。その結果として生じるアクシオン場 $a(r,t)$ は、陽子スピンとの線形（$H_l \propto S \cdot \nabla a$）および二次（$H_q \propto S \cdot \nabla a^2$）の勾配結合を通じて、擬似磁場 $B_p$ を生成する。信号は、2つの成分（太陽への空間勾配による径方向成分と、地球の公転運動による「アクシオン風」による横方向成分）からなる。これらの成分の相対的な強さは、アクシオンの質量（周波数）に依存する。
2. データ前処理: 各ステーションの生の時系列データは、商用電源ノイズ（50/60 Hz）および緩やかなドリフト（0.0003 Hzの高域通過フィルタ）を除去するためにフィルタリングされた。データは64 Hzのサンプリングレートにビン分割された。ゲインのドリフトを補正するために、1時間ごとの校正パルスが使用された。
3. タイムシフト・クロス相関: グローバルな信号と局所的なノイズを区別するために、本解析ではクロス相関パイプラインを採用している。信号の日常的な変調振幅は、各ステーションの地理的位置と方位に基づいて変化するため、異なるステーション間のデータは、予測される信号のピークを一致させるように時間シフト（$\Delta t_j$）される。
4. スペクトル解析: ペアごとのクロス相関は20分間のウィンドウ内で計算される。得られた信号はパワースペクトル密度（PSD）に変換され、正規化された後、重み付けされたネットワークレベルのスペクトルへと統合される。
5. 候補選択とフィッティング:
   • ステージ1: 実行中のすべての日で正規化PSDが閾値（$P_\star = 3.35$）を超える周波数を候補として特定する。
   • ステージ2: 候補となる周波数に対して、組み合わせられたPSDの日常的な変調パターンを、アクシオンハローのシミュレーションから導出された理論モデルに適合させる。テスト統計量 $Z$ を用いて、適合の質をアクシオンモデルと定数（ノイズのみ）モデルの間で比較する。
   • 有意性: 候補が有意であると見なされるのは、適合された振幅がゼロではなく、かつテスト統計量 $Z$ の $p$ 値が、$Z$ 値の経験的分布から導出されたグローバルな有意性閾値（$p_{crit}$）を下回る場合のみである。

主な貢献

• 初の専用探索: 本研究は、世界的な磁力計ネットワークを用いた、重力的に束縛された太陽アクシオンハローの初の専用探索を提示している。
• 信号モデリング: 著者らは、基底状態の太陽ハローにおける、線形および二次のアクシオン・陽子勾配結合によって生成される擬似磁場の包括的な信号モデルを開発した。これには、周波数スペクトル全体にわたる、径方向成分と横方向成分の相互作用および日常的な変調を明示的に考慮したモデルが含まれる。
• 解析パイプライン: 本研究は、世界規模のネットワークにおけるグローバルな信号を一致させ、相関のない局所ノイズを効果的に抑制しつつ、太陽束縛場のコヒーレンスを保持する、タイムシフトを利用した堅牢なクロス相関パイプラインを実装している。

結果

• 検出なし: 12箇所のGNOMEステーションによる69日間のデータの解析の結果、太陽アクシオンハローと一致する統計的に有意な候補信号は認められなかった。テスト統計量 $Z$ の観測された $p$ 値は、グローバルな有意性閾値を下回るものはなかった。
• 上限値: 検出が得られなかったため、本研究では、約 0.05 – 20 Hz の周波数範囲における、アクシオン誘起擬似磁場の振幅に対する95%信頼区間の上限値を設定した。
• 結合制約: これらの磁場の上限値は、線形（$1/f_l$）および二次（$1/f_q$）のアクシオン・陽子結合定数への制約へと変換された。制約は、以下の2つのベンチマークとなるハロー密度シナリオに対して導出された：
  1. 行星エフェメリスによって許容される最大過剰密度。
  2. 四次自己相互作用による重力捕捉モデルが予測する密度。
• 二次結合の感度: 二次結合の場合、導出された上限値は、広範なパラメータ空間において、既存の天体物理学的境界（例：超新星1987A）を2桁以上上回っている。線形結合については、想定されるハロー密度に応じて、既存の境界と同等または補完的な結果となっている。

意義と主張
本論文は、重力捕捉および熱化によって形成された太陽アクシオンハローの存在に対する、最初の実験的な制約を確立したと主張している。結果は、GNOMEが、太陽束縛場特有の日常的な変調シグネチャを標的とした、グローバルなネットワークアプローチを通じて、エキゾチックなスピン結合を探索できる能力を持っていることを示している。

著者らは、これらの制約が、アクシオンが基底状態にあるという仮定、および（四次自己相互作用による）形成メカニズムが効率的であるという仮定に基づいていることを控えめに述べている。また、二次結合については、その結果が従来の天体物理学的限界を大幅に改善するものであると述べている。本研究は、得られたデータが、トポロジカル欠陥や確率論的な場など、他のエキゾチックな物理学の対象に対しても再解析が可能であることを強調している。また、アルカリ金属・希ガス共磁力計を用いた将来のGNOMEの反復は、陽子および中性子のスピンの両方に対してさらなる感度向上をもたらすであろう。