Ground-based and JWST Observations of SN 2022pul. II. Evidence from Nebular Spectroscopy for a Violent Merger in a Peculiar Type Ia Supernova

We present an analysis of ground-based and JWST observations of SN 2022pul, a peculiar “03fg-like” (or “super-Chandrasekhar”) Type Ia supernova (SN Ia), in the nebular phase at 338 days postexplosion. Our combined spectrum continuously covers 0.4–14 μm and includes the first mid-infrared spectrum of a 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. A strong, broad, centrally peaked [Ne ii] line at 12.81 μm was previously predicted as a hallmark of “violent merger” SN Ia models, where dynamical interaction between two sub-M Ch white dwarfs (WDs) causes disruption of the lower-mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to reproduce the optical iron emission better, and add mass in the innermost region (<2000 km s−1) to account for the observed narrow [O i] λ λ6300, 6364 emission. A violent WD–WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SNe Ia.

1. INTRODUCTION Supernovae (SNe) imprint important information about the composition, geometry, kinematics, and ionization of the ejected material in their late-time (nebular) emission-line spectra.In the nebular phase (≳ 100 days past maximum light; Bowers et al. 1997;Branch et al. 2008;Silverman et al. 2013;Friesen et al. 2014;Black et al. 2016), the ejecta have expanded and the opacity has dropped enough to see directly through even the innermost regions and allow forbidden lines to form.In white-dwarf (WD) supernovae, decay of radioactive 56 Co (produced by the radioactive decay of 56 Ni synthesized in the explosion) provides most of the energy that excites the low-density nebular ejecta, producing forbidden emission lines throughout the ejecta (Axelrod 1980;Nadyozhin 1994;Jerkstrand 2017).These glowing ashes contain clues for determining the SN origin.Nebular-phase spectroscopy is therefore an important forensic tool for investigating the remains of exploded WDs, whose explosion channels are not well understood.
The majority of WD SNe are classified as Type Ia supernovae (SNe Ia); however, a growing number and variety of extreme or peculiar SN Ia subtypes have been discovered (for recent reviews, see Taubenberger 2017;Jha et al. 2019;Liu et al. 2023).These peculiar SNe result from the thermonuclear explosion of a WD, but they generally do not follow the Phillips (1993) relationship between luminosity and lightcurve shape that makes normal SNe Ia so useful for cosmological distance measurements (Phillips 1993;Benetti et al. 2005;Taubenberger 2017;Jha et al. 2019).No consensus has been reached concerning the possible and dominant pathways connecting progenitor systems, explosion mechanisms, and observed WD SN properties.
The origins of normal SNe Ia are particularly elusive in part because variations between objects are relatively small and many of the models predict similar observational properties.By contrast, extreme and peculiar WD SNe exhibit distinguishing features that can make identification of their origins clearer (Foley et al. 2013;Kromer et al. 2013;Hsiao et al. 2020).Understanding the scenarios that produce peculiar WD SNe helps us understand the unusual and extreme ways WDs can explode while simultaneously narrowing the model space for normal SNe Ia.
A rare type of peculiar WD SN characterized by high luminosities, broad light curves, C II absorption at early times, relatively low ejecta velocities, and low ionization state nebular-phase spectra are "03fg-like" (commonly called "super-Chandrasekhar") SNe Ia (Howell et al. 2006;Hicken et al. 2007;Scalzo et al. 2010;Silverman et al. 2011;Taubenberger et al. 2013;Hsiao et al. 2020;Ashall et al. 2021;Srivastav et al. 2023).Explaining the light curve of SN 2003fg, the class prototype, by decay of 56 Ni alone requires a total mass above 1.4 M ⊙ , the Chandrasekhar mass limit (M Ch ) for a nonrotating WD (Howell et al. 2006).A mass above this limit can be achieved either by additional support from differential rotation in a single WD (Yoon & Langer 2005) or a merger of two WDs with total mass greater than M Ch (Howell et al. 2006;Hicken et al. 2007).Alternative explanations that have been invoked to explain this diverse class of bright SNe Ia include the "core-degenerate" scenario, where a near-M Ch C/O WD explodes within the C-rich envelope of an asymptotic giant branch (AGB) star (Khokhlov et al. 1993;Hoeflich et al. 1996;Livio & Riess 2003;Quimby et al. 2007;Kashi & Soker 2011;Hsiao et al. 2020;Lu et al. 2021), and the violent merger of two sub-M Ch WDs (Dimitriadis et al. 2022;Srivastav et al. 2023;Siebert et al. 2023).These potential avenues account for the excess luminosity and relatively low velocity of these objects by interaction with the AGB star envelope or C/O-rich circumstellar material (CSM) from the disrupted smaller WD, respectively.For more comprehensive reviews of 03fg-like SNe Ia, see Taubenberger (2017), Ashall et al. (2021), and Paper I (Siebert et al. 2024).
In contrast to the classical WD merger scenario where the merger produces a single-object remnant that later explodes, the violent merger model proceeds via dynamical interaction between two sub-M Ch WDs that causes disruption of the secondary (lower mass) WD and subsequent detonation of the primary (higher mass) WD (for a review, see Pakmor 2017).The primary WD is burned first and its ejecta expand beyond the secondary by the time that the secondary is burned ∼ 1 s later and the ejecta of the secondary expand into the center of the ejecta of the primary.Thus, the intermediatemass elements (IMEs) produced by burning of the secondary will be centrally located within the ejecta.Depending on the mass ratio, mass of the primary WD, density of the secondary WD at the time of detonation, and composition of both WDs, a fairly wide range of SN properties can be produced (Pakmor et al. 2011(Pakmor et al. , 2012;;Kromer et al. 2013;Pakmor 2017).In general, violent mergers are asymmetric scenarios which can result in enhanced production of IMEs from burning of the secondary, lower levels of iron-group elements (IGEs) from burning of the sub-M Ch primary, centrally located C/O/Ne/Mg and IMEs, and a low ionization state at late times due to high central ejecta density.
Nebular spectra can reveal all of these properties.This information can be difficult to disentangle, however, owing to significant line crowding and blending at optical wavelengths.Spectroscopy with JWST in the near-infrared (NIR, which extends beyond the NIR wavelength range accessible from the ground) and mid-infrared (MIR) grants access to emission lines from several species of IMEs, and the lines are comparatively isolated so individual line profiles can be distinguished.These NIR and MIR spectral features help clarify the information contained in the optical and groundbased NIR spectra, aiding in discriminating between progenitor systems and explosion mechanisms.Observations of SN 2021aefx, the first SN Ia with JWST NIR+MIR nebular spectra, revealed strong, stable Ni indicative of a high-mass WD progenitor and evidence for stratified ejecta resulting from some type of detonation (Kwok et al. 2023).Further modeling found SN 2021aefx to be consistent with an offcenter delayed detonation of a near-M Ch mass WD (DerKacy et al. 2023).
We obtained a JWST NIR+MIR nebular spectrum of SN 2022pul, a peculiar 03fg-like SN Ia (JWST Cycle 1 GO 2072; PI S. W. Jha).Our JWST spectrum of SN 2022pul is the first MIR spectrum of an 03fg-like SN Ia, and it exhibits a strong thermal (T ≈ 500 K) dust continuum (analyzed in Pa- per III; Johansson et al. 2024).It also displays several differences from the optical-MIR spectrum of SN 2021aefx.Additional details about SN 2022pul and its unique photometric and spectroscopic evolution can be found in Paper I (Siebert et al. 2024).
Here, we present analysis of the dust-continuumsubtracted optical through MIR nebular-phase spectrum of SN 2022pul at 338 d post explosion (this assumes that the dust is external, see Johansson et al. 2024 for a discussion).In Section 2 we identify NIR and MIR emission lines and highlight distinctive spectral properties in comparison to SN 2021aefx.Section 3 presents emission-line-profile fits to the dominant spectral features and their implications on the distribution of material in the ejecta.In Section 4 we compare SN 2022pul to nebular spectroscopic predictions (Blondin et al. 2023) of the violent merger model from Pakmor et al. (2012).We discuss the implications of our results and conclude in Section 5 that SN 2022pul was most likely produced by the violent merger of two WDs.

SPECTRAL ANALYSIS
The continuous optical+NIR+MIR spectrum of SN 2022pul in the nebular phase at 338 rest-frame days post-explosion (estimated explosion date MJD 59785.3)exhibits many unique properties, including a clear NIR+MIR dust continuum that closely resembles a blackbody (see also Siebert et al. 2024).Here, we focus on analysis of the spectral lines, while a detailed analysis of the dust properties is given by Johansson et al. (2024).We remove the dust contribution to the spectrum by subtracting off a blackbody of T = 500 K from the 2.5-14 µm region, and use this dust-subtracted spectrum throughout this work.We do not observe other spectral signs of dust such as CO or SiO emission (although line overlap may complicate clear identification of dust ex-tinction in the optical).For details of the observations and data reduction of the spectra used in this work, see Siebert et al. (2024).
Shown in Figure 1, several highlights of SN 2022pul's combined nebular spectrum include a low mean ionization state, strong emission from IMEs, relatively weak emission from IGEs, highly asymmetric line profiles, narrow [O I] λλ6300, 6364 and remarkably strong [Ca II] λλ7291, 7324 in the optical spectrum, detection of [S IV] at 10.51 µm, and the first unambiguous detection of [Ne II] at 12.81 µm in an SN Ia.

Spectroscopic Comparison with SN 2021aefx
SN 2021aefx, a spectroscopically and photometrically regular SN Ia near maximum light, is currently the only other SN Ia with a published continuous spectrum from 0.3 µm to 14 µm in the nebular phase (270 d post-explosion; Kwok et al. 2023), making it a natural comparison object for the peculiar "super-Chandra" 03fg-like SN 2022pul.After subtracting the dust continuum, SN 2022pul has a greater fraction of total flux in the NIR and MIR, and a smaller fraction of optical flux than SN 2021aefx.The light curve of SN 2022pul shows that at late times, the I-band flux exceeds the flux in other bands (probably due to the strong [Ca II] λλ7291, 7324 doublet) and declines more slowly than typical SNe Ia (Siebert et al. 2024).This is consistent with findings from Ashall et al. (2021) that 03fg-like SNe Ia are brighter and decline more slowly in the NIR than normal SNe Ia.
This redder color may be due in part to the comparatively low overall ionization state of SN 2022pul, where additional flux from lines of lower ionization falls in the infrared.The flux from optical [Fe III] around 5000 Å is greatly diminished relative to most SNe Ia (Siebert et al. 2024 ilarly, the JWST NIRSpec spectrum reveals that the NIR [Fe III] lines around 2.2 and 2.9 µm are weak while the NIR [Fe II] lines around 1.3 and 1.6 µm are strong compared with SN 2021aefx.In the MIR, we see the same ionization trend in the Ni, Ar, and Co where the ionization is lower in SN 2022pul than in SN 2021aefx.In agreement with the ionization inferred from the MIR, the strength of the [Ca II] λλ7291, 7324 doublet in SN 2022pul is remarkable (Siebert et al. 2024), whereas in SN 2021aefx the spectral complex near 7300 Å was attributed to [Fe II] and [Ni II] with little to no contribution from [Ca II] (Kwok et al. 2023) -yet they have similar strengths of [Ca IV] near 3.2 µm.This shift to a lower mean ionization for SN 2022pul suggests a larger density or mass of the ejecta, but the presence of Ca IV also indicates that some higher ionization bubbles are present in both, discussed further in Section 4.
The NIR does not display lines other than those present in SN 2021aefx, though the [Fe III] is weaker and the three [Fe II] features in the range 4-5 µm are stronger.Blondin et al. (2023) showed that the broad feature near 3.2 µm is dominated by [Ca IV] 3.21 µm, with weaker contributions from [Fe III] 3.23 µm and [Ni I] 3.12 µm.The NIRSpec observation of SN 2022pul further confirms this broad line around 3.2 µm as [Ca IV] 3.21 µm because the [Fe III] and Ni emission are weaker in SN 2022pul than in SN 2021aefx, but the 3.2 µm feature has a similar strength.Several unique features in the MIRI/LRS spectrum of SN 2022pul are shown in Figure 2 and described below.

[Ne II] emission
In addition to all of the NIR and MIR lines present in SN 2021aefx (line identifications are given by Kwok et al. 2023), SN 2022pul exhibits a strong, broad, centrally peaked feature around 12.8 µm that we identify as [Ne II] 12.81 µm (see Figure 2).Predicted by Blondin et al. (2023), the presence of strong, centrally peaked [Ne II] 12.81 µm is a distinguishing feature of a violent WD-WD merger not present in any of the other models considered in their study.Ne is produced centrally in a violent WD-WD merger because it is a product of burning the less dense secondary WD, whose ejecta expands into the center of the primary WD's ejecta after the primary has burned and expanded.The observed [Ne II] 12.81 µm line is broad, asymmetric, and dissimilar to other line-profile shapes seen in SN 2022pul and other SNe Ia with MIR spectra.

Strong Intermediate-Mass Elements
One of the major differences between the MIR spectra of SN 2022pul and SN 2021aefx is the strength of the IME (hereafter defined as S, Ar, and Ca) emission lines.Whereas in SN 2021aefx, [Ar II] 6.98 µm was blended out by the neighboring Ni emission lines such that its shape was uncertain, [Ar II] 6.98 µm is the dominant emission line in the MIR spectrum of SN 2022pul.Its peak flux is blueshifted by ∼ 2000 km s −1 , and is over twice as strong in F ν as [Ne II] 12.81 µm, the next strongest MIR emission line.Similarly, [Ca II] λλ7291, 7324 dominates the optical spectrum.
The [Ar III] 8.99 µm emission line is also quite strong in SN 2022pul: roughly three times stronger in F ν than in SN 2021aefx at a very similar time post-explosion (see .The [Ar III] 8.99 µm line continues to exhibit a boxy, flat-topped shape, with a similar width, but the slanted top is much steeper.This steep slant can be produced by an asymmetric shell of [Ar III] emission, discussed further in Section 3. The feature at 10.5 µm in SN 2022pul results from a combination of [S IV] 10.51 µm and [Co II] 10.52 µm emission.The strength and flat shape of the feature at 10.5 µm cannot be explained by [Co II] 10.52 µm alone, as is further discussed in Section 3, constituting a firm detection of [S IV].Kwok et al. (2023) did not invoke a contribution from [S IV] 10.51 µm in SN 2021aefx; however, placed in context with SN 2022pul, the "shoulder" feature starting at 10 µm in SN 2021aefx might actually be [Si IV] with a similar boxy shape as the [Ar III] but with about half the strength.Reevaluation of [S IV] in SN 2021aefx will be included in future work.
All of the IME emission lines in SN 2022pul are strong relative to SN 2021aefx, and their peak emission is asymmetric and blueshifted.

Weak Iron-group Elements
In contrast to the strength of the IMEs, the IGEs (hereafter Fe, Co, and Ni) appear somewhat weaker in the MIR than in SN 2021aefx.Part or all of this apparent decrease in IGEs may be due to the shift to lower ionization.For example, the [Ni II] 6.64 µm line is much stronger in SN 2022pul, but the [Ni III] lines at 7.35 and 11.00 µm are several times weaker, and the [Ni IV] 8.41 µm line is so weak that it is difficult to see without an arcsinh or log scaling of the flux.This points toward some of the flux in the higher ionization IGE emission lines being shifted to the singly ionized emission lines.
The isolated [Co III] 11.89 µm line shows a distinctly asymmetric shape whose peak is redshifted by ∼ 3500 km s −1 .Compared with SN 2021aefx, the [Co III] 11.89 µm line is somewhat narrower.As will be shown and discussed in Section 3, the asymmetric redshifted shape of the [Co III] 11.89 µm line profile is consistent with all of the other IGE ions throughout the NIR and MIR, as well as the redder side of the optical.
3. EMISSION-LINE PROFILES Nebular emission-line profiles are important indicators of the geometry of the SN ejecta because they imprint the flux at all projected line-of-sight velocities (and therefore projected radius) through the ejecta.The shape of the lines arises from a combination of the density and excitation of the elements within the ejecta (for a review, see Jerkstrand 2017).In SN 2022pul, we find evidence in the highly asymmetric line profiles for distinct distributions of the IGEs (Fe, Co, Ni), IMEs (S, Ar, Ca), O, and Ne within the ejecta.Unfortunately, it is not possible to determine the full threedimensional (3D) geometry that gives rise to each profile owing to projection effects and the possibility of bubbles, mushroom shapes, clumps, or other such non-spherically symmetric distributions.However, by fitting line profiles throughout the optical, NIR, and MIR, we can determine the composi- Notably, the Ar and Ca (IMEs) have a blueshifted peak with a more extended red tail, while the Co and Fe (IGEs) display an opposite trend, with a redshifted peak and a more extended blue tail.
tion within slabs perpendicular to the line of sight, which resolves some of the geometrical layout.In this section we fit line profiles consistently across the optical+NIR+MIR and discuss their implications.
For our line-profile fitting, we need relatively high resolution to discern line-profile shapes, so we combine an optical MMT/Binospec spectrum at 332 d post-explosion, a NIR Keck/NIRES spectrum at 316 d post-explosion (the nearest-phase ground-based NIR spectrum available), and the JWST/NIRSpec+MIRI spectrum beyond 1.83 µm at 338 d post-explosion (Siebert et al. 2024).Despite the difference in phase, we use the Keck/NIRES spectrum (rebinned to increase the signal-to-noise ratio per pixel) instead of the JWST/NIRSpec spectrum in the range 1.0-1.83µm owing to the much higher resolution.In addition to removing the dust contribution, we also remove all remaining underlying MIR continuum in the JWST data (see Johansson et al. 2024, for discussion of the MIR continuum).
The optical, NIR, and MIR spectra of SN 2022pul are all essential to this work: the higher resolution of the groundbased optical and NIR spectra provides clarity for several key line-profile shapes, and the lower resolution JWST data contain isolated lines and emission from IMEs.These properties of the data are complementary for interpreting the line profiles of SN 2022pul.

Line-Fitting Procedure
Following the general approach of nebular line fitting from Maguire et al. (2018), Flörs et al. (2018Flörs et al. ( , 2020)), and Kwok et al. (2023), we model the observed spectrum as a superposition of line profiles for all reasonably strong contributing lines.Restrictions on the fitting are enforced such that all lines of the same ion have the same profile shape parameters and kinematic offset from the central wavelength.We note similarities between the nebular spectrum of SN 2022pul and the nebular spectrum of the violent WD-WD merger model of Pakmor et al. (2012) computed by Blondin et al. (2023) (their MERGER model).Thus, we compute this model out to the same phase as SN 2022pul (338 d post-explosion) and we fix the relative line strengths of each contributing line in our fits to the relative line strengths in the model.This allows us to base the relative line strengths on temperatures and densities similar to those in SN 2022pul.Details of and comparisons to the merger model are presented in Section 4.
While the relative line strengths within a given ion are fixed by the MERGER model at 338 d, we allow the relative amplitudes between each ion to vary as fitting parameters.We fit the optical through MIR and include the following ions: Fe I, Fe II, Fe III, Co II, Co III, Ni II, Ni III, Ni IV, S IV, Ca II, Ca IV, Ar II, Ar III, O I, and Ne II.Additionally, to account for asymmetry in SN 2022pul, we fit each ion with non-Gaussian line profiles, further described below.
The resolution of the NIRSpec/Prism and MIRI/LRS data introduces significant uncertainties into the fits for parameters that measure width.For features beyond 1.8 µm, these resolution uncertainties are on the order of 1200 km s −1 near 3 µm, 2000 km s −1 near 6 µm, and 500 km s −1 near 12 µm.Uncertainties from the fitting were calculated using a bootstrap resampling method.Uncertainties in width parameters such as the full width at half-maximum intensity (FWHM) and inner-shell velocity were obtained by adding the resolution and fitting uncertainties in quadrature.

Iron-Group Elements: Fe, Co, Ni
Starting with the isolated and well-resolved [Co III] 11.89 µm line, we approximate the shape of the line profile by the sum of two Gaussians with FWHM = 5400 ± 1900 km s −1 and kinematic offset v offset = 4600±2300 km s −1 , and FWHM = 8000 ± 2200 km s −1 and kinematic offset v offset = 1400 ± 1800 km s −1 .As shown in Figure 3, the [Fe II] 1.26 µm line has a shape very similar (taking into account the neighboring [Fe II] lines that create its red tail) to that of [Co III] 11.89 µm, so we fit all the [Fe II] lines across 1.20-1.83µm with a sum of two Gaussians and obtain a similar, but slightly more peaked, profile for [Fe II] with FWHM = 4700 ± 1200 km s −1 and kinematic offset v offset = 3800±500 km s −1 , and FWHM = 7400±2100 km s −1 and kinematic offset v offset = −600 ± 1300 km s −1 (see Figure 4 for a visualization of a double Gaussian profile).
A possible explanation for why Fe and Co do not exhibit exactly the same profiles is that while Co mostly traces the original 56 Ni profile, Fe is more widely distributed because it is partially primordial.In the regions dominated by the IMEs and O+Ne, Fe may be a significant coolant since primordial Fe is ∼ 400 times more abundant than Co.Additionally, or alternatively, the Co may be less sharply peaked as this may be the region corresponding to partial burning in quasi nuclear statistical equilibrium (qNSE) (DerKacy et al. 2023).
Although the MIR isolated Ni lines are weak and lowresolution, [Ni II], [Ni III], and [Ni IV] all exhibit a very similar line profile shape.Likewise, the [Co II] lines have a similar shape as [Co III] and the [Fe III] lines have a similar shape as [Fe II].Thus, to reduce degeneracy and make fitting many ions (and therefore parameters) more tractable, we use the same profile shape parameters (excluding amplitude) for all ionization states of the same element for Co, Fe, and Ni.We find that the Ni lines are slightly better fit by the Fe parameters than by Co, so we use the same profileshape parameters as [Fe II] for all Ni lines.In principle, we would not necessarily expect different ionization states of the same element to arise from the same emitting regions; however, Figure 4 shows that our constraints on these profile fits can reproduce the emission features very well from the optical through the MIR.This suggests that there is not significant ionization stratification in the IGEs (i.e., emission from different ionization states is coming from similar ejecta regions).Furthermore, the IGE lines all peak with a redshift of ∼ 3500 km s −1 ; the consistency of these profiles indicates that the IGEs are present in roughly the same ejecta regions.

Intermediate-Mass Elements: S, Ar, Ca
The emission from IMEs in SN 2022pul is very strong and displays more line-profile variation than the IGEs.Despite the variation, the IME emission-line profiles share similar shapes and all exhibit a blueshifted peak at ∼ −2000 km s −1 , the opposite direction as the IGEs.They also exhibit ionization stratification, shown in Figure 5.The higher ionization states of [S IV], [Ar III], and [Ca IV] all have a generally asymmetric boxy shape, whereas the lower ionization [Ar II] and [Ca II] appear to have an additional blueshifted, Gaussian-like component.We choose to fit the Ca II and Ar II features with the sum of a Gaussian and a slanted flat-top profile.Physically this could mean that the central, denser regions of the ejecta have lower overall ionization, requiring this additional Gaussian component (see Figure 6 for a visualization of an asymmetric shell and additional Gaussian component profile).
The right panel of Figure 4 shows our best-fit line profiles for the IMEs.The [Ar III] 8.99 µm line is isolated and exhibits a clear slanted flat-topped shape with Gaussian wings.This flat-topped Gaussian profile can arise from a thick shell of [Ar III] emission, where one side of the shell is thicker than the other, causing the slant.The viewing angle of this asymmetric shell also contributes to the observed slant (DerKacy et al. 2023).The parameters for our [Ar III] 8.99 µm fit are FWHM = 10, 300 ± 1100 km s −1 for the Gaussian wings, overall kinematic offset v offset = 700 ± 300 km s −1 , inner-shell velocity v min = 5400 ± 900 km s −1 , and offset of the shell thickness v offset, shell = −1000 ± 200 km s −1 .
As described above, we fix the fit parameters (excluding amplitude) for [Co II] 10.52 µm and [Ni III] 11.00 µm, and then we fit the remaining contribution to the feature at 10.5 µm as the [S IV] 10.51 µm line.This yields a slanted boxy profile with Gaussian wings for [S IV], with a similar width and slope as [Ar III], but a narrower inner-shell ra- dius.Its parameters are FWHM = 8200 ± 800 km s −1 , kinematic offset v offset = −1000 ± 200 km s −1 , inner-shell velocity v min = 2500 ± 700 km s −1 , and offset of the shell thickness v offset, shell = −1300 ± 300 km s −1 .
When we assume the peak of the [Co II] 10.52 µm line is redshifted (like all the relatively isolated IGE ions consistently display), we recover a blueshifted peak for [S IV] that agrees with the other IMEs.Because the IGEs are redshifted and the IMEs are blueshifted, we are able to distinguish the contributions to the 10.5 µm feature from both [Co II] 10.52 µm and [S IV] 10.51 µm.We do not at-tribute the full flat-topped, boxy shape of the feature at 10.5 µm to [S IV] 10.51 µm alone because the neighboring [Co III] 11.89 µm line is quite strong, and SN 2022pul has a low mean ionization state, so we do not expect the contribution from [Co II] 10.52 µm to be negligible.
The [Ca IV] 3.21 µm line can also be fit reasonably well by an asymmetric boxy profile with slightly larger width than [Ar III] and shallower slope.Its parameters are FWHM = 7500 ± 1300 km s −1 , kinematic offset v offset = −2400 ± 200 km s −1 , inner-shell velocity v min = 5900 ± 1100 km s −1 , and offset of the shell thickness v offset, shell = [Ar II] 6.98 µm and [Ca II] λλ7291, 7324 show similar line profiles (see Figure 3), and the [Ca II] line is helpful for interpreting the general shape of the low-resolution [Ar II] line.Figure 5 shows that the [Ar III] 8.99 µm line shape fits nicely with the [Ar II] 6.98 µm line if a blueshifted component is added.Thus, we fit [Ar II] 6.98 µm by summing a Gaussian with a forced contribution from the shape of the [Ar III].Since the [Ca II] shares a similar shape to [Ar II], with a Gaussian-like peak and a somewhat boxy red wing, we also fit it with the sum of a Gaussian and a forced contribution from the shape of the [Ar III].For [Ar II] the additional Gaussian component has FWHM = 9300 ± 1700 km s −1 and a kinematic offset v offset = −2300 ± 100 km s −1 and for [Ca II] the additional Gaussian component has FWHM = 4700 ± 200 km s −1 and a kinematic offset v offset = −1500 ± 100 km s −1 .Because the optical 0.73 µm region is complex with contribution from many lines, including [Fe II], [Ni II], and [Ar III], we fit all of these additional lines simultaneously with the [Ca II] doublet and find their contribution to be relatively small.While we cannot definitively rule out other line profile shapes for [Ca II] and [Ar II], it is clear that they exhibit an additional blueshifted peak of emission not seen in the higher IME ionization states.
The shapes of these IME line profiles and their ionization stratification point to the presence of broad, asymmetric shells of higher ionization IME material, and a somewhat narrower, blueshifted region of lower ionization IME material.

O and Ne
In the center-bottom panel of Figure 4, we display our bestfit line profiles for [O I] and [Ne II].The [O I] λλ6300, 6364 feature is remarkably narrow (FWHM = 2500 ± 100 km s −1 ), allowing us to resolve emission from both components of the doublet.We do not detect emission from either [O II] λ3727 or [O III] λλ4959, 5007; however, this region of the spectrum is dominated by low-ionization-state IGEs, complicating the interpretation of the blended features.The [O I] emission is best modeled by a single skewed Gaussian velocity component.Similar skewed (or "sawtoothed") emission has been observed in the [Fe II] and [Ca II] features of other 03fg-like SNe Ia (Taubenberger 2017;Siebert et al. 2023Siebert et al. , 2024)).This asymmetry could potentially be explained by either an intrinsically asymmetric ejecta distribution or dust extinction from within the SN ejecta.
If a large amount of dust were present in the ejecta, assuming it did not render the [O I] emission invisible, one would expect the redshifted emission to be more heavily extinguished than the blueshifted emission.This would lead to an emission profile whose total offset is preferentially blueshifted, appears steeper on the blueshifted side, and is more shallow-sloped on the redshifted side (Jerkstrand 2017).There is some evidence that the [Ca II] features of 03fg-like SNe are preferentially blueshifted (Siebert et al. 2024); however, one would expect more dust attenuation at shorter wavelengths.Furthermore, the [O I] in SN 2022pul is redshifted by 260 ± 10 km s −1 , so the systemic velocity of the ejecta (and dust) would need to be redshifted if dust extinction were the cause of the asymmetric line profile.This could be the case; however, we do not observe evidence for such a systemic shift in the other nebular features, and therefore favor an asymmetric ejecta distribution.
The [Ne II] 12.81 µm emission has a unique profile compared to the other elements.There is a narrow peak close to ∼ 0 km s −1 , indicating that it is present in the central regions of the ejecta along with [O I].However, most of the [Ne II] emission is likely not coincident with the [O I] ow- ing to its much larger ionization potential.The overall broad [Ne II] line suggests that it forms at large velocities similar to the IMEs and IGEs.We model the broad component of the [Ne II] line profile with a broad Gaussian on the red side (FWHM = 10000 ± 500 km s −1 ), and an asymmetric slanted flat-top on the blue side (FWHM = 7200 ± 600 km s −1 , kinematic offset v offset = 1400 ± 200 km s −1 , inner-shell velocity v min = 7000 ± 500 km s −1 , and offset of the shell thickness v offset, shell = 800 ± 100 km s −1 ).The one-sided asymmetric slant is tilted in the opposite direction as the tilt of the IMEs, and could arise from a half-shell or some type of missing bubble of [Ne II] on the blueshifted side (see Figure 4 for visualization of this profile).The overall shape of this profile does not appear as blueshifted as the other IMEs, and the differing locations of missing emission indicate that the distributions of Ne and the IMEs are distinctly different.
The highly asymmetric [Ne II] line shares some of the same projected velocities, but it cannot coexist with the IMEs and IGEs, indicating large-scale composition asymmetry.Additionally, the small-scale structure seen through-out the [Ne II] profile may be caused by the presence of a large number of randomly distributed blobs (e.g., see Jerkstrand 2017, Fig. 3); alternatively, it may be due to higher noise levels at the longest MIRI/LRS wavelengths.Other higher resolution optical and NIR emission lines, including the [O I] λλ6300, 6364 feature, do not show convincing evidence for small structure in the line profiles that would suggest clumping on large scales (see Section 4).In Section 5, we discuss the detection of [Ne II] 12.81 µm in relation to predictions from WD-WD violent merger models.

Geometric Interpretation: Distinct Ejecta Distributions
The difference in line profiles suggests differences in the geometric distributions of material in the ejecta.In order to more easily compare between the bulk emission coming from each group of elements (IGEs: Fe, Co, Ni; IMEs: S, Ar, Ca; O, Ne), rather than from particular elements and ions, we sum all of the line profiles across our full wavelength range within these three groups in velocity space and calculate the bolometric luminosity along the line of sight in bins of 500 km s −1 (dν L bol ).Each velocity bin represents a slab of emission on the plane of the sky, giving a 2D projected spatial distribution of the elements.Shown in the bottom right panel of Figure 6, the summed profiles of the IGEs, IMEs, and O+Ne are close to those of the individual profiles, but they erase the small variations seen from profile to profile.
Using the geometric interpretations of the line-profile components that we adopt in our fits, we construct cartoons for the projected 2D distribution of our three element groups to aid visualization of the distinct regions (see Figure 6).Assuming spherical symmetry, Gaussian emission distributions produce Gaussian line profiles, and slanted flat-topped (or boxy) profiles result from asymmetric shells of emission.We assume spherical symmetry for the individual components of each profile (resulting in overall axial symmetry for the 2D projected distributions) and stick to combinations of these two basic line-profile shapes, so our cartoons are oversimplifications.The real ejecta are probably asymmetric and complex, potentially including bubbles, mushroom shapes, and clumps.Pakmor et al. (2012) show the distribution of elements in their WD violent merger hydrodynamical simulation in their Fig. 2. Interestingly, our cartoons bear some similarities to these physically realistic distributions.
We stress that these cartoons are not to be taken literally and that they are 2D projections.Overlap in line-ofsight velocity (v LOS ) between the element groups does not necessarily imply that they are physically co-located in 3D space.Specifically, all of our distributions show emission at v LOS = 0 km s −1 ; however, this is a projection effect because, from theoretical models, there should be very little to no IGEs or IMEs in the most central regions where the [O I] emission is located.

Nickel-56 Mass
At the late phase of SN 2022pul (338 d post explosion), the 56 Ni (half-life of ∼ 6 d) has long since decayed away and the energy deposition that continues to power the spectral emission lines comes from the radioactive decay of 56 Co (halflife of ∼ 77 d).Locally-deposited gamma-rays and positrons from the decay can be absorbed by the IGEs and some IMEs in the IGE-rich regions, while the IMEs and unburned elements outside this region can absorb nonlocally deposited gamma-rays and secondary, low-energy photons.Thus, to estimate the initial mass of 56 Ni synthesized in the explosion, we account for the total flux by integrating the spectrum (F = F ν dν) in the range 0.4-14 µm.We then estimate the luminosity from the emission lines by correcting the flux to a distance of 16 Mpc (see Paper I, Siebert et al. 2024, for a discussion of distance to SN 2022pul).In the low-density nebular phase, however, most of the energy escapes.To account for this, we divide by the fraction of the total decay energy that is absorbed by the ejecta from the MERGER model (Blondin et al. 2023, and described further in Section 4), which has a value E absorbed /E total = 0.055.
Using the uncertainty in the distance (16 ± 2 Mpc) as the dominant source of uncertainty in our calculation, we estimate a bolometric luminosity (i.e., deposited 56 Co-decay power) at 338 d post explosion of 4.6 ± 1.2 × 10 41 erg s −1 , which corresponds to a mass of synthesized 56 Ni of 0.66 ± 0.17 M ⊙ (see equations from Nadyozhin 1994;Branch & Wheeler 2017, Section 5.4).This estimate may also suffer from uncertainties in reddening, absolute-flux calibration, gamma-ray escape fraction, and assumption of spherical symmetry, but these uncertainties should be smaller than the 25% uncertainty introduced by the distance.Our estimated 56 Ni mass is in rough agreement with the 0.6 M ⊙ of 56 Ni produced in the violent merger simulation by Pakmor et al. (2012) and the calculated 56 Ni masses for SN 2011fe, the canonical normal SN Ia (Bora et al. 2022, and references therein).and a violent merger of two sub-M Ch WDs (0.9 + 1.1 M ⊙ ) [MERGER; Pakmor et al. 2012]), which they compared to a spectrum of SN 2021aefx covering the wavelength range 0.3-14 µm from Kwok et al. (2023).These model spectra are generated assuming spherical symmetry, so the line profile shapes do not reflect ejecta asymmetries that might be seen in two-dimensional models.Only the MERGER model displayed a strong, centrally-peaked line due to [Ne II] 12.81 µm, resulting from the presence of Ne in the innermost ejecta layers of this model.None of the other models considered by Blondin et al. (2023) presented such a strong [Ne II] feature, and only the GCD model displayed a weak flat-topped line, whose intensity at line center was a factor of ∼ 3 smaller than observed in SN 2021aefx.Only by artificially setting a minimum Ne mass fraction of 10 −2 in the inner ejecta of the GCD model were they able to match the peak intensity and profile shape of the potential weak [Ne II] line in SN 2021aefx.Future analysis of subsequent observations of SN 2021aefx should clarify whether [Ne II] 12.81 µm is actually present in SN 2021aefx (Ashall et al., in prep.).
The violent merger model naturally predicts the presence of Ne in the inner ejecta layers, which is associated with the delayed and incomplete burning of the lower-mass secondary WD.Moreover, the higher density in layers ≲ 10, 000 km s −1 enhances the recombination rate there, which in turn results in a lower mean ionization state compared to the other explosion models.For comparing to the nebular spectrum of SN 2022pul, we recomputed the four models studied by Blondin et al. (2023) at a later time of 338 d post-explosion, using the same numerical setup.The full optical through MIR comparison of these models to the observed spectrum is displayed in Figure 7 and highlights of the results are shown in Figure 8.
Given the limited range of models from very different origins, we compare qualitatively rather than quantitatively.Of the four models, we favor the MERGER model for several    The lower ionization results in a better match to the numerous [Fe II] features in the 1-2 µm range, yet the optical [Fe III] complex around 0.5 µm is still too strong compared to SN 2022pul, as is the [Ni III] 7.35 µm line.This suggests the ionization in SN 2022pul is even lower than predicted in the MERGER model.We consider clumping to be a viable mechanism to further reduce the ionization state in this  model.Using the same approach as Blondin et al. (2023), we recomputed a clumped version of the MERGER model with a volume-filling factor f = 0.5 (i.e., resulting in a factor of two increase in density of the clumps).Clumpy ejecta can leave an impact on emission lines (e.g., nebular spectra of the Type Ib SN 1985F and of the Type IIb SN 1993J; Filippenko & Sargent 1989;Filippenko et al. 1994;Houck & Fransson 1996), but these would need to be large clumps with a large density contrast (i.e., dense shells in 1D or large clumps in 3D), whose physical origin in SN Ia ejecta is unclear.In our modeling approach the clumps are significantly smaller than the photon mean free path, so we see a global impact on the ionization state (through enhanced recombination) but no small structure in the line profiles.
As seen in Figure 9, this clumped model is able to better match the flux level in the [Fe III]-dominated complex, while maintaining a satisfactory match out to ∼ 3 µm.The lower ionization also results in a stronger [Ar II] 6.98 µm line and a weaker [Ni III] 7.35 µm line.However, it also results in an even weaker [Ca IV] 3.12 µm line compared to SN 2022pul.The same is true of [Ar III] 8.99 µm.Moreover, the [Ni II] 6.63 µm and [Ni III] 7.35 and 11.00 µm lines remain too strong.This suggests that the clumping is not uniform, affecting some regions more strongly than others.
The mismatch in lines from specific elements over several ionization stages suggests that abundance also plays a role.Here, a lower abundance of stable Ni and a higher abundance of Ar/Ca could help resolve some of these discrepancies.In particular, a lower abundance of stable Ni (∼ 0.03 M ⊙ in this MERGER model) seems warranted to reproduce the lines of Ni II and Ni III in the observed spectrum.Pakmor et al. (2013) show in their helium-ignited violent merger of a system of 0.9 + 0.76 M ⊙ WDs that the central density and the Ar/Ca abundance is higher in the innermost layers, while the stable Ni abundance is lower and more centrally concentrated than the merger from Pakmor et al. (2012), which the MERGER model is based on.Additionally, the amount of stable Ni in the merger model from Pakmor et al. (2012) scales with the metallicity of the progenitor system (Z = Z ⊙ in the model), so lower stable Ni abundance may imply the progenitor system had subsolar metallicity.There may be a variety of inner-ejecta properties that could be produced through violent merger scenarios that would better match the observations of SN 2022pul.
Despite the lower ionization of the clumped MERGER model, the [O I] λλ6300, 6364 emission remains very weak compared with SN 2022pul.This is not related to an abundance issue, since oxygen completely dominates the ejecta layers below ∼ 2000 km s −1 in the MERGER model (mass fraction ≳ 0.5).Radial compression does not change the optical depth, so the mass of inner material must be increased to boost the amount of gamma-ray energy intercepted by the inner, O-rich ejecta layers.We were able to match the intensity and width of the narrow [O I] feature seen in SN 2022pul by artificially increasing the density below ∼ 2000 km s −1 .In practice we followed the procedure from Dessart et al. (2020) and added a Gaussian component to the original density profile centered on the inner boundary of our spatial grid (v min ≈ 580 km s −1 ) with a width σ ≈ 600 km s −1 , resulting in an additional 0.1 M ⊙ in these layers (∼ 50% of which is O).The original mass of O (0.53 M ⊙ ) was all located in the inner ejecta and the total ejecta mass in the model is increased through this modification.The fractional O abundance in the innermost region is not particularly important, as long as [O I] λλ6300, 6364 remains the dominant coolant.
As seen from Figure 9, the narrow [O I] feature is well reproduced, while the rest of the spectrum remains largely unchanged compared to the original model.The physical origin of such an increase in density in the inner ejecta is likely the result of a higher-mass secondary WD (e.g., Pakmor et al. 2013).Other contributing effects could be how completely the secondary WD burned, a different C/O fraction of the initial secondary WD, and compression of the innermost layers through the 56 Ni bubble effect (e.g., Wang 2005;Dessart et al. 2021).The changes to the ejecta properties (density, temperature, mean ionization) in our modified merger model Fe i 1.187  where clumping and 0.1 M ⊙ central mass have been introduced are shown in Figure 10, and a detailed list of important emission lines in our modified merger model is given in Table 1.
While additional fine tuning of the abundances and densities of particular elements could help improve the match of any of the models, only the MERGER model can be modified to create the emergence of the narrow, central [O I] feature.This is because the other models predict essentially no O or Ne in the low-velocity central regions, while the MERGER model is actually dominated by O in the center (mass fraction ∼0.5).Fine tuning would not produce the emergence of the central O and Ne lines in the other models because they are not dominant coolants in this region, and adding a signif-icant central mass of O and Ne to the other models would be unphysical.(Taubenberger et al. 2019) and SN 2021zny (Dimitriadis et al. 2023)  the presence and diversity of oxygen and neon.Both turbulent pure deflagrations (Kozma et al. 2005;Fink et al. 2014) and violent WD mergers (Pakmor et al. 2012;Kromer et al. 2013;Pakmor et al. 2013) naturally predict neon and unburned oxygen present at low velocities.A violent WD-WD merger provides an explanation for the large-scale asymmetries in the nebular emission features described in Section 3, though a turbulent deflagration could also produce asymmetries.However, a pure deflagration would produce a lowluminosity explosion with different early-time spectroscopic properties, narrower nebular line profiles owing to the lower explosion energy, and weaker emission from IMEs.To check this, we ran the N100def model from Fink et al. (2014) with CMFGEN to 338 d post explosion (not shown here) and indeed found that the MERGER model is a significantly better match to SN 2022pul for the above reasons.We therefore favor a WD merger over a pure deflagration for SN 2022pul.
The large parameter space in WD mergers could help explain the diversity of 03fg-like SNe.The variation among objects may be explained by different WD mass ratios, by the specific explosion mechanism of one or both WDs, and by observations along different lines of sight.Recent studies of 03fg-like SNe Ia have favored WD merger models occurring within a dense C/O-rich CSM (Dimitriadis et al. 2022;Srivastav et al. 2023;Dimitriadis et al. 2023;Siebert et al. 2023) due to rapidly evolving light-curve bumps, lowionization nebular spectra sometimes with [O I], asymmetric and blueshifted nebular emission lines, and a larger fraction of flux in the NIR at late times than normal SNe Ia.Our observations of SN 2022pul enhance the evidence for this model for 03fg-like SNe Ia through clear dust detection, unusual light-curve rise shape, asymmetric line profiles, strong IMEs, and central O and Ne.
In the case of SN 2022pul, we also observe distinct ejecta components with different composition offset in their bulk velocities.The secondary WD may have been partially, but not completely, disrupted when the primary WD detonated.When burning front passes over the lower density secondary WD material, it too can burn, synthesizing O, Ne, and the IMEs (whereas the higher density primary WD would produce the IGEs) (Pakmor 2017).The degree of ejecta asymmetry will be smaller if the secondary is more compact when the primary is burned because the ashes of the primary will be less hindered by a less disrupted secondary (Pakmor 2017).
In SN 2022pul, we observe central O and Ne, very strong emission from IMEs, and a high degree of asymmetry in the nebular line features.Furthermore, our modified merger model demonstrates that additional, concentrated central mass (0.1 M ⊙ at ≲ 2000 km s −1 ) is needed to reproduce the [O I] λλ6300, 6364 emission.This very centrally peaked Notes: All wavelengths are given in air.Forbidden and semiforbidden transitions are noted using the appropriate brackets around the ion name.Wavelengths marked with a ' †' symbol denote transitions connected to the ground state.Ions with the same superscript correspond to transitions within the same multiplet.
and asymmetric [O I] emission is reminiscent of expectations of stripped material from a companion (Marietta et al. 2000;Botyánszki et al. 2018), and may indicate that the secondary WD left a significant amount of mass in the innermost region.The central, broad Ne and strong IMEs suggest that the secondary WD material was at least partially burned (the O may be a combination of burned C and unburned O, while unburned Ne is not expected for the lower-mass secondary).A higher mass, and therefore less disrupted secondary WD, could produce the strong IMEs and the centrally concentrated material.However, the high degree of asymmetry would suggest a more disrupted secondary, which perhaps could have left behind centrally concentrated material if it was only partially burned.This may be due in part to a viewing-angle effect where the ejecta appear more symmetric in the vertical direction, and more asymmetric in the plane of rotation (Bulla et al. 2016;Pakmor 2017).Previous measurements of 03fg-like SNe Ia have not shown strong continuum polarization (Tanaka et al. 2010;Cikota et al. 2019).More polarization data and construction of nebular spectra along different viewing angles from violent merger hydrodynamical simulations would help resolve this and improve our understanding of the 3D ejecta distributions giving rise to the asymmetric line profiles.
Another key part of understanding this progenitor scenario is its connection to the inferred presence of a dense C/O-rich CSM from the shape of the light-curve rise (Siebert et al. 2024).This requires considering scenarios where the explosion occurs either pre-merger during a dynamical interaction phase (commonly known as a "perimerger"), or scenarios where explosion occurs post-merger after the complete disruption of the secondary WD (Pakmor et al. 2010;van Kerkwijk et al. 2010;Pakmor et al. 2012;Moll & Woosley 2013;Guillochon et al. 2010;Raskin & Kasen 2013).
In both cases, simulations of the ejecta are highly asymmetric and the synthetic light curves cover a large range in luminosity depending on viewing angle (with the highest luminosity observed along the equator).This light-curve line-ofsight dependence might also vary over time because at early phases the interaction dominates over decay power.One beneficial consequence of the post-merger scenario is that the complete disruption of the secondary WD leads to the formation of a massive C/O-rich accretion disk.The early-time light curve of SN 2022pul is consistent with ejecta interaction with this kind of CSM (Noebauer et al. 2016;Siebert et al. 2024).More work is needed to understand if a violent merger is capable of producing enough CSM at the required distance to explain this kind of interaction signature.
The large amount of CSM (0.6 M ⊙ ) adopted in the Noebauer et al. ( 2016) interaction model is difficult to reconcile with the central oxygen and neon seen in SN 2022pul.If this much material from the secondary WD was pulled off into an exterior spherical shell, it seems unlikely that the secondary would have enough leftover mass to create the highdensity central O giving rise to the [O I] doublet and the strong emission from the IMEs.If the CSM is responsible for damping the kinetic energy of the explosion, leading to the boosted luminosity and low-velocity C II at early times, the CSM must cover a substantial fraction of the exploding system (Maeda et al. 2023).Howell et al. (2006) suggest that the extra binding energy associated with the additional mass may be enough to cause sufficiently slow ejecta velocity.However, Raskin & Kasen (2013) find that the amount of CSM stripped from the secondary WD in their models is small, less than 0.005 M ⊙ .If the dust emission arises from CSM, assuming a dust-to-gas ratio of 0.01, the dust mass for SN 2022pul given by Johansson et al. (2024) implies a CSM mass of 0.05-0.11M ⊙ , depending on composition.Further exploration is needed to address whether these lower CSM masses (or perhaps a different composition) could produce the interaction signature we observe.Hsiao et al. (2020) present evidence supporting the "coredegenerate" model for LSQ14fmg, where a near-M Ch WD interacts with the degenerate core of an AGB star within a common envelope stripped of its H and He through a superwind phase prior to explosion.We disfavor this scenario to explain SN 2022pul because the explosion is expected to be a spherical explosion of a single object after secular accretion.It therefore does not explain the observed central location of O and Ne or the large-scale ejecta asymmetries.It is tempting to speculate whether a modified version of this model, in which the AGB degenerate core material also undergoes thermonuclear fusion, could be developed for SN 2022pul.Perhaps a triple system could be invoked (Kushnir et al. 2013;Grishin & Perets 2022), where one WD or star is disrupted to become the surrounding CSM and precipitates the violent merger of the remaining two WDs.Rajamuthukumar et al. (2023) found that triple systems contribute significantly to the overall SN Ia rate.Again, more detailed modeling over a wide range of possible parameters and scenarios involving WD-WD violent mergers is needed to better understand the system giving rise to this unique SN.
We summarize our main observations and discussion points below.
• The presence of O and Ne in the central ejecta are hallmarks of the WD violent merger SN Ia model.
• Asymmetric nebular emission-line profiles show distinct distributions of IGEs, IMEs, and O and Ne in the ejecta that can be naturally explained by the violent merger model.
• A pure deflagration model is disfavored as it would produce lower luminosities, narrower nebular line widths, and weaker emission from IMEs than observed in SN 2022pul.
• Especially strong emission from IMEs could come from burning of a disrupted, lower-density secondary WD.
• Our spectrum of SN 2022pul is quite well matched by the MERGER model nebular spectrum computed by Blondin et al. (2023) from the violent WD-WD merger model of Pakmor et al. (2012).The agreement between the observed and model spectra is further improved by the addition of clumping and 0.1 M ⊙ to the central mass.Future observations of 03fg-like SNe Ia would benefit from JWST spectroscopy, especially in the higher resolution instrument modes which reach to even longer wavelengths and reveal additional lines.Complementary spectropolarimetry would give insight into line-profile asymmetries and constrain aspect angles (DerKacy et al. 2023).If this class of objects indeed arises from WD mergers, the observations should show a large degree of diversity in their asymmetric line profiles, and JWST observations could additionally establish a firmer connection to the dust surrounding or formed in these objects.This work is based on observations made with the NASA/ESA/CSA JWST as part of program #02072.We thank Shelly Meyett for her consistently excellent work scheduling the JWST observations, Sarah Kendrew for assistance with the MIRI observations, and Glenn Wahlgren for help with the NIRSpec observations.We also thank the anonymous referee whose suggestions significantly improved the paper.The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under National Aeronautics and Space Administration (NASA) contract NAS 5-03127 for JWST.Support for this program

Figure 3 .
Figure 3. [Ar II] and [Ca II] emission-line profiles in velocity space compared with [Co III] and [Fe II].The shapes of the [Ar II] and [Ca II] are quite consistent with each other, and the shapes of [Co III] and [Fe II] are nearly identical (taking into account the neighboring [Fe II] lines that create the red tail on the [Fe II] 1.26 µm feature).Notably, the Ar and Ca (IMEs) have a blueshifted peak with a more extended red tail, while the Co and Fe (IGEs) display an opposite trend, with a redshifted peak and a more extended blue tail.

Figure 4 .
Figure 4. Left: Fits to Co emission lines.The vertical dashed red line shows that the peak of the IGE distributions is highly redshifted at ∼ 3500 km s −1 .Neighboring IME lines are dashed blue lines while neighboring Fe lines are dashed gray lines.Center top: Similar to the Co fits, but for Fe (neighboring Co lines are dashed grey).Right: Fits to the IME emission lines.Neighboring IGE lines are shown in dashed red.The vertical dashed blue line shows the peak of the IME distributions, which is highly blueshifted at ∼ 2000 km s −1 .Center bottom: Fits to the optical [O I] λλ6300, 6364 doublet and the MIR [Ne II] 12.81 µm line.The [O I] is very central and narrow.The [Ne II] peaks at the same location as the [O I] (vertical dashed green line); however, it has an additional Gaussian-like red wing and a flat-topped-like blue wing slanted in the direction opposite that of the IMEs.

Figure 5 .
Figure 5.Comparison of the [Ar II] 6.98 µm and [Ar III] 8.99 µm lines, as well as the optical [Ca II] λλ7291, 7324 and NIR [Ca IV] 3.21 µm lines.The lower ionization states show a more prominent blueshifted peak (vertical dashed blue line), while the higher ionization states exhibit more of an asymmetric, slanted flat-topped shape.The bottom panel shows a comparison between the high ionization states of the IMEs.The [S IV] 10.51 µm appears more flat-topped owing to the neighboring [Co II] 10.52 µm line peaked toward the red.The different profile shapes between ionization states imply ionization stratification.−200 ± 200 km s −1 .Simultaneously fitting the nearby [Fe III] around 2.9 µm, we constrain the contribution to the [Ca IV] profile from [Fe III] 3.23 µm to be very minimal.Thus, most of the flux in this feature is from [Ca IV], though [Ni I] 3.12 µm may contribute to the narrow peak on the blue side of the [Ca IV] emission.[ArII] 6.98 µm and [Ca II] λλ7291, 7324 show similar line profiles (see Figure3), and the [Ca II] line is helpful for interpreting the general shape of the low-resolution [Ar II] line.Figure5shows that the [Ar III] 8.99 µm line shape fits nicely with the [Ar II] 6.98 µm line if a blueshifted component is added.Thus, we fit [Ar II] 6.98 µm by summing a Gaussian with a forced contribution from the shape of the [Ar III].Since the [Ca II] shares a similar shape to [Ar II], with a Gaussian-like peak and a somewhat boxy red wing, we also fit it with the sum of a Gaussian and

Figure 6 .
Figure 6.Top left, top right, bottom left: Cartoon projected 2D geometric interpretations of the IGE, IME, and O+Ne line profiles, assuming axial symmetry.Bottom right: Bolometric luminosity per projected velocity bin (500 km s −1 ) along the line of sight for the IGEs (Fe, Co, Ni), the IMEs (S, Ar, Ca), and O+Ne.There are clear differences in the overall line profiles from each element group, indicating distinct geometric distributions of material.This bolometric luminosity plot, derived from the spectral line profiles of the data (independent of the cartoon models), gives the true relative contributions of the different element groups without suggesting microscopic mixing.

Figure 7 .
Figure7.Spectra of our reference model set at 338 days post explosion compared to the dust and MIR-continuum-subtracted spectrum of SN 2022pul over the wavelength ranges (from top to bottom) optical (0.35-1 µm), NIR (1-5 µm), lower MIR (5-14 µm), and upper MIR (14-28 µm).The SN 2022pul spectrum has been corrected for redshift and MW extinction.The synthetic fluxes correspond to the assumed distance to SN 2022pul of 16 Mpc; they have not been rescaled or normalized in any way.The FX−Yµm label gives the fraction of the total optical to MIR flux (0.35-14 µm) for SN 2022pul emitted within the wavelength range of each plot.We include line identifications based on their Sobolev equivalent width, as inBlondin et al. (2023) (see Table1).Transitions connected to the ground state are marked with a " †" symbol.For consecutive lines within a multiplet (connected by a horizontal line), we give the g f -weighted mean wavelength of the transitions.

Figure 8 .
Figure 8. Similar to Figure 7 but highlighting the 0.7-2 µm and 10-14 µm regions where the models diverge substantially.The MERGER model matches the observed spectrum of SN 2022pul best, notably producing stronger [Ca II] near 0.73 µm, capturing the NIR [Fe II] lines extremely well, and predicting the strong observed [Ne II] 12.81 µm line that is not seen in the other models.However, there are also deficiencies to the MERGER model, such as the strength of the Ni lines being overpredicted (e.g., [Ni III] 11.00 µm).

Figure 9 .
Figure 9. Similar to Figure 7 but for our reference (blue) and modified (red) merger models.
5. DISCUSSION & CONCLUSIONS As discussed in Paper I (Siebert et al. 2024), the detection of [O I] emission in WD SNe is exceedingly rare.SN 2022pul is now the third 03fg-like SN Ia to exhibit [O I], along with SN 2012dn

Figure 10 .
Figure 10.Top left: Density at 338 d for MERGER model and our modified merger model with clumping and 0.1 M⊙ of mass added to the central (≲ 2000 km s −1 ) region.Bottom left: Temperature profiles.Top right: Mean Fe ionization profiles.The mean ionization is defined as Σiin i+ /Σin i+ , where n i+ is the number density of ionization stage i for Fe, such that a mean ionization of ∼ 2 indicates that Fe 2+ is the dominant stage.Middle right: Mean Ni ionization profiles.Bottom right: Mean O ionization profiles.

•
It is still somewhat unclear how the violent merger model can produce both the C/O-rich CSM described by Paper I (Siebert et al. 2024) and Paper III (Johansson et al. 2024), and the central O and Ne seen in the spectra.• Further work and more detailed models are needed to constrain the masses of the WDs involved in the merger.Our observations of SN 2022pul are a prime example of the importance of gathering spectra of WD SNe across the optical, NIR, and MIR wavelengths.The MIR observation of [Ne II] is complementary to the observation of [O I] in the optical; the optical [Ca II], NIR [Ca IV], and MIR [Ar II], [Ar III], and [S IV] mutually help interpret the IME distribution and ionization stratification; and the NIR [Fe II] and MIR [Co III] lines establish the distribution of the IGEs.Furthermore, the higher resolution ground-based optical and NIR observations were essential to interpret the lower resolution JWST observations, and isolated MIR spectral lines were key to unraveling the complex line blending in the optical and NIR.
Full optical+NIR+MIR comparison between the dust-continuum-subtracted spectrum of SN 2022pul at 338 rest-frame days postexplosion (texp = MJD 59785.3) at d = 16 Mpc, and SN 2021aefx at 270 d rest-frame days post-explosion at d = 18 Mpc from Kwok et al. (2023).The optical spectrum of SN 2022pul is from MMT/Binospec at a similar epoch of 332 d post-explosion and the NIR+MIR spectrum is from JWST/NIRSpec+MIRI.The prominent stable Ni emission lines in SN 2021aefx are weaker (relatively) in SN 2022pul, as is the [Co III]11.88 µm line.In SN 2022pul, the singly ionized emission lines are relatively strong compared with the higher ionization states.A strong [Ne II] feature at 12.81 µm is observed in SN 2022pul.The flux axis uses a non-linear (arcsinh) scale to better show all the features across a wide range of wavelength and Fν .
Figure 2. Comparison between the dust-continuum-subtracted MIR spectrum of SN 2022pul at 338 rest-frame days post-explosion (MJD 59785.3)and the MIR spectrum of SN 2021aefx at 340d rest-frame days post-explosion from DerKacy et al. (2023) scaled to the distance of SN 2022pul (16 Mpc).The prominent stable 58 Ni and Co emission lines in SN 2021aefx are somewhat weaker and have shifted to lower mean ionization state in SN 2022pul.The [Ar II] 6.98 µm line dominates the SN 2022pul spectrum, with strong [Ar III] 8.99 µm and [S IV] 10.51 µm as well.[Ne II] at 12.81 µm is strong, broad, and centrally peaked in SN 2022pul.All line profiles are considerably asymmetric.An inset is provided on a log scale to show the very weak [Ni IV] 8.41 µm and [Ni III] 11.00 µm lines.
Table1).Transitions connected to the ground state are marked with a " †" symbol.For consecutive lines within a multiplet (connected by a horizontal line), we give the g f -weighted mean wavelength of the transitions.
centrally located, whereas the other models predict very flat profiles for all ionization species, indicating that the IMEs are absent in the center; and (3) it is the only model consistent with the strong, centrally peaked [Ne II] 12.81 µm line.
. The velocity distribution of [O I] varies greatly among these SNe, while its offset is consistent with originating from the central regions of the ejecta in each case.Additionally, SN 2022pul exhibits [Ne II] peaked in the central regions and extending farther out as well.A viable explosion model should be able to explain

Table 1 .
Lines in the wavelength range 0.35-14 µm in our modified merger model (MERGER with clumping + 0.1 M⊙ central mass) at 338 days post explosion whose absolute Sobolev EW exceeds 5% of the largest absolute EW in the wavelength range 0.35-1 µm.