Cosmic Shear Field Research Articles

Synthetic light-cone catalogues of modern redshift and weak lensing surveys with abacussummit

ABSTRACT The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy, and dark matter. However, the development of high-fidelity theoretical models is a necessary stepping stone. Here, we present public high-resolution weak lensing maps on the light-cone, generated using the N-body simulation suite abacussummit, and accompanying weak lensing mock catalogues, tuned to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument. Available in this release are maps of the cosmic shear, deflection angle, and convergence fields at source redshifts ranging from z = 0.15 to 2.45 as well as cosmic microwave background convergence maps for each of the 25 base-resolution simulations ($L_{\rm box} = 2000\, h^{-1}\, {\rm Mpc}$ and Npart = 69123) as well as for the two huge simulations ($L_{\rm box} = 7500\, h^{-1}\, {\rm Mpc}$ and Npart = 86403) at the fiducial abacussummit cosmology. The pixel resolution of each map is 0.21 arcmin, corresponding to a healpix Nside of 16 384. The sky coverage of the base simulations is an octant until z ≈ 0.8 (decreasing to about 1800 deg2 at z ≈ 2.4), whereas the huge simulations offer full-sky coverage until z ≈ 2.2. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper Suprime-Cam data sets. The mock catalogues are validated against theoretical predictions for various clustering and lensing statistics, such as correlation multipoles, galaxy–shear, and shear–shear, showing excellent agreement. All products can be downloaded via a Globus endpoint (see Data Availability section).

Response approach to the integrated shear 3-point correlation function: the impact of baryonic effects on small scales

ABSTRACT The integrated shear 3-point correlation function ζ± is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function ξ± by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate ζ± that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale ζ± is dominated by the non-linear matter bispectrum in the squeezed limit, which can be evaluated accurately using the non-linear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale ζ± measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the non-linear matter power spectrum, thereby permitting to account for these astrophysical effects on ζ± similarly to how they are currently accounted for on ξ±. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to ξ±, a combined $\xi _{\pm }\ \&\ \zeta _{\pm }$ analysis can lead to improvements of order $20\!-\!40{{\ \rm per\ cent}}$ on the constraints of cosmological parameters such as σ8 or the dark energy equation of state parameter w0. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These encouraging results motivate future works on the integrated shear 3-point correlation function towards applications to real survey data.

Mapping Spatially Varying Additive Biases in Cosmic Shear Data

In this paper we address the challenge of extracting maps of spatially varying unknown additive biases from cosmic shear data. This is done by exploiting the isotropy of the cosmic shear field, and the anisotropy of a typical additive bias field, using an autocorrelation discrepancy map; which identifies significant non-Gaussian components of the map. We test this approach using simulations and find that the autocorrelation discrepancy map produces spatially varying features that are indicative of the additive bias field both in amplitude and spatial variation. We then apply this to the Dark Energy Survey Year 1 data, and find evidence for spatially varying additive biases of at most 0.002 on large-scales. The method can be used to empirically inform modelling of the spatially varying additive bias field in any cosmological parameter inference, and can act as a validation test for cosmic shear surveys.

The integrated three-point correlation function of cosmic shear

ABSTRACT We present the integrated three-point shear correlation function iζ± – a higher order statistic of the cosmic shear field – which can be directly estimated in wide-area weak lensing surveys without measuring the full three-point shear correlation function, making this a practical and complementary tool to two-point statistics for weak lensing cosmology. We define it as the one-point aperture mass statistic Map measured at different locations on the shear field correlated with the corresponding local two-point shear correlation function ξ±. Building upon existing work on the integrated bispectrum of the weak lensing convergence field, we present a theoretical framework for computing the integrated three-point function in real space for any projected field within the flat-sky approximation and apply it to cosmic shear. Using analytical formulae for the non-linear matter power spectrum and bispectrum, we model iζ± and validate it on N-body simulations within the uncertainties expected from the sixth year cosmic shear data of the Dark Energy Survey. We also explore the Fisher information content of iζ± and perform a joint analysis with ξ± for two tomographic source redshift bins with realistic shape noise to analyse its power in constraining cosmological parameters. We find that the joint analysis of ξ± and iζ± has the potential to considerably improve parameter constraints from ξ± alone, and can be particularly useful in improving the figure of merit of the dynamical dark energy equation of state parameters from cosmic shear data.

