Ejecta, Rings, and Dust in SN 1987A with JWST MIRI/MRS

Supernova ( SN ) 1987A is the nearest supernova in ∼ 400 yr. Using the JWST MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring ( ER ) , and outer rings in the mid-infrared 12,927 days ( 35.4 yr ) after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. The broad emission lines ( 280 – 380 km s − 1 FWHM ) that are seen from all singly-ionized species originate from the expanding ER, with properties consistent with dense post-shock cooling gas. Narrower emission lines ( 100 – 170 km s − 1 FWHM ) are seen from species originating from a more extended lower-density component whose high ionization may have been produced by shocks progressing through the ER or by the UV radiation pulse associated with the original supernova event. The asymmetric east – west dust emission in the ER has continued to fade, with constant temperature, signifying a reduction in dust mass. Small grains in the ER are preferentially destroyed, with larger grains from the progenitor surviving the transition from SN into SNR. The ER dust is ﬁ t with a single set of optical constants, eliminating the need for a secondary featureless hot dust component. We ﬁ nd several broad ejecta emission lines from [ Ne II ] , [ Ar II ] , [ Fe II ] , and [ Ni II ] . With the exception of [ Fe II ] 25.99 μ m, these all originate from the ejecta close to the ring and are likely to be excited by X-rays from the interaction. The [ Fe II ] 5.34 to 25.99 μ m line ratio indicates a temperature of only a few hundred K in the inner core, which is consistent with being powered by

