Weak Lensing Surveys Research Articles

Radio weak lensing shear measurement in the visibility domain – II. Source extraction

This paper extends the method introduced in Rivi et al. (2016b) to measure galaxy ellipticities in the visibility domain for radio weak lensing surveys. In that paper we focused on the development and testing of the method for the simple case of individual galaxies located at the phase centre, and proposed to extend it to the realistic case of many sources in the field of view by isolating visibilities of each source with a faceting technique. In this second paper we present a detailed algorithm for source extraction in the visibility domain and show its effectiveness as a function of the source number density by running simulations of SKA1-MID observations in the band 950-1150 MHz and comparing original and measured values of galaxies' ellipticities. Shear measurements from a realistic population of 10^4 galaxies randomly located in a field of view of 1 deg^2 (i.e. the source density expected for the current radio weak lensing survey proposal with SKA1) are also performed. At SNR >= 10, the multiplicative bias is only a factor 1.5 worse than what found when analysing individual sources, and is still comparable to the bias values reported for similar measurement methods at optical wavelengths. The additive bias is unchanged from the case of individual sources, but is significantly larger than typically found in optical surveys. This bias depends on the shape of the uv coverage and we suggest that a uv-plane weighting scheme to produce a more isotropic shape could reduce and control additive bias.

Probing supervoids with weak lensing

The cosmic microwave background (CMB) has non-Gaussian features in the temperature fluctuations. An anomalous cold spot surrounded with a hot ring, called the Cold Spot is one of such features. If a large underdence region (supervoid) resides towards the Cold Spot, we would be able to detect a systematic shape distortion in the images of background source galaxies via weak lensing effect. In order to estimate the detectability of such signals, we used the data of $N$-body simulations to simulate full-sky ray-tracing of source galaxies. We searched for a most prominent underdense region using the simulated convergence maps smoothed at a scale of 20 degree and obtained tangential shears around it. The lensing signal expected in a concordant $\Lambda$CDM model can be detected at a signal-to-noise ratio $S/N\sim3$. If a supervoid with a radius of $\sim 200\,h^{-1}\,\textrm{Mpc}$ and a density contrast $\delta_0 \sim -0.3$ at the centre resides at a redshift $z\sim 0.2$, on-going and near-future weak gravitational lensing surveys would detect a lensing signal with $S/N\gtrsim4$ without resorting to stacking. From the tangential shear profile, we can obtain a constraint on the projected mass distribution of the supervoid.

Increasing the lensing figure of merit through higher order convergence moments

The unprecedented quality, the increased dataset, and the wide area of ongoing and near future weak lensing surveys allows to move beyond the standard two points statistics thus making worthwhile to investigate higher order probes. As an interesting step towards this direction, we expolore the use of higher order moments (HOM) of the convergence field as a way to increase the lensing Figure of Merit (FoM). To this end, we rely on simulated convergence to first show that HOM can be measured and calibrated so that it is indeed possible to predict them for a given cosmological model provided suitable nuisance parameters are introduced and then marginalized over. We then forecast the accuracy on cosmological parameters from the use of HOM alone and in combination with standard shear power spectra tomography. It turns out that HOM allow to break some common degeneracies thus significantly boosting the overall FoM. We also qualitatively discuss possible systematics and how they can be dealt with.

Full-sky Ray-tracing Simulation of Weak Lensing Using ELUCID Simulations: Exploring Galaxy Intrinsic Alignment and Cosmic Shear Correlations

Abstract The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mock weak-lensing surveys. In this work, we use the full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID N-body simulation run with WMAP9 cosmology. In our model, we assume that the shape of the central elliptical galaxy follows that of the dark matter halo, and that of the spiral galaxy follows the halo spin. Using the mock galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of the Kilo-Degree Survey (KiDS) and Deep Lens Survey (DLS). We find that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy; otherwise, the shear correlations on small scales are too high. Most importantly, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of the elliptical galaxy. Our results explain the origin of the detected positive GI term in the weak-lensing surveys. We conclude that in future analyses, the GI model must include the dependence on galaxy types in more detail.

Constraining cosmology with the velocity function of low-mass galaxies

The number density of field galaxies per rotation velocity, referred to as the velocity function, is an intriguing statistical measure probing the smallest scales of structure formation. In this paper we point out that the velocity function is sensitive to small shifts in key cosmological parameters such as the amplitude of primordial perturbations ($\sigma_8$) or the total matter density ($\Omega_{\rm m}$). Using current data and applying conservative assumptions about baryonic effects, we show that the observed velocity function of the Local Volume favours cosmologies in tension with the measurements from Planck but in agreement with the latest findings from weak lensing surveys. While the current systematics regarding the relation between observed and true rotation velocities are potentially important, upcoming data from HI surveys as well as new insights from hydrodynamical simulations will dramatically improve the situation in the near future.

