[論文レビュー] Euclid: Validation of the MontePython forecasting tools

Santiago Casas, Julien Lesgourgues+127 more|arXiv (Cornell University)|Jan 1, 2023

Insurance, Mortality, Demography, Risk Management参考文献 2被引用数 2

ひとこと要約

この論文は、モンテカルロ・プロセスを用いたフィッシャー行列予測の妥当性を、Euclid宇宙ミッションのためのモンテピーロの予測パイプラインに対して検証する。本研究では、模擬尤度関数の数値的2階微分を用いて得られたフィッシャー行列予測と、既存のCosmicFish予測およびEuclidコラボレーションのEP:VII結果とを厳密に比較した。複数の設定において、2%以内の極めて良好な一致が確認され、HP精度設定が安定性を確保するために不可欠であることが示された。さらに、フィッシャー楕円とMCMCによる信頼領域が一致することから、ガウス近似が有効であることが裏付けられ、これによりΛCDMを超える拡張的宇宙論モデルに対する信頼性の高い予測が可能になった。

ABSTRACT

The Euclid mission of the European Space Agency will perform a survey of weak lensing cosmic shear and galaxy clustering in order to constrain cosmological models and fundamental physics. We expand and adjust the mock Euclid likelihoods of the MontePython software in order to match the exact recipes used in previous Euclid Fisher matrix forecasts for several probes: weak lensing cosmic shear, photometric galaxy clustering, the cross-correlation between the latter observables, and spectroscopic galaxy clustering. We also establish which precision settings are required when running the Einstein-Boltzmann solvers CLASS and CAMB in the context of Euclid. For the minimal cosmological model, extended to include dynamical dark energy, we perform Fisher matrix forecasts based directly on a numerical evaluation of second derivatives of the likelihood with respect to model parameters. We compare our results with those of other forecasting methods and tools. We show that such MontePython forecasts agree very well with previous Fisher forecasts published by the Euclid Collaboration, and also, with new forecasts produced by the CosmicFish code, now interfaced directly with the two Einstein-Boltzmann solvers CAMB and CLASS. Moreover, to establish the validity of the Gaussian approximation, we show that the Fisher matrix marginal error contours coincide with the credible regions obtained when running Monte Carlo Markov Chains with MontePython while using the exact same mock likelihoods. The new Euclid forecast pipelines presented here are ready for use with additional cosmological parameters, in order to explore extended cosmological models.

研究の動機と目的

• EuclidコラボレーションおよびCosmicFishから得られた既存のフィッシャー行列予測と比較し、MontePythonのEuclid予測パイプラインの妥当性を検証すること。
• Einstein-BoltzmannソルバーCAMBおよびCLASSの数値的安定性と、Euclid予測の文脈における必要精度設定を評価すること。
• 同一の模擬尤度関数を用いてMCMCで得られた信頼領域と比較することで、フィッシャー行列予測におけるガウス近似の妥当性を検証すること。
• MontePythonにおけるフィッシャーとMCMCモードの整合性を保証することで、非標準的宇宙論モデルに対する堅牢で拡張可能な予測を可能にすること。
• CAMBおよびCLASSの精度設定を文書化・標準化し、予測パイプラインにおける相互の整合性と数値的安定性を確保すること。

提案手法

• モンテピーロにおいて、有限差分法を用いて模擬Euclid尤度関数の最適適合点まわりの2階微分を数値的に計算し、フィッシャー行列を算出する。
• 同じ理論的予測を用いるが異なるアルゴリズム実装を採用するCosmicFishの予測と、モンテピーロの予測を直接比較する。
• CAMBおよびCLASSによるオンザフライおよび事前計算済みのEinstein-Boltzmannソルバー（EBS）呼び出しを併用し、一貫性と数値的ロバストネスを検証する。
• 同一の模擬尤度関数を用いてモンテピーロでMCMCチェインを走らせ、得られた信頼領域とフィッシャー楕円を比較することで、ガウス近似の妥当性を検証する。
• 異なるEBS精度設定（例：DP対HP）およびステップサイズを体系的にテストし、数値的ノイズと高次非ガウス性の影響を最小限に抑える最適な設定を同定する。
• CAMBおよびCLASSとの4通りのインターフェース（事前計算済み理論ファイル有無を含む）を実装・ベンチマークし、パイプライン間の一貫性を確保する。

実験結果

リサーチクエスチョン

• RQ1 ΛCDM+{w0, wa}モデルに対して、モンテピーロのフィッシャー行列予測は、CosmicFishおよびEuclidコラボレーションのEP:VII予測とどの程度一致するか？
• RQ2 CAMBおよびCLASSで、Euclid予測のための安定的かつ正確なフィッシャー行列予測を実現するために必要な精度設定（例：DP対HP）は何か？
• RQ3 フィッシャー楕円とMCMCで得られた信頼領域の重複度から、Euclid予測におけるガウス近似の妥当性はどの程度保証されるか？
• RQ4 同一の模擬尤度関数を用いた場合、モンテピーロパイプラインはフィッシャーモードとMCMCモードの両方で一貫した予測を信頼できるか？
• RQ5 ステップサイズおよび数値的ノイズが、Euclid予測の文脈における有限差分法に基づくフィッシャー行列計算の精度に及ぼす影響は何か？

主な発見

• モンテピーロのフィッシャー行列予測は、写真的および分光的プローブの両方において、悲観的設定下でもCosmicFishおよびEP:VII予測と2%以内の一致を示した。
• モンテピーロパイプラインにおける安定的かつ正確なフィッシャー行列計算には、CAMBおよびCLASSのHP精度設定が不可欠であるが、CosmicFishではDP設定で十分である。
• モンテピーロのフィッシャー行列のマージナル誤差等高線は、同一の模擬尤度関数を用いてMCMCチェインから得られた信頼領域と極めて良好に一致しており、ΛCDM+{w0, wa}モデルにおけるガウス近似の妥当性が裏付けられた。
• CAMBおよびCLASSの両方をPythonラッパー経由で呼び出す場合、予測間の差異が0.5%未満に抑えられ、EBS間の高い一貫性が確認された。
• MP/Fisherパイプラインは、小さなステップサイズでは数値的ノイズに敏感であり、大きなステップサイズでは非ガウス性の影響を受けるが、写真的プローブでは10%、分光的プローブでは5%のマージナルエラーに対するステップサイズとHP設定を組み合わせることで安定した結果が得られた。
• 本検証により、モンテピーロの模擬尤度関数が堅牢かつ一貫していることが確認され、フィッシャーおよびMCMCモードの両方で、ΛCDMを超える拡張的宇宙論モデルに対する信頼性の高い予測が可能になった。