Fourier Power Function Shapelets (FPFS) shear estimator: Performance on image simulations

Abstract

We reinterpret the shear estimator developed by Zhang & Komatsu within the framework of shapelets and propose the Fourier Power Function Shapelets (FPFS) shear estimator. Four shapelet modes are calculated from the power function of every galaxy's Fourier transform after deconvolving the point-spread function (PSF) in Fourier space. We propose a novel normalization scheme to construct dimensionless ellipticity and its corresponding shear responsivity using these shapelet modes. Shear is measured in a conventional way by averaging the ellipticities and responsivities over a large ensemble of galaxies.With the introduction and tuning of a weighting parameter, noise bias is reduced below one per cent of the shear signal. We also provide an iterative method to reduce selection bias. The FPFS estimator is developed without any assumptions regarding galaxy morphology or any approximations for PSF correction. Moreover, our method does not rely on heavy image manipulations or complicated statistical procedures.We test the FPFS shear estimator using several Subaru Hyper Suprime- Cam (HSC)-like image simulations and the main results are as follows. (i) For simulations that only contain isolated galaxies, the amplitude of the multiplicative bias is below 1 × 10-2. (ii) For more realistic simulations, which also contain blended galaxies, the blended galaxies are deblended by the first-generation HSC deblender before shear measurement. A multiplicative bias of (-5.71 ± 0.31) × 10-2 is found. The blending bias is calibrated by image simulations. Finally, we test the consistency and stability of this calibration.