J. J. Eldridge

The death of massive stars - I. Observational constraints on the progenitors of type II-P supernovae

The progenitors of many type II core-collapse supernovae have now been identified directly on pre-discovery imaging. Here we present an extensive search for the progenitors of type Ibc supernovae in all available pre-discovery imaging since 1998. There are 12 type Ibc supernovae with no detections of progenitors in either deep ground-based or Hubble Space Telescope archival imaging. The deepest absolute BV R magnitude limits are between 4 m and 5 m . We compare these limits with the observed Wolf-Rayet population in the Large Magellanic Cloud and estimate a 16 per cent probability we have failed to detect such a progenitor by chance. Alternatively the progenitors evolve significantly before core-collapse or we have underestimated the extinction towards the progenitors. Reviewing the relative rates and ejecta mass estimates from lightcurve modelling of Ibc SNe, we find both incompatible with Wolf-Rayet stars with initial masses > 25M⊙being the only progenitors. We present binary evolution models that fit these observational constraints. Stars in binaries with initial masses . 20M⊙lose their hydrogen envelopes in binary interactions to become low mass helium stars. They retain a low mass hydrogen envelope until � 10 4 years before core-collapse; hence it is not surprising that galactic analogues have been difficult to identify.

Spectral population synthesis including massive binaries

We have constructed a new code to produce synthetic spectra of stellar populations that includes massive binaries. We have tested this code against the broad-band colours of unresolved young massive stellar clusters in nearby galaxies, the equivalent widths of the Red and Blue Wolf–Rayet bumps in star-forming Sloan Digital Sky Survey galaxies and the ultraviolet and optical spectra of the star-forming regions Tol-A and B in NGC 5398. In each case, we find a good agreement between our models and observations. We find that in general binary populations are bluer and have fewer red supergiants, and thus significantly less flux in the I band and at longer wavelengths, than single star populations. Also we find that Wolf–Rayet stars occur over a wider range of ages up to 107 yr in a stellar population including binaries, increasing the ultraviolet flux and Wolf–Rayet spectral features at later times. In addition, we find that nebula emission contributes significantly to these observed properties and must be considered when comparing stellar models with observations of unresolved stellar populations. We conclude that incorporation of massive stellar binaries can improve the agreement between observations and synthetic spectral synthesis codes, particularly for systems with young stellar populations.

Stellar population effects on the inferred photon density at reionization

The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium (IGM). We quantify the ionizing flux arising from synthetic stellar populations which incorporate the evolution of interacting binary stars. We determine that these show ionizing flux boosted by 60 per cent at 0.05 ≤ Z ≤ 0.3 Z⊙ and a more modest 10–20 per cent at near-solar metallicities relative to star-forming populations in which stars evolve in isolation. The relation of ionizing flux to observables such as 1500 A continuum and ultraviolet spectral slope is sensitive to attributes of the stellar population including age, star formation history and initial mass function (IMF). For a galaxy forming 1 M⊙ yr−1, observed at >100 Myr after the onset of star formation, we predict a production rate of photons capable of ionizing hydrogen, Nion = 1.4 × 1053 s−1 at Z = Z⊙ and 3.5 × 1053 s−1 at 0.1 Z⊙, assuming a Salpeter-like IMF. We evaluate the impact of these issues on the ionization of the IGM, finding that the known galaxy populations can maintain the ionization state of the Universe back to z ∼ 9, assuming that their luminosity functions continue to MUV = −10, and that constraints on the IGM at z ∼ 2–5 can be satisfied with modest Lyman-continuum photon escape fractions of 4–24 per cent depending on assumed metallicity.

Runaway stars as progenitors of supernovae and gamma-ray bursts

