Open Access
Issue
A&A
Volume 686, June 2024
Article Number A227
Number of page(s) 7
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202449637
Published online 14 June 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model.

Open access funding provided by Max Planck Society.

1. Introduction

Despite decades of observational and theoretical efforts, we still do not know the progenitor systems and the explosion mechanism of Type Ia supernovae. These questions are deeply connected but need to be answered individually. A detailed overview of the current state of the field is available in recent reviews (e.g. Ruiter 2020; Liu et al. 2023).

Spectra in the nebular phase of Type Ia supernovae more than 100 d after the explosion directly probe the centre of the ejecta, when they become optically thin (Ruiz-Lapuente & Lucy 1992). They, as well as supernova remnants (Seitenzahl et al. 2019), are thereby in principle well suited to distinguish explosion models.

The arrival of JWST and its ability to obtain continuous spectra of Type Ia supernovae in the nebular phase all the way through the near- and mid-infrared enables us not only to observe isolated nickel, cobalt, and iron lines in nebular spectra, but also to measure and compare their detailed line shapes (Kwok et al. 2023; DerKacy et al. 2023; Blondin et al. 2023; Siebert et al. 2024; Kwok et al. 2024). This new window into the inner ejecta of Type Ia supernova explosions may allow us to finally understand their explosion mechanism and solve the progenitor question.

This paper revisits existing explosion simulations of the most popular Chandrasekhar-mass and sub-Chandrasekhar-mass explosion scenarios currently in the running to explain normal Type Ia supernovae. We show that the ejecta of all these models are inherently multidimensional, and that both the explosion and very likely the radiation transport have to be simulated in more than one dimension to correctly interpret recent JWST observations.

We also show that the structure of the ejecta, and in particular the distribution of the stable and radioactive nickel in them, are characteristic of the different explosion scenarios. Both can in principle be directly probed with nebular spectra of Type Ia supernovae. Therefore, new and future JWST nebular spectroscopic observations of Type Ia supernovae will likely be very valuable for constraining and understanding their physics.

We start with a brief summary of the explosion physics of Type Ia supernovae in Sect. 2. We then compare different explosion models of Chandrasekhar-mass white dwarfs in Sect. 3 and sub-Chandrasekhar-mass white dwarfs in Sect. 4, and quantify the similarities and differences of the ejecta of different explosion scenarios. In Sect. 5 we discuss implications for obtaining synthetic observables and for the interpretation of observations. We conclude with a summary of our main results and a brief outlook in Sect. 6.

2. The physics of Type Ia supernova explosions

To understand the connection between explosion mechanism and production of radioactive and stable nickel and the need for multidimensional models, we first took a step back and looked at what we know about Type Ia supernovae explosion physics.

There is a broad consensus that the progenitor systems of Type Ia supernovae are binary systems with at least one massive carbon-oxygen white dwarf. Since isolated white dwarfs are inert, we believe that interaction with a companion star leads to the explosion. Upon ignition, unstable nuclear burning disrupts the massive white dwarf and produces about a solar mass of ejecta with significant amounts of radioactive 56Ni, which decays via 56Co to 56Fe and powers the light curve of a Type Ia supernova (Diehl et al. 2015; Churazov et al. 2015). Here, we focused specifically on normal Type Ia supernovae that make up ∼70% of all Type Ia supernovae and produce a typical amount of 0.6 ± 0.3 M of radioactive 56Ni (Li et al. 2011; Taubenberger 2017).

To connect progenitor systems with observables of explosions and understand the production of radioactive and stable nickel, we needed to look at the explosion mechanism and the nuclear burning that powers the explosion. Unstable nuclear burning of degenerate matter proceeds as a thin burning front separating burning ashes and unburned fuel1. They come in two fundamentally distinct modes (see e.g. Röpke 2017, for a review).

Deflagration flames propagate subsonically via thermal conduction. Therefore, the unburned fuel reacts to the approaching flame before it is burned. Moreover, the whole white dwarf responds to the energy release from the deflagration and expands on a timescale comparable to the timescale on which the deflagration burns the white dwarf.

Detonation fronts are in essence shock waves sustained by energy released from nuclear burning immediately behind the shock front. They propagate supersonically with respect to the sound speed of the fuel. A detonation burns the fuel before it can react to the approaching burning front. Once started, a detonation burns the whole white dwarf essentially instantaneously, faster than the sound crossing time of the white dwarf.

An important consequence is that fuel, as well as the white dwarf as a whole, can react to the approaching deflagration, but not to the approaching detonation. When a deflagration flame is ignited in a white dwarf, the energy release from nuclear burning causes the white dwarf to expand on a timescale comparable to the timescale on which the deflagration flame propagates through the white dwarf. Therefore, the density of a parcel of fuel when it is burned is typically significantly lower than its initial density was at the instant when the deflagration ignited. In contrast, a detonation burns all of the white dwarf’s material essentially at the density it had at the moment the detonation formed.

The density of the fuel when it is burned by either a deflagration or a detonation is crucial because it principally determines the composition of the ashes. In particular, this includes the amount of neutronisation. Carbon-oxygen white dwarfs with solar metallicity are close to a Ye = 0.5, that is, they have an equal number of protons and electrons. At densities ρ ≳ 109 g cm−3 weak reactions become important enough to substantially neutronise the fuel when it is burned and decrease Ye (Seitenzahl et al. 2009). This systematically changes the typical isotopes created in the nuclear burning.

