N132D: Oxygen-Rich Supernova Remnant

• N132D is a young, oxygen-rich supernova remnant characterized by high-velocity, hydrogen-poor ejecta, a complex 3D structure, and an estimated age of ~2500–3100 years.
• Multiwavelength studies have reconstructed a torus-like ejecta (with a distinct runaway knot) and revealed intricate shock dynamics from optical, X-ray, and infrared observations.
• Its gamma-ray spectrum and shock behavior strongly support hadronic cosmic-ray acceleration, while debates over the progenitor mass and explosion geometry remain unresolved.

N132D is a young, oxygen-rich core-collapse supernova remnant in the Large Magellanic Cloud, located near the LMC bar and close to the galaxy’s rotation center. It is the X-ray-brightest remnant in the LMC and a luminous GeV–TeV source, and it has been used as a laboratory for three-dimensional ejecta reconstruction, reverse-shock evolution, shock–cloud interaction, nonequilibrium plasma physics, dust formation, and cosmic-ray acceleration (Vogt et al., 2010, Collaboration et al., 2017, Vink et al., 2021). Its observed structure is markedly multi-phase: an optical and infrared O-rich ejecta system occupies the interior, Fe-rich X-ray plasma is concentrated toward the southern interior, the outer shell is shaped by strong density contrasts, and molecular and atomic gas surround and in places penetrate the remnant.

1. Classification, age, and evolutionary state

N132D belongs to the rare class of oxygen-rich young supernova remnants, whose optical filaments are dominated by high-velocity, hydrogen- and helium-poor ejecta rich in $\alpha$-process elements such as O, Ne, S, Ar, and Ca (Vogt et al., 2010). It is generally interpreted as the product of a core-collapse explosion, although subtype assignments differ across the literature: prior optical/UV work is summarized as favoring a Type Ib origin, the infrared literature describes it as a Type Ib/Ibc object, and comparison with Cassiopeia A has been used to suggest a possible Type IIb-like event (Collaboration et al., 2017, Rho et al., 2023, Law et al., 2020).

Published age estimates are not identical but occupy a relatively narrow range. A three-dimensional optical reconstruction derived an explosion age of $2450 \pm 195$ yr from the mean expansion of the O-rich ejecta, HST proper motions yielded $2770 \pm 500$ yr, deep RGS work characterized the object as about $2500$–$3100$ yr old, and a molecular-environment study quoted about $\sim 3150$ yr (Law et al., 2020, Banovetz et al., 2023, Suzuki et al., 2020, Sano et al., 2015). These values collectively place N132D between the canonical very young ejecta-dominated remnants and the more evolved mixed-morphology systems.

Its evolutionary classification reflects that intermediate status. Optical kinematics have been interpreted as showing that the remnant is approaching the end of the reverse-shock phase before entering the fully thermalized Sedov stage, while hard X-ray and gamma-ray studies describe it as being in transition from a young remnant with strong particle acceleration toward a middle-aged remnant dominated by thermal X-rays and hadronic gamma rays (Vogt et al., 2010, Bamba et al., 2018). A plausible implication is that N132D is best understood not as a single-phase remnant, but as a system in which young ejecta signatures and more mature swept-up-shell signatures coexist.

2. Three-dimensional ejecta structure and multiwavelength morphology

Integral-field and long-slit spectroscopy show that the optically bright oxygen-rich ejecta are organized into a large coherent structure rather than a simple filled shell. WiFeS observations reconstructed most of the [O III] ejecta into a ring of about $\sim 12$ pc in diameter inclined by $\sim 25^\circ$ to the line of sight, while later Magellan/IMACS+GISMO data described the dominant structure as a broken torus tilted by about $28^\circ$ with respect to the plane of the sky and having a radius of $4.4$ pc at $2450 \pm 195$0 kpc (Vogt et al., 2010, Law et al., 2020). In the latter reconstruction, the ejecta exhibit a blue-shifted radial-velocity asymmetry from about $2450 \pm 195$1 to $2450 \pm 195$2, and the average expansion velocity of $2450 \pm 195$3 was converted to age using the ballistic relation

