An Empirical Calibration of the Tip of the Red Giant Branch Distance Method in the Near Infrared. I. HST WFC3/IR F110W and F160W Filters

The Tip of the Red Giant Branch (TRGB)-based distance method in the I band is one of the most efficient and precise techniques for measuring distances to nearby galaxies (D<= 15 Mpc). The TRGB in the near infrared (NIR) is 1 to 2 magnitudes brighter relative to the I band, and has the potential to expand the range over which distance measurements to nearby galaxies are feasible. Using Hubble Space Telescope (HST) imaging of 12 fields in 8 nearby galaxies, we determine color-based corrections and zero points of the TRGB in the Wide Field Camera 3 IR (WFC3/IR) F110W and F160W filters. First, we measure TRGB distances in the I band equivalent Advanced Camera System (ACS) F814W filter from resolved stellar populations with the HST. The TRGB in the ACS F814W filter is used for our distance anchor and to place the WFC3/IR magnitudes on an absolute scale. We then determine the color dependence (a proxy for metallicity/age) and zero point of the NIR TRGB from photometry of WFC3/IR fields which overlap with the ACS fields. The new calibration is accurate to ~1% in distance, relative to the F814W TRGB. Validating the accuracy of the calibrations, we find that the distance modulus for each field using the NIR TRGB calibration agrees with the distance modulus of the same fields as determined from the F814W TRGB. This is a JWST preparatory program and the work done here will directly inform our approach to calibrating the TRGB in JWST NIRCam and NIRISS photometric filters.

1. TIP OF THE RED GIANT BRANCH DISTANCES Distances to galaxies are fundamental measurements and are essential to placing myriad galaxy properties (e.g., mass, luminosity, size, etc.) on absolute scales.In addition, extragalactic distance measurements are critical rungs in the cosmic distance ladder enabling local measurements of the Hubble Constant (H 0 ; e.g., Freed- Beaton et al. 2016;McQuinn et al. 2017).It is one of the most efficient and precise methods, it can be employed in almost any galaxy type requiring only that they host a well-defined population of Red Giant Branch (RGB) stars, and the physical principles underlying the method are well understood (e.g., Salaris & Cassisi 1997;Madore et al. 1997;Salaris et al. 2002).
The TRGB is defined as the end of the RGB stellar evolutionary phase.A low-mass star (< 2M ⊙ ) ceases burning hydrogen in its core as it exits the mainsequence and begins burning hydrogen in a shell around its inert, degenerate helium core.The star evolves up the RGB by expanding, increasing in luminosity, and dumping helium ash onto a still degenerate helium core.The core pressure and temperature build until a critical temperature of ∼ 10 8 K is reached, at which point helium fusion ignites, lifting the degeneracy.This point is defined as the TRGB.The sudden lifting of the degeneracy causes a release of energy referred to as the Helium Flash, after which the star moves onto either the Horizontal Branch phase, decreasing in brightness and increasing in surface temperature, or onto the Red Clump.Because the Helium Flash occurs at a specific stellar core mass, the TRGB luminosity is stable and predictable across stars with a wide range of properties and can be used as a standard candle (e.g., Lee et al. 1993;Beaton et al. 2018).
While the bolometric luminosity of the TRGB is approximately constant, the apparent brightness of the TRGB in specific bandpasses will vary on a star-to-star basis due to differences in intrinsic stellar properties such as chemical content (i.e., metallicity), and, to a lesser extent, age.The impact of these variations in the apparent brightness of stars at the TRGB correlates with color and must be carefully quantified in order for the TRGB to be used as a precision distance indicator.
TRGB distances are typically measured in the I band where the brightness of the TRGB has been found to be approximately constant and metallicity/age-based corrections are modest, albeit with somewhat larger metallicity/age dependence for more metal-rich stellar populations (e.g., Lee et al. 1993;Madore & Freedman 1995;Sakai et al. 1996).This metallicity/age dependence is well understood and is accounted for in a number of TRGB calibrations using color-based corrections (e.g., Méndez et al. 2002;Makarov et al. 2006;Rizzi et al. 2007;Jang & Lee 2017;Freedman 2021).
The TRGB is identified as a sharp discontinuity in the stellar luminosity function (LF), something that can be measured more accurately if the color correction is applied first.Two methods are commonly used to identify the location of the TRGB: (1) the Sobel edge-detection filter based on a zero-sum Sobel kernel, which is itself a discrete approximation of a gradient (e.g., Lee et al. 1993;Madore & Freedman 1995;Sakai et al. 1996), and (2) a maximum-likelihood (ML) fitter which assumes a simple power law with a cutoff in the TRGB region plus a power law of a second slope for a stellar population brighter than the TRGB (e.g., asymptotic giant branch (AGB) stars; e.g., Méndez et al. 2002;Makarov et al. 2006).Once the TRGB luminosity has been identified and measured, the distance to the galaxy can be determined using a zero point calibrated from previous studies (e.g., Lee et al. 1993;Rizzi et al. 2007;Jang & Lee 2017).
As one moves from the I band to longer wavelengths in the near-infrared (NIR), the TRGB increases in brightness by ∼1-2 mag offering observational gains.If the NIR TRGB can be robustly calibrated, the TRGB-based distance method has the potential to be applied to galaxies at larger distances and over a greater volume.In addition, shorter observation times will be required to measure high-precision distances, further increasing the efficiency of the method.A TRGB calibration in the NIR is timely given the new and upcoming space-based telescopes designed to operate in the NIR, i.e., the James Webb Space Telescope (JWST) and the Nancy Grace Roman Observatory.
In contrast to the TRGB in the I band, the luminosity of the NIR TRGB exhibits a greater metallicity/age dependence in the sense that the TRGB stars have a larger range in brightness as a function of metallicity/age (e.g., Ferraro et al. 2000;Valenti et al. 2004;Salaris & Girardi 2005;McQuinn et al. 2019).There have been previous calibrations of the color dependence of the NIR TRGB through a combination of theoretical stellar evolutionary models and photometric data (e.g., Dalcanton et al. 2012a;Wu et al. 2014;Serenelli et al. 2017;Hoyt et al. 2018;Madore et al. 2018;Górski et al. 2018;Durbin et al. 2020;Freedman et al. 2020;Madore et al. 2023a).However, current theoretical models from stellar evolution libraries vary in their zero point and metallicity/age-dependent slopes of the NIR TRGB, and are not always self-consistent across different wavelengths (e.g., Serenelli et al. 2017;Durbin et al. 2020).
In this work, we calibrate the NIR TRGB in the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) F110W and F160W filters.Our calibration is entirely empirical as we do not use theoretical stellar evolution libraries, and is tied to the well understood F814W TRGB (Jang & Lee 2017;Freedman 2021).In addition, this program is a JWST preparatory program and the work done here will directly inform our approach to calibrating the TRGB in several JWST Near Infrared Camera (NIRCam) and Near Infrared Imager and Slitless Spectrograph (NIRISS) photometric filters (JWST-GO-1638; PI McQuinn).
The paper is organized as follows: We present the HST observations, data reduction, photometry, artificial star tests (ASTs), and foreground-extinction corrected color magnitude diagrams (CMDs) in § 2. We describe the TRGB-based distance method used to determine the distances to our sample, which anchor our calibration, and briefly compare our distance measurements to TRGB distances in the literature in § 3. Next, we discuss our methods for calibrating the NIR TRGB color-based correction and zero point and we present the NIR TRGB calibration in the F110W and F160W filters in § 4. Finally, we discuss our results in the context of previous calibrations and summarize our findings in § 5 and § 6, respectively.

