CAMB v2: cosmological power spectra for high-precision surveys

技術要約: CAMB 2: 高精度サーベイのための宇宙論的パワースペクトル

問題提起
Simons Observatoryやその発展形であるAdvanced SOといった、次世代の宇宙マイクロ波背景放射（CMB）および大規模構造サーベイは、観測の不確かさを大幅に下回る数値誤差を持つ理論的パワースペクトルを要求している。従来のCAMB（v1.x）やCLASSは、Planck時代の解析（$\ell \lesssim 2000$において$\sim 0.1\%$の精度が必要）には十分であったが、より広いスケール範囲および非平坦な宇宙論において、より厳格なパーミル（$10^{-3}$）精度の要件を満たすには課題があった。具体的には、曲がった時空における超球面ベッセル関数の評価、再結合履歴の取り扱い、質量を持つニュートリノのサンプリング、およびダークエネルギーの摂動に関する更新が、精度を犠牲にすることなく堅牢性と速度を確保するために不可欠であった。

手法および主要な開発事項
本論文は、レンズ効果を受けたCMBおよび物質パワースペクトルに焦点を当てた、CAMB（v2.0.0）の包括的なアップデートについて記述している。開発哲学は、汎用的な特殊関数ライブラリに依存して精度と速度のトレードオフを固定する手法から、物理的問題とサンプリングされるグリッドに明示的に較正されたカスタム精度の数値的手法へとシフトしている。この開発の大部分は、アルゴリズムの提案、デバッグ、および数値的変更の検証を行うために、人間の監督下にある大規模言語モデル（LLM）およびAIエージェントを用いて行われた。

主な技術的貢献は以下の通りである：

1. 超球面ベッセル関数（線形視積分）:

   • 革新性: 非平坦な宇宙論における線形視積分に使用される超球面ベッセル関数（$\phi^\nu_l$）の新しい取り扱い。以前のWKB近似のように曲がった動径方程式をエアリー関数へ写像するのではなく、これを平坦な球面ベッセル方程式へと写像する。
   • メカニズム: **オルバー・アクション写像（Olver action map）**を用い、曲がった方程式の作用積分を、転回点（turning point）を通じて平坦な方程式の作用積分と一致させる。これにより、座標変換 $z(\chi)$ と振幅因子 $A(\chi)$ が定義される。
   • 特性: この近似は平坦な極限において正確であり、転回点を通っても滑らかであり、近平坦な領域では平坦なベッセル引数と振幅の単純な再スケーリングへと還元される。
   • 実装: この写像は、ソースポイント間での効率的なヌメロフ（Numerov）伝播のための正確な初期値を提供し、各ステップでの高価な逐次評価を回避する。低マルチポールや無効なコーナーに対しては、フォールバックメカニズム（厳密な漸化式、エアリー近似）が使用される。

2. 再結合および背景進化:

   • 再結合モデル: デフォルトのRECFAST有効数レベル近似を、高精度なCosmoRecのリファレンスに対して再較正した。ヘリウムIのレート方程式に新しい現象論的な補正因子を追加し、特にヘリウム再結合の窓（$1500 \lesssim z \lesssim 3000$）におけるコヒーレントな残留物を除去した。
   • BBNの一貫性: コードは、古いPArthENoPEベースのテーブルに代わり、$\Omega_b h^2$ および $\Delta N_{eff}$ に基づいてPRIMAT (2024) から原始ヘリウム存在量（$Y_p$）を補間するようになった。
   • 積分: 再結合の硬い（stiff）初期フェーズは、アダプティブなDormand–Prince積分器に引き渡される前に、暗黙的Rosenbrock法（ROS2）によって処理される。

