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 empirical metallicity dependence of the mass-loss rate of O- and early B-type stars

We present a comprehensive study of the observational dependence of the mass-loss rate in stationary stellar winds of hot massive stars on the metal content of their atmospheres. The metal content of stars in the Magellanic Clouds is discussed, and a critical assessment is given of state-of-the-art mass-loss determinations of OB stars in these two satellite systems and the Milky-Way. Assuming a power-law dependence of mass loss on metal content,. M. Z(m), and adopting a theoretical relation between the terminal flow velocity and metal content, v(infinity). Z(0.13) (Leitherer et al. 1992, ApJ, 401, 596), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 105.2 L circle dot. Within the errors, this result is in agreement with the prediction m = 0.69 +/- 0.10 by Vink et al. (2001, A& A, 369, 574). Absolute empirical values for the mass loss, based on Ha and ultraviolet (UV) wind lines, are found to be a factor of two higher than predictions in this high luminosity regime. If this difference is attributed to inhomogeneities in the wind, and this clumping does not impact the predictions, this would imply that luminous O and early-B stars have clumping factors in their Ha and UV line forming regions of about a factor of four. For lower luminosity stars, the winds are so weak that their strengths can generally no longer be derived from optical spectral lines (essentially Ha) and one must currently rely on the analysis of UV lines. We confirm that in this low-luminosity domain the observed Galactic WLR is found to be much steeper than expected from theory (although the specific sample is rather small), leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L similar to 10(4.75) L circle dot, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars would turn out to be overestimated as a result of wind clumping, but the degree of clumping would be rather independent of metallicity, the scalings derived in this study are expected to remain correct.

Mass-loss rates for hot luminous stars:The influence of line branching

The effect of photon frequency redistribution by line branching on mass-loss rates for hot luminous stars is investigated. Monte Carlo simulations are carried out for a range of OB star models which show that previous mass-loss calculations which neglect non-resonance line scattering overestimate mass-loss rates for luminous O stars by ~20 per cent. For luminous B stars the effect is somewhat larger, typically ~50 per cent. A Wolf-Rayet star model is used to investigate line branching in the strong wind limit. In this case the effect of line branching is much greater, giving mass-loss rates that are smaller by a factor ~3 from computations which neglect branching.

PTF11iqb: cool supergiant mass-loss that bridges the gap between Type IIn and normal supernovae

PTF11iqb was initially classified as a TypeIIn event caught very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on day 2, but the narrow emission weakened quickly and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Halpha emission exhibited a complex, multipeaked profile reminiscent of SN1998S. In terms of spectroscopic evolution, we find that PTF11iqb was a near twin of SN~1998S, although with weaker interaction with circumstellar material (CSM) at early times, and stronger CSM interaction at late times. We interpret the spectral changes as caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped by the expanding SN ejecta photosphere, but then revealed again after the end of the plateau when the photosphere recedes. The light curve can be matched with a simple model for weak CSM interaction added to the light curve of a normal SN~II-P. This plateau requires that the progenitor had an extended H envelope like a red supergiant, consistent with the slow progenitor wind speed indicated by narrow emission. The cool supergiant progenitor is significant because PTF11iqb showed WR features in its early spectrum --- meaning that the presence of such WR features in an early SN spectrum does not necessarily indicate a WR-like progenitor. [abridged] Overall, PTF11iqb bridges SNe~IIn with weaker pre-SN mass loss seen in SNe II-L and II-P, implying a continuum between these types.

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.

SN 2011A: A low-luminosity interacting transient with a double plateau and strong sodium absorpton

We present optical photometry and spectroscopy of the optical transient SN 2011A. Our data span 140 days after discovery including BVRI u′g′r′i′z′ photometry and 11 epochs of optical spectroscopy. Originally classified as a type IIn supernova (SN IIn) due to the presence of narrow Hα emission, this object shows exceptional characteristics. First, the light curve shows a double plateau, a property only observed before in the impostor SN 1997bs. Second, SN 2011A has a very low luminosity (M_V=-15.72), placing it between normal luminous SNe IIn and SN impostors. Third, SN 2011A shows low velocity and high equivalent width absorption close to the sodium doublet, which increases with time and is most likely of circumstellar origin. This evolution is also accompanied by a change in line profile; when the absorption becomes stronger, a P Cygni profile appears. We discuss SN 2011A in the context of interacting SNe IIn and SN impostors, which appears to confirm the uniqueness of this transient. While we favor an impostor origin for SN 2011A, we highlight the difficulty in differentiating between terminal and non-terminal interacting transients.

A search for circumstellar material around B-type stars in the Galactic halo

