Pre-explosion Properties of Helium Star Donors to Thermonuclear Supernovae

Helium star–carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code MESA, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to near-UV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multiday orbit. However, extending our binary models to initial CO WD masses of up to 1.2 M ⊙, we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star–WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition.

Predicted spatial and velocity distributions of ejected companion stars of helium accretion-induced thermonuclear supernovae

Context. Thermonuclear supernovae (SNe), a subset of which are the highly important SNe Type Ia, remain one of the more poorly understood phenomena known to modern astrophysics. In recent years, the single degenerate helium (He) donor channel, where a white dwarf star accretes He-rich matter from a hydrogen-depleted companion, has emerged as a promising candidate progenitor scenario for these events. An unresolved question in this scenario is the fate of the companion star, which would be evident as a runaway hot subdwarf O/B stars (He sdO/B) in the aftermath of the SN event. Aims. Previous studies have shown that the kinematic properties of an ejected companion provide an opportunity to closer examine the properties of an SN progenitor system. However, with the number of observed objects not matching predictions by theory, the viability of this mechanism is called into question. In this study, we first synthesize a population of companion stars ejected by the aforementioned mechanism, taking into account predicted ejection velocities, the inferred population density in the Galactic mass distribution, and subsequent kinematics in the Galactic potential. We then discuss the astrometric properties of this population. Methods. We present 106 individual ejection trajectories, which were numerically computed with a newly developed, lightweight simulation framework. Initial conditions were randomly generated, but weighted according to the Galactic mass density and ejection velocity data. We then discuss the bulk properties (Galactic distribution and observational parameters) of our sample. Results. Our synthetic population reflects the Galactic mass distribution. A peak in the density distribution for close objects is expected in the direction of the Galactic centre. Higher mass runaways should outnumber lower mass ones. If the entire considered mass range is realised, the radial velocity distribution should show a peak at 500 km s−1. If only close US 708 analogues are considered, there should be a peak at (∼750 − 850) km s−1. In either case, US 708 should be a member of the high-velocity tail of the distribution. Conclusions. We show that the puzzling lack of confirmed surviving companion stars of thermonuclear SNe, though possibly an observation-related selection effect, may indicate a selection against high mass donors in the SD He donor channel.

An Elementary Model for a Self-Accelerating Outward Propagating Flame Subject to the Rayleigh–Taylor Instability: Transition to Detonation

Within the Boussinesq approximation, an elementary model for the deflagration-to-detonation transition triggered by self-acceleration of an expanding flame is formulated and explored. The self-acceleration is sustained by the intrinsic Rayleigh–Taylor instability until the Deshaies–Joulin deflagrability threshold is reached, followed by an abrupt transition to detonation. Emergence of the threshold is caused by positive feedback between the accelerating flame and the flame-driven pressure shock that results in the thermal runaway when the flame speed reaches a critical level. The model offers a simple mechanism that may be responsible for the transition to detonation in thermonuclear supernovae.

Heavy elements nucleosynthesis on accreting white dwarfs: building seeds for the p-process

ABSTRACT The origin of the proton-rich trans-iron isotopes in the Solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor’s rapidly accreting white dwarf (RAWD) phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf (WD) models with initial masses ≥ 0.85 $\, \mathrm{M}_\odot$ using the stellar code mesa The near-surface layers of the 1, 1.26, 1.32 and 1.38 $\, \mathrm{M}_\odot$ models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous RAWD models at lower mass, we find that the H-shell flashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 1015 cm−3 in the most massive WDs. Through the recurrence of the H-shell flashes these intermediate neutron densities can be sustained effectively for a long time leading to high-neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to supersolar abundances for p-nuclei with A > 96. Our models show that SNIa are a viable p-process production site.

SDSS J124043.01+671034.68: the partially burned remnant of a low-mass white dwarf that underwent thermonuclear ignition?