Results of gravitational lensing and primordial gravitational waves from the POLARBEAR experiment

POLARBEAR is a Cosmic Microwave Background radiation (CMB) polarization experiment that is located in the Atacama Desert in Chile. The scientific goals of the experiment are to characterize the B-mode signal from gravitational lensing, as well as to search for B-mode signals created by primordial gravitational waves (PGWs). Polarbear started observations in 2012 and has published a series of results. These include the first measurement of a nonzero B-mode angular auto-power spectrum at sub-degree scales where the dominant signal is gravitational lensing of the CMB. In addition, we have achieved the first measurement of crosscorrelation between the lensing potential, which was reconstructed from the CMB polarization data alone by Polarbear, and the cosmic shear field from galaxy shapes by the Subaru Hyper Suprime-Cam (HSC) survey. In 2014, we installed a continuously rotating half-wave plate (CRHWP) at the focus of the primary mirror to search for PGWs and demonstrated the control of low-frequency noise. We have found that the low-frequency B-mode power in the combined dataset with the Planck high-frequency maps is consistent with Galactic dust foreground, thus placing an upper limit on the tensor-to-scalar ratio of r < 0.90 at the 95% confidence level after marginalizing over the foregrounds.

Evidence for the Cross-correlation between Cosmic Microwave Background Polarization Lensing from Polarbear and Cosmic Shear from Subaru Hyper Suprime-Cam

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

3D cosmic shear: Numerical challenges, 3D lensing random fields generation, and Minkowski functionals for cosmological inference

Cosmic shear - the weak gravitational lensing effect generated by fluctuations of the gravitational tidal fields of the large-scale structure - is one of the most promising tools for current and future cosmological analyses. The spherical-Bessel decomposition of the cosmic shear field ("3D cosmic shear") is one way to maximise the amount of redshift information in a lensing analysis and therefore provides a powerful tool to investigate in particular the growth of cosmic structure that is crucial for dark energy studies. However, the computation of simulated 3D cosmic shear covariance matrices presents numerical difficulties, due to the required integrations over highly oscillatory functions. We present and compare two numerical methods and relative implementations to perform these integrations. We then show how to generate 3D Gaussian random fields on the sky in spherical coordinates, starting from the 3D cosmic shear covariances. To validate our field-generation procedure, we calculate the Minkowski functionals associated with our random fields, compare them with the known expectation values for the Gaussian case and demonstrate parameter inference from Minkowski functionals from a cosmic shear survey. This is a first step towards producing fully 3D Minkowski functionals for a lognormal field in 3D to extract Gaussian and non-Gaussian information from the cosmic shear field, as well as towards the use of Minkowski functionals as a probe of cosmology beyond the commonly used two-point statistics.

The alignment of galaxy spin with the shear field in observations

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

Hierarchical cosmic shear power spectrum inference

We develop a Bayesian hierarchical modelling approach for cosmic shear power spectrum inference, jointly sampling from the posterior distribution of the cosmic shear field and its (tomographic) power spectra. Inference of the shear power spectrum is a powerful intermediate product for a cosmic shear analysis, since it requires very few model assumptions and can be used to perform inference on a wide range of cosmological models \emph{a posteriori} without loss of information. We show that joint posterior for the shear map and power spectrum can be sampled effectively by Gibbs sampling, iteratively drawing samples from the map and power spectrum, each conditional on the other. This approach neatly circumvents difficulties associated with complicated survey geometry and masks that plague frequentist power spectrum estimators, since the power spectrum inference provides prior information about the field in masked regions at every sampling step. We demonstrate this approach for inference of tomographic shear $E$-mode, $B$-mode and $EB$-cross power spectra from a simulated galaxy shear catalogue with a number of important features; galaxies distributed on the sky and in redshift with photometric redshift uncertainties, realistic random ellipticity noise for every galaxy and a complicated survey mask. The obtained posterior distributions for the tomographic power spectrum coefficients recover the underlying simulated power spectra for both $E$- and $B$-modes.

The impact of supersurvey modes on cosmological constraints from cosmic shear fields