19 B-type stars, selected from the Palomar-Green Survey, have been observed at infrared wavelengths to search for possible infrared excesses, as part of an ongoing programme to investigate the nature of early-type stars at high Galactic latitudes. The resulting infrared fluxes, along with Stromgren photometry, are compared with theoretical flux profiles to determine whether any of the targets show evidence of circumstellar material, which may be indicative of post-main- sequence evolution. Eighteen of the targets have flux distributions in good agreement with theoretical predictions. However, one star, PG 2120 + 062, shows a small near-infrared excess, which may be due either to a cool companion of spectral type F5-F7, or to circumstellar material, indicating that it may be an evolved object such as a post-asymptotic giant branch star, in the transition region between the asymptotic giant branch and planetary nebula phase, with the infrared excess due to recent mass loss during giant branch evolution.

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.

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 x 10(-4) M-circle dot 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 similar to 8 x 10(-3) M-circle dot. 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.

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 la event with some resemblance to SN 1991bg.

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

Chemistry and Distribution of Daughter Species in the Circumstellar Envelopes of O-Rich AGB Stars

Context. Thanks to the advent of Herschel and ALMA, new high-quality observations of molecules present in the circumstellar envelopes of asymptotic giant branch (AGB) stars are being reported that reveal large differences from the existing chemical models. New molecular data and more comprehensive models of the chemistry in circumstellar envelopes are now available.
Aims: The aims are to determine and study the important formation and destruction pathways in the envelopes of O-rich AGB stars and to provide more reliable predictions of abundances, column densities, and radial distributions for potentially detectable species with physical conditions applicable to the envelope surrounding IK Tau.
Methods: We use a large gas-phase chemical model of an AGB envelope including the effects of CO and N2 self-shielding in a spherical geometry and a newly compiled list of inner-circumstellar envelope parent species derived from detailed modeling and observations. We trace the dominant chemistry in the expanding envelope and investigate the chemistry as a probe for the physics of the AGB phase by studying variations of abundances with mass-loss rates and expansion velocities.
Results: We find a pattern of daughter molecules forming from the photodissociation products of parent species with contributions from ion-neutral abstraction and dissociative recombination. The chemistry in the outer zones differs from that in traditional PDRs in that photoionization of daughter species plays a significant role. With the proper treatment of self-shielding, the N → N2 and C+→ CO transitions are shifted outward by factors of 7 and 2, respectively, compared with earlier models. An upper limit on the abundance of CH4 as a parent species of (≲2.5 × 10-6 with respect to H2) is found for IK Tau, and several potentially observable molecules with relatively simple chemical links to other parent species are determined. The assumed stellar mass-loss rate, in particular, has an impact on the calculated abundances of cations and the peak-abundance radius of both cations and neutrals: as the mass-loss rate increases, the peak abundance of cations generally decreases and the peak-abundance radius of all species moves outwards. The effects of varying the envelope expansion velocity and cosmic-ray ionization rate are not as significant.

The VLT-FLAMES tarantula survey:XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations - in terms of spectral type, luminosity class, and spatial location - in the field of view. Methods. We measured projected rotational velocities, 3e sin i, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(3e), of the equatorial rotational velocity, 3e. Results. The distribution of 3e sin i shows a two-component structure: a peak around 80 km s1 and a high-velocity tail extending up to 600 km s^-1 This structure is also present in the inferred distribution P(3e) with around 80% of the sample having 0 <3e ≤ 300 km s^-1 and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the highvelocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts. © 2013 Author(s).

Chemistry in AGB stars: successes and challenges

Emission and absorption line observations of molecules in late-type stars are a vital component in our understanding of stellar evolution, dust formation and mass loss in these objects. The molecular composition of the gas in the circumstellar envelopes of AGB stars reflects chemical processes in gas whose properties are strong functions of radius with density and temperature varying by more than ten and two orders of magnitude, respectively. In addition, the interstellar UV field plays a critical role in determining not only molecular abundances but also their radial distributions. In this article, I shall briefly review some recent successful approaches to describing chemistry in both the inner and outer envelopes and outline areas of challenge for the future.

Understanding B-type supergiants in the low metallicity environment of the SMC

Spectroscopic analyses of 7 SMC B-type supergiants and 1 giant have been undertaken using high resolution optical data obtained on the VLT with UVES. FASTWIND, a non-LTE, spherical, line-blanketed model atmosphere code was used to derive atmospheric and wind parameters of these stars as well as their absolute abundances. Mass-loss rates, derived from H-alpha profiles, are in poor agreement with metallicity dependent theoretical predictions. Indeed the wind-momenta of the SMC stars appear to be in good agreement with the wind-momentum luminosity relationship (WLR) of Galactic B-type stars, a puzzling result given that line-driven wind theory predicts a metallicity dependence. However the galactic stars were analysed using unblanketed model atmospheres which may mask any dependence on metallicity. A mean nitrogen enhancement of a factor of 14 is observed in the supergiants whilst only an enrichment of a factor of 4 is present in the giant, AV216. Similar excesses in nitrogen are observed in O-type dwarfs and supergiants in the same mass range, suggesting that the additional nitrogen is produced while the stars are still on the main-sequence. These nitrogen enrichments can be reproduced by current stellar evolution models, which include rotationally induced mixing, only if large initial rotational velocities of 300 kin s(-1) are invoked. Such large rotational velocities appear to be inconsistent with observed v sin i distributions for O-type stars and B-type supergiants. Hence it is suggested that the currently available stellar evolution models require more efficient mixing for lower rotational velocities.