Exploring non-equilibrium effects in sequential freeze-in

この論文は、多成分ダークセクターの凍結生成（freeze-in）において局所熱平衡からの逸脱が重要であることを示し、最小の二スカラーモデルを用いた研究で、従来の数密度アプローチと位相空間レベルの計算との間にダークマターの存在量で最大 1 桁の差異が生じることを明らかにしています。

要約

この論文は、宇宙の謎「ダークマター（暗黒物質）」がどのようにして生まれたか、そしてその「生まれ方」が私たちが観測する未来にどう影響するかを、新しい視点から探求した研究です。

専門用語を避け、わかりやすい比喩を使って解説しましょう。

1. 物語の舞台：「見えない国」と「熱いお風呂」

まず、宇宙の初期を想像してください。そこは超高温の「熱いお風呂（標準模型の粒子たち）」に満ちています。
一方、ダークマターは、このお風呂とはほとんど触れ合わない「見えない国」に住む住人です。

• 従来の考え（お湯に溶ける砂糖）：
  昔の研究者は、「見えない国の住人」も、お風呂の温度に合わせて温まり、均一に分布している（熱平衡状態）と仮定していました。つまり、お風呂の温度と同じ温度で、均一に溶けている砂糖のようなイメージです。

• この論文の発見（冷たい氷のかけら）：
  しかし、この論文の著者たちは、「待てよ！見えない国の住人は、お風呂の温度に追いつく前に、独自の動きをしているかもしれない」と考えました。彼らは「非平衡（非熱的）」な状態、つまりお風呂の温度とは異なる独自の温度や、偏った分布を持っている可能性があるのです。

2. 実験室：「2 人の双子」と「仲介者」

この研究では、ダークマターを生成する仕組みを、シンプルなモデルでシミュレーションしました。

• 登場人物：

  1. H（ヒッグス粒子）： お風呂の主人。
  2. ϕ（ファイ）： 仲介者（メッセンジャー）。H から生まれ、ダークマターへつなぐ役割。
  3. S（エス）： 真のダークマター（安定して残るもの）。

• 生成のプロセス（「連続した注文」）：

  1. H が ϕ を生み出す。
  2. ϕ が集まって、S（ダークマター）に変化する。
     この「H → ϕ → S」という連続したプロセスが、この研究の核心です。

3. 何が起きたのか？「お湯の温度」を信じるな！

ここで面白いことが起きました。研究者たちは、3 つの異なる方法でダークマターの量を計算しました。

1. 方法 A（従来の計算）： 「ϕ も S も、お風呂（H）と同じ温度で均一に動いている」と仮定する。
2. 方法 B（少し詳しい計算）： 「ϕ と S は、お風呂とは違う独自の温度を持っているかもしれない」と仮定する。
3. 方法 C（完全な計算）： 「個々の粒子がどう動いているか、すべて追跡する（非平衡な状態をそのまま計算する）。」

結果：

• 方法 A（従来の考え）と方法 C（完全な計算）の間には、最大で「10 倍」もの差が出ました！
• なぜ？
  • 比喩： お風呂から「熱い石（高エネルギーの粒子）」が飛び出してくるのを想像してください。
  • 従来の計算（方法 A）は、「お風呂全体が均一に熱いから、飛び出す石も平均的な熱さだ」と考えます。
  • しかし、実際の計算（方法 C）では、「飛び出す石は、実は非常に熱い石だけが選ばれて飛び出している」ことがわかりました。
  • この「熱い石」だけが、次の段階でダークマター（S）を作る鍵を握っているのです。
  • 従来の計算はこの「熱い石」の存在を見逃し、ダークマターの量を過大評価してしまう（または過小評価してしまう）ことがわかりました。

4. なぜこれが重要なのか？「未来へのメッセージ」

この発見は、ダークマターを探す未来のプロジェクトにとって非常に重要です。

• 探偵の道具：
  私たちは、ダークマターが星と星の間で衝突して光（ガンマ線）を放つ「間接検出」や、巨大な加速器で長生きする粒子を探す「前方物理実験」でダークマターを見つけようとしています。
• 誤った地図：
  もし、従来の「均一な温度」という間違った仮説に基づいて予測を立てると、「ここに来れば見つかる！」と地図を間違えてしまいます。
