SN 2023ixf: Dense CSM Type II Supernova

• SN 2023ixf is a nearby Type II supernova characterized by a dense circumstellar environment and flash-ionization features that serve as benchmarks for studying progenitor mass loss and explosion physics.
• Multi-wavelength observations from optical to X-ray reveal a rapid luminosity rise, a bright plateau, and complex ejecta–CSM interactions that enhance our understanding of shock dynamics and mass-loss history.
• Diverse analyses—including flash spectroscopy, hydrodynamical modeling, and dust evolution studies—yield conflicting progenitor mass estimates, highlighting the challenges of reconciling high-mass dusty solutions with low-mass constraints.

SN 2023ixf is a Type II supernova in Messier 101 whose proximity, early discovery, and unusually dense circumstellar environment made it one of the most intensively studied core-collapse supernovae of the decade. The literature treats it as a nearby benchmark for connecting red-supergiant progenitor physics, pre-supernova mass loss, shock breakout in circumstellar material, late-time ejecta–CSM interaction, dust formation, and multi-messenger constraints. Published studies adopt distances of 6.4 Mpc, $6.71\pm0.14$ Mpc, $6.85\pm0.15$ Mpc, and 6.9 Mpc, and they consistently place first light about one day before the 2023 May 19 discovery in clear light at 14.9 mag (Jacobson-Galán, 10 Jul 2025, Hiramatsu et al., 2023).

1. Discovery, classification, and observational setting

SN 2023ixf was discovered by Koichi Itagaki on 2023 May 19 17:27:15 UT and classified hours later as a Type II supernova with flash-ionization features. Multiple early-light analyses place the explosion epoch near MJD $60082.74$–$60082.79$, and the event occurred in the southeastern spiral arm of M101, near an H II region (Jacobson-Galán, 10 Jul 2025, Hiramatsu et al., 2023). Its brightness and distance made it the closest and brightest core-collapse supernova in nearly a decade in one IceCube study, which quotes 6.4 Mpc and $B$-band magnitude 10.8 (Mand et al., 9 Jul 2025).

The early optical evolution was unusually luminous for a Type II event. One densely sampled photometric study measured $M_V=-18.18\pm0.09$ mag at 2023-05-25 21:37 UTC, while another reported a rapid rise of about 5 days to $M_V\approx-18.2$ mag followed by a bright plateau at $M_V\approx-17.6$ mag through 30 days (Sgro et al., 2023, Hiramatsu et al., 2023). The review literature treats this luminosity excess as a defining signature of strong circumstellar interaction superposed on an otherwise Type II photospheric evolution (Jacobson-Galán, 10 Jul 2025).

A recurrent point in the literature is that the earliest spectrum was not representative of a quiescent, “clean-envelope” Type II supernova. Instead, SN 2023ixf exhibited a flash-ionization phase with narrow emission and electron-scattering wings, after which it transitioned into a more ordinary Type II photospheric state. This distinction matters because progenitor and explosion inferences depend sensitively on whether the dense inner CSM is modeled explicitly (Jacobson-Galan et al., 2023, Jacobson-Galán, 10 Jul 2025).

2. Early-time circumstellar interaction and the compact inner CSM

The first days after explosion established SN 2023ixf as a canonical flash-spectroscopy event. UV/optical spectra showed H I, He I/II, C IV, and N III/IV/V with narrow cores and broad symmetric wings attributed to photo-ionization and electron scattering in dense, close-in CSM. These flash features persisted for about 5–8 days, then faded as the ejecta overran the confined material (Jacobson-Galan et al., 2023, Hiramatsu et al., 2023).

Several studies derived a compact CSM scale from those timescales. Optical light-curve color evolution and flash-feature persistence imply a characteristic radial extent of $\sim(3-7)\times10^{14}$ cm in one study, while CMFGEN and HERACLES modeling favored dense CSM confined to $r=(0.5-1)\times10^{15}$ cm with $6.85\pm0.15$0 for $6.85\pm0.15$1, corresponding to enhanced mass loss during the last $6.85\pm0.15$2-6 years before explosion (Hiramatsu et al., 2023, Jacobson-Galan et al., 2023). The review literature synthesizes these results as evidence for a very dense inner zone inside $6.85\pm0.15$3 cm, embedded within a broader wind-like environment (Jacobson-Galán, 10 Jul 2025).

X-ray observations independently required a dense but spatially limited medium. NuSTAR detected hard X-rays at $6.85\pm0.15$4 d and $6.85\pm0.15$5 d, with spectra described by hot thermal bremsstrahlung through a thick neutral absorber whose column dropped from $6.85\pm0.15$6 to $6.85\pm0.15$7. The absorbed $6.85\pm0.15$8–79 keV luminosity stayed near $6.85\pm0.15$9, and the interpretation was forward-shock emission from a dense but confined CSM with $60082.74$0 at $60082.74$1 cm and $60082.74$2 for $60082.74$3 (Grefenstette et al., 2023).