1. INTRODUCTION Supernovae (SNe) are important contributors to the chemical evolution of the universe over cosmic time, however, their relative net contribution to dust production in the early universe is still debated.While extant studies of SNe have shown that these events are prodigious sources for dust production (e.g.Matsuura et al. 2011;Niculescu-Duvaz et al. 2022;Shahbandeh et al. 2023), they have found that their ability to also destroy dust is similarly (and potentially equally) proficient (e.g., Gall & Hjorth 2018;Slavin et al. 2020;Kirchschlager et al. 2023).Assessing the effects of SNe on their host galaxies, and understanding how their dust production and destruction processes occurs over months to years, requires time-series observations at a high angular resolution to build realistic models.
SN 1987A provides a once-in-a-lifetime opportunity to investigate the evolution of a supernova into a supernova remnant (SNR;McCray & Fransson 2016).On February 23, 1987, light from a massive blue supergiant (BSG) Sanduleak -69 202 in the Large Magellanic Cloud (LMC), ∼50 kpc from us (Pietrzyński et al. 2019), which collapsed and exploded as a supernova was visible to the naked eye (e.g., Arnett et al. 1989); producing an initial burst of neutrinos consistent with the formation of a neutron star (Burrows 1988).In the intervening years, SN 1987A has been extensively studied across the electromagnetic spectrum revealing important insights into the physics of supernovae, the properties of stellar explosions, and the evolution of stars.See Table 1 for a listing of some key events in the evolution of SN 1987A.
The SN 1987A morphology is comprised of asymmetrical inner ejecta, a dense equatorial ring (ER) approximately 2 arcsec in diameter inclined at 43 • from the line-of-sight, and two faint outer rings 5 arcsec in diameter.The outer rings, thought to be produced at the same time as the ER, have an expansion velocity of 26 km s −1 (Crotts & Heathcote 2000), whilst the ejecta is expanding at a rate of ∼ 2000-10000 km s −1 (McCray 1993;Fransson et al. 2013;Larsson et al. 2016;Kangas et al. 2022b).
Emission from the ER is the dominant component of SN 1987A observations from the X-ray to the midinfrared (IR).The ER was initially photoionized by the SN shock breakout, but that emission has faded away and the emission from the ER is now dominated by shock interaction.The ER was likely produced in a binary merger event during the progenitor's red supergiant phase of evolution (Morris & Podsiadlowski 2009) and was initially expanding at 10.3 km s −1 relative to the SN.As shocks from the blast wave interacted with the dense gas in the ER, ∼30 hot spots began to form and steadily brighten in the ring in optical wavelengths (Pun et al. 1995;Bouchet et al. 2000;Fransson et al. 2015;McCray & Fransson 2016), increasing the expansion velocity of the ring to ∼700 km s −1 (Larsson et al. 2019a;Kangas et al. 2022a).Once the initial shock from the blast wave passed (at ∼7,500 days, or 20.5 yrs) these hot spots began to fade and there is increasing interaction outside the ER as the forward shock expands outwards.Interactions with reverse shocks propagating into the SN ejecta also affect recent processing, gradually destroying the ring in the process (Fransson et al. 2015;Larsson et al. 2019a;Orlando et al. 2019;Arendt et al. 2020;Kangas et al. 2022b).
Using the T-ReCS instrument, Bouchet et al. (2004) resolved 10 µm emission from SN 1987A at day 6067 (year 16.6) tracing its warm dust morphology; this emission is dominated by the ER, with less than 10 −5 M ⊙ of warm progenitor dust.Spitzer InfraRed Spectrograph (IRS) observations revealed the ER is comprised of warm (T d ≈ 180 K) silicate dust (Bouchet et al. 2006;Dwek et al. 2010;Arendt et al. 2016).This emission brightened up to day 8708 (year 23.8) as the dust grains were collisionally heated by the shocked gas, however, the shape of the dust spectrum remained unchanged.
Cryogenic Spitzer IRAC fluxes (obtained from day 6130, year 16.8) appeared to be strengthening in agreement with the IRS spectra, however, a secondary hot T d > 350 K dust component of undetermined composition appeared to be necessary to fit the blue-end (3 -8 µm) of the spectral energy distribution (SED; Dwek et al. 2010).Imaging with the warm Spitzer mission continued to monitor the SN 1987A light curves until day 11,885 (year 32.5) (Arendt et al. 2020).These appear to peak between day 8500-9000 (year 23.2-24.6)before undergoing a gradual fading in luminosity since day 9810 (year 26.9) post-explosion, which is consistent with the passage of the forward shock through the ER and the gradual destruction of dust emitting at these shorter wavelengths.
In the IR, CO and SiO vibration bands from the ejecta became apparent in the integrated spectra of SN 1987A from day 200 post-explosion (Roche et al. 1991;Spyromilio et al. 1988), and dust formation in the ejecta occurred between 400 and 500 days, consequently obscuring some of the optical and near-IR emission and re-radiating this light in the mid-IR (Lucy et al. 1989;Bouchet et al. 1991;Wooden et al. 1993).This dust was characterised by a featureless single-temperature graybody emission spectrum.
The Herschel far-IR observations (Matsuura et al. 2011) unveiled a large reservoir of newly condensed cold dust with temperature T d = 17 − 23 K and mass M d = 0.3 − 0.8 M ⊙ , depending on whether its composition arises solely from amorphous carbon or a mix of amorphous carbon and silicates (Matsuura et al. 2015).The origin of this dust was unambiguously confirmed to be from the inner ejecta of SN 1987A by the Atacama Large Millimetre Array (ALMA; Indebetouw et al. 2014;Matsuura et al. 2017).This extremely efficient dust condensation in the ejecta (measured 24 years after the explosion) requires almost all the expected Mg and Si produced in the SN to condense into grains and is significantly higher than the dust mass measured by Bouchet et al. (1991) and Wooden et al. (1993) in the first few years post-explosion.This ejecta dust, 36 years after the explosion, has not yet experienced processing by the reverse shock (France et al. 2010;Chevalier & Fransson 2017); it remains unclear how much of this dust will survive and dissipate into the interstellar medium once this interaction occurs.
The SN 1987A ejecta has a complex asymmetric morphology.In the optical, the ejecta is elongated northsouth with a faint central region (Larsson et al. 2011), and is currently undergoing a free expansion interior to the ER (Larsson et al. 2016(Larsson et al. , 2019b)).ALMA observations show the CO and SiO molecules in the ejecta have different distributions, with cold dust emission coincident with the faint optical region (Kamenetzky et al. 2013;Indebetouw et al. 2014).These observations agree with theoretical models of dust formation in SN that show CO and SiO form in different regions of the ejecta a few hundred days after explosion (Sarangi & Cherchneff 2013;Sluder et al. 2018).
Several near-and mid-IR lines from the ejecta were observed during the first years after explosion with both ground-based instruments and the Kuiper Airborne Observatory.These have proved to be extremely useful as diagnostics of the temperature, density and abundances in the ejecta (Aitken et al. 1988;Rank et al. 1988;Spyromilio et al. 1990;Wooden et al. 1993;Colgan et al. 1994).In particular, the line ratios of the [Fe II] lines indicated a low filling factor for the iron in the core of the ejecta (Li et al. 1993;Kozma & Fransson 1998), indicative of 56 Ni heating, during the first weeks after explosion, creating low-density bubbles in the core.Also the [Fe II] line profiles were very useful for probing the asymmmetry and abundance mixing of the ejecta (Spyromilio et al. 1990;Haas et al. 1990).These observations have provided some of the most important constraints on the explosion mechanism.
The launch of JWST (Rigby et al. 2023;Gardner et al. 2023) with its superb sensitivity and high spatial resolution images and integral field spectrographs (IFS) has opened a new era in mid-IR observations, enabling us to obtain IR data on SN 1987A at a comparable resolution to the optical HST data and sub-mm ALMA data.For instance, NIRSpec observations reveal asymmetry in the ejecta, likely due to a binary merger prior to explosion and morphological 3D structures in the circumstellar medium (Larsson et al. 2023).
With the Mid-Infrared Instrument (MIRI; Wright et al. 2023) Medium Resolution Spectrometer (MRS) onboard JWST, we can finally observe and spatially resolve the ER and ejecta of SN 1987A at mid-IR wavelengths with a high S/N ratio.The paper is organised as follows.In Section 2 we describe the observations and data reduction.Our results are presented in Section 3. Finally, in Section 4 we discuss the implications of these data.The conclusions are summarised in Section 5.The optical emission from the ER starts fading, the blast wave has passed the ER.
The ER fades at most wavelengths.2013 9500 First detection of optical spots outside the ER.
provided complete spectral coverage from 4.9 to 27.9 µm at medium resolution (R ∼ 4000-1500; Jones et al. 2023).As no wavelength channel was prioritised for the SN 1987A MRS observations, a 4-point extended source dither pattern in the 'negative' direction optimised for 'ALL' channels was selected.The MRS fields of view range from 3.2 ′′ × 3.7 ′′ in channel 1 to 6.6 ′′ × 7.7 ′′ in channel 4 (Wells et al. 2015;Argyriou et al. 2023).No target acquisition was employed and no dedicated background observations were taken.The total MRS integration time corresponded to 9396s, or approximately 3.75 hr including all overheads.Simultaneous MIRI imaging in the F560W, F770W, and F1000W filters was obtained to improve the astrometric accuracy of the MRS data.These observations were taken using the slowr1 readout pattern (to reduce data volume), with 10 groups and 3 integrations, using the same dither pattern as the MRS data.
To reduce our MRS data we used a development version of the JWST Calibration Pipeline (Bushouse et al. 2023) retrieved on May-23-20231 , with versions 11.17.0 and 'jwst 1090.pmap' of the Calibration Reference Data System (CRDS) and CRDS context, respectively, and development versions of the FLT-5 MRS calibration files which include an improved photometric calibration derived from the Cycle 1 calibration program.We ran all raw data files through the Detector1Pipeline to convert the ramp files into rate images.Since inaccuracies in the spacecraft pointing information can be introduced by guide star catalogue errors and roll uncertainty (see Pontoppidan et al. 2022), we corrected the world coordinate system (WCS) reference keywords in our rate im-ages by determining the offset of point sources detected in the F560W simultaneous imaging field to their Gaia Data Release 32 counterparts.This resulted in an offset correction of −2.825 ′′ in RA and −0.630 ′′ in Dec.This correction is expected to be accurate to ≲0.1 ′′ (Patapis et al. 2023), with the upper limit due to uncertainties in the relative astrometry between the MRS and imager.The pointing correction must be applied before running stage 2 of the pipeline as the WCS solution is attached to the files in its first step.
We ran all files through Spec2Pipeline, with the optional residual fringe step switched on to produce flux calibrated rate images, from which we constructed spectral cubes for all of the 12 MRS sub-bands using Spec3Pipeline.MRS data usually contain large numbers of 'warm' pixels which are not flagged in the bad pixel mask as they evolve over time.For data where no dedicated background is available, these can have a significant impact on cubes and extracted spectra, where they manifest as spurious emission lines.To account for these when producing our cubes, we used the updated outlier detection step in Spec3Pipeline to capture and omit these pixels from the cube-building process.We show sample slices from each of the 12 sub-band cubes in Fig. 1 to highlight some features of our MRS data.
Channel 1 and 2 cubes suffer from a 'striping' in the x-dimension of the cube frame, which is residual dark current.Even though the pipeline subtracts a dark current reference from the data, drifts in the effective dark current over time cause such residuals.Unfortunately, without dedicated background observations, this resid- ual cannot be removed without risking the subtraction of real signal from SN1987A.The channel 4 cubes suffer from an increased signal towards the right edge in the cube frame.This likely results from a more uncertain flux calibration near the slice edges.We omitted these regions from our analysis.
Two final aspects of our data to contend with were the degradation in spatial resolution and increasing background from band 1A through 4C, both of which are evident in Fig. 1.The former presents challenges in creating representative full-band MRS spectra from small regions of the SN 1987A system.These are discussed in more detail below.As we had no dedicated background pointing we were required to use regions within the sub-band FOVs to create a background that could be subtracted from our spectra.We defined circular regions surrounding SN 1987A for channels 1-3, with the backgrounds extracted from all pixels exterior to the circle.For channel 4 we used a polygonal region surrounding SN 1987A but omitting the bright edge of the cubes mentioned previously.These background regions are shown in Fig. 2, along with the location of the outer rings defined using their morphology in the [Ne ii] 12.81354 µm line.It is clear from Fig. 2 that the outer rings take up a significant fraction of the channel 1 and 2 FOVs making it impractical to have a 'true' background.The sizes of background boundaries in these channels are therefore a trade-off between producing a representative background while mitigating contamination from the line emission from outer rings (see Sect. 3.1).
Due to the complex morphology of SN 1987A, we have extracted spectra from several locations in the MRS cube.These include 1) the total ER; 2) segments of equal size from the east and west sides of the ER; 3) four smaller cardinal point regions from the N, S, E and W portions of the ER, comparable in size to the hot-spots observed in the optical, 4); the central ejecta, and finally 5); the entire SN 1987A system to enable comparisons with past Spitzer observations3 .All regions are shown in Fig. 2. Spectral regions were defined using SAO Image DS9 (Joye & Mandel 2003) and extracted using the aperture module in the Astropy Photutils package4 (Bradley et al. 2022a) by parsing the regions using the Astropy Regions package5 (Bradley et al. 2022b).We applied an additional post-pipeline residual fringe correction to our spectra to account for the high-frequency fringing in channels 3 and 4, thought to originate in the MIRI/MRS dichroics, which are not effectively removed in the pipeline.This correction is now included in the JWST Calibration Pipeline package under the extract 1d step.Finally, to produce the total, full-band MRS spectrum of SN 1987A we scaled the MRS sub-band fluxes for minor offsets by determining the median flux difference in the overlapping wavelength regions, starting from the band 1A.These subband offsets were typically at the ∼ 1% level, with the worst-case offset of ∼ 8% being at the band 2C/3A overlap region where MRS fringing is most problematic.The scaled sub-band spectra were stitched using the combine 1d tool in the JWST Calibration Pipeline package.With the MRS we have detected and spatially resolved the 5-28 µm emission from the ejecta, ER and outer rings of SN 1987A at day 12,927 (year 35.4) after the explosion.The ER is well resolved and isolated from the nearby Stars 2 (which is not visible in the mid-IR) and 3 (Figure 1) which have affected lower spatial resolution observations of the system, especially at shorter wavelengths.The ER produces most of the IR emission, however, this emission has a complex morphology and is not evenly distributed over the ring.It is brightest on the west side.This is consistent with the evolution of the peak hot spot emission observed by HST in the ER from the northeast to the west side circa 2006 (Fransson et al. 2016;Larsson et al. 2019a), due to the asymmetric distribution of matter in the ER.Similar behaviours are also observed in the X-ray (Frank et al. 2016) and near-IR (Kangas et al. 2022a;Larsson et al. 2023).
The line emission from the ER originates from dense clumps of gas, observed as hotspots in the optical images, while we find that the continuum emission is more spatially extended.This is illustrated for the H I 7.460 µm Pfund-α line and nearby continuum in Figure 3.In order to quantify the differences in the spatial distribution, we fit the ER with an elliptical model with a Gaussian radial profile and sinusoidal intensity in the azimuthal direction.While this model does not account for all of the substructure in the ER, it captures the bulk of the emission and hence provides a good estimate of the overall size.The semi-major (minor) axes of the best-fit ellipses are 0. ′′ 850±0.′′ 006 (0. ′′ 632±0.′′ 002) for the line image, compared to 0. ′′ 988 ± 0. ′′ 004 (0. ′′ 721 ± 0. ′′ 002) for the continuum image.
In the collapsed MIRI/MRS data cubes the overall morphology is comparable to the T-ReCS 10 µm data observed by Bouchet et al. (2004Bouchet et al. ( , 2006) ) at days 7241, 7565, 8720 and VLT/VISIR 10-18 µm images at days 10950-10976 (year 30.0) (Matsuura et al. 2022).Higher spatial resolution JWST mid-IR images of SN 1987A obtained on the same day post-explosion as the MIRI/MRS data including the ER morphology are discussed in depth in Bouchet et al. (2024).Strikingly the resolved emission from the [Ne ii] lines reveals the outer ring structure (see Sect. Figure 4 shows the total background-subtracted MRS spectrum of SN 1987A, for the region defined by the white annulus shown in the upper part of Fig.   dust.These spectral features and their interpretation are discussed in detail in the following sections.

Emission line spectral properties
To investigate the emission line properties of the ER and its adjacent regions, we analyzed the spectrum corresponding to the whole ring, shown in Figure 4, as well as the spectra extracted for the four cardinal point subregions of the ER that are shown in the lower part of Fig. 2. Emission line radial velocities, full-width halfmaximum (FWHM) values and integrated line fluxes were measured using the Gaussian emission line fitting (ELF) routines written by P. J. Storey within the dipso spectral analysis package (Howarth et al. 2014).Results are presented in Table 2 for the whole ER spectrum and in Table 3 for the four cardinal point spectra.The solid angle corresponding to each spectral extraction is listed in the table headers.
Sixteen emission line detections are listed in each table.H i 8-6 7.502 µm was also detected.In addition, a very broad 26-µm feature with multiple peaks (FWZI∼8000 km s −1 ) is present in each of the spectra, corresponding to [Fe ii] 25.98839 µm, superposed on what is a narrow [O iv] 25.8903 µm line (see Fig. 4).A weaker blueshifted component to this feature is discussed in Section 3.4 and shown in Figure 15.For the whole ER spectrum the integrated [Fe ii] 25.99 µm line flux is measured to be (2.3 ± 0.1) × 10 −21 W cm −2 , similar to the flux levels measured for this line from a time sequence of Spitzer spectra obtained up to day 7954 (year 21.8) (Arendt et al. 2016).A broad redshifted emission line corresponding to [Fe ii] 5.340169 µm is also present in the whole ER spectrum, with a flux of (6.0 ± 0.25) × 10 −22 W cm −2 , a heliocentric line centre radial velocity of +2450±30 km s −1 and an FWHM of 1600±70 km s −1 .The only cardinal point spectrum in which the 5.34 µm line is detected is the Ring-South spectrum (flux=(6.2± 0.13) × 10 −23 W cm −2 , velocity= +2990 ± 11 km s −1 , FWHM= 1055 ± 27 km s −1 ).These lines are due to the ejecta-ring interaction and are discussed further in Section 3.4.
A comparison of the MRS line fluxes for the whole ER listed in Table 2 with  The 10-40 µm Spitzer-IRS spectra discussed by Arendt et al. (2016) had a resolving power of R=600.The ∼5 times greater resolving power of the 5-27 µm JWST-MRS spectra allows higher resolution FWHM and radial velocity measurements to be made.Inspection of the FWHM values for the lines in the whole ER spectrum listed in Table 2 reveals two groupings: (a) eight lines with FWHM values of between 280 and 380 km s −1 , namely all the lines that arise from singly ionized species (Ni + , Ar + , H + , Ne + and Fe + ), together with lines from doubly ionized Na 2+ and Ar 2+ ; and (b) seven lines arising mainly from more highly ionized species (Ne 5+ , Ne 4+ , Ne 2+ , S 3+ , S 2+ and O 3+ ), with FWHM values between 99 and 171 km s −1 .The [S iv] 10.51 µm, [Ne vi] 7.65 µm and [Ne v] 14.32 µm lines are particularly narrow (99-107 km s −1 FWHM)not much larger than the instrumental resolution measured at those wavelengths by Jones et al. (2023), while the larger line widths (123-145 km s −1 ) measured for the [S iii] 18.71 µm, [Ne v] 24.32 µm and [O iv] 24.89 µm lines are consistent with the instrumental resolutions measured by Jones et al. (2023) for this lower resolution part of the MRS range.
The general expansion of the ER may contribute to the broadening of some of the lines seen in the whole-ER spectrum, an effect which should be reduced in the individual cardinal point spectra.Comparison of the FWHM values measured for the whole ER (Table 2) with those measured for the cardinal point spectra (Table 3) shows that for the singly ionized species their FWHM values are appreciably lower in the cardinal point spectra.However, this is not the case for the much narrower lines that arise from the more highly ionized species, for which the whole-ER and cardinal point FWHM values are similar.This behaviour indicates that while the emission from the singly and doubly ionized species originates from the ER itself, that from the more highly ionized species may have a partly different origin.This is also shown by the line maps presented in Section 3.1.2.Support for this interpretation comes from a consideration of the radial velocity differences seen in Table 3 across the four cardinal point spectra.For lines from singly ionized species, the mean difference between the radial velocities measured for the Ring-South and Ring-North spectra was +286.7±45.5 km s −1 , while the East-West mean difference was +46.9±29.5 km s −1 .In contrast, for the more highly ionized species, the corresponding radial velocity differences were much smaller, +45.7±49.1 km s −1 (South-North) and −7.1±10.0km s −1 (East-West).We take this as confirmation that the lines from singly and doubly ionized species originate principally from the expanding ER and that the much narrower lines from more ionized species must originate from a different emission component.Gröningsson et al. (2008) presented 6 km s −1 resolution optical spectra of northern and southern parts of the ER, obtained in 2002 at days 5702-5705 (year 15.6) past explosion, at a time when the northern quadrant was the brightest part of the ER.They saw a mixture of narrow lines (FWHM=10-29 km s −1 ), from H i up to [Ne v], and broad lines (line widths larger than 200 km s −1 ), from H i up to [Fe xiv].They attributed the narrow lines to unshocked pre-ionized gas in the ER and the broad lines to recombining post-shock gas.By the time of our MIRI-MRS observations of the ER on day 12927 (year 35.4), the only broad emission lines seen were from the singly ionized species, consistent with the ongoing recombination of the post-shock gas and with Gröningsson et al. (2008)'s prediction that "As more and more gas cools we expect the width of especially the low ionization lines to increase." From optical emission line diagnostics and thermally broadened line widths, Gröningsson et al. (2008) estimated electron densities n e ∼ 1500-5000 cm −3 and electron temperatures T e ∼ 6500-24,000 K for their narrowline component.The narrow-line [Ne v] 14.32/24.32µm line flux ratio is the only electron density diagnostic available from our MRS spectra of the ER.For the Ring-West spectrum, the brightest of the four cardinal point spectra (see Table 3), this ratio is measured to be 1.151±0.115,while for the whole-ER spectrum (see Table 2) the ratio is found to be 1.021±0.059.Using the atomic data of Galavis et al. (1997) and Griffin & Badnell (2000), these ratios correspond to electron densities n e of between 700 and 4300 cm −3 for electron temperatures in the range T e = 5000 − 25000 K, similar to the parameters derived by Gröningsson et al. (2008) for their narrow line component.
The detection of two [Ni ii] lines in the Ring-West spectrum allows a constraint to be placed on the electron temperature of this part of the ER.Using the atomic data of Nussbaumer & Storey (1982), Quinet & Le Dourneuf (1996) and Bautista (2004), the observed [Ni ii] 6.64/10.68µm flux ratio of 9.71±0.35requires an electron temperature T e ≤ 6500 K.For temperatures higher than this the 10.68 µm line (which originates from a higher level than the upper level of the 6.64 µm Table 3. Detected emission lines in the four cardinal point spectra, their measured heliocentric radial velocities, full-width half maxima, and line fluxes.transition) would become too strong.This temperature limit is consistent with strong post-shock cooling having taken place in the dense ER clumps.

Line emission maps
Here, we discuss the spatially resolved line emission structure of SN 1987A over the MRS IFS for the different narrow-and broad-line-emitting species.For each selected emission line, we formed continuum-subtracted images by subtracting from each line-centred image a mean image from adjacent line-free regions on either side of the line, each off-image having the same spectral bandwidth as the line-centred image.In order to check for radial velocity variations across this structure that might support a very large expanding ring interpretation, we extracted spectra of the [S iii] 18.71 µm line at three positions in the ring-like structure: one from near its north-eastern edge; one near its north-western edge; and one just to the north of the ER.However, the measured [S iii] radial velocities showed a difference of only 4.9±0.7 km s −1 between the NE and NW positions, with the velocity of the position north of the ER measured to be in between those of the other two.We interpret these results as not providing strong support for a very large expanding ring interpretation.
We have identified spectral features of both small and large spatial scales by manually searching through the spectral axis MRS sub-band datacube using small and large apertures at varying locations around the system.We then identified the emitting species associated with each feature before slicing a sub-cube around the feature using the Spectral-Cube 6 Astropy affiliated package (Ginsburg et al. 2014).We carefully fitted and subtracted the continuum from each of the sub-cube spaxels using regions away from the spectral feature and 6 https://spectral-cube.readthedocs.io/en/latest/accounting for the variation in sampling and residual fringing.With the resulting sub-cubes, we created sets of moment maps using Spectral-Cube.From the 0th, 1st, and 2nd moments, we produced integrated line intensity maps, velocity maps which account for the systemic barycentric velocity of SN 1987A of 287 km s −1 (Gröningsson et al. 2008), and FWHM maps.We applied a cut across all maps based on an integrated intensity threshold to mask regions with little or no line emission.The moment maps for a selection of lines are shown in Figs. 6, 7 and 8. Lines from singly ionized species, the majority of those plotted in the three figures, all show a strong radial velocity difference between the blue-shifted northern part of the ER and the red-shifted southern part, consistent with the cardinal point radial velocity measurements for singly ionized species listed in Table 3, and with previous ground-based results (e.g.Kjaer et al. 2007).
The radial velocity map shown in Fig. 8 for [Ne v] 14.32 µm, a representative of the more highly ionized species, does not show any radial velocity changes across the ER, consistent with its velocity measurements in Table 3, which show differences of only 2.0±1.8 and 13.6±2.5 km s −1 between its East-West and South-North radial velocities, respectively.But Table 3 also shows that the fluxes measured for [Ne v] and other high ionization species peak in the Ring-West aperture, indicating a relationship between the narrow-line component and the ER material.However, the narrow-line component also extends to the outer rings and beyond, up to the edges of the MRS FoV -see Figs. 5 and 8 for e.g.[Ne v], with the mean surface brightness of the [Ne v] 14.32-µm line in the background aperture (Fig. 2) corresponding to 4 per cent of its surface brightness in the Ring-West aperture.

Dust Composition
The MRS spectrum of the ER, shown in Figure 4, is dominated by continuum emission from amorphous silicate dust, characterised by the broad ∼10 µm and ∼18 µm spectral features.No evidence for polycyclic aromatic hydrocarbons (PAH) emission is present in the spectrum; the improved resolution of the MRS has not revealed any new solid-state dust species in the mid-IR spectra compared to Spitzer/IRS observations.
In order to determine which dust properties best reproduce the shape of the dust continuum observed with JWST, we fitted the 5.3-28 µm MRS and the 0.9-5.3µm NIRSpec data (Larsson et al. 2023) with two-temperature dust models of varying grain compositions.We accounted for the synchrotron emission for  .day 12,927 using the model constructed from ALMA observations of the ER that is described by Equation 1 of Cendes et al. (2018).We used the parameters from the torus model listed in Table 6 of Cendes et al. (2018) and the updated value of the spectral index α of 0.7 from Cigan et al. (2019).The parameters for the synchrotron component were fixed and did not vary in our fit.We also included bound-free and free-free emission from H I, He I, and He II with a temperature of 10 4 K that was used by Larsson et al. (2023) to fit the NIRSpec data.
The models are shown in Figure 9 and the best-fit temperatures and corresponding dust masses are listed in Table 4.We find that the 0.9-28 µm spectral range is bestreproduced by a single grain composition of astrodust  2023), astrodust reproduces the observational properties of dust in the diffuse ISM and is made up of composite grains that consist of different compositions (primarily silicates and oxides) on small scales but are effectively one uniform composition on larger scales of >0.05 µm.We assumed that the grains are oblate spheroids with an axial ratio of 1.4:1, a porosity of 0.2, and a density of 2.74 g cm −3 (see Hensley & Draine 2023).Since the mass absorption coefficients are sizeindependent in the mid-IR wavelength range, the choice of grain size did not affect our best-fit parameters.
We also attempted to fit the spectra using a silicate dust grain composition (Draine & Lee 1984), but as was found in earlier studies (e.g.Dwek et al. 2010), we required an additional featureless dust component to account for the shorter wavelength emission.For this sec-ond component, we used carbon dust from Draine & Lee (1984) (middle panel of Figure 9) and pure iron dust from Ordal et al. (1985) (bottom panel).Both combinations fit the spectra fairly well, even though the fits are slightly poorer than the astrodust in the 3-5 µm range.The corresponding dust masses and temperatures are listed in Table 4.The total dust mass required from the pure Fe dust component is unreasonably high, so we disfavour this composition as a dominant component producing the shorter-wavelength dust continuum.
The appropriateness of the astrodust ISM dust model of Draine & Hensley (2021) for the ER dust around SN 1987A is a valid issue, given that its immediate progenitor was a B supergiant.While 'normal' B supergiants do not form dust in their outflows (Barlow & Cohen 1977), some luminous blue variables with B supergiant spectra have circumstellar dust shells that exhibit mid-IR silicate emission bands (e.g.Voors et al. 2000), but existing observations lack sufficient signal-to-noise for comparisons to be made to the astrodust model.Another possibility is that the ER dust might have been produced during an earlier M supergiant phase of the progenitor star.The mid-IR spectral survey of M supergiants by Speck et al. (2000) found a wide variety of 10-µm silicate feature shapes, including a 'broad feature' group whose 10-µm profiles exhibit similarities to that of SN 1987A's ER.The need for additional opacity shortwards of 7-µm in cool star dust envelope models, over and above that provided by pure silicates, has long been recognised, leading to the introduction of 'dirty silicate' models for such sources (e.g.Jones & Merrill 1976).
Out of the grain compositions that we used, only the astrodust composition could reproduce the observed spectrum without the need for an additional featureless grain composition.The reason for this can be seen in the comparison of the mass absorption coefficient (κ) values for each of the grain species, shown in Figure 10.It can be seen that the astrodust opacity is significantly larger in the 2.5-7.0 µm range, and the peak-to-continuum ratio for the broad 9.8 and 20 µm features lower, compared to silicates from Draine & Lee (1984) and the LMC carbonaceous-silicate grains from Weingartner & Draine (2001).