Galaxy Sample
The luminosity of the TRGB has a dependence on stellar metallicity/age, which changes as a function wavelength.In the F814W filter, the dependence is modest but is still important to consider not only in its calibration, but also in observing strategies.Specifically, the TRGB in F814W is relatively flat for metal poor stellar populations (e.g., Beaton et al. 2016), so the outer stellar fields of galaxies which are thought to be primarily metal poor are commonly targeted for TRGB distance measurements.Outer stellar fields offer additional advantages including (1) reduced impact from crowding effects where the apparent brightness of a star can be biased to higher or lower magnitudes in regions of relatively high stellar densities (e.g., in stellar disks), (2) there are fewer young stellar populations that can blur the TRGB, and (3) there is less internal extinction.
However, there is observational evidence that suggests galaxy halos can be populated with stars that span a wider range of metallicity.Recent studies of spiral halos have shown metallicity gradients (e.g., Ibata et al. 2014;Gilbert et al. 2014;Cohen et al. 2020), high median metallicities (i.e., [Fe/H] > −1.2 dex) as far as 100 kpc from galaxy centers (e.g., Conroy et al. 2019), significant field-to-field variations in the mean metallicity of RGB stars for an individual galaxy (including along the major vs. minor axes), and metallicity variations of 1 dex across a sample of 6 spiral halos (Monachesi et al. 2016).
In the NIR, because the luminosity of the TRGB has a more significant dependence on metallicity/age (e.g., Salaris & Girardi 2005;Serenelli et al. 2017;McQuinn et al. 2019), the impact of variations in metallicity/age on the TRGB luminosity must be carefully accounted for.We have crafted our calibration sample to explicitly account for these metallicity/age variations.Our calibration of the NIR TRGB is based on observations of galaxies with a wide range of metallicities.The observations include new imaging of outer stellar fields in four nearby galaxies obtained for this program, and archival imaging of four low-mass galaxies.We supplemented the new observations with the archival data to increase the metallicity baseline of the sample.
For these new observations, we made use of the separation of the HST Advanced Camera System/Wide Field Camera (ACS/WFC; Ford et al. 1998) and WFC3/IR on the sky to image two outer stellar fields per galaxy with WFC3 in the parallel observing mode.This observing strategy enables us to account for potential changes in the stellar content of outer stellar fields.Field one was selected with ACS, and the roll angle was set such that WFC3/IR imaged another region of the outer stellar field.A second visit was planned for each field such that orientation of ACS and WFC3/IR were flipped.Figure 1a presents GALEX images for each galaxy with footprint targets of ACS (green) and WFC3/IR (white).
We placed tight roll angle constraints on our observations to ensure that each field was imaged with both instruments.
The archival data consist of HST observations of 4 nearby dwarf galaxies with overlapping images from the same set of filters as in our new observations.The archival data include imaging of both the stellar disk and part of the halos.Figure 1b shows GALEX images of the 4 dwarf galaxies with footprint targets of ACS (green) and WFC3/IR (white).The combined observations provide a uniform data set for our calibration.
Table 1 presents a summary of the HST observations including galaxy coordinates and observation details.The new observations were obtained as part of the JWST preparatory program HST-GO-15917.The archival observations were obtained over three different programs: ACS imaging is from HST-GO-10210 and HST-GO-10915 (ACS Nearby Galaxy Survey Treasury; Dalcanton et al. 2009), and the WFC3/IR imaging is part of HST- GO-16162 (LUVIT;Boyer et al. in prep.).

Image Alignment
Precise, relative astrometric alignment between all images is required for high-fidelity photometry.We used TweakReg in the DrizzlePac 3.0 package (Hack et al. 2013;Avila et al. 2015) to align our ACS and WFC3/IR raw images to a common world coordinate system (WCS).TweakReg uses a point source detection algorithm to extract astrometric source catalogs for each input image.The alignment algorithm calculates the best fits for shifts, rotations, and scale factors between the images and a reference image based on matching the point sources across images.Images in one filter are first aligned to each other and are then combined into a single drizzled image using AstroDrizzle in DrizzlePac 3.0.The drizzled image is then used as a reference image to which images in all other filters are aligned.
We selected the drizzled F814W image as the reference image, as this image is a mid-wavelength image in our data set and the ACS data have a higher resolution than WFC3.TweakReg then updates the World Coordinate System (WCS) in the fits header for the non-reference images.
We required the TweakReg alignment to be ≤ 0.1 pixels for ACS images, and ≤ 0.15 pixels for WFC3/IR.We used TweakReg to align all the F814W images to each other and then used AstroDrizzle to generate a drizzled F814W reference image.Next, we aligned all remaining CTE-corrected flc.fitsACS images and the flt.fitsWFC3/IR images to the drizzled F814W image.Finally, we used AstroDrizzle to make drizzled images in the three other filters (F606W, F110W, and F160W) and visually inspected their alignment using DS9 (Joye & Mandel 2003).

