Constraining a $f(R, L_m)$ Gravity Cosmological Model with Observational Data

本研究では、多様な観測データセットを用いたベイズMCMC解析により、空間的に平坦な$f(R, L_m)$重力モデルを制約し、それが標準的な$\Lambda$CDM宇宙論に代わる実行可能な選択肢となり得る一方で、より高い推論された$H_0$の値を通じてハッブル・テンションを解決する可能性があることを見出した。

要約

宇宙を、巨大で膨張する風船だと想像してみてください。長い間、科学者たちはこの風船が一定のペースで膨らんでいるか、あるいはガソリンが切れかかった車のように、むしろ減速していると考えてきました。しかし、1990年代後半、奇妙なことが発見されました。風船はただ膨らんでいるだけでなく、加速しているのです。まるで、誰かが突然アクセルを踏み込んだかのようです。

この「キック（加速）」を説明するために、標準的な科学モデル（ΛCDMと呼ばれます）は、「ダークエネルギー」と呼ばれる目に見えない、謎めいた力が風船を押し広げていると説明します。これは、風船の中に隠れたモーターが入っていると言うようなものです。

この論文は、異なるアイデアを提案しています。隠れたモーターを追加する代わりに、著者たちは、**「道路のルール（重力の法則）」**自体が、私たちが考えているものとは少し異なっている可能性を示唆しています。彼らは、$f(R, L_m)$ 重力と呼ばれる理論を探求しています。

以下は、彼らが何を行い、何を見出したのかについての簡単な解説です。

1. 新しいルールブック：幾何学と物質の混合

アインシュタインの古い理論では、重力はトランポリンのようなものです。重い物体（物質）は布地を曲げ、他の物体はその方へと転がっていきます。

• 標準的な見方： 布地（幾何学）と物体（物質）は相互作用しますが、それぞれ自分のレーンに留まっています。
• 新しい見方 ($f(R, L_m)$)： 著者たちは、布地と物体がもっとしっかりと「手を繋いでいる」ことを示唆しています。それらは直接結合（カップリング）しています。もしトランポリンの布地がその上のボウリングの玉に囁きかけ、ボウリングの玉もまた布地の伸び方を変えるように囁き返せるとしたら、どうなるかを想像してみてください。

この「囁き合い」（結合）が新しい力を生み出します。これには、謎の「ダークエネルギー」というモーターを必要としません。加速は、この新しい重力の仕組みから自然に生じるのです。

2. 探偵の仕事：証拠の検証

著者たちは単に方程式を書いたのではありません。彼らは、自分たちの新しい理論が宇宙に残された手がかりに適合するかどうかをチェックする探偵として行動しました。彼らは4つの主要な「証拠」を用いました。

• 宇宙年代計（Cosmic Chronometers）： 古い星の年齢を測定することで、異なる時期に宇宙がどれくらいの速さで膨張していたかを確認しました。
• 超新星（Pantheon+）： 爆発する星を「標準光源」として使い、宇宙の端々までの距離を測定しました。
• BAO（バリオン音響振動）： 銀河の分布に見られる「化石化した」波紋（木の年輪のようなもの）を観察しました。
• CMB（宇宙マイクロ波背景放射）： ビッグバンの名残であり、初期宇宙のスナップショットを与えてくれます。

彼らは強力なコンピュータ手法（MCMC）を使用して、これらのデータに対して自分たちの理論の何百万ものバリエーションをテストし、どのバージョンが最もよく適合するかを確認しました。

3. 結果：より良い適合か？

彼らの「新しいルールブック」は、標準モデルと比較して何を予測したのでしょうか。

• 膨張速度 ($H_0$)： 科学界には「ハッブル・テンション」と呼ばれる有名な不一致があります。宇宙の膨張速度を測る方法の一つは低い数値を示し、もう一つは高い数値を示します。標準モデルは低い数値の方に傾いています。
  • 論文の主張： 彼らの新しいモデルは、自然に高い速度を予測し、これは「高い数値」の測定結果により良く一致します。これは、彼らの理論が他の部分を壊すことなく、この不一致を解決できる可能性があることを示唆しています。