Temporal Evolution
The comparison between the Spitzer/IRS spectra from previous epochs to the day 12,927 (year 35.4) total MIRI/MRS spectra extracted from SN 1987A is shown in the left panel of Figure 11.The Spitzer IRS data from days 6000 to 8000 (years 16.4 to 21.9) have been reduced Note-Dust temperatures and masses for the models shown in Figure 9.The listed uncertainties are the 3-σ uncertainties from the fit.  in a consistent manner using the optimal extraction from the Combined Atlas of Sources with Spitzer IRS Spectra (CASSIS; Lebouteiller et al. 2011), with the Short-Low scaled to match the Long-Low flux following Dwek et al. (2010).The MRS spectrum is obtained from within an elliptical region just inside the background boundaries (see Fig. 2) and includes all emission from the ER and inner region and is thus comparable to the Spitzer IRS data in which the SN 1987A system is not spatially resolved.It can be seen that since day ∼8000 (year 21.9), the brightness of the broad 20 µm silicate feature has continued to increase whereas the part below 15 µm has decreased.This effect is clearly illustrated in the right panel of Figure 11 where the spectra have been normal- ized to the flux density of the 20-21 µm range of the day 6805 (year 18.6) spectrum.As was shown by Dwek et al. (2010), the Spitzer spectra (grey curves) have an almost identical spectral shape.However, while the day 12927 (year 35.4) MRS spectrum shows the same approximate shape above 18 µm (the spectral region that samples the cold dust component), the relative flux densities at shorter wavelengths are significantly lower.
In order to determine the evolution of the dust parameters over time, we fitted all spectra shown in Figure 9 with a two-temperature dust model using the astrodust composition.For a consistent comparison between Spitzer IRS and MRS, we only fitted the 5-28 µm range common to both instruments.We included the previously discussed time-varying synchrotron component (Cendes et al. 2018) in the model.The resulting dust temperatures and masses are summarized in Table 5 and shown in Figure 12.Interestingly, the dust temperature has remained the same within the uncertainties across all epochs.The cold dust mass measured from the most recent MRS data has increased slightly since day ∼8000 (year 21.9), but the mass of the hot dust component dropped by more than 40% (right panel of Figure 12).This is likely indicative of significant destruction of the small grain component over time, as previously noted by Arendt et al. (2020) based on Spitzer IRAC imaging observations.The temporal evolution does not occur in a uniform way over the entire ER.As was shown with IRAC imaging at 3.6 and 4.5 µm (Arendt et al. 2020) and with mid-IR imaging (Matsuura et al. 2022) the brightening of the ER started in the east in the 2000's after which it faded again and the brightening started in the west by the late 2010's.Figure 1 confirms that the W is still brighter than other portions of the ring over the whole wavelength range of the MRS.A comparison of the spectra taken from the E and W segments of the ER (Figure 2) is shown in Figure 13.The W is approximately 50% brighter than the E part.Apart from the lower overall flux in the E, there are also some spectral changes that can be seen in the spectra normalised at 20 µm in the right panel of this figure.The eastern part has a relatively stronger long wavelength part and the 10 µm silicate band is weaker and its peak has slightly shifted to the red in comparison to the W.This spectral difference can be explained by a relatively higher dust temperature in the W. To model the dust parameters like the dust mass and temperature of the two segments requires a careful determination of the local underlying synchrotron continuum for the two ring segments and is planned for a future paper.

