[論文レビュー] High-accuracy emulators for observables in $Λ$CDM, $N_\mathrm{eff}$, $Σm_ν$, and $w$ cosmologies

Boris Bolliet, A. Spurio Mancini+5 more|arXiv (Cornell University)|Mar 2, 2023

Computational Physics and Python Applications被引用数 4

ひとこと要約

本論文では、CosmoPowerフレームワークを用いて、宇宙論的観測量—CMBパワースペクトル、物質パワースペクトル、赤方偏移-距離関係、大規模構造プローブ—の高精度ニューラルネットワークエミュレータを提示する。ΛCDM、$w$CDM、$N_{\mathrm{eff}}$、$\Sigma m_\nu$モデルにおける高精度のアインシュタイン＝ボルツマン計算のトレーニングを経て、CMBでは$\ell=10^4$まで、物質パワースペクトルでは$k=50\,\mathrm{Mpc}^{-1}$まで、0.5%未満の精度を達成し、MCMCを用いたラップトップ上でのステージIV CMBおよびステージIII LSSデータ解析を可能にする。

ABSTRACT

We use the emulation framework CosmoPower to construct and publicly release neural network emulators of cosmological observables, including the Cosmic Microwave Background (CMB) temperature and polarization power spectra, matter power spectrum, distance-redshift relation, baryon acoustic oscillation (BAO) and redshift-space distortion (RSD) observables, and derived parameters. We train our emulators on Einstein-Boltzmann calculations obtained with high-precision numerical convergence settings, for a wide range of cosmological models including $Λ$CDM, $w$CDM, $Λ$CDM+$N_\mathrm{eff}$, and $Λ$CDM+$Σm_ν$. Our CMB emulators are accurate to better than 0.5% out to $\ell=10^4$ which is sufficient for Stage-IV data analysis, and our $P(k)$ emulators reach the same accuracy level out to $k=50 \,\, \mathrm{Mpc}^{-1}$, which is sufficient for Stage-III data analysis. We release the emulators via an online repository CosmoPower Organisation, which will be continually updated with additional extended cosmological models. Our emulators accelerate cosmological data analysis by orders of magnitude, enabling cosmological parameter extraction analyses, using current survey data, to be performed on a laptop. We validate our emulators by comparing them to CLASS and CAMB and by reproducing cosmological parameter constraints derived from Planck TT, TE, EE, and CMB lensing data, as well as from the Atacama Cosmology Telescope Data Release 4 CMB data, Dark Energy Survey Year-1 galaxy lensing and clustering data, and Baryon Oscillation Spectroscopic Survey Data Release 12 BAO and RSD data.

研究の動機と目的

• 遅いボルツマンコードの呼び出しを高速なニューラルネットワークエミュレータに置き換えることで、宇宙論的パラメータ推論を加速すること。
• ΛCDM、$w$CDM、ΛCDM+$N_{\mathrm{eff}}$、ΛCDM+$\Sigma m_\nu$を含む多様なモデルにおいて、主要な宇宙論的観測量の高精度（0.5%未満）のエミュレーションを達成すること。
• 標準的なハードウェア（例：ラップトップ）上で、現在および将来のCMBおよび大規模構造データの高精度なMCMC解析を可能にすること。
• 宇宙論コミュニティが利用可能な公開リポジトリとして、継続的に更新可能なエミュレータのリポジトリを提供すること。
• エミュレータの妥当性を、既存のコード（class、camb）およびPlanck、ACT DR4、DES Y1、BOSS DR12の実データ尤度関数と比較して検証すること。

提案手法

• エミュレータは、宇宙論的観測量を対象としたTensorFlowベースのニューラルネットワークパイプラインであるCosmoPowerフレームワークを用いて構築された。
• トレーニングデータは、厳密な数値収束設定を用いた高精度なアインシュタイン＝ボルツマンソルバー（cambおよびclass）の出力から生成された。
• エミュレータは、$\Omega_{\mathrm{b}}$、$\Omega_{\mathrm{c}}$、$n_{\mathrm{s}}$、$A_{\mathrm{s}}$、$H_0$、$N_{\mathrm{eff}}$、$\Sigma m_\nu$、$w$といった多様な宇宙論的パラメータを含む、複数の宇宙論的モデルにおいてトレーニングされた。
• フレームワークは、CMB温度（TT）、偏光（TE、EE）、レンズ効果（PP）、物質パワースペクトル（PKL、PKNL）、距離-赤方偏移関係（DA、H）、$\sigma_8(z)$、および派生パラメータ（DER）のエミュレーションをサポートする。
• 検証は、エミュレータ出力とclassおよびcambの結果との比較、および実データ尤度関数（Planck、ACT DR4、DES Y1、BOSS DR12）からの事後分布の再現性によって実施された。
• エミュレータはオンラインリポジトリ（CosmoPower Organisation）を通じて公開されており、cobaya や cosmosis といった標準的なMCMCパイプラインへの統合を想定して設計されている。

実験結果

リサーチクエスチョン

• RQ1 ニューラルネットワークエミュレータは、$\ell=10^4$および$k=50\,\mathrm{Mpc}^{-1}$までCMBおよび物質パワースペクトルの予測において0.5%未満の精度を達成できるか？
• RQ2 これらのエミュレータは、実データ（例：Planck、ACT DR4、DES Y1、BOSS DR12）からの宇宙論的パラメータ制約を高い忠実度で再現できるか？
• RQ3 ステージIIIおよびステージIVの宇宙論的調査において、エミュレータベースのMCMC推論は、従来のボルツマンコード呼び出しに比べて、速度と精度の面でどの程度置き換え可能か？
• RQ4 $w$CDM や $\Sigma m_\nu$ などの拡張モデルを含む多様な宇宙論的モデルにおいて、エミュレータの頑健性はどの程度か？
• RQ5 エミュレータは、cobaya や cosmosis といった標準的な宇宙論的解析フレームワークに効果的に統合できるか？

主な発見

• CMBエミュレータは、$\ell=10^4$まで0.5%未満の精度を達成し、ステージIV CMBデータ解析に十分な水準に達している。
• 物質パワースペクトルエミュレータは、$k=50\,\mathrm{Mpc}^{-1}$まで0.5%未満の精度を維持しており、ステージIIIの大規模構造調査の要件を満たしている。
• エミュレータは、PlanckのTT、TE、EE、CMBレンズ効果尤度関数に加え、ACT DR4、DES Y1、BOSS DR12のデータから得られる周辺化事後分布を正確に再現した。
• エミュレータを用いたMCMCパラメータ推論は、標準的なラップトップでも実行可能であり、宇宙論的解析の速度を桁違いに向上させた。
• エミュレータはCosmoPower Organisationリポジトリを通じて公開されており、今後の新しい宇宙論的モデルへの継続的アップデートを想定して設計されている。
• フレームワークは拡張可能であり、将来のBモードパワースペクトルやその他の高度な観測量のエミュレータ対応も想定している。