Producing type Iax supernovae from a specific class of helium-ignited WD explosions?

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. However, their progenitors are still uncertain. We study whether the population properties of this class might be understood if the events originate from a subset of sub-Chandrasekhar mass explosions. In this potential progenitor population, a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star; ignition of that helium layer then leads to ignition of the CO WD. We incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain rates and delay times for such events. The predicted Galactic event rate of these explosions is ~1.5\times10^{-3}{yr}^{-1} according to our standard model, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are ~70Myr-800Myr, consistent with the fact that most of SNe Iax have been discovered in late-type galaxies. If the explosions are assumed to be double-detonations -- following current model expectations -- then based on the CO WD masses at explosion we also estimate the distribution of resulting SN brightness (-13 \gtrsim M_{bol} \gtrsim -19mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not. We suggest that the high mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might be necessary to produce SN 2002cx-like characteristics, perhaps even by changing the nature of the CO ignition.Comment: 8 pages, 10 figures, 1 table, accepted for publication in Astronomy and Astrophysic

On the role of recombination in common-envelope ejections

The energy budget in common-envelope events (CEEs) is not well understood, with substantial uncertainty even over to what extent the recombination energy stored in ionised hydrogen and helium might be used to help envelope ejection. We investigate the reaction of a red-giant envelope to heating which mimics limiting cases of energy input provided by the orbital decay of a binary during a CEE, specifically during the post-plunge-in phase during which the spiral-in has been argued to occur on a time-scale longer than dynamical. We show that the outcome of such a CEE depends less on the total amount of energy by which the envelope is heated than on how rapidly the energy was transferred to the envelope and on where the envelope was heated. The envelope always becomes dynamically unstable before receiving net heat energy equal to the envelope's initial binding energy. We find two types of outcome, both of which likely lead to at least partial envelope ejection: "runaway" solutions in which the expansion of the radius becomes undeniably dynamical, and superficially "self-regulated" solutions, in which the expansion of the stellar radius stops but a significant fraction of the envelope becomes formally dynamically unstable. Almost the entire reservoir of initial helium recombination energy is used for envelope expansion. Hydrogen recombination is less energetically useful, but is nonetheless important for the development of the dynamical instabilities. However, this result requires the companion to have already plunged deep into the envelope; therefore this release of recombination energy does not help to explain wide post-common-envelope orbits.Comment: 17 pages, 10 figures, submitted to MNRAS. Comments are welcom

Globular cluster formation efficiencies from black-hole X-ray binary feedback

We investigate a scenario in which feedback from black-hole X-ray binaries (BHXBs) sometimes begins inside young star clusters before strong supernova feedback. Those BHXBs could reduce the gas fraction inside embedded young clusters whilst maintaining virial equilibrium, which may help globular clusters (GCs) to stay bound when supernova-driven gas ejection subsequently occurs. Adopting a simple toy model with parameters guided by BHXB population models, we produce GC formation efficiencies consistent with empirically-inferred values. The metallicity dependence of BHXB formation could naturally explain why GC formation efficiency is higher at lower metallicity. For reasonable assumptions about that metallicity dependence, our toy model can produce a GC metallicity bimodality in some galaxies without a bimodality in the field-star metallicity distribution.Comment: Accepted to ApJ Letters on 19th July. 6 pages. The definitive version is available from: http://iopscience.iop.org/2041-8205/809/1/L16

Identification of the Long-Sought Common-Envelope Events

Common-envelope events (CEEs), during which two stars temporarily orbit within a shared envelope, are believed to be vital for the formation of a wide range of close binaries. For decades, the only evidence that CEEs actually occur has been indirect, based on the existence of systems that could not be otherwise explained. Here we propose a direct observational signature of CEE arising from a physical model where emission from matter ejected in a CEE is controlled by a recombination front as the matter cools. The natural range of timescales and energies from this model, as well the expected colors, light-curve shapes, ejection velocities and event rate, match those of a recently-recognized class of red transient outbursts.Comment: 6 main and 22 supplemental pages, 5 total figures, one table and 2 movies. This is the authors version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science, vol 339, 2013, http://www.sciencemag.org/content/339/6118/433.abstrac

Sub-Chandrasekhar White Dwarf Mergers as the Progenitors of Type Ia Supernovae

Type Ia supernovae (SNe Ia) are generally thought to be due to the thermonuclear explosions of carbon–oxygen white dwarfs (COWDs) with masses near the Chandrasekhar mass. This scenario, however, has two long-standing problems. First, the explosions do not naturally produce the correct mix of elements, but have to be finely tuned to proceed from subsonic deflagration to supersonic detonation. Second, population models and observations give formation rates of near-Chandrasekhar WDs that are far too small. Here, we suggest that SNe Ia instead result from mergers of roughly equal-mass CO WDs, including those that produce sub-Chandrasekhar mass remnants. Numerical studies of such mergers have shown that the remnants consist of rapidly rotating cores that contain most of the mass and are hottest in the center, surrounded by dense, small disks. We argue that the disks accrete quickly, and that the resulting compressional heating likely leads to central carbon ignition. This ignition occurs at densities for which pure detonations lead to events similar to SNe Ia. With this merger scenario, we can understand the type Ia rates and have plausible reasons for the observed range in luminosity and for the bias of more luminous supernovae toward younger populations. We speculate that explosions of WDs slowly brought to the Chandrasekhar limit—which should also occur—are responsible for some of the “atypical” SNe Ia

Aspherical supernova explosions and formation of compact black hole low-mass X-ray binaries

It has been suggested that black-hole low-mass X-ray binaries (BHLMXBs) with short orbital periods may have evolved from BH binaries with an intermediate-mass secondary, but the donor star seems to always have higher effective temperatures than measured in BHLMXBs (Justham, Rappaport & Podsiadlowski 2006). Here we suggest that the secondary star is originally an intermediate-mass (\sim 2-5 M_{\sun}) star, which loses a large fraction of its mass due to the ejecta impact during the aspherical SN explosion that produced the BH. The resulted secondary star could be of low-mass (\la 1 M_{\sun}). Magnetic braking would shrink the binary orbit, drive mass transfer between the donor and the BH, producing a compact BHLMXB.Comment: 4 pages, accepted for publication in MNRAS Letter

STROOPWAFEL: Simulating rare outcomes from astrophysical populations, with application to gravitational-wave sources

Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds $\sim$25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor $\sim$3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.Comment: Accepted. Data and scripts to reproduce main results is publicly available. The code for the STROOPWAFEL algorithm will be made publicly available. Early inquiries can be addressed to the lead autho