Nonlinear Matter Power Spectrum from relativistic $N$-body Simulations: $\Lambda_{\rm s}$CDM versus $\Lambda$CDM

この論文は、一般相対論に基づく相対論的 N 体シミュレーションを用いて、符号反転する宇宙定数を持つ$\Lambda_{\rm s}$CDM 模型と標準的な$\Lambda$CDM 模型を比較し、前者が非線形物質パワースペクトルに観測可能かつ$\sigma_8$や$S_8$のスケール不変な変化では説明できない特徴的な増幅パターンを生み出すことを示すことで、宇宙の星形成率のピークと整合的な重力学的な予測を提供している。

要約

この論文は、宇宙の進化を説明する現在の標準モデル（$\Lambda$CDM モデル）に、少しだけ「不思議な魔法」を加えた新しいシナリオ（$\Lambda_s$CDM モデル）が、実際の宇宙の構造形成にどのような影響を与えるかを、スーパーコンピュータを使った大規模なシミュレーションで検証したものです。

専門用語を避け、日常の比喩を使ってわかりやすく解説します。

1. 物語の舞台：宇宙という「巨大なクレープ」

まず、宇宙を想像してください。そこには「ダークマター」という見えない土台があり、その上に銀河や銀河団（星の集まり）が乗っています。
これまでの標準モデル（$\Lambda$CDM）では、宇宙の膨張を加速させる「ダークエネルギー」という力は、昔から一定の強さで働き続けてきたと考えられていました。まるで、一定の力で押され続けるクレープのようなものです。

しかし、この論文は**「実は、その力は過去に一度、方向を変えていたのではないか？」**と問いかけます。

2. 魔法のスイッチ：「 sign-switching（符号の切り替え）」

新しいモデル（$\Lambda_s$CDM）では、宇宙の歴史の途中（約 100 億年前、赤方偏移$z \approx 2$の頃）に、ダークエネルギーの正体が**「プラスからマイナスに、そしてまたプラスに戻る」**という劇的な変化を起こしたと仮定しています。

• マイナスの時代（AdS 相）：
  宇宙の膨張を「少しだけブレーキ」かけるような、あるいは重力を少し強くするような働きをしました。

  • 比喩： 宇宙というクレープが膨らむのを、一時的に「少しだけ押さえつけ」た状態です。
  • 結果： 押さえつけられると、中の物質（ダークマター）がより強く引き寄せられやすくなります。まるで、風船を少し潰すと中の空気が集まりやすくなるように、銀河の種（構造）が急成長しました。

• プラスの時代（dS 相）：
  その後、スイッチが戻り、いつものように宇宙の膨張を加速させました。

  • 比喩： 再び風船を膨らませる力を強くしました。
  • 結果： 成長が少し抑制されますが、すでに「マイナスの時代」に急成長した銀河の種は消えません。その**「急成長の痕跡」がそのまま残ります。**

3. シミュレーションの結果：「波のうねり」

研究者たちは、この「スイッチが入った宇宙」と「通常の宇宙」をスーパーコンピュータでシミュレーションし、銀河の分布（パワースペクトル）を比較しました。

そこで発見されたのは、**「特定のサイズの銀河団にだけ、目立つ盛り上がり（ Crest ）」**でした。

• 通常の宇宙（$\Lambda$CDM）： 銀河の分布は滑らかです。
• 新しい宇宙（$\Lambda_s$CDM）： 銀河団（グループや小さなクラスター）のサイズに相当する場所（波長）で、20% 近くも銀河が密集しているという「うねり」が見つかりました。

重要なポイント：
この「うねり」は、単に全体的に銀河が増えただけ（$\sigma_8$という値を変えるだけ）ではありません。

• 時限爆弾のような特徴： この盛り上がりは、宇宙の歴史の中で「いつ、どこで」スイッチが入ったかによって、その位置が移動します。
  • 昔（赤方偏移$z=2$頃）は、小さな銀河の集まりでピークでした。
  • 現在（$z=0$）になると、そのピークは少し大きな銀河団のサイズに「漂流」して現れます。
  • これは、**「宇宙の成長履歴そのものが、銀河の分布に刻み込まれている」**ことを意味します。

4. なぜこれが重要なのか？「宇宙の真昼」へのヒント

この研究で注目すべきは、その「盛り上がり」が起きるタイミングと場所です。

• タイミング： 宇宙の星が最も活発に作られていた時期（「宇宙の正午」と呼ばれる$z \approx 1.5 \sim 2$の時代）と一致します。
• 場所： 銀河団や小さな銀河グループのサイズです。

比喩：
もし、この新しいモデルが正しければ、宇宙の歴史書（銀河の分布）には、**「ある時期に、星を作る工場が突然フル稼働した」という証拠が、銀河団のサイズに刻まれていることになります。
現在の標準モデルでは、この「急激な成長」を自然に説明するのが難しいのですが、この「スイッチモデル」なら、「ブレーキが一時的にかかったことで、逆に星の材料がギュッと集まりやすくなった」**と説明がつくのです。

