Weak Lensing Power Spectrum Research Articles

Constraining holographic dark energy and analyzing cosmological tensions

We investigate cosmological constraints on the holographic dark energy (HDE) using the state-of-the-art cosmological datasets: Planck CMB angular power spectra and weak lensing power spectra, Atacama Cosmology Telescope (ACT) temperature power spectra, baryon acoustic oscillation (BAO) and redshift-space distortion (RSD) measurements from six-degree-field galaxy survey and Sloan Digital Sky Survey (DR12 & DR16) and the Cepheids-Supernovae measurement from SH0ES team (R22). We also examine the HDE model and ΛCDM with and without Neff (effective number of relativistic species) being treated as a free parameter. We find that the HDE model can relieve the tensions of H0 and S8 to certain degrees. With “Planck+ACT+BAO+RSD” datasets, the constraints are H0=69.70±1.39km s−1Mpc−1 and S8=0.823±0.011 in HDE model, which brings down the Hubble tension down to 1.92σ confidence level (C.L.) and the S8 tension to (1−2)σ C.L. By adding the R22 data, their values are improved as H0=71.86±0.93km s−1Mpc−1 and S8=0.813±0.010, which further brings the Hubble tension down to 0.85σ C.L. and relieves the S8 tension. We also quantify the goodness-of-fit of different models with Akaike information criterion (AIC) and Bayesian information criterion (BIC), and find that the HDE agrees with the observational data better than the ΛCDM and other extended models (treating Neff as free for fitting).

On the sensitivity of weak gravitational lensing to the cosmic expansion function

ABSTRACT We analyse the functional derivative of the cosmic-shear power spectrum $C_\ell ^\gamma$ with respect to the cosmic expansion function. Our interest in doing so is two-fold: (i) In view of attempts to detect minor changes of the cosmic expansion function that may be due to a possibly time-dependent dark-energy density, we wish to know how sensitive the weak-lensing power spectrum is to changes in the expansion function. (ii) In view of recent empirical determinations of the cosmic expansion function from distance measurements, independent of specific cosmological models, we wish to find out how uncertainties in the expansion function translate to uncertainties in the cosmic-shear power spectrum. We find the following answers: relative changes of the expansion function are amplified by the cosmic-shear power spectrum by a factor ≈2–6, weakly depending on the scale factor where the change is applied, and the current uncertainty of one example for an empirically determined expansion function translates to a relative uncertainty of the cosmic-shear power spectrum of $\approx 10~{{\ \rm per\ cent}}$.

Thermal Energy Census with the Sunyaev–Zel’dovich Effect of DESI Galaxy Clusters/Groups and Its Implication on the Weak-lensing Power Spectrum

We carry out a thermal energy census of hot baryons at z < 1, by cross correlating the Planck Modified Internal Linear Combination Algorithm (MILCA) y map with 0.8 million clusters/groups selected from the Yang et al. catalog. The thermal Sunyaev–Zel’dovich effect around these clusters/groups is reliably obtained, which enables us to make our model constraints based on one-halo (1h) and two-halo (2h) contributions, respectively. (1) The total measurement signal-to-noise (S/N) of the one-halo term is 63. We constrain the Y–M relation over the halo mass range of 1013–1015 M ⊙ h −1, and find Y ∝ M α with α = 1.8 at z = 0.14 (α = 2.1 at z = 0.75). The total thermal energy of gas bound to clusters/groups increases from 0.1 meV cm−3 at z = 0.14 to 0.22 meV cm−3 at z = 0.75. (2) The 2h term is used to constrain the bias-weighted electron pressure 〈b y P e 〉. We find that 〈b y P e 〉 (in units of meV cm−3) increases from 0.24 ± 0.02 at z = 0.14 to 0.45 ± 0.02 at z = 0.75. These results lead to several implications. (i) The hot gas fraction f gas in clusters/groups monotonically increase with the halo mass, where f gas of a 1014 M ⊙ h −1 halo is ∼50% (25%) of the cosmic mean at z = 0.14 (0.75). (ii) By comparing the 1h and 2h terms, we obtain a tentative constraint on the thermal energy of unbound gas. (iii) The above results lead to significant suppression of the matter and weak-lensing power spectrum at small scales. These implications are important for astrophysics and cosmology, and we will further investigate them with improved data and gas modeling.

Correction to: The matter fluctuation amplitude inferred from the weak lensing power spectrum and correlation function in CFHTLenS data

Correction to: The matter fluctuation amplitude inferred from the weak lensing power spectrum and correlation function in CFHTLenS data

