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.

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.

The death of massive stars - II. Observational constraints on the progenitors of Type Ibc supernovae

The progenitors of many Type II core-collapse supernovae (SNe) have now been identified directly on pre-discovery imaging. Here, we present an extensive search for the progenitors of Type Ibc SNe in all available pre-discovery imaging since 1998. There are 12 Type Ibc SNe with no detections of progenitors in either deep ground-based or Hubble Space Telescope archival imaging. The deepest absolute BVR magnitude limits are between -4 and - 5 mag. We compare these limits with the observed Wolf-Rayet population in the Large Magellanic Cloud and estimate a 16 per cent probability that 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 light-curve modelling of Ibc SNe, we find both incompatible with Wolf-Rayet stars with initial masses >25 M⊙ being the only progenitors. We present binary evolution models that fit these observational constraints. Stars in binaries with initial masses ≲ 20 M⊙ lose their hydrogen envelopes in binary interactions to become low-mass helium stars. They retain a low-mass hydrogen envelope until ≈104 yr before core-collapse; hence, it is not surprising that Galactic analogues have been difficult to identify.

Possible binary progenitors for the Type Ib supernova iPTF13bvn

Cao et al. reported a possible progenitor detection for the Type Ib supernovae iPTF13bvn for the first time. We find that the progenitor is in fact brighter than the magnitudes previously reported by approximately 0.7-0.2 mag with a larger error in the bluer filters. We compare our new magnitudes to our large set of binary evolution models and find that many binary models with initial masses in the range of 10-20M(circle dot) match this new photometry and other constraints suggested from analysing the supernova. In addition, these lower mass stars retain more helium at the end of the model evolution indicating that they are likely to be observed as Type Ib supernovae rather than their more massive, Wolf-Rayet counter parts. We are able to rule out typical Wolf-Rayet models as the progenitor because their ejecta masses are too high and they do not fit the observed SED unless they have a massive companion which is the observed source at the supernova location. Therefore only late-time observations of the location will truly confirm if the progenitor was a helium giant and not a Wolf-Rayet star.

Supersolar Ni/Fe production in the Type IIP SN 2012ec

SN 2012ec is a Type IIP supernova (SN) with a progenitor detection and comprehensive photospheric phase observational coverage. Here, we present Very Large Telescope and Public ESO Spectroscopic Survey of Transient Objects observations of this SN in the nebular phase. We model the nebular [O I] lambda lambda 6300, 6364 lines and find their strength to suggest a progenitor main-sequence mass of 13-15 M-circle dot. SN2012ec is unique among hydrogen-rich SNe in showing a distinct line of stable nickel [Ni II] lambda 7378. This line is produced by Ni-58, a nuclear burning ash whose abundance is a sensitive tracer of explosive burning conditions. Using spectral synthesis modelling, we use the relative strengths of [Ni II] lambda 7378 and [Fe II] lambda 7155 (the progenitor of which is Ni-56) to derive a Ni/Fe production ratio of 0.20 +/- 0.07 (by mass), which is a factor 3.4 +/- 1.2 times the solar value. High production of stable nickel is confirmed by a strong [Ni II] 1.939 mu m line. This is the third reported case of a core-collapse SN producing a Ni/Fe ratio far above the solar value, which has implications for core-collapse explosion theory and galactic chemical evolution models.

Deflagrations in hybrid CONe white dwarfs: a route to explain the faint Type Iax supernova 2008ha

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon-oxygen core surrounded by an oxygen-neon mantle. Being born with masses similar to 1.1 M-aS (TM), hybrid WDs in a binary system may easily approach the Chandrasekhar mass (M-Ch) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-M-Ch hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 M-aS (TM) of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near-M-Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-M-Ch bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching M-Ch to be of the order of 1 per cent of the Galactic SN Ia rate.

Sher 25: pulsating but apparently alone

The blue supergiant Sher 25 is surrounded by an asymmetric, hourglass-shaped circumstellar nebula, which shows similarities to the triple-ring structure seen around SN 1987A. From optical spectroscopy over six consecutive nights, we detect periodic radial velocity variations in the stellar spectrum of Sher 25 with a peak-to-peak amplitude of ~ 12 km s^-1 on a time-scale of about 6 d, confirming the tentative detection of similar variations by Hendry et al. From consideration of the amplitude and time-scale of the signal, coupled with observed line profile variations, we propose that the physical origin of these variations is related to pulsations in the stellar atmosphere, rejecting the previous hypothesis of a massive, short-period binary companion. The radial velocities of two other blue supergiants with similar bipolar nebulae, SBW1 and HD 168625, were also monitored over the course of six nights, but these did not display any significant radial velocity variations. 

The nebular spectra of SN 2012aw and constraints on stellar nucleosynthesis from oxygen emission lines

We present nebular-phase optical and near-infrared spectroscopy of the Type IIP supernova SN 2012aw combined with non-local thermodynamic equilibrium radiative transfer calculations applied to ejecta from stellar evolution/explosion models. Our spectral synthesis models generally show good agreement with the ejecta from a M_ZAMS = 15 M_⊙progenitor star. The emission lines of oxygen, sodium, and magnesium are all consistent with the nucleosynthesis in a progenitor in the 14-18 M_⊙ range.We also demonstrate how the evolution of the oxygen cooling lines of [O I] λ5577, [O I] λ6300, and [O I] λ6364 can be used to constrain the mass of oxygen in the non-molecularly cooled ashes to < 1 M_⊙, independent of the mixing in the ejecta. This constraint implies that any progenitor model of initial mass greater than 20 M_⊙ would be difficult to reconcile with the observed line strengths. A stellar progenitor of around M_ZAMS = 15 M_⊙ can consistently explain the directly measured luminosity of the progenitor star, the observed nebular spectra, and the inferred pre-supernova mass-loss rate.We conclude that there is still no convincing example of a Type IIP supernova showing the nucleosynthesis products expected from an M_ZAMS > 20 M_⊙ progenitor. © 2014 The Author. Published by Oxford University Press on behalf of the Royal Astronomical Society.