Other Dust Emission
In contrast to the inner ER, no dust emission was detected in the region corresponding to the two outer rings which are three times more extended than the ER and currently remain unaffected by the blast wave.
The ejecta dust was too faint to be detected by Spitzer.Dust emission may be present in the JWST/MRS ejecta spectra, however, this requires more mature calibration and extraction procedures to separate any potential contamination from the ring, as such the ejecta's solid-state component is subject to a future paper after follow-up Cycle 2 JWST/MRS observations have been obtained (PID #2763).

Emission from the ejecta and ejecta -ER interaction
During our search for small-scale emission line structures in the SN 1987A system we found that, in addition to the primary emission line peaks with radial velocities consistent with an origin in the ER, outer rings, or diffuse ISM, some lines had additional high-velocity blue-and/or red-shifted features in the northeast and southwest of the ER.Two of these, [Fe ii] 5.340169 µm and [Fe ii] 25.98839 µm have already been described in Sect.3.1.1.Further to these, high-velocity components were detected for [Ni ii] 6.636 µm, [Ar ii] 6.985 µm, and [Ne ii] 12.814 µm.Because the spectral and spatial information complements each other nicely, we discuss first the integrated line spectra and later the spatial maps.One may note that especially at the longer wavelengths the velocity resolution is considerably better than the spatial, although smoothed over a surface of constant line of sight velocity.From observations of the HST light curve and imaging, we know that the ejecta is powered by radioactive decay of 44 Ti and X-rays from the ejecta -ER interaction (Larsson et al. 2011).The X-ray-shielded inner ejecta is mainly powered by radioactive decay of 44 Ti, while the outer parts are heated by the photoelectric absorption of X-rays (Fransson et al. 2013).In the spectra, the ejecta are seen as broad lines expanding with velocity 2000 − 10, 000 km s −1 , where the dense, metal-rich core of the SN has a velocity of 2000 − 3000 km s −1 , mainly emitting in lines from Fe i-ii, O i, Mg i-ii, Si i and Ca ii, but also from H i, and He i, mixed into the core (Chugai et al. 1997;Kozma & Fransson 1998;Jerkstrand et al. 2011).The 44 Ti input is partly in gamma-rays and partly positrons, where the former mainly escapes at these late epochs while the positrons are deposited locally in the Fe and O-rich gas.Because of the very low temperature of the gas in the core, ≲ 200 K, it is expected that a major fraction of this energy will be emitted in the [Fe ii] 25.99 µm line (Kozma & Fransson 1998;Jerkstrand et al. 2011).This and other [Fe ii] lines are therefore of special interest for understanding the conditions in the core of the SN.
In Fig. 14 we show an extraction of the central part (see Fig. 2) of the ejecta in the [Fe ii] 25.99 µm line, together with the NIR [Fe i] 1.443 µm and [Fe ii] 1.534 µm lines from the JWST NIRSpec spectra (Larsson et al. 2023).The three lines have similar profiles, in spite of different ionization stages and excitation energies.The [Fe ii] 25.99 µm line, which has the best spectral resolution, extends from ∼ 3400 km s −1 on the blue side to ∼ 5300 km s −1 on the red side.The line is quite asymmetric, with a steep blue side and a more extended red side.
This central extraction showed only faint traces of the other lines mentioned above.However, an extraction also covering the ER (Fig. 2) had a clear signal for all of these.In Fig. 15 we show on the left panel a comparison of this extraction with the central extraction.In this plot we have scaled the total extraction by a factor of 0.14 to agree at ∼ 1500 km s −1 , mainly reflecting the small central extraction region compared to the total extraction (Fig. 2).The most notable difference is the strong peak at ∼ 2200 km s −1 , but also a fainter peak at ∼ 4100 km s −1 .The peak at ∼ 2200 km s −1 corresponds to the interaction region in the SW, as can be seen on the line maps below.The central peak, not seen in the central extraction, may come from the part of the dense ejecta above the ER plane, which has been bright also in earlier SINFONI observations (Larsson et al. 2016), when it was still located inside the ER.Because of the expansion, it is now projected against the ER.
In the middle panel, we show a comparison with the [Fe ii] 5.340 µm line.This comparison is interesting because it shows a major difference with the [Fe ii] 25.99 µm line in that the low-velocity emission is very low compared to the emission peaks at ∼ +2200 km s −1 and ∼ −1800 km s −1 .This indicates that there is little [Fe ii] 5.340 µm emission from the central parts of the ejecta, as was indicated from the central extraction.It also separates out the emission from the ER interaction on the blue side, which is not easily distinguished from the 25.99 µm line.
The right panel of Fig. 15 shows a corresponding comparison with the [Ne ii] 12.81 µm line.The [Ne ii] line exhibits a central component that originates from the ER and which is stronger and narrower than the central [Fe ii] 25.99-µm line emission.Outside of its central ER emission, the [Ne ii] profile exhibits a similar line shape as the [Fe ii] 5.340 µm line, with a stronger red component and a weaker blue component at similar velocities.
To determine the morphology and radial velocity of these structures spatially we again produced intensity and velocity maps using the 0th and 1st moments using Spectral-Cube, this time masking the contribution of the low-velocity component.The resulting maps are shown in Fig. 16.Significant, highly red-shifted emission centred at ∼ +2000 km s −1 and projected towards the southwestern ER is observed in all lines.In the cases of [Ar ii] 6.985 µm, [Ni ii] 6.636 µm, and [Fe ii] 25.988 µm, a blue-shifted component is also observed, though projected towards the northeast of the ER.The location and velocity of these features agree with the JWST NIRSpec study of Larsson et al. (2023) wards the northeast and southwest of the ER.Larsson et al. (2023) attributed these features to the interaction of the expanding ejecta with the reverse shock and also noted that northeast-southwest elongated morphology has been observed in all other atomic lines from the ejecta (Wang et al. 2002;Kjaer et al. 2010;Larsson et al. 2013Larsson et al. , 2016)).We suggest here that the high-velocity emission observed in our MIRI/MRS observations results from that same ejecta-reverse shock interaction.A more detailed analysis of the ejecta and interaction regions will be presented in future work.