• 「キック」の瞬間： 両方のモデルは、宇宙が減速から加速へと切り替わったのが、約70億年から80億年前（赤方偏移がおよそ0.7から0.8の時）であるという点で一致しています。彼らのモデルはこの切り替えが0.79で起きたと予測しており、これはデータに完璧にフィットします。
• 宇宙の年齢： 彼らは宇宙の年齢が約133.4億年であると計算しました。これは標準モデルの推定値である130.6億年に非常に近い数値です。

4. 判定：それはより優れているのか？

著者たちは、統計的な「スコアカード」（AICおよびBICと呼ばれます）を用いて、彼らのモデルを標準モデルと比較しました。

• これは、2つのレシピを比較することに似ています。標準的なレシピはシンプルです（材料が少ない）。新しいレシピには、もう一つの材料（新しい重力の結合）が追加されています。
• 通常、材料を追加すると、それが劇的に美味しくならない限り、レシピの評価は「悪化」します。
• 結果： 彼らの新しいレシピは、わずかに「美味しく」なりました。統計的なスコアは、標準モデルよりもほんの少しだけ低く（＝より良く）なりました。これは、彼らの理論が、謎の「ダークエネルギー」モーターを必要とせずに、データを同等か、あるいはわずかに優れたレベルで説明できる、実行可能で一貫性のある代替案であることを意味します。

まとめ

著者たちは、重力と物質が直接対話する新しい理論を構築しました。彼らはこれを、現在得られる最高のデータを用いて宇宙の歴史に対してテストしました。その結果、以下のことが分かりました。

1. 謎の「ダークエネルギー」を必要とせずに、なぜ宇宙が加速しているのかを説明できる。
2. 宇宙の膨張速度の測定値がなぜ食い違うのかというパズル（ハッブル・テンション）を解く助けとなる。
3. 彼らのモデルは、現在の標準モデルとほぼ同等、あるいはわずかに優れたレベルで観測データに適合する。

要するに、彼らは宇宙の車を走らせる新しい方法を見つけました。それは古い方法と同じくらいスムーズに感じられますが、古い地図では説明できなかった「道の上の凹凸」を説明できるかもしれないのです。

技術概要

技術要約：観測データを用いた $f(R, L_m)$ 重力宇宙論モデルの制約

問題提起
標準的な $\Lambda$CDM 宇宙論モデルおよび一般相対性理論（GR）は、宇宙定数や未知のダークエネルギー成分を導入することなく、宇宙の後期における加速膨張を説明する上で重大な課題に直面している。さらに、「ハッブル・テンション」として知られる持続的な不一致が存在する。これは、初期宇宙の測定値（例：Planck 2018 CMB データによる $H_0 \approx 67.4$ km/s/Mpc）と、局所的な後期の決定値（例：SH0ES による $H_0 \approx 73$ km/s/Mpc）との間の乖離である。曲率と物質の非最小結合を含む修正重力理論は、これらの問題を解決するための潜在的な理論的枠組みを提供する。本研究では、曲率 $R$ と物質ラグランジアン密度 $L_m$ の両方に依存する $f(R, L_m)$ 重力が、$\Lambda$CDM に代わる選択肢としてどの程度生存可能かを調査する。

手法
著者らは、$f(R, L_m)$ の枠組みの中で、空間的に平坦なフリードマン・ルメートル・ロバートソン・ウォーカー（FLRW）宇宙論モデルを構築した。

1. 理論的定式化: 本研究では、$f(R, L_m) = \alpha R + L_m^\beta + \gamma$ という特定の関数形式を採用した（ここで $\alpha, \beta, \gamma$ は自由パラメータである）。計量に関する作用を変化させることで、著者らは修正された場の方程式を導出した。状態方程式 $p = (1-n)\rho$ を持つ完全流体を仮定することで、$H(z)$ に関する解を解析的に求め、以下の式を得た：
   $$H(z) = H_0 \sqrt{(1-\lambda) + \lambda(1+z)^{3(1+w)}}$$
   ここで、$\lambda$ および $w$ は、モデル定数と状態方程式パラメータ $n$ から導かれる無次元パラメータである。