• 正しい地図：
  この論文は、「実は、粒子の動きはもっと複雑で、偏っている。だから、予測される信号の強さや場所も、10 倍も変わる可能性がある」と警告しています。

まとめ

この論文は、**「ダークマターの誕生は、お湯に溶ける砂糖のように均一ではない。むしろ、熱い石だけが飛び跳ねるような、激しく偏ったプロセスだったかもしれない」**と教えてくれます。

もしこの「非平衡な動き」を無視して計算を続ければ、私たちがダークマターを見つけるための実験計画が、大きく外れてしまう可能性があります。そのため、より精密な計算ツールを開発し、粒子の「個々の動き」まで追跡することが、今後の宇宙の謎を解く鍵だと主張しています。

一言で言えば：
「ダークマターの量や性質を正しく知るには、単なる『平均値』ではなく、粒子たちがどう『踊っているか』まで詳しく見る必要があるよ！」という、宇宙論における重要な気づきを提供した研究です。

技術概要

論文の技術的サマリー：「Exploring non-equilibrium effects in sequential freeze-in」

1. 研究の背景と問題提起

暗黒物質（DM）の正体とその生成メカニズムは、現代の素粒子物理学および宇宙論における未解決の重大な課題です。特に、DM の運動量分布 $f(p)$ が熱平衡分布 $f_{eq}(p)$ を追従しているという仮定は、従来の WIMP（Weakly Interacting Massive Particles）の凍結アウト（freeze-out）モデルや、標準模型（SM）プラズマとの相互作用が極めて弱い凍結イン（freeze-in）モデルにおいて一般的に用いられています。

しかし、多成分からなるダークセクター（暗黒セクター）において、内部ダイナミクスが複雑な場合、この「熱平衡分布」という近似は破綻する可能性があります。特に、中間粒子を介した**逐次凍結イン（sequential freeze-in）**シナリオでは、最終的な DM の残留密度が、中間粒子の非熱的な運動量分布に敏感に依存することが予想されます。これまでの研究では、非平衡効果の定量的な影響、特に位相空間（phase-space）レベルでの詳細な評価が体系的に行われておらず、従来の数密度（number-density）アプローチとの乖離がどの程度生じるかが不明確でした。

2. 提案されたモデルと手法

著者らは、この問題を解明するために、ヒッグスポータルを介して SM と結合する2 成分のスカラーダークセクターモデルを提案・解析しました。

モデルの概要

• 構成粒子: 2 つの実スカラー粒子 $\phi$（媒介粒子）と $S$（安定した DM 候補）。
• 対称性: $S$ は安定であり DM となりますが、$\phi$ は SM 粒子へ崩壊する不安定な粒子です。
• 生成メカニズム:
  1. 直接凍結イン: ヒッグス崩壊 $h \to SS$ による直接生成。
  2. 逐次凍結イン: ヒッグス崩壊 $h \to \phi\phi$ により生成された $\phi$ が、$SS$へ変換される過程（$\phi\phi \to SS$）を介した生成。
  3. ダーク凍結アウト: 変換率がハッブル率を超え、ダークセクター内で熱平衡が達成された後の凍結アウト。

計算手法の比較

DM の残留密度を計算する際、以下の 3 つの異なるボルツマン方程式アプローチを比較・適用しました。

1. 数密度ボルツマン方程式（nBE）:

   • 従来のアプローチ。すべてのダークセクター粒子が SM プラズマと運動学的平衡（Kinetic Equilibrium）にあると仮定し、数密度のみを追跡します。
   • 分布関数は SM 温度 $T_{SM}$ における平衡分布と仮定されます。

2. 結合ボルツマン方程式（cBE）:

   • 数密度だけでなく、有効温度（分布の 2 次モーメント）も追跡します。
   • 各ダークセクター粒子が自己平衡状態にある（分布関数の形状は平衡分布だが、温度 $T_i$ は $T_{SM}$ と異なる）と仮定します。

3. 完全ボルツマン方程式（fBE）:

   • 本論文の核心。分布関数 $f(p)$ の形状そのものに対する仮定を置かず、運動量空間全体での時間発展を数値的に解きます。
   • 非熱的な分布の歪み（特に高運動量テール）を直接考慮します。
   • 実装には、2 成分モデルに対応したコード「DRAKE」の改良版を使用しました。