Early IR data added another layer to the same picture. NEOWISE-R detected a clear IR excess at day 3.631, fit by a $60082.74$4 K blackbody with radius $60082.74$5 cm, interpreted as circumstellar dust heated by UV emission from shock–CSM interaction. At day 10.836, the SED was instead consistent with ejecta-dominated emission, suggesting that the initial shock-powered reprocessing phase was short-lived (Dyk et al., 2024). A plausible implication is that the immediate environment was not a single homogeneous wind, but a structured dusty CSM with strong radial gradients.

3. Progenitor star and the mass debate

Pre-explosion imaging established a red supergiant progenitor candidate, but not a single uncontested mass. Archival HST, Spitzer, and ground-based near-IR data consistently indicate a dusty, variable RSG with strong mid-IR excess and long-timescale variability of roughly 1000 days, yet progenitor-mass estimates differ sharply because they depend on dust treatment, pulsation phase, assumed circumstellar extinction, and stellar-evolution models (Kilpatrick et al., 2023, Niu et al., 2023, Qin et al., 2023).

Approach	Key inference	Representative studies
HST red-source identification with modest extra dust	$60082.74$6–10 $60082.74$7 dusty RSG	(Pledger et al., 2023)
HST/Spitzer/NIR SED with dusty shell	$60082.74$8, $60082.74$9	(Kilpatrick et al., 2023)
C-rich dusty SED + local SFH	$60082.79$0–$60082.79$1 from SED; $60082.79$2–$60082.79$3 from environment	(Niu et al., 2023)
Time-dependent dusty-RSG SED	$60082.79$4, $60082.79$5 mag	(Qin et al., 2023)
CFHTLS $60082.79$6-band + HETDEX environment	$60082.79$7	(Liu et al., 2023)
Hydrodynamical and nebular constraints	$60082.79$8, $60082.79$9, or 10–15 $B$0	(Bersten et al., 2023, Folatelli et al., 14 Feb 2025, Michel et al., 17 Mar 2025)

The lower-mass pre-explosion interpretation identifies a single progenitor candidate consistent with Source A in HST imaging, fits a dusty-shell SED with $B$1, $B$2 K, dust-shell radius $B$3, and total circumstellar mass $B$4, and infers an initial mass near $B$5 or, more broadly, $B$6–$B$7 (Kilpatrick et al., 2023). A still lower-mass HST-based interpretation argued for $B$8–10 $B$9, but explicitly emphasized blending and dust as major uncertainties (Pledger et al., 2023).

Higher-mass interpretations treat circumstellar obscuration and variability more aggressively. One SED study favored a C-rich dusty envelope, $M_V=-18.18\pm0.09$0, $M_V=-18.18\pm0.09$1 K, $M_V=-18.18\pm0.09$2, and $M_V=-18.18\pm0.09$3–$M_V=-18.18\pm0.09$4, with a mass-loss rate $M_V=-18.18\pm0.09$5 at least $M_V=-18.18\pm0.09$6 years before explosion (Niu et al., 2023). Another, using a time-dependent dusty-RSG model, derived $M_V=-18.18\pm0.09$7 K, $M_V=-18.18\pm0.09$8, $M_V=-18.18\pm0.09$9, $M_V\approx-18.2$0 mag, and $M_V\approx-18.2$1, yielding $M_V\approx-18.2$2 and periodic variability with $M_V\approx-18.2$3 days (Qin et al., 2023). HETDEX plus CFHTLS $M_V\approx-18.2$4-band imaging pushed the preferred mass even higher, to $M_V\approx-18.2$5, while also finding low local extinction and near-solar metallicity on the metal-rich $M_V\approx-18.2$6 branch (Liu et al., 2023).

Subsequent hydrodynamical and nebular analyses moved the consensus back toward the low-mass end. One-dimensional radiation-hydrodynamics modeling of the bolometric light curve and photospheric velocities favored $M_V\approx-18.2$7, $M_V\approx-18.2$8 erg, and $M_V\approx-18.2$9, with an upper limit $M_V\approx-17.6$0 (Bersten et al., 2023). Late-time nebular work using [O I], [Ca II], and H$M_V\approx-17.6$1 found 10–15 $M_V\approx-17.6$2 or $M_V\approx-17.6$3, reinforcing the view that SN 2023ixf may have arisen from a relatively low-mass red supergiant despite the high-luminosity dusty-SED solutions (Folatelli et al., 14 Feb 2025, Michel et al., 17 Mar 2025). The resulting controversy is not merely observational scatter; it reflects genuinely different treatments of circumstellar obscuration, pulsation, envelope stripping, and the degree to which CSM interaction biases both pre- and post-explosion inferences.