Cosmic Shear with Einstein Rings

Abstract We explore a new technique to measure cosmic shear using Einstein rings. In Birrer et al., we showed that the detailed modeling of Einstein rings can be used to measure external shear to high precision. In this Letter, we explore how a collection of Einstein rings can be used as a statistical probe of cosmic shear. We present a forecast of the cosmic shear information available in Einstein rings for different strong lensing survey configurations. We find that, assuming that the number density of Einstein rings in the COSMOS survey is representative, future strong lensing surveys should have a cosmological precision comparable to the current ground-based weak lensing surveys. We discuss how this technique is complementary to the standard cosmic shear analyses since it is sensitive to different systematics and can be used for cross-calibration.

Identification of SIPS J2045-6332 as a Partially Resolved Binary

We show that SIPS J2045-6332, a late M/early L object previously identified as a candidate spectral mix binary, shows an elongated image shape. Using shape measurement techniques originally developed for cosmological weak lensing surveys on VISTA VHS images we show that this likely blended binary has an implied position angle of ~290 degrees East of North with a secondary companion that is likely to be a late L dwarf. This object would be a good follow-up target for high resolution imaging studies

Weak Gravitational Lensing of Finite Beams.

The standard theory of weak gravitational lensing relies on the infinitesimal light beam approximation. In this context, images are distorted by convergence and shear, the respective sources of which unphysically depend on the resolution of the distribution of matter-the so-called Ricci-Weyl problem. In this Letter, we propose a strong-lensing-inspired formalism to describe the lensing of finite beams. We address the Ricci-Weyl problem by showing explicitly that convergence is caused by the matter enclosed by the beam, regardless of its distribution. Furthermore, shear turns out to be systematically enhanced by the finiteness of the beam. This implies, in particular, that the Kaiser-Squires relation between shear and convergence is violated, which could have profound consequences on the interpretation of weak-lensing surveys.

Partially acoustic dark matter cosmology and cosmological constraints

Observations of the cosmic microwave background (CMB) together with weak lensing measurements of the clustering of large scale cosmological structures and local measurements of the Hubble constant pose a challenge to the standard $\Lambda$CDM cosmological model. On one side CMB observations imply a Hubble constant that is lower than local measurements and an amplitude of the lensing signal that is higher than direct measurements from weak lensing surveys. We investigate a way of relieving these tensions by adding dark radiation tightly coupled to an acoustic part of the dark matter sector and compare it to massive neutrino solutions. While these models offer a way of separately relieving the Hubble and weak lensing tensions they are prevented from fully accommodating both at the same time since the CMB requires additional cold dark matter when adding acoustic dark matter or massive neutrinos to preserve the same sharpness of the acoustic peaks which counteracts the desired growth suppression.

Non-negative Matrix Factorization for Self-calibration of Photometric Redshift Scatter in Weak-lensing Surveys

Abstract Photo-z error is one of the major sources of systematics degrading the accuracy of weak-lensing cosmological inferences. Zhang et al. proposed a self-calibration method combining galaxy–galaxy correlations and galaxy–shear correlations between different photo-z bins. Fisher matrix analysis shows that it can determine the rate of photo-z outliers at a level of 0.01%–1% merely using photometric data and do not rely on any prior knowledge. In this paper, we develop a new algorithm to implement this method by solving a constrained nonlinear optimization problem arising in the self-calibration process. Based on the techniques of fixed-point iteration and non-negative matrix factorization, the proposed algorithm can efficiently and robustly reconstruct the scattering probabilities between the true-z and photo-z bins. The algorithm has been tested extensively by applying it to mock data from simulated stage IV weak-lensing projects. We find that the algorithm provides a successful recovery of the scatter rates at the level of 0.01%–1%, and the true mean redshifts of photo-z bins at the level of 0.001, which may satisfy the requirements in future lensing surveys.

Linear and non-linear Modified Gravity forecasts with future surveys

Abstract Modified Gravity theories generally affect the Poisson equation and the gravitational slip in an observable way, that can be parameterized by two generic functions ( η and μ ) of time and space. We bin their time dependence in redshift and present forecasts on each bin for future surveys like Euclid. We consider both Galaxy Clustering and Weak Lensing surveys, showing the impact of the non-linear regime, with two different semi-analytical approximations. In addition to these future observables, we use a prior covariance matrix derived from the Planck observations of the Cosmic Microwave Background. In this work we neglect the information from the cross correlation of these observables, and treat them as independent. Our results show that η and μ in different redshift bins are significantly correlated, but including non-linear scales reduces or even eliminates the correlation, breaking the degeneracy between Modified Gravity parameters and the overall amplitude of the matter power spectrum. We further apply a Zero-phase Component Analysis and identify which combinations of the Modified Gravity parameter amplitudes, in different redshift bins, are best constrained by future surveys. We extend the analysis to two particular parameterizations of μ and η and consider, in addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future surveys will be able to constrain the current values of η and μ at the 2–5% level when using only linear scales (wavevector k 0 . 15 h/Mpc), depending on the specific time parameterization; sensitivity improves to about 1 % when non-linearities are included.