ABSTRACT The white dwarf SDSS J124043.01+671034.68 (SDSS J1240+6710) was previously found to have an oxygen-dominated atmosphere with significant traces of neon, magnesium, and silicon. A possible origin via a violent late thermal pulse or binary interactions has been suggested to explain this very unusual photospheric composition. We report the additional detection of carbon, sodium, and aluminium in far-ultraviolet and optical follow-up spectroscopy. No iron-group elements are detected, with tight upper limits on titanium, iron, cobalt, and nickel, suggesting that the star underwent partial oxygen burning, but failed to ignite silicon burning. Modelling the spectral energy distribution and adopting the distance based on the Gaia parallax, we infer a low white dwarf mass, $M_\mathrm{wd}=0.41\pm 0.05\, \mathrm{M}_\odot$. The large space velocity of SDSS J1240+6710, computed from the Gaia proper motion and its radial velocity, is compatible with a Galactic rest-frame velocity of ≃ 250 km s−1 in the opposite direction with respect to the Galactic rotation, strongly supporting a binary origin of this star. We discuss the properties of SDSS J1240+6710 in the context of the recently identified survivors of thermonuclear supernovae, the D6 and LP 40−365 stars, and conclude that it is unlikely related to either of those two groups. We tentatively suggest that SDSS J1240+6710 is the partially burned remnant of a low-mass white dwarf that underwent a thermonuclear event.

44Ti ejecta in young supernova remnants

Context. Tracing unstable isotopes produced in supernova nucleosynthesis provides a direct diagnostic of supernova explosion physics. Theoretical models predict an extensive variety of scenarios, which can be constrained through observations of the abundant isotopes 56Ni and 44Ti. Direct evidence of the latter was previously found only in two core-collapse supernova events, and appears to be absent in thermonuclear supernovae. Aims. We aim to to constrain the supernova progenitor types of Cassiopeia A, SN 1987A, Vela Jr., G1.9+0.3, SN1572, and SN1604 through their 44Ti ejecta masses and explosion kinematics. Methods. We analyzed INTEGRAL/SPI observations of the candidate sources utilizing an empirically motivated high-precision background model. We analyzed the three dominant spectroscopically resolved de-excitation lines at 68, 78, and 1157 keV emitted in the decay chain of 44Ti→44Sc→44Ca. The fluxes allow the determination of the production yields of 44Ti. Remnant kinematics were obtained from the Doppler characteristics of the lines. Results. We find a significant signal for Cassiopeia A in all three lines with a combined significance of 5.4σ. The fluxes are (3.3 ± 0.9) × 10−5 ph cm−2 s−1, and (4.2 ± 1.0) × 10−5 ph cm−2 s−1 for the 44Ti and 44Sc decay, respectively. This corresponds to a mass of (2.4 ± 0.7) × 10−4 M⊙ and (3.1 ± 0.8) × 10−4 M⊙, respectively. We obtain higher fluxes for 44Ti with our analysis of Cassiopeia A than were obtained in previous analyses. We discuss potential differences. We interpret the line width from Doppler broadening as expansion velocity of (6400 ± 1900) km s−1. We do not find any significant signal for any other candidate sources. Conclusions. We obtain a high 44Ti ejecta mass for Cassiopeia A that is in disagreement with ejecta yields from symmetric 2D models. Upper limits for the other core-collapse supernovae are in agreement with model predictions and previous studies. The upper limits we find for the three thermonuclear supernovae (G1.9+0.3, SN1572 and SN1604) consistently exclude the double detonation and pure helium deflagration models as progenitors.

Determining the 56Ni distribution of type Ia supernovae from observations within days of explosion

Recent studies have shown how the distribution of 56Ni within the ejected material of type Ia supernovae can have profound consequences on the observed light curves. Observations at early times can therefore provide important details on the explosion physics in thermonuclear supernovae, which are poorly constrained. To this end, we present a series of radiative transfer calculations that explore variations in the 56Ni distribution. Our models also show the importance of the density profile in shaping the light curve, which is often neglected in the literature. Using our model set, we investigate the observations that are necessary to determine the 56Ni distribution as robustly as possible within the current model set. We find that this includes observations beginning at least ∼14 days before B-band maximum, extending to approximately maximum light with a relatively high (≲3 day) cadence, and in at least one blue and one red band (such as B and R, or g and r) are required. We compare a number of well-observed type Ia supernovae that meet these criteria to our models and find that the light curves of ∼70–80% of objects in our sample are consistent with being produced solely by variations in the 56Ni distributions. The remaining supernovae show an excess of flux at early times, indicating missing physics that is not accounted for within our model set, such as an interaction or the presence of short-lived radioactive isotopes. Comparing our model light curves and spectra to observations and delayed detonation models demonstrates that while a somewhat extended 56Ni distribution is necessary to reproduce the observed light curve shape, this does not negatively affect the spectra at maximum light. Investigating current explosion models shows that observations typically require a shallower decrease in the 56Ni mass towards the outer ejecta than is produced for models of a given 56Ni mass. Future models that test differences in the explosion physics and detonation criteria should be explored to determine the conditions necessary to reproduce the 56Ni distributions found here.