No Arabic abstract
It has long been recognized that the neutrinos detected from the next core-collapse supernova in the Galaxy have the potential to reveal important information about the dynamics of the explosion and the nucleosynthesis conditions as well as allowing us to probe the properties of the neutrino itself. The neutrinos emitted from thermonuclear - type Ia - supernovae also possess the same potential, although these supernovae are dimmer neutrino sources. For the first time, we calculate the time, energy, line of sight, and neutrino-flavor-dependent features of the neutrino signal expected from a three-dimensional delayed-detonation explosion simulation, where a deflagration-to-detonation transition triggers the complete disruption of a near-Chandrasekhar mass carbon-oxygen white dwarf. We also calculate the neutrino flavor evolution along eight lines of sight through the simulation as a function of time and energy using an exact three-flavor transformation code. We identify a characteristic spectral peak at $sim 10$ MeV as a signature of electron captures on copper. This peak is a potentially distinguishing feature of explosion models since it reflects the nucleosynthesis conditions early in the explosion. We simulate the event rates in the Super-K, Hyper-K, JUNO, and DUNE neutrino detectors with the SNOwGLoBES event rate calculation software and also compute the IceCube signal. Hyper-K will be able to detect neutrinos from our model out to a distance of $sim 10$ kpc. At 1 kpc, JUNO, Super-K, and DUNE would register a few events while IceCube and Hyper-K would register several tens of events.
Despite their use as cosmological distance indicators and their importance in the chemical evolution of galaxies, the unequivocal identification of the progenitor systems and explosion mechanism of normal type Ia supernova (SN Ia) remains elusive. The leading hypothesis is that such a supernova is a thermonuclear explosion of a carbon-oxygen white dwarf, but the exact explosion mechanism is still a matter of debate. Observation of a galactic SN Ia would be of immense value in answering the many open questions related to these events. One potentially useful source of information about the explosion mechanism and progenitor is the neutrino signal. In this paper we compute the expected neutrino signal from a gravitationally confined detonation (GCD) explosion scenario for a SN~Ia and show how the flux at Earth contains features in time and energy unique to this scenario. We then calculate the expected event rates in the Super-K, Hyper-K, JUNO, DUNE, and IceCube detectors and find both Hyper-K and IceCube would see a few events for a GCD supernova at 1 kpc or closer, while Super-K, JUNO, and DUNE would see a events if the supernova were closer than ${sim}0.3$ kpc. The distance and detector criteria needed to resolve the time and spectral features arising from the explosion mechanism, neutrino production, and neutrino oscillation processes are also discussed. The neutrino signal from the GCD is then compared with the signal from a deflagration-to-detonation transition (DDT) explosion model computed previously. We find the overall event rate is the most discriminating feature between the two scenarios followed by the event rate time structure. Using the event rate in the Hyper-K detector alone, the DDT can be distinguished from the GCD at 2$sigma$ if the distance to the supernova is less than $2.3;{rm kpc}$ for a normal mass ordering and $3.6;{rm kpc}$ for an inverted ordering.
The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.
We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled statistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions that can, with a particular choice of transition density, produce a similar average and range of Ni-56 masses to those inferred from observations. Within this framework, we utilize a more realistic simmered white dwarf progenitor model with a flame model and energetics scheme to calculate the amount of Ni-56 and NSE material synthesized for a suite of simulated explosions in which the transition density is varied in the range 1-3x10^7 g/cc. We find a quadratic dependence of the NSE yield on the log of the transition density, which is determined by the competition between plume rise and stellar expansion. By considering the effect of metallicity on the transition density, we find the NSE yield decreases by 0.055 +/- 0.004 solar masses for a 1 solar metallicity increase evaluated about solar metallicity. For the same change in metallicity, this result translates to a 0.067 +/- 0.004 solar mass decrease in the Ni-56 yield, slightly stronger than that due to the variation in electron fraction from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat ambiguous, but the dependence we find is comparable to that inferred from some studies.
The double-detonation explosion model has been considered a candidate for explaining astrophysical transients with a wide range of luminosities. In this model, a carbon-oxygen white dwarf star explodes following detonation of a surface layer of helium. One potential signature of this explosion mechanism is the presence of unburned helium in the outer ejecta, left over from the surface helium layer. In this paper we present simple approximations to estimate the optical depths of important He I lines in the ejecta of double-detonation models. We use these approximations to compute synthetic spectra, including the He I lines, for double-detonation models obtained from hydrodynamical explosion simulations. Specifically, we focus on photospheric-phase predictions for the near-infrared 10830 AA~and 2 $mu$m lines of He I. We first consider a double detonation model with a luminosity corresponding roughly to normal SNe Ia. This model has a post-explosion unburned He mass of 0.03 $M_{odot}$ and our calculations suggest that the 2 $mu$m feature is expected to be very weak but that the 10830 AA~feature may have modest opacity in the outer ejecta. Consequently, we suggest that a moderate-to-weak He I 10830 AA~feature may be expected to form in double-detonation explosions at epochs around maximum light. However, the high velocities of unburned helium predicted by the model ($sim 19,000$~km~s$^{-1}$) mean that the He I 10830 AA~feature may be confused or blended with the C I 10690~AA~line forming at lower velocities. We also present calculations for the He I 10830 AA~and 2 $mu$m lines for a lower mass (low luminosity) double detonation model, which has a post-explosion He mass of 0.077 $M_{odot}$. In this case, both the He I features we consider are strong and can provide a clear observational signature of the double-detonation mechanism.
We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that both our ADR and energy release methods do not generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well-behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia Supernovae.