Unfortunately, the ignition of white dwarfs is one of the unsolved questions in the puzzle of the explosion mechanism of Type Ia supernovae. Despite fundamental physical differences, different explosion models are not easily distinguished by their observables around peak brightness (Röpke et al. 2012), because the composition of the outer layers of the ejecta are typically similar between different explosion models and the observables are rather insensitive to most of the detailed structure of the inner, optically thick part of the ejecta.

Because neutronisation during nuclear burning strongly depends on density, the ratio between radioactive 56Ni (Ye = 0.5) and stable 58Ni (Ye < 0.5) becomes a sensitive probe of the burning conditions of the inner parts of the white dwarf. Moreover, because radioactive and stable nickel are both produced in substantial amounts in Type Ia supernovae they offer an important observational diagnostic (see e.g. Blondin et al. 2022).

Radioactive 56Ni powers the light curve of Type Ia supernovae around maximum light via radioactive decay to 56Co and further to 56Fe. In the nebular phase, 100 d or later after the explosion, it has already completely decayed. Therefore, nickel detected in nebular spectra or supernova remnants is completely dominated by stable 58Ni that was synthesised directly in the explosion. In contrast, cobalt and iron in the nebular phase are dominated by decay products of 56Ni, and therefore trace the radioactive 56Ni created in the explosion. Observing nickel, cobalt, and iron in nebular spectra or remnants of type Ia supernovae thus directly informs us about the nucleosynthesis conditions in the centre of the explosion itself. It is therefore crucial to understand what to expect for different explosion scenarios to interpret observations.

Here we analysed the structure of the inner ejecta of multidimensional explosion simulations of the most popular Chandrasekhar-mass and sub-Chandrasekhar-mass explosion scenarios currently in the running to explain normal Type Ia supernovae (Ruiter 2020; Liu et al. 2023) with a focus on the distribution of radioactive and stable nickel in the ejecta.

3. Chandrasekhar-mass models

Chandrasekhar-mass explosion models start from a carbon-oxygen white dwarf that is accreting material from a companion star and is approaching the Chandrasekhar mass. The high pressure leads to nuclear reactions in its centre (Cameron 1959), which slowly increase its temperature and drive large-scale convection (e.g. Höflich & Stein 2002). At some point, convection becomes insufficient to transport away the energy generated from nuclear burning (e.g. Woosley et al. 2004). At this time, a local runaway leads to unstable nuclear burning and proceeds either as a deflagration or a detonation. This ignition is hard to model in detail, but probably occurs in a single spot slightly off-centre (Zingale et al. 2011; Nonaka et al. 2012; Malone et al. 2014). Importantly, at this point the central density of the white dwarf is typically higher than 109 g cm−3, and almost all of its mass sits above 107 g cm−3 (e.g. Gasques et al. 2005; Lesaffre et al. 2006; Piersanti et al. 2022).

At the moment, we cannot make an ab initio statement about the nature of the initial burning. Ignition of either a deflagration or a detonation is in principle possible. If the initial ignition led to a detonation, all the mass of the white dwarf would burn at its density at ignition. Because most of the mass resides at a density larger than 107 g cm−3, the burning products are almost exclusively iron group elements, and significant production of intermediate-mass elements like silicon does not occur (Arnett 1969; Seitenzahl & Townsley 2017). Such an explosion has never been observed even though it would be very bright.

Therefore, the initial ignition has to proceed as a deflagration. It initially burns at densities close to 109 g cm−3 in the centre of the Chandrasekhar-mass white dwarf and its energy release causes the white dwarf to expand. Because the deflagration starts burning at high densities, its ashes always contain a high fraction of stable nickel. The deflagration ashes rise buoyantly and drive turbulence in the whole white dwarf. For the most likely off-centre ignition (Nonaka et al. 2012), only part of the white dwarf is burned and ejected, typically leaving a bound remnant with a mass of ∼1 M. The ejecta are dominated by iron group elements. Ejecta of such an explosion have been shown to provide a compelling match to the peculiar subclass of 02cx-like Type Ia supernovae (Kromer et al. 2013; Magee et al. 2016; Kawabata et al. 2021; Lach et al. 2022a; Dutta et al. 2022; Singh et al. 2024). For the energetically most optimistic central isotropic ignition the whole white dwarf does explode, but the ejecta do not at all resemble normal Type Ia supernovae (Röpke et al. 2006, 2007).

However, it seems possible to start a detonation after the energy release from the initial deflagration has expanded the white dwarf significantly. When a significant fraction of the white dwarf has dropped below a density of 107 g cm−3, a detonation that burns the white dwarf essentially instantaneously at the density it has then, produces substantial amounts of intermediate mass elements.

Two fundamentally different, physically plausible paths have been proposed to lead to a detonation. In the delayed detonation model (Khokhlov 1991), strong turbulence in the white dwarf leads to a transition of the deflagration flame into a detonation. In the gravitationally confined detonation model (Plewa et al. 2004), a bubble of deflagration ashes leaves the surface of the white dwarf and drives a flow around it. The flow of deflagration ashes converges on the opposite side and starts a detonation via compression.