DISCUSSION
The MIRI MRS observations presented above reveal the first spatially resolved spectral maps in the 5 to 28 µm wavelength region.The gas emission lines arise in the ER, the ejecta, the outer rings and the region between the ER and the outer rings.The continuum, which is dominated by dust emission, arises primarily in the ER.
The improved spectral resolution over Spitzer reveals the kinematic origin of the lines as well as enabling the discovery of new gas emission lines such as [Fe ii] 5.34 µm, [Na iii] 7.32 µm, [Ar iii] 8.99 µm, and [Ni ii] 10.68 µm (Sect.3.1).The origin of the narrow but centrally located [Ar ii] lines is discussed by Fransson et al. (2024).The [Ni ii] lines set an electron temperature of T e ≤ 6500 K in the western lobe.This relatively low temperature plus the low degree of ionization of the broad-line emitting species in the ER are consistent with dense post-shock cooling gas emission.The current higher density of the western side compared to the rest of the ER has likely also elevated the temperature and heating of the dust in the western lobe (Matsuura et al. 2022).
The origin of the narrow-line component, associated with the spectra of the more highly ionized species, is less clear.On the one hand, the surface brightnesses of the narrow lines peak at the equatorial ring, indicating some degree of association with the material in the ER.On the other hand, the narrow-line emission component shows no significant radial velocity variations and extends far beyond the ER, permeating the outer rings and reaching the edges of the MRS field with declining surface brightness.Candidate mechanisms to produce this narrow-line, higher degree of ionization, emission component include photoionization by X-ray and ultraviolet photons from shocks passing through the ER, and flash ionization by the original UV photon pulse of circumstellar or interstellar material of sufficiently low density that recombination has not yet completed.
Previously Dwek et al. (2010) and Arendt et al. (2016) used a two-component dust model to fit the lowresolution Spitzer/IRS spectra, consisting of a warm dust amorphous silicate component with a temperature of 180 K to model the 8-35 µm emission and a secondary featureless hot dust component fit to the 5-8 µm flux assumed to have an amorphous carbon (T = 460 K) composition.This secondary component may also be fit by iron grains as its optical constants are also featureless which would result in hot dust temperatures of ∼350 K. Arendt et al. (2016) determined temperatures of 190 K and 525 K for the silicate and carbonaceous components, respectively.Using the same sets of mass absorption coefficients as Dwek et al. (2010) and Arendt et al. (2016) we obtain slightly cooler temperatures for both components in our dust model fits to the MRS data, as we account for the contribution from the extrapolated long-wavelength synchrotron component and the H and He continuum emission in the near-IR.Our multi- In our analysis, astrodust (Hensley & Draine 2023) optical constants best fit the NIRSpec and MIRI spectra.Astrodust is comprised of C, O, Mg, Si, and Fe in a single-grain material, including amorphous silicates in the form Mg 1.3 (Fe,Ni) 0.3 SiO 3.6 and no large carbonaceous grains.This eliminates the need for two-grain compositions (silicates plus a featureless component typically assumed to be carbon) spatially coexisting, which form from completely divergent condensation pathways in the RGB phase of evolution.Grains formed during this phase are primarily oxygen-rich, due to CNO processes in the core resulting in low carbon abundances, with any remaining C atoms locked up in CO molecules.The additional benefit of our model is that the cooler temperature we obtain for our hot dust component can be reached by collisional heating under thermal equilibrium, eliminating the need for the more problematic stochastic heating which very small amorphous hydrocarbons from previous fits would require (Arendt et al. 2016;Dwek et al. 2010).
While beyond the scope of the present study, detailed modelling of the various components of the SN 1987A MRS spectra with laboratory data is an important step toward understanding the evolution of dust in SNe.Using Dorschner et al. (1995) optical constants for amorphous silicates which Jones et al. (2014) found best fit the spectra of RSGs in the LMC results in fits to the ER which are too sharp to reproduce the 10 and 20-micron features.Likely, some additional components at the few percent level including SiO 2 , Al 2 O 3 , CaCO 3 , Fe oxides, carbides, and sulphur compounds, which are all included in the artificial astrodust composite grains, would be required.Pure Fe is ruled out as a viable candidate due to the large dust mass that would be required (Table 4).
The evolution of the hot dust <15 µm component is consistent with the steady decline seen in the 3.6 and 4.5 µm flux measured during Spitzer's warm mission.The Spitzer cryogenic data pre-day 8000 was found to follow the same trends as the X-ray emission prior to diverging post-day 8000 (post-year 21.9) (Arendt et al. 2016(Arendt et al. , 2020)).This corresponds to the collisional heating of the dust in the ER before a transition to radiative heating or a change in the physical conditions in the ring, for instance, grain destruction or a change in the dust grain size distribution.Thermal sputtering timescale suggests that grain destruction may become prominent from day 9200 (year 25.2) after the explosion (Dwek et al. 2010).This dust processing and destruction are conclusively shown for the hot dust in the MRS spectra obtained on day 12927 (year 35.4) , however, the effects are negligible for the cold dust emission which may even be growing in mass (though not significantly) due to the expanding propagation of the forward shock wave sweeping up remnant circumstellar material.
Dust grains in the ring are expected to have a distribution of grain sizes that are heated to different temperatures, with the resulting emission spectrum reflecting a continuous temperature gradient as a function of grain size.In our fits, this gradient is approximated by a twotemperature-component model.The temperature differential between our hot and cold components in the astrodust model is, therefore, likely a consequence of the dust grain size distribution, with the hot dust component corresponding to smaller grains and the cold component to the larger ones.The trend seen in Figure 12 can thus be interpreted as circumstellar material from the progenitor being swept up by the blast wave expanding through the ring.As the shock propagates through the ring, the mass of the shock-heated dust grows in proportion with the total swept-up gas mass, as seen on days 7000-8000 (Dwek et al. 2010).As they are replenished, the smaller grains are being simultaneously destroyed in the ER on timescales of a few years (Arendt et al. 2016).Once the grains stop being replenished in sufficient amounts, the destruction of the small grains dominates, whereas the large grains are more likely to survive the passage of the high-velocity shock, resulting in the divergence in dust masses between the hot and cold components seen in the MRS data.This scenario requires the hot grains to be non-carbonaceous, as the binding energy makes carbonaceous grains relatively stable against thermal and non-thermal sputtering (Fischera et al. 2002).
SOFIA observations (Matsuura et al. 2019) identify a 31.5-70µm excess which may be attributed to dust reformation after the passage of the forward shock, warm dust emission from the ejecta or the survival of large grains.Matsuura et al. (2019) prefer the reformation of dust grains after the shock front.There is however no indication in the MRS data of a warm dust component from the ejecta nor a re-formed dust component or substantial grain growth in the ring, given the consistency between the long-wavelength MIRI/MRS and Spitzer/IRS data.This suggests that dust reformation in the post-shock region of this system does not easily occur or occurs on much longer time scales given the low densities (Biscaro & Cherchneff 2014).Our results support the survival of large grains in the shocks.
The different line shapes between the [Fe ii] 25.99 µm line on the one hand and the [Fe ii] 5.340 µm and [Ne ii] 12.81 µm lines, on the other hand, show that the excitation of the outer parts of the ejecta, close to the ER, and the central region are different.The excitation energy of the upper level of the [Fe ii] 25.99 µm line corresponds to 550 K, while it is 2676 K for the [Fe ii] 5.340 µm line, 1115 K for the [Ne ii] 12.81 µm line and 2045 K for the [Ar ii] 6.985 µm line.From the modelling of the 8-year spectrum in Jerkstrand et al. (2011) the temperature in the 44 Ti-powered inner core is 70 − 170 K.The excitation of the [Fe ii] 25.99 µm line can therefore occur at low temperatures, while the much higher excitation energies of the other lines require a higher electron temperature.This is consistent with the observed dip in the central line profile of the [Fe ii] 5.340 µm line.The temperature at 35 years may be different from the above, but these line profiles indicate a very low temperature in the core also at the present epoch.Models of the X-ray transfer in the ejecta from the ejecta -ER interaction also show that most of the X-rays have not yet penetrated to the central parts of the core (Fransson et al. 2013;Kangas et al. 2022b).The high-velocity Fe-rich regions which have expanded to close to the reverse shock are, however, not shielded from the X-rays due to a much smaller column density between the reverse shock and ejecta.They will therefore be heated efficiently, as is directly seen from the line profiles and the maps, and will have a much higher temperature.More detailed modelling of the shock interaction is required for a more quantitative estimate.