4. Explosion properties, plateau behavior, and geometry

Hydrodynamical modeling of the photospheric phase treats SN 2023ixf as a broadly normal Type II explosion once the strongest early CSM-interaction phase is excluded. The preferred low-mass model gives $M_V\approx-17.6$4, pre-SN mass $M_V\approx-17.6$5, radius $M_V\approx-17.6$6, explosion energy $M_V\approx-17.6$7 erg, and nickel mass $M_V\approx-17.6$8, with Fe II $M_V\approx-17.6$9 used as the primary photospheric velocity tracer (Bersten et al., 2023). Independent nebular analyses infer $\sim(3-7)\times10^{14}$0Ni masses near $\sim(3-7)\times10^{14}$1–$\sim(3-7)\times10^{14}$2, explosion energies spanning $\sim(3-7)\times10^{14}$3 to $\sim(3-7)\times10^{14}$4 erg or, in plateau-scaling work, $\sim(3-7)\times10^{14}$5–$\sim(3-7)\times10^{14}$6 erg for alternative models (Michel et al., 17 Mar 2025, Forde et al., 16 Apr 2025).

The plateau itself was short and morphologically informative. One nebular study measured a plateau of $\sim(3-7)\times10^{14}$7 days, with midpoint of the drop at $\sim(3-7)\times10^{14}$8 days, and interpreted incomplete $\sim(3-7)\times10^{14}$9-ray trapping with $r=(0.5-1)\times10^{15}$0 days as evidence for low ejecta mass (Michel et al., 17 Mar 2025). Another light-curve study argued that plateau luminosity and duration can be matched by a low-mass, relatively unstripped progenitor, but that the steep drop at the end of the plateau is better reproduced by a higher-initial-mass progenitor with substantial prior mass loss and a reduced H-rich envelope (Forde et al., 16 Apr 2025). This suggests that the envelope structure, rather than the total stellar mass alone, is the key quantity controlling the plateau morphology.

Spectropolarimetry added a geometric constraint not captured by spherically symmetric light-curve fits. Six optical spectropolarimetric epochs from $r=(0.5-1)\times10^{15}$1 to $r=(0.5-1)\times10^{15}$2 days measured continuum polarization $r=(0.5-1)\times10^{15}$3 at days $r=(0.5-1)\times10^{15}$4 and $r=(0.5-1)\times10^{15}$5, followed by a drop to $r=(0.5-1)\times10^{15}$6 by day $r=(0.5-1)\times10^{15}$7, together with a $r=(0.5-1)\times10^{15}$8 rotation of the polarization position angle. The interpretation is that the ejecta emerged from a dense asymmetric CSM between days $r=(0.5-1)\times10^{15}$9 and $6.85\pm0.15$00, and that the outer ejecta and the circumstellar material do not share the same symmetry axis (Vasylyev et al., 2023). This result is frequently invoked to explain why optical, X-ray, radio, and progenitor inferences are not all simultaneously satisfied by a single spherical wind model.

5. Nebular evolution, distant CSM, and the onset of shock-powered late phases

By the nebular phase, SN 2023ixf had moved beyond the compact inner CSM probed by flash spectroscopy and began sampling larger-radius structures. A GMOS spectrum at 445 days showed a dramatic transformation relative to 259 days: H$6.85\pm0.15$01 developed a complex profile consisting of a central peaked component plus a broad boxy component, interpreted as a mixture of radioactive ejecta emission and ejecta–CSM interaction (Folatelli et al., 14 Feb 2025). In that picture, the boxy H$6.85\pm0.15$02 traces a shell at $6.85\pm0.15$03 cm, implying mass loss about 500–1000 years before explosion for a $6.85\pm0.15$04 wind, while [O I] and [Ca II] remained comparatively unaffected and continued to support a progenitor mass of 10–15 $6.85\pm0.15$05 (Folatelli et al., 14 Feb 2025).

A related nebular study tracked four spectra from day $6.85\pm0.15$06 to $6.85\pm0.15$07 and found narrow, blueshifted oxygen lines corresponding to an ejected oxygen mass $6.85\pm0.15$08, again favoring $6.85\pm0.15$09. The same work emphasized unusually broad, multi-component H$6.85\pm0.15$10, with high-velocity hydrogen features emerging around day $6.85\pm0.15$11 at $6.85\pm0.15$12, possibly from shock interaction with a dense hydrogen-rich region at $6.85\pm0.15$13 cm (Michel et al., 17 Mar 2025). A plausible implication is that SN 2023ixf did not encounter a single recent outflow, but a radially stratified CSM that records mass loss on timescales from years to centuries.

