A forward-modelling approach to overcome PSF smearing and fit flexible models to the chemical structure of galaxies

Benjamin Metha,Simon Birrer,Tommaso Treu,Michele Trenti,Xuheng Ding,Xin Wang

2024-03-13

Abstract:Historically, metallicity profiles of galaxies have been modelled using a radially symmetric, two-parameter linear model, which reveals that most galaxies are more metal-rich in their central regions than their outskirts. However, this model is known to yield inaccurate results when the point-spread function (PSF) of a telescope is large. Furthermore, a radially symmetric model cannot capture asymmetric structures within a galaxy. In this work, we present an extension of the popular forward-modelling python package LENSTRONOMY, which allows the user to overcome both of these obstacles. We demonstrate the new features of this code base through two illustrative examples on simulated data. First, we show that through forward modelling, LENSTRONOMY is able to recover accurately the metallicity gradients of galaxies, even when the PSF is comparable to the size of a galaxy, as long as the data is observed with a sufficient number of pixels. Additionally, we demonstrate how LENSTRONOMY is able to fit irregular metallicity profiles to galaxies that are not well-described by a simple surface brightness profile. This opens up pathways for detailed investigations into the connections between morphology and chemical structure for galaxies at cosmological distances using the transformative capabilities of JWST. Our code is publicly available and open source, and can also be used to model spatial distributions of other galaxy properties that are traced by its surface brightness profile.

Instrumentation and Methods for Astrophysics,Astrophysics of Galaxies

What problem does this paper attempt to address?

The problems that this paper attempts to solve mainly focus on two aspects:

1. **Overcoming the blurring effect of the point spread function (PSF)**: When dealing with the metallicity distribution of high - redshift galaxies, traditional methods lead to inaccurate measurement results due to the large point spread function (PSF) of telescopes. Especially when the size of the PSF is equal to or larger than the size of the galaxy, the traditional radially symmetric model will underestimate the metallicity gradient, thus causing deviation.

2. **Fitting non - radially symmetric chemical structure models**: Traditional radially symmetric models cannot capture the asymmetric structures inside galaxies, such as the metallicity changes caused by spiral arms and merger events. These asymmetric structures are more common in high - redshift galaxies, so a more flexible modeling method is required to describe these complex chemical structures.

To solve these problems, the author proposes a forward - modeling - based method and extends the functions of the popular Python package `lenstronomy`. This method can:

- **Accurately recover the metallicity gradient**: Even when the PSF is large or the pixelation effect is significant, the metallicity gradient of the galaxy can be accurately recovered through forward modeling.
- **Fit complex metallicity distributions**: It can fit non - radially symmetric metallicity distributions, thereby better describing the complex structures inside galaxies, such as spiral arms and merger events.

Specifically, the author achieves this goal through the following steps:

1. **Introducing light - weighted convolution**: Considering that metallicity is a light - weighted quantity, the author introduces the light - weighted convolution formula to take into account the effects of PSF blurring and pixelation. The formula is as follows:
\[
\tilde{Z}(\mathbf{x})=\frac{\int Z(\mathbf{x}')I(\mathbf{x}')K(\mathbf{x}' - \mathbf{x})d\mathbf{x}'}{\int I(\mathbf{x}')K(\mathbf{x}' - \mathbf{x})d\mathbf{x}'}
\]
where \( Z(\mathbf{x}) \) is the intrinsic metallicity of the galaxy, \( I(\mathbf{x}) \) is the surface brightness distribution of the galaxy, and \( K(\mathbf{x}) \) is the kernel function of the PSF.

2. **Multi - component model**: For galaxies composed of multiple surface brightness distributions, the author proposes a method for calculating the total metallicity distribution:
\[
Z(\mathbf{x})=\frac{\sum_{i = 1}^{n}I_i(\mathbf{x})Z_i(\mathbf{x})}{\sum_{i = 1}^{n}I_i(\mathbf{x})}
\]
This formula allows the metallicities of different components to be mixed together to simulate complex galaxy structures.

3. **Forward - modeling framework**: By using the forward - modeling framework, the author can verify the effectiveness of this method in simulated data. The results show that even when the PSF is large or the pixelation effect is significant, this method can still accurately recover the metallicity gradient of the galaxy.

In conclusion, the method proposed in this paper aims to overcome the limitations of existing models in dealing with the metallicity distribution of high - redshift galaxies, providing a more accurate and flexible modeling tool, which is helpful for in - depth understanding of the chemical evolution history of galaxies.