In both cases, the ashes of the deflagration have already risen substantially when the detonation starts, and the central part of white dwarf at this time is dominated by unburned material. The detonation then burns the rest of the white dwarf essentially instantaneously and produces some radioactive nickel with little stable nickel, because the lower density of the fuel avoids additional neutronisation. The detonation ashes then occupy the central parts of the final ejecta.

We show the structure of stable and radioactive nickel in the ejecta of Chandrasekhar-mass explosion models quantitatively in Fig. 1. The figure shows the mass distribution of radioactive and stable nickel in velocity–density space 100 s after the explosion when the ejecta are fully in homologous expansion. We focused on three representative, well-studied, Chandrasekhar-mass explosion models that all produce about 0.6 M of radioactive 56Ni, typical of a normal type Ia supernova. They are N100, a delayed detonation model (Seitenzahl et al. 2013), where turbulence leads to a detonation after 0.9 s, 82_d1.0, a gravitationally confined detonation (Lach et al. 2022b), where a detonation forms from a converging shock on the surface of the white dwarf 6.4 s after the ignition of the initial deflagration, and DDC10, a 1D model for a delayed detonation (Blondin et al. 2013) where the deflagration is artificially turned into a detonation 1.7 s after the ignition of the deflagration when the density at the flame drops below 2.7 × 107 g cm−3.

thumbnail Fig. 1.

Distribution of radioactive 56Ni (top row) and stable 58Ni (bottom row) in density–velocity space 100 s after the explosion. Each histogram is normalised individually to the total mass of the type of nickel it shows. The columns show different Chandrasekhar-mass explosion models: a 3D delayed detonation (N100; Seitenzahl et al. 2013), a 3D gravitationally confined detonation (r82.d1.0; Lach et al. 2022b), and a 1D delayed detonation (DDC10; Blondin et al. 2013). The solid line in the left and middle panels shows the spherically averaged density profile of the total mass of the ejecta. The dashed lines indicate if nickel was synthesised in the deflagration or in the detonation. In the 3D models the deflagration ashes are located at higher velocities, i.e. further out in the ejecta, than the detonation ashes. In the 1D model this hierarchy is inverted. In the 3D models nickel is distributed over a wide range of densities at a fixed velocity, so spherical averaging 3D ejecta significantly changes their physical properties.

3.1. The separate location of radioactive and stable nickel in the ejecta

In all Chandrasekhar-mass models shown in Fig. 1, radioactive and stable nickel are clearly separated. In the 3D models, the nickel created in the deflagration ends up at larger velocities than the nickel created in the detonation. This is a direct consequence of the buoyant rise in the deflagration ashes that drive Rayleigh–Taylor instabilities, leading to volume-filling turbulence in white dwarfs. The centre of the white dwarf is then filled with fuel when the detonation forms, and the central ejecta are made up exclusively of detonation ashes. The 1D model has a fundamentally different ejecta structure with an inverted hierarchy of detonation and deflagration ashes. This inversion is a direct result of the 1D nature of the model. It does not allow the model to describe the buoyant rise of deflagration ashes. Instead, they remain inside the detonation ashes that are created at larger radii after the deflagration has turned into a detonation. However, the buoyant rise of the deflagration ashes in 3D is unavoidable for all physically plausible ignition conditions.

It has been suggested that a large-scale smooth initial deflagration front combined with an early ignition of a detonation could technically prevent the inversion of the ashes (Leung & Nomoto 2018; Kobayashi et al. 2020). Crucially though, these conditions seem impossible to reach from physical ignition scenarios, because of buoyancy and Rayleigh–Taylor instabilities (Zingale et al. 2011; Seitenzahl et al. 2013). Even implausibly strong magnetic fields cannot suppress these instabilities on large scales (Hristov et al. 2021). Moreover, the lack of turbulence in those models removes any physically plausible reason for the turbulent ignition of a detonation.

3.2. Density variations in the ejecta at a fixed velocity

As shown in Fig. 1, in all 3D Chandrasekhar-mass explosion models, nickel is spread over a significant range of densities at fixed velocity. The deflagration ashes in the N100 delayed detonation model are located in dense blobs, at a significantly higher density (up to factor of two) than the spherically averaged density of the total ejecta at the same velocity. Most of the detonation ashes in contrast sit at slightly lower densities than the spherical average of the ejecta, because of a tail to higher densities at low velocities. The 82_d1.0 model shows a similar tail to nickel at higher densities that is caused by large-scale asymmetries in the white dwarf after the deflagration bubble has risen to the surface. We note that microscopic mixing on unresolved scales will not change any of these conclusions, because the scales on which the deflagration ashes and detonation ashes are separated are well resolved in the simulations. The detonation ashes occupy comparable densities in the 3D models and the 1D model, but the 1D model has much higher densities in the inner deflagration ashes. A possible reason is that the deflagration ashes have to lift all of the unburned material of the white dwarf to expand. This challenge arises because, in the 1D model, the deflagration ashes cannot rise buoyantly and pass beyond the unburned material.

4. Sub-Chandrasekhar-mass models