3. 質量を持つニュートリノ:

   • 背景密度と圧力は、事前計算されたテーブルではなく、直接的な低次フィットを用いて計算される。
   • 運動量サンプリング: ニュートリノの質量に基づき、Fermi–Diracカーネルの正確な積分を保証するために、サンプリングされる運動量ポイントの数（3から4または5へ）を自動的に調整する。
   • 階層の切り替え: 高エネルギーの相対論的種については摂動論的な相対論的展開を用い、後期進化については非相対論的な速度積分階層を用い、質量と赤方偏移に基づいてこれらを切り替える。

4. ダークエネルギー (PPF):

   • ファントム・ディバイド（$w = -1$）を横切るダークエネルギー流体に対して、**安定化されたパラメータ化されたポスト・フリードマン（PPF）**進化を実装した。
   • 正則化: 高い $k/H$（緩和率が $k^2$ に比例して増大する領域）における数値的な硬さを回避するため、緩和率を閾値 $\beta_\star = 30$ でキャップしている。これにより、連続性を維持しつつハードカットオフによるアーティファクトを避け、数値的な硬さと積分許容誤差への感度を取り除く。

5. レンズ効果と精度:

   • レンズ計算には、切断によるリンギングを抑制するための新しい $C^2$ 連続（五次）テーパーを用いた、曲がった空の相関関数法を使用している。
   • コードは、要求された $\ell_{max}$ まで収束を保証するために必要な最大波数（$k_{max}$）を自動的に決定する。

結果と検証
著者らは、曲率、ニュートリノ質量、ダークエネルギー、および再結合モデルの変動を含む100以上の宇宙論的構成のグリッドを用いて、精度パラメータを倍増させた「ブーストされた」リファレンス計算に対して、新しいデフォルト設定を検証した。

• 収束性: アンブーストのデフォルト設定は、主要なレンズ効果を受けたCMBの範囲（$600 \le \ell \le 3500$）および準線形物質パワースペクトルにおいて、保守的な$10^{-3}$の点別収束目標を満たしている。
• CLASSとの比較: 一致した物理的仮定と高精度設定で実行した場合、CAMB v2とCLASS v3.3.4は、平坦な$\Lambda$CDMにおいて$\sim 10^{-4}$のレベルで一致する。
• パフォーマンス:
  • 平坦なモデル: 新しいデフォルトは、高$k$サポートに対して一致させた場合、以前のバージョンと同等か、わずかに高速であるが、より高い精度を提供している（高$\ell$において誤差を$\sim 1\%$から$< 0.1\%$に低減）。
  • 非平坦なモデル: 新しい超球面ベッセルスキームは、一致した精度レベルにおいて、非平坦モデルの計算コストをCAMB v1.6.6と比較して17分の1に削減した。近平坦なモデル（$|\Omega_K| \le 10^{-3}$）では、平坦なケースに対するオーバーヘッドはわずかである。
• 尤度への影響: グリッドレベルの一致により、中央値の $\Delta\chi^2 \approx 0.4$ がブーストされたリファレンスに対して得られ、これは許容範囲内であり、数値誤差が次世代サーベイのパラメータ推論において無視できることを示している。

意義と主張
本論文は、CAMB v2.0.0が、ユーザーが精度パラメータを手動で調整することなく、次世代の精密なCMBおよび大規模構造サーベイに必要な速度と数値精度を提供することを主張している。中心的な貢献は、超球面関数のための平坦なベッセル比較法であり、これが曲がった空間の積分における技術的ボトルネックを解決した。

さらに、本論文は方法論的な転換を強調している。すなわち、汎用ライブラリに頼るのではなく、物理的問題のニーズに合わせて特別に較正されたカスタム精度の数値ルーチンの開発である。著者らは、複雑な数値アルゴリズムを開発、テスト、および検証するために、人間の監督下にあるAIエージェントを使用することの証明としてこの研究を提示しており、このアプローチが、多数の小さな数値的変更を適合させ検証するという労力の大きい作業を実用的なものにできることを示している。本論文は物理的な謎を解明すると主張するものではなく、次世代の宇宙論的データ解析のための堅牢で高忠実度な計算ツールを提供することを目的としている。