3. 主要な結果

3.1 ベンチマーク点における非平衡効果の定量化

6 つのベンチマーク点（BM1-BM6）を選定し、異なる生成レジーム（直接、逐次、ダーク凍結アウトなど）における結果を比較しました。

• 直接凍結イン（BM1）:
  • DM が直接生成される場合、nBE、cBE、fBE の結果は一致します。熱平衡分布の仮定が有効であることを確認しました。
• 逐次凍結イン（BM2, BM3）:
  • 大きな乖離: 中間粒子 $\phi$ を介して DM が生成される場合、nBE と fBE の間で最大で 1 桁（オーダー）の残留密度の差が生じました。
  • メカニズム:
    1. 運動量分布の歪み: 変換過程 $\phi\phi \to SS$ は運動量保存則により、$\phi$ の高運動量テール（high-momentum tail）に強く依存します。nBE は SM 温度での平衡分布を仮定するため、この高運動量成分を過大評価します。一方、fBE は非平衡的な分布の歪み（高運動量テールの枯渇）を正しく捉え、変換率を抑制します。
    2. 温度の分離: cBE および fBE では、$\phi$ の温度 $T_\phi$ が $T_{SM}$ よりも低くなることを示し、これがさらに変換率を抑制しました。
  • BM3（より重い DM）では、運動学的抑制がさらに強まり、nBE と fBE の差が顕著になりました。
• ダーク凍結アウト（BM5, BM6）:
  • 変換結合定数 $\lambda_{S\phi}$ が大きい場合、ダークセクター内部で熱平衡が迅速に達成されます。この場合、cBE は fBE とよく一致しますが、nBE とは依然として差が生じます（ダークセクターの温度が SM 温度と異なるため）。

3.2 相互作用強度への依存性

結合定数 $\lambda_{S\phi}$（変換）および $\lambda_{h\phi}$（ポータル）を変化させたスキャンを行いました。

• 変換率がハッブル率と同程度になると、逐次凍結インからダーク凍結アウトへ遷移します。
• 非熱効果（fBE と nBE の乖離）は、変換過程が運動量に敏感な領域（サブスレッショルド領域）で最大となり、結合定数の変化に対して非単調な振る舞いを示すことが分かりました。

3.3 間接検出と前方物理実験への影響

• 間接検出（CTA, CMB）: 残留密度の計算誤差は、間接検出信号（$\gamma$線など）の予測強度に直接影響します。nBE に基づく予測は、実際の信号強度を過大評価する可能性があります。
• 前方物理実験（FASER, SHiP, MATHUSLA）: 媒介粒子 $\phi$ の寿命と混合角 $\sin\theta$ の関係は、実験の感度領域を決定します。非平衡効果を考慮した正確な残留密度の算出は、これらの実験で探索されるパラメータ空間の解釈に不可欠です。

4. 結論と学術的意義

1. 非平衡効果の重要性の再確認:
   凍結インシナリオ、特に逐次生成や多成分ダークセクターにおいて、運動量分布の非熱的歪みを無視することは許されません。従来の数密度アプローチ（nBE）は、残留密度を最大で 1 桁以上過大評価する可能性があり、DM の質量や結合定数の制限を誤ったものにする恐れがあります。

2. 計算手法の必要性:
   正確な予測を行うためには、位相空間レベルでの完全なボルツマン方程式（fBE）の解法、あるいは少なくとも有効温度を考慮した cBE の使用が必須であることが示されました。

3. ツール開発への貢献:
   本研究で開発された計算手法は、次期公開版の「DRAKE」コードに実装される予定であり、将来の多成分ダークセクターモデルの研究基盤を提供します。

4. 現象論的意義:
   提案された 2 成分スカラーモデルは、間接検出（銀河中心過剰など）や前方物理実験（長寿命粒子探索）の両方で検出可能なシグナルを持ちます。非平衡効果を正しく扱うことで、これらの実験結果の解釈や、DM 候補としてのモデルの妥当性評価がより精緻に行えるようになります。

総じて、本論文は「凍結インにおける非平衡ダイナミクス」が単なる理論的な興味ではなく、観測可能な物理量（残留密度、検出信号）に決定的な影響を与える重要な要素であることを実証的に示した点に大きな意義があります。