HybPSF: Hybrid PSF reconstruction for the observed JWST NIRCam image

Lin Nie,Huanyuan Shan,Guoliang Li,Lei Wang,Charling Tao,Qifan Cui,Yushan Xie,Dezi Liu,Zekang Zhang

2023-08-31

Abstract:The James Webb Space Telescope (JWST) ushers in a new era of astronomical observation and discovery, offering unprecedented precision in a variety of measurements such as photometry, astrometry, morphology, and shear measurement. Accurate point spread function (PSF) models are crucial for many of these measurements. In this paper, we introduce a hybrid PSF construction method called HybPSF for JWST NIRCam imaging data. HybPSF combines the WebbPSF software, which simulates the PSF for JWST, with observed data to produce more accurate and reliable PSF models. We apply this method to the SMACS J0723 imaging data and construct supplementary structures from residuals obtained by subtracting the WebbPSF PSF model from the data. Our results show that HybPSF significantly reduces discrepancies between the PSF model and the data compared to WebbPSF. Specifically, the PSF shape parameter ellipticity and size comparisons indicate that HybPSF improves precision by a factor of approximately 10 for \$R^2\$ and \$50\%\$ for \$e\$. This improvement has important implications for astronomical measurements using JWST NIRCam imaging data.

Instrumentation and Methods for Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to improve the accuracy of the point - spread function (PSF) model in the images of the James Webb Space Telescope (JWST) Near - Infrared Camera (NIRCam). Specifically, the authors introduce a hybrid PSF reconstruction method (HybPSF), aiming to combine the PSF simulated by the WebbPSF software and the observational data to generate a more accurate and reliable PSF model.

### Background and Problem Description

1. **Importance of PSF**:
- An accurate PSF model is crucial for a variety of astronomical measurements, such as photometry, astrometry, morphology, and shear measurements.
- Inaccurate PSF modeling can lead to systematic errors, affecting the precision of studies such as the weak gravitational lensing effect and galaxy cluster mass reconstruction.

2. **Existing Challenges**:
- For space telescopes (such as JWST), star images are usually under - sampled, and the number of PSF stars in a single exposure may not be sufficient to construct a PSF model using classical methods.
- Forward - modeling methods (such as WebbPSF) can simulate PSF, but there are differences between the actual observational data and the simulation results, especially in the inner part of the PSF.

### The Solution Proposed in the Paper

To address the above challenges, the authors propose the HybPSF method. The specific steps are as follows:

1. **Extract Star Images**:
- Extract PSF star images from JWST NIRCam observational data, construct mask regions and weight images, and then normalize each image.

2. **Construct Residual Images**:
- Use WebbPSF to generate the corresponding PSF model and subtract it from the star image to obtain the residual image.

3. **Calculate the χ² Equation**:
- Align WebbPSF to the center of each data and use the data, WebbPSF, and mask image to calculate the χ² equation.

4. **Construct the HybPSF Model**:
- Apply polynomial fitting to the coefficients Ci and mi of each chip and combine the WebbPSF and Resi models at the target position. The highest order of the polynomial depends on the number of PSF stars used for PSF reconstruction.

Through these steps, the HybPSF method significantly reduces the difference between the PSF model and the observational data and improves the accuracy of PSF shape parameters (such as ellipticity and size). Specifically, HybPSF improves the accuracy of R² by about 10 times and the accuracy of ellipticity e by 50%.

### Conclusion

The HybPSF method significantly improves the accuracy of the PSF model in JWST NIRCam images by combining WebbPSF simulation and observational data, which is of great significance for high - precision astronomical measurements.