Exponential $f(R)$ cosmology with massive neutrinos as a dynamical dark energy framework

技術要約：質量を持つニュートリノを伴う指数関数型 $f(R)$ 宇宙論

問題提起
標準的な $\Lambda$CDM 宇宙論モデルは成功を収めているものの、持続的な理論的および観測的課題に直面している。理論的な問題には、宇宙定数 $\Lambda$ に関連する微調整問題や同時性問題、および確認されたダークマター粒子の欠如が含まれる。観測的には、顕著なテンションが存在する。最も注目すべきは「ハッブル・テンション」であり、これは局所的な距離梯子による測定値（$H_0 \approx 73.04$ km s$^{-1}$ Mpc$^{-1}$）と、$\Lambda$CDM フレームワーク内での宇宙マイクロ波背景放射（CMB）からの初期宇宙の推論（$H_0 \approx 67.36$ km s$^{-1}$ Mpc$^{-1}$）との間の $\sim 5\sigma$ の不一致である。さらに、DESI コラボレーションによる最近のバリオン音響振動（BAO）データは、$\Lambda$CDM からの逸脱を示唆しており、動的なダークエネルギーの存在を暗示している。

二次的なテンションは、全ニュートリノ質量（$\sum m_\nu$）に関するものである。ニュートリノ振動実験は $\sum m_\nu \gtrsim 0.06$ eV という下限を確立しているが、$\Lambda$CDM モデル内の宇宙論的データは、しば_しばこの閾値を下回る上限値（例：DESI による $\sum m_\nu < 0.064$ eV）を導き出す。現象論的な $w_0w_a$CDM モデルはこれらの境界を緩和できる可能性があるが、強固な理論的基礎を欠いている。本研究では、理論的に動機付けられた修正重力、具体的には指数関数型 $f(R)$ 重力が、質量を持つニュートリノと組み合わさることで、標準モデルとの高赤方偏移での整合性を維持しつつ、ハッブル・テンションとニュートリノ質量制約を同時に解決できるかどうかを調査する。

手法
著者らは、以下のラグランジアンで定義される一般化された指数関数型 $f(R)$ 重力モデルを用いている：
$$f(R) = R - 2\Lambda \left\{ 1 - \exp\left[ -\beta \left( \frac{R}{2\Lambda} \right)^\alpha \right] \right\}$$
ここで、$\alpha$ と $\beta$ は、高曲率における一般相対性理論（GR）から低曲率における修正重力への遷移を制御する無次元パラメータである。本モデルは、冷たいダークマター、バリオン、放射、および質量を持つニュートリノからなる物質セクターを含む。ニュートリノのエネルギー密度は、相対論的領域から非相対論的領域への遷移を考慮し、フェルミ・ディラック分布に対する位相空間積分を通じて計算される。

宇宙論的ダイナミクスは、空間的に平坦な FLRW メトリックから導かれる修正フリードマン方程式から導出される。数値解析を容易にするため、著者らは無次元のハッブルパラメータ $E(z) = H(z)/H_0^*$ を用いて、$f(R)$ パラメータを等価な $\Lambda$CDM パラメータ（$H_0^*, \Omega_i^*$）にマッピングしている。$E(z)$ と無次元のリッチスカラー $\mathcal{R}(z)$ に関する結合微分方程式系は、高赤方偏移の初期条件（モデルが $\Lambda$CDM を回収する点）から $z=0$ まで数値的に積分される。

ベイズ的マルコフ連鎖モンテカルロ（MCMC）解析が、emcee パッケージを用いて、パラメータ空間 $(H_0^*, \alpha, \beta, \omega_b^*, \omega_{bc}^*, \omega_\nu^*)$ を制約するために実行される。解析には、以下の4つの観測データセットが用いられる：

1. 宇宙年代計 (CC): $H(z)$ を測定した31個のデータ点。
2. DESI DR2 BAO: 6つの赤方偏移ビンにおける $D_M(z)/r_d$ および $D_H(z)/r_d$ の測定値。
3. CMB: 音響スケール $\theta_*$ および物理密度 $\omega_b, \omega_{bc}$ に関する圧縮された Planck 2018 尤度。
4. Pantheon+ (Pan+): 1550 個の Type Ia 超新星の距離モジュラス。