Multi-Camera Photometry
To robustly calibrate the TRGB in the NIR we required a catalog of matched stars from all of the data for each galaxy.We performed simultaneous cross-camera point spread function (PSF) photometry on the ACS and WFC3/IR images.PSF photometry on the ACS CTE-corrected and WFC3/IR flat-field calibrated images was executed using DOLPHOT, a modified version of the WFPC2-specific HSTphot package that includes instrument-specific support for ACS and WFC3 (Dolphin 2000(Dolphin , 2016)).DOLPHOT uses a deep, cosmic-ray cleaned and drizzled image, produced using AstroDrizzle, as a reference image to identify point sources and alignment solutions.The drizzled F814W was selected as the reference image as ACS has a larger field of view, higher spatial resolution, and a smaller PSF relative to WFC3/IR.DOLPHOT then locates and iteratively measures the flux of the same point sources in each raw image.
We adopt the DOLPHOT input parameters used in the Panchromatic Hubble Andromeda Treasury (PHAT) and Panchromatic Hubble Andromeda: Triangulum Extended Region (PHATTER) surveys (Williams et al. 2014(Williams et al. , 2021)).We conducted extensive input parameter tests by varying several parameters away from those presented in the PHAT survey.We find in particular that for the parameter FitSky, which sets the method DOLPHOT uses to estimate the local background around each source, a value of 2 produces the best results in both the number and quality of photometered sources in fields with low or high stellar densities.This is in agreement the with the results from the PHATTER survey (Williams et al. 2021).

Artificial Star Tests
Figure 1a.Clockwise from the top left: M81, NGC253, NGC300, NGC2403.Shown are color images using GALEX nearultraviolet (NUV; blue) and far-ultraviolet (FUV; red) filters, and an average of the NUV and FUV filters (green), with the HST ACS and WFC3/IR instrument footprints indicated in green and white, respectively.Images are north up, east left.GALEX images were obtained courtesy of NASA/JPL-Caltech.
We performed ASTs using DOLPHOT in order to quantify photometric completeness and photometric errors due to blending in our stellar catalogs.We injected 200k − 500k ASTs in each field that spanned the CMD parameter space of all four filters (F606W, F814W, F110W, and F160W).We generated the artificial star input lists required to run ASTs using the Bayesian Extinction And Stellar Tool (BEAST; Gordon et al. 2016).While DOLPHOT has the ability to generate artificial star input lists for more than two bands, the artifical star input lists that the BEAST generates should be higherfidelity, since DOLPHOT uses a simple scaling relation between the bands.Briefly, the BEAST uses an observed photometry catalog and a grid of stellar evolution models to predict sets of spectral energy distributions (SEDs).A BEAST artificial stars input list is then cre-ated by sampling the SEDs ensuring that the full range of the model fluxes in each filter is spanned.

Photometric Catalog
We created high-fidelity stellar catalogs by applying quality cuts to the DOLPHOT photometric outputs.The initial output from DOLPHOT includes parameters which provide information about the quality of the recovered points sources.We use a number of these parameters to cull sources that we define as good stars.First, we remove sources with the error flag parameter ≥ 4. The error flag tells us how well a source was recovered in the image.We also remove sources with the object type parameter < 2, which tells us about the confidence that a measured object is a star and is not elongated or too faint.Next, we exclude sources with a signal-to-noise ratio (SNR) ≤ 5 in each filter.We then apply cuts based on the crowding parameter, which is a measure of the difference between the initial brightness of a star and its value after subtracting the flux from neighboring sources.We use the sharpness parameter to remove sources that are not point like (i.e., cosmic rays or extended sources such as background galaxies).
We make our quality cuts empirically and uniformly across all fields, on a per-filter basis.We iteratively adjust the quality parameters until the diffraction spikes from bright foreground stars are minimized in the photometry catalog.Shown in Figure 2 is an example of the criteria used for NGC2403-2, and is representative of our data.Here we show the stars that pass our quality cuts (black points), and those that fail to meet our criteria (gray points).In F606W and F814W, the maximum crowding value is 0.1 mag, and in F110W and F160W the maximum crowding value allowed is 0.2 mag.In all filters, we chose to include only sources with a value of sharpness in the range −0.2 ≤ sharpness ≤ 0.2 mag.

Color-Magnitude Diagrams
Figure 3 shows extinction corrected optical CMDs from the spatially-clipped ACS FOV with for each galaxy (as shown in green in Figures 1a and 1b).Extinction corrections for all the galaxies in the sample were taken from the dust maps of Schlegel et al. (1998) with a re-calibration from Schlafly & Finkbeiner (2011).A V values are listed in Table 1.Extinction corrections are small with the maximum values of A F606W = 0.19, A F814W = 0.12, A F110W = 0.069, and A F160W = 0.041 for Antlia Dwarf, and average values for the entire sample of A F606W = 0.052, A F814W = 0.032, A F110W = 0.018, and A F160W = 0.011.Table 2 summarizes the  20 21 22 23 24 25 26 27 28   F606W   0 filter coefficients and final filter specific total extinction values.
Each galaxy has a well-populated RGB and the TRGB is identifiable by eye as a distinct discontinuity at the upper part of the RGB.We also note the presence of mainsequence, blue helium-burning (BHeB), and red heliumburning (RHeB) stars, as well as some AGB stars in all fields except for NGC253.We select spatial regions in either the stellar halos or outside the dense stellar disks to minimize contamination from these young, non-RGB stellar populations in ACS, as described in Section 2.2.6.
Figure 4 presents extinction corrected NIR CMDs from the full WFC3/IR FOV for each galaxy (as shown in white in Figures 1a and 1b).All NIR CMDs display a similar RGB sequence but the RGB has a narrower color baseline and the TRGB is more steeply sloped compared to the optical CMDs as predicted from stellar isochrones (e.g., McQuinn et al. 2019).

Spatial Cuts
The TRGB is best identified in a CMD with a wellpopulated RGB sequence and minimal contamination from non-RGB stars (i.e., AGB or RHeB stars).Because the TRGB is based on measuring the discontinuity in the LF, retention of non-RGB stars in the photometry with both colors and apparent brightnesses near the RGB can blur the break in the LF.
Thus, we apply spatial cuts to our observations to reduce contamination from non-RGB populations.Our spatial cuts are based on a 2-step approach and use information from the ACS stellar catalogs and the location of the ACS FOV relative to the main stellar disks of the galaxies.First, we exclude regions where we find evidence of significant young stellar populations based on GALEX images (see Figures 1a and 1b).Second, we refine the spatial cuts by making CMD selections of any remaining young stellar populations, plotting their spatial distributions, and excluding any region with high fractions of young stars.
This 2-step approach is modified for the four dwarf galaxies with archival observations in which the ACS FOV covers the central region and outer stellar field.We adopt ellipse structural parameters from McConnachie (2012) for Antlia Dwarf, Sloan Digital Sky Survey Data Release 16 (SDSS DR16; Ahumada et al. 2020) for UGC 8833, and SDSS DR12 (Alam et al. 2015) for UGC 9128 and for UGC 9240.We then generate elliptical annuli starting from the center and moving outward, and count the number of MS stars and RGB stars in each annuli.We use the ratio of MS stars to RGB stars to optimize our final FOV selection.
The left panel of Figure 5 presents an example of the spatial cuts applied to M81-1.The point sources in the upper left corner of the image (green) are separated from the rest of the sources (orange).The right panel presents the resulting CMD of the two regions color-coded in the same manner.There is a clear separation of the stellar populations in the these two regions.