$2450 \pm 195$4

with

$2450 \pm 195$5

The ballistic interpretation is supported independently by HST proper motions. Using 120 individual O-rich knots over a 16 yr baseline, the center of expansion was measured at $2450 \pm 195$6, $2450 \pm 195$7 (J2000), with a $2450 \pm 195$8 uncertainty of $2450 \pm 195$9, and the proper motions were found to remain consistent with homologous expansion (Banovetz et al., 2023). The same study reported a median proper motion of $2770 \pm 500$0, corresponding to about $2770 \pm 500$1 at $2770 \pm 500$2 kpc, reinforcing the picture of long-lived, only weakly decelerated ejecta moving through a low-density cavity.

A recurrent feature in optical reconstructions is the “runaway knot” (RK), spatially separated from the main torus and nearly perpendicular to its plane. The Magellan reconstruction assigned it a total space velocity of $2770 \pm 500$3 and noted its coincidence with X-ray emission enhanced in Si relative to both the LMC and N132D’s bulk ejecta, motivating interpretation as material originating deep within the progenitor (Law et al., 2020). HST proper-motion work revised this kinematic picture downward, combining a typical motion of about $2770 \pm 500$4 in the RK region with a Doppler speed of $2770 \pm 500$5 to obtain a three-dimensional velocity of about $2770 \pm 500$6 (Banovetz et al., 2023). The existence of a peculiar high-velocity feature is therefore robust, whereas its exact dynamical extremity depends on the reconstruction method.

Infrared observations extend the ejecta morphology beyond the optical. Spitzer, WISE, and Herschel data show high-velocity [Ne II] $2770 \pm 500$7m, [Ne III] $2770 \pm 500$8m, and [O IV] $2770 \pm 500$9m concentrated in a central ring that coincides with the optical and X-ray ejecta, and Herschel imaging reveals a cold dust component spatially associated with the same central ejecta (Rho et al., 2023). The fitted dust mass of that central ejecta region is $2500$0, and the $2500$1m dust feature requires forsterite grains (Rho et al., 2023). This establishes that the interior structure is simultaneously kinematic, chemical, and dusty.

3. Circumstellar and interstellar environment

N132D is expanding into a highly structured environment rather than a uniform medium. Mopra CO data identified three molecular clouds apparently interacting with the remnant, two of which lie to the west and east and form a cavity-like CO structure along the X-ray shell, while a third lies to the northeast (Sano et al., 2015). The total molecular mass of those three clouds was estimated as $2500$2 using the conversion factor $2500$3, and this mass was argued to be about an order of magnitude smaller than earlier estimates because only a subset of the larger giant molecular cloud is actually interacting with the remnant (Sano et al., 2015).

ALMA and ASTE later resolved the molecular environment at much finer scale. Diffuse CO emission was found not only at the southern edge but also inside the X-ray shell, where ALMA resolved nine molecular clouds, labeled A–I, at an angular resolution of $2500$4, corresponding to about $2500$5 pc at the distance of the LMC (Suzuki et al., 2020). Their typical sizes and masses are about $2500$6 pc and $2500$7, with individual cloud masses of about $2500$8–$2500$9 and densities of about $3100$0–$3100$1 (Suzuki et al., 2020). Inside the X-ray shell the ratio

$3100$2

reaches $3100$3–$3100$4, in contrast to $3100$5 for the southern giant molecular cloud outside the remnant, and this was interpreted as evidence that shock heating has occurred (Suzuki et al., 2020).

The geometry of the interaction is also constrained spectroscopically. For cloud F near the center of the remnant, Chandra spectroscopy found additional LMC absorption of $3100$6, a forward-shock component whose normalization went to zero, and a hotter thermal component with $3100$7 keV and $3100$8, consistent with recent shock heating (Suzuki et al., 2020). Combined with the mismatch between the CO-derived proton column and the X-ray absorbing column, this led to the interpretation that cloud F is shock-engulfed and located on the near side of the remnant (Suzuki et al., 2020).

