[論文レビュー] Euclid preparation - LXXI. Simulations and nonlinearities beyond ΛCDM. 3. Constraints on f(R) models from the photometric primary probes
本稿は、Euclidミッションの光度測定プローブを用いてf(R)重力の制約を評価し、非線形物質パワースペクトルモデル化における不確実性がパラメータ回復に与える影響を検討する。f(R)パラメータ制約におけるバイアスは、バリオン物理学とのデゲネラシーに起因するが、理論的誤差またはバリオンパラメータの事前分布を組み込むことで軽減可能であり、悲観的設定では事前分布に依存しない上限としてlog₁₀|fR₀| < −5.6が得られる。
We study the constraint on $f(R)$ gravity that can be obtained by photometric primary probes of the Euclid mission. Our focus is the dependence of the constraint on the theoretical modelling of the nonlinear matter power spectrum. In the Hu-Sawicki $f(R)$ gravity model, we consider four different predictions for the ratio between the power spectrum in $f(R)$ and that in $Λ$CDM: a fitting formula, the halo model reaction approach, ReACT and two emulators based on dark matter only $N$-body simulations, FORGE and e-Mantis. These predictions are added to the MontePython implementation to predict the angular power spectra for weak lensing (WL), photometric galaxy clustering and their cross-correlation. By running Markov Chain Monte Carlo, we compare constraints on parameters and investigate the bias of the recovered $f(R)$ parameter if the data are created by a different model. For the pessimistic setting of WL, one dimensional bias for the $f(R)$ parameter, $\log_{10}|f_{R0}|$, is found to be $0.5 σ$ when FORGE is used to create the synthetic data with $\log_{10}|f_{R0}| =-5.301$ and fitted by e-Mantis. The impact of baryonic physics on WL is studied by using a baryonification emulator BCemu. For the optimistic setting, the $f(R)$ parameter and two main baryon parameters are well constrained despite the degeneracies among these parameters. However, the difference in the nonlinear dark matter prediction can be compensated by the adjustment of baryon parameters, and the one-dimensional marginalised constraint on $\log_{10}|f_{R0}|$ is biased. This bias can be avoided in the pessimistic setting at the expense of weaker constraints. For the pessimistic setting, using the $Λ$CDM synthetic data for WL, we obtain the prior-independent upper limit of $\log_{10}|f_{R0}|< -5.6$. Finally, we implement a method to include theoretical errors to avoid the bias.
研究の動機と目的
- Euclidの光度測定主プローブを用いて、非線形物質パワー スペクトルモデル化がf(R)重力の制約に与える影響を評価すること。
- FORGE、e-Mantis、ReACTなどのN体シミュレーションに基づく異なる予測の間の差異が、f(R)パラメータの回復に与える影響を調査すること。
- バリオン物理学が、バリオンパラメータとのデゲネラシーを通じてf(R)パラメータ制約をバイアス化する役割を検討すること。
- 不確実な非線形パワー スペクトル予測に起因するバイアスを低減するため、尤度関数に理論的誤差を組み込む手法を開発・検証すること。
- 理論的モデル化の不確実性に強く耐性を持つ、事前分布に依存しないf(R)重力の堅牢な制約を確立すること。
提案手法
- 非線形物質パワー スペクトル比の4種類の異なる予測を用いる：フィッティング式、ハロモデル反応アプローチ（ReACT）、およびダークマターのみのN体シミュレーションに基づく2つのエミュレータ（FORGEとe-Mantis）。
- これらの予測をMontePython宇宙論的パラメータ推定フレームワークに統合し、弱引力レンズ（WL）、光度測定銀河クラスタリング、およびそれらの相互相関の角度スペクトルを計算する。
- 異なる理論的モデル化仮定下でのf(R)パラメータ制約を比較するために、マルコフ連鎖モンテカルロ（MCMC）計算を実行する。
- 合成データを1つのモデル（例：FORGE）で生成し、別のモデル（例：e-Mantis）でフィッティングすることで、回復されたf(R)パラメータのバイアスを評価する。
- バリオンフィケーションエミュレータ（BCemu）を組み込み、バリオン物理学がWL制約に与える影響とデゲネラシーを検討する。
- FORGEとReACTの保守的な差異を用いて理論的誤差推定を実施し、非線形パワー スペクトル予測の不確実性を考慮する。
実験結果
リサーチクエスチョン
- RQ1異なる非線形物質パワー スペクトルの理論的モデルが、Euclidの光度測定データにおけるf(R)重力パラメータの制約にどのように影響を与えるか？
- RQ2f(R)パラメータとバリオン物理学パラメータの間のデゲネラシーが、log₁₀|fR₀|の1次元マージナル化制約にどの程度バイアスを引き起こすか？
- RQ3理論的誤差推定は、非線形パワー スペクトルモデル化の不正確さに起因するf(R)パラメータ回復のバイアスを軽減できるか？
- RQ4弱引力レンズデータのみを用いた悲観的観測設定下で、log₁₀|fR₀|の事前分布に依存しない上限は何か？
- RQ5異なる非線形構造形成モデルを用いる場合、バリオン物理学の組み込みがf(R)重力制約の堅牢性にどのように影響するか？
主な発見
- 弱引力レンズデータのみを用いた悲観的設定では、合成データをFORGEで生成しe-Mantisでフィッティングした場合、log₁₀|fR₀|パラメータの1次元バイアスが0.5σに達する。
- バリオン物理学を組み込むと、非線形ダークマター予測の差異がバリオンパラメータの調整で相殺されるため、log₁₀|fR₀|のマージナル化制約におけるバイアスが悪化する。
- デゲネラシーが存在するにもかかわらず、楽観的設定ではf(R)パラメータと2つの主要なバリオンパラメータが制約可能であるが、log₁₀|fR₀|の1次元マージナル化制約は依然としてバイアスを示す。
- 弱引力レンズ用のΛCDM合成データを用いた悲観的設定では、事前分布に依存しない上限としてlog₁₀|fR₀| < −5.6が得られる。
- FORGEとReACTの差異を用いた理論的誤差推定はバイアスを低減し、堅牢性を向上させる。これは、信頼できる制約を得るためにはこのような誤差モデル化が不可欠であることを示唆する。
- 本研究は、エミュレータの精度が極めて重要であることを確認した。FORGEはe-Mantisに比べてより大きなエミュレーション誤差を示しており、これがモデル依存のバイアスを引き起こしており、より優れたエミュレーション技術によりこれを軽減できる。
より良い研究を、今すぐ始めましょう
論文設計から論文執筆まで、研究時間を劇的に削減しましょう。
クレジットカード登録不要
このレビューはAIが作成し、人間の編集者が確認しました。