MEASURING THE DISTANCES IN THE F814W LUMINOSITY FUNCTION
In this section, first we introduce the calibration of the F814W TRGB magnitude as a function of the F606W−F814W color.Second, we discuss two of the methods developed to measure the magnitude of the TRGB, and justify the method we employ to measure the TRGB in our sample.Third, we adopt our F814W TRGB zero point and we present our final F814W-based distance measurements.Finally, we compare our results to those reported in the literature.

Color Correction of the TRGB in F814W
The first step taken to identify the TRGB luminosity is to correct for the well-known, modest dependence of the TRGB on metallicity/age.We account for this dependence by using color-based corrections.We adopt the only the color correction (not the zero point) from Jang & Lee (2017), named the quadratic transformation (QT) system, derived for the HST filters which has been recently shown as most accurately characterizing the color-based dependence of the TRGB (Hoyt 2021).We reproduce the calibration here, over the applicable color range F606W-F814W> 1.5, for convenience: where F814W is the original magnitude, F814W rect is the color-corrected magnitude, and α = 0.159 ± 0.010 and β = −0.047± 0.020 are the best-fit parameters for the quadratic and linear terms, respectively.For F606W-F814W≤ 1.5, the TRGB is consistent with being flat (i.e., at a constant magnitude) and no correction is applied.In Eq. 3 of Jang & Lee (2017, Eq. 1 in this text) the authors identify the pivot point for the color-based correction as the average color of the TRGB (F606W−F814W = 1.1 mag) in their sample.The use of the average color as the pivot point accounts for differences between the color of the stars in the observed galaxy and the color used to anchor the zero point in the calibration.We apply this correction on a per-star basis prior to measuring the location of TRGB (i.e., we rectify the TRGB).Fitting for the TRGB after applying the correction has the effect of increasing the sharpness of the discontinuity in the LF.Note that this color-correction only applies to 2 of the 12 fields as most of the data had colors F606W-F814W< 1.5.

TRGB Measurements Methods
Several techniques have been developed to measure the F814W TRGB magnitude.Of them, the technique first applied to identify the discontinuity in the LF was  the Sobel filter edge detection technique (Lee et al. 1993;Sakai et al. 1996Sakai et al. , 1997)).This technique is based on applying a Sobel Kernel in the form [−1, 0, +1], an approximation of a derivative, to a binned LF.Convolving the Sobel Kernel with the binned LF results in a response function indicating the greatest change in number of stars from one magnitude bin to the next.
We employ a more sophisticated Bayesian Maximum Likelihood (ML) technique to measure the F814W TRGB magnitude in each galaxy.The ML technique uses a parametric form of the RGB LF to fit the observed LF.The probability distribution takes into account the photometric uncertainty distribution and completeness function from ASTs (see Makarov et al. 2006;Madore et al. 2023b, for a full discussion of the technique, and for a detailed discussion of quantifying tip uncertainties respectively).We adopt the same theoretical LF form used in Makarov et al. (2006): where A is the slope of the RGB with a normal prior of µ A = 0. uncertainty is based on the range of solutions returning the log of the probability, log (P ), within 0.5 of the maximum.
When fitting for the TRGB in the F814W filter, we select stars consistent with the luminosity and colors of RGB stars on a per-field basis as a well constrained RGB sequence will yield the most robust TRGB measurements.If stars blueward of the RGB sequence (e.g., BHeB, RHeB, or main sequence stars) remain in the CMD, they can bias the stellar LF and reduce the strength of the break in the LF associated with the TRGB.The color ranges we selected to constrain the RGB stars are shown in Figure 3 as orange boxes with stars included as black points, while the stars excluded are shown as gray points.Our color selections are empirically determined such that the blue edge of the box follows the blue edge of the RGB, while the red edge is set to encompass the reddest sources.Even though the ML approach takes into account the possible presence of a secondary population, we nonetheless explored restricting our fits to stars within a red color limit (thereby reducing the potential number of redder AGB stars).When we compared the TRGB values measured without applying this additional color restriction, we find the residuals between the fits are consistent with zero within the statistical uncertainties.We chose to use the full population of redder stars as this selection method avoids artificially biasing the AST-based completeness function used in the ML fitting.
We run the ML technique on the foregroundextinction corrected photometry, after applying the color-based correction described in Eq. 1.The location of the TRGB magnitude identified with the ML technique (indicated as a horizontal red line in Figure 3) is shown in each CMD with uncertainties from the ML fit shaded in blue.Table 3 presents the best-fit TRGB luminosities from the ML technique for the sample.The uncertainties include statistical uncertainties from the ML probability distribution function.
This ML technique is preferred over the Sobel filter technique for measuring the F814W TRGB magnitude as it does not require binning the LF, and considers the photometric measurement uncertainties and completeness.The ML technique also fits the F814W LF while accounting for the AGB population.Compared with the Sobel filter technique, the ML uncertainties are typically lower, primarily due to the Sobel filters dependence on bin width.In Section 4, however, we employ the Sobel filter technique to measure the slope and zero point of the NIR TRGB for two reasons.First, we perform our calibration based on a stacked stellar catalog of every field in the sample.As of this writing, the ML method does not have a way to simultaneously handle the results of ASTs from multiple fields.It assumes any input photometry is from one field, and that the ASTs characterize the photometric noise properties of that field only.Second, the Sobel filter proves to be an excellent tool to optimize and ultimately determine our best-fit slope to the TRGB.

Distances
We use the zero point calibration from Freedman (2021) to convert TRGB luminosities to distance moduli (µ).Their final F814W zero point is based on four independent calibrations of the TRGB absolute magnitude: the megamaser galaxy NGC 4258, Galactic Globular clusters (i.e., using Gaia DR3), the Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC).The independent calibrations are found to be internally selfconsistent at the 1% level.We therefore choose to adopt this F814W TRGB zero point measured.We note that, since our NIR TRGB calibration is tied to the F814W TRGB calibration, the results of our calibration can be updated to account for improvements to the F814W TRGB zero point and uncertainties.
Statistical uncertainties from the ML fit are reported as the total uncertainty on the TRGB measurement (see §4.3 for details).Final distance moduli and distances are also reported in Table 3 with statistical uncertainties on each fit.For the four galaxies with two observed fields, the best-fitting distances agree within the uncertainties.

