SN 2023ixf in Messier 101: A Variable Red Supergiant as the Progenitor Candidate to a Type II Supernova

All hydrogen-rich supernovae (SN II) with directly identified progenitor stars have been interpreted to come from systems with initial mass <20 M_⊙ (Smartt 2015). With the exception of the blue supergiant progenitor of the peculiar SN II 1987A (Hillebrandt et al. 1987; Arnett et al. 1989), the yellow supergiant progenitor stars of hydrogen-poor SNe IIb (e.g., Aldering et al. 1994), and the luminous blue variable (LBV) progenitor stars to SNe IIn (e.g., Gal-Yam & Leonard 2009), all of these systems are red supergiants (RSGs). These stars have massive, extended, hydrogen envelopes and make up the majority of directly identified progenitor stars to core-collapse SNe (SNe 2003gd, 2004A, 2004et, 2005cs, 2006my, 2008bk, 2009hd, 2009kr, 2009md, 2012A, 2012aw, 2012ec, 2016cok, 2017eaw, 2018aoq, 2020jfo, and 2022acko; Smartt et al. 2004; Maund & Smartt 2009; Elias-Rosa et al. 2010; Fraser et al. 2010; Crockett et al. 2011; Fraser et al. 2011; Kochanek et al. 2012; Maund et al. 2013; Tomasella et al. 2013; Fraser et al. 2014; Maund et al. 2014; Kochanek et al. 2017; Kilpatrick & Foley 2018; O'Neill et al. 2019; Rui et al. 2019; Van Dyk et al. 2019; Sollerman et al. 2021; Van Dyk et al. 2023a, 2023b). The lack of >20 M_⊙ RSG progenitor stars of SNe II despite the fact that they make up ≈15% of RSGs following a Salpeter initial mass function and RSGs with $\mathrm{log}(L/{L}_{\odot })\gt 5.2$ are observed in the LMC, M31, and M33 (Drout et al. 2012; Neugent et al. 2020; Neugent 2021a, 2021b) has been noted as the "RSG problem" (Smartt et al. 2009; Smartt 2009, although see also Davies & Beasor 2018).

Theoretically, massive RSGs are predicted to have compact oxygen cores, and many of them may collapse directly to black holes as "failed SNe," leading to a paucity of high-mass counterparts to SNe II (Sukhbold et al. 2016). This scenario broadly agrees with the light curves and nucleosynthetic yields of SNe II, which also favor lower-mass progenitor stars (Brown et al. 2013; Müller et al. 2017; Morozova et al. 2018) as well as direct evidence for a high-mass RSG in NGC 6946 whose optical counterpart disappeared (Adams et al. 2017; see also Neustadt et al. 2021a and Byrne & Fraser 2022). This source also left behind a weak infrared (IR) transient consistent with expectations for mass ejection in failed SNe (Lovegrove & Woosley 2013; Piro 2013; Fernández et al. 2018). Long time baseline follow up of nearby galaxies with deep, high-resolution imaging can constrain the fraction of disappearing stars, such as the estimate by the "Survey for Nothing" that ≈16% of massive stars produce failed SNe in Neustadt et al. (2021b), close to the value required by an upper mass threshold for successful explosions of >20 M_⊙ stars. Simultaneously constraining the fractions and mass distributions of failed SNe and SN II progenitor stars is therefore a powerful tool for probing massive star structure and the latest stages of stellar evolution.

SN II progenitor stars also exhibit a wide range of circumstellar densities in their immediate vicinity (<10¹⁵ cm) as implied by flash spectroscopy (Gal-Yam et al. 2014; Khazov et al. 2016; Yaron et al. 2017; Terreran et al. 2022; Tinyanont et al. 2022), early photometric evolution (Morozova et al. 2017, 2018), as well as evidence for preexplosion variability and eruptions (Kilpatrick & Foley 2018; Jacobson-Galán et al. 2022). These features may have significant implications for the interpretation of their preexplosion counterparts in the absence of multiband, multiepoch imaging. For example, the vast majority of SN II preexplosion counterparts are identified in F814W imaging from the Hubble Space Telescope (HST; e.g., Smartt 2009; Davies & Beasor 2018). This filter is blueward of the peak of RSG spectral energy distributions (SEDs) and may be significantly impacted by circumstellar extinction in the presence of a dusty shell. Moreover, many RSGs exhibit well-known modes of variability (Stothers 1969; Jurcevic et al. 2000; Guo & Li 2002; Yang & Jiang 2011; Soraisam et al. 2018) that may become even more extreme as they approach core collapse (Yoon & Cantiello 2010; Fuller 2017; Davies et al. 2022), although SN 2016cok is a counterexample whose progenitor star exhibited very little variability (Kochanek et al. 2017). Without multiepoch imaging in which their average luminosities can be estimated, interpretation of photometry for SN II preexplosion counterparts is complicated by large systematic uncertainties.

Here we present preexplosion imaging to the nearby SN II 2023ixf discovered in Messier 101 (M101) on 2023 May 19 (Itagaki 2023). These data cover ultraviolet to mid-IR bands from 1999–2019. We demonstrate that there is a single credible progenitor candidate to SN 2023ixf and estimate its luminosity, temperature, and initial stellar mass as well as its variability and total circumstellar material (CSM) inferred from a significant mid-IR excess. We find it was significantly variable in the mid-IR and compare that timescale with well-observed RSGs. We summarize the total data set and our reduction procedure in Section 2 and the analysis and modeling of those data in Section 3. Finally, we discuss the broader implications of this progenitor candidate in Section 4 and conclude in Section 5.

We assume a line-of-sight extinction through the Milky Way of A_V = 0.025 mag (assuming R_V = 3.1, this is E(B − V) = 0.008 mag) from Schlafly & Finkbeiner (2011). We also adopt the latest Cepheid distance to M101 of 6.85 ± 0.15 Mpc from Riess et al. (2022). Finally, throughout this paper we reference W. V. Jacobson-Galán et al. (2023, in preparation), who demonstrate that SN 2023ixf appears to be a normal type II SN with broad lines of hydrogen. We also assume a host reddening to SN 2023ixf of E(B − V) = 0.033 mag from W. V. Jacobson-Galán et al. (2023, in preparation), derived from Na i D line absorption in optical spectra of this event. Given the small value for this line-of-sight reddening, it does not significantly impact our results and we adopt a total-to-selective extinction ratio in the host galaxy of R_V = 3.1 (implying A_V = 0.10 mag), however we acknowledge that this could range from R_V = 2–6 (implying A_V = 0.07–0.20 mag).

2.1. HST

The site of SN 2023ixf was observed with HST's WFPC2, Advanced Camera for Surveys (ACS), and WFC3 over seven epochs from 1999 March 23 to 2018 March 30, or 24.2 to 5.1 yr before discovery (Table 1). Following methods described in Kilpatrick et al. (2022) and Kilpatrick et al. (2021), we used the custom python-based pipeline hst123 ²⁰ to download, align, and drizzle all HST imaging (for details, see Hack et al. 2021), and perform photometry in dolphot (Dolphin 2016). We used recommended dolphot settings for each imager as described in the respective manual. ²¹

Table 1. HST Photometry of the SN 2023ixf Progenitor Candidate

MJD	Instrument	Filter	m	σm
			(mag)	(mag)
51260.9786	WFPC2	F656N	>23.433	⋯
51261.0390	WFPC2	F675W	26.422	0.230
51261.1120	WFPC2	F547M	>26.273	⋯
51345.9897	WFPC2	F656N	>23.776	⋯
51346.0529	WFPC2	F547M	>26.416	⋯
52593.9933	ACS/WFC	F435W	>27.393	⋯
52594.0096	ACS/WFC	F555W	>27.099	⋯
52594.0215	ACS/WFC	F814W	24.881	0.059
52878.3224	WFPC2	F336W	>27.025	⋯
53045.0069	ACS/WFC	F658N	25.332	0.284
56735.8571	WFC3/UVIS	F673N	>24.629	⋯
58207.4384	ACS/WFC	F658N	>25.488	⋯
58207.4561	ACS/WFC	F435W	>27.799	⋯

Note. All magnitudes are in the AB system.

A machine-readable version of the table is available.

Download table as:  DataTypeset image

The final stacked imaging of M101 observed in 2002 by ACS is shown in Figure 1 as an RGB image (F814W, F555W, and F435W). We also show each image in which we obtain a detection of a counterpart at the explosion site of SN 2023ixf, which includes WFPC2 F675W and ACS F658N imaging. In addition, we have deep constraints in the bluer bands F336, F435W, and F555W, which we consider in the context of a binary companion below.

Zoom In Zoom Out Reset image size

Within 02 of the reported position of SN 2023ixf, there are two sources detected in F814W. This is clearly seen in Figure 1 where a counterpart is located at the site of SN 2023ixf in the ACS imaging and appears as a blended source in the F814W panel. The brighter source has m_F814W = 24.881 ± 0.059 mag ²² that we refer to as "Source A" (we note that this is the same source identified in Pledger & Shara 2023) and is blended with the fainter "Source B" approximately 01 (2.0 ACS/WFC pixels) to the northeast with m_F814W = 25.955 ± 0.125 mag. Below we consider which, if either, of these sources may be the preexplosion counterpart to SN 2023ixf and the extent to which any blended emission from other sources may contaminate the photometry of that source in other bands with poorer resolution.

2.2. Spitzer Space Telescope/IRAC

2.3. Ground-based IR Imaging

2.4. GMOS Imaging of SN 2023ixf

Gemini-North/GMOS obtained a series of 10 × 1 s images and 3 × 60 s images in the gri bands on 2023 June 3 and 5. We obtained all such imaging from the Gemini Data Archive, ²⁶ and following standard procedures in astropy, we removed the bias from these frames using the overscan correction. We then calibrated each frame using DoPhot photometry (Schechter et al. 1993) and Pan-STARRS standard stars in each frame of the GMOS images (Flewelling et al. 2020). To obtain the deepest possible image of the field surrounding SN 2023ixf without saturating the SN position, we masked saturated pixels from the SN and stacked all images in swarp (Bertin 2010) into a single g+r+i frame, weighting the individual frames by the inverse variance of the sky pixels across the individual bands. The final image centered on the position of SN 2023ixf is shown in Figure 2.

Zoom In Zoom Out Reset image size

3.1. Aligning Pre- and Postexplosion Imaging

3.2. The SED of the SN 2023ixf Progenitor System

3.2.1. Single Blackbody Fit

Assuming that the SN 2023ixf preexplosion counterpart is dominated by a single SED from its progenitor star and with no variability between each epoch (though see the Spitzer analysis in Section 3.3), we can model the nature of this source from the ultraviolet to mid-IR. We initially adopt a simple blackbody spectrum and derive its temperature and luminosity using a Markov Chain Monte Carlo (MCMC) approach in the python-based package emcee. Assuming the distance, Milky Way extinction, and host extinction given above, we derive the in-band magnitudes for a blackbody spectrum of a given temperature and luminosity using the filter transmission functions for each space- and ground-based bandpass and using pysynphot. Following methods in Kilpatrick et al. (2021), we fit a blackbody model by sampling the posterior distribution over the range of model parameters, and report their posterior means and standard deviations.

Following this method, we find that the SN 2023ixf progenitor candidate is consistent with a $\mathrm{log}(L/{L}_{\odot })={4.73}_{-0.08}^{+0.07}$ and T_eff = 1640 ± 20 K blackbody as shown in Figure 3. This implied photospheric temperature is extremely low, even for the latest M supergiant spectral types (i.e., the coolest RSGs have T_eff = 3400–3500 K; Levesque et al. 2006; Davies et al. 2013; Davies & Beasor 2018). The photospheric radius implied by T_eff = 1640 K (≈2900 R_⊙) suggests that we are seeing material extended well beyond the envelope of an RSG such as a shell of CSM in the local environment around the SN 2023ixf progenitor star. We infer that the photosphere in the IR is dominated by a component of host dust, and instead turn to a multicomponent SED fit below.

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

3.2.2. MARCS RSG and Circumstellar Dust Fit

To fit the counterpart with a more realistic optical to mid-IR SED, we use a combined RSG and dusty CSM spectrum initially presented in Kilpatrick & Foley (2018) and based on DUSTY radiative transfer models (from Kochanek et al. 2012; see also Ivezic & Elitzur 1997; Ivezic et al. 1999; Elitzur & Ivezić 2001). This model uses a MARCS RSG spectrum of an arbitrary temperature and luminosity (see Gustafsson et al. 2008 for details), which is reprocessed through a shell of graphitic dust at a given temperature and mass (i.e., similar to carbon-bearing species around massive RSGs; e.g., Royer et al. 2010). In general, we fit for the RSG luminosity, stellar temperature, dust temperature, and V-band optical depth through the mass of dust. These assumptions yield a mass, luminosity, and radius for the dust assuming an r⁻² density profile (see Kochanek et al. 2012; Kilpatrick & Foley 2018 for more details). Assuming a dust-to-gas ratio and wind speed, we can then derive the total mass of CSM and mass-loss rate, which we give below.

From this model, we find that the effective temperature of the counterpart is more in line with known RSGs at ${T}_{\mathrm{eff}}={3920}_{-160}^{+200}$ K while the overall luminosity remains the same at $\mathrm{log}(L/{L}_{\odot })=4.74\pm 0.07$. In order to fit the IR excess observed in the Spitzer bands, we require a circumstellar shell of dust with a V-band optical depth of τ_V = 5.8 ± 0.2 (corresponding to A_V = 4.6 ± 0.2 mag) and an effective temperature of T_dust = 880 ± 40 K. This shell would have an effective radius of ${8600}_{-800}^{+900}$ R_⊙ and a total dust mass of ${5.0}_{-0.8}^{+1.1}\times {10}^{-7}$ M_⊙, or a total mass of 5 × 10⁻⁵ M_⊙ assuming a dust-to-gas ratio of 0.01 in the circumstellar environment (consistent with ratios in the environments of SNe from Fox et al. 2010).

Assuming it was being produced by a constant wind with an r⁻² profile with a velocity of v_wind, the implied mass-loss rate is $\dot{M}/({v}_{\mathrm{wind}})=1.3\pm 0.1\times {10}^{-6}\,{M}_{\odot }\,{\mathrm{yr}}^{-1}/(50\,\mathrm{km}\,{{\rm{s}}}^{-1})$. We assume a velocity of 50 km s⁻¹ for consistency with W. V. Jacobson-Galán et al. (2023, in preparation), noting that this value is similar to other RSGs (e.g., 30–50 km s⁻¹ for VY CMa or NML Cyg; Knapp et al. 1982; Decin et al. 2006), however high-resolution spectroscopy of unshocked material from early in the evolution of SN 2023ixf can more precisely constrain this value.

In the fits above, we do not include the ACS F658N detection despite it being spatially coincident with Source A, as this counterpart is likely dominated by Hα emission that we do not include in our model. However, assuming that the emission in this filter contains Hα and continuum emission from an RSG, we estimate that the total Hα flux density is 1.3 × 10⁻¹⁹ erg s⁻¹ cm⁻² Å⁻¹ (corrected for host and Milky Way extinction). This corresponds to a total Hα luminosity of 4.9 × 10³⁹ erg s⁻¹ or 1300 L_⊙. This value far exceeds the expectations for Hα emission in massive RSGs (e.g., 1 L_⊙ for VY CMa in Smith et al. 2001), and so may be unassociated with the progenitor star.

RSG binaries such as VV Cephei exhibit high Hα luminosities due to mass transfer from the RSG to an early-type binary companion (Wright 1977). The extreme variability in this line suggests a change in the mass transfer rate that is correlated with the properties of the wide binary orbit for this system. However, for VV Cephei, the equivalent width of ≈7–25 Å during the most active phases of this source implies a luminosity of L_Hα ≈ 10–50 L_⊙ (Pollmann 2005). We cannot rule out a wide binary with our optical HST limits (Section 3.2.3), which could also explain the Hα detection for an even more extreme mass transfer rate.

Assuming that the underlying star implied by our RSG model is a single source with $\mathrm{log}(L/{L}_{\odot })=4.74\pm 0.07$, we consider its initial mass by comparing to MIST (Choi et al. 2016), STARS (Eldridge & Tout 2004), Geneva (Ekström et al. 2012), and KEPLER (Woosley & Heger 2007) models. We show the derived luminosity in comparison to the MIST models in Figure 4. All models assume solar metallicity, and either a rotating or nonrotating star. In general, the final luminosity of a model SN progenitor depends on the He core luminosity, which is higher in models including rotation and overshooting. For a MIST model assuming a star at solar metallicity, we find an initial mass of 11 ± 1 M_⊙. The STARS models indicate M_ZAMS = 11 ± 1 M_⊙, while the the Geneva rotating models produce a star of similar final luminosity at M_ZAMS = 11−12 M_⊙. Finally, the KEPLER (nonrotating) models indicate M_ZAMS = 12 ± 1. Thus the star could feasibly come from a system ranging from 10–13 M_⊙. Second dredge-up in 6–9 M_⊙ stars can increase the final luminosity substantially (e.g., Eldridge et al. 2007; Jones et al. 2013), with the stars along the asymptotic giant branch (AGB) ending up more luminous and with cooler photospheric temperatures than their more massive M-type counterparts. While our DUSTY SED models favor a higher T_eff than those of AGB stars (typically T_eff ∼ 3100 K) it is possible the progenitor is a cooler AGB-type star. Measurements of the nickel mass created in the explosion and the oxygen mass ejected (both after 100–200 days) will provide interesting constraints on the core mass and explosion mechanism.

Finally, we examine the consistency of our derived mass-loss rate and initial mass with prescriptions from Beasor et al. (2020). Applying their luminosity-dependent mass-loss rates, we find that for RSGs of $\mathrm{log}(L/{L}_{\odot })=4.74$ they predict $\dot{M}=0.7$–4.2 × 10⁻⁶ M_⊙ yr⁻¹, ²⁸ which is in close agreement with our inferred value. Similarly, applying their initial mass (for 10–12 M_⊙) and luminosity-dependent parameterization, ²⁹ we derive 0.4–1.1 × 10⁻⁶ M_⊙ yr⁻¹. These values are close to our derived mass-loss rate of 1.3 × 10⁻⁶ M_⊙ yr⁻¹, especially considering the model-dependent nature of these relations, the uncertain wind velocity, and the uncertainties in the model fitting parameters.

3.2.3. Constraints on a Binary Companion from Optical Limits

We also consider the possibility that the progenitor star evolved in a binary and exploded as the primary star in that system. Comparing our photometry to BPASS v2.2.1 binary star models (Eldridge et al. 2017), we examine all systems for which the combined flux from the primary and secondary at the time the primary explodes is fainter than our limiting magnitudes. We emphasize that these models do not include circumstellar extinction or predictions for the mid-IR luminosity, and we only use the limits from our bluer bands where we predict the primary star to be faint in order to constrain the presence of a companion star.

Although we examine all bands contained in the BPASS models for which we have deep limits (F336W, F435W, and F555W), our most constraining limit comes from ACS F555W with m_F555W > 27.1 mag, corresponding to M_F555W > −2.2 mag with no additional extinction from the CSM. This could be the case for a companion at wide separations with minimal additional extinction. For BPASS models with solar metallicity, this limit excludes any systems with a secondary star with >6.4 M_⊙. In scenarios where the V-band circumstellar extinction (A_V = 4.4 mag, implying M_F555W > −6.6 mag) is taken into account, virtually all BPASS models are consistent with our limits. A close binary therefore remains a possibility for SN 2023ixf, whose presence could better be constrained with deep optical imaging after the SN fades.

3.3. Preexplosion Variability of the Progenitor Candidate and Implications for Mass Loss

The SN 2023ixf progenitor star exhibited extreme variability in the Spitzer/IRAC bands several years before explosion (Figure 5). We show the light curve of the counterpart at 3.6 and 4.5 μm from ∼2012–2020 in Figure 5. For comparison, we also show the 3.9 μm light curve of α Orionis (Betelgeuse) from 2017 January 1 to 2022 December 12 (Taniguchi et al. 2022), but shifted in time by 4 yr so it overlaps with that of the SN 2023ixf counterpart and scaled to its average flux at ∼23 μJy. In the SED analysis above, we emphasize that we accounted for intrinsic variations in the Spitzer bands by including the standard deviation across the light curves in the average Channel 1 and 2 fluxes.

Zoom In Zoom Out Reset image size

The SN 2023ixf progenitor star exhibited significant mid-IR variability with an average of 22.44 and 23.99 μJy and peak-to-peak variability of 15.64 and 17.86 μJy at 3.6 and 4.5 μm, respectively (roughly 0.8 mag or 70% in both bands). These extreme variations appear correlated in the two Spitzer bands with approximately the same overall magnitude, which would be consistent with a mode of variability where the visible photosphere expands and contracts with at most small variations in effective temperature. We also note that these variations are similar in amplitude to the high-luminosity end of large-amplitude, cool pulsators observed in the LMC (O'Grady et al. 2020, 2023).

Moreover, the light curve from 2012–2020 exhibits as quasi-sinusoidal variation with a timescale of roughly 2.8 yr (1000 day). We infer this timescale via a Fourier transform of the Spitzer/IRAC photometry, although the sampling of the light curve and the fact that we only observe peak-to-peak variations over ∼2.5 cycles in Figure 5 suggests it is uncertain and could vary anywhere from 2.6–3.0 yr. Combined, this evidence is similar to the κ-mechanism pulsations in Betelgeuse, which are the primary mode of variability in that star and are driven by changes in the atmospheric opacity (see, e.g., Li & Gong 1994; Heger et al. 1997). Paxton et al. (2013) observe these modes with timescales of 1–8 yr directly in MESA models of RSGs where the structure of the star is resolved with sufficiently high time resolution, also in close agreement with the simulations in Yoon & Cantiello (2010) and our inferred timescale. These pulsations drive expansion and contraction in the atmosphere at a nearly constant temperature (e.g., Levesque & Massey 2020), resulting in overall changes to the luminosity from Betelgeuse.

In addition, the significant variability correlated across both bands supports the conclusion that the Spitzer counterpart is dominated by a single source. As this variability is extreme in the IR even for a single RSG, the flux is unlikely to contain significant emission from a second source (e.g., Source B) compared with the minimum flux level of our light curve in Figure 5.

Assuming this mechanism is responsible for the variability in the SN 2023ixf counterpart, the required changes in luminosity are ≈50% larger than in Betelgeuse. Our derived mass-loss rate is comparable with Betelgeuse (which has a rate from 0.2–2 × 10⁻⁶ M_⊙ yr⁻¹; Dolan et al. 2016), which matches expectations for the κ-mechanism driving strong mass loss with large variations in luminosity. However, our estimate precludes a "superwind" generated in the CSM (e.g., Yoon & Cantiello 2010; Davies et al. 2022) up to the point where our data cut off 3.6 yr before explosion. The preexplosion mass-loss rate could be enhanced if the star was significantly more active during the final 3 yr before core collapse, which is predicted by W. V. Jacobson-Galán et al. (2023, in preparation).

Current estimates of the maximum luminosity of SN II progenitor stars are dominated by a handful of direct counterpart detections in the literature (e.g., Davies & Beasor 2018; Kochanek 2020), and the highest-luminosity stars in those samples contribute significant weight to statistical analyses of the distribution from which they are drawn.

Assuming the SN 2023ixf candidate counterpart is dominated by its progenitor star emission in the HST, Spitzer, and ground-based imaging, SN 2023ixf provides the best example to date of the optical to mid-IR SED of an SN II progenitor star (comparable to SN 2017eaw; Kilpatrick & Foley 2018; Rui et al. 2019; Van Dyk et al. 2019), but with a low foreground host-galaxy extinction and precise distance. The significant variability and large implied circumstellar extinction support the presence of such features in other SN II progenitor stars given the lack of multiepoch, multiband imaging in which they can be analyzed. The vast majority of directly detected SN II progenitor stars in recent analyses (Smartt 2015; Davies & Beasor 2018; Kochanek 2020) were identified from F814W imaging. Our analysis of SN 2023ixf demonstrate that it has extreme bolometric corrections when accounting for its mid-IR luminosity. This has been explored in the past as a possible explanation to the RSG problem (Walmswell & Eldridge 2012), but our dust mass estimates imply a much larger extinction than has been previously explored for SN II progenitor stars. However even with the significant stellar flux in the near- and mid-IR, the total preexplosion luminosity indicates an initial mass significantly below 20 M_⊙. Future systems with Spitzer or JWST detections similar to SN 2023ixf can constrain the distribution of that extinction from their IR excess.

We have presented direct imaging showing that there is a credible progenitor candidate to the type II SN 2023ixf in M101 at 6.85 Mpc. In summary, this imaging demonstrates:

• 1.  
  The candidate progenitor star to SN 2023ixf is most consistent with a $\mathrm{log}(L/{L}_{\odot })=4.74\pm 0.07$ RSG with an effective temperature of ${T}_{\mathrm{eff}}={3920}_{-160}^{+200}$ K. Following single-star evolutionary tracks, this would place the progenitor star's initial mass at 11 ± 2 M_⊙, placing it within the range of other low- to moderate-mass RSG progenitors to SNe II (e.g., Smartt 2015).
• 2.  
  Modeling of the mid-IR SED from this counterpart suggests that it was enshrouded in a dusty shell of CSM similar to SN 2017eaw (Kilpatrick & Foley 2018). The implied mass-loss rate for a wind that could produce this shell divided by its wind speed is $\dot{M}/{v}_{\mathrm{wind}}=1.3\pm 0.1\times {10}^{-6}\,{M}_{\odot }\,{\mathrm{yr}}^{-1}/(50\,\mathrm{km}\,{{\rm{s}}}^{-1})$. This is comparable to more massive stars such as Betelgeuse but low compared to a "superwind" or mass-loss rates from immediately before explosion inferred in W. V. Jacobson-Galán et al. (2023, in preparation).
• 3.  
  Spitzer/IRAC photometry exhibits significant preexplosion variability that is correlated in both bands. We also see evidence for a 2.8 yr (1000 day) timescale in this variability, similar to but generally stronger in amplitude than pulsations in other RSGs driven by opacity changes in their atmosphere (i.e., the κ mechanism; Li & Gong 1994; Heger et al. 1997; Yoon & Cantiello 2010; Paxton et al. 2013).

Future studies of the progenitor stars to SNe II will greatly benefit from multiband, multiepoch imaging of resolved stellar populations using deep optical and IR surveys of nearby galaxies, such as those by the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope (Ivezić et al. 2019; Gezari et al. 2022). To understand these stars better, their preexplosion evolution, and the exact cause of the RSG problem, only detailed SEDs and light curves from the optical to mid-IR, such as those that we present for the SN 2023ixf counterpart, can shed light on the pathways through which SN II progenitor stars evolve and explode.

We thank Divakara Mayya for providing helpful comments on our manuscript and Daisuke Taniguchi for providing comments and the Betelgeuse mid-IR light curve.

C.D.K. is partly supported by a CIERA postdoctoral fellowship and acknowledges support from a NASA grant for HST-GO-16136. M.R.D. acknowledges support from the NSERC through grant RGPIN-2019-06186, the Canada Research Chairs Program, the Canadian Institute for Advanced Research (CIFAR), and the Dunlap Institute at the University of Toronto. A.G. is supported by the National Science Foundation Graduate Research Fellowship Program under grant No. DGE1746047. A.G. also acknowledges funding from the Center for Astrophysical Surveys Fellowship at UIUC/NCSA and the Illinois Distinguished Fellowship. K.S.M. acknowledges funding from the European Unions Horizon 2020 research and innovation program under ERC Grant Agreement No. 101002652 and Marie Sklodowska-Curie Grant Agreement No. 873089. C.G. is supported by a VILLUM FONDEN Young Investigator Grant (project number 25501). J.H. was supported by a VILLUM FONDEN Investigator grant (project number 16599). S.J.S. acknowledges funding from STFC grants ST/X006506/1 and ST/T000198/1.

The UCSC team is supported in part by NASA grants NNG17PX03C and 80NSSC22K1518, NSF grant AST-1815935, the Gordon & Betty Moore Foundation, the Heising-Simons Foundation, and by a fellowship from the David and Lucile Packard Foundation to R.J.F.

The Young Supernova Experiment (YSE) and its research infrastructure is supported by the European Research Council under the European Union's Horizon 2020 research and innovation program (ERC Grant Agreement 101002652, PI: K. Mandel), the Heising-Simons Foundation (2018-0913, PI: R. Foley; 2018-0911, PI: R. Margutti), NASA (NNG17PX03C, PI: R. Foley), NSF (AST-1720756 and AST-1815935, PI: R. Foley; AST-1909796 and AST-1944985, PI: R. Margutti), the David & Lucille Packard Foundation (PI: R. Foley), VILLUM FONDEN (project 16599, PI: J. Hjorth), and the Center for AstroPhysical Surveys (CAPS) at the National Center for Supercomputing Applications (NCSA) and the University of Illinois Urbana-Champaign.

The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen's University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, STScI, NASA under grant NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, NSF grant AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation.

YSE-PZ (Coulter et al. 2023) was developed by the UC Santa Cruz Transients Team with support from NASA grants NNG17PX03C, 80NSSC19K1386, and 80NSSC20K0953; NSF grants AST-1518052, AST-1815935, and AST-1911206; the Gordon & Betty Moore Foundation; the Heising-Simons Foundation; a fellowship from the David and Lucile Packard Foundation to R.J.F.; Gordon and Betty Moore Foundation postdoctoral fellowships and a NASA Einstein fellowship, as administered through the NASA Hubble Fellowship program, and grant HST-HF2-51462.001 to D.O.J.; and a National Science Foundation Graduate Research Fellowship, administered through grant No. DGE-1339067, to D.A.C.

This research is based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 526555. These observations are associated with programs 6829 (PI: Chu), 9490 (PI: Kuntz), 9720 (PI: Chandar), and 13361 (PI: Blair). This work is based in part on archival data obtained with the Spitzer Space Telescope, which was operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by an award issued by JPL/Caltech. Based on observations obtained at the international Gemini Observatory, a program of NSFs NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini Observatory partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovaçőes e Comunicaçőes (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). This work was enabled by observations made from the Gemini-North telescope, located within the Maunakea Science Reserve and adjacent to the summit of Maunakea. We are grateful for the privilege of observing the Universe from a place that is unique in both its astronomical quality and its cultural significance.

Facilities: Gemini (GMOS - , NIRI) - , HST (WFPC2 - , ACS - , WFC3) - , Mayall (NEWFIRM) - Kitt Peak National Observatory's 4 meter Mayall Telescope, and Spitzer (IRAC) - .

Software: astropy (Astropy Collaboration et al. 2022), dolphot (Dolphin 2016), DoPhot (Schechter et al. 1993), drizzlepac (Hack et al. 2021), emcee (Foreman-Mackey et al. 2013), forwardmodel (Kilpatrick & Rubin 2023), hst123 (Kilpatrick 2021), IRAF (Tody 1986), MOPEX (Makovoz & Khan 2005), progenitors (Kilpatrick 2023), pyraf (Science Software Branch at STScI 2012), swarp (Bertin 2010), and YSE-PZ (Coulter et al. 2022, 2023).

All data and analysis products presented in this article are available upon request. Analysis code and photometry used in this paper are available on GitHub: https://github.com/charliekilpatrick/progenitors; a copy of this material has been deposited to Zenodo: 10.5281/zenodo.8060641. The HST data used in this paper can be found in MAST: 10.17909/dqc4-yx93. The Spitzer Space Telescope data used in this paper can be found in the Spitzer Heritage Archive (IRSA 2022) (AORS: 44605696, 44605952, 45237504, 48186624, 50627072, 50627328, 50627840, 52776448, 52776704, 52776960, 52777216, 52777472, 52777728, 52777984, 52778240, 60830720, 60830976, 60831232, 60831488, 60831744, 60832000, 66022400, 66022656, 66022912, 66023168, 66023424, and 66023680). Gemini data (GN-2010A-Q-27) are publicly available on the Gemini data archive at https://archive.gemini.edu/ and NEWFIRM data (Prop IDs: 2012A-0540) are from the NOIRLab Astro Data Archive at https://astroarchive.noirlab.edu/.