Optical integral-field spectroscopy resolves the dense cloud population that the blast wave is crushing. From 20 shocked clouds, shock models yielded a mean cloud shock velocity of about $3100$9, a mean pre-shock density of about $\sim 3150$0, a mean shock ram pressure of about $\sim 3150$1 dynes cm$\sim 3150$2, and a mean post-shock cooling time to $\sim 3150$3 K of about $\sim 3150$4 yr (Dopita et al., 2018). The same study concluded that [Fe X] and [Fe XIV] arise in faster, partially radiative shocks around $\sim 3150$5 moving through lower-density gas near cloud peripheries, and that dust destruction is nearly complete in those faster shocks but incomplete in the slower cloud shocks until the recombination zone is reached (Dopita et al., 2018).

A later XRISM/Resolve analysis supplied direct microphysical support for ongoing ion–neutral interaction. By adding a charge-exchange component, the fit improved for weak features at $\sim 3150$6 and $\sim 3150$7 keV; the combined significance was reported as $\sim 3150$8, the charge-exchange flux was constrained to no more than $\sim 3150$9 of the total source flux in the $\sim 12$0–$\sim 12$1 keV band, and the interaction velocity was limited to $\sim 12$2 (Gu et al., 4 Apr 2025). This is consistent with the broader picture of active shock–cloud interfaces rather than a purely collisionally ionized shell.

4. X-ray plasma structure, Fe-rich ejecta, and shock dynamics

High-resolution soft X-ray spectroscopy shows that the bulk of the soft emission is multi-temperature, ionizing, and dominated by swept-up LMC gas. Deep XMM-Newton/RGS data required at least three thermal NEI components with $\sim 12$3, $\sim 12$4, and $\sim 12$5 keV, all sharing an ionization timescale of $\sim 12$6 (Suzuki et al., 2020). The inferred abundances are broadly consistent with LMC ISM values, supporting a swept-up origin for the soft component, although the forbidden-to-resonance ratios of O VII and Ne IX are higher than expected for a typical NEI plasma and were attributed to resonance scattering, charge exchange, or unresolved superposition of temperature components (Suzuki et al., 2020).

Spatially resolved Chandra spectroscopy similarly found that the rim is dominated by shocked circumstellar/interstellar material, with mean rim abundances of O, Ne, Mg, Si, S, and Fe near local LMC values, but with localized O enhancement on the north-west rim and S enhancement on the north-east blast wave (Sharda et al., 2020). Interior spectroscopy identified a relatively hot Si-rich plasma with $\sim 12$7 keV that is also responsible for the Fe K emission, while continuum-subtracted imaging showed that the Fe K emission is spread through much of the southern interior and does not reach the blast wave (Sharda et al., 2020). Deep XMM-Newton analysis later refined the hard component to a hot plasma of $\sim 12$8 keV in the CIE case, with Ca, Fe, and Ni abundances significantly enhanced over typical LMC ISM abundances, and concluded that the Fe K emission arises from ejecta heated by the reverse shock (Foster et al., 28 Apr 2025).

Microcalorimeter spectroscopy has been decisive for the kinematics of the X-ray-emitting plasma. Hitomi first showed that the Fe K complex is redshifted by about $\sim 12$9 relative to the local LMC ISM, while S XV He$\sim 25^\circ$0 is consistent with the local systemic velocity of $\sim 25^\circ$1, implying that Fe traces asymmetric ejecta whereas S traces swept-up shell gas (Collaboration et al., 2017). XRISM/Resolve then measured Si/S He$\sim 25^\circ$2 with $\sim 25^\circ$3 and $\sim 25^\circ$4, Fe He$\sim 25^\circ$5 with $\sim 25^\circ$6 and $\sim 25^\circ$7, and Fe Ly$\sim 25^\circ$8 with $\sim 25^\circ$9 and $28^\circ$0 (Collaboration, 2024). The width decomposition was written as

$28^\circ$1

