Non-linear evolution in $f(R)$ gravity: perturbative modelling of the Chameleon mechanism

以下は、提示された論文「Non-linear evolution in f(R) gravity: perturbative modelling of the Chameleon mechanism（f(R) 重力における非線形進化：カメレオン機構の摂動論的モデリング）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

標準的な$\Lambda$CDM 宇宙モデルは、ダークエネルギーの起源やハッブル定数・$\sigma_8$の緊張関係（tensions）など、理論的・観測的な課題に直面しています。これに対抗する有力な候補として、一般相対性理論（GR）を修正する$f(R)$重力理論が挙げられます。
$f(R)$重力では、アインシュタイン・ヒルベルト作用におけるリッチスカラー$R$を一般関数$f(R)$に置き換えることで、追加の自由度（スカラー場）を導入します。しかし、太陽系内の高精度実験と整合させるためには、高密度領域では GR が回復し、低密度の宇宙論的スケールでは修正が現れる「スクリーニング機構」が必要です。この$f(R)$モデルにおける代表的な機構がカメレオン機構です。

課題:

• カメレオン機構を考慮すると、スカラー場$\chi$（$\chi = \Phi - \Psi$、$\Phi$と$\Psi$は計量摂動のポテンシャル）の方程式が非線形になります。
• 従来の数値シミュレーション（N-body シミュレーション）は計算コストが高く、特に非線形領域でのスカラー場を解くために FFT 以外の sophisticated なソルバー（マルチグリッド法など）が必要となり、計算が重くなります。
• 解析的な摂動論的アプローチは、非線形項を無視するか、あるいは複雑な境界値問題を数値的に厳密に解く必要があり、効率的な解析解の構築が困難でした。

2. 手法 (Methodology)

本研究では、球対称な摂動（滑らかな補償トッパット・プロファイル）における非線形進化を解析的に近似しつつ、カメレオン機構を摂動論的に取り込む新しい手法を提案しました。

• モデル: Hu-Sawicki モデル（$f(R)$の具体的な形式）を採用し、背景宇宙の進化は$\Lambda$CDM と一致するようにパラメータを設定しました。
• 変数の定義: 重力ポテンシャルの和$\Phi_+ = (\Phi+\Psi)/2$と差$\chi = \Phi-\Psi$を定義します。$\Phi_+$はポアソン方程式を満たしますが、$\chi$は非線形微分方程式を満たします。
• 摂動論的解法:
  • $\chi$の方程式を、カメレオン機構を無視した場合の解（1 次近似）を基準として、密度揺らぎ$\delta$やリッチスカラーの摂動$\delta R$の関数として展開します。
  • $\chi = \sum \epsilon^n \chi^{(n)}$と仮定し、$\chi^{(1)}$（線形解）を求めた後、その解を用いて$\delta R$を計算し、それをソース項として高次項（2 次〜4 次）を逐次的に解くアルゴリズムを構築しました。
  • 各反復ステップで、線形微分演算子の逆演算（解析解または FFT）を使用するため、計算コストが低く抑えられます。
• 進化シミュレーション:
  • 混合オイラー・ラグランジュ法を採用。密度と速度の時間進化はラグランジュ座標（殻の位置追跡）で、ポテンシャルの空間分布はオイラー座標で計算します。
  • 各時間ステップで、密度分布からリッチスカラー$R$を再計算し、それに基づいてコンプトン波長（$\bar{x}_C$）を更新することで、カメレオン機構の空間・時間依存性を追跡します。

3. 主要な貢献 (Key Contributions)

• 非線形方程式の効率的な摂動解: カメレオン機構を含む$f(R)$重力の非線形スカラー場方程式を、数値ソルバーに依存せず、解析的な反復法で近似解く手法を確立しました。
• コンプトンスケールの時空間進化の可視化: 密度揺らぎの進化に伴うコンプトンスケール（$\bar{x}_C$）の空間的・時間的変化を詳細に追跡し、それが密度分布にどのようにフィードバックするかを明らかにしました。
• 密度 - 速度発散関係（DVDR）の非一対一性の発見: カメレオン機構が働く領域では、同じ密度値に対して複数の速度発散値が存在しうる（非一対一性）ことを初めて示しました。

4. 結果 (Results)

研究では、コンプトンスケールと摂動の長さスケールの比$Q = \bar{x}_C / x_{top}$に基づいて、3 つの進化領域を定義し、それぞれで異なる振る舞いを示しました。

1. 強場領域 ($Q \gg 1$):
   • 摂動スケールがコンプトンスケールより十分小さい場合。
   • カメレオン機構は機能せず、重力定数が$G_{eff} = \frac{4}{3}G$に増幅された GR と同じように進化します。密度成長は GR よりも速くなります。
2. 弱場領域 ($Q \ll 1$):
   • 摂動スケールがコンプトンスケールより十分大きい場合。
   • 修正重力の効果は抑制され、進化は標準的な GR とほぼ区別がつかなくなります。
3. 中間領域 ($Q \sim 1$):
   • 摂動スケールがコンプトンスケールと同程度の場合。
   • トッパットの縁での密度増大: 既知の現象として、トッパットの縁（密度勾配が急峻な部分）で密度が増大し、トッパット形状が崩れることが確認されました。
   • 原点付近の密度増大（新規発見）: 本研究で初めて報告された現象として、トッパットの中心（原点）付近でもわずかな密度増大が観測されました。
     • メカニズム: 外殻が内側に落下する過程で、中心付近の密度が一旦低下し、それによりコンプトンスケールが局所的に増大します。増大したコンプトンスケールは中心付近での引力を強め、さらに物質が中心に集まるという正のフィードバックループが生じ、中心付近に密度のピーク（バンプ）が形成されます。

密度 - 速度発散関係（DVDR）:

• 標準 GR やカメレオン機構を無視した$f(R)$では、非線形領域でも密度と速度発散の関係は単一の曲線（Zeldovich 曲線）で記述されます。
• しかし、カメレオン機構が有効な中間領域では、同じ密度値に対して、位置（半径）によって異なる速度発散値を持つようになり、DVDR が多価関数となります。これは、密度がコンプトンスケールを決定し、それが再び密度成長に影響を与えるという非線形フィードバック（ヒステリシス的な挙動）に起因します。

5. 意義 (Significance)

• 計算効率の向上: 従来の N-body シミュレーションや境界値問題ソルバーに依存せず、解析的な摂動展開を用いることで、f(R) 重力の非線形進化を低コストでシミュレートできる手法を提供しました。
• 物理的洞察の深化: カメレオン機構が宇宙大規模構造の形成に与える影響を、単なる「成長率の変化」だけでなく、「密度プロファイルの形状変化（縁と中心の増大）」や「DVDR の多価性」といった新しい特徴として捉えました。
• 将来の観測への示唆: 将来の宇宙大規模構造サーベイ（Euclid, DESI, Rubin 観測所など）において、密度分布の詳細な形状や、密度 - 速度関係の非単調性を観測することで、標準重力と修正重力を区別する新たな手がかりが得られる可能性があります。特に、DVDR のヒステリシス的な挙動は、修正重力の強力なシグネチャとなり得ます。

この研究は、球対称という特殊なケースに基づいていますが、得られた洞察は 3 次元ランダム初期条件を持つ N-body シミュレーションの解析や、修正重力理論の観測的検証において重要な指針となると期待されます。