5. 今後の展望：「証拠」を探す旅

この論文は、単なる理論的な話ではありません。
「もしこのモデルが正しければ、特定のサイズの銀河団の観測データに、この『20% ほどの盛り上がり』が見えるはずだ」と予測しています。

• 検証方法： 今後の望遠鏡（Euclid や LSST など）を使って、弱い重力レンズ効果や銀河団の数を精密に測れば、この「盛り上がり」があるかどうか、すぐにチェックできます。
• 意義： もしこの「盛り上がり」が見つければ、宇宙のダークエネルギーは一定ではなく、過去に劇的な変化（スイッチ）を起こしていたという、画期的な発見になります。逆に、見つからなければ、このモデルは否定されます。

まとめ

この論文は、**「宇宙のダークエネルギーは、過去に一度、正負を逆転させるスイッチを踏んだのではないか？」という大胆な仮説を、スーパーコンピュータで検証し、「その結果、特定のサイズの銀河団に『成長の傷跡』が残っているはずだ」**と予測したものです。

それはまるで、**「宇宙というクレープの表面に、過去に誰かが指で強く押した跡（盛り上がり）が残っている」**と発見し、その跡を調べることで、宇宙の秘密（ダークエネルギーの正体）を解き明かそうとする挑戦です。

技術概要

この論文は、標準的な宇宙モデル（$\Lambda$CDM）と、宇宙定数が符号を反転する（負から正へ）「$\Lambda_s$CDM」モデルを比較し、相対論的 N 体シミュレーションを用いて非線形領域における物質パワースペクトルの違いを明らかにした研究です。以下に、問題設定、手法、主要な貢献、結果、そして意義について詳細な技術的サマリーを記述します。

1. 問題設定 (Problem)

現在の標準宇宙モデルである$\Lambda$CDM は、宇宙マイクロ波背景放射（CMB）などの初期宇宙の観測と非常に良く一致していますが、ハッブル定数（$H_0$）や構造形成の振幅（$S_8$）に関する「宇宙論的緊張（tensions）」が近年顕在化しています。特に、CMB からの推定値と局所距離梯子による測定値の間の不一致（$H_0$ 緊張）や、弱い重力レンズ観測と CMB 予測の間の不一致（$S_8$ 緊張）は、標準モデルの不完全性を示唆しています。

これらの緊張を解決するための候補として、$\Lambda_s$CDM モデル（符号反転する宇宙定モデル）が提案されています。このモデルでは、赤方偏移 $z^\dagger \sim 2$ 付近で宇宙定数 $\Lambda_s$ が負（Anti-de Sitter, AdS 的）から正（de Sitter, dS 的）へ急激に反転すると仮定します。

• AdS 相 ($z > z^\dagger$): 宇宙定数が負のため、ハッブル摩擦が減少し、密度揺らぎの成長が促進されます。
• dS 相 ($z < z^\dagger$): 宇宙定数が正となり、通常の加速膨張に戻りますが、転移後の膨張率は$\Lambda$CDM よりも大きくなります。

しかし、このモデルが**非線形領域（構造形成の最終段階）**でどのような物質パワースペクトルの特徴を生み出すか、特に観測可能なスケールで$\Lambda$CDM とどのように区別できるかについては、相対論的 N 体シミュレーションを用いた詳細な検証が不足していました。

2. 手法 (Methodology)

本研究では、一般相対性理論に基づいた相対論的 N 体シミュレーションコード**「gevolution」**を使用しました。

• シミュレーション設定:
  • モデル: 急激な符号反転を近似した「瞬間的$\Lambda_s$CDM」モデルと、標準的な$\Lambda$CDM モデルを比較。
  • パラメータセット: 2 つの異なるデータセットに基づく最良適合パラメータを使用。
    1. Planck-only: Planck CMB データのみの解析。
    2. Full: Planck + BAO + Pantheon+SH0ES + KiDS-1000 の結合データセット。
  • 物理成分: バリオンと冷たい暗黒物質（CDM）を単一の粒子集団として共進化させ、質量を持つニュートリノ（$\sum m_\nu = 0.06$ eV）は線形摂動として取り扱いました。
  • 計算条件: 体積 $(2080 \text{ Mpc}/h)^3$、空間分解能 $0.5 \text{ Mpc}/h$、4160 個のグリッド点。
• 検証:
  • 線形領域では、コード「CLASS」のニュートンゲージ出力と gevolution の出力を比較し、背景進化の実装が正しいことを確認しました。
  • 非線形領域では、絶対パワースペクトルだけでなく、両モデルの比 $P_{\Lambda_s\text{CDM}} / P_{\Lambda\text{CDM}}$ を解析することで、解像度の影響を低減し、特徴的なシグナルを抽出しました。

3. 主要な貢献と結果 (Key Contributions & Results)

A. 非線形パワースペクトルにおける特徴的な「山（Crest）」の発見