Constraints on the dark matter and dark energy interactions from weak lensing bispectrum tomography

We estimate uncertainties of cosmological parameters for phenomenological interacting dark energy models using weak lensing convergence power spectrum and bispectrum. We focus on the bispectrum tomography and examine how well the weak lensing bispectrum with tomography can constrain the interactions between dark sectors, as well as other cosmological parameters. Employing the Fisher matrix analysis, we forecast parameter uncertainties derived from weak lensing bispectra with a two-bin tomography and place upper bounds on strength of the interactions between the dark sectors. The cosmic shear will be measured from upcoming weak lensing surveys with high sensitivity, thus it enables us to use the higher order correlation functions of weak lensing to constrain the interaction between dark sectors and will potentially provide more stringent results with other observations combined.

On Using a Space Telescope to Detect Weak-lensing Shear

Ignoring redshift dependence, the statistical performance of a weak-lensing survey is set by two numbers: the effective shape noise of the sources, which includes the intrinsic ellipticity dispersion and the measurement noise, and the density of sources that are useful for weak-lensing measurements. In this paper, we provide some general guidance for weak-lensing shear measurements from a “generic” space telescope by looking for the optimum wavelength bands to maximize the galaxy flux signal-to-noise ratio (S/N) and minimize ellipticity measurement error. We also calculate an effective galaxy number per square degree across different wavelength bands, taking into account the density of sources that are useful for weak-lensing measurements and the effective shape noise of sources. Galaxy data collected from the ultra-deep UltraVISTA Ks-selected and R-selected photometric catalogs (Muzzin et al. ) are fitted to radially symmetric Sérsic galaxy light profiles. The Sérsic galaxy profiles are then stretched to impose an artificial weak-lensing shear, and then convolved with a pure Airy Disk PSF to simulate imaging of weak gravitationally lensed galaxies from a hypothetical diffraction-limited space telescope. For our model calculations and sets of galaxies, our results show that the peak in the average galaxy flux S/N, the minimum average ellipticity measurement error, and the highest effective galaxy number counts all lie around the K-band near 2.2 μm.

Massive data compression for parameter-dependent covariance matrices

Abstract We show how the massive data compression algorithm MOPED can be used to reduce, by orders of magnitude, the number of simulated data sets which are required to estimate the covariance matrix required for the analysis of Gaussian-distributed data. This is relevant when the covariance matrix cannot be calculated directly. The compression is especially valuable when the covariance matrix varies with the model parameters. In this case, it may be prohibitively expensive to run enough simulations to estimate the full covariance matrix throughout the parameter space. This compression may be particularly valuable for the next generation of weak lensing surveys, such as proposed for Euclid and Large Synoptic Survey Telescope, for which the number of summary data (such as band power or shear correlation estimates) is very large, ∼104, due to the large number of tomographic redshift bins which the data will be divided into. In the pessimistic case where the covariance matrix is estimated separately for all points in an Monte Carlo Markov Chain analysis, this may require an unfeasible 109 simulations. We show here that MOPED can reduce this number by a factor of 1000, or a factor of ∼106 if some regularity in the covariance matrix is assumed, reducing the number of simulations required to a manageable 103, making an otherwise intractable analysis feasible.

Multipolar moments of weak lensing signal around clusters

Context. Upcoming weak lensing surveys such as Euclid will provide an unprecedented opportunity to quantify the geometry and topology of the cosmic web, in particular in the vicinity of lensing clusters. Aims. Understanding the connectivity of the cosmic web with unbiased mass tracers, such as weak lensing, is of prime importance to probe the underlying cosmology, seek dynamical signatures of dark matter, and quantify environmental effects on galaxy formation. Methods. Mock catalogues of galaxy clusters are extracted from the N-body PLUS simulation. For each cluster, the aperture multipolar moments of the convergence are calculated in two annuli (inside and outside the virial radius). By stacking their modulus, a statistical estimator is built to characterise the angular mass distribution around clusters. The moments are compared to predictions from perturbation theory and spherical collapse. Results. The main weakly chromatic excess of multipolar power on large scales is understood as arising from the contraction of the primordial cosmic web driven by the growing potential well of the cluster. Besides this boost, the quadrupole prevails in the cluster (ellipsoidal) core, while at the outskirts, harmonic distortions are spread on small angular modes, and trace the non-linear sharpening of the filamentary structures. Predictions for the signal amplitude as a function of the cluster-centric distance, mass, and redshift are presented. The prospects of measuring this signal are estimated for current and future lensing data sets. Conclusions. The Euclid mission should provide all the necessary information for studying the cosmic evolution of the connectivity of the cosmic web around lensing clusters using multipolar moments and probing unique signatures of, for example, baryons and warm dark matter.

