skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Baryonic effects in cosmic shear tomography: PCA parametrization and importance of extreme baryonic models

Abstract

Baryonic effects are amongst the most severe systematics to the tomographic analysis of weak lensing data which is the principal probe in many future generations of cosmological surveys like LSST, Euclid etc.. Modeling or parameterizing these effects is essential in order to extract valuable constraints on cosmological parameters. In a recent paper, Eifler et al. (2015) suggested a reduction technique for baryonic effects by conducting a principal component analysis (PCA) and removing the largest baryonic eigenmodes from the data. In this article, we conducted the investigation further and addressed two critical aspects. Firstly, we performed the analysis by separating the simulations into training and test sets, computing a minimal set of principle components from the training set and examining the fits on the test set. We found that using only four parameters, corresponding to the four largest eigenmodes of the training set, the test sets can be fitted thoroughly with an RMS $$\sim 0.0011$$. Secondly, we explored the significance of outliers, the most exotic/extreme baryonic scenarios, in this method. We found that excluding the outliers from the training set results in a relatively bad fit and degraded the RMS by nearly a factor of 3. Therefore, for a direct employment of this method to the tomographic analysis of the weak lensing data, the principle components should be derived from a training set that comprises adequately exotic but reasonable models such that the reality is included inside the parameter domain sampled by the training set. The baryonic effects can be parameterized as the coefficients of these principle components and should be marginalized over the cosmological parameter space.

Authors:
ORCiD logo [1];  [1]
  1. Fermilab
Publication Date:
Research Org.:
Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States)
Sponsoring Org.:
USDOE Office of Science (SC), High Energy Physics (HEP) (SC-25)
OSTI Identifier:
1409075
Report Number(s):
FERMILAB-PUB-17-223-A; arXiv:1707.02332
1609189
DOE Contract Number:
AC02-07CH11359
Resource Type:
Journal Article
Resource Relation:
Journal Name: Astrophys.J.
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS

Citation Formats