CONCLUSIONS
We have presented mid-IR spectra of SN 1987A obtained with the MRS instrument onboard JWST on day 12,927 post-explosion.The 5-28 µm MRS spectroscopic observations spatially resolve the ER, ejecta, outer rings and surrounding environment.We find: 1. Infrared line emission from the ER coincides with the dense clumps of gas observed as hotspots in the optical images, while we find that the dust continuum emission is more spatially extended.4. The velocity structures of the ER lines indicate that they arise from two regions: The singly ionized species show broad lines (280-380 km s −1 FWHM) and their radial velocities indicate that they originate from the expanding ER.Their low degree of ionization, together with the upper limit of 6500 K on the ER electron temperature provided by the ratio of the [Ni ii] 6.64-and 10.68-µm lines, are consistent with recombination of postshock gas.Lines from the more highly ionized species are narrower (99-171 km s −1 ) and their lack of any radial velocity variations across the ER and beyond indicates that they arise from a different emission component.The ratio of the [Ne v] 14.32-and 24.32-µm lines at the ER implies electron densities for this component of between 700-4300 cm −3 for T e = 5000−25000 K, while a [Ne v] 14.32-µm continuum-subtracted image shows the emission to also be associated with material located in the outer rings and beyond.The higher degree of ionization of these species may have been produced by either X-ray and ultraviolet photons from shocks progressing through the ER, or else by flash ionization by the UV photon pulse associated with the original supernova explosion.

