Natural Saturation Of The Sterile Neutrino Dark Matter Resonant Production By a High Lepton Flavor Asymmetry In Primordial Plasma

技術要約：ステライル・ニュートリノ・ダークマターの共鳴生成における自然な飽和

問題提起
keV質量領域のステライル・ニュートリノは、初期宇宙におけるアクティブ・ステライル振動を通じて生成される、有力なダークマター候補である。非共鳴的な熱的生成（ドドルソン＝ウィドロウ機構）はX線観測によって大部分が排除されているが、原始レプトン非対称性によって駆動される共鳴生成（シ・フラー機構）は依然として有効な経路である。標準的な取り扱いでは、初期レプトン非対称性（$L$）を増大させると、物質ポテンシャルが増強され、共鳴変換が促進されて最終的なステライル・ニュートリノの存在量が増加すると示唆されている。

しかし、この線形関係は、以下の2つの無視されていた物理効果によって崩壊する：

1. 共鳴の遅延： 初期非対称性が大きくなると、共鳴生成の時期がより後期の宇宙論的時刻（より低い温度）へとシフトする。
2. アクティブ・ニュートリノ振動による枯渇： このより低い温度（$T \lesssim 20$ MeV）において、アクティブ・ニュートリノの振動が効率的になり、3つのアクティブ・フレーバー間のフレーバー非対称性を再分配する。これにより、共鳴を駆動するために必要な特定のフレーバー非対称性が、生成期の前または最中に大幅に減少、あるいはほぼ消失する可能性がある。

さらに、非常に大きな非対称性においては、電気電荷化学ポテンシャルが荷電パイ中間子の質量に近づき、パイ中間子凝縮が発生して、電荷中性条件および熱力学的背景を変化させる可能性がある。先行研究（例：文献[18]）は、$L < 0.1$ の非対称性に限定されており、非対称性の動的な枯渇やパイ中間子凝縮を十分に考慮していなかった。

手法
著者らは、宇宙論的背景とステライル・ニュートリノ生成を自己整合的に進化させるための、結合された数値フレームワークを開発した：

• 熱力学的背景： システムは、電荷中性とバリオン非対称性の制約下で、化学ポテンシャル（$\mu_{\nu_e}, \mu_{\nu_\mu}, \mu_{\nu_\tau}, \mu_Q$）を解く。ハドロン部門は核子と軽メソンを用いてモデル化されている。
• パイ中間子凝縮： 非対称性が大きい場合（$|\mu_Q| \approx m_\pi$）、著者らはクォーク・メソンモデル（文献[19]に従う）を用いてパイ中間子凝縮相を組み込んでいる。フェルミオンの熱力学的ポテンシャルの寄与を無視した簡略化されたメソン処方を採用しており、これにより凝縮振幅の解析解を得ている。これは、物質ポテンシャルに寄与する電荷密度および化学ポテンシャルを修正する。
• アクティブ・ニュートリノ振動： 著者らは、アクティブ・ニュートリノの密度行列に関する量子キネティック方程式を解くために、COFLASYコードの修正版を利用している。このコードは、Pythonによるバックグラウンド・ソルバーと結合されている。
  • 結合スキーム： 進化は減少する温度グリッド上で進行する。各ステップ $T_i$ において、フレーバー非対称性がバックグラウンド・ソルバーに渡され、化学ポテンシャルおよび熱力学的量が更新される。これらはその後、COFLASYに渡され、アクティブ・ニュートリノ部門を $T_{i+1}$ まで進化させる。結果として生じるフレーバー非対称性の変化（$\delta L_\alpha$）が、次のステップの状態を更新する。
  • 初期化： 2つの領域がテストされている：「継続（continuation）」（完全な動的状態ベクトルを保持する）および「フレッシュステップ（fresh-step）」（更新された非対称性に基づいて密度行列を再初期化する）。フレッシュステップ法は、大きな非対称性において（継続法で数値的不安定性が生じるため）より安定していることが示された。
• ステライリ・ニュートリノ生成： ステライル・ニュートリノの存在量は、大きな非対称性の領域における量子キネティック方程式を再現する、一般化された有効確率処方（文献[18]）を用いて計算される。生成は、共鳴の狭幅近似を用いて解析的に扱われ、共鳴通過に関するボルツマン方程式を積分する。

主な貢献