Owing to the mass-sheet degeneracy, cosmic shear maps do not probe directly the Fourier modes of the underlying mass distribution on scales comparable to the survey size and larger. To assess the corresponding effect on attainable cosmological parameter constraints, we quantify the information on super-survey modes in a lognormal model and, when interpreted as nuisance parameters, their degeneracies to cosmological parameters. Our analytical and numerical calculations clarify the central role of super-sample covariance (SSC) in shaping the statistical power of cosmological observables. Reconstructing the background modes from their non-Gaussian statistical dependence to small scales modes yields the renormalized convergence. This diagonalizes the spectrum covariance matrix, and the information content of the corresponding power spectrum is increased by a factor of two over standard methods. Unfortunately, careful calculation of the Cramer-Rao bound shows that the information recovery can never be made complete, any observable built from shear fields, including optimal sufficient statistics, are subject to severe information loss, typically $80\%$ to $90\%$ below $\ell \sim 3000$ for generic cosmological parameters. The lost information can only be recovered from additional, non-shear based data. Our predictions hold just as well for a tomographic analysis, and/or full sky surveys.

Investigating clustering dark energy with 3D weak cosmic shear

As observational evidence increasingly consolidates the case for a cosmological constant being the source of the Universe's accelerated expansion, the question whether, and if so, how well, future experiments could detect deviations from this standard scenario is raised with urgency. Assuming a dark energy component different from a cosmological constant, the observable effects in general include gravitational clustering described by the fluid's (rest-frame) speed of sound. We employ 3d weak cosmic shear, a proposed method to take advantage of the full three-dimensional information inherent to the cosmic shear field, to explore the capability of future surveys to detect dark energy clustering and the signature of an enhanced amplitude of the matter power spectrum on large scales. For this purpose, we present adequate numerical methods facilitating 3d weak cosmic shear calculations. We find that the possible constraints heavily depend on the dark energy equation of state w. If w is not very close to -1, constraining the squared sound speed within an order of magnitude seems possible with a combination of Euclid and Planck data.

Magnification as a tool in weak lensing

Weak lensing surveys exploit measurements of galaxy ellipticities. These measurements are subject to errors, which degrade the cosmological information that can be extracted from the surveys. Here, we propose a way of using the galaxy data themselves to calibrate the measurement errors. In particular, the cosmic shear field, which causes the galaxies to appear elliptical, also changes their sizes and fluxes. Information about the sizes and fluxes of the galaxies can be added to the shape information to obtain more robust information about the cosmic shear field. The net result will be tighter constraints on cosmological parameters such as those which describe dark energy.

3D photometric cosmic shear

Here we present a number of improvements to weak lensing 3D power spectrum analysis, 3D cosmic shear, that uses the shape and redshift information of every galaxy to constrain cosmological parameters. We show how photometric redshift probability distributions for individual galaxies can be directly included in this statistic with no averaging. We also include the Limber approximation, considerably simplifying full 3D cosmic shear analysis, and we investigate its range of applicability. Finally we show the relationship between weak lensing tomography and the 3D cosmic shear field itself; the steps connecting them being the Limber approximation, a harmonic-space transform and a discretisation in wavenumber. Each method has its advantages: 3D cosmic shear analysis allows straightforward inclusion of all relevant modes, thus ensuring minimum error bars, and direct control of the range of physical wavenumbers probed, to avoid the uncertain highly nonlinear regime. On the other hand, tomography is more convenient for checking systematics through direct investigation of the redshift dependence of the signal. Finally, for tomography, we suggest that the angular modes probed should be redshift-dependent, to recover some of the 3D advantages.

Weak lensing from space: first cosmological constraints from three-point shear statistics★

We use weak lensing data from the Hubble Space Telescope COSMOS survey to measure the second- and third-moments of the cosmic shear field, estimated from about 450,000 galaxies with average redshift <z> ~ 1.3. We measure two- and three-point shear statistics using a tree-code, dividing the signal in E, B and mixed components. We present a detection of the third-order moment of the aperture mass statistic and verify that the measurement is robust against systematic errors caused by point spread function (PSF) residuals and by the intrinsic alignments between galaxies. The amplitude of the measured three-point cosmic shear signal is in very good agreement with the predictions for a WMAP7 best-fit model, whereas the amplitudes of potential systematics are consistent with zero. We make use of three sets of large Lambda CDM simulations to test the accuracy of the cosmological predictions and to estimate the influence of the cosmology-dependent covariance. We perform a likelihood analysis using the measurement and find that the Omega_m-sigma_8 degeneracy direction is well fitted by the relation: sigma_8 (Omega_m/0.30)^(0.49)=0.78+0.11/-0.26. We present the first measurement of a more generalised three-point shear statistic and find a very good agreement with the WMAP7 best-fit cosmology. The cosmological interpretation of this measurement gives sigma_8 (Omega_m/0.30)^(0.46)=0.69 +0.08/-0.14. Furthermore, the combined likelihood analysis of this measurement with the measurement of the second order moment of the aperture mass improves the accuracy of the cosmological constraints, showing the high potential of this combination of measurements to infer cosmological constraints.