Comparison with TRGB Distances from the Literature
Since we require robust F814W TRGB-based distances to anchor the NIR TRGB calibration, we compared distance moduli between our measurements and TRGB-based measurements in the literature.Based on the values reported in the NASA/IPAC Extragalactic Database (NED)1 our measurements agree within 1σ of the mean for each target.Note that these methods include slightly different approaches to measuring and calibrating the TRGB.
We included a more direct test and focus on TRGB measurements from the literature that match our approach.Specifically, we selected distance moduli values from the literature sources to compare with our data based on the following criteria: TRGB-based distance measurements were made using HST F814W data and color-based corrections, and zero points used in the measurement were empirical.We compared the individual measurements to our measurements.Again, we found that the measurements agree within their uncertainties.These literature values are tabulated in Table 4.

IDENTIFYING AND CALIBRATING THE TRGB IN THE NIR
While the TRGB magnitude in the F814W filter is observed to have a modest dependence on metallicity, the NIR TRGB is known to have a stronger dependence.Specifically, it has been shown that the brightness of the TRGB increases with increasing metallicity in the NIR (Valenti et al. 2004;Salaris & Girardi 2005;McQuinn et al. 2019).Observationally, this effect manifests as the TRGB luminosity increasing at redder colors.
To calibrate the slope and zero point of the NIR TRGB in the F110W and F160W filters we made a combined stellar catalog from the entire sample to provide a broad metallicity, hence color, baseline.We found that the stars in a CMD of any individual field did not cover a wide enough color range to robustly characterize the slope of the NIR TRGB for calibration purposes.Additionally, by combining the individual catalogs we increased the total number of stars in the CMD.Shown in the top panels of Figure 6 are the IR CMDs generated from the combined stellar catalogs.The combined stellar catalog is created after converting all foregroundextinction corrected apparent magnitudes to absolute magnitudes using the distance moduli measured from the F814W TRGB listed in Table 3.
Here, we provide a mathematical basis for the TRGB calibration, introduce our method for identifying the TRGB in the F110W and F160W filters, and then describe our methods for measuring and calibrating its slope as a function of color, and then for measuring the zero point of the NIR TRGB in both filters.

Mathematical Background on the TRGB Characterization
The magnitudes of TRGB stars in the NIR exhibit an approximately linear relationship with color such that their brightness increases from bluer to redder colors.Like the F814W TRGB, the color dependence of the NIR TRGB can be calibrated via a color-based correction (e.g., Rizzi et al. 2007;Wu et al. 2014;Jang & Lee 2017).By characterizing its color dependence, the NIR TRGB can be rectified to a single magnitude and subsequently employed as a robust distance indicator.
We characterize the color-based correction via the parameter β (i.e., the slope of the NIR TRGB), the median F110W-F160W color of the TRGB stars, which sets the pivot point at which the TRGB is rectified to a single magnitude (i.e., the zero point), and the absolute magnitude of the TRGB at the pivot point.We adopt for our color-based correction the following equation: where β is the slope of the TRGB in the chosen NIR filter, F110W-F160W is the color of the individual stars,   (Rizzi et al. 2007, ; the most current calibration at the time of their studies).While their adopted calibration differs from ours, the methodology remains the similar.
and γ is the median F110W-F160W color of the TRGB stars, M NIR is the color-corrected TRGB magnitude, and M TRGB NIR,γ is the magnitude of the TRGB at the pivot point in the NIR filters.When applied to the apparent magnitude of stars the equation is: where m NIR,rect and m NIR,orig are the extinction corrected rectified and original magnitudes of the stars in a NIR filter.
Mathematically, one can transform directly from the calibration in one bandpass, F110W for example, directly to another, F160W, for the same color F110W-F160W (e.g., Wu et al. 2014;Serenelli et al. 2017;Madore & Freedman 2020).The transformation can also be applied between different bandpasses (see Madore & Freedman 2020).Briefly, for the filters F110W and F160W and using Eq.3: where M TRGB F110W,γ is the absolute magnitude in F110W at the pivot point γ, and β is the slope of the TRGB in F110W.We can then determine the transformation from F110W to F160W.First we define F160W as: Next we substitute in for F110W with Eq. 5.
then expand the terms on the right hand side, combine like terms and expand βγ, and combine like terms once more to reach the calibration equation for F160W in the form of Eq. 5.
F160W = M TRGB F110W,γ − γ (10) The resulting transformation from F110W to F160W can then be characterized by the coefficients While it a straightforward process to transform the parameters α, β, and M TRGB NIR from F110W to the F160W, we chose to instead calibrate each filter independently and only then use the transformation as a check on our empirical results.The final values we report for M TRGB F110W,γ , M TRGB F160W,γ , and β are derived independently and empirically from individual fits to the data (see §4.2).Running the calibrations separately for each filter enables us to consider the photometric uncertainties and their impact on the calibration independently.We then apply the mathematical transformation to the F110W best-fit calibration values as a consistency check to determine the theoretical calibration parameter values in F160W (see §4.3).