2. 観測データ: モデルパラメータ（$H_0, \lambda, w$）は、以下の包括的なデータセットを用いて制約された：
   • コズミック・クロノメーターからの 35 個のハッブルパラメータ（$H(z)$）の測定値。
   • Pantheon+ コンピレーションを含む 1701 個の Type Ia 超新星ライトカーブ（統計的および系統的誤差のための全共分散行列を含む）。
   • DESI 初期リリースからのバリオン音響振動（BAO）データ。
   • 宇宙マイクロ波背景放射（CMB）のシフトパラメータ。
3. 統計解析: 著者らは、パラメータ空間を探索するために emcee サンプラーを用いたベイズ・マルコフ連鎖モンテカルロ（MCMC）法を採用した。総カイ二乗統計量（$\chi^2_{total}$）を最小化して最良適合値を求め、68% 信頼区間を算出した。モデルの性能は、赤池情報量基準（AIC）およびベイズ情報量基準（BIC）を用いて標準的な $\Lambda$CDM モデルと比較評価された。

主な貢献と結果

• 解析的導出: 従来のアンザッツに基づいたアプローチに依存していた先行研究とは異なり、本研究では $f(R, L_m)$ 重力の基礎的な場の方程式から直接かつ厳密に、宇宙論的解および $H(z)$ の具体的な形式を導出した。
• パラメータ制約: MCMC 解析により、$f(R, L_m)$ モデルにおける以下の 1$\sigma$ 信頼レベルでの最良適合値が得られた：
  • $H_0 = 73.750^{+0.159}_{-0.156}$ km/s/Mpc
  • $\lambda = 0.262^{+0.007}_{-0.007}$
  • $w = -0.005^{+0.001}_{-0.001}$
• ハッブル・テンションの解決: 推定された $H_0$ の値は、Planck 2018 $\Lambda$CDM 値（$67.4$ km/s/Mpc）よりも有意に高く、局所的な距離梯子による測定値（例：SH0ES）と密接に一致している。著者らは、$f(R, L_m)$ 重力における追加の自由度が、後期の観測と矛盾することなくハッブル・テンションを緩和する潜在的なメカニズムを提供していることを示唆している。
• 宇宙の進化: モデルは、減速から加速への遷移赤方偏移 $z_t \approx 0.79$ を予測しており、これは観測データと一致している。推定される宇宙の年齢は $13.34$ Gyr であり、$\Lambda$CDM の推定値である $13.06$ Gyr と同程度である。
• モデル比較: AIC および BIC を用いた統計解析は、$f(R, L_m)$ モデルが追加の自由パラメータを持っているにもかかわらず、$\Lambda$CDM よりも結合データセットに対してわずかに優れた適合を示すことを示している。差は緩やかであり（$\Delta$BIC $\approx -5.78, \Delta$AIC $\approx -11.30$）、修正重力シナリオに対する中程度の統計的選好を示唆している。
• 物理定数: 現在の物質密度 $\rho_0 = 0.534 \times 10^{-30}$ g/cm$^3$ および $n=1$（ダスト時代）を仮定することで、固有のモデル定数が決定された：$\beta \approx 1.005, \alpha \approx 4.51 \times 10^5, \gamma \approx -1.14 \times 10^{-29}$。

意義と主張
本論文は、$f(R, L_m)$ 重力シナリオが、標準的な $\Lambda$CDM モデルに代わる実現可能かつ一貫した選択肢であることを主張している。その主な意義は、以下の点にある：

1. $\Lambda$ なしでの加速の説明: 本モデルは、曲率と物質の結合によって駆動される後期宇宙の加速を再現し、宇宙定数の必要性を排除する。
2. ハッブル・テンションへの対処: この枠組みは自然に高いハッブル定数の値をサポートし、初期宇宙と後期宇宙の測定値の間の乖離に対する潜在的な解決策を提供する。
3. 一貫性の維持: モデルは現在の観測制約（CMB, BAO, SNe Ia）と整合しており、減速から加速へのシームレスな遷移を予測する。

著者らは、$\Lambda$CDM に対する統計的選好は控えめであるものの、$f(R, L_m)$ の枠組みは、特にダークエネルギーの性質やハッブル・テンションに関する未解決の問いを解決するための有望な道筋を提示していると結論付けている。今後の課題としては、宇宙論的摂動、重力波伝播の研究、および今後の高精度サーベイ（Euclid, Roman, Rubin）への本モデルの適用が挙げられている。