COSEBIs: Extracting the full E-/B-mode information from cosmic shear correlation functions

Cosmic shear is considered one of the most powerful methods for studying the properties of Dark Energy in the Universe. As a standard method, the two-point correlation functions $xi_\pm(theta)$ of the cosmic shear field are used as statistical measures for the shear field. In order to separate the observed shear into E- and B-modes, the latter being most likely produced by remaining systematics in the data set and/or intrinsic alignment effects, several statistics have been defined before. Here we aim at a complete E-/B-mode decomposition of the cosmic shear information contained in the $xi_\pm$ on a finite angular interval. We construct two sets of such E-/B-mode measures, namely Complete Orthogonal Sets of E-/B-mode Integrals (COSEBIs), characterized by weight functions between the $xi_\pm$ and the COSEBIs which are polynomials in $theta$ or polynomials in $ln(theta)$, respectively. Considering the likelihood in cosmological parameter space, constructed from the COSEBIs, we study their information contents. We show that the information grows with the number of COSEBI modes taken into account, and that an asymptotic limit is reached which defines the maximum available information in the E-mode component of the $xi_\pm$. We show that this limit is reached the earlier (i.e., for a smaller number of modes considered) the narrower the angular range is over which $xi_\pm$ are measured, and it is reached much earlier for logarithmic weight functions. For example, for $xi_\pm$ on the interval $1'\le \theta\le 400'$, the asymptotic limit for the parameter pair $(Omega_m, sigma_8)$ is reached for $\sim 25$ modes in the linear case, but already for 5 modes in the logarithmic case. The COSEBIs form a natural discrete set of quantities, which we suggest as method of choice in future cosmic shear likelihood analyses.

A new cosmic shear function: optimized E-/B-mode decomposition on a finite interval

The decomposition of the cosmic shear field into E- and B-mode is an important diagnostic in weak gravitational lensing. However, commonly used techniques to perform this separation suffer from mode-mixing on very small or very large scales. We introduce a new E-/B-mode decomposition of the cosmic shear two-point correlation on a finite interval. This new statistic is optimised for cosmological applications, by maximising the signal-to-noise ratio (S/N) and a figure of merit (FoM) based on the Fisher matrix of the cosmological parameters Omega_m and sigma_8. We improve both S/N and FoM results substantially with respect to the recently introduced ring statistic, which also provides E-/B-mode separation on a finite angular range. The S/N (FoM) is larger by a factor of three (two) on angular scales between 1 and 220 arc minutes. In addition, it yields better results than for the aperture-mass dispersion <M_ap>^2, with improvements of 20% (10%) for S/N (FoM). Our results depend on the survey parameters, most importantly on the covariance of the two-point shear correlation function. Although we assume parameters according to the CFHTLS-Wide survey, our method and optimisation scheme can be applied easily to any given survey settings and observing parameters. Arbitrary quantities, with respect to which the E-/B-mode filter is optimised, can be defined, therefore generalising the aim and context of the new shear statistic.

Cosmological model discrimination with weak lensing

Weak gravitational lensing provides a unique way of mapping directly the dark matter in the Universe. The majority of lensing analyses use the two-point statistics of the cosmic shear field to constrain the cosmological model, a method that is affected by degeneracies, such as that between σ8 and Ωm which are respectively the rms of the mass fluctuations on a scale of 8 Mpc/h and the matter density parameter, both at z = 0. However, the two-point statistics only measure the Gaussian properties of the field, and the weak lensing field is non-Gaussian. It has been shown that the estimation of non-Gaussian statistics for weak lensing data can improve the constraints on cosmological parameters. In this paper, we systematically compare a wide range of non-Gaussian estimators to determine which one provides tighter constraints on the cosmological parameters. These statistical methods include skewness, kurtosis, and the higher criticism test, in several sparse representations such as wavelet and curvelet; as well as the bispectrum, peak counting, and a newly introduced statistic called wavelet peak counting (WPC). Comparisons based on sparse representations indicate that the wavelet transform is the most sensitive to non-Gaussian cosmological structures. It also appears that the most helpful statistic for non-Gaussian characterization in weak lensing mass maps is the WPC. Finally, we show that the σ8–Ωm degeneracy could be even better broken if the WPC estimation is performed on weak lensing mass maps filtered by the wavelet method, MRLens.