$\Lambda_s$CDM モデルでは、$\Lambda$CDM に対して物質パワースペクトルの比に、赤方偏移に依存する明確な「山（crest）」構造が現れることが確認されました。

• AdS 相での増幅: 転移前（$z > z^\dagger$）の AdS 相では、ハッブル摩擦の減少により摂動成長が促進され、高波数側（$k \sim 1-3 \, h \text{ Mpc}^{-1}$）でパワースペクトルが増大します。
• 転移後の抑制と残存: 転移後（$z < z^\dagger$）は dS 相となり、より高い膨張率によって成長が抑制されますが、AdS 相で増幅された揺らぎは完全に消去されず、非線形結合を通じて物理的スケールが大きい方へ移動します。
• 特徴的な形状:
  • 転移直後 ($z \sim 2$): 比 $P_{\Lambda_s\text{CDM}} / P_{\Lambda\text{CDM}}$ は $k \simeq 1-3 \, h \text{ Mpc}^{-1}$ 付近でピークに達し、増幅度は 20-25% に達します。
  • 現在 ($z=0$): このピークは時間とともに大きな物理的スケール（小さな $k$）へ移動し、$k \simeq 0.6-1.0 \, h \text{ Mpc}^{-1}$ 付近で 15-20% の持続的な増幅（uplift）として残ります。
  • データセット依存性: Planck-only データでは転移が早い ($z^\dagger \approx 1.9$) ためピークが早期に現れ、Full データでは転移が遅い ($z^\dagger \approx 1.7$) ためピークがやや遅れて強く現れます。

B. 観測可能なスケールとの対応

この増幅が生じる波数範囲は、観測的に非常に重要なスケールに対応しています。

• 物理的スケール: $k \simeq 0.6-1.0 \, h \text{ Mpc}^{-1}$ は、銀河団（poor-cluster）や銀河群（group）のスケール（質量 $M \sim 10^{13}-10^{14} M_\odot$）に対応します。
• 観測手法: このスケールは、弱い重力レンズ（weak lensing）、銀河 - 銀河レンズ、銀河団の数え上げ、および熱的サニャエフ・ゼルドビッチ（tSZ）効果の感度が最も高い領域です。

C. 従来のパラメータ調整との区別

この特徴は、単なる $\sigma_8$ や $S_8$ のスケール非依存な再スケーリング（振幅の全体的な変更）では説明できません。

• 動的な特徴: 増幅のピーク位置（$k$）が赤方偏移とともに移動し、転移のタイミング（$z^\dagger$）と強く相関しています。これは、モデルの動的な歴史（AdS 相の持続時間と強度、dS 相での減衰）を反映した「形状特徴（shape feature）」です。

D. 「宇宙の正午（Cosmic Noon）」との関連

観測的に恒星形成率密度（SFRD）がピークに達する「宇宙の正午（$z \sim 1.5-2$）」と、$\Lambda_s$CDM モデルでパワースペクトルの増幅が最大となる時期が一致します。

• このモデルは、恒星形成が活発な時期に、重力背景が強化されることを示唆しており、銀河の形成やハローの成長効率に直接的な影響を与える可能性があります。

4. 意義 (Significance)

1. $\Lambda_s$CDM モデルの検証可能性の確立:
   本研究は、$\Lambda_s$CDM モデルが単なるパラメータの調整ではなく、物質パワースペクトルに**「赤方偏移に依存する局所的な増幅の山」**という明確で検証可能なシグナルを予言することを示しました。これは、将来の高精度観測（Euclid, LSST, CMB-S4 など）で直接検出、あるいは排除可能なターゲットとなります。

2. 宇宙論的緊張の解決への示唆:
   このモデルは、$H_0$ 緊張と $S_8$ 緊張を同時に緩和するメカニズムを提供します。AdS 相での成長促進と、その後の抑制が、観測される構造形成の振幅とハッブル定数の値を自然に調整する可能性があります。

3. 非線形領域における一般相対論的効果の理解:
   相対論的 N 体シミュレーションを用いることで、背景の宇宙論的進化（宇宙定数の符号反転）が非線形構造形成に与える影響を、ニュートン近似を超えて正確に捉えることができました。

4. 今後の研究方向:
   本研究は「急激な転移（瞬間的近似）」に基づいていますが、滑らかな転移（スカラー場モデルなど）や、修正重力理論への拡張、バリオンフィードバックの効果を組み込んだ研究が次のステップとして期待されます。また、この「山」のシグナルを、tSZ や弱い重力レンズのデータと組み合わせることで、宇宙の加速膨張のメカニズム（宇宙定数の符号反転の有無）を決定づけることが可能になるでしょう。

要約すると、この論文は、$\Lambda_s$CDM モデルが標準モデルと区別される決定的な「指紋」を非線形物質パワースペクトルに刻み込むことをシミュレーションで証明し、次世代の宇宙論観測による新物理の発見への道筋を示した重要な研究です。