Phase space analysis of an exponential model in $f(Q)$ gravity including linear dark-sector interactions

技術要約：線形ダークセクター相互作用を含む $f(Q)$ 重力における指数型モデルの位相空間解析

問題提起
標準的な $\Lambda$CDM 宇宙モデルは、観測的には成功しているものの、宇宙定数問題や偶然性問題（現在のエポックにおけるダークエネルギーとダークマターの密度の見かけ上の微調整）といった理論的な課題に直面している。DESI（Dark Energy Spectroscopic Instrument）コラボレーションによるデータを含む近年の観測結果は、厳密な宇宙定数ではなく、動的なダークエネルギー流体の可能性を示唆している。これは、修正重力理論、特に非計量スカラー $Q$ の関数である重力作用に基づく $f(Q)$ 重力（拡張対称テレパラレル重力）の探索を動機付けている。様々な $f(Q)$ モデルが提案されてきたが、放射、物質、およびダークエネルギー優占のエポックを含む宇宙膨張の全履歴を再現する上で、$\Lambda$CDM の滑らかな摂動として機能する特定の指数型モデルの妥当性を分析する必要がある。さらに、ダークマター（DM）とダークエネルギー（DE）の間の相互作用を導入することは、偶然性問題を解決するための潜在的なメカニズムを提供するため、$f(Q)$ の枠組み内でのそのような結合系の動力学的解析が必要となる。

手法
著者らは、関数形式 $F(Q) = 2\Lambda e^{-(b\Lambda/Q)^n}$ で定義される指数型 $f(Q)$ モデルを解析するために、動力学系形式を採用している。本研究は、修正フリードマン方程式を自律的な微分方程式系に還元するという Böhm er らによって提案された戦略に従っている。

1. 変数定義: ハッブルパラメータ $H$ に依存する標準的な正規化手法の非自律性を克服するため、著者らは一連の無次元動力学変数を導入している：
   • $x = 8\pi\rho_m / 3H^2$ （物質密度パラメータ）
   • $y = 8\pi\rho_r / 3H^2$ （放射密度パラメータ）
   • $m = H^2 / (\Lambda + H^2)$ （ハッブル率および宇宙定数に関連する変数）
   • 変数 $Z$ （幾何学的修正に関連）は、$m$ を用いて表現され、系を閉じる。
2. 自律系の構築: 展開パラメータ $\eta = \ln a$ に関してフリードマン方程式を書き換えることで、$\{x, m\}$ のペアに関する閉じた方程式系を導出している（$y$ は制約式によって決定される）。
3. 臨界点解析: 著者らは、$\dot{x} = 0$ および $\dot{m} = 0$ となる平衡点を見つけることで系を解析している。これらの点の安定性は、ヤコビ行列と固有値を用いた線形安定性解析によって決定される。
4. 相互作用の拡張: ダークマターとダークエネルギーの間の線形相互作用カーネル $Q = \xi H \rho_{DE}$ を含むようにモデルを拡張している。連続の方程式は、この相互作用項を含めるために修正されている。そして、結合パラメータ $\xi$ に依存する新しい臨界点と安定条件を見つけるために、動力学系を再解析している。

主要な貢献および結果

• 非相互作用モデルの解析:

  • 自律系は、標準的な宇宙論的エポックに対応する3つの臨界点を与える：
    1. 放射優占期 ($R$): 座標 $\{0, 1\}$。固有値 $\{1, 4\}$ は、これが不安定なリペラー（反発点）であることを示す。
    2. 物質優占期 ($M$): 座標 $\{1, 1\}$。固有値 $\{-1, 3\}$ は、これがサドルポイント（鞍点）であり、放射優占からダークエネルギー優占への遷移として機能することを示す。
    3. ド・ジッター・アトラクター ($P_{dS}$): 座標 $\{0, m_{DE}\}$。$m_{DE}$ の値は、パラメータ $b$ に依存する超越方程式の解である。線形安定性解析により、$b > 0$ および $b < 0$ の両方において固有値が負であることが示され、$P_{dS}$ が安定なアトラクターであることが確認された。
  • 宇宙論的パラメータ: 本モデルは、各エポックにおける有効状態方程式 ($w_{eff}$) と減速パラメータ ($q$) を正常に再現している：放射については $w_{eff} = 1/3, q=1$、物質については $w_{eff} = 0, q=1/2$、そしてド・ジッター相については $w_{eff} = -1, q=-1$ である。
  • パラメータ $b$ の影響: パラメータ $b$ は、ド・ジッター相への遷移のタイミングを制御する。$b < 0$ の場合、モデルはファントム的な振る舞い ($w_{eff} < -1$) を示し、$b > 0$ はクインテッセンス的な振る舞いに対応する。ド・ジッター点は $b < 1/6$ の場合にのみ存在し、より大きな値に対しては、点は複素数となり位相空間から消失する。

• 相互作用モデルの解析:

  • 線形相互作用 $Q = \xi H \rho_{DE}$ を含めることで、放射および物質の点は安定性において変化しない（それぞれ不安定およびサドル）。
  • 新しい相互作用ド・ジッター点 ($P_{int}$) が出現する。その座標は $b$ と結合 $\xi$ の両方に依存する。
  • 安定性: 解析されたパラメータ範囲において、$P_{int}$ は安定なアトラクターであることが判明した。
  • 相互作用の物理的意味合い:
    • $\xi > 0$（DE から DM へのエネルギー転送）の場合、アトラクターにおいて物質密度は完全に消失せず ($x_{int} \neq 0$)、ダークエネルギーと共にダークマターが存続する宇宙を意味する。
    • $\xi < 0$（DM から DE へのエネルギー転送）の場合、モデルはアトラクターにおいて負の物質密度パラメータ ($\Omega_m < 0$) を予測し、これは非物理的である。
    • 相互作用は「ファントム・クロッシング」を誘発し、有効状態方程式 $w_{eff}$ が $-1$の障壁を越え、遠い未来において$w_{eff} < -1$ となる。これは、晩期宇宙における零エネルギー条件（NEC）の破れを示唆している。

意義および結論
本論文は、提案された指数型 $f(Q)$ モデルが、厳密な宇宙定数を必要とすることなく、幾何学的な重力の修正から生じる、標準的な宇宙優占シーケンス（放射 $\to$ 物質 $\to$ ダークエネルギー）を定性的に再現する実行可能な宇宙論的シナリオであることを示している。本モデルは $\Lambda$CDM の摂動展開として振る舞い、$H(z)$ データセットを用いた従来の統計的制約と一致している。

線形ダークセクター相互作用の導入は、そのようなモデルが晩期の進化を変化させ、潜在的に偶然性問題を解決できる一方で、結合の符号と大きさによっては、負の密度という非物理的な結果やエネルギー条件（NEC）の破れをもたらす可能性があることを明らかにしている。著者らは、このような動的なダークエネルギーモデルは（DESIなどの）最近の観測的ヒントによって動機付けられているものの、パラメータ空間の拡大（$b$ と $\xi$ の導入）が統計的制約をより困難にしていると結論付けている。本研究は、これらの拡張重力モデルの妥当性を判断するために、今後のデータリリースを用いたさらなる制約の必要性を強調している。