Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, 15$^{+9}_{−8}$$$^{+15}_{−14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that φ Persei is a likely progenitor and that eccentric neutron star – white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.E.Z. is supported by the Netherlands Research School for Astronomy (NOVA). S.D.M. acknowledges support by a Marie Sklodowska-Curie Action Incoming Fellowship (H2020 MSCA-IF-2014, project id 661502). R.G.I. thanks STFC for his Rutherford fellowship (ST/L003910/1), the DAAD for funding TS, and Churchill college for funding his bi-fellowship and for access to their library. C.B. acknowledges NASA ADAP grant NNX15AM03G S01 and NSF/AST-1412980. S.C.Y. was supported by the Korea Astronomy and Space Science Institute under the R&D program (Project No. 3348-20160002) supervised by the Ministry of Science, ICT and Future Planning. The authors further acknowledge the Leiden Lorentz Center workshop “The Impact of Massive Binaries Throughout the Universe” and the Munich Institute for Astro- and Particle Physics (MIAPP) of the DFG cluster of excellence “Origin and Structure of the Universe” for supporting the “Physics of Supernovae”

Predicting the Presence of Companions for Stripped-envelope Supernovae: The Case of the Broad-lined Type Ic SN 2002ap

Many young, massive stars are found in close binaries. Using population synthesis simulations we predict the likelihood of a companion star being present when these massive stars end their lives as core-collapse supernovae (SNe). We focus on stripped-envelope SNe, whose progenitors have lost their outer hydrogen and possibly helium layers before explosion. We use these results to interpret new Hubble Space Telescope observations of the site of the broad-lined Type Ic SN 2002ap, 14 years post-explosion. For a subsolar metallicity consistent with SN 2002ap, we expect a main-sequence (MS) companion present in about two thirds of all stripped-envelope SNe and a compact companion (likely a stripped helium star or a white dwarf/neutron star/black hole) in about 5% of cases. About a quarter of progenitors are single at explosion (originating from initially single stars, mergers, or disrupted systems). All of the latter scenarios require a massive progenitor, inconsistent with earlier studies of SN 2002ap. Our new, deeper upper limits exclude the presence of an MS companion star >8–10 $\,{M}_{\odot }$, ruling out about 40% of all stripped-envelope SN channels. The most likely scenario for SN 2002ap includes nonconservative binary interaction of a primary star initially $\lesssim 23\,{M}_{\odot }$. Although unlikely (<1% of the scenarios), we also discuss the possibility of an exotic reverse merger channel for broad-lined Type Ic events. Finally, we explore how our results depend on the metallicity and the model assumptions and discuss how additional searches for companions can constrain the physics that govern the evolution of SN progenitors.This work is based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. It is also based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory (JPL), California Institute of Technology, under a contract with NASA. Support was provided by NASA through grants GO-14075 and AR-14295 from STScI. We thank Andrew Dolphin for his patient advice on how best to implement artificial star tests in Dolphot. E.Z. is supported by a grant of the Netherlands Research School for Astronomy (NOVA). S.d.M. acknowledges support by a Marie Sklodowska-Curie Action (H2020 MSCA-IF-2014, project BinCosmos, ID 661502). A.V.F.'s group is also grateful for generous financial assistance from the Christopher R. Redlich Fund, the TABASGO Foundation, and NSF grant AST-1211916. N.S. is grateful for support from NSF grants AST-1312221 and AST-1515559. R.G.I. thanks the STFC for funding his Rutherford Fellowship under grant ST/L003910/1 and Churchill College, Cambridge, for his fellowship and access to their library

On the formation history of Galactic double neutron stars

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including DNS formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova natal kick distributions, remnant mass prescriptions, and common-envelope energetics.We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models.We find that mass transfer from a stripped post-heliumburning secondary (case BB) on to a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian σ = 30 km s-1) and high (σ = 265 km s-1) components is preferred over a single high-kick component. We conclude that the observed DNS mass distribution can place strong constraints on model assumptions

The Tarantula Massive Binary Monitoring. I. Observational campaign and OB-type spectroscopic binaries

ontext. Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way. Aims. The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region. Methods. In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries. Results. Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus’s core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2–O7) tend to have shorter orbital periods than later spectral types (O9.2–O9.7). Conclusions. Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z⊙) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe’s peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z ~ 0.5 Z⊙