NIR TRGB Color Dependence and Zero Point
Our method for calibrating the color-based correction for the NIR TRGB simultaneously measures the slope and zero point of the TRGB.First, we identify a fiducial slope by eye in the NIR CMDs and use that value as a prior in the analysis.We make CMD cuts in the range of the RGB to reduce non-RGB contamination blueward of the RGB, similar to our method for constraining the RGB in the F814W filter.We identify the median color of the TRGB within a color range 0.73 ≤ F110W-F160W ≤ 1.0 mag at γ = 0.92 mag.
Second, we employ an algorithm based on the Nelder-Mead formalism via Scipy's minimization routine.
We construct a binned LF from the rectified CMD within the defined TRGB color range using a bin width of 0.01 mag and smooth over the LF with the Gaussian-weighted locally estimated scatterplot smoothing (GLOESS) method (Persson et al. 2004;Monson et al. 2017).GLOESS is a non-parametric weighted least squares method in which the fitting is done locally and the weights are described by a Gaussian profile with the bandwidth of the kernel described by τ .We choose a smoothing scale of τ = 0.05 (see § 4.3 for the systematic uncertainty associate with the choice of smoothing scale).We then apply the GLOESSsmoothed LF with a Sobel filter kernel of width [-1,0,1] convolved with a weighting that suppresses Poisson noise (e.g., Madore et al. 2009;Hatt et al. 2017).The minimization routine explores the slope parameter space to maximize the Sobel response at the TRGB magnitude under the assumption that the largest Sobel response corresponds to the slope best characterizing the TRGB color dependence.We bound the slope values to ranges of (−1.7, −0.6) and (−2.8, −1.5) for F110W and F160W, respectively, use a tolerance of 1 × 10 −6 for the Nedler-Mead algorithm, and search for the Sobel peak response and corresponding TRGB magnitude over a range of ±0.2 mag about the expected location of the TRGB.
Third, we quantify the best-fit slope, zero point, and their uncertainties by running 1500 Monte Carlo (MC) simulations.Each MC simulation repeats step two for resampled magnitudes of the stars based on the uncertainties in the µ F814W measurements and on the photometric uncertainties determined from ASTs.The latter uncertainties are determined by finding the 2000 recovered artificial stars closest (e.g., d = magnitude 2 + color 2 ) to each observed star and calculating the standard deviation and median (i.e., the bias) of the difference in their recovered and input magnitudes.The bias ±0.1 mag about the TRGB magnitudes in F110W and F160W is < 0.001 mag.In addition, we resample the colors of the stars based on the photomet-ric uncertainties in the two NIR filters.We measure the final slope and zero point values at the 50th percentile of the MC output, and measure asymmetric uncertainties by taking the differences in the 16th and 84th percentiles from the median.We provide the following relationships in the WFC3/IR Vegamag system for the absolute magnitude calibrations in the F110W and F160W filters: Figure 6 presents the resulting best-fit slopes (top), zero points (middle), and corner plots (bottom) with uncertainties on the best-fit values from the MC resampling procedure overplotted on the stacked F110W (left) and F160W (right) CMDs.The best-fit slopes and zero points are indicated by black solid lines and statistical uncertainties from the MC simulations are shown as red shaded regions about the fits.The CMDs in the top panels are the original, unrectified data to demonstrate how the luminosity of the TRGB depends on the color of the stars, while the CMDs in the middle two panels have been rectified using the best-fit slopes.The corner plots show the joint probability of the fits from the 1500 MC trials for each filter and the 16 th , 50 th , and 84 th percentile fits to the 1d histograms as black dashed lines.We find that there is a modest correlation between the slopes and zero points in both filters.The slope and zero point values are shown in the bottom panels and are described in Eqs. 13 and 14.We again emphasize that these zero points are tied to the TRGB color of F110W-F160W = 0.92 mag.
Finally, we use the mathematical transformations derived in § 4. 1 (Eqs. 11 & 12) to check the empirically derived calibrations.In an ideal scenario it is guaranteed that the mathematical transformation and the empirically derived slopes and zero points are identical.However, if the slopes and zero points are not well fit, they can differ from each other.Our empirically fit M TRGB F160W,γ agrees exactly with the transformed zero point, while the empirically measured slope differs from the transformed slope by 0.01 mag/mag which is less than the uncertainty on the slope.Thus, the agreement is within uncertainties as expected.