The most important difference between explosion models of sub-Chandrasekhar white dwarfs and Chandrasekhar-mass white dwarfs is the density profile of the white dwarf at ignition (Seitenzahl & Townsley 2017). The central density of sub-Chandrasekhar-mass white dwarfs is substantially lower (a 1.2 M carbon-oxygen white dwarf has a central density of only ∼2 × 108 g cm−3). Deflagrations in these sub-Chandrasekhar-mass white dwarfs would cause the white dwarf to expand quickly, and only produce tiny amounts of 56Ni. Therefore, only immediate detonations are relevant for sub-Chandrasekhar-mass white dwarfs. The detonation at low densities avoids any significant neutronisation during nuclear burning. Thus, the amount of stable nickel produced in sub-Chandrasekhar-mass models is typically comparable to the amount created by the detonation in delayed detonation Chandrasekhar-mass models, but significantly lower than in their deflagration.

Sub-Chandrasekhar-mass explosion scenarios can be differentiated by their ignition mechanism. In the classic double detonation scenario, a sub-Chandrasekhar-mass carbon-oxygen white dwarf accretes helium until its helium shell reaches sufficient conditions in the densest regions to lead to unstable nuclear burning that starts a detonation. The detonation in the helium shell then wraps around the white dwarf, sending a shockwave into it from all directions. When the shockwave converges in its core (or potentially when it interacts with itself on the opposite site of the ignition point or when it first hits the carbon-oxygen core) it starts a carbon detonation that burns the whole white dwarf (Bildsten et al. 2007; Fink et al. 2007, 2010; Shen & Moore 2014; Gronow et al. 2021).

If the companion is another carbon-oxygen white dwarf that donates helium from a thin shell, the helium detonation can form dynamically (Guillochon et al. 2010; Pakmor et al. 2013, 2021, 2022). This happens only when the secondary white dwarf is so close that it is about to be disrupted. At this time, the material in the accretion stream can be dense enough to be degenerate. Its impact on the surface of the primary white dwarf starts a detonation in the helium shell of the primary white dwarf like a hammer hitting an anvil. This scenario works for significantly smaller helium shells on the primary white dwarf, but for the purpose of this paper, it is essentially identical to the traditional double detonation scenario if the secondary white dwarf survives (Boos et al. 2021; Shen et al. 2021b; Pakmor et al. 2022). However, the secondary white dwarf can probably ignite as well via a similar mechanism, where the impact of the ejecta of the exploding primary white dwarf starts a detonation in the helium shell of the secondary. In this case it also explodes, but with a delay of a few seconds, significantly changing the inner ejecta of the explosion (Pakmor et al. 2022; Boos et al. 2024). Finally, if there is little to no helium present in a binary system of two carbon-oxygen white dwarfs, the disruption of the secondary white dwarf when the system merges and the impact this has on the primary white dwarf might directly lead to a carbon detonation on the surface of the primary white dwarf in the violent merger scenario (Pakmor et al. 2010, 2011, 2012b). In this case, both white dwarfs are burned by a detonation with only a short delay.

Figure 2 shows the mass distribution of radioactive and stable nickel in velocity–density space 100 s after the explosion for three typical sub-Chandrasekhar-mass explosion models. At this time the ejecta are fully in homologous expansion.

thumbnail Fig. 2.

Distribution of radioactive 56Ni (top row) and stable 58Ni (bottom row) in density–velocity space 100 s after the explosion. Each histogram is normalised individually to the total mass of the type of nickel it shows. The columns show different sub-Chandrasekhar-mass merger explosion models: a violent merger (1109; Pakmor et al. 2012a) and a helium ignited merger in which only the primary white dwarf (oneExp; Pakmor et al. 2022) explodes or both white dwarfs explode (twoExp; Pakmor et al. 2022). The solid line shows the spherically averaged density profile of the total mass of the ejecta. The distributions of radioactive and stable nickel in the ejecta are highly correlated. Only the oneExp model is close to spherical symmetry. In both the 1109 and twoExp models the secondary white dwarf explodes as well, compressing the ejecta of the primary white dwarf that contain primarily nickel. In the 1109 model the secondary explosion compresses the nickel to much higher densities than in the twoExp models because it happens with a delay of only 1 s compared to 4 s in the twoExp model.

The 1109 model is a violent merger of a 1.1 M and a 0.9 M white dwarf (Pakmor et al. 2012a). The carbon detonation starts on the surface of the primary white dwarf when it is hit by the secondary white dwarf. It first burns the primary, then, with a delay of ∼1 s, also the secondary white dwarf. Only the detonation that burns the primary white dwarf produces nickel; the central density of the secondary is too low to burn to iron group elements in the detonation in this case (but we note that the production of some nickel in the secondary would occur for an initially more massive secondary white dwarf). The ejecta of the secondary white dwarf expand into the central ejecta of the primary and compress them. Because of the small delay of only ∼1 s, the inner ejecta of the explosion of the primary white dwarf, which are nickel-rich, are significantly compressed. Most nickel in the final ejecta ends up at densities significantly larger than the average density of the ejecta at the same velocity, because a significant amount of material of the ejecta of the secondary ends up at similar velocities. The final ejecta of 1109 at velocities smaller than 2000 km s−1 contain essentially no nickel and are made up purely by the ashes of the secondary white dwarf.