Emission
5. The overall shape and positions of the features in the ER spectrum are generally consistent with the Spitzer/IRS spectra obtained between days 6000 and 8000.No resonances from other solidstate features are apparent in the higher-resolution MRS spectra.However, the relative flux densities < 15 µm are significantly lower in the MRS spectra.
6.The 0.9-28 µm NIRSpec and MIRI spectral range is best-reproduced by a single set of astrodust optical constants emitting at 157±4 and 334±12 K.This is contrary to previous models of dust emission from the ER which require a secondary pop-ulation of very small, carbon or iron grains of unknown origin radiating at significantly higher temperatures to spatially coexist with the cool astronomical silicate grains.The inclusion of the continuum from synchrotron and bound-free emission in our models significantly lowers the temperature of the 'hot' component even when this secondary population is represented by amorphous carbon or iron grains.
7. Fitting both the MRS and Spitzer spectra with the same two-temperature model produces a consistent dust temperature in the ring between day ∼ 7000 and 12927 (year 19.2 and 35.4).However, there is a significant difference in the mass of the hot dust at day 12927 compared to 7955, suggesting the hot, smaller dust grains are being destroyed.
8. The large dust grains from the progenitor preferentially survive the SN explosion and processing by the forward shock as it evolves into an SNR.
9. There is a shift in the amorphous silicate feature position between the East and West ER components indicating a difference in the shock processing.
10. Our spectra show a number of lines from the ejecta which, with the exception of the [Fe ii] 25.99-µm line, only have been seen during the first years after explosion.This includes [Fe ii] 5.34 µm, [Ar ii] 6.99 µm and [Ne ii] 12.81 µm.The [Fe ii] 25.99 µm line is in the centre still powered by the radioactive 44 Ti decay, while the emission at ∼ 2200 km s −1 is powered by the ER -ejecta interaction.This is also consistent with the line profiles and spatial distribution of the other lines above.The central [Fe ii] 5.34 µm to 25.99 µm line ratio is consistent with a very cool central ejecta, shielded from the X-rays.