Assessing theoretical uncertainties for cosmological constraints from weak lensing surveys

ABSTRACT Weak gravitational lensing is a powerful probe, which is used to constrain the standard cosmological model and its extensions. With the enhanced statistical precision of current and upcoming surveys, high-accuracy predictions for weak lensing statistics are needed to limit the impact of theoretical uncertainties on cosmological parameter constraints. For this purpose, we present a comparison of the theoretical predictions for the non-linear matter and weak lensing power spectra, based on the widely used fitting functions ($\texttt {mead}$ and $\texttt {rev-halofit}$ ), emulators ($\texttt {EuclidEmulator}$ , $\texttt {EuclidEmulator2}$ , $\texttt {BaccoEmulator}$ , and $\texttt {CosmicEmulator}$ ), and N-body simulations (pkdgrav3). We consider the forecasted constraints on the $\Lambda \texttt {CDM}$ and $\texttt {wCDM}$ models from weak lensing for stage III and stage IV surveys. We study the relative bias on the constraints and their dependence on the assumed prescriptions. Assuming a $\Lambda \texttt {CDM}$ cosmology, we find that the relative agreement on the S8 parameter is between 0.2 and 0.3σ for a stage III-like survey between the above predictors. For a stage IV-like survey the agreement becomes 1.4–3.0σ. In the $\texttt {wCDM}$ scenario, we find broader S8 constraints, and agreements of 0.18–0.26σ and 0.7–1.7σ for stage III and stage IV surveys, respectively. The accuracies of the above predictors therefore appear adequate for stage III surveys, whereas the fitting functions would need improvements for future stage IV surveys. Furthermore, we find that, of the fitting functions, $\texttt {mead}$ provides the best agreement with the emulators. We discuss the implication of these findings for the preparation of future weak lensing surveys, and the relative impact of theoretical uncertainties to other systematics.

Almanac: Weak Lensing power spectra and map inference on the masked sphere

We present a field-based signal extraction of weak lensing from noisy observations on the curved and masked sky. We test the analysis on a simulated Euclid-like survey, using a Euclid-like mask and noise level. To make optimal use of the information available in such a galaxy survey, we present a Bayesian method for inferring the angular power spectra of the weak lensing fields, together with an inference of the noise-cleaned tomographic weak lensing shear and convergence (projected mass) maps. The latter can be used for field-level inference with the aim of extracting cosmological parameter information including non-gaussianity of cosmic fields. We jointly infer all-sky $E$-mode and $B$-mode tomographic auto- and cross-power spectra from the masked sky, and potentially parity-violating $EB$-mode power spectra, up to a maximum multipole of $\ell_{\rm max}=2048$. We use Hamiltonian Monte Carlo sampling, inferring simultaneously the power spectra and denoised maps with a total of $\sim 16.8$ million free parameters. The main output and natural outcome is the set of samples of the posterior, which does not suffer from leakage of power from $E$ to $B$ unless reduced to point estimates. However, such point estimates of the power spectra, the mean and most likely maps, and their variances and covariances, can be computed if desired.

Non-Gaussian likelihood of weak lensing power spectra

The power spectrum of weak lensing fluctuations has a non-Gaussian distribution due to its quadratic nature. On small scales the Central Limit Theorem acts to Gaussianize this distribution but non-Gaussianity in the signal due to gravitational collapse is increasing and the functional form of the likelihood is unclear. Analyses have traditionally assumed a Gaussian likelihood with non-linearity incorporated into the covariance matrix; here we provide the theory underpinning this assumption. We calculate, for the first time, the leading-order correction to the distribution of angular power spectra from non-Gaussianity in the underlying signal and study the transition to Gaussianity. Our expressions are valid for an arbitrary number of correlated maps and correct the Wishart distribution in the presence of weak (but otherwise arbitrary) non-Gaussianity in the signal. Surprisingly, the resulting distribution is not equivalent to an Edgeworth expansion. The leading-order effect is to broaden the covariance matrix by the usual trispectrum term, with residual skewness sourced by the trispectrum and the square of the bispectrum. Using lognormal lensing maps we demonstrate that our likelihood is uniquely able to model both large and mildly non-linear scales. We provide easy-to-compute statistics to quantify the size of the non-Gaussian corrections. We show that the full non-Gaussian likelihood can be accurately modelled as a Gaussian on small, non-linear scales. On large angular scales non-linearity in the lensing signal imparts a negligible correction to the likelihood, which takes the Wishart form in the full-sky case. Our formalism is equally applicable to any kind of projected field.