The oneExp model is a double detonation explosion in a binary system of a 1.05 M and a 0.7 M white dwarf with thin helium shells of 0.03 M on each of them (Pakmor et al. 2022). Helium accretion on the primary white dwarf dynamically starts a helium detonation just a few orbits before the secondary would be disrupted and otherwise merge with the primary white dwarf. The helium detonation triggers a secondary carbon detonation in the core of the primary. In the oneExp model, the secondary white dwarf survives the explosion, by assumption (cf. twoExp below). The ejecta of the explosion are close to spherical, uniquely so among all 3D explosion models we consider here. The inner ejecta of the oneExp model are practically identical to classic double detonation models and simplified models of artificially centrally ignited sub-Chandrasekhar-mass carbon-oxygen white dwarfs (Sim et al. 2010; Shen et al. 2018, 2021a).

The twoExp model starts out identical to the oneExp model, but after the ejecta of the explosion of the primary white dwarf hit the secondary, a similar double detonation mechanism is assumed to occur there as well: the impact of the ejecta of the primary detonates its helium shell and causes a carbon core detonation, which destroys the secondary (Pakmor et al. 2022; Boos et al. 2024). Similar to the 1109 violent merger model, the ejecta of the explosion of the secondary white dwarf do not produce any iron group elements and expand into the centre of the ejecta of the primary white dwarf, albeit with a much longer delay of ∼4 s (compared to ∼1 s for the violent merger model). They still significantly compress the nickel-rich inner ejecta of the explosion of the primary white dwarf, but much less so than in the 1109 model. They also dominate the ejecta out to larger velocities (about 5000 km s−1 compared to the 2000 km s−1 in 1109).

5. Implications for Type Ia supernova nebular spectra and comparison to observations

The result that the ejecta of Type Ia supernovae are asymmetric is not new, and has been discussed in detail already in the context of Chandrasekhar-mass delayed detonation models (Gerardy et al. 2007; Kasen et al. 2009; Maeda et al. 2010) and violent mergers (Bulla et al. 2016). Our results emphasise that this is true in general for essentially all currently discussed explosion models of Type Ia supernovae. These asymmetries are fundamental to all 3D explosion models and faithfully simulating Type Ia supernova explosions crucially requires multidimensional simulations. In addition, we directly show that spherically averaging the ejecta of multidimensional explosion simulations before computing synthetic observables fundamentally changes the ejecta structure of the explosion and that 1D Chandrasekhar-mass explosions models are inherently unphysical. Therefore, extreme care should be taken when observational results are based on 1D explosion models (e.g. DerKacy et al. 2023), or for 1D radiation transfer simulation of multidimensional explosion simulations.

The detonation ashes of all 3D Chandrasekhar-mass models and the oneExp model end up at similar densities in the final ejecta. Regions that have both higher density and large mass fractions of radioactive nickel are only present in the 1109 model (mass density much higher than the Chandrasekhar-mass models) and the twoExp model (factors of a few higher). In both sub-Chandrasekhar-mass models, the ejecta of the primary white dwarf are compressed to higher densities by the impact of the ejecta of the secondary white dwarf that explodes slightly later. The densities directly influence the ionisation states of nickel and iron, which are crucial for the formation and intensity of specific lines in the nebular spectra (Flörs et al. 2020; Blondin et al. 2023).

Different timings between the explosions of the primary and the secondary white dwarf are likely possible for merger models, at least in the whole range between 1 s (1109) and 4 s (twoExp). Nebular spectra covering a wide range of wavelengths through the infrared might allow us to infer the density of nickel and iron in the ejecta and therefore constrain properties of the progenitor system.

In sub-Chandrasekhar-mass models, the radioactive and stable nickel are co-located. Therefore, we naively expect similar line shapes for nickel and iron. In Chandrasekhar-mass models the radioactive and stable nickel are globally much less correlated. Still, the detonation ashes, which dominate the central ejecta for multidimensional models, contain about the same ratio of stable to radioactive nickel as in sub-Chandrasekhar-mass models. The deflagration ashes of Chandrasekhar-mass models, however, which contain a much higher ratio of stable to radioactive nickel, are located at much higher velocities in the ejecta. In the most extreme model, the gravitationally confined detonation model r82_d1.0, the deflagration ashes sit at velocities larger than 10 000 km s−1. If the emission in nebular spectra of Chandrasekhar-mass models is dominated by the detonation ashes they might look very much like sub-Chandrasekhar-mass models. Detecting high-velocity stable Nickel that is fundamental to multidimensional Chandrasekhar-mass models, but seems impossible to get in sub-Chandrasekhar-mass models, could be one of the strongest possible links from observations to explosion scenario.

Crucially, this also means that the global ratio of stable to radioactive nickel in the ejecta that is currently often used to compare models and observations (Maguire et al. 2018; Dhawan et al. 2018; Flörs et al. 2020; Blondin et al. 2022; Liu et al. 2023), might be a poor indicator for nebular spectra properties of Type Ia supernovae explosion models. Type Ia supernova remnants might be different in this aspect, because all radioactive nickel and cobalt has long decayed and their luminosity comes from the reverse shock (Seitenzahl et al. 2019). Nevertheless, we expect that for Chandrasekhar-mass explosions the nickel and iron lines in remnants will predominantly originate from different regions. Therefore, it is important to take the full multidimensional density and composition structure of the ejecta into account for any analysis of late-time Type Ia supernova observables.