Error Budget
We provide a summary of our total error budget and tabulate the sources of uncertainty in Table 5.The sources of uncertainty that we consider are: • Foreground extinction (sys) • Smoothing Scale (sys) • F814W QT color correction (sys) • NIR Zero Point (stat) • F814W Zero Point (stat and sys) For the foreground extinction systematic uncertainty we adopt a systemic uncertainty of 0.01 mag equal to 10% of the maximum extinction value in F814W.Based on the choice of GLOESS smoothing scale τ we determined a systematic uncertainty of 0.02 mag by re-running the MC simulation with smoothing scales of τ = 0.02 and τ = 0.08.We adopt an estimate of 0.01 mag systematic uncertainty for applying the F814W color correction applied to measure the TRGB in F814W.The statistical uncertainty on the NIR zero points is 0.04 mag in both filters.
Finally, we adopt the statistical uncertainty of 0.015 mag and the systematic uncertainty of 0.035 mag reported in Freedman (2021) for the F814W TRGB zero point adopted in anchoring distance scale for our NIR TRGB calibration.As discussed in §3.3, their systematic uncertainty is based on a weighted mean of F814W TRGB zero point calibrations tied to geometric anchors in the LMC, SMC, Galactic Globular Clusters, and the megamaser galaxy NGC 4258.The geometric distances to the LMC and SMC are determined from detached eclipsing binaries (Pietrzyński et al. 2013(Pietrzyński et al. , 2019;;Graczyk et al. 2020).The geometric distance to NGC 4258 is determined from its nuclear water megamaser sources (Humphreys et al. 2013;Riess et al. 2016).From Table 5, we see that the dominant contribution to the statistical uncertainty is the NIR zero point and the dominant contribution to the systematic uncertainty is the F814W magnitude of the TRGB.

Comparison of F814W and NIR Distances
In order to validate our NIR TRGB calibration, we re-measure the distances to each field in the individual NIR CMDs using our new calibration and compare them to the F814W distances.To measure the TRGB magnitudes in the F110W and F160W filters, we apply the same method used in the original F814W TRGB measurements.First, we apply spatial cuts to the photometry catalogs as described in §2.2.6 in order to reduce contamination from non-RGB populations.Second, we apply our color-based correction to each star within the color range 0.73 ≤ F110W-F160W ≤ 1.0 mag with the appropriate equation for F110W (see Eq. 13) or F160W (see Eq. 14).Third, we use the ML technique described in § 3 to the foreground-extinction corrected NIR photometry.
We use the zero point to determine the distance moduli in the NIR filters and compare the results to the distance moduli as determined in the F814W filter.Figure 7 presents the distance modulus comparison for the F160W and F814W results.The top panel shows the F160W distance moduli plotted as a function of the F814W distance moduli, with the top panel showing the distance moduli and the bottom panel showing the residuals between the distance moduli.We fit a linear model to the data and derive posterior probability distributions and the Bayesian evidence with the nested sampling Monte Carlo algorithm MLFriends (Buchner, 2014;2019) using the UltraNest 2 package (Buchner 2021).The right side of each panel shows the corner plots the UltraNest regression fits to the orange points in the top of the distance modulus panel.Table 6 presents both the F110W and F160W TRGB fits and resulting distance moduli and distances.In both comparisons, the UltraNest fits to the distance moduli result in approximately 1-to-1 relationships (regression slope fit = 1.01±0.04;see the right panel in Figure 7).Similarly, the standard deviation of the weighted mean based on distance modulus residuals is µ = 0.02 ± 0.02, confirming that the distance moduli agree.The F160W and F814W distance moduli are in excellent agreement for every field indicating a very robust calibration.

COMPARISON WITH PREVIOUS NIR TRGB CALIBRATIONS
Calibrations for the NIR TRGB have been reported in a number of studies based on HST photometry in WFC3/IR F110W and F160W, as well as based on theoretical stellar evolution models.These calibrations include color-based corrections, zero points, or both.In Figure 8, we compare the results of our F160W TRGB calibration with those from previous studies and pro-2 https://johannesbuchner.github.io/UltraNest/vide a summary of the calibrations in Table 7.Each calibration is plotted within the applicable color ranges (see Table 7 for the color ranges explored in each study).The individual calibrations are color-coded by the study from which they originate.The figure demonstrates the our calibration is in excellent agreement with the calibrations from the literature.We discuss individual differences between calibrations in greater detail below.Dalcanton et al. (2012a) used data from 23 nearby galaxies (D ≲ 4 Mpc) to determine a calibration from the TRGB only in the F160W filter.They anchored their calibration to TRGB distances that they measured in the F814W filter calibrated from theoretical PAR-SEC isochrones (Girardi et al. 2010;Bressan et al. 2012;Marigo et al. 2017).To measure the NIR TRGB, Dalcanton et al. (2012a) applied the edge-detection filter described in Méndez et al. (2002) to a Gaussian-smoothed luminosity function (Sakai et al. 1997) to each field.They first measured the apparent magnitude in each of their fields, then converted to absolute magnitudes by adopting their F814W TRGB-based distances as truth.However, there are known systematic differences in the predicated TRGB absolute magnitude depending on the adopted theoretical stellar evolutionary library which change both the shape and location of the TRGB as a function of color (e.g., Beaton et al. 2018;Durbin et al. 2020).Additionally, Dalcanton et al. (2012a) reported concerns that their analysis was impacted by the processing and calibration of WFC3/IR images, which at the time was in its early stages and has since improved.Durbin et al. (2020) re-reduced the data from Dalcanton et al. (2012a), including updating the anchoring F814W TRGB-based distances, and applied a mutliwavelength-covariance method for measuring the TRGB across several HST filters.Durbin et al. (2020) leverage a combination of observed photometry and synthetic photometry to determine their multiwavelength TRGB measurements.However, they report a systematic uncertainty on the order of ∼ 0.1 mag depending on the adopted stellar evolutionary library.They caution a conservative 10% systematic uncertainty included in any model-based absolute TRGB calibration.Wu et al. (2014) based their calibration on HST WFC3/IR data in the F110W and F160W filters for 32 fields; they combined the 23 galaxies covering 26 fields from Dalcanton et al. (2012a) with several M31 fields from the Panchromatic Hubble Andromeda Treasury (PHAT; Dalcanton et al. 2012b).They adopted ACS/WFC F814W filter TRGB-based distance moduli from the Extragalactic Distance Database (EDD; Tully et al. 2009) as their anchoring distances for all galaxies excluding M31.The EDD uses the extinc-  2009) to derive their own M31 distance modulus.They found a break in the TRGB magnitude color dependence at F110W-F160W ≈ 0.95.The change in slope at redder colors is predicted by stellar libraries and may be partially due to line-blanketing in the atmospheres of RGB stars where the metals preferentially absorb light at optical wavelengths and re-emit at NIR wavelengths.The line-blanketing effect is non-linear and may have a stronger effect at higher stellar metallicities leading to a non-linear increase in the color dependence of the TRGB luminosity function (Barker et al. 2004;Mager et al. 2008;Wu et al. 2014).
We find hints of a potential change in the slope in both F110W and F160W, albeit at a slightly bluer color than Wu et al. (2014).However, our data are not suited to find the slope change to a high degree of confidence, and we therefore assume a linear dependence across our full color range.In addition, our calibration does not extend to the colors and magnitudes for M31 where Wu et al. (2014) find that the TRGB magnitude in the NIR may no longer be linear with color (Barker et al. 2004;Mager et al. 2008).Future data sets which include more stars across a wider range of colors can help constrain this effect.An empirical calibration of the slope of the TRGB at NIR wavelengths that has a wider metallicity range has wide applications to both observers and to stellar modelers.
Over the color range in our calibration, the method that Wu et al. (2014) used to calibrate the TRGB in the F110W and F160W filters differs from the method developed in this study.They use the ML method described in Section 3.3 to measure the TRGB in each WFC3/IR filter prior to correcting for any color dependence, and derive a TRGB color from the median color of stars within a magnitude interval of 0.05 mag fainter than the measured TRGB.In contrast, we first determine our color-based correction from our combined photometry catalog, then measure the TRGB magnitudes from the rectified CMDs.These differences between our calibration and the Wu et al. (2014) calibration may help explain the offset as shown in Figure 8. Serenelli et al. (2017) report a theoretical calibration for the NIR TRGB-based on two different stellar evolutionary codes, BaSTI (Pietrinferni et al. 2004(Pietrinferni et al. , 2006) ) and GARSTEC (Weiss & Schlattl 2008).Serenelli et al. (2017) use bolometric corrections determined from theoretical and empirical calibrations to measure the color-magnitude relationship of the F110W and F160W TRGB.They report that uncertainties in the bolometric corrections make it challenging to generate an accurate reference set for TRGB absolute magnitudes.We find that our new calibration is in overall agreement with previous calibrations.

CONCLUSION
The TRGB in the NIR has the potential to expand the range over which efficient and precise distance measurements are possible.However, for the NIR TRGB to be useful as a precision distance indicator it must be carefully calibrated.We presented an entirely empirical calibration, relative to the F814W TRGB calibration, for the TRGB in the WFC3/IR F110W and F160W filters including color-based corrections and zero point calibrations accurate to ±0.043 (stat) ±0.043 (sys) mag in distance modulus based on the standard deviation of the mean of the residuals weighted by the individual distance modulus uncertainties, shown in Figure 7 as the gray shaded region in the residual plot.We derived high-fidelity cross-matched photometric catalogs based on observations of 8 galaxies and 12 fields in F606W and F814W filters, and WFC3/IR F110W and F160W filters.From the photometric catalogs we mea-sured TRGB-based distances in the F814W filter using the empirical QT calibration from (Jang & Lee 2017) and zero point from (Freedman 2021) to anchor our calibration on a uniform distance scale.
To measure the color-magnitude relationship in the F110W and F160W filters we generated a combined stellar catalog from the entire sample to provide the broad-0.6 0.7 0.8 0.9 1.0 1.1 F110W-F160W est metallicity baseline.From the combined stellar catalog we measure the slope, color pivot point, and zero point of the TRGB in the F110W and F160W filters.
We provide the results of our calibration with uncertainties in §4.2 & §4.3 (Eqs.13 & 14).Finally we verified the internal precision of the calibration to ∼ 1% in D based on a fit to the mean and weighted standard deviation of the mean in the distance modulus residuals µ F160W − µ F814W .These new calibrations enable distance measurements over a larger volume than is accessible with the ACS F814W TRGB method.
The calibration presented in this study is robustly measured in the WFC3/IR F110W and F160W filters within the color range 0.73 ≤ F110W-F160W ≤ 1.0 mag.We emphasize that, similarly to in the F814W filter, the TRGB in the F110W and F160W filters is most easily identified when minimal non-RGB stellar populations are present in the CMD.These include AGB and RHeB stars that can mask the location of the TRGB in a CMD, especially since the color baseline of the RGB is narrow.Careful choice of observing fields can help mitigate the masking effect from non-RGB populations.In particular, imaging the outer stellar fields of galaxies may help to reduce the impact of non-RGB populations, since we expect older RGB stars to dominate in this spatial regime.The next phase of the NIR TRGB cal-ibration is to measure color-based corrections and zero points in the JWST filters on the NIRCam and NIRISS instruments as part of JWST-GO-1638 (PI McQuinn).

Figure 1b .
Figure 1b.Clockwise from the top left: Antlia Dwarf, UGC8833, UGC9128, and UGC9240.Shown are color images using GALEX near-ultraviolet (NUV; blue) and far-ultraviolet (FUV; red) filters, and an average of the NUV and FUV filters (green), with the HST ACS and WFC3/IR instrument footprints indicated in green and white, respectively.Images are north up, east left.GALEX images were obtained courtesy of NASA/JPL-Caltech.

Figure 2 .
Figure 2. Example of crowding (top) and sharpness (bottom) values as a function of magnitude for NGC2403-2 from the photometry output of DOLPHOT.From left to right: F606W, F814W, F110W, and F160W filters.Point sources included in our final photometry catalogs are in black and excluded point sources are in gray.Dashed horizontal lines bracket the accepted range of values.
Figure 3. Foreground-extinction corrected optical CMDs with spatial cuts applied to the photometry.From top to bottom in the first and second rows: M81, NGC253, NGC2403, and NGC300.The first and third columns show field 1, while the second and fourth columns show field 2. The third row shows the 4 dwarf galaxies: Antlia Dwarf, UGC9128, UGC8833, and UGC9240.Each panel consists of a CMD using the F606W and F814W filters.Representative photometric uncertainties per magnitude bin are shown.Boxed regions (orange) show the CMD selections to isolate the RGB sequence.We empirically determine the blue edge of box such that it follows the blue edge of the RGB, while the red edge encompasses all sources redward of the blue edge.Point sources inside (outside) the selection box are shown as black (gray) points.The locations of the TRGB magnitudes as identified from the F814W filter with the ML TRGB fit method are shown as horizontal lines (red) with uncertainties indicated as shaded light blue regions around the fit (see § 3 for more details).

Figure 4 .
Figure 4. Foreground-extinction corrected NIR CMDs with the F110W and F160W filters.The ordering of CMDs is the same as in Figure 3.The star counts in each overlapping ACS-WFC3/IR region are provided in the legends.Representative photometric uncertainties per magnitude bin are shown.

Figure 5 .
Figure 5. Example of spatial cuts that were carefully applied to specific fields.Left: M81-1 ACS color image with point sources from our final photometry catalog overlaid.Stars in the upper left corner of the image (green points) are excluded from the final photometry catalog used to fit the TRGB.The remaining stars (orange points) are kept after our spatial cuts.Right: Extinction corrected optical CMD of the stars from the spatially excluded region (green) and spatially included region (orange).The spatial cut removes most main-sequence stars and a significant fraction of the HeB sequences, while retaining the majority of the RGB stars.

Figure 6 .
Figure 6.The final results from the NIR TRGB calibration MC simulations.The left and right sides of the figure show the results for F110W and F160W, respectively.(Top and middle panels) NIR CMDs with the best-fit calibration values (black solid lines) and uncertainties (red shaded region) overplotted.The CMDs include photometry from all galaxies after placing the catalogs on an absolute magnitude scale.The top two panels panels present the unrectified CMDs with slope fits of βF110W = −1.15mag/mag and βF160W = −2.16mag/mag.The middle two panels present the rectified CMDs with zero point fits of M W F C3 F110W = −4.87mag and M W F C3 F160W = −5.79mag .Uncertainties are measured through 1500 MC simulations by simultaneously resampling the statistical uncertainties from µF814W and the photometric errors based on ASTs (see Section 4.2 for more details).(Bottom panel) The join probability distributions of the MC simulations for the slope and zero point fits.The 16 th , 50 th and 84 th percentile fits to the 1d histograms are shown as black dashed lines.The corner plot demonstrates a modest correlation between the slope and zero point fits and is propagated to the statistical uncertainties on the zero point.

Figure 7 .
Figure7.We present a comparison of the distance moduli based on the F814W TRGB and NIR F160W TRGB.Left: The top panel shows the distance modulus comparison (filled green circles).The orange dashed line and grey shaded region shows the linear regression fit to the green points, and the solid black line is the 1-to-1 line.We fit the data with UltraNest (Buchner 2021).The bottom panel shows the residuals in DM ∆µ = µ0,F160W −µ0,F814W.The positive residuals indicate that the NIR modulus is fainter, thus the measured distance is further away, than the F814W-based distance.The black dashed line marks ∆µ = 0, and the orange dashed line and gray shaded region are the mean and standard deviation of the mean weighted by the uncertainties on the data.The residual plot demonstrates that the distance moduli are in excellent agreement within the uncertainties based on both the weighted mean and standard deviation of the weighted mean of the residuals (μF160W ≈ −0.02 ± 0.02 (stat)).As a fractional uncertainty in distance, this represents a precision of ∼ 1%.Right: Corner plot for the fit parameters (slope, offset, and intrinsic scatter) from the ML fit to the data in the distance moduli panel.The result of the fit is consistent with a 1-to-1 relationship between µ0,F160W and µ0,F814W.

DurbinFigure 8 .
Figure 8. F160W TRGB color-magnitude relationships.Left: Calibrations from this work (gray) as well as other studies from the literature with uncertainties where available.Note that the uncertainties about the Wu et al. (2014) calibration are shown only for the zero point.The mismatch in the Wu et al. (2014) calibration where the slope changes at a color F110W-F160W= 0.95 is due to differences in the uncertainties as a function of color.Right: The same color-magnitude relationships in the left panel now presented relative to the Durbin et al. (2020) calibration.The calibrations agree within the uncertainties over a narrow range in color with the Wu et al. (2014) calibration deviating the most significantly at the reddest colors.

Table 1 .
Galaxy Sample and Summary of the HST Observations

Table 3 .
Summary of Distances Measured from F814W LF

Table 5 .
Summary of the Error Budget for the NIR TRGB Calibration

Table 6 .
Summary of Distances Measured from NIR LF