Kernel-based emulator for the 3D matter power spectrum from CLASS

The 3D matter power spectrum, $P_{\delta}(k,z)$ is a fundamental quantity in the analysis of cosmological data such as large-scale structure, 21cm observations, and weak lensing. Existing computer models (Boltzmann codes) such as CLASS can provide it at the expense of immoderate computational cost. In this paper, we propose a fast Bayesian method to generate the 3D matter power spectrum, for a given set of wavenumbers, $k$ and redshifts, $z$. Our code allows one to calculate the following quantities: the linear matter power spectrum at a given redshift (the default is set to 0); the non-linear 3D matter power spectrum with/without baryon feedback; the weak lensing power spectrum. The gradient of the 3D matter power spectrum with respect to the input cosmological parameters is also returned and this is useful for Hamiltonian Monte Carlo samplers. The derivatives are also useful for Fisher matrix calculations. In our application, the emulator is accurate when evaluated at a set of cosmological parameters, drawn from the prior, with the fractional uncertainty, $\Delta P_{\delta}/P_{\delta}$ centred on 0. It is also $\sim 300$ times faster compared to CLASS, hence making the emulator amenable to sampling cosmological and nuisance parameters in a Monte Carlo routine. In addition, once the 3D matter power spectrum is calculated, it can be used with a specific redshift distribution, $n(z)$ to calculate the weak lensing and intrinsic alignment power spectra, which can then be used to derive constraints on cosmological parameters in a weak lensing data analysis problem. The software ($\texttt{emuPK}$) can be trained with any set of points and is distributed on Github, and comes with a pre-trained set of Gaussian Process (GP) models, based on 1000 Latin Hypercube (LH) samples, which follow roughly the current priors for current weak lensing analyses.

First Evidence of Intrinsic Alignments of Red Galaxies at z &gt; 1: Cross Correlation between CFHTLenS and FastSound Samples

Abstract We report the first evidence for intrinsic alignment (IA) of red galaxies at z &gt; 1. We measure the gravitational shear-intrinsic ellipticity cross correlation function at z ∼ 1.3 using galaxy positions from the FastSound spectroscopic survey and galaxy shapes from the Canada France Hawaii telescope lensing survey data. Adopting the nonlinear alignment model, we obtain a 2.4σ level detection of the IA amplitude A LA = 27.48 − 11.54 + 11.53 (and 2.6σ with A LA = 29.43 − 11.49 + 11.48 when weak lensing contaminations are taken into account), which is larger than the value extrapolated from the constraints obtained at lower redshifts. Our measured IA is translated into a ∼20% contamination of the weak-lensing power spectrum for the red galaxies. This marginal detection of IA for red galaxies at z &gt; 1 motivates the continuing investigation of the nature of IA for weak lensing studies. Furthermore, our result provides the first step to utilize IA measurements in future high-z surveys as a cosmological probe, complementary to galaxy clustering and lensing.

Mixed dark matter: matter power spectrum and halo mass function

We investigate and quantify the impact of mixed (cold and warm) dark matter models on large-scale structure observables. In this scenario, dark matter comes in two phases, a cold one (CDM) and a warm one (WDM): the presence of the latter causes a suppression in the matter power spectrum which is allowed by current constraints and may be detected in present-day and upcoming surveys. We run a large set of N-body simulations in order to build an efficient and accurate emulator to predict the aforementioned suppression with percent precision over a wide range of values for the WDM mass, Mwdm, and its fraction with respect to the totality of dark matter, fwdm. The suppression in the matter power spectrum is found to be independent of changes in the cosmological parameters at the 2% level for k≲ 10 h/Mpc and z≤ 3.5. In the same ranges, by applying a baryonification procedure on both ΛCDM and CWDM simulations to account for the effect of feedback, we find a similar level of agreement between the two scenarios. We examine the impact that such suppression has on weak lensing and angular galaxy clustering power spectra. Finally, we discuss the impact of mixed dark matter on the shape of the halo mass function and which analytical prescription yields the best agreement with simulations. We provide the reader with an application to galaxy cluster number counts.

Primordial Power Spectrum reconstruction from CMB Weak Lensing Power Spectrum

