Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of Weak Gravitational Lensing - II : Including Finite Support Effects in Cosmic Shear Covariance Matrices
405
0
0.0
(
0
)
Added by
Joachim Harnois-Deraps
Publication date
2014
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
(Abridged) We investigate and quantify the impact of finite simulation volume on weak lensing two- and four-point statistics. These {it finite support} (FS) effects are modelled for several estimators, simulation box sizes and source redshifts, and validated against a new large suite of 500 $N$-body simulations. The comparison reveals that our theoretical model is accurate to better than 5 per cent for the shear correlation function $xi_{+}(theta)$ and its error. We find that the most important quantities for FS modelling is the ratio between the measured angle $theta$ and the angular size of the simulation box at the source redshift, $theta_{box}(z_s)$, or the multipole equivalent $ell / ell_{box}(z_s)$. When this ratio reaches 0.1, independently of the source redshift, the shear correlation function $xi_+$ is suppressed by 5, 10, 20 and 25 percent for $L_{box}= 1000$, $500$, $250$ and $147mbox{Mpc}/h$ respectively. When it reaches 0.2, the suppression exceeds 25 percent even for the largest box. The same effect is observed in $xi_{-}(theta)$, but at much larger angles. This has important consequences for cosmological analyses using $N$-body simulations to calibrate the impact of non-linear gravitational clustering or to estimate errors and systematics effects, and should not be overlooked. We propose simple semi-analytic solutions to correct for these finite box effects with and without the presence of survey masks, and the method can be generalized to any weak lensing estimator. This offers a graceful solution to the important problem of estimating accurate covariance matrices for weak lensing studies: there is no need to run extra large simulation volumes, as long as the box effects are corrected.
rate research
Read More
Gravitational lensing surveys have now become large and precise enough that the interpretation of the lensing signal has to take into account an increasing number of theoretical limitations and observational biases. Since the lensing signal is the strongest at small angular scales, only numerical simulations can reproduce faithfully the non-linear dynamics and secondary effects at play. This work is the first of a series in which all gravitational lensing corrections known so far will be implemented in the same set of simulations, using realistic mock catalogues and non-Gaussian statistics. In this first paper, we present the TCS simulation suite and compute basic statistics such as the second and third order convergence and shear correlation functions. These simple tests set the range of validity of our simulations, which are resolving most of the signals at the sub-arc minute level (or $ell sim 10^4$). We also compute the non-Gaussian covariance matrix of several statistical estimators, including many that are used in the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). From the same realizations, we construct halo catalogues, computing a series of properties that are required by most galaxy population algorithms. These simulation products are publicly available for download.
Upcoming weak lensing surveys will probe large fractions of the sky with unprecedented accuracy. To infer cosmological constraints, a large ensemble of survey simulations are required to accurately model cosmological observables and their covariances. We develop a parallelized multi-lens-plane pipeline called UFalcon, designed to generate full-sky weak lensing maps from lightcones within a minimal runtime. It makes use of L-PICOLA, an approximate numerical code, which provides a fast and accurate alternative to cosmological $N$-Body simulations. The UFalcon maps are constructed by nesting 2 simulations covering a redshift-range from $z=0.1$ to $1.5$ without replicating the simulation volume. We compute the convergence and projected overdensity maps for L-PICOLA in the lightcone or snapshot mode. The generation of such a map, including the L-PICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps to calculate the spherical harmonic power spectra, which we compare to theoretical predictions and to UFalcon results generated using the full $N$-Body code GADGET-2. We then compute the covariance matrix of the full-sky spherical harmonic power spectra using 150 UFalcon maps based on L-PICOLA in lightcone mode. We consider the PDF, the higher-order moments and the variance of the smoothed field variance to quantify the accuracy of the covariance matrix, which we find to be a few percent for scales $ell sim 10^2$ to $10^3$. We test the impact of this level of accuracy on cosmological constraints using an optimistic survey configuration, and find that the final results are robust to this level of uncertainty. The speed and accuracy of our developed pipeline provides a basis to also include further important features such as masking, varying noise and will allow us to compute covariance matrices for models beyond $Lambda$CDM. [abridged]
We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample of galaxies, with well-calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $6.17$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal, and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1sigma$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = sigma_8sqrt{Omega_{rm m}/0.3}$. Through a series of two-point null-tests, we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. The PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is, however, shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test, we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses.
We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58 in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level.
The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is proposed to obtain deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). At present, the observations for the CDFS field have been completed, and comprise in total about 4.9 deg$^2$ down to $r_mathrm{AB}$$sim$26 mag. In the companion paper by Fu et al. (2018), we present the weak lensing shear measurements for $r$-band images with seeing $le$ 0.9 arcsec. In this paper, we perform image simulations to calibrate possible biases of the measured shear signals. Statistically, the properties of the simulated point spread function (PSF) and galaxies show good agreements with those of observations. The multiplicative bias is calibrated to reach an accuracy of $sim$3.0%. We study the bias sensitivities to the undetected faint galaxies and to the neighboring galaxies. We find that undetected galaxies contribute to the multiplicative bias at the level of $sim$0.3%. Further analysis shows that galaxies with lower signal-to-noise ratio (SNR) are impacted more significantly because the undetected galaxies skew the background noise distribution. For the neighboring galaxies, we find that although most have been rejected in the shape measurement procedure, about one third of them still remain in the final shear sample. They show a larger ellipticity dispersion and contribute to $sim$0.2% of the multiplicative bias. Such a bias can be removed by further eliminating these neighboring galaxies. But the effective number density of the galaxies can be reduced considerably. Therefore efficient methods should be developed for future weak lensing deep surveys.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