6. Conclusions

We have shown that we really need 3D explosion simulations and likely 3D radiative transfer simulations to obtain faithful synthetic observables of the ejecta of Type Ia supernova explosion models.

We have focused on various models under consideration to explain normal Type Ia supernovae. Of all the models we included, only the double detonation scenario with a surviving companion has inner ejecta that are anywhere close to spherical symmetry.

We have shown that the central ejecta of 3D Chandrasekhar-mass delayed detonation explosion models are detonation ashes. The deflagration ashes have velocities higher than 5000 km s−1. In contrast, 1D Chandrasekhar-mass explosion models inherently lead to an unphysical inversion of the detonation and deflagration ashes in the ejecta (Fig. 1). In 3D Chandrasekhar-mass explosion models, both radioactive and stable nickel are distributed over a range of densities that varies by up to a factor of two in the ejecta at the same velocity.

In sub-Chandrasekhar-mass merger models in which the secondary white dwarf also explodes (violent merger or helium ignited merger), the nickel-rich ashes of the primary white dwarf can be substantially compressed by the ejecta of the secondary white dwarf. The details of the progenitor and explosion mechanism imprint themselves on the timing difference between the explosion of the primary white dwarf and the secondary white dwarf (1 s for the violent merger, 4 s for the helium-ignited merger) as well as the magnitude of this compression. Since the nebular emission is highly sensitive to the range of densities that contribute to the emitting region of the ejecta, this could be used to learn about the progenitor systems of Type Ia supernovae in the future.

Recent and future JWST nebular spectra of Type Ia supernovae allow us to significantly improve our interpretation of observed Type Ia supernovae and match them with theoretical explosion models (Kwok et al. 2023, 2024; DerKacy et al. 2023; Blondin et al. 2023). However, to do so we need faithful multidimensional models that simulate the explosion and obtain synthetic observables. Until we can compute nebular spectra and remnant models fully in 3D with all necessary physics involved, we will need to develop new, better approximations than spherically symmetric averaging (e.g. Blondin et al. 2023) and show that they represent full 3D models at least reasonably well.


1

Note, however, that at low very densities deflagration flames enter the distributed burning regime and cannot be considered as thin any longer (Röpke & Hillebrandt 2005).

Acknowledgments

AJR acknowledges support from the Australian Research Council under award number FT170100243. This research was supported in part by the National Science Foundation grant PHY-1748958 to the Kavli Institute for Theoretical Physics (KITP). ST acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (LENSNOVA: grant agreement no. 771776). The work of FKR is supported by the Klaus-Tschira Foundation, by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – 537700965, and by the European Union (ERC, ExCEED, project number 1010962). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. SAS acknowledges funding from STFC grant ST/X00094X/1. This work was supported by the ‘Programme National de Physique Stellaire’ (PNPS) of CNRS/INSU cofunded by CEA and CNES.