1. 飽和メカニズムの特定： 本論文は、原始レプモン非対称性を増大させても、ステライル・ニュートリノの存在量が無限に増加するわけではないことを示している。ある閾値を超えると、共鳴が、アクティブ振動が駆動となる非対称性を枯渇させる温度へとシフトし、その結果、生成率が低下または飽和する。
2. パイ中間子凝縮の導入： 著者らは、大きな非対称性におけるパイ中間子凝縮の効果を組み込み、それが電荷中性条件および物質ポテンシャルをどのように修正するかを示している。これは、$L \gtrsim 0.4$ における正確な共鳴追跡において極めて重要である。
3. 結合進化フレームワーク： 著者らは、宇宙論的背景（パイ中間子凝縮を含む）とアクティブ・ニュートリノのフレーバー振動を同時に進化させる、自己整合的な数値実装を提示している。これは、高い非対称性における共鳴生成をモデリングするために必要不可欠なステップである。
4. 混合角への制約： 最大限に達成可能な存在量を初期非対称性に対してマッピングすることにより、観測されたダークマター密度を説明するために必要な最小のアクティブ・ステライル混合角（$\sin^2 2\theta$）の更新された制約を導出している。

結果

• 存在量の飽和： 代表的なフレーバー構成（逆階層および順階層シナリオ）について、ステライル・ニュートリノの存在量 $\Omega_{\nu_s}$ は、初期非対称性 $L$ とともに最初は増加する。しかし、$L$ が増加する（通常 $L \gtrsim 0.5-1.0$）と、共鳴がより低い温度へとシフトし、そこでアクティブ振動が関連するフレーバー非対称性を急速に減少させる。これにより、生成効率が急激に低下する。さらに大きな非対称性では、共鳴は、すべての非対称性がゼロへと駆動される電子フレーバー振動の時期に入り、存在量のさらなる成長を妨げるプラトー（停滞状態）に入る可能性がある。
• 最小混合角： この飽和効果は、要求される混合角の下限を意味する。非対称性が最大であっても、混合角は任意に小さくすることはできない。著者らは、共鳴メカニズムおよび現在の観測限界と一致する、$\sin^2 2\theta$ のおよそ2桁（例：$m_s \sim 40$ keV に対して $10^{-18}$ から $10^{-16}$）にわたるパラメータ空間を特定している。
• 運動量スペクトル： 得られる運動量分布は非熱的である。非対称性が増加するにつれ、平均運動量は増加し（より高温のダークマター）、$\langle x \rangle \approx 3.2$ に近づく。生成がアクティブ振動と重なる場合、複雑なスペクトル特徴（例：二峰性ピーク）が現れる。しかし、著者らは、これらの重なり合う領域における数値的な不確実性が、詳細なスペクトル形状の完全なロバストな決定を妨げていると述べている（特にケースI：タウ・フレーバー混合の場合）。
• 数値的安定性： 「フレッシュステップ」初期化法は、「継続」法よりも堅牢であり、アクティブ振動が急速に進む大きな非対称性において、後者に観察される数値的不安定性を回避できることが証明された。

意義と主張
著者らは、本研究が、共鳴的なステライル・ニュートリノ生成に対する「自然な飽和メカニズム」を提供していると主張しており、これは大きなレプトン非対称性に依存するモデルに対する追加の制約として機能する。遅延共鳴とアクティブ振動による枯訳の相互作用が、このメカニズムを通じて達成可能な最大ダークマター密度を制限する。

本論文の結論は以下の通りである：

• 現在のX線望遠鏡の感度は、このメカニズムを完全に探索するには不十分であり、予測されるスペクトルのピークを検出するためには、混合角に対する一桁上の感度向上が必要である。
• 最大の存在量および運動量スペクトルの正確な決定は、現在、結合進化の数値的な複雑さによって制限されている。背景進化と量子キネティック方程式を統合した統一的な数値フレームワークが、アクティブ振動とステライル生成の重なり合う時期を解決するために必要である。
• 大きな原始非対称性は、パイ中間子凝縮を介した一次相転移を引き起こし、確率的な重力波背景を生成する可能性がある。そのような信号の検出は、共鳴的なステライル・ニュートリノ・ダークマターに必要な大きな非対称性の存在を示す、独立した証拠となり得る。

著者らは、最も複雑な領域（ケースI）における詳細な運動量スペクトルに関して慎重な姿勢を維持しており、現在の数値結果に見られる小規模な不規則性は、さらなる収束チェックなしに物理的な解釈を与えるべきではないと認めている。