[논문 리뷰] Euclid: Validation of the MontePython forecasting tools

Santiago Casas, Julien Lesgourgues외 127명|arXiv (Cornell University)|2023. 01. 01.

Insurance, Mortality, Demography, Risk Management참고 문헌 2인용 수 2

한 줄 요약

이 논문은 몽테피리안 예측 파이프라인의 유클리드 우주선 임무에 대한 정밀한 검증을 위해, 몽테피리안에서 수치적 두 번째 도함수를 이용해 계산한 피셔 행렬 예측값을 가짜 우주론적 우도 함수 기반으로 산출한 결과와 기존의 코스믹피시 예측값 및 유클리드 협회(EP:VII)의 피셔 행렬 결과를 철저히 비교한다. 다양한 설정에서 거의 완벽한 일치(2% 이내)를 확인하였고, 고정밀도 설정(High Precision, 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.

연구 동기 및 목표

• 유클리드 협회 및 코스믹피시에서 제공하는 기존의 피셔 행렬 예측값과 대비하여 몽테피리안의 유클리드 예측 파이프라인을 검증하는 것.
• 유클리드 예측의 맥락에서 앙트로피-보츠만 해석기(CAMB) 및 CLASS의 수치적 안정성과 필요 정밀도 설정을 평가하는 것.
• 동일한 가짜 우도 함수를 사용하여 MCMC로 유도한 신뢰 영역과 비교함으로써 피셔 행렬 예측에서의 정규 근사의 타당성을 검증하는 것.
• 몽테피리안에서 피셔 모드와 MCMC 모드 간의 일관성을 확보함으로써 비표준 우주론 모델에 대한 견고하고 확장 가능한 예측을 가능하게 하는 것.
• CAMB 및 CLASS에 대한 정확도 설정을 문서화하고 표준화하여 상호 일치성과 수치적 안정성을 확보하는 것.

제안 방법

• 몽테피리안에서 유한 차분 기법을 사용해 가짜 유클리드 우도 함수의 최적 적합점 주변에서 두 번째 도함수를 수치적으로 계산하여 피셔 행렬을 산출한다.
• 동일한 이론 예측을 사용하지만 알고리즘 구현 방식이 다른 코스믹피시의 결과와 몽테피리안의 피셔 예측 결과를 직접 비교한다.
• CAMB 및 CLASS를 통한 실시간 및 사전 계산된 앙트로피-보츠만 해석기(EBS) 호출을 모두 사용하여 일관성과 수치적 탄력성을 시험한다.
• 동일한 가짜 우도 함수를 사용해 몽테피리안에서 MCMC 체인을 실행하고, 유도된 신뢰 영역과 피셔 타원체를 비교함으로써 정규 근사의 타당성을 검증한다.
• 다양한 EBS 정밀도 설정(DP 대비 HP)과 스텝 크기를 체계적으로 테스트하여 수치 노이즈와 고차 비정규성 영향을 최소화하는 최적 설정을 규명한다.
• CAMB 및 CLASS와 코스믹피시 간의 네 가지 인터페이스(사전 계산된 이론 파일 유무 포함)를 구현하고 벤치마킹하여 파이프라인 간 일관성을 확보한다.

실험 결과

연구 질문

• RQ1 ΛCDM+{w0, wa} 모델에 대해 몽테피리안의 피셔 행렬 예측값은 코스믹피시 및 유클리드 협회(EP:VII)의 예측값과 얼마나 일치하는가?
• RQ2 유클리드 예측에서 안정적이고 정확하며 상호 일치하는 피셔 행렬 예측을 확보하기 위해 CAMB 및 CLASS에서 필요한 정밀도 설정(DP 대비 HP)은 무엇인가?
• RQ3 피셔 타원체와 MCMC로 유도된 신뢰 영역 간의 겹침을 통해 유클리드 예측에서 정규 근사의 타당성은 어느 정도로 확인되는가?
• RQ4 동일한 가짜 우도 함수를 사용할 때 몽테피리안 파이프라인은 피셔 모드와 MCMC 모드에서 신뢰할 만한 일관성을 유지하는가?
• RQ5 스텝 크기와 수치적 노이즈가 유클리드 예측의 맥락에서 유한 차분 기반 피셔 행렬 계산의 정확도에 미치는 영향은 무엇인가?

주요 결과

• 피isher 행렬 예측값이 악조건 설정에서도 광학 및 분광 측정 기반 프로브에서 코스믹피시 및 EP:VII 예측값과 2% 이내로 일치함을 확인함.
• 몽테피리안 파이프라인에서 안정적이고 정확한 피셔 행렬 계산을 위해서는 CAMB 및 CLASS에서 HP 정밀도 설정이 필수적임을 규명하였고, 코스믹피시에서는 DP 설정으로도 충분함.
• 동일한 가짜 우도 함수를 사용해 MCMC 체인을 실행한 결과 유도된 신뢰 영역과 몽테피리안의 피셔 행렬 마진별 오차 윤곽선 간 뛰어난 일치를 보였으며, 이는 ΛCDM+{w0, wa} 모델에 대해 정규 근사가 타당함을 확인함.
• CAMB 및 CLASS를 모두 파이썬 래퍼를 통해 호출할 경우, 두 해석기 간 예측값 간 오차가 0.5% 미만으로 나타나, EBS 간 높은 일관성이 있음을 확인함.
• MP/Fisher 파이프라인은 작은 스텝 크기에서는 수치 노이즈에 매우 민감하고, 큰 스텝 크기에서는 고차 비정규성 영향을 받지만, 광학 프로브에 대해 10%, 분광 프로브에 대해 5%의 마진화된 오차를 기준으로 HP 설정과 함께 안정된 결과를 도출함.
• 검증 결과 몽테피리안의 가짜 우도 함수가 견고하고 일관됨을 확인하였으며, 이는 피셔 모드와 MCMC 모드에서 모두 ΛCDM를 초월한 확장된 우주론 모델에 대한 신뢰할 수 있는 예측을 가능하게 함.