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.

The massive binary companion star to the progenitor of supernova 1993J

The massive star that underwent a collapse of its core to produce supernova (SN)1993J was subsequently identified as a non-variable red supergiant star in images of the galaxy M81 taken before explosion(1, 2). It showed an excess in ultraviolet and B-band colours, suggesting either the presence of a hot, massive companion star or that it was embedded in an unresolved young stellar association1. The spectra of SN1993J underwent a remarkable transformation from the signature of a hydrogen-rich type II supernova to one of a helium-rich (hydrogen-deficient) type Ib(3, 4). The spectral and photometric peculiarities were best explained by models in which the 13�20 solar mass supergiant had lost almost its entire hydrogen envelope to a close binary companion(5, 6, 7), producing a 'type IIb' supernova, but the hypothetical massive companion stars for this class of supernovae have so far eluded discovery. Here we report photometric and spectroscopic observations of SN1993J ten years after the explosion. At the position of the fading supernova we detect the unambiguous signature of a massive star: the binary companion to the progenitor.

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.

Australia Telescope Compact Array 1.2cm Observations of the Massive Star Forming Region G305.2+0.2

We report on Australia Telescope Compact Array observations of the massive star-forming region G305.2+0.2 at 1.2 cm. We detected emission in five molecules towards G305A, confirming its hot core nature. We determined a rotational temperature of 26 K for methanol. A non-local thermodynamic equilibrium excitation calculation suggests a kinematic temperature of the order of 200 K. A time-dependent chemical model is also used to model the gas-phase chemistry of the hot core associated with G305A. A comparison with the observations suggest an age of between 2 × 104 and 1.5 × 105 yr. We also report on a feature to the south-east of G305A which may show weak Class I methanol maser emission in the line at 24.933 GHz. The more evolved source G305B does not show emission in any of the line tracers, but strong Class I methanol maser emission at 24.933 GHz is found 3 arcsec to the east. Radio continuum emission at 18.496 GHz is detected towards two H ii regions. The implications of the non-detection of radio continuum emission towards G305A and G305B are also discussed.

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 massive neutron star or low-mass black hole in 2S 0921-630

We report on the optical spectroscopy of the eclipsing halo low-mass X-ray binary 2S 0921-630, which reveals the absorption-line radial velocity curve of the K0 III secondary star with a semiamplitude K-2=92.89+/-3.84 km s(-1), a systemic velocity gamma=34.9+/-3.3 km s(-1), and an orbital period P-orb of 9.0035+/-0.0029 days (1 sigma). Given the quality of the data, we find no evidence for the effects of X-ray irradiation. Using the previously determined rotational broadening of the mass donor and applying conservative limits on the orbital inclination, we constrain the compact object mass to be 2.0-4.3 M-circle dot (1 sigma), ruling out a canonical neutron star at the 99% level. Since the nature of the compact object is unclear, this mass range implies that the compact object is either a low-mass black hole with a mass slightly higher than the maximum possible neutron star mass (2.9 M-circle dot) or a massive neutron star. If the compact object is a black hole, it confirms the prediction of the existence of low-mass black holes, while if the object is a massive neutron star, its high mass severely constrains the equation of state of nuclear matter.

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.

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

We present the results of a 10.5-yr, volume-limited (28-Mpc) search for supernova (SN) progenitor stars. In doing so we compile all SNe discovered within this volume (132, of which 27 per cent are Type Ia) and determine the relative rates of each subtype from literature studies. The core-collapse SNe break down into 59 per cent II-P and 29 per cent Ib/c, with the remainder being IIb (5 per cent), IIn (4 per cent) and II-L (3 per cent). There have been 20 II-P SNe with high-quality optical or near-infrared pre-explosion images that allow a meaningful search for the progenitor stars. In five cases they are clearly red supergiants, one case is unconstrained, two fall on compact coeval star clusters and the other twelve have no progenitor detected. We review and update all the available data for the host galaxies and SN environments (distance, metallicity and extinction) and determine masses and upper mass estimates for these 20 progenitor stars using the STARS stellar evolutionary code and a single consistent homogeneous method. A maximum likelihood calculation suggests that the minimum stellar mass for a Type II-P to form is m(min) = 8.5(-1.5)(+1) M-circle dot and the maximum mass for II-P progenitors is m(max) = 16.5 +/- 1.5 M-circle dot, assuming a Salpeter initial mass function holds for the progenitor population (in the range Gamma = -1.35(-0.7)(+0.3)). The minimum mass is consistent with current estimates for the upper limit to white dwarf progenitor masses, but the maximum mass does not appear consistent with massive star populations in Local Group galaxies. Red supergiants in the Local Group have masses up to 25 M-circle dot and the minimum mass to produce a Wolf-Rayet star in single star evolution (between solar and LMC metallicity) is similarly 25-30 M-circle dot. The reason we have not detected any high-mass red supergiant progenitors above 17 M-circle dot is unclear, but we estimate that it is statistically significant at 2.4 sigma confidence. Two simple reasons for this could be that we have systematically underestimated the progenitor masses due to dust extinction or that stars between 17-25 M-circle dot produce other kinds of SNe which are not II-P. We discuss these possibilities and find that neither provides a satisfactory solution. We term this discrepancy the 'red supergiant problem' and speculate that these stars could have core masses high enough to form black holes and SNe which are too faint to have been detected. We compare the Ni-56 masses ejected in the SNe to the progenitor mass estimates and find that low-luminosity SNe with low Ni-56 production are most likely to arise from explosions of low-mass progenitors near the mass threshold that can produce a core-collapse.

Rotating massive main-sequence stars II. Simulating a population of LMC early B-type stars as a test of rotational mixing

Context. Rotational mixing in massive stars is a widely applied concept, with far-reaching consequences for stellar evolution, nucleosynthesis, and stellar explosions.