Observational Constraints on the Progenitors of Core-Collapse Supernovae: The Case for Missing High-Mass Stars

Over the last 15 years, the supernova community has endeavoured to directly identify progenitor stars for core-collapse supernovae discovered in nearby galaxies. These precursors are often visible as resolved stars in high-resolution images from space-and ground-based telescopes. The discovery rate of progenitor stars is limited by the local supernova rate and the availability and depth of archive images of galaxies, with 18 detections of precursor objects and 27 upper limits. This review compiles these results (from 1999 to 2013) in a distance-limited sample and discusses the implications of the findings. The vast majority of the detections of progenitor stars are of type II-P, II-L, or IIb with one type Ib progenitor system detected and many more upper limits for progenitors of Ibc supernovae (14 in all). The data for these 45 supernovae progenitors illustrate a remarkable deficit of high-luminosity stars above an apparent limit of log L/L-circle dot similar or equal to 5.1 dex. For a typical Salpeter initial mass function, one would expect to have found 13 high-luminosity and high-mass progenitors by now. There is, possibly, only one object in this time-and volume-limited sample that is unambiguously high-mass (the progenitor of SN2009ip) although the nature of that supernovae is still debated. The possible biases due to the influence of circumstellar dust, the luminosity analysis, and sample selection methods are reviewed. It does not appear likely that these can explain the missing high-mass progenitor stars. This review concludes that the community's work to date shows that the observed populations of supernovae in the local Universe are not, on the whole, produced by high-mass (M greater than or similar to 18 M-circle dot) stars. Theoretical explosions of model stars also predict that black hole formation and failed supernovae tend to occur above an initial mass of M similar or equal to 18 M-circle dot. The models also suggest there is no simple single mass division for neutron star or black-hole formation and that there are islands of explodability for stars in the 8-120 M-circle dot range. The observational constraints are quite consistent with the bulk of stars above M similar or equal to 18 M-circle dot collapsing to form black holes with no visible supernovae.

Self-regulated Shocks in Massive Star Binary Systems

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using steady-state models that include a parameterized line force with X-ray ionization dependence. Our primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic X-ray luminosity of the shocks, L X, remains largely unaltered, with the exception of a modest increase at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks and the wind driving, we term this scenario as self-regulated shocks. This effect is found to greatly increase the range of binary separations at which a wind-photosphere collision is likely to occur in systems where the momenta of the two winds are significantly different. Furthermore, the excessive levels of X-ray ionization close to the shocks completely suppress the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveal reasonable agreement in terms of log (L X/L bol). The inclusion of self-regulated shocks improves the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for systems with unequal wind momenta (if considered to be a wind-photosphere collision).

The VLT-FLAMES tarantula survey:X. Evidence for a bimodal distribution of rotational velocities for the single early B-type stars

Aims. Projected rotational velocities (ve sin i) have been estimated for 334 targets in the VLT-FLAMES Tarantula Survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods. Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results. Projected rotational velocities range up to approximately 450 kms-1 and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having 0 <ve <100 km s^-1 and the high velocity component having v_e ∼ 250 km s^-1. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions. The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars.While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations. © ESO 2013.

On the progenitor of the type IIP SN 2013ej in M74