2.
OBSERVATIONS AND DATA REDUCTION SN 1987A was observed with the MIRI/MRS on 2022-07-16 as part of the guaranteed time programme #1232 (PI: Wright).The observations consist of 94 groups of three integrations using the fastr1 readout pattern for all three bands (SHORT, MEDIUM, LONG), which

Figure 1 .
Figure 1.Sample slices from our 12 MRS sub-band cubes, with the band labels shown in the top-left.The white lines in each pane represent 1.0 arcsec to highlight the increasing size of the FOV, as well as the decreasing spatial resolution from channels 1 to 4. 'Star 3' is visible to the lower left of the ER in bands 1A-2A.North is up, east is left.

Figure 2 .
Figure 2. Top row: The spectra extraction regions for the ER spectra plotted on band A in each channel.The total ER spectrum was extracted from the white elliptical annulus, with the east and west ER extracted from the segments on the left and right.Bottom row: Same as top but for the cardinal point and ejecta spectra.In both rows, the location of the outer rings is shown by the two faded ellipses, the channels 1-3 background inner boundaries are shown by the dashed cyan circles, and the channel 4 background region is shown by the dashed cyan polygon.The 'entire' SN 1987A system spectrum was extracted from a region slightly smaller than the inner boundary of the background regions and excluding Star 3. North is up, east is left.
3.1) and the central compact ejecta is seen in the [Ar ii] and [Fe ii] emission lines (see Sect. 3.4). 2

Figure 3 .
Figure 3. Illustration of the different spatial distributions of the line emission and continuum.The left panel shows an integrated image of the H I 7.460 µm line after subtraction of the continuum, the middle panel shows an image of the continuum on each side of the line, and the right panel show the two images together.The black contours in the left and middle panels show the hotspots in a HST/F657N image obtained on 2022-09-06 (dominated by Hα emission, Rosu et al. in preparation).Note the good agreement between the H I 7.460 µm line and the HST image, while the continuum is more spatially extended.

Figure 4 .
Figure 4.The total MIRI MRS background subtracted spectrum of SN 1987A 12,927 days (35.4 years) after the explosion.The mid-IR spectrum at this epoch is dominated by dust emission from amorphous silicates in the inner ER.The line emission is from a combination of the ionised ER and the central ejecta (with a broad emission profile, e.g., [Fe ii] at 26 µm).The most significant lines are highlighted.
those listed by Arendt et al. (2016) from their Spitzer 10-40 µm spectral sequence shows the latter fluxes to be generally an order of magnitude larger.This might be due to the Spitzer line flux measurements having been made with larger apertures.A comparison of the relative line intensity patterns seen in the JWST-MRS and Spitzer-IRS datasets indicates that, relative to [S iv] 10.51 µm, the low degree of ion-ization [Fe ii] and [Ne ii] lines have brightened by factors of 5-12 in the MRS spectra.Arendt et al. (2016) suspected the presence of [Ni ii] 6.64 µm, H i 7.46 µm and [Ne vi] 7.65 µm in their R=90 5-10 µm Spitzer spectra -these lines are confirmed by the MRS spectra.The MRS detection of [Ar iii] 8.99 µm appears to be new.The other new line detections are [Fe ii] 5.340169 µm, [Ni ii] 10.6822 µm and [Na iii] 7.31775 µm, both seen in the whole-Ring and in the Ring-West spectra.

aFigure 5 .
Figure 5. Continuum subtracted emission line images for those species with emission from outside the ER.The grey circles overlaid on the images are the locations of the outer rings determined from the [Ne ii] 12.8135 µm morphology, with the lightgrey contours showing the ER continuum morphology determined from an 'off' position adjacent to the line.The colour bars are in units of MJy/sr.

Fig. 5
shows continuum-subtracted images for nine emission lines, using display ranges designed to bring out fainter levels of emission.All of the images show strong ER emission, peaking at the western quadrant.The outer rings are prominent in the [Ne ii], [Ne iii], [Ne v], [Ne vi] and [O iv] images, with the observed distributions of the latter three ions predominantly following the morphology of the outer rings.All of the line images shown in Fig. 5 display some level of extended emission beyond the extent of the inner and outer rings.For [Ne v] and [O iv] this lower surface brightness emission appears to closely envelop the outer rings but for [Ne ii], [Ne iii] and, in particular, [S iii] 18.71 µm, part of a large ring-like structure fills the north-eastern half of the IFS FoV, passing through the northern part of the ER.

Figure 6 .
Figure 6.Moment maps for a selection of species prominent in the ring ordered row-by-row.The left panes show the integrated line intensity.The middle panes show the radial velocity of the line emission, corrected for the systemic barycentric velocity of SN 1987A of +287 km s −1 (Gröningsson et al. 2008).The right panes show the line FWHMs.The grey circles overlaid on the images are the locations of the outer rings derived determined from the [Ne II] 12.8135 µm morphology, with the grey contours showing the ER continuum morphology determined from an 'off' position adjacent to the line.

Figure 7 .Figure 8 .
Figure 7. Same as Fig.6 but for a selection of H I lines.

Figure 9 .Figure 10 .
Figure9.Fits to the MRS and NIRspec data using the two-temperature dust compositions of astrodust(Hensley & Draine 2023) in the top panel, carbon and silicates(Draine & Lee 1984) in the middle panel, and silicates(Draine & Lee 1984) and iron dust in the bottom panel.The synchrotron emission (orange line) has been extrapolated from ALMA observations based on the models described inCendes et al. (2018) andCigan et al. (2019) (see text for details).The emission from the H and He continuum and lines (purple line) is described inLarsson et al. (2023).The best-fit parameters for the dust components are summarized in Table4.

Figure 15 .
Figure 15.Left: Comparison of the [Fe ii] 25.99 µm line profiles from the inner ejecta and the total line profile including the ring.The spike at -1150 km s −1 is the narrow [O iv] 25.91 µm line from the ER.Middle: Comparison of the line profiles of the [Fe ii] 5.340 µm and [Fe ii] 25.99 µm lines.Note the lack of emission between ±1500 km s −1 in the 5.34 µm profile.Right: The same for the [Ne ii] 12.81 µm line.

Figure 16 .
Figure 16.Same as Fig.6 but for high-velocity line components and without the FWHM maps..

Table 2 .
Detected emission lines, their measured heliocentric radial velocities, full-width half maxima, and line fluxes, for the whole ER extraction (subtending a solid angle of 1.072200 × 10 −10 sr).
a Vacuum wavelengths from the compilation by van Hoof (2018) b H i 8-6 7.502493 µm was also detected

Table 4 .
Dust Model Fit Parameters

Table 5 .
Temporal Evolution in Dust Properties line radial velocities, FWHM values and integrated line fluxes were measured at four cardinal points in the ER and for the whole ring.Nineteen emission lines, from thirteen species, were detected.The majority of the lines identified have previously been detected in Spitzer spectra of SN 1987A, including [Ni ii] 6.64 µm, H i 7.46 µm and [Ne vi] 7.65 µm which were previously only tentatively identified.However, [Fe ii] 5.34 µm, [Na iii] 7.32 µm, [Ar iii] 8.99 µm and [Ni ii] 10.68 µm are new detections.3. The MRS lines are approximately an order of magnitude fainter compared to the Spitzer detections on day 7954 (year 21.8), which could be partly attributable to a smaller extraction region for the MRS lines.However, relative to [S iv] 10.51 µm, the low degree of ionization [Fe ii] and [Ne ii] lines have brightened by factors of 5-12 in the MRS spectra.