When a core-collapse supernova occurs in a binary system, the surviving star as well as the compact remnant emerging from the supernova may reach a substantial space velocity. With binary population synthesis modelling at solar and one-fifth of solar metallicity, we predict the velocities of such runaway stars or binaries. We compile predictions for runaway OB stars, red supergiants and Wolf–Rayet stars, either isolated or with a compact companion. For those stars or binaries which undergo a second stellar explosion we compute their further evolution and the distance travelled until a Type II or Type Ibc supernova or a long or short gamma-ray burst (GRB) occurs. We find our predicted population of OB runaway stars broadly matches the observed population of stars but, to match the fastest observed Wolf–Rayet runaway stars, we require that black holes receive an asymmetric kick upon formation. We find that at solar metallicity Type Ic supernova progenitors travel shorter distances than the progenitors of other supernova types because they are typically more massive and thus have shorter lifetimes. Those of Type IIP supernovae can fly farthest about 48 pc on average at solar metallicity, with about 8 per cent of them reaching 100 pc. In considering the consequences of assuming that the progenitors of long GRBs are spun-up secondary stars that experience quasi-homogeneous evolution, we find that such evolution has a dramatic effect on the population of runaway Wolf– Rayet stars and that some 30 per cent of GRBs could occur a hundred parsec or more from their initial positions. We also consider mergers of double compact object binaries consisting of neutron stars and/or black holes. We find the most common type of visible mergers are neutron star–black hole mergers that are roughly 10 times more common than neutron star–neutron star mergers. All compact mergers have a wide range of merger times from years to Gyr and are predicted to occur 300 times less often than supernovae in the Milky Way. We also find that there may be a population of low-velocity neutron stars that are ejected from a binary rather than by their own natal kick. These neutron stars need to be included when the distribution of neutron star kicks is deduced from observations.

The great escape: how exoplanets and smaller bodies desert dying stars

Mounting discoveries of extrasolar planets orbiting post-main sequence stars motivate studies aimed at understanding the fate of these planets. In the traditional “adiabatic” approximation, a secondary’s eccentricity remains constant during stellar mass loss. Here, we remove this approximation, investigate the full twobody point-mass problem with isotropic mass loss, and illustrate the resulting dynamical evolution. The magnitude and duration of a star’s mass loss combined with a secondary’s initial orbital characteristics might provoke ejection, modest eccentricity pumping, or even circularisation of the orbit. We conclude that Oort clouds and wide-separation planets may be dynamically ejected from 1M⊙ 7M⊙ parent stars during AGB evolution. The vast majority of planetary material which survives a supernova from a 7M⊙ 20M⊙ progenitor will be dynamically ejected from the system, placing limits on the existence of firstgeneration pulsar planets. Planets around > 20M⊙ black hole progenitors may easily survive or readily be ejected depending on the core collapse and superwind models applied. Material ejected during stellar evolution might contribute significantly to the free-floating planetary population.

Modelling the binary progenitor of Supernova 1993J

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code, we investigate potential progenitors for the Type IIb Supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al.) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient, the observations can be reproduced by a system consisting of a 15 M ⊙ primary and a 14 M ⊙ secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a 17 M ⊙ primary and a 16 M ⊙ secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.

Bolometric light curves and explosion parameters of 38 stripped-envelope core-collapse supernovae

Literature data are collated for 38 stripped-envelope core-collapse supernovae (SE SNe; i.e. SNe IIb, Ib, Ic and Ic-BL) that have good light curve coverage in more than one optical band. Using bolometric corrections derived in previous work, the bolometric light curve of each SN is recovered and template bolometric light curves provided. Peak light distributions and decay rates are investigated; SNe subtypes are not cleanly distinguished in this parameter space, although some grouping of types does occur and there is a suggestion of a Phillips-like relation for most SNe Ic-BL. The bolometric light curves are modelled with a simple analytical prescription and compared to results from more detailed modelling. Distributions of the explosion parameters shows the extreme nature of SNe Ic-BL in terms of their 56Ni mass and the kinetic energy, however ejected masses are similar to other subtypes. SNe Ib and Ic have very similar distributions of explosion parameters, indicating a similarity in progenitors. SNe~IIb are the most homogeneous subtype and have the lowest average values for 56Ni mass, ejected mass, and kinetic energy. Ejecta masses for each subtype and SE SNe as a whole are inconsistent with those expected from very massive stars. The majority of the ejecta mass distribution is well described by more moderately massive progenitors in binaries, indicating these are the dominant progenitor channel for SE SNe.

Massive Stars Exploding in a He-Rich Circumstellar Medium -- I. Type Ibn (SN 2006jc-Like) Events

