pith. sign in

arxiv: 2605.29038 · v1 · pith:GOEPSTFLnew · submitted 2026-05-27 · 🌌 astro-ph.HE

Nebular Fingerprints of a Violent White Dwarf Merger: 3D NLTE Modelling of Type Ia Supernovae

J. M. Pollin , R. Pakmor , L. J. Shingles , F. P. Callan , S. A. Sim , L. A. Kwok , C. E. Collins , F. K. Roepke
show 1 more author
A. L. McGarrity
This is my paper

Pith reviewed 2026-06-29 10:11 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords Type Ia supernovaewhite dwarf mergersnebular spectra3D radiative transferNLTE modelingviolent mergeroxygen emissionviewing angle effects
0
0 comments X

The pith

3D NLTE calculations of a violent white dwarf merger reveal [O I] emission from unburned secondary material and substantial viewing-angle effects.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper performs the first three-dimensional non-local thermodynamic equilibrium radiative transfer calculations on a violent merger of 1.1 and 0.7 solar mass white dwarfs. It demonstrates that multidimensional structure improves the ionization balance relative to one-dimensional models and produces nebular emission lines such as [O I] from unburned material that are absent in 1D. These results are compared to observed supernovae including SN 2021aefx and SN 2022pul to explore how nebular spectra can test progenitor channels. The calculations also show that stable nickel line strengths depend on both abundance and ionization state.

Core claim

By applying full NLTE radiative transfer including non-thermal electrons to the 3D hydrodynamical output of a violent merger, the authors show that multidimensional effects improve the ionization state and reveal [O I] features from unburned material associated with the secondary, while the model reproduces much of the spectrum of SN 2021aefx yet underpredicts [Ni II] and produces strong high-ionization stable-Ni lines; viewing-angle variations are large and distinct from D6-like double-detonation scenarios.

What carries the argument

The 3D hydrodynamical structure and composition of the 1.1 + 0.7 solar mass violent merger combined with full NLTE treatment of excitation, ionization, and non-thermal electrons.

If this is right

  • Multidimensional modelling improves the ionisation state and reveals features absent from 1D calculations, most notably [O I] from unburned material associated with the secondary.
  • The model reproduces much of the panchromatic spectrum of the normal SN 2021aefx but underpredicts [Ni II] while producing strong high-ionisation stable-Ni features.
  • Viewing-angle variation is substantial, with signatures distinct from D^6-like scenarios.
  • JWST nebular samples combined with multidimensional modelling can discriminate between progenitor channels.
  • The calculations suggest that SN 2022pul may require a similar merger configuration involving full disruption of the secondary or a more massive companion.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the input 3D structure is realistic, then nebular spectra can be used to map the distribution of unburned material left by the secondary.
  • The dependence of nickel lines on ionization state implies that abundance measurements from 1D models may need systematic corrections when applied to merger events.
  • Different viewing angles could produce observable spectral diversity that future multi-epoch observations might detect.

Load-bearing premise

The hydrodynamical structure and composition taken from the input violent-merger simulation accurately represent the physical conditions required for the subsequent NLTE radiative transfer to produce reliable spectra.

What would settle it

Absence of [O I] emission or lack of strong viewing-angle variation in nebular spectra of events similar to SN 2021aefx would contradict the predictions of this 3D violent-merger model.

Figures

Figures reproduced from arXiv: 2605.29038 by A. L. McGarrity, C. E. Collins, F. K. Roepke, F. P. Callan, J. M. Pollin, L. A. Kwok, L. J. Shingles, R. Pakmor, S. A. Sim.

