No Arabic abstract
In this paper we derive a full expression for the propagation of weak lensing shape measurement biases into cosmic shear power spectra including the effect of missing data. We show using simulations that terms higher than first order in bias parameters can be ignored and the impact of biases can be captured by terms dependent only on the mean of the multiplicative bias field. We identify that the B-mode power contains information on the multiplicative bias. We find that without priors on the residual multiplicative bias $delta m$ and stochastic ellipticity variance $sigma_e$ that constraints on the amplitude of the cosmic shear power spectrum are completely degenerate, and that when applying priors the constrained amplitude $A$ is slightly biased low via a classic marginalisation paradox. Using all-sky Gaussian random field simulations we find that the combination of $(1+2delta m)A$ is unbiased for a joint EE and BB power spectrum likelihood if the error and mean (precision and accuracy) of the stochastic ellipticity variance is known to better than $sigma(sigma_e)leq 0.05$ and $Deltasigma_eleq 0.01$, or the multiplicative bias is known to better than $sigma(m)leq 0.07$ and $Delta mleq 0.01$.
We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3leq z leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 leq ell leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $Lambda$ cold dark matter ($Lambda$CDM) model, we find $S_8equiv sigma_8(Omega_{rm m}/0.3)^alpha=0.800^{+0.029}_{-0.028}$ for $alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $sim 0.6sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.
Upcoming measurements of the small-scale primary cosmic microwave background (CMB) temperature and polarization power spectra ($TT$/$TE$/$EE$) are anticipated to yield transformative constraints on new physics, including the effective number of relativistic species in the early universe ($N_{rm eff}$). However, at multipoles $ell gtrsim 3000$, the primary CMB power spectra receive significant contributions from gravitational lensing. While these modes still carry primordial information, their theoretical modeling requires knowledge of the CMB lensing convergence power spectrum, $C_L^{kappakappa}$, including on small scales where it is affected by nonlinear gravitational evolution and baryonic feedback processes. Thus, the high-$ell$ primary CMB is sensitive to these late-time, nonlinear effects. Here, we show that inaccuracies in the modeling of $C_L^{kappakappa}$ can yield surprisingly large biases on cosmological parameters inferred from the primary CMB power spectra measured by the upcoming Simons Observatory and CMB-S4 experiments. For CMB-S4, the biases can be as large as $1.6sigma$ on the Hubble constant $H_0$ in a fit to $Lambda$CDM and $1.2sigma$ on $N_{rm eff}$ in a fit to $Lambda$CDM+$N_{rm eff}$. We show that these biases can be mitigated by explicitly discarding all $TT$ data at $ell>3000$ or by marginalizing over parameters describing baryonic feedback processes, both at the cost of slightly larger error bars. We also discuss an alternative, data-driven mitigation strategy based on delensing the CMB $T$ and $E$-mode maps. Finally, we show that analyses of upcoming data will require Einstein-Boltzmann codes to be run with much higher numerical precision settings than is currently standard, so as to avoid similar -- or larger -- parameter biases due to inaccurate theoretical predictions.
We compute the spherical-sky weak-lensing power spectrum of the shear and convergence. We discuss various approximations, such as flat-sky, and first- and second- order Limber equations for the projection. We find that the impact of adopting these approximations is negligible when constraining cosmological parameters from current weak lensing surveys. This is demonstrated using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that the reported tension with Planck Cosmic Microwave Background (CMB) temperature anisotropy results cannot be alleviated. For future large-scale surveys with unprecedented precision, we show that the spherical second-order Limber approximation will provide sufficient accuracy. In this case, the cosmic-shear power spectrum is shown to be in agreement with the full projection at the sub-percent level for l > 3, with the corresponding errors an order of magnitude below cosmic variance for all l. When computing the two-point shear correlation function, we show that the flat-sky fast Hankel transformation results in errors below two percent compared to the full spherical transformation. In the spirit of reproducible research, our numerical implementation of all approximations and the full projection are publicly available within the package nicaea at http://www.cosmostat.org/software/nicaea.
With the advent of large-scale weak lensing surveys there is a need to understand how realistic, scale-dependent systematics bias cosmic shear and dark energy measurements, and how they can be removed. Here we describe how spatial variations in the amplitude and orientation of realistic image distortions convolve with the measured shear field, mixing the even-parity convergence and odd-parity modes, and bias the shear power spectrum. Many of these biases can be removed by calibration to external data, the survey itself, or by modelling in simulations. The uncertainty in the calibration must be marginalised over and we calculate how this propagates into parameter estimation, degrading the dark energy Figure-of-Merit. We find that noise-like biases affect dark energy measurements the most, while spikes in the bias power have the least impact, reflecting their correlation with the effect of cosmological parameters. We argue that in order to remove systematic biases in cosmic shear surveys and maintain statistical power effort should be put into improving the accuracy of the bias calibration rather than minimising the size of the bias. In general, this appears to be a weaker condition for bias removal. We also investigate how to minimise the size of the calibration set for a fixed reduction in the Figure-of-Merit. These results can be used to model the effect of biases and calibration on a cosmic shear survey accurately, assess their impact on the measurement of modified gravity and dark energy models, and to optimise surveys and calibration requirements.
We introduce the Generalised Lensing and Shear Spectra GLaSS code which is available for download from https://github.com/astro-informatics/GLaSS It is a fast and flexible public code, written in Python, that computes generalized spherical cosmic shear spectra. The commonly used tomographic and spherical Bessel lensing spectra come as built-in run-mode options. GLaSS is integrated into the Cosmosis modular cosmological pipeline package. We outline several computational choices that accelerate the computation of cosmic shear power spectra. Using GLaSS, we test whether the assumption that using the lensing and projection kernels for a spatially-flat universe -- in a universe with a small amount of spatial curvature -- negligibly impacts the lensing spectrum. We refer to this assumption as The Spatially-Flat Universe Approximation, that has been implicitly assumed in all cosmic shear studies to date. We confirm that The Spatially-Flat Universe Approximation has a negligible impact on Stage IV cosmic shear experiments.