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).

Massive open star clusters using the VVV survey III. A young massive cluster at the far edge of the Galactic bar

Context. Young massive clusters are key to map the Milky Way’s structure, and near-infrared large area sky surveys have contributed strongly to the discovery of new obscured massive stellar clusters. Aims. We present the third article in a series of papers focused on young and massive clusters discovered in the VVV survey. This article is dedicated to the physical characterization of VVV CL086, using part of its OB-stellar population. Methods. We physically characterized the cluster using JHKS near-infrared photometry from ESO public survey VVV images, using the VVV-SkZ pipeline, and near-infrared K-band spectroscopy, following the methodology presented in the first article of the series. Results. Individual distances for two observed stars indicate that the cluster is located at the far edge of the Galactic bar. These stars, which are probable cluster members from the statistically field-star decontaminated CMD, have spectral types between O9 and B0 V. According to our analysis, this young cluster (1.0 Myr < age < 5.0 Myr) is located at a distance of 11+5-6 kpc, and we estimate a lower limit for the cluster total mass of (2.8+1.6-1.4) · 103 M⊙. It is likely that the cluster contains even earlier and more massive stars.

MY Camelopardalis, a very massive merger progenitor

Context. The early-type binary MY Cam belongs to the young open cluster Alicante 1, embedded in Cam OB3. Aims. MY Cam consists of two early-O type main-sequence stars and shows a photometric modulation suggesting an orbital period slightly above one day. We intend to confirm this orbital period and derive orbital and stellar parameters. Methods. Timing analysis of a very exhaustive (4607 points) light curve indicates a period of 1.1754514 ± 0.0000015 d. High-resolution spectra and the cross-correlation technique implemented in the todcor program were used to derive radial velocities and obtain the corresponding radial velocity curves for MY Cam. Modelling with the stellar atmosphere code fastwind was used to obtain stellar parameters and create templates for cross-correlation. Stellar and orbital parameters were derived using the Wilson-Devinney code, such that a complete solution to the binary system could be described. Results. The determined masses of the primary and secondary stars in MY Cam are 37.7 ± 1.6 and 31.6 ± 1.4M⊙, respectively. The corresponding temperatures, derived from the model atmosphere fit, are 42 000 and 39 000 K, with the more massive component being hotter. Both stars are overfilling their Roche lobes, sharing a common envelope. Conclusions. MY Cam contains the most massive dwarf O-type stars found so far in an eclipsing binary. Both components are still on the main sequence, and probably not far from the zero-age main sequence. The system is a likely merger progenitor, owing to its very short period.

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.

A study of the large-scale structure of the wind of WR 134

Diverses méthodes ont été utilisées pour étudier les étoiles Wolf-Rayet (WR) dans le but de comprendre les phénomènes physiques variés qui prennent place dans leur vent dense. Pour étudier la variabilité qui n'est pas strictement périodique et ayant des caractéristiques différentes d'une époque à l'autre, il faut observer pendant des périodes de temps suffisamment longues en adopter un échantillonnage temporel élevé pour être en mesure d'identifier les phénomènes physiques sous-jacents. À l'été 2013, des astronomes professionnels et amateurs du monde entier ont contribué à une campagne d'observation de 4 mois, principalement en spectroscopie, mais aussi en photométrie, polarimétrie et en interférométrie, pour observer les 3 premières étoiles Wolf-Rayet découvertes: WR 134 (WN6b), WR 135 (WC8) et WR 137 (WC7pd + O9). Chacune de ces étoiles est intéressante à sa manière, chacune présentant une variété différente de structures dans son vent. Les données spectroscopiques de cette campagne ont été réduites et analysées pour l'étoile présumée simple WR 134 pour mieux comprendre le comportement de sa variabilité périodique à long terme dans le cadre d'une étude des régions d'interactions en corotation (CIRs) qui se retrouvent dans son vent. Les résultats de cette étude sont présentés dans ce mémoire.

He II lambda 4686 IN eta CARINAE: COLLAPSE OF THE WIND-WIND COLLISION REGION DURING PERIASTRON PASSAGE

The periodic spectroscopic events in eta Carinae are now well established and occur near the periastron passage of two massive stars in a very eccentric orbit. Several mechanisms have been proposed to explain the variations of different spectral features, such as an eclipse by the wind-wind collision (WWC) boundary, a shell ejection from the primary star or accretion of its wind onto the secondary. All of them have problems explaining all the observed phenomena. To better understand the nature of the cyclic events, we performed a dense monitoring of eta Carinae with five Southern telescopes during the 2009 low-excitation event, resulting in a set of data of unprecedented quality and sampling. The intrinsic luminosity of the He II lambda 4686 emission line (L similar to 310 L-circle dot) just before periastron reveals the presence of a very luminous transient source of extreme UV radiation emitted in the WWC region. Clumps in the primary's wind probably explain the flare-like behavior of both the X-ray and He II lambda 4686 light curves. After a short-lived minimum, He II lambda 4686 emission rises again to a new maximum, when X-rays are still absent or very weak. We interpret this as a collapse of the WWC onto the "surface" of the secondary star, switching off the hard X-ray source and diminishing the WWC shock cone. The recovery from this state is controlled by the momentum balance between the secondary's wind and the clumps in the primary's wind.

The masses of the neutron and donor star in the high-mass X-ray binary IGR J18027-2016