We use natural seeing imaging of SN 2013ej in M74 to identify a progenitor candidate in archival Hubble Space Telescope (HST) + Advanced Camera for Survey images. We find a source coincident with the supernova (SN) in the F814W filter within the total 75 mas (~3 pc astrometric uncertainty; however, the position of the progenitor candidate in contemporaneous F435W and F555W filters is significantly offset. We conclude that the 'progenitor candidate' is in fact two physically unrelated sources; a blue source which is likely unrelated to the SN, and a red source which we suggest exploded as SN 2013ej. Deep images with the same instrument on board HST taken when the SN has faded (in approximately two year's time) will allow us to accurately characterize the unrelated neighbouring source and hence determine the intrinsic flux of the progenitor in three filters.We suggest that the F814W flux is dominated by the progenitor of SN 2013ej, and assuming a bolometric correction appropriate to an M-type supergiant, we estimate that the mass of the progenitor of SN 2013ej was between 8 and 15.5M⊙. 

On the triple peaks of SNHunt248 in NGC 5806

We present our findings on a supernova (SN) impostor, SNHunt248, based on optical and near-IR data spanning ~15 yr before discovery, to ~1 yr post-discovery. The light curve displays three distinct peaks, the brightest of which is at MR ~ −15.0 mag. The post-discovery evolution is consistent with the ejecta from the outburst interacting with two distinct regions of circumstellar material. The 0.5–2.2 μm spectral energy distribution at −740 d is well-matched by a single 6700 K blackbody with log (L/L⊙) ~ 6.1. This temperature and luminosity support previous suggestions of a yellow hypergiant progenitor; however, we find it to be brighter than the brightest and most massive Galactic late-F to early-G spectral type hypergiants. Overall the historical light curve displays variability of up to ~ ± 1 mag. At current epochs (~1 yr post-outburst), the absolute magnitude (MR ~ − 9 mag) is just below the faintest observed historical absolute magnitude ~10 yr before discovery.

Broad balmer wings in BA Hyper/Supergiants distorted by diffuse interstellar bands: five examples in the 30 Doradus Region from the VLT-Flames Tarantula Survey

Extremely broad emission wings at Hβ and Hα have been found in VLT-FLAMES Tarantula Survey data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of Hβ and the shortward wing of Hα. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables (LBVs). No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.

The VLT-FLAMES survey of massive stars: Observations in the Galactic clusters NGC 3293, NGC 4755 and NGC 6611

We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi- Element Spectrograph ( FLAMES) instrument at the Very Large Telescope ( VLT). Here we present observations of 269 Galactic stars with the FLAMES- Giraffe Spectrograph ( R similar or equal to 25 000), in fields centered on the open clusters NGC3293, NGC4755 and NGC6611. These data are supplemented by a further 50 targets observed with the Fibre- Fed Extended Range Optical Spectrograph ( FEROS, R = 48 000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC3293 and NGC4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared.

Chemical abundances and winds of massive stars in M31: a B-type supergiant and a WC star in OB 10

We present high quality spectroscopic data for two massive stars in the OB 10 association of M31, OB 10-64 (B0 la) and OB 10-WRI (WC6). Medium resolution spectra of both stars were obtained using the ISIS spectrograph on the William Herschel Telescope. This is supplemented with Hubble Space Telescope STIS UV spectroscopy and Keck I HIRES data for OB 10-64. A non- local thermodynamic equilibrium (LTE) model atmosphere and abundance analysis for OB 10-64 is presented, indicating that this star has similar photospheric CNO, Mg and Si abundances to solar neighbourhood massive stars. A wind analysis of this early B-type supergiant reveals a mass-loss rate of (M)over dot = 1.6 x 10(-6) M-circle dot yr(-1), and v(infinity) = 1650 km s(-1). The corresponding wind momentum is in good agreement with the wind momentum-luminosity relationship found for Galactic early-B supergiants. Observations of OB 10-WRI are analysed using a non-LTE, line-blanketed code, to reveal approximate stellar parameters of log L/L-circle dot similar to 5.7, T-* - 75 kK, v(infinity) similar to 3000 km s(-1), (M)over dot/(M-circle dot yr(-1)) similar to 10(-4.3) adopting a clumped wind with a filling factor of 10 per cent. Quantitative comparisons are made with the Galactic WC6 star HD 92809 (WR23) revealing that OB 10-WR1 is 0.4 dex more luminous, though it has a much lower C/He ratio (similar to0.1 versus 0.3 for HD 92809). Our study represents the first detailed, chemical model atmosphere analysis for either a B-type supergiant or a Wolf- Rayet (WR) star in Andromeda, and shows the potential of how such studies can provide new information on the chemical evolution of galaxies and the evolution of massive stars in the local Universe.

The VLT-FLAMES survey of massive stars: stellar parameters and rotational velocities in NGC 3293, NGC 4755 and NGC 6611

An analysis is presented of VLT-FLAMES spectroscopy for three Galactic clusters, NGC3293, NGC4755 and NGC6611. Non-LTE model atmosphere calculations have been used to estimate effective temperatures (from either the helium spectrum or the silicon ionization equilibrium) and gravities (from the hydrogen spectrum). Projected rotational velocities have been deduced from the helium spectrum (for fast and moderate rotators) or the metal line spectrum (for slow rotators). The origin of the low gravity estimates for apparently near main sequence objects is discussed and is related to the stellar rotational velocity. The atmospheric parameters have been used to estimate cluster distances (which are generally in good agreement with previous determinations) and these have been used to estimate stellar luminosities and evolutionary masses. The observed Hertzsprung-Russell diagrams are compared with theoretical predictions and some discrepancies including differences in the main sequence luminosities are discussed. Cluster ages have been deduced and evidence for non-coeval star formation is found for all three of the clusters. Projected rotational velocities for targets in the older clusters, NGC3293 and NGC4755, have been found to be systematically larger than those for the field, confirming recent results in other similar age clusters. The distribution of projected rotational velocities are consistent with a Gaussian distribution of intrinsic rotational velocities. For the relatively unevolved targets in the older clusters, NGC3293 and NGC4755, the peak of the velocity distribution would be 250 km s(-1) with a full-width-half-maximum of approximately 180 km s(-1). For NGC6611, the sample size is relatively small but implies a lower mean rotational velocity. This may be evidence for the spin-down effect due to angular momentum loss through stellar winds, although our results are consistent with those found for very young high mass stars. For all three clusters we deduce present day mass functions with Gamma-values in the range of -1.5 to -1.8, which are similar to other young stellar clusters in the Milky Way.

The VLT-FLAMES survey of massive stars: mass loss and rotation of early-type stars in the SMC

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method (Mokiem et al. 2005, A&A, 441, 711), which combines the stellar atmosphere code FASTWIND (Puls et al. 2005, A&A, 435, 669) with the genetic algorithm based optimisation routine PIKAIA (Charbonneau 1995, ApJS, 101, 309). Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The automated method allows for a detailed determination of the projected rotational velocities. The present day v(r) sin i distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v(r)) distribution that can be characterised by a mean velocity of about 160-190 km s(-1) and an effective half width of 100-150 km s(-1). The vr distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct (Maeder & Meynet 2001, A&A, 373, 555), the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring vr. The mass loss rates of the SMC objects having luminosities of log L-star/L-circle dot greater than or similar to 5.4 are in excellent agreement with predictions by Vink et al. (2001, A&A, 369, 574). However, for lower luminosity stars the winds are too weak to determine. M accurately from the optical spectrum. Three targets were classified as Vz stars, two of which are located close to the theoretical zero-age main sequence. Three lower luminosity targets that were not classified as Vz stars are also found to lie near the ZAMS. We argue that this is related to a temperature effect inhibiting cooler from displaying the spectral features required for the Vz luminosity class.

The VLT-FLAMES survey of massive stars: evolution of surface N abundances and effective temperature scales in the Galaxy and Magellanic Clouds

We present an analysis of high resolution VLT-FLAMES spectra of 61 B-type stars with relatively narrow-lined spectra located in 4 fields centered on the Milky Way clusters; NGC 3293 and NGC 4755 and the Large and Small Magellanic cloud clusters; NGC 2004 and NGC 330. For each object a quantitative analysis was carried out using the non-LTE model atmosphere code TLUSTY; resulting in the determination of their atmospheric parameters and photospheric abundances of the dominant metal species (C, N, O, Mg, Si, Fe). The results are discussed in relation to our earlier work on 3 younger clusters in these galaxies; NGC 6611, N11 and NGC 346 paying particular attention to the nitrogen abundances which are an important probe of the role of rotation in the evolution of stars. This work along with that of the younger clusters provides a consistent dataset of abundances and atmospheric parameters for over 100 B-type stars in the three galaxies. We provide effective temperature scales for B-type dwarfs in all three galaxies and for giants and supergiants in the SMC and LMC. In each galaxy a dependence on luminosity is found between the three classes with the unevolved dwarf objects having significantly higher effective temperatures. A metallicity dependence is present between the SMC and Galactic dwarf objects, and whilst the LMC stars are only slightly cooler than the SMC stars, they are significantly hotter than their Galactic counterparts.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds

Aims.We aim to provide the atmospheric parameters and rotational velocities for a large sample of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Methods: Atmospheric parameters, stellar masses, and rotational velocities have been estimated for approximately 250 early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006; Trundle et al. 2007) from the same survey, providing a sample of ~400 early-type objects. Results: Comparison of the rotational velocities with evolutionary tracks suggests that the end of core hydrogen burning occurs later than currently predicted and we argue for an extension of the evolutionary tracks. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities nor their position on the H-R diagram are predicted. We suggest that blue loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample, we find that (at 3s confidence level) massive stars (above 8 M?) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution that is a Gaussian peaking at 175 km s-1 (SMC) and 100 km s-1 (LMC) with a 1/e half width of 150 km s-1. We find that in NGC 346 in the SMC, the 10-25 M? main-sequence stars appear to rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through angular momentum loss via stellar winds at SMC metallicity, hence this could be a reflection of mass dependent birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25 M? stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 M? at low metallicities.

The VLT-FLAMES Survey of massive stars: Rotation and nitrogen enrichment as the key to understanding massive star evolution

Rotation has become an important element in evolutionary models of massive stars, specifically via the prediction of rotational mixing. Here we study a sample of stars, including rapid rotators, to constrain such models and use nitrogen enrichments as a probe of the mixing process. Chemical compositions (C, N, O, Mg, and Si) have been estimated for 135 early B-type stars in the Large Magellanic Cloud with projected rotational velocities up to similar to 300 km s(-1) using a non-LTE TLUSTY model atmosphere grid. Evolutionary models, including rotational mixing, have been generated attempting to reproduce these observations by adjusting the overshooting and rotational mixing parameters and produce reasonable agreement with 60% of our core hydrogen burning sample. We find (excluding known binaries) a significant population of highly nitrogen-enriched intrinsic slow rotators (nu sin i less than or similar to 50 km s(-1)) incompatible with our models (similar to 20% of the sample). Furthermore, while we find fast rotators with enrichments in agreement with the models, the observation of evolved (dex) fast rotators (log g < 3.7 dex) that are relatively unenriched (a further similar to 20% of the sample) challenges the concept of rotational mixing. We also find that 70% of our blue supergiant sample cannot have evolved directly from the hydrogen-burning main sequence. We are left with a picture where invoking binarity and perhaps fossil magnetic fields is required to understand the surface properties of a population of massive main- sequence stars.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds (vol 457, pg 265, 2006)

We correct the estimates of the dispersions in the rotational velocities for early-type stars in our Galaxy (Dufton et al. 2006, A&A, 457, 265) and the Magellanic Clouds (Hunter et al. 2008, A&A, 479, 541). The corrected values are pi(1/4) (i.e. approximately 33%) larger than those published in the original papers.

The VLT-FLAMES survey of massive stars: constraints on stellar evolution from the chemical compositions of rapidly rotating Galactic and Magellanic Cloud B-type stars

Aims. We have previously analysed the spectra of 135 early B-type stars in the Large Magellanic Cloud (LMC) and found several groups of stars that have chemical compositions that conflict with the theory of rotational mixing. Here we extend this study to Galactic and Small Magellanic Cloud (SMC) metallicities.