-
Observational signatures of thermonuclear electron-capture supernovae -- Ne II line strengthening and color evolution as traces of the explosion mechanism
Authors:
Alexander Holas,
Fionntan P. Callan,
Samuel W. Jones,
Friedrich K. Roepke,
Rüdiger Pakmor,
Alexandra Kozyreva,
Christine E. Collins,
Luke J. Shingles,
Stuart A. Sim,
Joshua M. Pollin
Abstract:
Thermonuclear electron-capture supernovae (tECSNe) are a potential fate of certain intermediate mass stars forming ONe cores at the end of their evolution. Simulations suggest that these explosions are a viable alternative to collapse, yet no synthetic observables exist that allow for their identification among observed transients. We present first of their kind synthetic observables of a tECSN si…
▽ More
Thermonuclear electron-capture supernovae (tECSNe) are a potential fate of certain intermediate mass stars forming ONe cores at the end of their evolution. Simulations suggest that these explosions are a viable alternative to collapse, yet no synthetic observables exist that allow for their identification among observed transients. We present first of their kind synthetic observables of a tECSN simulation, aiming to establish whether these explosions can occur in nature, and investigate potential observational signatures to separate them from similar transients such as pure deflagrations in CO white dwarfs. We carry out 3D photospheric phase and 1D late phase simulations using the radiative transfer code Artis. As input, we use a tECSN explosion simulation and a CO deflagration simulation with comparable $^{56}$Ni production, both computed with the Leafs code. Both models have similar observational characteristics, akin to SNe~Iax-like events. The tECSN ejecta model are characterized by a $M(^{56}\mathrm{Ni})/M_\mathrm{ej}$ ratio $25\%$ lower than that of comparable CO deflagration models. At early times, the tECSN model shows a slower decline in the red colors compared to the CO deflagration due the greater amount of Ti and Cr synthesized in the tECSN explosion. At late times, the tECSN model exhibits an exceptionally strong $12.8\,μ$m Ne II emission line, that strengthens substantially over time, whereas its strength remains largely unchanged in the CO deflagration. Our results suggest tECSNe could potentially result in SN~Iax-like transients. Importantly, we find no features that are in tension with existing observables. So far, there are no indicators that unambiguously and robustly separate tECSNe from deflagrations in CO white dwarfs. Nonetheless, our work highlights the potential importance of the mid-infrared wavelength range for distinguishing possible explosion mechanisms.
△ Less
Submitted 3 June, 2026;
originally announced June 2026.
-
Angular momentum drives proton-rich nucleosynthesis in hyperaccreting neutron stars in common envelopes
Authors:
Alexander D. Hall-Smith,
Sophie E. D. Abrahams,
Alison M. Laird,
Christian Aa. Diget,
Christopher Fryer,
Samuel W. Jones
Abstract:
Interacting binaries can produce a wide range of exotic systems, including X-ray binaries and merging neutron stars, through a mass transfer phase called Common Envelope (CE) evolution. A CE phase can occur during rapid expansion as a star as it moves off the main sequence. If the engulfed star is a compact object (e.g. neutron star), a CE phase can lead to hyperaccretion onto the neutron star. Pr…
▽ More
Interacting binaries can produce a wide range of exotic systems, including X-ray binaries and merging neutron stars, through a mass transfer phase called Common Envelope (CE) evolution. A CE phase can occur during rapid expansion as a star as it moves off the main sequence. If the engulfed star is a compact object (e.g. neutron star), a CE phase can lead to hyperaccretion onto the neutron star. Previous work focused on systems in which the accreting material has low angular momentum, studying turbulent outflows. This study investigates the impact of angular momentum on accreting material leading to the formation of an accretion disk. Disk accretion systems lead to very different nuclear burning conditions.
This paper presents the results of nucleosynthesis modelling of material ejected from an accretion disk surrounding a 1.5 M$_{\odot}$ neutron star in a CE with a 15 M$_{\odot}$ companion. As material is accreted towards the neutron star, sufficient heating will occur to eject a fracton of the material back into the surrounding envelope, producing a nucleosynthetic yield signature that differs from other explosions. We find that significant mass fractions of rp-process products are synthesised, thereby providing another mechanism for rp-process contribution to galactic chemical evolution, following ejection of the CE. Furthermore, later stages of the CE evolution the accrete helium leading to alpha-rich, supernova-like nucleosynthesis, producing $^{44}$Ti and $^{56}$Ni. Further work on modelling both the accretion disk wind, and the companion envelope ejection, is vital to understand the contributions of these scenarios to chemical evolution.
△ Less
Submitted 6 March, 2026;
originally announced March 2026.
-
Localized $^{18}$O production in white dwarf mergers
Authors:
Alexander Holas,
Veronica Agaeva,
Friedrich K. Roepke,
Samuel W. Jones,
Javier Moran-Fraile,
Marco Vetter,
Rüdiger Pakmor,
Philipp Podsiadlowski
Abstract:
The merger of a He white dwarf (WD) and a CO WD is the favored formation channel for R Coronae Borealis (RCB) stars. These stars exhibit ${^{16}}\mathrm{O}/{^{18}}\mathrm{O}$ ratios that are orders of magnitude lower than the solar value. However, it is not fully understood whether such low ${^{16}}\mathrm{O}/{^{18}}\mathrm{O}$ ratios can be achieved in WD merger remnants for the predicted lifetim…
▽ More
The merger of a He white dwarf (WD) and a CO WD is the favored formation channel for R Coronae Borealis (RCB) stars. These stars exhibit ${^{16}}\mathrm{O}/{^{18}}\mathrm{O}$ ratios that are orders of magnitude lower than the solar value. However, it is not fully understood whether such low ${^{16}}\mathrm{O}/{^{18}}\mathrm{O}$ ratios can be achieved in WD merger remnants for the predicted lifetime of RCB stars of around $10^4\,\mathrm{years}$. In this work, we perform detailed nucleosynthesis calculations of a 3D magnetohydrodynamical simulation of a merger of a $0.3\,M_\odot$ He WD and a $0.6\,M_\odot$ CO WD for $4000\,\mathrm{s}$ at which point a steady state in temperature and density is reached. From this point, we follow several radial zones to study the long-term production of ${^{18}}\mathrm{O}$ and its variability throughout the burning region.
We find that the asymmetric merger process leaves an imprint on the distribution of the abundances at the end of our hydrodynamic simulation. During the long-term evolution up to $100\,\mathrm{years}$, we observe ${^{16}}\mathrm{O}/{^{18}}\mathrm{O}$ ratios of order of unity, although the timescale on which ${^{18}}\mathrm{O}$ is destroyed again is highly location dependent. Importantly, our calculations suggest that in the outer layers of the burning shell, the dominant production channel is $^{14}\mathrm{C}(α,γ)^{18}\mathrm{O}$ instead of the commonly considered $^{14}\mathrm{N}(α,γ)^{18}\mathrm{F}(β^+)^{18}\mathrm{O}$ reaction, whereby the former can be sustained for longer periods of time. Furthermore, these outer regions do not reach the conditions necessary for fast $α$-captures in ${^{18}}\mathrm{O}$ to ${^{22}}\mathrm{Ne}$, thus being favorable to maintaining a low ${^{16}}\mathrm{O}/{^{18}}\mathrm{O}$ ratio.
△ Less
Submitted 18 February, 2026;
originally announced February 2026.
-
Drawing the line between explosion and collapse in electron-capture supernovae -- I. Impact of conductive flame speeds and ignition conditions on the explosion mechanism
Authors:
Alexander Holas,
Samuel W. Jones,
Friedrich K. Roepke,
Rüdiger Pakmor,
Christina Fakiola,
Giovanni Leidi,
Raphael Hirschi,
Ken J. Shen
Abstract:
Electron-capture supernovae (ECSNe) are commonly thought to result in a collapse to a neutron star. Recent work has shown that a thermonuclear explosion is also a possible outcome. The division between the two regimes has not yet been mapped out. In this study, we investigate the conditions under which the transition from thermonuclear explosion to collapse occurs, and what physical mechanisms dri…
▽ More
Electron-capture supernovae (ECSNe) are commonly thought to result in a collapse to a neutron star. Recent work has shown that a thermonuclear explosion is also a possible outcome. The division between the two regimes has not yet been mapped out. In this study, we investigate the conditions under which the transition from thermonuclear explosion to collapse occurs, and what physical mechanisms drive each outcome. We conducted a parameter study of 56 3D hydrodynamic simulations of ECSN in ONe white dwarfs using a level set based flame model implemented in the Leafs code. We varied both the ignition location and the central density at ignition to determine the conditions of the transition regime. Additionally, we explored two different laminar flame parameterizations and how they impact the simulation outcome. From our parameter study, we find a transition density in the range of $\logρ_c^{ini}=10.0$ and $10.15$ g cm$^{-3}$, depending on the ignition location and utilized laminar flame speed parameterization. Importantly, we find that for sufficiently high central densities, the burned ashes can sink into the core and trap large amounts of neutron-rich material in the bound remnant. In the transition regime between explosion and collapse, we find that the laminar flame speed plays a critical role by suppressing the formation of instabilities and thereby reducing the nuclear energy generation needed to overcome the collapse. We find that a thermonuclear explosion is possible for a wide range of parameters, whereby a more off-center ignition allows for higher central densities to still result in an explosion. Both the conditions at ignition and the flame physics are critical in determining the outcome. Detailed 3D hydrodynamic simulations of the preceding stellar evolution and the ignition process of the thermonuclear flame are necessary to accurately predict the outcome of ECSNe.
△ Less
Submitted 23 January, 2026;
originally announced January 2026.
-
A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics
Authors:
Samuel W. Jones,
Colin P. McNally,
Meritt Reynolds
Abstract:
Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on…
▽ More
Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. We present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.
△ Less
Submitted 7 January, 2026;
originally announced January 2026.
-
First models of s process in AGB stars of solar metallicity for the stellar evolutionary code ATON with a novel stable explicit numerical solver
Authors:
A. Yagüe López,
D. A. García-Hernández,
P. Ventura,
C. L. Doherty,
J. W. den Hartogh,
S. W. Jones,
M. Lugaro
Abstract:
Aims. We describe the first s-process post-processing models for asymptotic giant branch (AGB) stars of masses 3, 4 and 5 M at solar metallicity (Z=0.018) computed using the input from the stellar evolutionary code aton. Methods. The models are computed with the new code snuppat(S-process NUcleosynthesis Post-Processing code for aton), including an advective scheme for the convective overshoot tha…
▽ More
Aims. We describe the first s-process post-processing models for asymptotic giant branch (AGB) stars of masses 3, 4 and 5 M at solar metallicity (Z=0.018) computed using the input from the stellar evolutionary code aton. Methods. The models are computed with the new code snuppat(S-process NUcleosynthesis Post-Processing code for aton), including an advective scheme for the convective overshoot that leads to the formation of the main neutron source, 13C. Each model is post-processed with 3 different values of the free overshoot parameter. Included in the code snuppat is the novel Patankar-Euler-Deflhard explicit numerical solver, that we use to solve the nuclear network system of differential equations. Results. The results are compared to those from other s-process nucleosynthesis codes (Monash,fruity, and NuGrid), as well as observations of s-process enhancement in AGB stars, planetary nebulae, and barium stars. This comparison shows that the relatively high abundance of12C in the He-rich intershell in aton results in as-process abundance pattern that favours the second over the first s-process peak for all the masses explored. Also, our choice of an advective as opposed to diffusive numerical scheme for the convective overshoot results in significants-process nucleosynthesis also for the 5 M models, which may be in contradiction with observations.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
The Nucleosynthetic Yields of Core-Collapse Supernovae, prospects for the Next Generation of Gamma-Ray Astronomy
Authors:
Sydney Andrews,
Chris L. Fryer,
Samuel W. Jones,
Wesley P. Even,
Marco Pignatari
Abstract:
Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key…
▽ More
Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants.
△ Less
Submitted 22 December, 2019;
originally announced December 2019.
-
Nucleosynthetic Yields from Neutron Stars Accreting in Binary Common Envelopes
Authors:
James D. Keegans,
Chris L. Fryer,
Samuel W. Jones,
Benoit Cote,
Krzysztof Belczynski,
Falk Herwig,
Marco Pignatari,
Alison M. Laird,
Christian Aa. Diget
Abstract:
Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario f…
▽ More
Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.
△ Less
Submitted 5 February, 2019;
originally announced February 2019.