Two- and three-dimensional wide-field weak lensing mass maps from the Hyper Suprime-Cam Subaru Strategic Program S16A data

Abstract We present wide-field (167 deg2) weak lensing mass maps from the Hyper Supreme-Cam Subaru Strategic Program (HSC-SSP). We compare these weak lensing based dark matter maps with maps of the distribution of the stellar mass associated with luminous red galaxies. We find a strong correlation between these two maps with a correlation coefficient of ρ = 0.54 ± 0.03 (for a smoothing size of 8΄). This correlation is detected even with a smaller smoothing scale of 2΄ (ρ = 0.34 ± 0.01). This detection is made uniquely possible because of the high source density of the HSC-SSP weak lensing survey ($\bar{n}\sim 25\:$arcmin−2). We also present a variety of tests to demonstrate that our maps are not significantly affected by systematic effects. By using the photometric redshift information associated with source galaxies, we reconstruct a three-dimensional mass map. This three-dimensional mass map is also found to correlate with the three-dimensional galaxy mass map. Cross-correlation tests presented in this paper demonstrate that the HSC-SSP weak lensing mass maps are ready for further science analyses.

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.

Estimating the weak-lensing rotation signal in radio cosmic shear surveys

Weak lensing has become an increasingly important tool in cosmology and the use of galaxy shapes to measure cosmic shear has become routine. The weak-lensing distortion tensor contains two other effects in addition to the two components of shear: the convergence and rotation. The rotation mode is not measurable using the standard cosmic shear estimators based on galaxy shapes, as there is no information on the original shapes of the images before they were lensed. Due to this, no estimator has been proposed for the rotation mode in cosmological weak-lensing surveys, and the rotation mode has never been constrained. Here, we derive an estimator for this quantity, which is based on the use of radio polarisation measurements of the intrinsic position angles of galaxies. The rotation mode can be sourced by physics beyond $\Lambda$CDM, and also offers the chance to perform consistency checks of $\Lambda$CDM and of weak-lensing surveys themselves. We present simulations of this estimator and show that, for the pedagogical example of cosmic string spectra, this estimator could detect a signal that is consistent with the constraints from Planck. We examine the connection between the rotation mode and the shear $B$-modes and thus how this estimator could help control systematics in future radio weak-lensing surveys.

Looking through the same lens: Shear calibration for LSST, Euclid, and WFIRST with stage 4 CMB lensing

The next generation weak lensing surveys (i.e., LSST, Euclid and WFIRST) will require exquisite control over systematic effects. In this paper, we address shear calibration and present the most realistic forecast to date for LSST/Euclid/WFIRST and CMB lensing from a stage 4 CMB experiment (CMB S4). We use the CosmoLike code to simulate a joint analysis of all the two-point functions of galaxy density, galaxy shear and CMB lensing convergence. We include the full Gaussian and non-Gaussian covariances and explore the resulting joint likelihood with Monte Carlo Markov Chains. We constrain shear calibration biases while simultaneously varying cosmological parameters, galaxy biases and photometric redshift uncertainties. We find that CMB lensing from CMB S4 enables the calibration of the shear biases down to 0.2% - 3% in 10 tomographic bins for LSST (below the ~0.5% requirements in most tomographic bins), down to 0.4% - 2.4% in 10 bins for Euclid and 0.6% - 3.2% in 10 bins for WFIRST. For a given lensing survey, the method works best at high redshift where shear calibration is otherwise most challenging. This self-calibration is robust to Gaussian photometric redshift uncertainties and to a reasonable level of intrinsic alignment. It is also robust to changes in the beam and the effectiveness of the component separation of the CMB experiment, and slowly dependent on its depth, making it possible with third generation CMB experiments such as AdvACT and SPT-3G, as well as the Simons Observatory.

The effect of Limber and flat-sky approximations on galaxy weak lensing

We review the effect of the commonly-used Limber and flat-sky approximations on the calculation of shear power spectra and correlation functions for galaxy weak lensing. These approximations are accurate at small scales, but it has been claimed recently that their impact on low multipoles could lead to an increase in the amplitude of the mass fluctuations inferred from surveys such as CFHTLenS, reducing the tension between galaxy weak lensing and the amplitude determined by Planck from observations of the cosmic microwave background. Here, we explore the impact of these approximations on cosmological parameters derived from weak lensing surveys, using the CFHTLenS data as a test case. We conclude that the use of small-angle approximations for cosmological parameter estimation is negligible for current data, and does not contribute to the tension between current weak lensing surveys and Planck.