Aims. We report near-infrared observations of the supergiant donor to the eclipsing high mass X-ray binary pulsar IGR J18027-2016. We aim to determine its spectral type and measure its radial velocity curve and hence determine the stellar masses of the components. Methods. ESO/VLT observations of the donor utilising the NIR spectrograph ISAAC were obtained in the H and K bands. The multi-epoch H band spectra were cross-correlated with RV templates in order to determine a radial solution for the system. Results. The spectral type of the donor was confirmed as B0-1 I. The radial velocity curve constructed has a semi-amplitude of 23.8 ± 3.1 km s-1. Combined with other measured system parameters, a dynamically determined neutron star mass of 1.4  ±  0.2–1.6  ±  0.3 M⊙ is found. The mass range of the B0-B1 I donor was 18.6  ±  0.8–21.8  ±  2.4 M⊙. These lower and upper limits were obtained under the assumption that the system is viewed edge-on (i = 90° with β = 0.89) for the lower limit and the donor fills its Roche lobe (β = 1 with i = 73.1°) for the upper limit respectively.

The evolution and masses of the neutron star and donor star in the high mass X-ray binary OAO 1657−415

We report near-infrared radial velocity (RV) measurements of the recently identified donor star in the high mass X-ray binary (HMXB) system OAO 1657−415 obtained in the H band using ISAAC on the Very Large Telescope. Cross-correlation methods were employed to construct a RV curve with a semi-amplitude of 22.1 ± 3.5 km s−1. Combined with other measured parameters of this system it provides a dynamically determined neutron star (NS) mass of 1.42 ± 0.26 M⊙ and a mass of 14.3 ± 0.8 M⊙ for the Ofpe/WN9 highly evolved donor star. OAO 1657−415 is an eclipsing HMXB pulsar with the largest eccentricity and orbital period of any within its class. Of the 10 known eclipsing X-ray binary pulsars OAO 1657−415 becomes the ninth with a dynamically determined NS mass solution and only the second in an eccentric system. Furthermore, the donor star in OAO 1657−415 is much more highly evolved than the majority of the supergiant donors in other HMXBs, joining a small but growing list of HMXBs donors with extensive hydrogen depleted atmospheres. Considering the evolutionary development of OAO 1657−415, we have estimated the binding energy of the envelope of the mass donor and find that there is insufficient energy for the removal of the donor’s envelope via spiral-in, ruling out a common envelope evolutionary scenario. With its non-zero eccentricity and relatively large orbital period the identification of a definitive evolutionary pathway for OAO 1657−415 remains problematic, we conclude by proposing two scenarios which may account for OAO 1657−415 current orbital configuration.

Orbital and Physical Properties of the σ Ori Aa, Ab, B Triple System

We provide a complete characterization of the astrophysical properties of the σ Ori Aa, Ab, B hierarchical triple system and an improved set of orbital parameters for the highly eccentric σ Ori Aa, Ab spectroscopic binary. We compiled a spectroscopic data set comprising 90 high-resolution spectra covering a total time span of 1963 days. We applied the Lehman-Filhés method for a detailed orbital analysis of the radial velocity curves and performed a combined quantitative spectroscopic analysis of the σ Ori Aa, Ab, B system by means of the stellar atmosphere code FASTWIND. We used our own plus other available information on photometry and distance to the system for measuring the radii, luminosities, and spectroscopic masses of the three components. We also inferred evolutionary masses and stellar ages using the Bayesian code BONNSAI. The orbital analysis of the new radial velocity curves led to a very accurate orbital solution of the σ Ori Aa, Ab pair. We provided indirect arguments indicating that σ Ori B is a fast-rotating early B dwarf. The FASTWIND+BONNSAI analysis showed that the Aa, Ab pair contains the hottest and most massive components of the triple system while σ Ori B is a bit cooler and less massive. The derived stellar ages of the inner pair are intriguingly younger than the one widely accepted for the σ Orionis cluster, at 3 ± 1 Ma. The outcome of this study will be of key importance for a precise determination of the distance to the σ Orionis cluster, the interpretation of the strong X-ray emission detected for σ Ori Aa, Ab, B, and the investigation of the formation and evolution of multiple massive stellar systems and substellar objects.

A new survey of cool supergiants in the Magellanic Clouds

Aims. In this study we conduct a pilot program aimed at the red supergiant population of the Magellanic Clouds. We intend to extend the current known sample to the unexplored low end of the brightness distribution of these stars, building a more representative dataset with which to extrapolate their behaviour to other Galactic and extra-galactic environments. Methods. We select candidates using only near infrared photometry, and with medium resolution multi-object spectroscopy, we perform spectral classification and derive their line-of-sight velocities, confirming the nature of the candidates and their membership in the clouds. Results. Around two hundred new red supergiants have been detected, hinting at a yet to be observed large population. Using near- and mid-infrared photometry we study the brightness distribution of these stars, the onset of mass-loss, and the effect of dust in their atmospheres. Based on this sample, new a priori classification criteria are investigated, combining mid- and near-infrared photometry to improve the observational efficiency of similar programs to this.

The donor star of the X-ray pulsar X1908+075

High-mass X-ray binaries consist of a massive donor star and a compact object. While several of those systems have been well studied in X-rays, little is known for most of the donor stars as they are often heavily obscured in the optical and ultraviolet regime. There is an opportunity to observe them at infrared wavelengths, however. The goal of this study is to obtain the stellar and wind parameters of the donor star in the X1908+075 high-mass X-ray binary system with a stellar atmosphere model to check whether previous studies from X-ray observations and spectral morphology lead to a sufficient description of the donor star. We obtained H- and K-band spectra of X1908+075 and analysed them with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. For the first time, we calculated a stellar atmosphere model for the donor star, whose main parameters are: Mspec = 15 ± 6 M⊙, T∗ = 23-3+6 kK, log geff = 3.0 ± 0.2 and log L/L⊙ = 4.81 ± 0.25. The obtained parameters point towards an early B-type (B0–B3) star, probably in a supergiant phase. Moreover we determined a more accurate distance to the system of 4.85 ± 0.50 kpc than the previously reported value.

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.