References

  1. Arnett, W. D. 1969, Ap&SS, 5, 180 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bildsten, L., Shen, K. J., Weinberg, N. N., & Nelemans, G. 2007, ApJ, 662, L95 [Google Scholar]
  3. Blondin, S., Dessart, L., Hillier, D. J., & Khokhlov, A. M. 2013, MNRAS, 429, 2127 [NASA ADS] [CrossRef] [Google Scholar]
  4. Blondin, S., Bravo, E., Timmes, F. X., Dessart, L., & Hillier, D. J. 2022, A&A, 660, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Blondin, S., Dessart, L., Hillier, D. J., Ramsbottom, C. A., & Storey, P. J. 2023, A&A, 678, A170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Boos, S. J., Townsley, D. M., Shen, K. J., Caldwell, S., & Miles, B. J. 2021, ApJ, 919, 126 [NASA ADS] [CrossRef] [Google Scholar]
  7. Boos, S. J., Townsley, D. M., & Shen, K. J. 2024, ApJ. submitted, [arXiv:2401.08011] [Google Scholar]
  8. Bulla, M., Sim, S. A., Pakmor, R., et al. 2016, MNRAS, 455, 1060 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cameron, A. G. W. 1959, ApJ, 130, 916 [CrossRef] [Google Scholar]
  10. Churazov, E., Sunyaev, R., Isern, J., et al. 2015, ApJ, 812, 62 [NASA ADS] [CrossRef] [Google Scholar]
  11. DerKacy, J. M., Ashall, C., Hoeflich, P., et al. 2023, ApJ, 945, L2 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dhawan, S., Flörs, A., Leibundgut, B., et al. 2018, A&A, 619, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Diehl, R., Siegert, T., Hillebrandt, W., et al. 2015, A&A, 574, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Dutta, A., Sahu, D. K., Anupama, G. C., et al. 2022, ApJ, 925, 217 [NASA ADS] [CrossRef] [Google Scholar]
  15. Fink, M., Hillebrandt, W., & Röpke, F. K. 2007, A&A, 476, 1133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Fink, M., Röpke, F. K., Hillebrandt, W., et al. 2010, A&A, 514, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Flörs, A., Spyromilio, J., Taubenberger, S., et al. 2020, MNRAS, 491, 2902 [NASA ADS] [Google Scholar]
  18. Gasques, L. R., Afanasjev, A. V., Aguilera, E. F., et al. 2005, Phys. Rev. C, 72, 025806 [NASA ADS] [CrossRef] [Google Scholar]
  19. Gerardy, C. L., Meikle, W. P. S., Kotak, R., et al. 2007, ApJ, 661, 995 [NASA ADS] [CrossRef] [Google Scholar]
  20. Gronow, S., Collins, C. E., Sim, S. A., & Röpke, F. K. 2021, A&A, 649, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Guillochon, J., Dan, M., Ramirez-Ruiz, E., & Rosswog, S. 2010, ApJ, 709, L64 [Google Scholar]
  22. Höflich, P., & Stein, J. 2002, ApJ, 568, 779 [CrossRef] [Google Scholar]
  23. Hristov, B., Hoeflich, P., & Collins, D. C. 2021, ApJ, 923, 210 [NASA ADS] [CrossRef] [Google Scholar]
  24. Kasen, D., Röpke, F. K., & Woosley, S. E. 2009, Nature, 460, 869 [Google Scholar]
  25. Kawabata, M., Maeda, K., Yamanaka, M., et al. 2021, PASJ, 73, 1295 [NASA ADS] [CrossRef] [Google Scholar]
  26. Khokhlov, A. M. 1991, A&A, 245, 114 [NASA ADS] [Google Scholar]
  27. Kobayashi, C., Leung, S.-C., & Nomoto, K. 2020, ApJ, 895, 138 [CrossRef] [Google Scholar]
  28. Kromer, M., Fink, M., Stanishev, V., et al. 2013, MNRAS, 429, 2287 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kwok, L. A., Jha, S. W., Temim, T., et al. 2023, ApJ, 944, L3 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kwok, L. A., Siebert, M. R., Johansson, J., et al. 2024, ApJ, 966, 135 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lach, F., Callan, F. P., Bubeck, D., et al. 2022a, A&A, 658, A179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Lach, F., Callan, F. P., Sim, S. A., & Röpke, F. K. 2022b, A&A, 659, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Lesaffre, P., Han, Z., Tout, C. A., Podsiadlowski, P., & Martin, R. G. 2006, MNRAS, 368, 187 [NASA ADS] [CrossRef] [Google Scholar]
  34. Leung, S.-C., & Nomoto, K. 2018, ApJ, 861, 143 [Google Scholar]
  35. Li, W., Chornock, R., Leaman, J., et al. 2011, MNRAS, 412, 1473 [NASA ADS] [CrossRef] [Google Scholar]
  36. Liu, Z.-W., Röpke, F. K., & Han, Z. 2023, Res. Astron. Astrophys., 23, 082001 [CrossRef] [Google Scholar]
  37. Maeda, K., Benetti, S., Stritzinger, M., et al. 2010, Nature, 466, 82 [NASA ADS] [CrossRef] [Google Scholar]
  38. Magee, M. R., Kotak, R., Sim, S. A., et al. 2016, A&A, 589, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Maguire, K., Sim, S. A., Shingles, L., et al. 2018, MNRAS, 477, 3567 [NASA ADS] [CrossRef] [Google Scholar]
  40. Malone, C. M., Nonaka, A., Woosley, S. E., et al. 2014, ApJ, 782, 11 [NASA ADS] [CrossRef] [Google Scholar]
  41. Nonaka, A., Aspden, A. J., Zingale, M., et al. 2012, ApJ, 745, 73 [NASA ADS] [CrossRef] [Google Scholar]
  42. Pakmor, R., Kromer, M., Röpke, F. K., et al. 2010, Nature, 463, 61 [Google Scholar]
  43. Pakmor, R., Hachinger, S., Röpke, F. K., & Hillebrandt, W., 2011, A&A, 528, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Pakmor, R., Edelmann, P., Röpke, F. K., & Hillebrandt, W. 2012a, MNRAS, 424, 2222 [Google Scholar]
  45. Pakmor, R., Kromer, M., Taubenberger, S., et al. 2012b, ApJ, 747, L10 [NASA ADS] [CrossRef] [Google Scholar]
  46. Pakmor, R., Kromer, M., Taubenberger, S., & Springel, V. 2013, ApJ, 770, L8 [Google Scholar]
  47. Pakmor, R., Zenati, Y., Perets, H. B., & Toonen, S. 2021, MNRAS, 503, 4734 [Google Scholar]
  48. Pakmor, R., Callan, F. P., Collins, C. E., et al. 2022, MNRAS, 517, 5260 [NASA ADS] [CrossRef] [Google Scholar]
  49. Piersanti, L., Bravo, E., Straniero, O., Cristallo, S., & Domínguez, I. 2022, ApJ, 926, 103 [NASA ADS] [CrossRef] [Google Scholar]
  50. Plewa, T., Calder, A. C., & Lamb, D. Q. 2004, ApJ, 612, L37 [Google Scholar]
  51. Röpke, F. K. 2017, in Handbook of Supernovae, eds. A. W. Alsabti, & P. Murdin, 1185 [CrossRef] [Google Scholar]
  52. Röpke, F. K., & Hillebrandt, W. 2005, A&A, 429, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Röpke, F. K., Gieseler, M., Reinecke, M., Travaglio, C., & Hillebrandt, W. 2006, A&A, 453, 203 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Röpke, F. K., Hillebrandt, W., Schmidt, W., et al. 2007, ApJ, 668, 1132 [NASA ADS] [CrossRef] [Google Scholar]
  55. Röpke, F. K., Kromer, M., Seitenzahl, I. R., et al. 2012, ApJ, 750, L19 [CrossRef] [Google Scholar]
  56. Ruiter, A. J. 2020, IAU Symp., 357, 1 [Google Scholar]
  57. Ruiz-Lapuente, P., & Lucy, L. B. 1992, ApJ, 400, 127 [NASA ADS] [CrossRef] [Google Scholar]
  58. Seitenzahl, I. R., & Townsley, D. M. 2017, in Handbook of Supernovae, eds. A. W. Alsabti, & P. Murdin, 1955 [CrossRef] [Google Scholar]
  59. Seitenzahl, I. R., Townsley, D. M., Peng, F., & Truran, J. W. 2009, At. Data Nucl. Data Tables, 95, 96 [NASA ADS] [CrossRef] [Google Scholar]
  60. Seitenzahl, I. R., Ciaraldi-Schoolmann, F., Röpke, F. K., et al. 2013, MNRAS, 429, 1156 [NASA ADS] [CrossRef] [Google Scholar]
  61. Seitenzahl, I. R., Ghavamian, P., Laming, J. M., & Vogt, F. P. A. 2019, Phys. Rev. Lett., 123, 041101 [CrossRef] [Google Scholar]
  62. Shen, K. J., & Moore, K. 2014, ApJ, 797, 46 [Google Scholar]
  63. Shen, K. J., Kasen, D., Miles, B. J., & Townsley, D. M. 2018, ApJ, 854, 52 [Google Scholar]
  64. Shen, K. J., Blondin, S., Kasen, D., et al. 2021a, ApJ, 909, L18 [NASA ADS] [CrossRef] [Google Scholar]
  65. Shen, K. J., Boos, S. J., Townsley, D. M., & Kasen, D. 2021b, ApJ, 922, 68 [NASA ADS] [CrossRef] [Google Scholar]
  66. Siebert, M. R., Kwok, L. A., Johansson, J., et al. 2024, ApJ, 960, 88 [NASA ADS] [CrossRef] [Google Scholar]
  67. Sim, S. A., Röpke, F. K., Hillebrandt, W., et al. 2010, ApJ, 714, L52 [Google Scholar]
  68. Singh, M., Sahu, D. K., Barna, B., et al. 2024, ApJ, 965, 73 [NASA ADS] [CrossRef] [Google Scholar]
  69. Taubenberger, S. 2017, in Handbook of Supernovae, eds. A. W. Alsabti, & P. Murdin, 317 [CrossRef] [Google Scholar]
  70. Woosley, S. E., Wunsch, S., & Kuhlen, M. 2004, ApJ, 607, 921 [NASA ADS] [CrossRef] [Google Scholar]
  71. Zingale, M., Nonaka, A., Almgren, A. S., et al. 2011, ApJ, 740, 8 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1.