We use the modified and improved Richardson-Lucy (IRL) deconvolutionalgorithm to reconstruct the Primordial Power Spectrum (PPS) from the Weak Lensing Power Spectrum CL ϕϕ reconstructed from CMB anisotropies. This provides an independent window to observe and constrain the PPS PR (k) along different k scales as compared to CMB Temperature Power Spectrum. The Weak Lensing Power Spectrum does not contain secondary variations in power and hence is cleaner, unlike the Temperature Power Spectrum which suffers from lensing which must be addressed during PPS reconstructions. We demonstrate that the physical behaviour of the weak lensing kernel is unique and reconstructs broad features over k. We provide an in-depth analysis of the error propagation using simulated data and Monte-Carlo sampling, using Planck best-fit cosmological parameters to simulate the data with cosmic variance limited error bars. The error and initial condition analyses provide a clear picture of the optimal reconstruction region for the estimator while providing a detailed statistical insight of the results. We also provide an algorithm for PR (k) sampling sparsity to be used based on the given data and errors, to optimize statistical significance. Eventually we plan to use this method on actual mission data and provide a cross reference to PPS reconstructed from other sectors and any possible features in them.

The unequal-time matter power spectrum: impact on weak lensing observables

We investigate the impact of a common approximation of weak lensing power spectra: the use of single-epoch matter power spectra in integrals over redshift. We disentangle this from the closely connected Limber's approximation. We derive the unequal-time matter power spectrum at one-loop in standard perturbation theory and effective field theory to deal with non-linear physics. We compare these formalisms and conclude that the unequal-time power spectrum using effective field theory breaks for larger scales. As an alternative we introduce the midpoint approximation. We also provide, for the first time, a fitting function for the time evolution of the effective field theory counterterms based on the Quijote simulations. Then we compute the angular power spectrum using a range of approaches: the Limber approximation, and the geometric and midpoint approximations. We compare our results with the exact calculation at all angular scales using the unequal-time power spectrum. We use DES Y1 and LSST-like redshift distributions for our analysis. We find that the use of the Limber's approximation in weak lensing diverges from the exact calculation of the angular power spectrum on large-angle separations, ℓ < 10. Even though this deviation is of order 2% maximum for cosmic lensing, we find the biggest effect for galaxy clustering and galaxy-galaxy lensing. We show that not only is this true for upcoming galaxy surveys, but also for current data such as DES Y1. Finally, we make our pipeline and analysis publicly available as a Python package called unequalpy.

Amending the halo model to satisfy cosmological conservation laws

One of the most powerful tools in the arsenal of theoretical cosmologists is the halo model of large scale structure, which provides a phenomenological description of nonlinear structure in our universe. However, it is well known that there is no simple way to impose conservation laws in the halo model. This can severely impair the predictions on large scales for observables such as weak lensing or the kinematic Sunyaev-Zel'dovich effect, which should satisfy mass and momentum conservations, respectively. For example, the standard halo model overpredicts weak lensing power spectrum by $> 8\%$ on scales $>20$ degrees. To address this problem, we present an $Amended ~ Halo ~ Model$, explicitly separating the linear perturbations from $compensated$ halo profiles. This is guaranteed to respect conservation laws, as well as linear theory predictions on large scales. We then provide a simple fitting function for the compensated halo profiles, and discuss the modified predictions for 1-halo and 2-halo terms, as well as other cosmological observations such weak lensing power spectrum. Furthermore, we argue that the amended halo model provides a more accurate framework to capture physical effects that happen in the process of cosmological structure formation.

The matter fluctuation amplitude inferred from the weak lensing power spectrum and correlation function in CFHTLenS data

ABSTRACT Based on the cosmic shear data from the Canada–France–Hawaii Telescope Lensing Survey (CFHTLenS), Kilbinger et al. obtained a constraint on the amplitude of matter fluctuations of σ8(Ωm/0.27)0.6 = 0.79 ± 0.03 from the two-point correlation function (2PCF). This is ≈3σ lower than the value 0.89 ± 0.01 derived from Planck data on cosmic microwave background (CMB) anisotropies. On the other hand, based on the same CFHTLenS data, but using the power spectrum, and performing a different analysis, Liu et al. obtained the higher value of $\sigma _8(\Omega _\mathrm{m}/0.27)^{0.64}=0.87^{+0.05}_{-0.06}$. We here investigate the origin of this difference, by performing a fair side-by-side comparison of the 2PCF and power spectrum analyses on CFHTLenS data. We find that these two statistics indeed deliver different results, even when applied to the same data in an otherwise identical procedure. We identify excess power in the data on small scales (ℓ &amp;gt; 5000) driving the larger values inferred from the power spectrum. We speculate on the possible origin of this excess small-scale power. More generally, our results highlight the utility of analysing the 2PCF and the power spectrum in tandem, to discover (and to help control) systematic errors.

