We present here the cosmo-SLICS, a new suite of simulations specially designed for the analysis of current and upcoming weak lensing data beyond the standard two-point cosmic shear. We sampled the ...Ωm, σ8, h, w0 parameter space at 25 points organised in a Latin hyper-cube, spanning a range that contains most of the 2σ posterior distribution from ongoing lensing surveys. At each of these nodes we evolved a pair of N-body simulations in which the sampling variance is highly suppressed, and ray-traced the volumes 800 times to further increase the effective sky coverage. We extracted a lensing covariance matrix from these pseudo-independent light-cones and show that it closely matches a brute-force construction based on an ensemble of 800 truly independent N-body runs. More precisely, a Fisher analysis reveals that both methods yield marginalized two-dimensional constraints that vary by less than 6% in area, a result that holds under different survey specifications and that matches to within 15% the area obtained from an analytical covariance calculation. Extending this comparison with our 25 wCDM models, we probed the cosmology dependence of the lensing covariance directly from numerical simulations, reproducing remarkably well the Fisher results from the analytical models at most cosmologies. We demonstrate that varying the cosmology at which the covariance matrix is evaluated in the first place might have an order of magnitude greater impact on the parameter constraints than varying the choice of covariance estimation technique. We present a test case in which we generate fast predictions for both the lensing signal and its associated variance with a flexible Gaussian process regression emulator, achieving an accuracy of a few percent on the former and 10% on the latter.
Abstract
We present measurements of the weak gravitational lensing shear power spectrum based on
$450 \deg ^2$
of imaging data from the Kilo Degree Survey. We employ a quadratic estimator in two and ...three redshift bins and extract band powers of redshift autocorrelation and cross-correlation spectra in the multipole range 76 ≤ ℓ ≤ 1310. The cosmological interpretation of the measured shear power spectra is performed in a Bayesian framework assuming a ΛCDM model with spatially flat geometry, while accounting for small residual uncertainties in the shear calibration and redshift distributions as well as marginalizing over intrinsic alignments, baryon feedback and an excess-noise power model. Moreover, massive neutrinos are included in the modelling. The cosmological main result is expressed in terms of the parameter combination
$S_8 \equiv \sigma _8 \sqrt{\Omega _{\rm m}/0.3}$
yielding S
8 = 0.651 ± 0.058 (three z-bins), confirming the recently reported tension in this parameter with constraints from Planck at 3.2σ (three z-bins). We cross-check the results of the three z-bin analysis with the weaker constraints from the two z-bin analysis and find them to be consistent. The high-level data products of this analysis, such as the band power measurements, covariance matrices, redshift distributions and likelihood evaluation chains are available at http://kids.strw.leidenuniv.nl.
We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450 deg super( 2) of imaging data from the Kilo Degree Survey (KiDS). For a flat ... cold dark ...matter (...CDM) cosmology with a prior on H sub( 0) that encompasses the most recent direct measurements, we find S sub( 8) ... = 0.745 plus or minus 0.039. This result is in good agreement with other low-redshift probes of large-scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A 2.3... tension in S sub( 8) and 'substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved 'self-calibrating' version of lensfit validated using an extensive suite of image simulations. Four-band ugri photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov chains are available at http://kids.strw.leidenuniv.nl. (ProQuest: ... denotes formulae/symbols omitted.)
Correlations between the intrinsic shapes of galaxies and the large-scale galaxy density field provide an important tool to investigate galaxy intrinsic alignments, which constitute the major ...potential astrophysical systematic in cosmological weak lensing (cosmic shear) surveys, but also yield insight into the formation and evolution of galaxies. We measure galaxy position-shape correlations in the MegaZ-LRG sample for more than 800 000 luminous red galaxies for comoving transverse separations of 0.3 < rp < 60 h-1 Mpc, making the first such measurement with a photometric redshift sample. In combination with a re-analysis of several spectroscopic SDSS samples, we constrain an intrinsic alignment model for early-type galaxies over long baselines in redshift (z ≲ 0.7) and luminosity (4 mag) with high statistical precision. We develop and test the formalism to incorporate photometric redshift scatter in the modelling of these observations. For rp > 6 h-1 Mpc, the fits to galaxy position-shape correlation functions are consistent with the scaling with rp and redshift of a revised, nonlinear version of the linear alignment model (Hirata & Seljak 2004) for all samples. An extra redshift dependence ∝ (1 + z)ηother is constrained to ηother = −0.3 ± 0.8 (1σ). To obtain consistent amplitudes for all data, an additional dependence on galaxy luminosity ∝ Lβ with $\beta=1.1^{+0.3}_{-0.2}$β=1.1-0.2+0.3 is required. The normalisation of the intrinsic alignment power spectrum is found to be $(0.077 \pm 0.008)\, \rho_{\rm cr}^{-1}$(0.077±0.008) ρcr-1 for galaxies at redshift 0.3 and r band magnitude of − 22 (k- and evolution-corrected to z = 0). Assuming zero intrinsic alignments for blue galaxies, we assess the bias on cosmological parameters for a tomographic CFHTLS-like lensing survey given our new constraints on the intrinsic alignment model parameter space. Both the resulting mean bias and its uncertainty are smaller than the 1σ statistical errors when using the constraints from all samples combined. The addition of MegaZ-LRG data is critical to achieving constraints this strong, reducing the uncertainty in intrinsic alignment bias on cosmological parameters by factors of three to seven.
Aims. Cosmic shear is a powerful method to constrain cosmology, provided that any systematic effects are under control. The intrinsic alignment of galaxies is expected to severely bias parameter ...estimates if not taken into account. We explore the potential of a joint analysis of tomographic galaxy ellipticity, galaxy number density, and ellipticity-number density cross-correlations to simultaneously constrain cosmology and self-calibrate unknown intrinsic alignment and galaxy bias contributions. Methods. We treat intrinsic alignments and galaxy biasing as free functions of scale and redshift and marginalise over the resulting parameter sets. Constraints on cosmology are calculated by combining the likelihoods from all two-point correlations between galaxy ellipticity and galaxy number density. The information required for these calculations is already available in a standard cosmic shear data set. We include contributions to these functions from cosmic shear, intrinsic alignments, galaxy clustering and magnification effects. Results. In a Fisher matrix analysis we compare our constraints with those from cosmic shear alone in the absence of intrinsic alignments. For a potential future large area survey, such as Euclid, the extra information from the additional correlation functions can make up for the additional free parameters in the intrinsic alignment and galaxy bias terms, depending on the flexibility in the models. For example, the dark energy task force figure of merit is recovered even when more than 100 free parameters are marginalised over. We find that the redshift quality requirements are similar to those calculated in the absence of intrinsic alignments.
ABSTRACT
We present constraints on Horndeski gravity from a combined analysis of cosmic shear, galaxy–galaxy lensing and galaxy clustering from $450\, \mathrm{deg}^2$ of the Kilo-Degree Survey and ...the Galaxy And Mass Assembly survey.The Horndeski class of dark energy/modified gravity models includes the majority of universally coupled extensions to ΛCDM with one scalar field in addition to the metric. We study the functions of time that fully describe the evolution of linear perturbations in Horndeski gravity. Our results are compatible throughout with a ΛCDM model. By imposing gravitational wave constraints, we fix the tensor speed excess to zero and consider a subset of models including, e.g. quintessence and f(R) theories. Assuming proportionality of the Horndeski functions αB and αM (kinetic braiding and the Planck mass run rate, respectively) to the dark energy density fraction ΩDE(a) = 1 − Ωm(a), we find for the proportionality coefficients $\hat{\alpha }_\mathrm{ B} = 0.20_{-0.33}^{+0.20} \,$ and $\, \hat{\alpha }_\mathrm{ M} = 0.25_{-0.29}^{+0.19}$. Our value of $S_8 \equiv \sigma _8 \sqrt{\Omega _{\mathrm{m}}/0.3}$ is in better agreement with the Planck estimate when measured in the enlarged Horndeski parameter space than in a pure ΛCDM scenario. In our joint three-probe analysis, we report a downward shift of the S8 best-fitting value from the Planck measurement of $\Delta S_8 = 0.016_{-0.046}^{+0.048}$ in Horndeski gravity, compared to $\Delta S_8 = 0.059_{-0.039}^{+0.040}$ in ΛCDM. Our constraints are robust to the modelling uncertainty of the non-linear matter power spectrum in Horndeski gravity. Our likelihood code for multiprobe analysis in both ΛCDM and Horndeski gravity is publicly available at https://github.com/alessiospuriomancini/KiDSHorndeski.
We present a combined tomographic weak gravitational lensing analysis of the Kilo Degree Survey (KV450) and the Dark Energy Survey (DES-Y1). We homogenize the analysis of these two public cosmic ...shear datasets by adopting consistent priors and modeling of nonlinear scales, and determine new redshift distributions for DES-Y1 based on deep public spectroscopic surveys. Adopting these revised redshifts results in a 0.8
σ
reduction in the DES-inferred value for
S
8
, which decreases to a 0.5
σ
reduction when including a systematic redshift calibration error model from mock DES data based on the MICE2 simulation. The combined KV450+DES-Y1 constraint on
S
8
= 0.762
−0.024
+0.025
is in tension with the
Planck
2018 constraint from the cosmic microwave background at the level of 2.5
σ
. This result highlights the importance of developing methods to provide accurate redshift calibration for current and future weak-lensing surveys.
We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey. This is the first time that a full optical to near-infrared ...data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning 450 deg
2
, allows us to significantly improve the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and – most importantly – to solidify our knowledge of the redshift distributions of the sources. Based on a flat ΛCDM model we find
S
8
≡ σ
8
Ω
m
/0.3 = 0.737
+0.040
−0.036
in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with
S
8
differing by 2.3
σ
. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with
Planck
is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.
We present redshift distribution estimates of galaxies selected from the fourth data release of the Kilo-Degree Survey over an area of ∼1000 deg
2
(KiDS-1000). These redshift distributions represent ...one of the crucial ingredients for weak gravitational lensing measurements with the KiDS-1000 data. The primary estimate is based on deep spectroscopic reference catalogues that are re-weighted with the help of a self-organising map (SOM) to closely resemble the KiDS-1000 sources, split into five tomographic redshift bins in the photometric redshift range 0.1 <
z
B
≤ 1.2. Sources are selected such that they only occupy that volume of nine-dimensional magnitude-space that is also covered by the reference samples (‘gold’ selection). Residual biases in the mean redshifts determined from this calibration are estimated from mock catalogues to be ≲0.01 for all five bins with uncertainties of ∼0.01. This primary SOM estimate of the KiDS-1000 redshift distributions is complemented with an independent clustering redshift approach. After validation of the clustering-
z
on the same mock catalogues and a careful assessment of systematic errors, we find no significant bias of the SOM redshift distributions with respect to the clustering-
z
measurements. The SOM redshift distributions re-calibrated by the clustering-
z
represent an alternative calibration of the redshift distributions with only slightly larger uncertainties in the mean redshifts of ∼0.01 − 0.02 to be used in KiDS-1000 cosmological weak lensing analyses. As this includes the SOM uncertainty, clustering-
z
are shown to be fully competitive on KiDS-1000 data.