Distribution of radioactive 56Ni (top row) and stable 58Ni (bottom row) in density–velocity space 100 s after the explosion. Each histogram is normalised individually to the total mass of the type of nickel it shows. The columns show different Chandrasekhar-mass explosion models: a 3D delayed detonation (N100; Seitenzahl et al. 2013), a 3D gravitationally confined detonation (r82.d1.0; Lach et al. 2022b), and a 1D delayed detonation (DDC10; Blondin et al. 2013). The solid line in the left and middle panels shows the spherically averaged density profile of the total mass of the ejecta. The dashed lines indicate if nickel was synthesised in the deflagration or in the detonation. In the 3D models the deflagration ashes are located at higher velocities, i.e. further out in the ejecta, than the detonation ashes. In the 1D model this hierarchy is inverted. In the 3D models nickel is distributed over a wide range of densities at a fixed velocity, so spherical averaging 3D ejecta significantly changes their physical properties.

In the text
thumbnail Fig. 2.

Distribution of radioactive 56Ni (top row) and stable 58Ni (bottom row) in density–velocity space 100 s after the explosion. Each histogram is normalised individually to the total mass of the type of nickel it shows. The columns show different sub-Chandrasekhar-mass merger explosion models: a violent merger (1109; Pakmor et al. 2012a) and a helium ignited merger in which only the primary white dwarf (oneExp; Pakmor et al. 2022) explodes or both white dwarfs explode (twoExp; Pakmor et al. 2022). The solid line shows the spherically averaged density profile of the total mass of the ejecta. The distributions of radioactive and stable nickel in the ejecta are highly correlated. Only the oneExp model is close to spherical symmetry. In both the 1109 and twoExp models the secondary white dwarf explodes as well, compressing the ejecta of the primary white dwarf that contain primarily nickel. In the 1109 model the secondary explosion compresses the nickel to much higher densities than in the twoExp models because it happens with a delay of only 1 s compared to 4 s in the twoExp model.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.