and the broad Fe He$28^\circ$2 line was interpreted as the signature of very hot, recently shocked ejecta, whereas the narrow Si/S lines are consistent with shocked ISM (Collaboration, 2024). Assuming the Fe He$28^\circ$3 emission originates predominantly from ejecta, the observer-frame reverse-shock velocity at the time the bulk of the Fe ejecta were shock heated was estimated to lie in the range $28^\circ$4 (Collaboration, 2024).

Forward-shock dynamics have been inferred by two complementary methods. A thermal-diagnostic model applied to deep Chandra data yielded forward-shock velocities of $28^\circ$5–$28^\circ$6 around the rim, with moderate azimuthal dependence; these values are broadly consistent with southern proper-motion measurements but fall short of northern proper-motion speeds by as much as a factor of about $28^\circ$7, motivating interpretation in terms of highly efficient particle acceleration in the north (Okada et al., 21 Feb 2025). Direct Chandra proper-motion measurements subsequently found $28^\circ$8 for the bright southern rim and $28^\circ$9 for two features in the northeast blowout region (Long et al., 15 Sep 2025). In the south, the measured electron temperature of $4.4$0 keV is consistent with Coulomb equilibration plus adiabatic expansion after the shock, whereas the northeast temperature of $4.4$1 keV remains significantly below the value implied by the proper-motion speed (Long et al., 15 Sep 2025). This reinforces the view that N132D’s shock structure is strongly nonuniform.

5. Gamma rays and cosmic-ray acceleration

N132D is an unusually luminous extragalactic gamma-ray remnant. H.E.S.S. observations totaling $4.4$2 hr yielded a $4.4$3 detection above $4.4$4 TeV, and the joint H.E.S.S.+Fermi-LAT spectral energy distribution from $4.4$5 GeV to $4.4$6 TeV is well fit by a single power law with photon index $4.4$7 and normalization $4.4$8 at $4.4$9 TeV (Vink et al., 2021, Collaboration et al., 2021). A cutoff is not required; the best-fit exponential cutoff energy is $2450 \pm 195$00 TeV, and the quoted lower limit is $2450 \pm 195$01 TeV at $2450 \pm 195$02 confidence (Vink et al., 2021, Collaboration et al., 2021). The $2450 \pm 195$03–$2450 \pm 195$04 TeV luminosity is $2450 \pm 195$05 in one analysis and $2450 \pm 195$06 in another, making N132D more than an order of magnitude brighter than Cassiopeia A in the TeV band (Vink et al., 2021, Collaboration et al., 2021).

The multiwavelength interpretation strongly favors a hadronic origin. Inverse-Compton models are disfavored because the radio spectral index is $2450 \pm 195$07, which would imply a gamma-ray photon index near $2450 \pm 195$08, significantly harder than observed, and because N132D lacks strong nonthermal synchrotron X-rays (Vink et al., 2021). A leptonic fit requires a magnetic field of about $2450 \pm 195$09 and an electron energy budget of order $2450 \pm 195$10 erg, which was judged implausibly large (Vink et al., 2021). By contrast, hadronic modeling gives

$2450 \pm 195$11

corresponding to roughly $2450 \pm 195$12 of the explosion energy for $2450 \pm 195$13, and the unbroken gamma-ray spectrum implies parent protons extending to at least $2450 \pm 195$14 TeV (Vink et al., 2021). A separate H.E.S.S.-based hadronic fit used a proton cutoff energy of $2450 \pm 195$15 TeV and derived a lower limit of $2450 \pm 195$16 TeV for the parent proton distribution (Collaboration et al., 2021).

The gamma-ray phenomenology has also been linked directly to the cloud environment. ALMA/ASTE CO work argued that if the gamma rays are hadronic, the shock-engulfed molecular clouds inside the remnant provide target material for cosmic-ray protons, and using molecular-cloud densities of $2450 \pm 195$17–$2450 \pm 195$18 it estimated a conservative lower limit of $2450 \pm 195$19–$2450 \pm 195$20 erg for the accelerated proton population (Suzuki et al., 2020). A broader H.E.S.S. study considered three possible sites for the gamma rays: the shell interacting with dense H I, shocked dense cloudlets inside the remnant, and escaped cosmic rays illuminating nearby CO clouds to the southwest; the measured southwest offset of the gamma-ray centroid was noted but deemed not statistically significant (Vink et al., 2021). This suggests that the hadronic interpretation is robust, whereas the exact spatial partition between shell and cloud targets remains unresolved.

