Relativistic N-body simulations with mass-varying neutrinos

技術要約：質量変化ニュートリノを伴う相対論的N体シミュレーション

問題提起
近年の宇宙論データ、特にDESIとPlanckおよびACTを組み合わせた解析は、ニュートリノ質量の総和（$\Sigma m_\nu < 0.0642$ eV）に関する上限を示唆しており、これは局所的なニュートリノ振動実験から得られる下限（$\Sigma m_\nu > 0.06$ eV）と緊張関係にあります。この不一致は、ニュートリノの質量が宇宙時間とともに変化するモデルの調査を動機付けています。具体的には、本論文では、ダークエネルギーとして機能するクインテッセンス場（「コスモン」）に結合したニュートリノを扱うモデルを検証しています。このようなモデルはダークエネルギーの同時性問題を解決できる可能性がありますが、従来の数値研究では、宇宙論的スケールにおける非線形構造である「ニュートリノ・ランプ（塊）」が一般的に形成されることが示されており、これが大規模構造の観測結果と矛盾する可能性があります。さらに、これまでの解析は線形的な取り扱いや非相対論的な近似に依存していることが多く、重要なバックリアクション効果や、異方性ストレスおよび有効ダークエネルギー状態方程式の進化といった相対論的現象を見逃している可能性があります。

手法
著者らは、新しいGPU加速N体および格子シミュレーションコードであるGAsGRDを開発し、活用しました。このコードは、粒子と場の完全な相対論的力学を維持しながら、弱場展開（計量ポテンシャル $\Phi, \Psi \sim \epsilon$）においてアインシュタイン方程式を解くgevolutionフレームワークの拡張版です。

主な技術的実装は以下の通りです：

• 相対論的フレームワーク： コードはアインシュタイン・フレームで動作しますが、結合種（ニュートリノ）については、質量が $m_\nu = A(\phi)\tilde{m}_\nu$ として変化する種依存のジョルダン・フレームで扱います。ジョルダン・フレームにおける計量ポテンシャルは、大きな共形因子変化（$\Delta A$）に対応するために非摂動的に扱われますが、アインシュタイン・フレームのポテンシャルは効率化のために線形化されたままとなります。
• 力学： コードは、スカラー媒介力を含む粒子の測地線方程式を解き、4次ルンゲ＝クッタ積分器を用いてスカラー場（コスモン）を更新します。また、CDM背景に対して急速に進化するニュートリノを解像するために、ニュートリノに対する適応型タイムステップ・スキームを採用しています。
• バックリアクション： シミュレーションは、宇宙膨張率へのバックリアクション効果を考慮するため、体積平均エネルギー密度を用いてフリードマン方程式をオンザフライで解きます。これは線形ソルバーでしばしば無視される機能です。
• 初期条件： 初期条件は、ゼルドビッチ近似を用い、ニュートリノの速度分散のフェルミ・ディラック分布をサンプリングするCLASS-MaVaN（質量変化ニュートリノのためのCLASSの修正版）を使用して生成されます。
• モデルパラメータ： 指数関数的ポテンシャル $V(\phi) \propto e^{-\alpha \phi}$ および指数関数的結合 $A(\phi) = e^{\beta \phi}$ に焦点を当てています。2つのシミュレーション・スイートが実行されました：
  • 小規模ボックス・スイート (Ai)： ニュートリノ・ランプのサブストラクチャの解像度要件をテストするために設計されました（$L_{box} = 50-500 \, h^{-1}\text{Mpc}$）。
  • 大規模ボックス・スイート (Ei)： 宇宙論的な生存可能性とバックリアクション効果を探索するために設計されました（$L_{box} = 1-10 \, h^{-1}\text{Gpc}$）。パラメータは、$z=0$ における局所的なニュートリノ質量境界を満たしつつ、非線形補正を考慮するように調整されています。

主な結果

1. ニュートリノ・ランプ： シミュレーションは、非線形なニュートリノ・ランプの形成を裏付けています。しかし、従来の静的な記述とは異なり、これらのランプは非常に動的で、回転しており、コスモン場の振動によって周期的な集積と希薄化を起こしていることが判明しました。
2. 解像度と収束： ニュートリノ・ランプのサブストラクチャを解像するには、高い格子解像度（$N_{grid} \gtrsim 128$）と十分な粒子数が必要です。粗いシミュレーションでは、ランプの正しい方位や位相空間を捉えることができません。
3. バックリアクションと膨張史： 非線形バックリアクション効果は、背景の膨張史を著しく変化させます。線形的な取り扱いにおいて生存可能に見えるモデルも、非線形効果を含めると排除されるか、修正を余儀なくされます。シミュレーションによれば、これらのバックリアクション効果により、有効ハッブルパラメータ $H_0$ は $\nu\Lambda$CDMよりも4〜6%大きくなることがあります。
4. ファントム交差と状態方程式： 有効ダークエネルギー状態方程式（$w_{DE,eff}$）は振動挙動を示し、ファントム境界（$w = -1$）を何度も横切ります。この挙動は、スカラー場とニュートリノ間のエネルギー転移に起因します。著者らは、標準的な $w(a) \text{CDM}$ モデルを仮定している観測者は、仮定された物質密度 $\Omega_m$ に依存して、たとえ基礎となる物理が異なっていても、ファントム挙動を推論してしまう可能性があると指摘しています。
5. 構造形成： クインテッセンス場の存在とニュートリノ質量の変化は、$\nu\Lambda$CDMと比較して冷たい物質の集積を抑制しますが、その抑制の程度はモデルパラメータ（特に $\alpha$）によって異なります。
6. 相対論的効果： シミュレーションは、顕著な異方性ストレス（$X = \Phi - \Psi$）を明らかにしました。これは、特定のモデルにおいて小スケールでバルディアン・ポテンシャルの最大90%に達します。これは、非線形なコスモン場と半相対論的なニュートリノによって引き起こされます。

意義と主張
本論文は、GAsGRDが、線形ソルバーが見落とす非線形および相対論的効果を捉えることにより、質量変化ニュートリノ・モデルを観測に対して厳密にテストするための枠組みを提供すると主張しています。著者らは以下を強調しています：

• 元の一定結合モデルは、初期の散逸を通じて、観測上の問題となるニュートリノ・ランプの形成を遅らせるか、あるいは弱めることができます。
• ニュートリノの特性、ダークエネルギーの動力学、および重力ポテンシャルの間のスリップ（異方性ストレス）の観測を通じて、これらのモデルを $\Lambda$CDMと区別するためには、相対論的シミュレーションが不可欠です。
• ファントム交差や大規模なニュートリノ構造を含む豊かな現象論は、パラメータ空間のさらなる探索を動機付けるものです。
• 現在のパラメータ選択は冷たい物質の集積における抑制を示していますが、これらのモデルは、将来の観測プロジェクト（Euclid、LSST、Romanなど）が、このコードが予測可能な特定のシグネチャ（異方性ストレスや、重力レンズ効果および積分ザックス・ヴォルフェ効果などの二次CMB効果）を制約できれば、依然として有力な候補であり続けます。

著者らは、これらのモデルを標準的な宇宙論から決定的に識別するためには、このような非線形シミュレーションから構築されたエミュレータを用いたパラメータ空間の系統的な研究が必要であると結論付けています。