結果は、赤池情報量基準（AIC）およびベイズ情報量基準（BIC）を用いて、統計的選好度を評価するために、標準的な $\Lambda$CDM モデルおよび現象論的な $w_0w_a$CDM モデルと比較される。

主な貢献と結果
本研究は、最新の DESI DR2 データを用いて、一般化された指数関数型 $f(R)$ モデルのパラメータ（$\alpha, \beta$）および全ニュートリノ質量（$\sum m_\nu$）に対する初の結合制約を提供する。

• モデルの生存可能性: 指数関数型 $f(R)$ モデルは現在の観測と矛盾しない。安定性条件 $f_R > 0$ および $f_{RR} > 0$ は、全事後パラメータ空間および赤方偏移範囲にわたって満たされており、ゴーストやタキオン不安定性の不在を裏付けている。
• ハッブル・テンション: 本モデルはハッブル・テンションを部分的に緩和する。Pantheon+ 超新星を含めた場合、$f(R)$ モデルから推論される $H_0$ は $70.89 \pm 0.19$ km s$^{-1}$ Mpc$^{-1}$ である。この値は、$\Lambda$CDM ($71.43 \pm 0.13$) や $w_0w_a$CDM ($71.13 \pm 1.20$) と比較して、局所的な測定値とのテンションを軽減しているが、フルデータセットの組み合わせにおいては、$f(R)$ と $\Lambda$CDM の間に $2.3\sigma$ の不一致が依然として残っている。
• ニュートリノ質量制約: $f(R)$ フレームワークに質量を持つニュートリノを組み込むことで、$\sum m_\nu$ の上限値が $\Lambda$CDM と比較して緩やかに緩和される。Pantheon+ データを用いた場合、$1\sigma$ 上限は $\Lambda$CDM の $< 0.0126$ eV や $w_0w_a$CDM の $< 0.0613$ eV に対し、$f(R)$ では $\sum m_\nu < 0.0290$ eV となる。これは緩和ではあるものの、上限値は依然として振動実験の下限（$\sim 0.06$ eV）を下回っており、テンションは完全には解決されていない。
• 統計的選好度:
  • Pantheon+ データを含まない場合、$\Lambda$CDM が統計的に好まれる。
  • Pantheon+ データの追加により、AIC および BIC メトリックに従うと、$w_0w_a$CDM ($\Delta$AIC = 43.75) および指数関数型 $f(R)$ ($\Delta$AIC = 20.46) は、共に $\Lambda$CDM に対して「非常に強い」および「強い」証拠を示す。
• パラメータの縮退: 解析は、修正パラメータ（$\alpha, \beta$）と宇宙論的パラメータ（$H_0, \Omega_m$）の間の強い縮退を明らかにしており、これは修正重力が物質密度や膨張史の変化を部分的に補償できることを示している。$\beta$ の事後分布は非ガウス的であり、これは指数関数的補正の非線形な役割を反映している。

意義と主張
本論文は、修正重力とニュートリノ物理学の相互作用が、宇宙論的テンションを調査するための実行可能かつ理論的に堅牢な枠組みを提供することを主張している。現象論的なダークエネルギーのパラメータ化とは異なり、指数関数型 $f(R)$ モデルは、重力のラグランジアン修正に根ざした強固な理論的基礎を提供している。

著者らは、本モデルがハッブル・テンションやニュートリノ質量の不一致を完全に解決するものではないものの、重力セクターを修正しつつ質量を持つニュートリノを含めるという複合的なアプローチが、$\Lambda$CDM 単独よりも観測データへの統計的な適合性を向上させ、ニュートリノ質量の制約をより効果的に緩和できることを示していると強調している。本研究は、標準的な宇宙論とモデルを判別する上で、低赤方偏移のデータセット（DESI 等の SNe Ia）の重要性を浮き彫りにしている。著者らは、現在の解析が背景宇宙論に限定されているため、モデルの妥当性を完全に評価するには、将来的な摂動効果（弱重力レンズ効果、赤方偏移空間歪みなど）を含める必要があると結論づけている。