On the progenitor of the Type Ic SN 2013dk in the Antennae galaxies

We report the results of our search for the progenitor candidate of SN 2013dk, a Type Ic supernova (SN) that exploded in the Antennae galaxy system. We compare pre-explosion Hubble Space Telescope (HST) archival images with SN images obtained using adaptive optics at the ESO Very Large Telescope. We isolate the SN position to within 3σ uncertainty radius of 0.02 arcsec and show that there is no detectable point source in any of the HST filter images within the error circle. We set an upper limit to the absolute magnitude of the progenitor to be MF555W ≳ -5.7, which does not allow Wolf-Rayet (WR) star progenitors to be ruled out. A bright source appears 0.17 arcsec away, which is either a single bright supergiant or compact cluster, given its absolute magnitude of MF555W = -9.02 ± 0.28 extended wings and complex environment. However, even if this is a cluster, the spatial displacement of SN 2013dk means that its membership is not assured. The strongest statement that we can make is that in the immediate environment of SN 2013dk (within 10 pc or so), we find no clear evidence of either a point source coincident with the SN or a young stellar cluster that could host a massive WR progenitor.

Type Ia supernovae from exploding oxygen-neon white dwarfs

Context. The progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M-Ch, WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor.

Aims. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia.

Methods. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M-Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels.

Results. The ejecta structures of our simulated detonations in sub-M-Ch, ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models.

Conclusions. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (similar to 1.2 M-circle dot) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.

YSOVAR: Mid-infrared Variability Among YSOs in the Star Formation Region GGD12-15

We present an IR-monitoring survey with the Spitzer Space Telescope of the star-forming region GGD 12-15. More than 1000 objects were monitored, including about 350 objects within the central 5′, which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability project. The region was also the subject of a contemporaneous 67 ks Chandra observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5′ are young stellar objects (YSOs) based on a combination of their spectral energy distribution, IR variability, and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-infrared (mid-IR) variability. There are 39 periodic sources, and all but one is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded YSOs. Most of the Class I/ II objects exhibit redder colors in a fainter state, which is compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR light curves. However, we find that X-ray detected Class II sources have longer timescales for change in the MIR than a similar set of non-X-ray detected Class IIs.

Boron abundances in B-type stars: A test of rotational depletion during main-sequence evolution

Boron abundances have been derived for seven main-sequence B- type stars from Hubble Space Telescope STIS spectra around the B III lambda2066 line. In two stars, boron appears to be undepleted with respect to the presumed initial abundance. In one star, boron is detectable but is clearly depleted. In the other four stars, boron is undetectable, implying depletions of 1-2 dex. Three of these four stars are nitrogen enriched, but the fourth shows no enrichment of nitrogen. Only rotationally induced mixing predicts that boron depletions are unaccompanied by nitrogen enrichments. The inferred rate of boron depletion from our observations is in good agreement with these predictions. Other boron-depleted nitrogen-normal stars are identified from the literature. In addition, several boron- depleted nitrogen-rich stars are identified, and while all fall on the boron-nitrogen trend predicted by rotationally induced mixing, a majority have nitrogen enrichments that are not uniquely explained by rotation. The spectra have also been used to determine iron group (Cr, Mn, Fe, and Ni) abundances. The seven B-type stars have near-solar iron group abundances, as expected for young stars in the solar neighborhood. We have also analyzed the halo B-type star PG 0832 + 676. We find [Fe/H] = -0.88 +/- 0.10, and the absence of the B III line gives the upper limit [B/H] <-2.5. These and other published abundances are used to infer the star's evolutionary status as a post-asymptotic giant branch star.

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

The VLT-FLAMES survey of massive stars: wind properties and evolution of hot massive stars in the Large Magellanic Cloud

We have studied the optical spectra of a sample of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star forming region N11. Our observations sample the central associations of LH9 and LH10, and the surrounding regions. Stellar parameters are determined using an automated fitting method ( Mokiem et al. 2005), which combines the stellar atmosphere code fastwind ( Puls et al. 2005) with the genetic algorithm based optimisation routine PIKAIA ( Charbonneau 1995). We derive an age of 7.0 +/- 1.0 and 3.0 +/- 1.0 Myr for LH9 and LH10, respectively. The age difference and relative distance of the associations are consistent with a sequential star formation scenario in which stellar activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2. From helium and hydrogen line fitting we find the hottest three of these stars to be similar to 49- 54 kK ( compared to similar to 45- 46 kK for O3 stars). Detailed determination of the helium mass fraction reveals that the masses of helium enriched dwarfs and giants derived in our spectroscopic analysis are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally enhanced mixing leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. ( 2006). For bright giants and supergiants no such mass discrepancy is found; these stars therefore appear to follow tracks of modestly or non-rotating objects. The set of programme stars was sufficiently large to establish the mass loss rates of OB stars in this Z similar to 1/2 Z(circle dot) environment sufficiently accurate to allow for a quantitative comparison with similar objects in the Galaxy and the SMC. The mass loss properties are found to be intermediate to massive stars in the Galaxy and SMC. Comparing the derived modified wind momenta D-mom as a function of luminosity with predictions for LMC metallicities by Vink et al. ( 2001) yields good agreement in the entire luminosity range that was investigated, i.e. 5.0