SN 2023ixf was discovered in M101 within a day of the explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet UV spectra obtained with the Hubble Space Telescope 14 19 24 and 66 days after the explosion. Interaction between the supernova ejecta and circumstellar material CSM is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric Mg II emission on day 66. We compare our observations to CMFGEN supernova models that include CSM interaction inlineformula mmlmath overflowscrollmmlmover accenttruemmlmrowmmlmiMmmlmimmlmrowmmlmrowmmlmommlmommlmrowmmlmovermmlmoltmmlmommlmsupmmlmrowmmlmn10mmlmnmmlmrowmmlmrowmmlmommlmommlmn3mmlmnmmlmrowmmlmsupmmlmath inlineformula M SUBSUB yrSUP1SUP and find that the power from CSM interaction is decreasing with time from L SUBshSUB  5  10SUP42SUP erg sSUP1SUP to L SUBshSUB  1  10SUP40SUP erg sSUP1SUP between days 14 and 66. We examine the contribution of individual atomic species to the spectra on days 14 and 19 showing that the majority of the features are dominated by iron nickel magnesium and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae which show no definitive signs of CSM interaction and those with persistent signatures assuming the same progenitor radius and metallicity. Finally we show that the evolution and asymmetric shape of the Mg II   2796 2802 emission are not unique to SN 2023ixf. These observations add to the early measurements of dense confined CSM interaction tracing the massloss history of SN 2023ixf to 33 yr prior to the explosion and the density profile to a radius of 5.7  10SUP15SUP cm. They show the relatively short evolution from a quiescent red supergiant wind to high mass loss.
['Type II supernovae', 'Red supergiant stars', 'Stellar mass loss', 'Ultraviolet transient sources', 'Ultraviolet spectroscopy', 'Hubble Space Telescope', '1731', '1375', '1613', '1854', '2284', '761', 'Astrophysics - High Energy Astrophysical Phenomena', 'Astrophysics - Solar and Stellar Astrophysics']
Circumstellar Interaction in the Ultraviolet Spectra of SN 2023ixf 1466 Days After Explosion
{'Circumstellar Interaction in the Ultraviolet Spectra of SN 2023ixf 14-66 Days After Explosion': "K. Azalee Bostroem, 1, ∗ David J. Sand, 1 Luc Dessart, 2 Nathan Smith, 1 Saurabh W. Jha, 3 Stefano Valenti, 4 Jennifer E. Andrews, 5 Yize Dong ( 董 - 泽 ) , 4 Alexei V. Filippenko, 6 Sebastian Gomez, 7 Daichi Hiramatsu, 8, 9 Emily T. Hoang, 4 Griffin Hosseinzadeh, 10 D. Andrew Howell, 11, 12 Jacob E. Jencson, 13 Michael Lundquist, 14 Curtis McCully, 11, 12 Darshana Mehta, 4 Nicolas E. Meza Retamal, 4 Jeniveve Pearson, 1 Aravind P. Ravi, 4 Manisha Shrestha, 1 and Samuel Wyatt 15  \n1 Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA  \n2 Institut d'Astrophysique de Paris, CNRS-Sorbonne Universit'e, 98 bis boulevard Arago, 75014 Paris, France  \n3 Department of Physics and Astronomy, Rutgers, the State University of New Jersey,  \n136 Frelinghuysen Road, Piscataway, NJ 08854-8019, USA  \n4 Department of Physics and Astronomy, University of California, Davis, 1 Shields Avenue, Davis, CA 95616-5270, USA 5 Gemini Observatory, 670 North A'ohoku Place, Hilo, HI 96720-2700, USA  \n6 Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA  \n7 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218-2410, USA  \n8 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138-1516, USA  \n9 The NSF AI Institute for Artificial Intelligence and Fundamental Interactions, USA  \n10 Department of Astronomy & Astrophysics, University of California, San Diego, 9500 Gilman Drive, MC 0424, La Jolla, CA 92093-0424, USA  \n11 Las Cumbres Observatory, 6740 Cortona Drive, Suite 102, Goleta, CA 93117-5575, USA  \n12 Department of Physics, University of California, Santa Barbara, CA 93106-9530, USA  \n13 IPAC, Mail Code 100-22, Caltech, 1200 E. California Blvd., Pasadena, CA 91125  \n14 W. M. Keck Observatory, 65-1120 M¯amalahoa Highway, Kamuela, HI 96743-8431, USA  \n15 Astrophysics Science Division, NASA Goddard Space Flight Center, Mail Code 661, Greenbelt, MD 20771, USA", 'ABSTRACT': 'SN 2023ixf was discovered in M101 within a day of explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet (UV) spectra obtained with the Hubble Space Telescope 14, 19, 24, and 66 days after explosion. Interaction between the supernova ejecta and circumstellar material (CSM) is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric Mg II emission on day 66. We compare our observations to CMFGEN supernova models that include CSM interaction ( ˙ M < 10 -3 M ⊙ yr -1 ) and find that the power from CSM interaction is decreasing with time, from L sh ≈ 5 × 10 42 erg s -1 to L sh ≈ 1 × 10 40 erg s -1 between days 14 and 66. We examine the contribution of individual atomic species to the spectra on days 14 and 19, showing that the majority of the features are dominated by iron, nickel, magnesium, and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae which show no definitive signs of CSM interaction and those with persistent signatures assuming the same progenitor radius and metallicity. Finally, we show that the evolution and asymmetric shape of the Mg II λλ 2796, 2802 emission are not unique to SN 2023ixf. These observations add to the early measurements of dense, confined CSM interaction, tracing the mass-loss history of SN 2023ixf to ∼ 33 yr prior to the explosion and the density profile to a radius of ∼ 5 . 7 × 10 15 cm. They show the relatively short evolution from a quiescent red supergiant wind to high mass loss.  \nKeywords: Type II supernovae; Red supergiant stars; Stellar mass loss; Ultraviolet transient source; Ultraviolet spectroscopy; Hubble Space Telescope  \nbostroem@arizona.edu', '1. INTRODUCTION': "Hydrogen-rich supernovae (Type II) are thought to come from red supergiants (RSGs) with zero-age main sequence masses between ∼ 8 and ∼ 25 M ⊙ . While it is accepted that all RSGs lose mass, there are no quantitative mass-loss rates derived from stellar evolution theory. Adding to this, there is no consensus about the empirical mass-loss rates for RSGs; Massey et al. (2023) infer mass-loss rates higher than the canonical rates of de Jager et al. (1988) from the RSG luminosity function, while other studies (e.g. Beasor et al. 2020), measure lower mass-loss rates from the infrared (IR) excess of RSGs. Different mass-loss rates affect the final star that explodes, impacting the final hydrogen envelope mass as well as the density structure of the circumstellar material (CSM) into which the supernova explodes.  \nSignatures of mass loss can also be observed in supernovae, allowing us to constrain late-stage RSG mass loss using post-explosion observations. One definitive sign of CSM is the presence of narrow emission lines. As the number of young Type II supernovae discovered has increased with the prevalence of high-cadence wide-field surveys, it has become clear that a significant fraction of them have confined CSM (Khazov et al. 2016; Bruch et al. 2023; Jacobson-Gal'an et al. 2024). Early-time spectra of these supernovae display flash features: narrow emission lines with Lorentzian wings formed by the ionization of surrounding CSM by shock breakout and the interaction of the ejecta with CSM. These narrow features disappear a few days to weeks after explosion and the evolution then proceeds as a typical Type IIP/L supernovae.  \nWhile these lines disappear when the CSM is swept up, the influence of this CSM on the remaining photometric and spectroscopic evolution of the supernova is uncertain. Observations and recent theoretical modeling of the effects of the overall dynamics and CSM properties on the light curve and spectral evolution suggest that, even in the absence of narrow emission lines, the CSMejecta interaction will convert some of the kinetic energy of the supernova into thermal energy, boosting the supernova luminosity (Smith et al. 2015; Dessart et al. 2017; Dessart & Hillier 2022). Furthermore, this interaction will produce a bluer spectral energy distribution (SED) after a few days. The discrepancy between the U -V color curves of a supernova with and without CSM grows with time as the noninteracting supernova ejecta expand and cool (Dessart & Hillier 2022). Spectroscop-  \nically, high-velocity absorption features, from the cool dense shell that forms at the interface between ejecta and CSM, are predicted to be visible in H α , H β (Chugai et al. 2007), Na I D λλ 5890, 5896, and the Ca II nearinfrared (NIR) triplet λλλ 8498, 8542, 8662, and produce strong emission lines in the ultraviolet (Dessart & Hillier 2022). In particular, as the ejecta cool, Mg II λλ 2796, 2802 emerges in emission with increasing strength as the level of CSM interaction increases. Thus, the characterization of the long-term evolution of supernovae weeks to years after explosion probes RSG mass loss centuries to millennia prior to explosion.  \nObservational evidence supporting these predictions is unclear. Sample analyses found contradictory results about whether supernovae with flash features are systematically bluer and more luminous than noninteracitng supernovae (Khazov et al. 2016; Bruch et al. 2023). Additional contradictions are found in individual objects (e.g. SN 2018zd, Hiramatsu et al. 2021; SN 2020pni, Terreran et al. 2022; various, Jacobson-Gal'an et al. 2024; PTF11iqb, Smith et al. 2015; SN 2013fs, Bullivant et al. 2018). Spectroscopically, while an absorption feature is often observed blueward of H α (Guti'errez et al. 2017) and has been interpreted as high-velocity hydrogen due to CSM interaction, the relationship between the sample of supernovae with this feature and those with flash features is unclear. UV observations of Type II supernovae are especially sparse, with early-time data limited by the availability of sensitive rapid-response UV telescopes, and with late-time observations limited by the fading of the UV when CSM interaction is not present (Panagia et al. 1980; Pettini et al. 1982; Cassatella 1987; Baron et al. 2000; Fransson et al. 2005; Vasylyev et al. 2022, 2023a).  \nThe study of Type II supernovae with dense, confined CSM was recently enriched by the explosion of SN 2023ixf in the nearby galaxy M101. The proximity, early discovery and classification (Itagaki 2023; Perley & Gal-Yam 2023), and identification of narrow emission lines have led to this being the best-observed supernova with flash features to date. Of particular note is the periodic pre-explosion light curve of the RSG progenitor in the optical and IR (Jencson et al. 2023; Kilpatrick et al. 2023; Niu et al. 2023; Qin et al. 2023; Soraisam et al. 2023; Xiang et al. 2024), the incredibly deep multiband pre-discovery data which characterize the shock breakout and light-curve rise starting at an absolute magnitude of V ≈ -11 together with high-cadence photometry through the end of the plateau phase (JacobsonGal'an et al. 2023; Hiramatsu et al. 2023; Hosseinzadeh et al. 2023; Sgro et al. 2023; Teja et al. 2023; Yamanaka  \net al. 2023; Zimmerman et al. 2024; Bersten et al. 2024; Li et al. 2024; Martinez et al. 2024), intranight optical spectroscopy for more than a week after discovery (Bostroem et al. 2023a; Jacobson-Gal'an et al. 2023; Teja et al. 2023; Hiramatsu et al. 2023; Yamanaka et al. 2023; Zhang et al. 2023), early-time optical spectropolarimetry (Vasylyev et al. 2023b; Singh et al. 2024), highresolution spectroscopy of the first 17 days of evolution (Smith et al. 2023; Hiramatsu et al. 2023; Zimmerman et al. 2024), X-ray and radio detections (Berger et al. 2023; Grefenstette et al. 2023; Matthews et al. 2023; Chandra et al. 2024), and the first early-time UV spectra of a supernova with flash features (Zimmerman et al. 2024). While a coherent picture of these observations is still being explored, they point toward a confined, dense CSM with a radius R ≲ 10 15 cm and a mass-loss rate ˙ M ≈ 10 -2 -10 -4 M ⊙ yr -1 .  \nIn this paper, we describe UV spectroscopy of SN 2023ixf obtained with the Hubble Space Telescope (HST) at 14, 19, 24, and 66 days past explosion, all after the flash features had faded. Throughout the paper we use UTC 2023-05-18 18:00:00 (MJD 60082.75 ± 0.10; Hosseinzadeh et al. 2023) as the explosion epoch, E ( B -V ) = 0 . 0077 ± 0 . 0002 mag (Schlafly & Finkbeiner 2011) as the Milky Way (MW) extinction, and E ( B -V ) = 0 . 031 mag (Smith et al. 2023) as the host-galaxy extinction. In Section 2 we describe the data reduction and in Section 3 the evolution of the spectra, focusing on the UV. SN 2023ixf is compared to models with CSM interaction in Section 4 and to other supernovae with UV observations in Section 5. We close with a discussion of CSM interaction through day 66 in Section 6 and summarize our results in Section 7.", '2. OBSERVATIONS AND DATA REDUCTION': "HST observations were obtained on days 14, 19, 24, and 66 with the CCD detector of the Space Telescope Imaging Spectrograph (STIS) using the 52 '' × 0. '' 2 slit at the E1 position to mitigate flux loss to charge transfer efficiency loss and the G230LB, G430L, and G750L gratings (DD GO-17313; Bostroem et al. 2023b). Unfortunately, our planned visit on day 50 failed to acquire guide stars. However, the observation was successfully repeated on day 66. We obtained at least four exposures with each grating to automatically remove cosmic rays. A list of observations is given in Table 1.  \nBias, flat-fielding, and cosmic-ray rejection were automatically performed prior to download from the Mikulski Archive for Space Telescopes. Additionally, onedimensional (1D) spectra were automatically extracted, wavelength calibrated, and flux calibrated, prior to download for all G230LB and G430L observations ex-  \ncept the G230LB observation on day 66. At this phase, the UV flux was very low at the blue end, preventing the automatic identification of the spectrum location by the pipeline. These observations were manually extracted with the stistools.calstis pipeline, using the red end of the spectrum to identify the location of the trace.  \nThe G750L grating suffers fringing at the reddest wavelengths, which can be corrected using a contemporaneous fringe flat. We obtained a fringe flat with each visit using the 0 . '' 3 × 0 . '' 09 aperture, as recommended for the E1 aperture position. The fringe flat was applied to the G750L observations using the module stistools.defringe and 1D spectra were extracted using the stistools.x1d routine.  \nWe utilized data-quality flags 16 (high dark rate) and 512 (bad reference pixel) to eliminate bad pixels. In particular, we note that narrow emission features that appear in all spectra are flagged as high dark rate pixels and therefore removed from the analysis. We confirm that these are isolated pixels with high count rates in the dark files by examining the extraction location in the two-dimensional dark images themselves. Nevertheless, it is possible that there is a real unresolved emission line at the same pixel as at least one of these features. However, this would be contaminated by the high dark rate and therefore unable to provide any constraint on the CSM. Thus, given the characteristics described above, we treat these features as artifacts and do not consider them further.  \nWe add to this dataset the HST observations taken between days 3 and 11 as part of GO-17205 (Zimmerman et al. 2022), which have been presented by Zimmerman et al. (2024). 1D-extracted and flux-calibrated files were downloaded from the Mikulski Archive for Space Telescopes. As noted by Zimmerman et al. (2024), some of these exposures suffered from pointing, acquisition, and/or saturation issues, which led to unreliable flux due to the supernova not being aligned in the slit and the nonlinearity of the CCD response near the saturation point. To mitigate these issues, we scale all observations with each grating in a given visit to the spectrum with the most flux (assuming this is the best-centered observation) using a first-degree polynomial fit to the ratio of the spectra as the scale factor. A constant multiple of a few percent was then applied to align the G430L and G750L observations to the G230LB observation, and then the fluxes were combined by taking the median value at each wavelength, excluding pixels with data-quality flags of 16 and 512. We note that the majority of pixels are saturated in the range ∼ 3200-5000 ˚ A  \nTable 1. HST /STIS observations of SN 2023ixf from GO-17313  \nin visit 1, ∼ 2100-7600 ˚ A in visits 2 and 3, and ∼ 29007900 ˚ A in visit 4, making these wavelength ranges the most unreliable in terms of flux calibration. However, in this paper, we are primarily analyzing the later epochs (14 -66 days) which are not saturated and only use the early epochs ( < 10 days) to identify the presence and timing of relevant spectral features. Thus, we do not attempt to further correct for this effect.", '3. SPECTRAL EVOLUTION': "The near-UV (NUV) through NIR evolution of SN 2023ixf from day 3 to 66 is shown in Figure 1. No Type II supernovae (excluding Type IIn supernovae) since SN 1987A have been observed this early in the UV and with this cadence, revealing unprecedented details of the evolution. As has been noted by other authors, narrow emission lines are present in the early-time spectra, including notably a narrow P Cygni line from N IV λ 1718 which is present in the first spectrum and fades by day 8, and C III λ 2297 which is present from day 4 to day 8 (Yamanaka et al. 2023; Jacobson-Gal'an et al. 2023; Bostroem et al. 2023a; Teja et al. 2023; Hiramatsu et al. 2023; Zimmerman et al. 2024). These early narrow lines indicate dense CSM surrounding the supernova that has not been shocked. By day 14, these lines disappear and the optical spectrum is that of a young Type II supernova. Although slow to develop as a result of early interaction, a typical Type II supernova optical spectrum emerges with prominent Doppler-broadened P Cygni profiles in H I , He I , and other metals.  \nA more detailed view of the evolution from day 14 through day 66 is shown in Figure 2. Although still dominated by iron line blanketing (the broad absorption of flux due to a forest of iron lines that blend together), the UV continuum flux is clearly present in the day 14 spectrum. This is unusual for classical Type II supernovae in which the UV flux fades rapidly (Baron et al. 2000; Brown et al. 2007; Vasylyev et al. 2023a; Bostroem et al. 2023c).  \nWhile the flux blueward of 2500 ˚ A fades significantly by day 20, prominent absorption features can still be seen in the range 2500-3500 ˚ A, including Mg II , and just blueward Fe II /Fe III . In contrast with the fading UV flux, a prominent 'emission' line arises around 1925 ˚ A which reaches peak strength at 25 days before decreasing in strength.  \nBy day 66, there is very little UV flux as the supernova ejecta cool causing the SED to shift to longer wavelengths and iron absorption to increase in strength. However, the Mg II λλ 2796, 2802 doublet is clearly present as a boxy broad emission line with a slanted top. Also, at this time H α has fully developed a broad P Cygni profile, and other metal lines typical of Type II supernovae are present in the optical, including the Ca II NIR triplet λλλ 8498, 8542, 8662.  \nAcomparison of the Mg II and H α profiles provides insight into the origin of the Mg II emission (see Figure 3). When CSM is present, a cool dense shell forms at the interface between the ejecta and CSM. While the H α emission originates from the ejecta and the cool dense shell, the Mg II is thought to originate only from the  \nFigure 1. Evolution of the NUV to NIR spectra of SN 2023ixf from 3.5 to 66 days after explosion. Prominent features (N IV : 1718, 7109, 7122 ˚ A; C III : 2297 ˚ A; Mg II : 2796, 2802 ˚ A; He II : 3203, 4685.5, 4860 ˚ A; C IV : 5801, 5812 ˚ A; H α : 6563 ˚ A; He I : 6678.1 ˚ A; Ca II : 8498, 8542, 8662 ˚ A) are marked with vertical lines at their rest wavelength and labeled at the top of the figure. Narrow ISM absorption lines have been removed to highlight the supernova spectra.  \n<!-- image -->  \nshell and should thus be at a similar velocity as the outer edges of the ejecta. We compare Mg II and H α in Figure 3, finding that the profiles extend to similar velocities on both the red and blue sides, confirming that the Mg II is emitted from the cool dense shell at the outer edge of the ejecta. However, H α sits in the optical where there is continuum and forms throughout the ejecta and cool dense shell, leading to a P Cygni profile, while Mg II , without UV continuum and forming in the outer ejecta, is seen solely in emission. There is a shallow, narrow dip at the blue edge of the H α profile that could be absorption from the cool dense shell at ∼ 9200 km s -1 , although it is not clear if this feature is real given the signal-to-noise ratio (S/N) of the spectrum.  \nTo characterize the Mg II feature further, we fit three different profiles in Figure 3, representing different physical configurations. A spherical distribution of material with a Gaussian emitting profile will produce a Gaussian profile. Although the profile is clearly not Gaussian, it could be coming from a spherical distribution of material that is being attenuated. To explore the nonattenuated spherical distribution, we fit a Gaussian to the blue side of the Mg II profile from -10 , 000 to -7000 km s -1 . This fit is shown as a pink line in Figure 3; it has a mean of µ = 0 km s -1 and a FWHM of 7780 km s -1 . The observed blue side of Mg II is well described by a Gaussian profile. An alternate emitting region could  \nbe a sphere with a Gaussian profile that has had the center removed to form a thick shell. This configuration will have a flat-topped profile with Gaussian wings (Jerkstrand 2017) if the cool dense shell and ejecta are optically thin. Given our slanted top profile, we use the maximum of our observed flux as the flat-topped maximum of the model profile. This defines the inner radius velocity of the shell as v in = -6800 km s -1 , and the profile is shown as a dotted-blue line in the right panel of Figure 3.  \nTo model the slated top of the profile, we model the line as if it were generated by an asymmetric shell of material (i.e., a shell in which there is more flux from one side than the other reaching the observer). Following Kwok et al. (2023, 2024), we model this asymmetric shell as a Gaussian profile with an off-center hole cut out, which results in a slanted top instead of a flat-topped profile. We parameterize this model by the FWHM and mean of the Gaussian, the velocity of the center of the hole ( v c ), and the velocity of the radius of the hole ( v in ). We find FWHM = 7580 km s -1 , µ = -2995 km s -1 , v in = 5000 km s -1 , and v c = 1750 km s -1 , where v c and v in are measured relative to µ . This model fits the data exceptionally well and is shown as a green dashed line in the right panel of Figure 3. While the CSM is likely asymmetric, we caution that this parameterization is more likely to indicate different optical depths along the  \nFe II  \nFigure 2. Acomparison of the spectra of SN 2023ixf from 14.25 to 66.25 days post-explosion (solid lines) with the best-matching CSM interaction model (semi-transparent). The zero flux levels for each epoch are drawn as a dot-dashed line in the same color as the observed spectrum at that epoch. Over time, the UV emission declines and P Cygni features develop in the optical, until day 66, when only Mg II λλ 2796, 2802 is visible in the UV. This corresponds to a sequence of models with decreasing interaction power. Even though the model spectra do not match the SED of the observations, especially in the first two spectra where they appear to be cooler (resulting in lower-ionization species and more UV line blanketing), the match is impressive given the 1D nature and the ad-hoc implementation of the mass, density, and clumping of the cool dense shell. Prominent spectral features have been marked with dotted lines at 7000 km s -1 , the approximate location of the absorption throughout the spectroscopic sequence. The 'emission' feature at 1925 ˚ A is marked at rest with a dashed line. Narrow ISM absorption lines have been removed to better highlight the supernova spectra. Observed and model spectra at each epoch have been placed at the same distance with the same offset applied.  \n<!-- image -->  \nline of sight through the ejecta and cool dense shell, and should not be taken as a literal description of the ejecta. It does, however, demonstrate that this profile can be created by a shell of material in which a higher flux of photons from the blue side of the ejecta is reaching us, than photons from the red side. Additionally, it defines the inner edge of the blue side of the shell to be at velocity -6245 km s -1 (where the profile peaks), close to our earlier estimate of -6800 km s -1 . From these two values, we adopt v in = -6500 km s -1 . From these fits we conclude that the emitting region of the Mg II is well described by a Gaussian profile with a FWHM ≈ 7700 km s -1 and an inner radius of v in = -6500 km s -1 that has been attenuated on the red side due to the opacity of the ejecta and cool dense shell. We note that as the profile approaches the peak flux from the blue side, the rise is shallower than the Gaussian profile, possibly indicating asymmetry. Finally, we note that if the CSM is asymmetric, inner and outer shell velocities would represent projected velocities.  \nNarrow absorption lines are present in the spectra from Fe II λ 2344, λ 2374, λ 2383, λ 2587, λ 2600 and Mg II λλ 2796 , 2802, λ 2852 . 96. These lines are centered at the host redshift and there is no evolution in their equivalent widths with time. From this, we conclude that these lines are either from the interstellar medium (ISM) in the host galaxy or distant CSM and remove them from the remainder of the analysis.", '4. COMPARISON WITH LITERATURE MODELS': "Although by day 14 the narrow emission lines typically associated with strong CSM interaction have faded, it is possible that some energy from the interaction of the supernova ejecta with CSM is still contributing to the observed supernova spectrum. When supernova ejecta encounter CSM, a cool dense shell is formed between the forward and reverse shocks. The cool dense shell absorbs and reemits some of the shock power from CSM interaction. Dessart & Hillier (2022) published 1D (spherically symmetric) models of the interaction of supernova ejecta and cool dense shell with a low-density wind ( ˙ M < 10 -3  \nFigure 3. Left: H α and Mg II line profiles on day 66 in velocity space of SN 2023ixf compared to the CSM interaction models of Dessart & Hillier (2022) with L sh = 1 × 10 40 erg s -1 , scaled to the distance of SN 2023ixf. The shapes of both lines are well captured by the models, given that the cool dense shell is at a higher velocity in the models than in SN 2023ixf. Right: Mg II profile (black) fit with a pure Gaussian (pink) modeling spherical ejecta, a symmetric shell (blue dotted) modeling emission from the cool dense shell in the absence of optical-depth effects, and an asymmetric shell (green dashed) modeling the cool dense shell with optical depth.  \n<!-- image -->  \nM ⊙ yr -1 ). Unlike those of Dessart et al. (2017), which model interaction in a spherically symmetric dense, confined CSM, these models ignore the hydrodynamics of the CSM interaction in exchange for more accurate radiative transfer. While Dessart & Hillier (2022) compared their models to optical observations, they show that the signatures of this interaction are strongest in the UV, and our spectral series allows us to test these model predictions.  \nBriefly, these models start with a 15 M ⊙ progenitor of solar metallicity evolved with Modules for Experiments in Stellar Astrophysics (Paxton et al. 2011, 2013) which has been exploded with the radiation-hydrodynamics code V1D (Livne 1993; Dessart et al. 2010) and evolved to 10 days post explosion. On day 10, these observations are mapped to the nonlocal thermodynamic equilibrium radiation transport code CMFGEN (Hillier & Miller 1998; Hillier & Dessart 2012; Dessart et al. 2013; Hillier & Dessart 2019) assuming homologous expansion after this point. A cool dense shell is simulated by placing 0.1 M ⊙ of material at 11,700 km s -1 into the density profile. This velocity was chosen as the shock and cool dense shell are not expected to slow down much in the presence of a low-density CSM over a time scale of weeks (Dessart et al. 2017). The power from the CSM-ejecta  \ninteraction is then deposited at a constant rate into this region, mimicking the conversion of kinetic energy from the ejecta into radiative energy. Implicit in this implementation is the assumption that the cool dense shell is fully formed during the first few days after the explosion and does not grow in mass.  \nWith this setup, they produce models with no power deposited and continuous depositions of 1 × 10 40 , 1 × 10 41 , 5 × 10 41 , 1 × 10 42 , 2 . 5 × 10 42 , 5 × 10 42 , and 1 × 10 43 erg s -1 from day 15 through at least day 120 (some models are run longer). We emphasize that these models were not created for this specific supernova (or any supernova with early narrow emission lines). One difference between the model setup and SN 2023ixf is the presence of confined dense CSM, which produces the flash features observed in the early spectra. Dessart et al. (2017) find that in the presence of dense CSM, the maximum ejecta velocity decreased from 11,500 km s -1 in the low-density scenario to 7200 km s -1 in the highdensity scenario (see also Dessart & Jacobson-Gal'an 2023). We thus expect the velocity of the outer ejecta of SN 2023ixf to be slower than the models of Dessart & Hillier (2022). Another consequence of the dense CSM is that the initial conditions of the interaction region of the model may not be representative of the physi-  \nitions of the outer ejecta of SN 2023ixf, leading to a larger mismatch between the observations and the models in the earliest epoch than at later times.  \nWe visually compare each model at a given epoch to our observed flux, starting from the day 14.25 spectrum, and identify the model that best matches the observed flux. Although each of the models represents a steadystate wind, the cool, dense shell is optically thin to continuum radiation and radiates efficiently. For this reason, it will react promptly to a change in shock power. Thus, we expect that a decreasing mass-loss rate with radius would effectively have the same result as taking a different steady-state model for different epochs. The spectra of SN 2023ixf and best-matched models are shown in Figure 2 and the model properties are listed in Table 2. As expected, the model velocity is faster than the observed ejecta velocity leading to blueshifted and broader features. We find that the observations of SN 2023ixf indicate a progressive drop in interaction power over time, which translates into a drop in progenitor mass loss with increasing distance.  \nConsidering that these models were not produced for this specific supernova, it is impressive that no scaling beyond correcting for distance is required to match the observations and many of the features are reproduced, although at varying strengths. The first two epochs of models are too faint in the NUV and too bright in the optical, a discrepancy that improves with time. Nevertheless, we consider the best-matched model to be indicative of the CSM interaction level as the selected model most closely matches the observed flux: the model with more CSM interaction is too bright at all wavelengths and the model with less CSM interaction is even more discrepant for the majority of the spectrum. This consideration is further validated by the custom modeling in Section 6. One possible reason for the difference between the models and observations is that the temperature is too low in the emitting region in the early models. This is corroborated by the lack of He I λ 5875 in the models, whose formation requires higher temperatures (Dessart & Hillier 2005)  \nThe UV is dominated by a forest of iron-groupelement absorption lines that blend together, making it challenging to associate features with the emission or P Cygni profiles of individual species (Bostroem et al. 2023c; Dessart & Hillier 2005). To facilitate the identification of regions with contributions from different elements, we post-process the CMFGEN models to calculate the spectra on day 14 and day 19, omitting the bound-bound transitions of individual species. The full  \nspectrum is then divided by these spectra to obtain the relative flux contribution of each species. We perform this calculation for H, He, C, N, O, Ne, Na, Mg, Al, Si, S, K, Ca, Sc, Ti, Cr, Fe, Co, and Ni. Comparing the omitted spectrum to the full spectrum for each species, regions were identified as at least partially resulting from an individual species. These are marked in the top panel of each column of Figure 4, which shows the full model spectrum compared to the observed spectrum at the same epoch. Elements that have absorption features with a depth of > 2% of the total flux are shown in the bottom panel of Figure 4. In these panels, the strong influence of the iron-group elements (specifically Fe and Ni on day 14 and Fe, Ni, and Cr on day 19) on the overall spectrum is clear. In particular, we point out that the apparent developing 'emission line' at 1924 ˚ A is in fact a window of lower absorption between two prominent iron absorption regions, similar to the 'emission line' at 2970 ˚ A identified in the NUV spectra of SN 2022acko on day ∼ 20 (Bostroem et al. 2023c). Additionally, we show the continuum emission if no absorption is present in the top panel of each column of Figure 4, demonstrating how strongly the spectrum deviates from blackbody emission.  \nAnother emission line of particular note is Mg II in the day 66 spectrum, which is compared to the CMFGEN models in Figure 5. To account for the higher velocity of the models, we scale the models by the ratio of the blue edge of the feature in the observation to the blue edge in the model which preserves the morphology of the line. The Mg II feature is only present in models with CSM interaction after day ∼ 50 and rapidly increases in strength for greater interaction power. Our observation on day 66 falls between the models with a shock power of L sh = 1 × 10 40 erg s -1 and L sh = 1 × 10 41 erg s -1 and looks clearly distinct from models with no interaction.  \nIn the left panel of Figure 3, we compare the Mg II and H α profiles of our observation and the L sh = 1 × 10 40 erg s -1 model (without any shift in wavelength) on day 66. We find that both profiles are similar, although the model extends to bluer wavelengths due to the lack of dense CSM in the early supernova evolution, which caused the cool dense shell to form at lower velocities in SN 2023ixf. While the models have a large narrow absorption feature on the blue side of H α from the blue side of the cool dense shell, no such prominent feature is present in the observed H α profile. There is a hint of an absorption feature around 9200 km s -1 , although this is near the noise threshold of the spectrum and thus we cannot confidently associate it with this feature. On the red edge of H α , the models show a subtle red shelf.  \n)  \n1  \nÅ  \n1  \ns  \n2  \nm  \nc  \ng  \nr  \ne  \n(  \ny  \nt  \ni  \ns  \nn  \ne  \nD  \nx  \nu  \nl  \nF  \n<!-- image -->  \nFigure 4. Top: UV spectra of SN 2023ixf on day 14.25 (left) and 19.25 (right) compared to radiative-transfer models. Absorption lines from individual species are marked below the spectra. In the day 19.25 spectrum, the asterisk on Fe at 1924 ˚ A indicates that this is a window of minimal iron absorption rather than an emission line due to any atomic transition. The individual features are reproduced quite well, albeit at higher velocities in the models. Bottom: decomposition of the model spectrum at each epoch into species that contribute to the spectrum. Although the individual spectra are offset, the scales of the spectra are consistent, showing the relative intensity of each species. The dominance of iron-group elements (and in particular iron) in the spectra at both epochs is clear.  \n<!-- image -->  \nTable 2. Characteristics of the CSM around the progenitor of SN 2023ixf assuming v shock = 10000 km s -1 and v wind = 55 km s -1 .  \nNote -The last two rows bracket the true shock luminosity as the Mg II λλ 2796, 2802 emission lies between L sh = 10 40 and L sh = 10 41 erg s -1 .  \nFigure 5. Mg II line profiles of SN 2023ixf on day 66 in velocity space (black) compared to models with no power from CSM interaction (cyan) and the two lowest levels of CSM interaction modeled ( L sh = 1 × 10 40 erg s -1 ; blue, L sh = 1 × 10 41 erg s -1 ; magenta). This feature is present only in those models with CSM interaction, indicating that SN 2023ixf is interacting with CSM on day 66. The velocity scale of the models has been scaled to match the observations.  \n<!-- image -->  \nWith the S/N of our data, we cannot verify whether this feature is present.  \nThe model Mg II is very boxy, rising sharply at the blue edge of the H α profile to a plateau before falling off just after 0 km s -1 . This shape arises from a shell of material with the back side obscured by the cool dense shell and ejecta, leading to what appears to be blueshifted emission. Compared to our observations, the Mg II feature in the models is boxier and more blueshifted. Additionally, the Mg II emission in the observations rises more gradually to a peak and begins to fall off immediately, rather than plateau. Nevertheless, the observations have the same basic structure as the models, indicating a similar emitting mechanism and structure.", '5. COMPARISON WITH OTHER SUPERNOVAE': "We focus our comparison of SN 2023ixf to other supernovae on the small number of Type II supernovae with UV observations, defining a sample of objects with a range of CSM properties and supernova types (see Table 3). Although many of these objects continued to exhibit signs of CSM interaction later in their evolution (e.g., Pooley et al. 2002; Mauerhan & Smith 2012; Smith et al. 2017), we limit the discussion to their plateauphase evolution, which corresponds to our observations of SN 2023ixf.  \nSN 1993J was a Type IIb supernova, with the hydrogen lines fading and helium lines increasing in strength  \nsoon after explosion (Filippenko et al. 1993). Many lines of evidence indicated interaction with CSM including narrow emission lines (Benetti et al. 1994) and X-ray (Zimmermann et al. 1994) and radio (Pooley & Green 1993) detections. SN 1979C and SN 1998S are both Type II supernovae with steeply declining light curves and above-average photospheric phase luminosities. They also showed narrow emission lines for weeks after explosion; however, these did eventually disappear (Panagia et al. 1980; Leonard et al. 2000). SN 1979C was also detected at X-ray and radio wavelengths indicating interaction with CSM beyond the nominal RSG wind (Panagia et al. 1980). SN 1980K had a decline rate similar to that of SN 1979C, although it was about a magnitude fainter in the V band. While it did not exhibit narrow emission lines, it was detected in X-rays (Canizares et al. 1982) and radio (Weiler et al. 1992). SN 1987A was a peculiar Type II supernova, originating from a compact, blue supergiant progenitor. No CSM interaction was initially detected; however, months after explosion narrow emission lines appeared as high-energy photons ionized a ring of CSM. SN 1999em, SN 2021yja, and SN 2022wsp are considered fairly normal Type II supernovae (Leonard et al. 2000; Baron et al. 2000; GalYam et al. 2008; Hosseinzadeh et al. 2022; Vasylyev et al. 2022, 2023a), although Hosseinzadeh et al. (2022) note that SN 2021yja was exceptionally blue at early times and Vasylyev et al. (2023b) find that the absorption of H α and H β in SN 2022wsp was suppressed in the P Cygni profiles which they ascribe to CSM interaction. Finally, SN 2005cs and SN 2022acko were low-luminosity supernovae (Brown et al. 2007; Bostroem et al. 2023c). The only possible sign of CSM interaction in these supernovae is a 'ledge' feature in the early spectra, which has been attributed to CSM, although other possibilities exist (see Pearson et al. 2023, for a detailed discussion). SN 2021yja also showed the 'ledge' feature. Details on the explosion epochs, extinction values, and spectra considered in this analysis are given in Table 3. We note that UV data also exist for SN 2010jl, a superluminous Type IIn supernova. Throughout its evolution, the UV spectrum is dominated by the cool dense shell with narrow spectral features and enhanced flux, in contrast to the other supernovae in this sample where the spectra are dominated by the supernova photosphere. We therefore exclude it from our analysis.  \nIn Figure 6 and Figure 7 we show a comparison of the spectra of SN 2023ixf to this sample of UV spectra, scaled to the distance of SN 2023ixf, on day ∼ 10 and day ∼ 20, respectively. For this comparison, we exclude SN 1979C and SN 1998S, as their UV flux is significantly higher than that of SN 2023ixf and all mod-  \nconsidered in this paper. On day 10, SN 2023ixf is similar in brightness to SN 2021yja. Although it has a lower flux, SN 2022wsp is spectroscopically similar. The Mg II λλ 2796, 2802 doublet is present in absorption in SN 1980K, SN 2021yja, SN 2022wsp, and SN 2022acko, although it is blueshifted to higher velocities in SN 2021yja and SN 2022acko. Additionally, blueward of Mg II , the Fe absorption feature identified by Vasylyev et al. (2023a) in SN 2022wsp is clearly present in SN 2023ixf and SN 1980K. Although not a pronounced absorption trough, it is also possibly present in SN 2021yja and SN 2022acko and blended with Mg II in SN 1999em. Like SN 2023ixf, the UV spectra of SN 1980K, SN 2021yja, and SN 2022wsp all show an absorption complex of Fe and Ni between 1750 ˚ A and 2050 ˚ A. The Fe absorption complex between 2050 ˚ A and 2550 ˚ A, which is notably shallower in SN 2023ixf and SN 1980K than models predict, is fairly well reproduced in SN 2021yja, although all of the observed profiles have less prominent individual troughs than the model.  \nOne of the most distinctive features of CSM interaction is excess UV flux. On day 10, we see that SN 1980K has a higher flux than SN 2023ixf, SN 2021yja has a very similar flux while SN 2022wsp is a bit lower, SN 2022acko is significantly lower, and the remaining supernovae have virtually no flux on day 10. We note, however, that both excess UV flux and the slope of the UV SED are highly dependent on the assumed distance, extinction, and explosion epoch of the supernovae. In particular, the redshift-independent distance of NGC 4269 ranges from 4 to 7.8 Mpc (see Van Dyk et al. 2019, for a thorough discussion), significantly affecting the luminosity of SN 1980K. While the distances and thus overall luminosities are uncertain, individual features can indicate different levels of CSM interaction. For example, the development of an Fe absorption complex redward of 3000 ˚ A is not present in SN 1980K, SN 2021yja, and SN 2023ixf. This feature only appears in models with lower levels of CSM interaction. This line does appear to be developing in SN 2022wsp, consistent with the interpretation that the lower flux level is indicative of less CSM interaction. Comparing with the models, this suggests that SN 2021yja and SN 2022wsp have close to 1 × 10 43 erg s -1 interaction power added to the ejecta luminosity on day 10.  \nBostroem et al. (2023c) noted that SN 1999em has an extremely low UV flux relative to comparable sequences. It is interesting to consider both the anomalous SED and features of SN 1999em relative to the other supernovae in this sample, as SN 1999em is a relatively prototypical Type IIP supernova (although a high-velocity feature was detected in both H α and H β during the photo-  \nheric phase, interpreted as a sign of CSM interaction; Chugai et al. 2007). The explosion epoch used is halfway between the last nondetection and first detection. It is possible that, at discovery, SN 1999em was up to 5 days older, as suggested by Dessart & Hillier (2006), which would make the UV flux less unusual; however, even at a later phase, the shape of the UV flux is different from that of the other supernovae - either more suppressed below 2750 ˚ A or with increasing flux redward of that. The presence of an'emission' feature owing to the window of absorption at 1924 ˚ A further supports a later explosion epoch. Given this uncertainty, SN 1999em is plotted in both Figure 6 and Figure 7.  \nOn day 20, the UV flux has decreased in all supernovae. As was observed in SN 2022wsp (Vasylyev et al. 2023a), the Fe feature blueward of Mg II has faded in SN 2023ixf. Interestingly, SN 2021yja shows an extended feature that encompasses both the Fe and Mg II features, possibly indicating that these are blended. The slope and features of SN 2021yja and SN 2022wsp are very similar to the model with L sh = 2 . 5 × 10 42 erg s -1 , although they all have a larger flux. The slope of SN 2023ixf and SN 1999em is not reproduced by any model, even if an overall flux offset is applied. It is possible that models customized for each supernova with more sophisticated physics would provide a better fit. For example, a compelling time-dependent model for SN 1999em is presented in Dessart & Hillier (2006). The flux of SN 2021yja has fallen to a level similar to that of SN 2022wsp, indicating a lower level of power from CSM interaction. The fact that the flux of SN 2021yja fell relative to SN 2023ixf and SN 2022wsp from the flux levels on day 10 gives us confidence that this is not a distance effect. The window of absorption at 1924 ˚ A creates a clear 'emission line' in SN 2021yja, SN 2022wsp, and SN 2023ixf.  \nIn interacting supernovae, the Mg II emission is attributed to emission from the cool dense shell, which has been excited by radiation from the reverse shock (Chevalier & Fransson 1994; Dessart & Hillier 2022). Originating from a shell of material, the emission-line profile is expected to be broad, boxy, and blue shifted (Dessart & Hillier 2022). It has been noted in SN 1979C 14 yr after explosion (Fesen et al. 1999), in SN 1993J on day 670 (Fransson et al. 2005), in SN 1995N on day 943 (Fransson et al. 2002), in SN 1998S on days 28-485 (Fransson et al. 2005), and in SN 2010jl on days 34-543 (Fransson et al. 2014). However, these supernovae represent the more strongly interacting objects in our sample. It remains to be seen if this feature is present in supernovae that do not have any period of strong interaction with dense CSM.  \nFigure 6. Top : comparison of the UV spectrum of SN 2023ixf on day 9 with that of other Type II supernovae at a similar epoch (phase for each supernova given in the panel legend; all spectra scaled to the distance of SN 2023ixf): SN 1980K, SN 1987A, SN 1999em, SN 2005cs, SN 2021yja, SN 2022wsp, and SN 2022acko. At NUV wavelengths, SN 2023ixf is among the brightest in the sample. Bottom : day 9 (black) and 14 (gray) spectra of SN 2023ixf compared with CMFGEN models having varying degrees of shock power on day 14. Although the majority of UV observations have phases around day 10, day 14 is the first epoch modeled by Dessart & Hillier (2022). We show both SN 2023ixf spectra to demonstrate possible evolution from the spectra in the top panel and the models in the bottom. The day 14 spectrum best matches the model with L sh = 5 × 10 42 erg s -1 below ∼ 3000 ˚ A but is brighter than the model at shorter wavelengths. Elements are marked in both panels in black at -7000 km s -1 to match the SN 2023ixf absorption.  \n<!-- image -->  \nWe searched all available NUV spectra of our comparison sample of Type II supernovae for the broad, boxy Mg II emission that was observed on day 66 in SN 2023ixf. As Mg II is present in absorption on day 25 and then in emission on day 66, we do not know the exact timing of this transition. Additionally, this timing is related to the ejecta and CSM properties and thus may not be a universal property of Type II supernovae. We therefore consider all spectra with phases greater than 30 days. One complicating factor is the presence of Mg II absorption from both the MW and the host galaxy, which contaminates the supernova emission and can make it challenging to assess the shape of the emission profile. Whenever possible, we remove this narrow absorption from our data.  \nWe do not detect this feature in SN 1987A, although this is perhaps not surprising given its blue supergiant progenitor and lack of UV flux. The feature is clearly detected in the remaining supernovae as a broad, boxy, blue-shifted profile. In SN 1980K, we marginally detect it. However, the S/N of the spectra is fairly low and the feature is not as strong as in other supernovae. In SN 1979C we detect it developing already on day 24, although it could be as old as 32 days if the supernova explosion occurred immediately after the last nondetection. From this date, Mg II evolves from a blue-shifted, flat-topped profile into one with a slanted profile (higher blue edge) starting around day 60. Similarly, Mg II evolves in SN 1998S from a boxy, blue-shifted, symmetric profile on day 72 to an asymmetric profile on days 237 and 485. In SN 1993J, the Mg II emission is fairly constant with a boxy, asymmetric profile from day 173 to 669. This is the earliest epoch of data for SN 1993J after the UV continuum has faded. While the feature is clearly present as broad emission in SN 1995N and SN 2010jl (days 30-572), the ISM absorption from both the Milky Way and the host galaxy obscure the shape of the profile and we do not include these in our analysis.  \nThe Mg II emission features of SN 1979C, SN 1980K, SN 1993J, and SN 1998S are compared to that of SN 2023ixf in Figure 8. We choose the closest available epoch to day 66, and for SN 1998S we also include the first epoch at which the profile is asymmetric. From these comparisons, we find that if a change in the shape of the Mg II emission is observed, the emission begins as a symmetric, boxy, blue-shifted profile and always evolves to the same asymmetric, boxy profile, with more flux on the blue side. For supernovae that, like SN 2023ixf, always show an asymmetric profile, the asymmetry is always higher on the blue side. The consistency of this asymmetry tells us that this is not an effect of asymmetric CSM that we would expect  \nto observe from a different viewing angle, but rather an optical-depth effect. This is corroborated by the Dessart &Hillier (2022) models, which show a blue-shifted, boxy emission profile with no emission from ejecta at velocities below 7000 km s -1 (i.e. the emission is isolated to the cool dense shell). Additional validation is derived from a custom model that reproduces the models of Dessart & Hillier (2022) but with the cool dense shell at 8500 km s -1 and shows the same slanted profile as we observed in the Mg II feature. While optical-depth effects are responsible for the blue-shifted and slated profile at early times, at later times dust formation likely attenuates the red side of the feature. This sample shows a variety in maximum velocity of the blue side of these features, indicating diversity in the location of the cool dense shell, as well as other supernova parameters such as the ejecta mass, hydrogen envelope mass, explosion energy, and CSM mass. However, they uniformly rise to maximum more slowly than the models, indicating a different cool dense shell density profile than is implemented in the models.  \nOverall, this comparison confirms that UV spectra are sensitive to different levels of CSM interaction, and the Mg II feature can reveal weak CSM interaction even when it is not clearly seen in the optical. The sparse time series of the other supernovae that we have assembled indicates a diversity in features and their strength, velocity, and evolution, but a larger sample is required to evaluate further. Additionally, while the models reproduce many features, they are not able to reproduce the slope and possibly the overall flux of the observations, indicating room for further detailed modeling.", '6.1. Testing Model Predictions': "With our sequence of observations, we can test the following predictions of signatures of CSM interaction made by Dessart & Hillier (2022) through day 66.  \n- 1. Higher luminosity and flux in all filters for supernovae having more interaction, with this effect being more prominent in the UV and at later times when the ejecta luminosity has faded.\n- 2. Emergence of the Mg II doublet as broad, boxy emission ∼ 50 days after explosion, increasing in luminosity until day ∼ 175.\n- 3. Weakened absorption features as the photospheric spectrum of the supernova is filled in by emission from the cool dense shell.\n- 4. A sharp-edged red shelf extending beyond the photospheric emission in H α originating from the cool  \nFigure 7. Top : comparison of the spectrum of SN 2023ixf on day 19 with that of other Type II supernovae at a comparable epoch (phase for each spectrum given in the panel legend; all spectra scaled to the distance of SN 2023ixf): SN 1980K, SN 1987A, SN 1993J, SN 2021yja, SN 2022wsp, and SN 2022acko. SN 2023ixf remains UV bright, especially when compared with SN 2021yja, which has faded to a flux similar to that of SN 2022wsp. The explosion epoch of SN 1999em is uncertain and thus we show the spectrum here, too, to demonstrate its similarity to the other spectra at these epochs. Bottom : UV spectrum of SN 2023ixf on 19 days (black) compared with the CMFGEN models having varying degrees of shock power on day 20 (colors). Although at optical wavelengths the spectrum best matches the model with L sh = 2 . 5 × 10 42 , below ∼ 2750 ˚ A the spectrum best matches the model with L sh = 5 × 10 42 . Elements are marked in both panels in black at -7000 km s -1 to match the SN 2023ixf absorption. The asterisk denotes the feature caused by a window between two strong iron absorption troughs.  \n<!-- image -->  \nFigure 8. Comparison of the Mg II emission line in SN 2023ixf (day 66; black) to that of other supernovae that show this feature (SN 1979C, brown; SN 1993J, yellow; and SN 1998S, magenta). All lines have a blueshifted peak. In all but the first SN 1998S spectrum, the profiles are asymmetric indicating that the red side of the line is attenuated. The gap in some spectra is due to the removal of narrow Galactic and host ISM absorption.  \n<!-- image -->  \ndense shell, although this feature is only starting to become visible on day 50  \n- 5. A strong, narrow absorption feature at the velocity of the dense shell in H α , H β , Na I D λλ 5890, 5896, and the Ca II NIR triplet λλλ 8498, 8542, 8662 increasing in strength with time starting around day 30.  \nFor the first point, we find that the UV luminosity is consistent with models with CSM interaction at all epochs, and not just limited to the early epochs when the narrow emission features are present. The fact that the shock power in the best-matched model decreases with the lookback time indicates that the ejecta are sweeping through CSM created from decreasing massloss rates. This is consistent with the light-curve modeling results of Singh et al. (2024) which require extended CSM ( ˙ M=10 -4 M ⊙ yr -1 to R = 10 16 cm) to match the late-time UV excess. We detect a boxy Mg II emission feature at 66 days, indicating CSM interaction even without the presence of a UV continuum, confirming the second expectation. As described in the third point, the absorption components in the P Cygni profiles of the observed spectra are very shallow and slow to develop, especially in H α which is only clearly identified in the final epoch. This is also notes by Singh et al. (2024). The fourth prediction of a sharp red shelf in H α is subtle on day 66 and not visible prior to this date. At the S/N of our data, we cannot confidently identify it nor  \nrule out its presence. Similarly, the narrow high-velocity absorption component, if present, is significantly weaker than predicted in the models. This could be due to a lower density CSM, asymmetric CSM where the bulk of the CSM is offset from the line of sight, and/or the fragmenting of the cool dense shell due to Rayleigh Taylor instabilities which would spread this feature out in the velocity space.  \nTo explore whether clumping and lowering the velocity of the cool dense shell would provide a better fit to our observations, we calculate a custom CMFGEN model. We place the cool dense shell at 8500 km s -1 . At this velocity, the ejecta are optically thick. This has two effects: first, the material has not been able to cool, so the temperature is higher; and second, it responds much more quickly to the power from the CSM interaction. However, we found that this change alone was not enough to reproduce our observations and continued to further tune the model. To better approximate the physical conditions that create the emission from the cool dense shell (i.e., the reprocessing of X-rays from the forward and reverse shock), we add a uniform field of X-rays starting at 8000 km s -1 with a total power of 5 × 10 42 erg s -1 instead of inserting power into the cool dense shell directly. Given the energy of these X-rays, we add higher energy ions to our model, namely, C III , N III , O III , Mg III , and Ne III . We also increase the clumping from 1% to 10% volume filling factor outside of 8000 km s -1 , keeping the total mass the same.  \nIn Figure 9 we compare this model on day 14.6 to our day 14.2 observed spectrum and the best-matching model from Dessart & Hillier (2022) on day 14.6. Reducing the velocity of the cool dense shell aligns the absorption and emission. However, the resulting increase in the temperature of the cool dense shell shifts flux from the optical to the UV, causing the UV to overshoot the observed spectrum and the optical to underestimate the true flux. The change in power deposition from direct injection into the cool dense shell to an Xray field reduces the UV flux to align with the observations but does not resolve the low optical flux issue. The increase in clumping from a volume filling factor of 1%10% increases the optical flux to match the full SED. The observed spectrum is very well matched by our custom model, demonstrating the success of this approach. However, further improvements could be made to better match the strength of features in the UV, most notably the Fe II absorption complex in the range 2250-2760 ˚ A, which is too strong in the models.  \n6.2. CSM Density Profile and Progenitor Mass-loss History  \nFigure 9. Observed spectrum on day 14 compared to the L sh = 5 × 10 42 erg s -1 model from Dessart & Hillier (2022, blue) and our new model with the cool dense shell at lower velocity, increased clumping, and a more complete atomic model (pink). The new model is able to match the velocity of most of the features and provides a significantly better match to the SED.  \n<!-- image -->  \nUsing the best models at each epoch, we derive the density profile of the RSG progenitor of SN 2023ixf, which is given in Table 2 and plotted in Figure 10. We use the shock luminosity of each model, assuming a shock velocity of v sh = 10 , 000 km s -1 to find the CSM density at each epoch from ρ ( r ) = L sh / (2 πv 3 sh r 2 ) g cm -3 . Note that this is a simplified model that assumes both spherical symmetry and a shock velocity. From this, we find that the density decreases more rapidly than a constant wind, going from ρ = 5 . 2 × 10 -16 g cm -3 at 14 days after explosion to ρ = 4 . 8 × 10 -20 g cm -3 on day 66.  \nIf we further use the RSG wind velocity of v wind = 55 km s -1 measured by Zhang et al. (2023), we can derive a mass-loss rate from the density and a time of mass loss from the radius (see Figure 11). We find ˙ M = 8 . 7 × 10 -4 M ⊙ yr -1 14 days after the explosion, decreasing to 4 . 4 × 10 -4 M ⊙ yr -1 19 and 24 days after the explosion, and finally arriving at a canonical RSG wind of 1 . 7 × 10 -6 M ⊙ yr -1 at 66 days after the explosion. This traces the mass loss from elevated rates ∼ 7 yr before the explosion down to nominal RSG rates ∼ 35 yr before explosion. The clear presence of Mg II emission on day 66 indicates that we can detect ejectaCSM interaction, even with canonical RSG winds of 10 -6 M ⊙ yr -1 . This mass-loss rate and timing are dependent on the assumed wind velocity. If Zhang et al. (2023) measured a wind that had already been radiatively accelerated, the velocity could be lower, which would increase the period of time prior to the explosion  \nFigure 10. CSM density of SN 2023ixf as a function of the radius assuming a shock velocity of 10,000 km s -1 . The CSM density constraints from this paper are plotted with star symbols. The filled symbols are constrained by the overall UV and optical flux. The open stars represent the value from the Mg II line, which lies between L sh = 10 40 erg s -1 (black) and L sh = 10 41 erg s -1 (pink). Density constraints from flash spectroscopy are shown at the epoch of the first spectrum, day 4, when the majority of the narrow features disappear, and day 7, when all of the narrow features disappear (Bostroem et al. 2023a; Jacobson-Gal'an et al. 2023; Zhang et al. 2023). Mustard triangles represent a mass-loss rate of 10 -2 M ⊙ yr -1 and green inverted triangles represent a mass-loss rate of 10 -3 M ⊙ yr -1 . Values derived from the H α luminosity from Zhang et al. (2023) are shown as light-blue solid squares and the value from Bostroem et al. (2023a) is a light-blue open square. CSM density from NuSTAR X-ray observations (Grefenstette et al. 2023) are plotted as filled blue circles and X-ray observations from Chandra (Chandra et al. 2024) are shown as open blue circles. Lines of constant density corresponding to mass-loss rates of 10 -1 , 10 -2 , 10 -3 , 10 -4 , 10 -5 , and 10 -6 M ⊙ yr -1 are shown as solid gray lines. While there is significant scatter in the density, it is clear that density is decreasing more rapidly than would be expected from constant mass-loss rates down to the density expected for canonical RSG winds.  \n<!-- image -->  \nprobed by these observations and decrease the mass-loss rate.  \nThe CSM density and mass-loss rates derived in this paper can be combined with other published values to paint a more complete picture of the CSM around the progenitor of SN 2023ixf. We plot these densities and mass-loss rates in Figure 10 and Figure 11. Grefenstette et al. (2023) use NuSTAR X-ray obser-  \nvations on day 4.4 and day 11 to measure a column density of N H = 26 +5 -7 × 10 22 atoms cm -2 and N H = 5 . 6 ± 2 . 7 × 10 22 atoms cm -2 , respectively. Additionally, from Chandra observations on days 13 and 86, Chandra et al. (2024) find a column density of N H = 2 . 5 +0 . 40 -0 . 34 × 10 22 atoms cm -2 and N H = 0 . 36 +0 . 22 -0 . 17 × 10 22 atoms cm -2 , respectively. From the column density, the CSM density can be found assuming a density profile from a constant wind and that N H = ∫ ∞ R sh ρ ( r ) /m p dr . With these assumptions, the density at each epoch is ρ = ( N H m P ) /R sh . From the best-matching models of the observed flash features between 1 and 7 days, Bostroem et al. (2023a), JacobsonGal'an et al. (2023), and Zhang et al. 2023 infer a massloss rate of 10 -3 -10 -2 M ⊙ yr -1 from the CMFGEN models of Dessart et al. (2017) and Jacobson-Gal'an et al. (2023). As the line emission in these models is dependent on the density and mass loss is a parameterization of this, we use the model wind velocity of v wind = 50 km s -1 to find the density rather than the measured wind velocity. We plot the density derived from flash features at the time of the first spectrum (day 1.1), day 4 when many of the narrow features disappear, and day 7 when the H α narrow component disappears. Bostroem et al. (2023a) and Zhang et al. (2023) use the H α luminosity to calculate the CSM radius and density. Bostroem et al. (2023a) find R CSM ≳ 8 . 7 × 10 13 cm and ρ CSM = 3 . 4 × 10 -14 g cm -3 . Zhang et al. (2023) determine R CSM = 2 . 33, 2.59, and 2 . 74 × 10 14 cm and ρ CSM = 7 . 41, 7.26, and 7 . 31 × 10 -15 g cm -3 .  \nThese combined measurements show that the massloss rate of the progenitor of SN 2023ixf was relatively low until quite close to the time of core collapse ( R ≲ 5 × 10 15 cm or t ≲ 33 yr prior to explosion), at which point it increased to the large densities inferred from flash-feature observations. However, this figure also highlights that the systematics, uncertainties, and simplifying assumptions in these different techniques lead to 1 order of magnitude uncertainties in mass-loss rates and densities.", '7. SUMMARY': "We present UV observations of SN 2023ixf 14-66 days after the explosion. SN 2023ixf is UV bright, with continuum emission through day 24. Its UV spectrum shows weaker metal absorption than other Type II supernovae that do not exhibit definitive signs of CSM interaction. In the optical, the H α P Cygni profile is slow to develop and has a shallow absorption component, which we infer is due to the additional emission from CSM interaction. On day 66, the Mg II profile and H α emission component have a similar extent in veloc-  \nFigure 11. Same as Figure 10 for the mass-loss history of SN 2023ixf assuming a shock velocity ( v sh = 10 , 000 km s -1 ) and the CSM wind velocity ( v wind = 55 km s -1 ). In its last ∼ 30 yr of life, the mass-loss rate increased from a nominal RSG wind of ∼ 10 -6 M ⊙ yr -1 to ∼ 10 -2 M ⊙ yr -1 .  \n<!-- image -->  \ny space, but extremely different profile shapes, resulting from the lack of UV continuum from the supernova ejecta. We examine the contribution of different species to the UV spectrum of SN 2023ixf on days 14 and 19, finding that the majority of the absorption in the UV spectrum is due to iron, nickel, and magnesium at both epochs, with additional contributions by chromium at the later epoch.  \nWe compare SN 2023ixf to CMFGEN models of supernova ejecta with power from CSM interaction (Dessart & Hillier 2022). Impressively, even though these models were not made for a supernova like SN 2023ixf, which interacted with dense CSM during its early evolution, they match the overall evolution and the presence of individual features. It is likely that the cool dense shell in SN 2023ixf is broader and centered at a lower velocity (8500 vs. 11,700 km s -1 ) and higher temperature than the models. This is confirmed by a custom model that places the cool dense shell at 8500 km s -1 . Additionally, the predicted presence of the Mg II emission after 50 days is confirmed, demonstrating that SN 2023ixf is interacting with CSM at all epochs through day 66.  \nWe visually identify the model that best matches our observed spectrum at each epoch, finding that the power from CSM interaction decreases with time. From this, we conclude that at larger radii, there is lower density CSM, suggesting that the mass-loss rate increased  \nrapidly just prior to the explosion. While the flash features favor a confined CSM, these disappear after a week, limiting the radius out to which the mass-loss history can be traced by this method. Our UV observations are sensitive to much lower CSM densities, and we use them to trace the CSM density to significantly larger radii and the mass-loss history to earlier epochs. With a shock velocity of 10,000 km s -1 , our observations probe the CSM density between a radius of 1 . 2 × 10 15 cm to 5 . 2 × 10 15 cm and show the density evolving from 5 . 2 × 10 -16 g cm -3 to 4 . 8 × 10 -20 g cm -3 over this range. Using the wind velocity of 55 km s -1 , we find that the mass-loss rate decreases from 8 . 7 × 10 -4 M ⊙ yr -1 in our first observation to 1 . 7 × 10 -6 M ⊙ yr -1 in our last observation, showing the dramatic change in mass-loss rates 33-7 yr before the explosion. Additionally, the mass-loss rate derived from the day 66 spectrum shows that 33 yr before the explosion, the RSG progenitor's mass-loss rate was consistent with a quiescent RSG wind (Beasor et al. 2020).  \nWe also compare SN 2023ixf with other supernovae having UV spectra around days 10 and 20. These objects cover a range of UV fluxes, similar to the range spanned by the CMFGEN models with varying levels of CSM interaction. This suggests that we can use the UV to sensitively distinguish varying levels of CSM interaction. Around day 10, we see that Mg II is clearly present in absorption in the spectra with the most CSM interaction. Often appearing next to it is Fe II ; this can be seen as a distinct feature in SN 2023ixf and SN 2022wsp, and it may be present in the other supernovae either blended with Mg II owing to its formation at higher velocities or as a weaker feature. Mg II and Fe II are weaker around day 20 in all supernovae. Additionally, the apparent emission line at 1924 ˚ A is caused by a window of lower iron absorption, rather than emission from an individual species, demonstrating the complexity of associating a single species with an individual feature in UV spectra.  \nBased on the shape and width of the Mg II emission in the day 66 spectrum, we conclude that this emission is from the cool dense shell, which is well described by a thick shell with a Gaussian density profile: FWHM ≈ 7700 km s -1 and an inner radius velocity of ∼ -6500 km s -1 . We model the emission with an asymmetric Gaussian shell, finding that the profile is consistent with more of the observed emission coming from the blue side than the red side. Additionally, in other supernovae, this feature is either symmetric or asymmetric with a higher blue side, leading us to conclude that the profile shape is due to attenuation of the far side of the shell, likely from opacity in the shell and ejecta at early times and dust  \nattenuation in the later-time spectra (e.g. SN 1998S at 240d).  \nFinally, we compile the density and mass-loss measurements of SN 2023ixf from the literature using a consistent shock velocity of 10,000 km s -1 and a wind velocity of 55 km s -1 . We show that the density profile from 9 . 5 × 10 13 cm to 7 . 5 × 10 15 cm (0.5-42 yr before explosion) decreases from 1 . 1 × 10 -12 g cm -3 to 4 . 9 × 10 -20 g cm -3 . With a wind velocity of 55 km s -1 , this corresponds to a change in mass-loss rate of ∼ 10 -2 M ⊙ yr -1 to ∼ 10 -6 M ⊙ yr -1 . This shows that over the final ∼ 40 yr of the progenitor's life, the mass-loss rate went from that of a quiescent RSG to an extremely high mass-loss rate. The densities and mass-loss rates of these different techniques span about an order of magnitude, highlighting the uncertainties and systematic errors in our methods, even when some consistent assumptions are made.  \nSN 2023ixf is the best-observed Type II supernova since SN 1987A, providing us with a Rosetta Stone to interpret other observations and against which to benchmark theoretical predictions. Although relatively few UV observations of Type II supernovae exist (and even fewer time series), it is clear that this wavelength range is rich in information about the temperature, composition, and dynamics of the CSM and the ejecta. It also probes the pre-supernova mass-loss history, tracing it back to quiescent RSG winds 40 years before explosion and can provide insights into the role of CSM interaction in the diversity of optical properties that we observe in Type II supernovae. In particular, time series which look for the presence of Mg II emission for all Type II supernovae and trace its evolution to even later epochs are warranted inform our understanding of mass loss in RSGs, especially in the years leading up to explosion. The observations of SN 2023ixf provide a critical link between early- and late-time UV observations, and also between 'normal' Type II supernovae and stronglyinteracting Type II supernovae.", 'T able 3 . Prop erties of Sup erno v ae with UV Sp ectra Analyzed in This P ap er': 'Note -+ SN 1980K do es not h a v e a reliable explosion ep o c h. W e adopt an explosion ep o c h 7 da ys prior to an observ ed sp ectrum on 1980 No v em b er 1 whic h sho ws broad H α emission. † See also Dessart & H illier 2006 , ∗ See also Dessart et al. 2008 . (A): Sc hlafly & Finkb einer 2011 . (B): Gall et al. 2015 . (C): F errarese et al. 1996 . (D): Barb on et al. 1982b . (E): P anagia et al. 1980 . (F): P a lum b o 1982 . (G): Ben v en uti et al. 1982 . (H): F esen et al. 1999 . (I): V an Dyk et al. 2019 . (J): Barb on et al. 1982a . (K): P ettini et al. 1982 . (L): Roma niello et al. 2000 . (M): Castagnoli et al. 1987 . (N): Scuderi et al. 1996 . (O): W amstek er et al. 1987 . (P): Cassatella et al. 1987 . (Q): Cassatella 1987 . (R): Kirshner et al. 1987 . (S): P anagia 1988 . (T): Pun et al. 1995 . (U): F reedman et al. 1994 . (V): Le wis et al. 1994 . (W): Ergon et al. 2014 . (X): Jeffery et a l. 1994 . (Y): F ransson et al. 2005 . (Z): Willic k et al. 1997 . (AA): B os tro em et al. 2023a . (BB ): Leonard et al. 2000 . (CC): Len tz e t al. 2001 . (DD): F ransson et al. 2005 . (EE): Leonard et al. 2003 . (FF): Elmhamdi et al. 2003 . (GG): Silv erman et al. 2017 . (HH): Baron et al. ( 2000 ). (I I): T ully 1988 . (JJ): Tsv etk o v et al. 2006 . (KK): Gal-Y am et al. 2008 . (LL): Sabbi et al. 2018 . (MM): P astorello et al. 2006 . (NN): Baron et al. 2007 . (OO): Bro wn et al. 2007 . (PP): Zhang et al. 2012 . (QQ): Stoll et al. 2011 . (RR): F ransson et al. 2014 . (SS): Hosseinzadeh et al. 2022 . (TT): V asyly ev et al. 2022 . (UU): V asyly ev et al. 2023a . (VV): Anand et al. 2021 . (WW): Bostro em et al. 2023c .  \n- We thank the referee for their thoughtful comments. 1\n- This research is based on observations made with the 2\n- NASA/ESA Hubble Space Telescope obtained from the 3\n- Space Telescope Science Institute, which is operated by 4\n- the Association of Universities for Research in Astron5\n- omy, Inc., under NASA contract NAS 5-26555. These 6  \n7  \nobservations are associated with programs GO-17313  \n- and GO-17205. Huge thanks to our HST program coor8\n- dinator Alison Vick and our contact scientist Dan Welty 9\n- for assistance in planning and executing the observations. The data described here may be obtained from the MAST archive at doi:10.17909/787k-c897. 10 11 12  \n13  \n14  \nK.A.B. is supported by an LSST-DA Catalyst Fellow-  \nship; this publication was thus made possible through  \n- the support of grant 62192 from the John Templeton 15\n- Foundation to LSST-DA. A.V.F. is grateful for financial 16\n- assistance from the Christopher R. Redlich Fund and 17\n- many other donors. Time domain research by D.J.S. 18  \n19  \nand team is supported by NSF grants AST-1908972  \n- and 2108032, and by the Heising-Simons Foundation 20\n- under grant No. 20201864. Research by Y.D., S.V., 21\n- N.M.R, E.H., and D.M. is supported by NSF grant AST22\n- 2008108. 23', 'REFERENCES': 'Facilities:  \nHST(STIS)  \nSoftware: astropy (Astropy Collaboration 2013; 2018; 2022), CMFGEN (Hillier & Miller 1998; Dessart et al. 2013; Hillier & Dessart 2019), MatPLOTLIB (Hunter 2007), NumPy (Harris et al. 2020), Scipy (Virtanen et al. 2020), stistools (Sohn 2019)  \nZhang, J., Lin, H., Wang, X., et al. 2023, Science Bulletin, 68, 2548, doi: 10.1016/j.scib.2023.09.015  \nZhang, T., Wang, X., Wu, C., et al. 2012, AJ, 144, 131, doi: 10.1088/0004-6256/144/5/131  \nZimmerman, E. A., Gal-Yam, A., Bruch, R. J., et al. 2022, Explosions in Real-Time: Rapid UV Supernova Flash Spectroscopy, HST Proposal. Cycle 30, ID. #17205 Zimmerman, E. A., Irani, I., Chen, P., et al. 2024, Nature, 627, 759, doi: 10.1038/s41586-024-07116-6 Zimmermann, H. U., Lewin, W., Predehl, P., et al. 1994, Nature, 367, 621, doi: 10.1038/367621a0'}