We present new spectroscopic and photometric data of the Type Ibn supernovae 2006jc, 2000er and 2002ao. We discuss the general properties of this recently proposed supernova family, which also includes SN 1999cq. The early-time monitoring of SN 2000er traces the evolution of this class of objects during the first few days after the shock breakout. An overall similarity in the photometric and spectroscopic evolution is found among the members of this group, which would be unexpected if the energy in these core-collapse events was dominated by the interaction between supernova ejecta and circumstellar medium. Type Ibn supernovae appear to be rather normal Type Ib/c supernova explosions which occur within a He-rich circumstellar environment. SNe Ibn are therefore likely produced by the explosion of Wolf–Rayet progenitors still embedded in the He-rich material lost by the star in recent mass-loss episodes, which resemble known luminous blue variable eruptions. The evolved Wolf–Rayet star could either result from the evolution of a very massive star or be the more evolved member of a massive binary system.We also suggest that there are a number of arguments in favour of a Type Ibn classification for the historical SN 1885A (S-Andromedae), previously considered as an anomalous Type Ia event with some resemblance to SN 1991bg.

SN 2009md: another faint supernova from a low-mass progenitor

We present adaptive optics imaging of the core-collapse supernova (SN) 2009md, which we use together with archival Hubble Space Telescope data to identify a coincident progenitor candidate. We find the progenitor to have an absolute magnitude of V=-4.63+0.3-0.4 mag and a colour of V-I= 2.29+0.25-0.39 mag, corresponding to a progenitor luminosity of log L/Ls˜ 4.54 ± 0.19 dex. Using the stellar evolution code STARS, we find this to be consistent with a red supergiant progenitor with M= 8.5+6.5-1.5 Ms. The photometric and spectroscopic evolution of SN 2009md is similar to that of the class of sub-luminous Type IIP SNe; in this paper we compare the evolution of SN 2009md primarily to that of the sub-luminous SN 2005cs. We estimate the mass of 56Ni ejected in the explosion to be (5.4 ± 1.3) × 10-3 Ms from the luminosity on the radioactive tail, which is in agreement with the low 56Ni masses estimated for other sub-luminous Type IIP SNe. From the light curve and spectra, we show the SN explosion had a lower energy and ejecta mass than the normal Type IIP SN 1999em. We discuss problems with stellar evolutionary models, and the discrepancy between low observed progenitor luminosities (log L/Ls˜4.3-5 dex) and model luminosities after the second dredge-up for stars in this mass range, and consider an enhanced carbon burning rate as a possible solution. In conclusion, SN 2009md is a faint SN arising from the collapse of a progenitor close to the lower mass limit for core collapse. This is now the third discovery of a low-mass progenitor star producing a low-energy explosion and low 56Ni ejected mass, which indicates that such events arise from the lowest end of the mass range that produces a core-collapse SN (7-8 Ms).

Massive stars exploding in a He-rich circumstellar medium – III. SN 2006jc: infrared echoes from new and old dust in the progenitor CSM

We present near- (NIR) and mid-infrared (MIR) photometric data of the Type Ibn supernova (SN) 2006jc obtained with the United Kingdom Infrared Telescope (UKIRT), the Gemini North Telescope and the Spitzer Space Telescope between days 86 and 493 post-explosion. We find that the IR behaviour of SN 2006jc can be explained as a combination of IR echoes from two manifestations of circumstellar material. The bulk of the NIR emission arises from an IR echo from newly condensed dust in a cool dense shell (CDS) produced by the interaction of the ejecta outward shock with a dense shell of circumstellar material ejected by the progenitor in a luminous blue variable (LBV)-like outburst about two years prior to the SN explosion. The CDS dust mass reaches a modest 3.0 × 10^(−4) M_⊙ by day 230. While dust condensation within a CDS formed behind the ejecta inward shock has been proposed before for one event (SN 1998S), SN 2006jc is the first one showing evidence for dust condensation in a CDS formed behind the ejecta outward shock in the circumstellar material. At later epochs, a substantial and growing contribution to the IR fluxes arises from an IR echo from pre-existing dust in the progenitor wind. The mass of the pre-existing circumstellar medium (CSM) dust is at least ∼8 × 10^(−3) M_⊙. This paper therefore adds to the evidence that mass-loss from the progenitors of core-collapse SNe could be a major source of dust in the Universe. However, yet again, we see no direct evidence that the explosion of an SN produces anything other than a very modest amount of dust.