Forecasting super-sample covariance in future weak lensing surveys with SuperSCRAM

The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to super-sample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0.

Weak-lensing Power Spectrum Reconstruction by Counting Galaxies. II. Improving the ABS Method with the Shift Parameter

Abstract In Paper I of this series, we proposed an analytical method of blind separation (ABS) to extract the cosmic magnification signal in galaxy number distribution and reconstruct the weak-lensing power spectrum. Here, we report a new version of the ABS method with significantly improved performance. This version is characterized by a shift parameter, , with the special case of corresponding to the original ABS method. We have tested this new version, compared it with the previous one, and confirmed its superior performance in all investigated situations. Therefore, it supercedes the previous version. The proof of concept studies presented in this paper demonstrate that it may enable surveys such as LSST and SKA to reconstruct the lensing power spectrum at z ≃ 1 with 1% accuracy. We will test the new ABS method in more realistic simulations to verify its applicability to real data.

Flat-Sky Pseudo-Cls Analysis for Weak Gravitational Lensing

We investigate the use of estimators of weak lensing power spectra based on a flat-sky implementation of the Pseudo-Cl (PCl) technique, where the masked shear field is transformed without regard for masked regions of sky. This masking mixes power, and E-convergence and B-modes. To study the accuracy of forward-modelling and full-sky power spectrum recovery we consider both large-area survey geometries, and small-scale masking due to stars and a checkerboard model for field-of-view gaps. The power spectrum for the large-area survey geometry is sparsely-sampled and highly oscillatory, which makes modelling problematic. Instead, we derive an overall calibration for large-area mask bias using simulated fields. The effects of small-area star masks can be accurately corrected for, while the checkerboard mask has oscillatory and spiky behaviour which leads to percent biases. Apodisation of the masked fields leads to increased biases and a loss of information. We find that we can construct an unbiased forward-model of the raw PCls, and recover the full-sky convergence power to within a few percent accuracy for both Gaussian and lognormal-distributed shear fields. Propagating this through to cosmological parameters using a Fisher-Matrix formalism, we find we can make unbiased estimates of parameters for surveys up to 1,200 deg$^2$ with 30 galaxies per arcmin$^2$, beyond which the percent biases become larger than the statistical accuracy. This implies a flat-sky PCl analysis is accurate for current surveys but a Euclid-like survey will require higher accuracy.

Noise estimates for measurements of weak lensing from the Ly α forest

We have proposed a method for measuring weak lensing using the Lyman-alpha forest. Here we estimate the noise expected in weak lensing maps and power spectra for different sets of observational parameters. We find that surveys of the size and quality of the ones being done today and ones planned for the future will be able to measure the lensing power spectrum at a source redshift of z~2.5 with high precision and even be able to image the distribution of foreground matter with high fidelity on degree scales. For example, we predict that Lyman-alpha forest lensing measurement from the Dark Energy Spectroscopic Instrument survey should yield the mass fluctuation amplitude with statistical errors of 1.5%. By dividing the redshift range into multiple bins some tomographic lensing information should be accessible as well. This would allow for cosmological lensing measurements at higher redshift than are accessible with galaxy shear surveys and correspondingly better constraints on the evolution of dark energy at relatively early times.

Precision calculations of the cosmic shear power spectrum projection

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.

Weak-lensing Power Spectrum Reconstruction by Counting Galaxies. I. The ABS Method

Abstract We propose an analytical method of blind separation (ABS) of cosmic magnification from the intrinsic fluctuations of galaxy number density in the observed galaxy number density distribution. The ABS method utilizes the different dependences of the signal (cosmic magnification) and contamination (galaxy intrinsic clustering) on galaxy flux to separate the two. It works directly on the measured cross-galaxy angular power spectra between different flux bins. It determines/reconstructs the lensing power spectrum analytically, without assumptions of galaxy intrinsic clustering and cosmology. It is unbiased in the limit of an infinite number of galaxies. In reality, the lensing reconstruction accuracy depends on survey configurations, galaxy biases, and other complexities due to a finite number of galaxies and the resulting shot noise fluctuations in the cross-galaxy power spectra. We estimate its performance (systematic and statistical errors) in various cases. We find that stage IV dark energy surveys such as Square Kilometre Array and Large Synoptic Survey Telescope are capable of reconstructing the lensing power spectrum at and accurately. This lensing reconstruction only requires counting galaxies and is therefore highly complementary to cosmic shear measurement by the same surveys.