Mohammed, Irshad, and Gnedin, Nickolay Y. Baryonic effects in cosmic shear tomography: PCA parametrization and importance of extreme baryonic models. United States: N. p., 2017. Web.
Mohammed, Irshad, & Gnedin, Nickolay Y. Baryonic effects in cosmic shear tomography: PCA parametrization and importance of extreme baryonic models. United States.
Mohammed, Irshad, and Gnedin, Nickolay Y. 2017. "Baryonic effects in cosmic shear tomography: PCA parametrization and importance of extreme baryonic models". United States. doi:. https://www.osti.gov/servlets/purl/1409075.
@article{osti_1409075,
title = {Baryonic effects in cosmic shear tomography: PCA parametrization and importance of extreme baryonic models},
author = {Mohammed, Irshad and Gnedin, Nickolay Y.},
abstractNote = {Baryonic effects are amongst the most severe systematics to the tomographic analysis of weak lensing data which is the principal probe in many future generations of cosmological surveys like LSST, Euclid etc.. Modeling or parameterizing these effects is essential in order to extract valuable constraints on cosmological parameters. In a recent paper, Eifler et al. (2015) suggested a reduction technique for baryonic effects by conducting a principal component analysis (PCA) and removing the largest baryonic eigenmodes from the data. In this article, we conducted the investigation further and addressed two critical aspects. Firstly, we performed the analysis by separating the simulations into training and test sets, computing a minimal set of principle components from the training set and examining the fits on the test set. We found that using only four parameters, corresponding to the four largest eigenmodes of the training set, the test sets can be fitted thoroughly with an RMS $\sim 0.0011$. Secondly, we explored the significance of outliers, the most exotic/extreme baryonic scenarios, in this method. We found that excluding the outliers from the training set results in a relatively bad fit and degraded the RMS by nearly a factor of 3. Therefore, for a direct employment of this method to the tomographic analysis of the weak lensing data, the principle components should be derived from a training set that comprises adequately exotic but reasonable models such that the reality is included inside the parameter domain sampled by the training set. The baryonic effects can be parameterized as the coefficients of these principle components and should be marginalized over the cosmological parameter space.},
doi = {},
journal = {Astrophys.J.},
number = ,
volume = ,
place = {United States},
year = 2017,
month = 7
}
  • Systematic uncertainties that have been subdominant in past large-scale structure (LSS) surveys are likely to exceed statistical uncertainties of current and future LSS data sets, potentially limiting the extraction of cosmological information. Here we present a general framework (PCA marginalization) to consistently incorporate systematic effects into a likelihood analysis. This technique naturally accounts for degeneracies between nuisance parameters and can substantially reduce the dimension of the parameter space that needs to be sampled. As a practical application, we apply PCA marginalization to account for baryonic physics as an uncertainty in cosmic shear tomography. Specifically, we use CosmoLike to run simulatedmore » likelihood analyses on three independent sets of numerical simulations, each covering a wide range of baryonic scenarios differing in cooling, star formation, and feedback mechanisms. We simulate a Stage III (Dark Energy Survey) and Stage IV (Large Synoptic Survey Telescope/Euclid) survey and find a substantial bias in cosmological constraints if baryonic physics is not accounted for. We then show that PCA marginalization (employing at most 3 to 4 nuisance parameters) removes this bias. Our study demonstrates that it is possible to obtain robust, precise constraints on the dark energy equation of state even in the presence of large levels of systematic uncertainty in astrophysical processes. We conclude that the PCA marginalization technique is a powerful, general tool for addressing many of the challenges facing the precision cosmology program.« less
  • We present the cosmic shear signal predicted by two viable cosmological models in the framework of modified-action f(R) theories. We use f(R) models where the current accelerated expansion of the Universe is a direct consequence of the modified gravitational Lagrangian rather than Dark Energy (DE), either in the form of vacuum energy/cosmological constant or of a dynamical scalar field (e.g. quintessence). We choose Starobinsky's (St) and Hu and Sawicki's (HS) f(R) models, which are carefully designed to pass the Solar System gravity tests. In order to further support — or rule out — f(R) theories as alternative candidates to themore » DE hypothesis, we exploit the power of weak gravitational lensing, specifically of cosmic shear. We calculate the tomographic shear matrix as it would be measured by the upcoming ESA Cosmic Vision Euclid satellite. We find that in the St model the cosmic shear signal is almost completely degenerate with ΛCDM, but it is easily distinguishable in the HS model. Moreover, we compute the corresponding Fisher matrix for both the St and HS models, thus obtaining forecasts for their cosmological parameters. Finally, we show that the Bayes factor for cosmic shear will definitely favour the HS model over ΛCDM if Euclid measures a value larger than ∼ 0.02 for the extra HS parameter n{sub HS}.« less
  • The Euclid space mission, designed to probe evolution of the Dark Energy (DE), will map a large area of the sky at three adjacent near-IR filters, Y, J, and H. This coverage will also enable mapping source-subtracted cosmic infrared background (CIB) fluctuations with unprecedented accuracy on sub-degree angular scales. Here, we propose methodology, using the Lyman-break tomography applied to the Euclid -based CIB maps, to accurately isolate the history of CIB emissions as a function of redshift from 10 ≲ z ≲ 20 and to identify the baryonic acoustic oscillations (BAOs) at those epochs. To identify the BAO signature, wemore » would assemble individual CIB maps over conservatively large contiguous areas of ≳400 deg{sup 2}. The method can isolate the CIB spatial spectrum by z to sub-percent statistical accuracy. We illustrate this with a specific model of CIB production at high z normalized to reproduce the measured Spitzer -based CIB fluctuation. We show that even if the latter contains only a small component from high- z sources, the amplitude of that component can be accurately isolated with the methodology proposed here and the BAO signatures at z ≳ 10 are recovered well from the CIB fluctuation spatial spectrum. Probing the BAO at those redshifts will be an important test of the underlying cosmological paradigm and would narrow the overall uncertainties on the evolution of cosmological parameters, including the DE. Similar methodology is applicable to the planned WFIRST mission, where we show that a possible fourth near-IR channel at ≥2 μ m would be beneficial.« less
  • Here, we present a tomographic cosmic shear study from the Deep Lens Survey (DLS), which, providing a limiting magnitudemore » $${r}_{\mathrm{lim}}\sim 27$$ ($$5\sigma $$), is designed as a precursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent nonlinear model to account for the astrophysical systematics originating from intrinsic alignments of galaxy shapes. We find that the cosmological leverage of the DLS is among the highest among existing $$\gt 10$$ deg2 cosmic shear surveys. Combining the DLS tomography with the 9 yr results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives $${{\rm{\Omega }}}_{m}={0.293}_{-0.014}^{+0.012}$$, $${\sigma }_{8}={0.833}_{-0.018}^{+0.011}$$, $${H}_{0}={68.6}_{-1.2}^{+1.4}\;{\text{km s}}^{-1}\;{{\rm{Mpc}}}^{-1}$$, and $${{\rm{\Omega }}}_{b}=0.0475\pm 0.0012$$ for ΛCDM, reducing the uncertainties of the WMAP9-only constraints by ~50%. When we do not assume flatness for ΛCDM, we obtain the curvature constraint $${{\rm{\Omega }}}_{k}=-{0.010}_{-0.015}^{+0.013}$$ from the DLS+WMAP9 combination, which, however, is not well constrained when WMAP9 is used alone. The dark energy equation-of-state parameter w is tightly constrained when baryonic acoustic oscillation (BAO) data are added, yielding $$w=-{1.02}_{-0.09}^{+0.10}$$ with the DLS+WMAP9+BAO joint probe. The addition of supernova constraints further tightens the parameter to $$w=-1.03\pm 0.03$$. Our joint constraints are fully consistent with the final Planck results and also with the predictions of a ΛCDM universe.« less
  • A goal of forthcoming imaging surveys is to use weak gravitational lensing shear measurements to constrain dark energy. A challenge to this program is that redshifts to the lensed, source galaxies must be determined using photometric, rather than spectroscopic, information. We quantify the importance of uncalibrated photometric redshift outliers to the dark energy goals of forthcoming imaging surveys in a manner that does not assume any particular photometric redshift technique or template. In so doing, we provide an approximate blueprint for computing the influence of specific outlier populations on dark energy constraints. We find that outlier populations whose photo-z distributionsmore » are tightly localized about a significantly biased redshift must be controlled to a per-galaxy rate of (1-3) x 10{sup -3} to insure that systematic errors on dark energy parameters are rendered negligible. In the complementary limit, a subset of imaged galaxies with uncalibrated photometric redshifts distributed over a broad range must be limited to fewer than a per-galaxy error rate of F {sub cat} {approx}< (2-4) x 10{sup -4}. Additionally, we explore the relative importance of calibrating the photo-z's of a core set of relatively well-understood galaxies as compared to the need to identify potential catastrophic photo-z outliers. We discuss the degradation of the statistical constraints on dark energy parameters induced by excising source galaxies at high- and low-photometric redshifts, concluding that removing galaxies with photometric redshifts z {sup ph} {approx}> 2.4 and z {sup ph} {approx}< 0.3 may mitigate damaging catastrophic redshift outliers at a relatively small ({approx}<20%) cost in statistical error. In an Appendix, we show that forecasts for the degradation in dark energy parameter constraints due to uncertain photometric redshifts depend sensitively on the treatment of the nonlinear matter power spectrum. In particular, previous work using Peacock and Dodds may have overestimated the photo-z calibration requirements of future surveys.« less