The most explicit late-time synthesis comes from the panchromatic study extending through the first two years. It constructed UV-to-IR spectral energy distributions at $6.85\pm0.15$14 and $6.85\pm0.15$15 days over $6.85\pm0.15$16–30 $6.85\pm0.15$17, found that the multiband light curve followed a standard radioactive decay decline rate after the plateau until $6.85\pm0.15$18 days, and concluded that shock-powered emission then began to dominate. The transition was temporally consistent with 0.3–10 keV X-ray detections and broad boxy line emission reprocessing shock luminosity in a cold dense shell between forward and reverse shocks (Jacobson-Galán et al., 15 Aug 2025). Using absorbed radioactive decay power plus the detected X-ray luminosity, that work inferred total shock-powered emission consistent with nearly complete thermalization of the reverse-shock luminosity in a continuous wind-like CSM with $6.85\pm0.15$19 for $6.85\pm0.15$20, and showed that CMFGEN models with $6.85\pm0.15$21–$6.85\pm0.15$22 and $6.85\pm0.15$23–$6.85\pm0.15$24 of silicate dust in the cold dense shell and/or inner ejecta can reproduce the global properties of the late-time $6.85\pm0.15$25-day UV-to-IR spectra (Jacobson-Galán et al., 15 Aug 2025).

6. Infrared echoes, internal dust, and multi-messenger legacy

Dust-related diagnostics in SN 2023ixf evolved from early circumstellar reprocessing to later internal formation. The earliest NEOWISE-R and optical/UV/NIR SED at day 3.631 required a hot $6.85\pm0.15$26 K component plus a cooler $6.85\pm0.15$27 K IR component at radius $6.85\pm0.15$28 cm, consistent with dust in the progenitor’s circumstellar shell heated by shock-interaction radiation (Dyk et al., 2024). Much later, optical/NIR/JWST work spanning 150–750 days separated a radiative-flash circumstellar infrared echo from in-situ dust formation: the 1.8–33.6 day IR excess was attributed to echo-dominated emission from pre-existing dust, while progressive H$6.85\pm0.15$29 red-wing attenuation from about 132 days, CO detected by about 217 days, and a growing late silicate-bearing component were taken as the onset of internal dust formation (Singh et al., 14 Mar 2026).

In that nebular-phase dust analysis, H$6.85\pm0.15$30 profile modeling from 141 to 418 days yielded a silicate-equivalent internal dust mass rising from $6.85\pm0.15$31 to $6.85\pm0.15$32, whereas late JWST SED fitting required a colder silicate-bearing component growing to $6.85\pm0.15$33 by 601 days and $6.85\pm0.15$34 by 723 days (Singh et al., 14 Mar 2026). This suggests that dust signatures in SN 2023ixf cannot be reduced to a single mechanism: early IR luminosity was dominated by circumstellar echo physics, while late asymmetries and far-IR continuum increasingly required newly formed CDS/ejecta dust.

Multi-messenger follow-up did not detect a non-electromagnetic counterpart, but the null results are themselves constraining. IceCube found no significant neutrino excess in either a 4-day or 32-day search window, placing 90% confidence upper limits of $6.85\pm0.15$35 and $6.85\pm0.15$36, corresponding to isotropic-equivalent energies of $6.85\pm0.15$37 erg and $6.85\pm0.15$38 erg for an $6.85\pm0.15$39 spectrum (Mand et al., 9 Jul 2025). Review work also notes non-detections by Fermi-LAT and LIGO-Virgo-KAGRA, alongside strong X-ray and radio detections, reinforcing the picture that SN 2023ixf was an electromagnetic benchmark for dense-CSM interaction rather than a detected high-energy neutrino or gravitational-wave source (Jacobson-Galán, 10 Jul 2025).

The scientific legacy of SN 2023ixf lies in the unusually complete linkage it provides between pre-explosion progenitor imaging, flash spectroscopy, UV/X-ray/radio diagnostics of dense compact CSM, polarimetric evidence for asymmetry, nebular evidence for larger-radius shells, and late-time shock power and dust formation. The main unresolved issue is not whether circumstellar interaction mattered—it did—but how to reconcile the high-mass dusty-progenitor solutions with the lower-mass hydrodynamical and nebular results, and how to map the apparently multi-epoch, asymmetric mass-loss history onto a single pre-supernova evolutionary channel. That tension is precisely why SN 2023ixf remains a reference event for Type II supernovae with dense, confined, and structured circumstellar material (Jacobson-Galán, 10 Jul 2025).