Introduction
The petrofit (PetroFit) is a package for calculating the Petrosian profiles of galaxies and fitting light profiles. It offers end-to-end tools necessary to measure and fit images of galaxies. In this section, we will be defining the morphological properties that can be measured using this package as well as the methods petrofit uses to compute them. This package is based on Crawford (2006) and motivated by the need for an open-source package to compute Petrosian radii and magnitudes.
Sersic Profiles
The Sersic profile is a mathematical function that describes how the intensity of a galaxy varies with radius r from its center. The profile can be described using the following formulation:
\(I\left( r(x,y) \right) = I_e\exp\left\{-b_n\left[\left(\frac{r}{r_{e}}\right)^{(1/n)}-1\right]\right\}\)
The total flux within a radius r is giving by the following integral:
\(L(\leq r) = 2 \pi \int_{0}^{r} I(r^\prime) r^\prime dr^\prime\)
Where:
\(I\) is the intensity or flux at position
(x, y).\(r\) is the radius from the center that corresponds to
(x, y).\(r_{e}\) is the effective radius which is equal to the half-light radius.
\(I_e\) is the flux at the half light radius or \(I_e = I(r_{e})\).
\(n\) is the Sersic index which determines the “steepness” of the profile.
n = 4corresponds to a de Vaucouleurs’ profile whilen = 1gives the exponential profile.
The constant \(b_n\) is defined such that \(r_e\) contains half the total flux, and can be solved for numerically using:
\(\Gamma(2n) = 2\gamma (2n,b_n)\)
In python, \(b_n\) can be calculated as follows:
[1]:
from scipy.special import gammaincinv
n = 4 # Sersic index
b_n = gammaincinv(2. * n, 0.5)
Ellipticity and Elongation
Galaxies that can be described with a Sersic profile are usually not perfectly symmetrical and often display an elliptical distribution. We define two quantities that can be used to describe the elliptical components of galaxies. Ellipticity and elongation are defined as follows:
\(ellip = 1 - \frac{b}{a} = \frac{elong - 1}{elong}\)
\(elong = \frac{a}{b} = \frac{1}{1 - ellip}\)
Where \(a\) is the unit length of the semi-major axis and \(b\) is the corresponding semi-minor axis. A circle corresponds to \(ellip = 0\) and \(elong = 1\). Ellipticity ranges from 0 to 1 while elongation ranges from 1 to infinity. The SExtractor and Photutils packages use elongations for apertures while the Astropy-modeling sub-module uses ellipticity for Sérsic models.
Input Raidus
Elliptical Sersic profiles can be described using \(ellip\) (ellipticity) and \(\theta\) (rotation angle). In such cases, the radius is given by the following expression:
\(r(x, y) = \sqrt{\left(r_{maj}(x, y)\right)^2 - \left(\frac{r_{min}(x, y)}{(1 - ellip)}\right)^2}\)
Where:
\(r_{maj}(x, y) = (x - x_0)\cdot{cos(\theta)} + {(y - y_0)}\cdot sin(\theta)\)
\(r_{min}(x, y) = -(x - x_0)\cdot{sin(\theta)} + {(y - y_0)}\cdot cos(\theta)\)
Sersic Index
A Sersic index determines the degree of curvature of the profile. The larger values of \(n\) get, the more concentrated the profile is at smaller radii. As \(n\) gets larger, the Sersic model tends to a power-law with a slope equal to 5. n = 4 corresponds to a de Vaucouleurs’ profile which describes elliptical galaxies and the cores of spirals well; while n = 1 gives the exponential profile, which models the spiral disks well. n = 0.5 gives a low concentration profile
(Gaussian) that can be used to model inter-cluster glow in galaxy clusters (assuming an elliptical shape), and can be used as an alternative to a Ferrer Profile.
Effective Radius and Intensity
\(r_{e}\) is the effective radius or half-light radius which encloses half of the total flux of the galaxy and can be defined by the following expression:
\(L(\leq r_e) = \frac{1}{2}L(\leq \infty)\).
\(I_e\) is the intensity or flux exactly at \(r_{e}\) and can be defined as \(I_e = I(r_{e})\) and determines the amplitude of the profile. \(I_e\) can be related to the intensity at the center of the Sersic profile (\(I_0 = I(r_0)\)) as follows:
\(I_e = \frac{I_0}{\exp\left\{b_n\right\}}\)
In this plot we see a one dimensional cross-section of a Sersic Profile with \(r_e = 25\) and \(I_e = 1\). Notice how \(I_0\) is equal to 5.3 and not 1. The dotted blue area under the curve contains half of the total flux and is equal to the striped green area in value (the green area extends to infinity).
Curve of Growth
The curve of growth (COG) of a galaxy is the relationship between radius (from the galaxy center) and the total flux within that radius. It represents the cumulative flux enclosed at a given radius. The COG is very sensitive to background noise and fluctuations from the galaxy’s stochastic components. For a galaxy well described by a Sersic profile, the curve of growth is given by \(L(\leq r)\).
Below is a plot of the COG of a Sersic profile with \((I_e=1, r_e=25, n=1)\):
Petrosian Profiles
In petrofit, the Petrosian profile or Petrosian index \(\eta\) is defined as the inverse of the original formulation by Petrosian (1976), and has the property where \(\eta(0) = 1\) at the center of the galaxy and drops to zero at the edge:
\(\eta (r) = \frac{I(r)}{\langle I(r) \rangle} = I(r) \frac{A(r)}{L(\leq r)}\)
Where:
\(\eta (r)\) is the Petrosian at \(r\).
\(I(r)\) is the surface brightness at \(r\).
\(\langle I (r) \rangle\) is the average surface brightness within \(r\).
\(L(\leq r)\) is the total flux within \(r\).
\(A(r)\) is the aperture area. For elliptical apertures: \(A(r) = \pi \cdot (1 - ellip) \cdot r^2\)
Below is a plot Petrosian profile of a Sersic model with \((I_e=1, r_e=25, n=1)\):
Discrete Computation
Given an array of fluxes (“\(L\)”) and corresponding aperture areas (“\(A\)”), the Petrosian profile can be computed discretely as follows:
Estimate the surface brightness by finding the average flux between the current index (\(i\)) and the last index (\(i-1\)). Note that the gap between apertures affects the accuracy of the surface brightness at \(i\), thus it is recommended to use apertures with radii that are incremented by a small number of pixels:
\(I[i] \approx \frac {L[i] - L[i-1]} {A[i] - A[i-1]}\)
Estimate the average surface brightness by taking the flux at \(i\) and dividing it by the corresponding area:
\({\langle I[i] \rangle} = \frac {L[i]} {A[i]}\)
Compute the Petrosian by plugging in the estimated values:
\(\eta [i] = \frac{I[i]}{\langle I[i] \rangle} = \left(\frac {L[i] - L[i-1]} {A[i] - A[i-1]}\right) \cdot \frac {A[i]} {L[i]}\)
In discrete computations, the Petrosian can not be computed at the first index even if it corresponds to the center pixel (\(r[i_0]=0\)). This is because, in real images, the surface brightness of a galaxy is binned into pixels and to get \(I_0\), one would need to infinitely oversample (infinite resolution). In other words, each pixel corresponds to a total surface brightness integrated within the area of the pixel as opposed to the surface brightness at that “location”. As such,
petrofit sets the first Petrosian value to np.nan when returning a discretely computed array of Petrosian indices. Nevertheless, it is known by definition that the Petrosian index at the center of a galaxy is 1 and petrofit takes advantage of this fact and uses interpolation when computing radii internally.
Petrosian Radius and Eta
The Petrosian radius is the radius at which the Petrosian profile is equal to a special Petrosian index \(\eta_{petro}\). In petrofit, the Petrosian radius is defined as the radius where \(\eta = 0.2\). Thus the Petrosian radius and \(\eta_{petro}\) (“eta”) can be defined as:
\(\eta( r_{petro} ) = \eta_{petro} = 0.2\)
Here is a plot of the Petrosian radius and \(\eta=0.2\) for a Sersic profile with \((I_e=1, r_e=25, n=1)\):
Petrosian Total Flux Radius and Epsilon
The Petrosian total flux radius is, as the name implies, the radius which encloses the total flux of the galaxy. Since Sersic profiles extend to infinity, it is not possible to define a radius that contains 100% of the total flux (that would equate to a radius at infinity). It is also not practical to consider minute flux values at very large radii because real images of galaxies contain noise. For these reasons petrofit considers the radius that contains 99% of the galaxy’s light to be the
total flux radius. To attain this radius, we must multiply the Petrosian radius with a constant. This constant is called epsilon and is set to 2 by default. Thus we can define the Petrosian total flux radius as follows:
\(r_{total} = \epsilon \cdot r_{petro}\)
Petrosian Half Light Radius
The Petrosian half-light radius is the radius which contains half of the total flux. This quaintly is especially important because the half-light radius is equal to the Sersic effective radius. The half-light radius can be computed by numerically finding the radius that encloses half of the total flux. Consequently, the Petrosian half-light radius can be defined as the radius which satisfies the following expression:
\(L(\leq r_{hl}) = \frac{1}{2} L(\leq r_{total} )\)
Concentration Index
The concertation index is the ratio of the radii containing a specified fraction of the total flux and can be generalized as follows:
\(C_{r_{n_1}, r_{n_2}} = 5 \cdot log\left(\frac{r_{n_1}}{r_{n_2}}\right)\)
where \(n\) is the radius in which \(n\%\) of the light is enclosed. The concentration index is related to the profile shape and can be used as a proxy for morphology. Some commonly used concertation indices are the following:
\(C_{2080} = 5 \cdot log\left(\frac{r_{20}}{r_{80}}\right)\)
\(C_{5090} = 5 \cdot log\left(\frac{r_{50}}{r_{90}}\right)\)
Petrosian Corrections
To first order, the Petrosian magnitude is independent of redshift as it depends on the ratio of two surface brightnesses, but because it is based on profile shape, changes in the profile due to morphological K-corrections still apply. However, the real strength of Petrosian magnitudes vs. other variable aperture photometry such as Kron magnitudes (Kron 1980) is that the Petrosian radius only depends on light interior to this radius. This aspect leads to small random errors in the measurement of \(r_{petro}\). Nonetheless, the magnitude within \(2 \cdot r_{petro}\), relative to a true, total magnitude, is still profile dependent. Although 99% of the light from an exponential profile (\(n=1\)) is recovered within \(2 \cdot r_{petro}\), only 82% of the light for an \(r^{1/4}\) profile is measured within \(2 \cdot r_{petro}\) (Graham et al. 2005).
Because of this, we need to adjust \(r_{total}\) by finding the appropriate epsilon (\(\epsilon\)) value. Both the consternation index and “correct” \(r_{total}\) are related, thus a relationship between the consternation index and \(\epsilon\) can be derived by simulating Sersic profiles. For ideal images of Sersic models (no noise and no PSF), we can make a plot of \(C2080\) vs Epsilon of simulated profiles with different half-light radii:
The scatter that is seen in higher concentration is due to sampling of the model onto a pixel grid. Profiles with small half-light radii and high concentrations are especially affected because they require a large oversampling factor.
Real images of galaxies contain noise and light from a source is smeared according to the PSF of the instrument. The smearing of light results in the reduction of concentration, causing a skew in the relationship we derived above. To illustrate this, the plot below shows concentration indices for the same profiles as above, but convolved with a Hubble Space Telescope F105W PSF, sampled at 60 mas. It also shows the same polynomial as above for reference:
Notice how profiles with smaller half-light radii are affected (skewed) the most by the PSF convolution. For this reason, we must make corrections that account for the half-light radius and the shape of the PSF.
List of Galaxy Properties
The following galaxy properties can be computed using petrofit as described in this section:
Sersic Parameters
Petrosian Parameters
References
Bershady, M. A., Jangren, A., and Conselice, C. J., “Structural and Photometric Classification of Galaxies. I. Calibration Based on a Nearby Galaxy Sample”, The Astronomical Journal, vol. 119, no. 6, pp. 2645–2663, 2000. doi:10.1086/301386.
Crawford, S. M., “The distribution of star-forming galaxies in intermediate redshift galaxy clusters”, PhDT, 2006.
Graham, A. W., “Total Galaxy Magnitudes and Effective Radii from Petrosian Magnitudes and Radii”, The Astronomical Journal, vol. 130, no. 4, pp. 1535–1544, 2005. doi:10.1086/444475.
Graham, A. W. and Driver, S. P., “A Concise Reference to (Projected) Sérsic R1/n Quantities, Including Concentration, Profile Slopes, Petrosian Indices, and Kron Magnitudes”, Publications of the Astronomical Society of Australia, vol. 22, no. 2, pp. 118–127, 2005. doi:10.1071/AS05001.
Kent, S. M., “CCD Surface photometry of field galaxies.I. Observations.”, The Astrophysical Journal Supplement Series, vol. 56, pp. 105–141, 1984. doi:10.1086/190978.
Kron, R. G., “Photometry of a complete sample of faint galaxies.”, The Astrophysical Journal Supplement Series, vol. 43, pp. 305–325, 1980. doi:10.1086/190669.
Petrosian, V., “Surface Brightness and Evolution of Galaxies”, The Astrophysical Journal, vol. 210, p. L53, 1976. doi:10.1086/182301.