Earlier NuSTAR and Suzaku work had already framed N132D as a remnant in transition from young to middle-aged. That study found a tight upper limit on nonthermal X-rays, a likely recombining or very hot thermal plasma, and a proton cutoff energy of about $2450 \pm 195$21 TeV in the broadband spectral energy distribution, again favoring hadronic gamma rays over a leptonic origin (Bamba et al., 2018). Subsequent H.E.S.S. results extended that picture by showing that the gamma-ray spectrum remains a single power law beyond $2450 \pm 195$22 TeV, giving the remnant a “youthful” high-energy spectrum despite its mature age (Vink et al., 2021).

6. Progenitor, explosion geometry, and outstanding controversies

The progenitor mass of N132D remains actively debated because different analyses emphasize different tracers and different parts of the ejecta. Far-ultraviolet HST/COS spectroscopy of an O-rich knot found helium and silicon associated with the oxygen-rich material, no nitrogen, and a preferred abundance ratio $2450 \pm 195$23 for a shock velocity of about $2450 \pm 195$24; comparison with nucleosynthesis models then suggested a progenitor around $2450 \pm 195$25 (0909.4110). Infrared ejecta and dust analysis argued that the dust mass and ejecta properties are more consistent with a higher progenitor mass of about $2450 \pm 195$26 (Rho et al., 2023). By contrast, Chandra abundance ratios and cavity arguments yielded $2450 \pm 195$27, and deep XMM-Newton analysis of the hot Fe-rich ejecta component used the Ca/Fe and Ni/Fe mass ratios to infer $2450 \pm 195$28 (Sharda et al., 2020, Foster et al., 28 Apr 2025). The divergence of these estimates indicates that the progenitor problem is not settled.

The explosion geometry is similarly contested. One interpretive line treats the main O-rich structure as an equatorial ring or torus, perhaps produced by a bipolar explosion from a rotating progenitor, with the overall remnant shape strongly influenced by pre-supernova mass loss and with the runaway knot representing a possible polar jet-like feature (Vogt et al., 2010). A later reconstruction emphasized the similarity to Cassiopeia A: both remnants show a torus-like main shell, a rupture or break in that structure, and a high-velocity Si-enriched component, leading to the more cautious conclusion that both explosion dynamics and environmental interaction likely matter and that the balance between them remains unresolved (Law et al., 2020).

A more specific 2025 proposal attributed the remnant’s morphology to the jittering jets explosion mechanism. That work identified a point-symmetric morphology with two approximately perpendicular symmetry axes, one running from the northwest ear through the center of the iron-rich emission and the other following the long axis of the remnant and the northeastern extension of the central oxygen-rich gas (Soker, 1 Jul 2025). Within that framework, one jet pair is invoked to shape the O-rich torus and runaway knot, and a second pair is invoked to shape the large-scale elongation (Soker, 1 Jul 2025). This is a model-specific interpretation rather than a consensus result, but it illustrates how strongly N132D’s asymmetries constrain explosion models.

Despite these controversies, several points are comparatively secure. N132D preserves a large O-rich ejecta structure with long-lived ballistic motion, contains asymmetric Fe-rich ejecta that are redshifted with respect to the local ISM, expands into a cavity or wind-blown bubble with strong density contrasts, and is interacting with dense clouds that materially affect its morphology and high-energy output (Banovetz et al., 2023, Collaboration, 2024, Sano et al., 2015, Long et al., 15 Sep 2025). The unresolved issues concern weighting rather than existence: how much of the present structure is fossil explosion geometry, how much is reverse-shock processing, how much is imposed by the surrounding medium, and which ejecta component best samples the progenitor’s true nucleosynthetic yield.