We investigate the performance of the Metacalibration shear calibration framework using simulated imaging data for the Nancy Grace Roman Space Telescope (Roman) reference High-Latitude Imaging Survey ...(HLIS). The weak lensing program of the Roman mission requires the mean weak lensing shear estimate to be calibrated within about 0.03%. To reach this goal, we can test our calibration process with various simulations and ultimately isolate the sources of residual shear biases in order to improve our methods. In this work, we build on the Roman HLIS image simulation pipeline in Troxel et al. 2021 to incorporate several new realistic processing-pipeline updates necessary to more accurately process the imaging data and calibrate the shear. We show the first results of this calibration for six deg\(^2\) of the simulated reference HLIS using Metacalibration and compare these results to measurements on more simple, faster Roman-like image simulations. In both cases, we neglect the impact of blending of objects. We find that in the simplified simulations, Metacalibration can calibrate shapes to be within \(m=(-0.01\pm 0.10)\)%. When applied to the current most-realistic version of the simulations, the precision is much lower, with estimates of \(m=(-1.34\pm 0.67)\)% for joint multi-band single-epoch measurements and \(m=(-1.13\pm 0.60)\)% for multi-band coadd measurements. These results are all consistent with zero within 1-2\(\sigma\), indicating we are currently limited by our simulated survey volume. Further work on testing the shear calibration methodology is necessary at higher precision to reach the level of the Roman requirements, in particular in the presence of blending. Current results demonstrate, however, that the Metacalibration method can work on undersampled space-based Roman imaging data at levels comparable to the requirements of current weak lensing surveys.
PSFs of coadded images Mandelbaum, Rachel; Jarvis, Mike; Lupton, Robert H ...
arXiv (Cornell University),
02/2023
Paper, Journal Article
Odprti dostop
We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what ...properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows: 1. We provide an illustration of how linear coaddition schemes can produce a coadd that lacks a well-defined PSF even for relatively simple scenarios and choices of weight functions. 2. We provide a more formal demonstration of the fact that a linear coadd only has a well-defined PSF in the case that either (a) each input image has the same PSF or (b) the coadd is produced with weights that are independent of the signal. 3. We discuss some reasons that two plausible nonlinear coaddition algorithms (median and clipped-mean) fail to produce a consistent PSF profile for stars. 4. We demonstrate that all nonlinear coaddition procedures fail to produce a well-defined PSF for extended objects. In the end, we conclude that, for any purpose where a well-defined PSF is desired, one should use a linear coaddition scheme with weights that do not correlate with the signal and are approximately uniform across typical objects of interest.
Weak gravitational lensing, or weak lensing, is one of the most powerful probes for dark matter and dark energy science, although it faces increasing challenges in controlling systematic ...uncertainties as \edit{the statistical errors become smaller}. The Point Spread Function (PSF) needs to be precisely modeled to avoid systematic error on the weak lensing measurements. The weak lensing biases induced by errors in the PSF model second moments, i.e., its size and shape, are well-studied. However, Zhang et al. (2021) showed that errors in the higher moments of the PSF may also be a significant source of systematics for upcoming weak lensing surveys. Therefore, the goal of this work is to comprehensively investigate the modeling quality of PSF moments from the \(3^{\text{rd}}\) to \(6^{\text{th}}\) order, and estimate their impact on cosmological parameter inference. We propagate the \textsc{PSFEx} higher moments modeling error in the HSC survey dataset to the weak lensing \edit{shear-shear correlation functions} and their cosmological analyses. We find that the overall multiplicative shear bias associated with errors in PSF higher moments can cause a \(\sim 0.1 \sigma\) shift on the cosmological parameters for LSST Y10. PSF higher moment errors also cause additive biases in the weak lensing shear, which, if not accounted for in the cosmological parameter analysis, can induce cosmological parameter biases comparable to their \(1\sigma\) uncertainties for LSST Y10. We compare the \textsc{PSFEx} model with PSF in Full FOV (\textsc{Piff}), and find similar performance in modeling the PSF higher moments. We conclude that PSF higher moment errors of the future PSF models should be reduced from those in current methods to avoid a need to explicitly model these effects in the weak lensing analysis.
The High Latitude Survey of the Nancy Grace Roman Space Telescope is expected to measure the positions and shapes of hundreds of millions of galaxies in an area of 2220 deg\(^2\). This survey will ...provide high-quality weak lensing data with unprecedented systematics control. The Roman Space Telescope will survey the sky in near infrared (NIR) bands using Teledyne H4RG HgCdTe photodiode arrays. These NIR arrays exhibit an effect called persistence: charges that are trapped in the photodiodes during earlier exposures are gradually released into later exposures, leading to contamination of the images and potentially to errors in measured galaxy properties such as fluxes and shapes. In this work, we use image simulations that incorporate the persistence effect to study its impact on galaxy shape measurements and weak lensing signals. No significant spatial correlations are found between the galaxy shape changes induced by persistence. On the scales of interest for weak lensing cosmology, the effect of persistence on the weak lensing correlation function is about two orders of magnitude lower than the Roman Space Telescope additive shear error budget, indicating that the persistence effect is expected to be a subdominant contributor to the systematic error budget for weak lensing with the Roman Space Telescope given its current design.
The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. ...Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired "target" point spread function (PSF). We present here an initial application of Imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that Imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the "fidelity"), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final fidelity and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ("Paper II") explores the implications for weak lensing analyses.
One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial ...application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their 1- and 2-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1 - 6) x 10-4 per component depending on the band; the correlation functions are >= 2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman 2nd moment error budget, although the 4th moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all 4 filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman.
The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing ...astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.