FAst STatistics for weak Lensing (FASTLens): fast method for weak lensing statistics and map making

With increasingly large data sets, weak lensing measurements are able to measure cosmological parameters with ever greater precision. However this increased accuracy also places greater demands on the statistical tools used to extract the available information. To date, the majority of lensing analyses use the two point-statistics of the cosmic shear field. These can either be studied directly using the two-point correlation function, or in Fourier space, using the power spectrum. But analyzing weak lensing data inevitably involves the masking out of regions or example to remove bright stars from the field. Masking out the stars is common practice but the gaps in the data need proper handling. In this paper, we show how an inpainting technique allows us to properly fill in these gaps with only $N \log N$ operations, leading to a new image from which we can compute straight forwardly and with a very good accuracy both the pow er spectrum and the bispectrum. We propose then a new method to compute the bispectrum with a polar FFT algorithm, which has the main advantage of avoiding any interpolation in the Fourier domain. Finally we propose a new method for dark matter mass map reconstruction from shear observations which integrates this new inpainting concept. A range of examples based on 3D N-body simulations illustrates the results.

Galaxy-galaxy-galaxy lensing: Third-order correlations between the galaxy and mass distributions in the Universe

Galaxy-galaxy lensing (GGL) measures the 2-point cross-correlation between galaxies and mass in the Universe. In this work we seek to generalise this effect by considering the third-order correlations between galaxies and mass: galaxy-galaxy- galaxy lensing. Third-order correlations in the cosmic shear field have recently been reported in the VIRMOS-DESCART and CTIO surveys. Such data should also be ideal for measuring galaxy-galaxy-galaxy lensing. Indeed, the effects of these higher- order correlations may have already been detected in recent studies of galaxy-galaxy lensing. Higher-order cross-correlation functions contain invaluable information about the relationship between galaxies and their mass environments that GGL studies alone cannot detect. In this paper we lay out the basic relations for third-order cross correlations and their projections and introduce a new set of scale dependent third-order bias parameters. We define three new observables: two galaxy-shear-shear correlation functions, G±, and a galaxy-galaxy-shear correlation, G. We relate these to the various projected cross-bispectra and give practical estimators for their measurement. We note that the observational signature of these correlators is simply the excess shear-shear correlation measured about foreground galaxies (for G±) and the average tangential shear around foreground galaxy pairs (for G). These quantities are no more than second order in the shear and so should be more easily measurable than the shear 3-point correlation. Finally we derive expressions for the third order aperture mass statistics in terms of both the cross-bispectra and the real-space correlation functions. Such statistics provide a very localized measurement of the bispectra, thus encapsulating essentially all of the available third-order information, while remaining easily obtainable from observations of 3-point cross-correlation functions. In addition we find that utilising aperture statistics has the further benefit that they measure only the connected part of the third order correlation.

The consequences of parity symmetry for higher-order statistics of cosmic shear and other polar fields

We consider the parity transformation and the consequences of parity invariance on the n-point correlation function of the cosmic shear caused by gravitational lensing by the large-scale structure, or any other polar (or spin-2) field. The decomposition of the shear field into E -a ndB-modes then yields the result that any correlation function which contains an odd number of B-mode shear components vanishes for parity-invariant random shear fields. In particular, this result implies that the expectation value of the third-order cross aperture statistics, D M 3 E , vanishes for parity-invariant shear fields. Therefore, a significant detection of a non-zero value of D M 3 E in a cosmic shear survey does not indicate the presence of B-modes, but of an underestimate of the statistical uncertainty or cosmic variance, or a remaining systematic eect in the shear measurement. We argue that the parity invariance provides a very specific diagnostic of systematic eects in shear data. Our results apply as well to the linear polarization of the cosmic microwave background.