Figure 1
Figure 1. Figure 1: Density of key species in three planes: 𝑌-𝑍 (top), 𝑋-𝑍 (middle), and 𝑋-𝑌 (bottom; cos( 𝜃 ) = 0, i.e., the merger plane). The arrow in the bottom left panel indicates the direction of 𝜙 = 0°. The colour scale indicates the logarithmic density of each species. Note this snapshot is at 270 days after explosion, apart from 56Ni, which is shown at 0.002 days. global asymmetry that is expected to produce observe… view at source ↗
Figure 2
Figure 2. Figure 2: 1D and direction-averaged 3D optical (top; ∼0.35–1 µm), NIR (middle; ∼1–3.5 µm), and MIR (bottom; ∼6–27 µm) spectra for the 3DViolentMerger and 1DViolentMerger models at 270 days post-explosion, compared to SN 2021aefx (Kwok et al. 2023). The observational data have been corrected for a redshift of 𝑧 = 0.005017 and for reddening due to host-galaxy extinction, 𝐸 (𝐵 − 𝑉)host = 0.097 mag (Hosseinzadeh et al. … view at source ↗
Figure 3
Figure 3. Figure 3: Nebular emission and absorption spectra for the 3DViolentMerger model at 270 days across all wavelength ranges, from top to bottom: optical, NIR, and MIR. The positive axis is colour-coded to indicate the emitting ions, based on each Monte Carlo packet’s thermal emission type. The negative axis shows the corresponding absorption contributions from each ion, which is only present in the optical region. The … view at source ↗
Figure 4
Figure 4. Figure 4: Same as [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Panchromatic optical, NIR, and MIR spectral time series of the direction-averaged spectrum (top panel) and viewing-angle orientations ranked by median absolute deviation (bottom panel) from the 3DViolentMerger model. Observational spectra of the normal SN 2022aaiq, SN 2021aefx, and SN 2024gy, and the 03fg-like SN 2022pul are shown; all have been corrected for redshift and reddening. The flux is normalised … view at source ↗
Figure 6
Figure 6. Figure 6: Ion populations for a slice through the 3DViolentMerger model (cos( 𝜃 ) = 0, i.e. the merger plane) at 270 days post explosion. Each panel shows a 2D slice of the 3D ejecta mapped into velocity space, illustrating the ion populations of key species: O i, Ne ii–iii, S ii–iv, Ar ii–iii, Ca ii, Ca iv, Fe i–v, Co i–v, and Ni i–v. Dashed circles indicate radial velocities of 5,000 km s−1 and 10,000 km s−1 . The… view at source ↗
Figure 7
Figure 7. Figure 7: Ion populations for the 1DViolentMerger model at 270 days post explosion. Each panel shows a 2D projection of the 1D ion population, assuming spherical symmetry, as a function of radial velocity for key species: O i, Ne ii–iii, S ii–iv, Ar ii–iii, Ca ii, Ca iv, Fe i–v, Co i–v, and Ni i–v. All panels share a common radial velocity scale, with inner and outer dashed circles marking velocities of 5,000 km s−1… view at source ↗
Figure 8
Figure 8. Figure 8: Spectra of the 3DViolentMerger model for different viewing angles at ∼300 days post-explosion, shown for the optical (top), NIR (middle), and MIR (bottom). The lines of sight are oriented around the merger plane (i.e., cos( 𝜃 ) = 0.0), where significant variation in the synthetic observables occurs. Observed spectra of SN 2021aefx (270 days; Kwok et al. 2023) and SN 2022pul (338 days; Kwok et al. 2024) are… view at source ↗
Figure 9
Figure 9. Figure 9: Same as [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Viewing-angle spectra for the 3DViolentMerger model at 270 days post-explosion in the merger plane (cos( 𝜃 ) = 0), with the corresponding orientations (𝜙 = 0–360°) indicated in the first panel of each row. The spectra are compared to the spectrum of SN 2021aefx (black; Kwok et al. 2023), which has been corrected for redshift and reddening. Prominent IGE and IME features are shown: [Fe iii] 0.466 µm, [Co i… view at source ↗
Figure 11
Figure 11. Figure 11: Same as [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Velocity shifts (top) of the [Fe ii] 1.257 µm, [Fe iii] 0.470 µm, [Co iii] 0.589 µm, [Ni iii] 7.348 µm and [Fe iii] 22.925 µm features for differ￾ent observer orientations at cos( 𝜃 ) = 0, 𝜙 = 0–360° (i.e., the merger plane), and the corresponding FWHM (bottom) for the 3DViolentMerger model. The [Ca iv] feature therefore provides a clear example of how asymmetric ion structure maps directly to the emergen… view at source ↗
read the original abstract

Binary systems composed of two carbon-oxygen white dwarfs (WDs) are a leading progenitor candidate for Type Ia supernovae. One widely discussed scenario is the dynamically driven double-degenerate double-detonation (D$^6$) of a sub-Chandrasekhar-mass WD binary, where detonations are triggered by dynamical interaction. However, some systems are expected to undergo violent mergers, in which the primary ignites through direct carbon ignition as the secondary strikes its surface. We present the first 3D nebular-phase radiative-transfer calculations of a violent merger, using a $1.1 M_\odot$ and $0.7 M_\odot$ sub-Chandrasekhar binary. Our simulations employ a full NLTE (non local thermodynamic equilibrium) treatment of excitation and ionisation, including non-thermal electron contributions. By comparing 1D and 3D realisations, we show that multidimensional modelling improves the ionisation state and reveals features absent from 1D calculations, most notably [O I] from unburned material associated with the secondary. The model reproduces much of the panchromatic spectrum of the normal SN 2021aefx, but underpredicts [Ni II] while producing strong high-ionisation stable-Ni features, illustrating that stable-Ni signatures depend not only on abundance, but also on ionisation state. Although the model does not reproduce the strong [Ar II] and [Ne II] emission observed in the 03fg-like SN 2022pul, our calculations suggest that this event may require a similar merger configuration, involving full disruption of the secondary or a more massive companion with more extensive burning. Finally, viewing-angle variation is substantial, with signatures distinct from D$^6$-like scenarios, suggesting that JWST nebular samples, combined with multidimensional modelling, can discriminate between channels.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper presents the first 3D NLTE nebular-phase radiative-transfer calculations for a violent white-dwarf merger using a 1.1 + 0.7 solar-mass binary. It compares 1D and 3D realisations to show that multidimensional modelling improves the ionisation state and reveals [O I] emission from unburned secondary material absent in 1D. The model reproduces much of the panchromatic spectrum of SN 2021aefx (underpredicting [Ni II] while producing strong high-ionisation stable-Ni features) but does not match the strong [Ar II] and [Ne II] in the 03fg-like SN 2022pul, while finding substantial viewing-angle variation distinct from D^6 scenarios.

Significance. If the hydrodynamical input is representative, this is a significant pioneering calculation demonstrating the necessity of 3D NLTE modelling (including non-thermal electrons) for accurate nebular spectra of violent mergers. The forward-modelling approach, explicit 1D–3D comparison, and discussion of viewing-angle diagnostics provide a concrete framework for using future JWST samples to discriminate progenitor channels.

major comments (1)
  1. [Abstract] Abstract: the claim that 3D modelling reveals [O I] from unburned secondary material and produces viewing-angle signatures distinct from D^6 scenarios is load-bearing on the adopted 1.1 + 0.7 Msun violent-merger hydrodynamical simulation. Systematic offsets in oxygen distribution, density profile or ignition geometry would directly alter the predicted line strengths and angular dependence independent of the NLTE solver; the manuscript should quantify this sensitivity.
minor comments (1)
  1. The abstract refers to 'panchromatic spectrum' comparisons; specifying the exact wavelength coverage and which ions dominate the fit would improve clarity.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for their positive assessment and constructive feedback. We address the major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that 3D modelling reveals [O I] from unburned secondary material and produces viewing-angle signatures distinct from D^6 scenarios is load-bearing on the adopted 1.1 + 0.7 Msun violent-merger hydrodynamical simulation. Systematic offsets in oxygen distribution, density profile or ignition geometry would directly alter the predicted line strengths and angular dependence independent of the NLTE solver; the manuscript should quantify this sensitivity.

    Authors: We agree that the specific predictions for [O I] emission and viewing-angle signatures depend on the details of the adopted 1.1 + 0.7 M⊙ hydrodynamical simulation. Our study isolates the effects of 3D NLTE radiative transfer by comparing 1D and 3D realisations of the same input structure, demonstrating that multidimensional effects improve ionisation and produce [O I] from secondary material. The viewing-angle variations reflect the asymmetric 3D geometry of this violent merger, which differs from D^6 scenarios. However, a quantitative assessment of sensitivity to changes in oxygen distribution, density profile or ignition geometry would require additional hydrodynamical models, which lies beyond the scope of this first 3D NLTE calculation. We will revise the abstract (and add a corresponding statement in the discussion) to explicitly note that the results are for this specific binary configuration and that systematic exploration of hydrodynamical variations is an important topic for future work. This will appropriately qualify the claims without altering the core demonstration of multidimensional effects. revision: partial

standing simulated objections not resolved
  • Quantification of the sensitivity of the [O I] line strengths and viewing-angle signatures to systematic variations in the hydrodynamical input (oxygen distribution, density profile or ignition geometry)

Circularity Check

0 steps flagged

No significant circularity; forward modeling from fixed hydro input

full rationale

The paper executes a one-way forward chain: a pre-existing 3D hydrodynamical merger simulation (1.1 + 0.7 M⊙) supplies density, composition and velocity structure; NLTE radiative transfer is then applied without any parameter adjustment or fitting to the target supernovae spectra. Comparisons to SN 2021aefx and SN 2022pul are post-hoc evaluations only. No equations or claims reduce the output spectra to the input abundances by construction, no fitted parameters are relabeled as predictions, and no self-citation chain is invoked to justify the central NLTE results. The 3D-versus-1D ionisation differences arise directly from the geometry of the supplied hydro model rather than from any circular redefinition.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim depends on the accuracy of the input 3D hydro merger model and standard assumptions in the NLTE radiative transfer code; no new entities are introduced.

free parameters (1)
  • Binary component masses = 1.1 M_sun primary, 0.7 M_sun secondary
    1.1 and 0.7 solar masses chosen to represent a sub-Chandrasekhar violent merger; selection affects composition and is not derived from the spectral data.
axioms (1)
  • domain assumption Hydrodynamical merger simulation supplies correct 3D density, velocity, and abundance structure for radiative transfer input.
    Radiative transfer calculations are performed on top of this prior hydro output.

pith-pipeline@v0.9.1-grok · 5932 in / 1335 out tokens · 43603 ms · 2026-06-29T10:11:42.531139+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

  1. [1]

    J., Collins C

    ARTIS Collaboration Shingles L. J., Collins C. E., Holas A., Callan F., Sim S.,2024,artistools(v2024.12.9),doi:10.5281/zenodo.14337284,https: //doi.org/10.5281/zenodo.14337284 ARTIS Collaboration et al., 2025, artis (v2025.08.01), doi:10.5281/zenodo.16684298,https://doi.org/10.5281/ zenodo.16684298 6 https://github.com/artis-mcrt/artistools/ 3D Nebular Si...

    work page doi:10.5281/zenodo.14337284 2024