Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore of high importance that their behaviour is properly accounted for in synthetic model populations. Here we present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations, of stars with 9-500 Msun and with metallicities ranging from Solar down to 1/250 Solar. The stellar models were computed with the 'Bonn' evolutionary code (covering core-hydrogen- and core-helium-burning phases, both complete). Post-processing for publication has been done using optimized methods, developed by our team, for massive and very massive stars. Interpolation and population synthesis were also performed on the models by our newly developed routine synStars. Eight of the grids represent slowly rotating massive stars with normal/classical evolution, while one grid represents fast rotating, chemically-homogeneously evolving models. Apart from the common stellar parameters such as mass, radius, surface temperature, luminosity and mass loss rate, we present stellar wind properties such as estimated wind velocity and kinetic energy of the wind. Additionally, we provide complete chemical yields of 34 isotopes, and estimates for the masses of the compact object remnants. The 'Bonn' Optimized Stellar Tracks (BoOST) project is published as simple tables - including stellar models, interpolated tracks and synthetic populations - thus ideal for further scientific applications. For example, star-formation studies could be done with BoOST to cover broad metallicity ranges, and so could be simulations of high-redshift galaxies. Additionally, gravitational-wave event rate predictions could be refined using BoOST by accounting for very massive stars at low-metallicity.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \hspace{1mm} \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO's measured rate of ${53.2}_{-28.2}^{+55.8} \hspace{1mm} \text{Gpc}^{-3}\text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \hspace{1mm} \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \hspace{1mm} \textrm{M}_{\odot}$ above the PISN gap. We further investigate the effects of momentum kicks during black hole formation, calculate cosmological and magnitude limited PISN rates and investigate the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations.

Type Ia supernovae (SNe Ia) are manifestations of stars deficient of hydrogen and helium disrupting in a thermonuclear runaway. While explosions of carbon-oxygen white dwarfs are thought to account for the majority of events, part of the observed diversity may be due to varied progenitor channels. We demonstrate that helium stars with masses between $\sim$1.8 and 2.5 M$_{\odot}$ may evolve into highly degenerate, near-Chandrasekhar mass cores with helium-free envelopes that subsequently ignite carbon and oxygen explosively at densities $\sim(1.8-5.9)\times 10^{9}$g cm$^{-3}$. This happens either due to core growth from shell burning (when the core has a hybrid CO/NeO composition), or following ignition of residual carbon triggered by exothermic electron captures on $^{24}$Mg (for a NeOMg-dominated composition). We argue that the resulting thermonuclear runaways is likely to prevent core collapse, leading to the complete disruption of the star. The available nuclear energy at the onset of explosive oxygen burning suffices to create ejecta with a kinetic energy of $\sim$10$^{51}$ erg, as in typical SNe Ia. Conversely, if these runaways result in partial disruptions, the corresponding transients would resemble SN Iax events similar to SN 2002cx. If helium stars in this mass range indeed explode as SNe Ia, then the frequency of events would be comparable to the observed SN Ib/c rates, thereby sufficing to account for the majority of SNe Ia in star-forming galaxies.

Star clusters are the building blocks of galaxies. They are composed of stars of nearly equal age and chemical composition, allowing us to use them as chronometers and as testbeds for gauging stellar evolution. It has become clear recently that massive stars are formed preferentially in close binaries, in which mass transfer will drastically change the evolution of the stars. This is expected to leave a significant imprint in the distribution of cluster stars in the Hertzsprung-Russell diagram. Our results, based on a dense model grid of more than 50,000 detailed binary-evolution calculations, indeed show several distinct, coeval main-sequence (MS) components, most notably an extended MS turnoff region, and a group of near-critical rotating stars that is spread over a large luminosity range on the red side of the classical MS. We comprehensively demonstrate the time evolution of the features in an animation, and we derive analytic expressions to describe these features. We find quantitative agreement with results based on recent photometric and spectroscopic observations. We conclude that while other factors may also be at play, binary evolution has a major impact on the MS morphology of young star clusters.

We report the discovery of a handful of optical hydrogen-poor knots in the central part of an extended infrared nebula centred on the [WO1] star WR 72, obtained by spectroscopic and imaging observations with the Southern African Large Telescope (SALT). Wide-field Infrared Survey Explorer (WISE) images of the nebula show that it is composed of an extended almost circular halo (of $\approx6$ arcmin or $\approx2.4$ pc in diameter) and an elongated and apparently bipolar inner shell (of a factor of six smaller size), within which the knots are concentrated. Our findings indicate that WR 72 is a new member of the rare group of hydrogen-poor planetary nebulae, which may be explained through a very late thermal pulse of a post-AGB star, or by a merger of two white dwarfs.

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars.

Be stars are rapidly rotating B main sequence stars, which show line emission due to an outflowing disc. By studying the evolution of rotating single star models, we can assess their contribution to the observed Be star populations. We identify the main effects which are responsible for single stars to approach critical rotation as functions of initial mass and metallicity, and predict the properties of populations of rotating single stars. We perform population synthesis with single star models of initial masses ranging between 3 and 30 solar masses, initial equatorial rotation velocities between 0 and 600 kms$^{-1}$ at compositions representing the Milky Way, Large and Small Magellanic Clouds. These models include efficient core-envelope coupling mediated by internal magnetic fields and correspond to the maximum efficiency of Be star production. We predict Be star fractions and the positions of fast rotating stars in the colour-magnitude diagram. We identify stellar wind mass-loss and the convective core mass fraction as the key parameters which determine the time dependance of the stellar rotation rates. Using empirical distributions of initial rotational velocities,our single star models can reproduce the trends observed in Be star fractions with mass and metallicity. However,they fail to produce a significant number of stars rotating very close to critical. We also find that rapidly rotating Be stars in the Magellanic Clouds should have significant surface nitrogen enrichments, which may be in conflict with abundance determinations of Be stars. Single star evolution may explain the high number of Be stars if 70 to 80% of critical rotationwould be sufficient to produce the Be phenomenon. However even in this case, the unexplained presence of many Be stars far below the cluster turn-off indicates the importance of the binary channel for Be star production.

$\gamma$ Cas stars are a $\sim$1% minority among classical Be stars with hard but only moderately strong continuous thermal X-ray flux and mostly very early-B spectral type. The X-ray flux has been suggested to originate from matter accelerated via magnetic disk-star interaction, by a rapidly rotating neutron star (NS) companion via the propeller effect, or by accretion onto a white dwarf (WD) companion. In view of the growing number of identified $\gamma$ Cas stars and the only imperfect matches between these suggestions and the observations, alternative models should be pursued. Two of the three best-observed $\gamma$ Cas stars, $\gamma$ Cas itself and $\pi$ Aqr, have a low-mass companion with low optical flux; interferometry of BZ Cru is inconclusive. Binary-evolution models are examined for their ability to produce such systems. The OB+He-star stage of post-mass transfer binaries, which is otherwise observationally unaccounted, can potentially reproduce many observed properties of $\gamma$ Cas stars. The interaction of the fast wind of helium stars with the disk and/or with the wind of Be stars may give rise to the production of hard X-rays. While not modelling this process, it is shown that the energy budget is favourable, and that the wind velocities may lead to hard X-rays as observed in $\gamma$ Cas stars. Furthermore, their observed number appears to be consistent with the evolutionary models. Within the Be+He-star binary model, the Be stars in $\gamma$ Cas stars are conventional classical Be stars. They are encompassed by O-star+Wolf-Rayet systems towards higher mass, where no stable Be decretion disks exist, and by Be+sdO systems at lower mass where the sdO winds may be too weak to cause the $\gamma$ Cas phenomenon. In decreasing order of the helium-star mass, the descendants could be Be+black-hole, Be+NS or Be+WD binaries.

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

We report the discovery of a spiral-like nebula with the Wide-field Infrared Survey Explorer (WISE) and the results of optical spectroscopy of its associated star TYC 8606-2025-1 with the Southern African Large Telescope (SALT). We find that TYC 8606-2025-1 is a G8 III star of $\approx3$ M$_\odot$, showing a carbon depletion by a factor of two and a nitrogen enhancement by a factor of three. We also derived an excess of s-process elements, most strongly for barium, which is a factor of three overabundant, indicating that TYC 8606-2025-1 is a mild barium star. We thereby add a new member to the small group of barium stars with circumstellar nebulae. Our radial velocity measurements indicate that TYC 8606-2025-1 has an unseen binary companion. The advanced evolutionary stage of TYC 8606-2025-1, together with the presence of a circumstellar nebula, implies an initial mass of the companion of also about 3 M$_\odot$. We conclude that the infrared nebula, due to its spiral shape, and because it has no optical counterpart, was ejected by the companion as a consequence of a very late thermal pulse, during about one orbital rotation.

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.

Gravitational wave emission can lead to the coalescence of close pairs of compact objects orbiting each other. For the case of neutron stars such mergers may yield masses above the Tolman-Oppenheimer-Volkoff limit, leading to the formation of black holes. For the case of white dwarfs the merger product may exceed the Chandrasekhar limit, leading either to a thermonuclear explosion as Type Ia supernova, or to a collapse forming a neutron star. If a Type Ia supernova explosion is avoided, the merger of two massive white dwarfs is expected to form a hydrogen- and helium-free circumstellar nebula with a hot and luminous, rapidly rotating and highly magnetized central star for several 10,000 yr before its final collapse. Here we report the discovery of a hot star with an emission line dominated spectrum in the centre of a circular mid-infrared nebula. Both the star and the nebula appear to be free of hydrogen and helium. Our tailored stellar atmosphere and wind models indicate a stellar surface temperature of about 200,000 K, and a record outflow velocity of 16,000 km/s. This extreme velocity, together with the derived mass outflow rate, imply rapid stellar rotation and a strong magnetic field aiding the wind acceleration. The Gaia distance of the star leads to a luminosity of 10^{4.5} Lsun, which matches models of the post-merger evolution of super-Chandrasekhar mass white dwarfs. The high stellar temperature and the nebular size argue for a short remaining lifetime of the star, which will produce a bright optical and high-energy transient upon collapse. Our observations indicate that super-Chandrasekhar mass white dwarf mergers can indeed avoid a thermonuclear explosion as Type Ia supernova, and provide empirical evidence for magnetic field generation in stellar mergers.

Helium accretion induced explosions in CO white dwarfs (WDs) are considered promising candidates for a number of observed types of stellar transients, including supernovae (SNe) of Type Ia and Type Iax. However, a clear favorite outcome has not yet emerged. We explore the conditions of helium ignition in the white dwarf and the final fates of helium star-WD binaries as function of their initial orbital periods and component masses. We compute 274 model binary systems with the Binary Evolution Code (BEC), where both components are fully resolved. Stellar and orbital evolution is computed simultaneously, including mass and angular momentum transfer, tides, and gravitational wave emission, as well as differential rotation and internal hydrodynamic and magnetic angular momentum transport. We find that helium detonations are expected only in systems with the shortest initial orbital periods, and for initially massive white dwarfs (MWD > 1.0 MSun ) and lower mass donors (Mdonor < 0.8 MSun), with accumulated helium layers mostly exceeding 0.1 MSun. Upon detonation, these systems would release the donor as a hypervelocity pre-WD runaway star, for which we predict the expected range of kinematic and stellar properties. Systems with more massive donors or initial periods exceeding 1.5 h will likely undergo helium deflagrations after accumulating 0.1 - 0.001 MSun of helium. Helium ignition in the white dwarf is avoided in systems with helium donor stars below - 0.6 MSun, and lead to three distinctly different groups of double white dwarf systems. The size of the parameter space open to helium detonation corresponds to only about 3 % of the galactic SN Ia rate, and to 10 % of the SN Iax rate, while the predicted large amounts of helium (>0.1 MSun) in progenitors cannot easily be reconciled with observations of archetypical SN Ia. ...

Context. The evolution of massive stars is strongly influenced by internal mixing processes such as semiconvection, convective core overshooting, and rotationally induced mixing. None of these is currently well constrained. Aims. We investigate models for massive stars in the Small Magellanic Cloud (SMC). We aim to constrain the various mixing efficiencies by comparing model results to observations. Methods. We use the stellar evolution code MESA to compute more than 60 grids of detailed evolutionary models for stars with initial masses of 9 to 100 Msol, in each grid assuming different combinations of mixing efficiencies. Results. We find that for most of the combinations of the mixing efficiencies, models in a wide mass range spend core-helium burning either only as blue supergiants, or only as red supergiants. The latter case corresponds to models that maintain a shallow slope of the hydrogen/helium (H/He) gradient separating the core and the envelope of the models. Only a small part of the mixing parameter space leads to models that produce a significant number of blue and red supergiants, which both exist abundantly in the SMC. Interestingly, these models contain steep H/He gradients, as is required to understand the hot, hydrogen-rich Wolf-Rayet stars in the SMC. We find that unless it is very fast, rotation has a limited effect on the H/He profiles in our models. Conclusions. While we use specific implementations of the considered mixing processes, they comprehensively probe the two first order structural parameters: the core mass and the H/He gradient in the core-envelope interface. Our results imply that, in massive stars, the H/He gradients above the helium cores become very steep. Our model grids can be tested with future observational surveys of the massive stars in the SMC, and thereby help to considerably reduce the uncertainties in models of massive star evolution.

We investigate how the proximity of supergiant donor stars to the Eddington-limit, and their advanced evolutionary stage, may influence the evolution of massive and ultra-luminous X-ray binaries with supergiant donor stars (SGXBs and ULXs). We construct models of massive stars with different internal hydrogen/helium gradients and different hydrogen-rich envelope masses, and expose them to slow mass loss to probe the response of the stellar radius. In addition, we compute the corresponding Roche-lobe overflow mass-transfer evolution with our detailed binary stellar evolution code, approximating the compact objects as point masses. We find that a hydrogen/helium gradient in the layers beneath the surface, as it is likely present in the well-studied donor stars of observed SGBXs, can enable nuclear timescale mass-transfer in SGXBs with a BH or a NS accretor, even for mass ratios in excess of 20. In our binary evolution models, the donor stars rapidly decrease their thermal equilibrium radius and can therefore cope with the inevitably strong orbital contraction imposed by the high mass ratio. Our results open a new perspective for understanding the large number of Galactic SGXBs, and their almost complete absence in the SMC. They may also offer a way to obtain more ULX systems, to find nuclear timescale mass-transfer in ULX systems even with neutron star accretors, and shed new light on the origin of the strong B-field in these neutron stars.

The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of $\log (M_{\rm dyn}/M_{\odot})=$3.8$\pm$0.3 for Hodge 301, and an upper limit of $\log (M_{\rm dyn}/M_{\odot})$<3.1$\pm$0.8 for SL 639, assuming Virial equilibrium. Analysis of the multi-epoch data reveals one RV-variable, potential binary candidate (VFTS744), which is likely a semi-regular variable asymptotic giant branch star. We estimate an upper limit on the observed binary fraction for our sample of 0.3, where we are sensitive to maximum periods for individual objects in the range of 1 to 10 000 days and mass-ratios above 0.3 depending on the data quality. From simulations of the RV measurements from binary systems given the current data we conclude that systems within the parameter range q>0.3, $\log$P[days]<3.5, would be detected by our variability criteria, at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters where orbital periods are uniformly distributed in the 3.3<$\log$P[days]<4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. [Abridged]

We estimate physical parameters for the late-type massive stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud (LMC). The observational sample comprises 20 candidate red supergiants (RSGs) which are the reddest (($B-V$) $>$ 1 mag) and brightest ($V$ $<$ 16 mag) objects in the VFTS. We use optical and near-IR photometry to estimate their temperatures and luminosities, and introduce the luminosity-age diagram to estimate their ages. We derive physical parameters for our targets, including temperatures from a new calibration of $(J-K_{\rm s})_{0}$ colour for luminous cool stars in the LMC, luminosities from their $J$-band magnitudes (thence radii), and ages from comparisons with state-of-the-art evolutionary models. We show that interstellar extinction is a significant factor for our targets, highlighting the need to take it into account in analysis of the physical parameters of RSGs. We find that some of the candidate RSGs could be massive AGB stars. The apparent ages of the RSGs in the Hodge 301 and SL 639 clusters show a significant spread (12-24 Myr). We also apply our approach to the RSG population of the relatively nearby NGC 2100 cluster, finding a similarly large spread. We argue that the effects of mass-transfer in binaries may lead to more massive and luminous RSGs (which we call `red stragglers') than expected from single-star evolution, and that the true cluster ages correspond to the upper limit of the estimated RSG ages. In this way, the RSGs can serve as a new and potentially reliable age tracer in young star clusters. The corresponding analysis yields ages of 24$^{+5}_{-3}$ Myr for Hodge 301, 22$^{+6}_{-5}$ Myr for SL 639, and 23$^{+4}_{-2}$ Myr for NGC 2100.

We present a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport. Our simulations cover single star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and $12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising by the end of the simulations. Although less energetic than typical events, our models are compatible with observations of less energetic explosions of low-mass progenitors. In six of our models, the mass outflow rate already exceeds the accretion rate onto the proto-neutron star, and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback. While the proto-neutron star is still accelerated by the gravitational tug of the asymmetric ejecta, the acceleration can be extrapolated to obtain estimates for the final kick velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$ and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\, \mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical neutron star birth properties can be naturally obtained in the neutrino-driven paradigm. We find a loose correlation between the explosion energy and the kick velocity. There is no indication of spin-kick alignment, but a correlation between the kick velocity and the neutron star angular momentum, which needs to be investigated further as a potential point of tension between models and observations.

We report the results of our study of the high-velocity (\approx160 km/s) runaway O star CPD-64 2731 and its associated horseshoe-shaped nebula discovered with the Wide-field Infrared Survey Explorer. Spectroscopic observations with the Southern African Large Telescope and spectral analysis indicate that CPD-64 2731 is a fast-rotating main-sequence O5.5 star with enhanced surface nitrogen abundance. We derive a projected rotational velocity of \approx300 km/s which is extremely high for this spectral type. Its kinematic age of \approx6 Myr, assuming it was born near the Galactic plane, exceeds its age derived from single star models by a factor of two. These properties suggest that CPD-64 2731 is a rejuvenated and spun-up binary product. The geometry of the nebula and the almost central location of the star within it argue against a pure bow shock interpretation for the nebula. Instead, we suggest that the binary interaction happened recently, thereby creating the nebula, with a cavity blown by the current fast stellar wind. This inference is supported by our results of 2D numerical hydrodynamic modelling.

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of $38 \pm 17 \mathrm{km \ s^{-1}}$ (1$\sigma$ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbor other stars more massive than 150 $M_\odot$ of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date.

We present a comprehensive stellar atmosphere analysis of 329 O- and B-type stars in the Small Magellanic Cloud (SMC) from the RIOTS4 survey. Using spectroscopically derived effective temperature (Teff) and surface gravities, we find that classical Be stars appear misplaced to low Teff and high luminosity in the spectroscopic Hertzsprung-Russell diagram (sHRD). Together with the most luminous stars in our sample, the stellar masses derived from the sHRD for these objects are systematically larger than those obtained from the conventional HRD. This suggests that the well-known, spectroscopic mass-discrepancy problem may be linked to the fact that both groups of stars have outer envelopes that are nearly gravitationally unbound. The non-emission-line stars in our sample mainly appear on the main-sequence, allowing a first estimate of the terminal-age main-sequence (TAMS) in the SMC, which matches the predicted TAMS between 12 and 40$\,$M$_{\odot}$ at SMC metallicity. We further find a large underabundance of stars above $\sim 25\,$M$_{\odot}$ near the ZAMS, reminiscent of such earlier findings in the Milky Way and LMC.

The first gravitational wave detections of mergers between black holes and neutron stars represent a remarkable new regime of high-energy transient astrophysics. The signals observed with LIGO-Virgo detectors come from mergers of extreme physical objects which are the end products of stellar evolution in close binary systems. To better understand their origin and merger rates, we have performed binary population syntheses at different metallicities using the new grid-based binary population synthesis code ComBinE. Starting from newborn pairs of stars, we follow their evolution including mass loss, mass transfer and accretion, common envelopes and supernova explosions. We apply the binding energies of common envelopes based on dense grids of detailed stellar structure models, make use of improved investigations of the subsequent Case BB Roche-lobe overflow and scale supernova kicks according to the stripping of the exploding stars. We demonstrate that all the double black hole mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as well as the double neutron star merger GW170817, are accounted for in our models in the appropriate metallicity regime. Our binary interaction parameters are calibrated to match the accurately determined properties of Galactic double neutron star systems, and we discuss their masses and types of supernova origin. Using our default values for the input physics parameters, we find a double neutron star merger rate of about 3.0 Myr^-1 for Milky-Way equivalent galaxies. Our upper limit to the merger-rate density of double neutron stars is R=400 yr^-1 Gpc^-3 in the local Universe (z=0).

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

A previous spectroscopic study identified the very massive O2 III star VFTS 16 in the Tarantula Nebula as a runaway star based on its peculiar line-of-sight velocity. We use the Gaia DR2 catalog to measure the relative proper motion of VFTS 16 and nearby bright stars to test if this star might have been ejected from the central cluster, R136, via dynamical ejection. We find that the position angle and magnitude of the relative proper motion (0.338 +/- 0.046 mas/yr, or approximately 80 +\- 11 km/s) of VFTS 16 are consistent with ejection from R136 approximately 1.5 +/- 0.2 Myr ago, very soon after the cluster was formed. There is some tension with the presumed age of VFTS 16 that, from published stellar parameters, cannot be greater than 0.9 +0.3/-0.2 Myr. Older ages for this star would appear to be prohibited due to the absence of He I lines in its optical spectrum, since this sets a firm lower limit on its effective temperature. The dynamical constraints may imply an unusual evolutionary history for this object, perhaps indicating it is a merger product. Gaia DR2 also confirms that another very massive star in the Tarantula Nebula, VFTS 72 (alias BI253; O2 III-V(n)((f*)), is also a runaway on the basis of its proper motion as measured by Gaia. While its tangential proper motion (0.392 +/-0.062 mas/yr or 93 +/-15 km/s) would be consistent with dynamical ejection from R136 approximately 1 Myr ago, its position angle is discrepant with this direction at the 2$\sigma$ level. From their Gaia DR2 proper motions we conclude that the two ~100 solar mass O2 stars, VFTS 16 and VFTS72, are fast runaway stars, with space velocities of around 100 km/s relative to R136 and the local massive star population. The dynamics of VFTS16 are consistent with it having been ejected from R136, and this star therefore sets a robust lower limit on the age of the central cluster of ~1.3 Myr.

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged]

Recent studies show that more than 70% of massive stars do not evolve as effectively single stars, but as members of interacting binary systems. The evolution of these stars is thus strongly altered compared to similar but isolated objects. We investigate the occurrence of parsec-scale mid-infrared nebulae around early-type stars. If they exist over a wide range of stellar properties, one possible overarching explanation is non-conservative mass transfer in binary interactions, or stellar mergers. For ~3850 stars (all OBA stars in the Bright Star Catalogue [BSC], Be stars, BeXRBs, and Be+sdO systems), we visually inspect WISE 22 $\mu$m images. Based on nebular shape and relative position, we distinguish five categories: offset bow shocks structurally aligned with the stellar space velocity, unaligned offset bow shocks, and centered, unresolved, and not classified nebulae. In the BSC, we find that 28%, 13%, and 0.4% of all O, B, and A stars, respectively, possess associated infrared (IR) nebulae. Additionally, 34/234 Be stars, 4/72 BeXRBs, and 3/17 Be+sdO systems are associated with IR nebulae. Aligned or unaligned bow shocks result from high relative velocities between star and interstellar medium (ISM) that are dominated by the star or the ISM, respectively. About 13% of the centered nebulae could be bow shocks seen head- or tail-on. For the rest, the data disfavor explanations as remains of parental disks, supernova remnants of a previous companion, and dust production in stellar winds. The existence of centered nebulae also at high Galactic latitudes strongly limits the global risk of coincidental alignments with condensations in the ISM. Mass loss during binary evolution seems a viable mechanism for the formation of at least some of these nebulae. In total, about 29% of the IR nebulae (2% of all OBA stars in the BSC) may find their explanation in the context of binary evolution.

We present a catalog of relative proper motions for 368,787 stars in the 30 Doradus region of the Large Magellanic Cloud (LMC), based on a dedicated two-epoch survey with the Hubble Space Telescope (HST) and supplemented with proper motions from our pilot archival study. We demonstrate that a relatively short epoch difference of 3 years is sufficient to reach a $\sim$0.1 mas yr$^{-1}$ level of precision or better. A number of stars have relative proper motions exceeding a 3-sigma error threshold, representing a mixture of Milky Way denizens and 17 potential LMC runaway stars. Based upon 183 VFTS OB-stars with the best proper motions, we conclude that none of them move faster than $\sim$0.3 mas yr$^{-1}$ in each coordinate -- equivalent to $\sim$70 km s$^{-1}$. Among the remaining 351 VFTS stars with less accurate proper motions, only one candidate OB runaway can be identified. We rule out any OB star in our sample moving at a tangential velocity exceeding $\sim$120 km s$^{-1}$. The most significant result of this study is finding 10 stars over wide range of masses, which appear to be ejected from the massive star cluster R136 in the tangential plane to angular distances from $35^{\prime\prime}$ out to $407^{\prime\prime}$, equivalent to 8-98 pc. The tangential velocities of these runaways appear to be correlated with apparent magnitude, indicating a possible dependence on the stellar mass.

We model the late evolution and mass loss history of rapidly rotating Wolf-Rayet stars in the mass range $5\,\rm{M}_{\odot}\dots 100\,\rm{M}_{\odot}$. We find that quasi-chemically homogeneously evolving single stars computed with enhanced mixing retain very little or no helium and are compatible with Type\,Ic supernovae. The more efficient removal of core angular momentum and the expected smaller compact object mass in our lower mass models lead to core spins in the range suggested for magnetar driven superluminous supernovae. Our more massive models retain larger specific core angular momenta, expected for long-duration gamma-ray bursts in the collapsar scenario. Due to the absence of a significant He envelope, the rapidly increasing neutrino emission after core helium exhaustion leads to an accelerated contraction of the whole star, inducing a strong spin-up, and centrifugally driven mass loss at rates of up to $10^{-2}\,\rm{M}_{\odot}~\rm{yr^{-1}}$ in the last years to decades before core collapse. Since the angular momentum transport in our lower mass models enhances the envelope spin-up, they show the largest relative amounts of centrifugally enforced mass loss, i.e., up to 25\% of the expected ejecta mass. Our most massive models evolve into the pulsational pair-instability regime. We would thus expect signatures of interaction with a C/O-rich circumstellar medium for Type~Ic superluminous supernovae with ejecta masses below $\sim 10\,\rm{M}_{\odot}$ and for the most massive engine-driven explosions with ejecta masses above $\sim 30\,\rm{M}_{\odot}$. Signs of such interaction should be observable at early epochs of the supernova explosion, and may be related to bumps observed in the light curves of superluminous supernovae, or to the massive circumstellar CO-shell proposed for Type~Ic superluminous supernova Gaia16apd.

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances. The VFTS obtained spectroscopy for B-type 434 stars. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities <80 \kms and were not classified as supergiants. In addition, nitrogen abundances for 34 similar stars from a previous survey have been re-evaluated. Approximately 75-80% of the targets have nitrogen enhancements of less than 0.3 dex, consistent with them having experienced only small amounts of mixing. However, stars with low projected rotational velocities <40 \kms and significant nitrogen enrichments are found in both our samples and simulations imply that these cannot all be rapidly rotating objects observed near pole-on. Adopting an enhancement threshold of 0.6 dex, we observed five stars, yet stellar evolution models with rotation predict only 1.25$\pm$1.11 for random stellar viewing inclinations. The excess of such objects is estimated to be 20-30% of all stars with current rotational velocities of less than 40 \kms and 2-4% of the total single B-type sample. These estimates constitute lower limits for stars that appear inconsistent with current grids of stellar evolutionary models. Including targets with smaller nitrogen enhancements of greater than 0.2 dex implies larger percentages of targets that are inconsistent with current evolutionary models, viz. 70% of the stars with rotational velocities less than 40 \kms and 6-8% of the total single stellar population. We consider possible explanations of which the most promising would appear to be breaking due to magnetic fields or stellar mergers with subsequent magnetic braking.

Massive stars play a key role in the evolution of the Universe. Our goal is to compare observed and predicted properties of single Galactic O stars to identify and constrain uncertain physical parameters and processes in stellar evolution and atmosphere models. We used a sample of 53 objects with spectral types from O3 to O9.7. For 30 of these, we determined the main photospheric and wind parameters, using optical spectroscopy and applying the FASTWIND code. For the remaining objects, literature data, obtained by means of the CMFGEN code, were used instead. The properties of our sample were compared to published predictions based on two grids evolution models that include rotationally induced mixing. Within each luminosity class, we find a close correlation of N surface abundance and luminosity, and a stronger N enrichment in more massive and evolved O stars. Additionally, a correlation of the surface nitrogen and helium abundances is observed. The large number of nitrogen-enriched stars above ~30 solar masses argues for rotationally induced mixing as the most likely explanation. However, none of the considered models can match the observed trends correctly, especially in the high mass regime. We confirm mass discrepancy for objects in the low mass O-star regime. We conclude that the rotationally induced mixing of helium to the stellar surface is too strong in some of the models. We also suggest that present inadequacies of the models to represent the N enrichment in more massive stars with relatively slow rotation might be related to problematic efficiencies of rotational mixing. We are left with a picture in which invoking binarity and magnetic fields is required to achieve a more complete agreement of the observed surface properties of a population of massive main- sequence stars with corresponding evolutionary models.

Wolf-Rayet star's winds can be so dense and so optically thick that the photosphere appears in the highly supersonic part of the outflow, veiling the underlying subsonic part of the star, and leaving the initial acceleration of the wind inaccessible to observations. We investigate the conditions and the structure of the subsonic part of the outflow of Galactic WR stars, in particular of the WNE subclass; our focus is on the conditions at the sonic point. We compute 1D hydrodynamic stellar structure models for massive helium stars adopting outer boundaries at the sonic point. We find that the outflows of our models are accelerated to supersonic velocities by the radiative force from opacity bumps either at temperatures of the order of 200kK by the Fe opacity bump or of the order of 50kK by the HeII opacity bump. For a given mass-loss rate, the conditions in the subsonic part of the outflow are independent from the detailed physical conditions in the supersonic part. The close proximity to the Eddington limit at the sonic point allows us to construct a Sonic HR diagram, relating the sonic point temperature to the L/M ratio and the stellar mass-loss rate, thereby constraining the sonic point conditions, the subsonic structure, and the stellar wind mass-loss rates from observations. The minimum mass-loss rate necessary to have the flow accelerated to supersonic velocities by the Fe opacity bump is derived. A comparison of the observed parameters of Galactic WNE stars to this minimum mass-loss rate indicates that their winds are launched to supersonic velocities by the radiation pressure arising from the Fe-bump. Conversely, models which do not show transonic flows from the Fe opacity bump form inflated envelopes. We derive an analytic criterion for the appearance of envelope inflation in the subphotospheric layers.

The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {{\sc fastwind}} with the genetic fitting algorithm {{\sc pikaia}} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and K\"{o}hler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.]

The majority of Wolf-Rayet (WR) stars represent the stripped cores of evolved massive stars who lost most of their hydrogen envelope. In low metallicity environments, such as the Small Magellanic Cloud (SMC), stellar winds are weaker and binary interaction is expected to dominate WR-star formation. However, the WR binary fraction appears to be ~40% at any metallicity. We use the recently determined physical properties of the twelve known SMC WR stars to explore their possible formation channels through comparisons with grids of SMC models, simulated with the detailed stellar evolution code MESA. These include models of rapidly rotating single stars, which experience (partial) chemically homogeneous evolution (CHE). We find that CHE is not able to account for the majority of the SMC WR stars. However, the apparently single WN star SMC AB12 and the double WR system SMC AB5 (HD 5980) appear consistent with this channel. We also analyze core helium burning stellar models assuming constant hydrogen gradients in their envelopes. We find a dichotomy in the envelope hydrogen gradients required to explain the observed temperatures of the SMC WR stars. Shallow gradients are found for the WR stars with O star companions, consistent with binary models where mass transfer occurs early, which is in agreement with their binary properties. On the other hand, much steeper hydrogen gradients are inferred for the group of hot apparently single WR stars. Since the hydrogen profiles in post main sequence models of massive stars become steeper with time, we conclude that these stars have likely been stripped by a companion during a phase of common envelope evolution. The companions, either main sequence stars or compact objects, are expected to still be present. A corresponding search might identify the first immediate double black hole binary progenitor with masses as high as those detected in GW150914.

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35.

Anomalous surface abundances are observed in a fraction of the low-mass stars of Galactic globular clusters, that may originate from hot-hydrogen-burning products ejected by a previous generation of massive stars. We present and investigate a scenario in which the second generation of polluted low-mass stars can form in shells around cool supergiant stars within a young globular cluster. Simulations of low-Z massive stars (M$_{\rm i}\sim$ 150$-$600 M$_{\odot}$) show that both core-hydrogen-burning cool supergiants and hot ionizing stellar sources are expected to be present simulaneously in young globular clusters. Under these conditions, photoionization-confined shells form around the supergiants. We find that the shell is gravitationally unstable on a timescale that is shorter than the lifetime of the supergiant, and the Bonnor-Ebert mass of the overdense regions is low enough to allow star formation. Since the low-mass stellar generation formed in this shell is made up of the material lost from the supergiant, its composition necessarily reflects the composition of the supergiant wind. We show that the wind contains hot-hydrogen-burning products, and that the shell-stars therefore have very similar abundance anomalies that are observed in the second generation stars of globular clusters. Considering the mass-budget required for the second generation star-formation, we offer two solutions. Either a top-heavy initial mass function is needed with an index of $-$1.71..$-$2.07. Alternatively, we suggest the shell-stars to have a truncated mass distribution, and solve the mass budget problem by justifiably accounting for only a fraction of the first generation. Even without forming a photoionizaton-confined shell, the cool supergiant stars predicted at low-Z could contribute to the pollution of the interstellar medium of the cluster from which the second generation was born.

The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ~36Msun and ~29Msun. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (PoWR), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. We provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.

(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models for optically-thick winds and assess how both types of models can simultaneously fulfill their respective sonic-point conditions. Fixing the outer wind law and terminal wind velocity, we obtain unique solutions for the mass-loss rates of optically-thick, radiatively-driven winds of WR stars in the phase of core He-burning. The resulting mass-loss relations as a function of stellar parameters, agree well with previous empirical relations. Furthermore, we encounter stellar mass limits below which no continuous solutions exist. While these mass limits agree with observations of WR stars in the Galaxy, they are in conflict with observations in the LMC. While our results confirm in particular the slope of oft-used empirical mass-loss relations, they imply that only part of the observed WN population can be understood in the framework of the standard assumptions of a smooth transonic flow and compact stellar core. This means that alternative approaches, such as a clumped and inflated wind structure, or deviations from the diffusion limit at the sonic point may have to be invoked. Qualitatively, the existence of mass limits for the formation of WR-type winds may be relevant for the non-detection of low-mass WR stars in binary systems, which are believed to be progenitors of Type Ib/c supernovae. The sonic-point conditions derived in this work may provide a possibility to include optically-thick winds in stellar evolution models in a more physically motivated form than in current models.

Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $\gamma$-ray burst.

Almost since the beginning, massive stars and their resultant supernovae have played a crucial role in the Universe. These objects produce tremendous amounts of energy and new, heavy elements that enrich galaxies, encourage new stars to form and sculpt the shapes of galaxies we see today. The end of millions of years of massive star evolution and the beginning of hundreds or thousands of years of supernova evolution are separated by a matter of a few seconds, in which some of the most extreme physics found in the Universe causes the explosive and terminal disruption of the star. Key questions remain unanswered in both the studies of how massive stars evolve and the behaviour of supernovae, and it appears the solutions may not lie on just one side of the explosion or the other or in just the domain of the stellar evolution or the supernova astrophysics communities. The need to view massive star evolution and supernovae as continuous phases in a single narrative motivated the Theo Murphy international scientific meeting "Bridging the gap: from massive stars to supernovae" at Chicheley Hall in June 2016, with the specific purpose to simultaneously address the scientific connections between theoretical and observational studies of massive stars and their supernovae, through engaging astronomers from both communities.

Aims: Past observations of fast-rotating massive stars exhibiting normal nitrogen abundances at their surface have raised questions about the rotational mixing paradigm. We revisit this question thanks to a spectroscopic analysis of a sample of bright fast-rotating OB stars, with the goal of quantifying the efficiency of rotational mixing at high rotation rates. Method: Our sample consists of 40 fast rotators on the main sequence, with spectral types comprised between B0.5 and O4. We compare the abundances of some key element indicators of mixing (He, CNO) with the predictions of evolutionary models for single objects and for stars in interacting binary systems. Results: The properties of half of the sample stars can be reproduced by single evolutionary models, even in the case of probable or confirmed binaries that can therefore be true single stars in a pre-interaction configuration. The main problem for the rest of the sample is a mismatch for the [N/O] abundance ratio (we confirm the existence of fast rotators with a lack of nitrogen enrichment) and/or a high helium abundance that cannot be accounted for by models. Modifying the diffusion coefficient implemented in single- star models does not solve the problem as it cannot simultaneously reproduce the helium abundances and [N/O] abundance ratios of our targets. Since part of them actually are binaries, we also compared their chemical properties with predictions for post-mass transfer systems. We found that these models can explain the abundances measured for a majority of our targets, including some of the most helium-enriched, but fail to reproduce them in other cases. Our study thus reveals that some physical ingredients are still missing in current models.

Tight binaries of helium white dwarfs (He WDs) orbiting millisecond pulsars (MSPs) will eventually "merge" due to gravitational damping of the orbit. The outcome has been predicted to be the production of long-lived ultra-compact X-ray binaries (UCXBs), in which the WD transfers material to the accreting neutron star (NS). Here we present complete numerical computations, for the first time, of such stable mass transfer from a He WD to a NS. We have calculated a number of complete binary stellar evolution tracks, starting from pre-LMXB systems, and evolved these to detached MSP+WD systems and further on to UCXBs. The minimum orbital period is found to be as short as 5.6 minutes. We followed the subsequent widening of the systems until the donor stars become planets with a mass of ~0.005 Msun after roughly a Hubble time. Our models are able to explain the properties of observed UCXBs with high helium abundances and we can identify these sources on the ascending or descending branch in a diagram displaying mass-transfer rate vs. orbital period.

Gravitational waves from the binary black hole (BH) merger GW150914 may enlighten our understanding of ultra-luminous X-ray sources (ULXs), as BHs>30Msun can reach luminosities>4x10^39 erg s^-1 without exceeding their Eddington limit. It is then important to study variations of evolutionary channels for merging BHs, which might instead form accreting BHs and become ULXs. It was recently shown that massive binaries with mass ratios close to unity and tight orbits can undergo efficient rotational mixing and evolve chemically homogeneously, resulting in a compact BH binary. We study similar systems by computing ~120000 detailed binary models with the MESA code covering a wide range of initial parameters. For initial mass ratios M2/M1~0.1-0.4, primaries >40Msun can evolve chemically homogeneously, remaining compact and forming a BH without undergoing Roche-lobe overflow. The secondary then expands and transfers mass to the BH, initiating a ULX phase. We predict that ~1 out of 10^4 massive stars evolves this way, and that in the local universe 0.13 ULXs per Msun yr^-1 of star-formation rate are observable, with a strong preference for low-metallicities. At metallicities log Z>-3, BH masses in ULXs are limited to 60Msun due to the occurrence of pair-instability supernovae which leave no remnant, resulting in an X-ray luminosity cut-off. At lower metallicities, very massive stars can avoid exploding as pair-instability supernovae and instead form BHs with masses above 130Msun, producing a gap in the ULX luminosity distribution. After the ULX phase, neutron-star-BH binaries that merge in less than a Hubble time are produced with a formation rate <0.2 Gpc^-3 yr^-1. We expect that upcoming X-ray observatories will test these predictions, which together with additional gravitational wave detections will provide strict constraints on the origin of the most massive BHs that can be produced by stars.

When a massive star in a binary system explodes as a supernova its companion star may be polluted with heavy elements from the supernova ejecta. Such a pollution had been detected in a handful of post-supernova binaries (Gonzalez Hernandez et al. 2011), but none of them is associated with a supernova remnant. We report the discovery of a solar-type star in a close, eccentric binary system with a neutron star within the young Galactic supernova remnant RCW 86. Our discovery implies that the supernova progenitor was a moving star, which exploded near the edge of its wind bubble and lost most of its initial mass due to common-envelope evolution shortly before core collapse. We find that the solar-type star is strongly polluted with calcium and other elements, which places the explosion within the class of calcium-rich supernovae -- faint and fast transients (Filippenko et al. 2003; Kasliwal et al. 2012}, whose origin is strongly debated (Kawabata et al. 2010; Waldman et al. 2011), and provides the first observational evidence that supernovae of this type can arise from core-collapse explosions (Kawabata et al. 2010; Moriya et al. 2010).

A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$.

Based mostly on stellar models which do not include rotation, CO white dwarfs which accrete helium at rates of about $\sim 10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$ have been put forward as candidate progenitors for a number of transient astrophysical phenomena, including supernovae of Type Ia, and the peculiar and fainter Type Iax supernovae. Here we study the impact of accretion-induced spin-up including the subsequent magnetic field generation, angular momentum transport, and viscous heating on the white dwarf evolution up to the point of helium ignition. We resolve the structure of the helium accreting white dwarf models with a one dimensional Langrangian hydrodynamic code, modified to include rotational and magnetic effects. We find magnetic angular momentum transport, which leads to quasi solid-body rotation, profoundly impacts the evolution of the white dwarf models. Our rotating lower mass ($0.54$ and $0.82~\mathrm{M}_\odot$) models accrete up to 50\% more mass up to ignition compared to the non-rotating case and the opposite for our more massive models. Furthermore, we find that rotation leads to up to 10-times smaller helium ignition densities, except for the lowest adopted initial white dwarf mass. Ignition densities of order $10^6\,\mathrm{g/cm}^3$ are only found for the lowest accretion rates and for large amounts of accreted helium $\gtrsim 0.4\,\mathrm{M}_\odot$. However, correspondingly massive donor stars would transfer mass at much higher rates. We expect explosive He-shell burning to mostly occur as deflagrations and at $\dot{M}>2\cdot10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$, regardless of white dwarf mass. Our results imply that helium accretion onto CO white dwarfs at the considered rates is unlikely to lead to explosions of the CO core or to classical Type,Ia supernovae, but may instead produce events like faint and fast hydrogen-free supernovae.

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo (z = 0.256), consisting of a multi-colour lightcurve from -30 d to +70 d in the rest-frame and a series of 6 spectra from PESSTO covering -7 d to +50 d. This is among the densest spectroscopic coverage, and best-constrained rising lightcurve, for a fast-declining hydrogen-poor SLSN. The bolometric lightcurve can be reproduced with a millisecond magnetar model with ~ 4 M_sol ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected ~ 6 M_sol of CO-rich material with a kinetic energy of ~ 7 x 10^51 erg, and suggests a partially thermalised additional source of luminosity between -2 d and +22 d. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS spectrum at +300 days shows no trace of SN emission lines and we place limits on the strength of [O I] from comparisons with other Ic SNe. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with 12 + log(O/H) = 8.2 in both the R23 and N2 scales (corresponding to ~ 0.3 Z_sol). We propose that the three bright regions in the host system are interacting, which thus triggers star-formation and forms young stellar populations.

The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs.

HD 164492C is a spectroscopic triple stellar system that has been recently detected to possess a strong magnetic field. We have obtained high-resolution spectroscopic and spectropolarimetric observations over a timespan of two years and derived physical, chemical, and magnetic properties for this object. The system is formed by an eccentric close spectroscopic binary (Ca1-Ca2) with a period of 12.5 days, and a massive tertiary Cb. We calculated the orbital parameters of the close pair, reconstructed the spectra of the three components, and determined atmospheric parameters and chemical abundances by spectral synthesis. From spectropolarimetric observations, multi-epoch measurements of the longitudinal magnetic fields were obtained. The magnetic field is strongly variable on timescales of a few days, with a most probable period in the range of 1.4-1.6 days. Star Cb with Teff~25000 K is the apparently fastest rotator and the most massive star of this triple system and has anomalous chemical abundances with a marked overabundance of helium, 0.35+/-0.04 by number. We identified this star as being responsible for the observed magnetic field, although the presence of magnetic fields in the components of the Ca pair cannot be excluded. Star Ca1 with a temperature of about 24000 K presents a normal chemical pattern, while the least massive star Ca2 is a mid-B type star (Teff~15000 K) with an overabundance of silicon. The obtained stellar parameters of the system components suggest a distance of 1.5 kpc and an age of 10-15 Myr.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

We present an ongoing, systematic search for extragalactic infrared transients, dubbed SPIRITS --- SPitzer InfraRed Intensive Transients Survey. In the first year, using Spitzer/IRAC, we searched 190 nearby galaxies with cadence baselines of one month and six months. We discovered over 1958 variables and 43 transients. Here, we describe the survey design and highlight 14 unusual infrared transients with no optical counterparts to deep limits, which we refer to as SPRITEs (eSPecially Red Intermediate Luminosity Transient Events). SPRITEs are in the infrared luminosity gap between novae and supernovae, with [4.5] absolute magnitudes between -11 and -14 (Vega-mag) and [3.6]-[4.5] colors between 0.3 mag and 1.6 mag. The photometric evolution of SPRITEs is diverse, ranging from < 0.1 mag/yr to > 7 mag/yr. SPRITEs occur in star-forming galaxies. We present an in-depth study of one of them, SPIRITS 14ajc in Messier 83, which shows shock-excited molecular hydrogen emission. This shock may have been triggered by the dynamic decay of a non-hierarchical system of massive stars that led to either the formation of a binary or a proto-stellar merger.

Many Type Ic superluminous supernovae have light-curve decline rates after their luminosity peak which are close to the nuclear decay rate of 56Co, consistent with the interpretation that they are powered by 56Ni and possibly pair-instability supernovae. However, their rise times are typically shorter than those expected from pair-instability supernovae, and Type Ic superluminous supernovae are often suggested to be powered by magnetar spin-down. If magnetar spin-down is actually a major mechanism to power Type Ic superluminous supernovae, it should be able to produce decline rates similar to the 56Co decay rate rather easily. In this study, we investigate the conditions for magnetars under which their spin-down energy input can behave like the 56Ni nuclear decay energy input. We find that an initial magnetic field strength within a certain range is sufficient to keep the magnetar energy deposition within a factor of a few of the 56Co decay energy for several hundreds of days. Magnetar spin-down needs to be by almost pure dipole radiation with the braking index close to 3 to mimic 56Ni in a wide parameter range. Not only late-phase 56Co-decay-like light curves, but also rise time and peak luminosity of most 56Ni-powered light curves can be reproduced by magnetars. Bolometric light curves for more than 700 days are required to distinguish the two energy sources solely by them. We expect that more slowly-declining superluminous supernovae with short rise times should be found if they are mainly powered by magnetar spin-down.

We study supernova ejecta-companion interactions in a sample of realistic semidetached binary systems representative of Type Ia supernova progenitor binaries in a single-degenerate scenario. We model the interaction process with the help of a high-resolution hydrodynamic code assuming cylindrical symmetry. We find that the ejecta hole has a half-opening angle of 40--50$^\circ$ with the density by a factor of 2-4 lower, in good agreement with the previous studies. Quantitative differences from the past results in the amounts and kinematics of the stripped companion material and levels of contamination of the companion with the ejecta material can be explained by different model assumptions and effects due to numerical diffusion.We analyse and, for the first time, provide simulation-based estimates of the amounts and of the thermal characteristics of the shock-heated material responsible for producing a prompt, soft X-ray emission. Besides the shocked ejecta material, considered in the original model by Kasen, we also account for the stripped, shock-heated envelope material of stellar companions, which we predict partially contributes to the prompt emission. The amount of the energy deposited in the envelope is comparable to the energy stored in the ejecta. The total energy budget available for the prompt emission is by a factor of about 2-4 smaller than originally predicted by Kasen. Although the shocked envelope has a higher characteristic temperature than the shocked ejecta, the temperature estimates of the shocked material are in good agreement with the Kasen's model. The hottest shocked plasma is produced in the subgiant companion case.

We report the results of spectroscopic observations and numerical modelling of the H II region IRAS 18153-1651. Our study was motivated by the discovery of an optical arc and two main-sequence stars of spectral type B1 and B3 near the centre of IRAS 18153-1651. We interpret the arc as the edge of the wind bubble (blown by the B1 star), whose brightness is enhanced by the interaction with a photoevaporation flow from a nearby molecular cloud. This interpretation implies that we deal with a unique case of a young massive star (the most massive member of a recently formed low-mass star cluster) caught just tens of thousands of years after its stellar wind has begun to blow a bubble into the surrounding dense medium. Our two-dimensional, radiation-hydrodynamics simulations of the wind bubble and the H II region around the B1 star provide a reasonable match to observations, both in terms of morphology and absolute brightness of the optical and mid-infrared emission, and verify the young age of IRAS 18153-1651. Taken together our results strongly suggest that we have revealed the first example of a wind bubble blown by a main-sequence B star.

We investigate light-curve and spectral properties of ultra-stripped core-collapse supernovae. Ultra-stripped supernovae are the explosions of heavily stripped massive stars which lost their envelopes via binary interactions with a compact companion star. They eject only ~ 0.1 Msun and may be the main way to form double neutron-star systems which eventually merge emitting strong gravitational waves. We follow the evolution of an ultra-stripped supernova progenitor until iron core collapse and perform explosive nucleosynthesis calculations. We then synthesize light curves and spectra of ultra-stripped supernovae using the nucleosynthesis results and present their expected properties. Ultra-stripped supernovae synthesize ~ 0.01 Msun of radioactive 56Ni, and their typical peak luminosity is around 1e42 erg/s or -16 mag. Their typical rise time is 5 - 10 days. Comparing synthesized and observed spectra, we find that SN 2005ek, some of the so-called calcium-rich gap transients, and SN 2010X may be related to ultra-stripped supernovae. If these supernovae are actually ultra-stripped supernovae, their event rate is expected to be about 1 per cent of core-collapse supernovae. Comparing the double neutron-star merger rate obtained by future gravitational-wave observations and the ultra-stripped supernova rate obtained by optical transient surveys identified with our synthesized light-curve and spectral models, we will be able to judge whether ultra-stripped supernovae are actually a major contributor to the binary neutron star population and provide constraints on binary stellar evolution.

The origin and evolution of magnetism in OB stars is far from being well understood. With approximately 70 magnetic OB stars known, any new object with unusual characteristics may turn out to be a key piece of the puzzle. We report the detection of an exceptionally strong magnetic field in the He-strong B2IV star CPD-62 2124. Spectropolarimetric FORS2 and HARPSpol observations were analysed by two independent teams and procedures, concluding on a strong longitudinal magnetic field of approximately 5.2 kG. The quantitative characterisation of the stellar atmosphere yields an effective temperature of 23650$\pm$250 K, a surface gravity of 3.95$\pm$0.10 dex and a surface helium fraction of 0.35$\pm$0.02 by number. The metal composition is in agreement with the cosmic abundance standard, except for Mg, Si and S, which are slightly non-solar. The strong and broad ($\sim$300 km/s) disc-like emission displayed by the H$\alpha$ line suggests a centrifugal magnetosphere supported by the strong magnetic field. Our results imply that CPD-62 2124 is an early B-type star hosting one of the strongest magnetic fields discovered to date, and one of the most evolved He-strong stars known, with a fractional main-sequence lifetime of approximately 0.6.

Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.

The "B fields in OB stars" (BOB) collaboration is based on an ESO Large Programme, to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. In the framework of this programme, we carried out low-resolution spectropolarimetric observations of a large sample of massive stars using FORS2 installed at the ESO VLT 8-m telescope. We determined the magnetic field values with two completely independent reduction and analysis pipelines. Our in-depth study of the magnetic field measurements shows that differences between our two pipelines are usually well within 3sigma errors. From the 32 observations of 28 OB stars, we were able to monitor the magnetic fields in CPD-57 3509 and HD164492C, confirm the magnetic field in HD54879, and detect a magnetic field in CPD-62 2124. We obtain a magnetic field detection rate of 6+-3% for the full sample of 69 OB stars observed with FORS2 within the BOB programme. For the pre-selected objects with a v sin i below 60 km/s, we obtain a magnetic field detection rate of 5+-5%. We also discuss X-ray properties and multiplicity of the objects in our FORS2 sample with respect to the magnetic field detections.

We review our current understanding of the progenitors of both long and short duration gamma-ray bursts (GRBs). Constraints can be derived from multiple directions, and we use three distinct strands; i) direct observations of GRBs and their host galaxies, ii) parameters derived from modeling, both via population synthesis and direct numerical simulation and iii) our understanding of plausible analog progenitor systems observed in the local Universe. From these joint constraints, we describe the likely routes that can drive massive stars to the creation of long GRBs, and our best estimates of the scenarios that can create compact object binaries which will ultimately form short GRBs, as well as the associated rates of both long and short GRBs. We further discuss how different the progenitors may be in the case of black hole engine or millisecond-magnetar models for the production of GRBs, and how central engines may provide a unifying theme between many classes of extremely luminous transient, from luminous and super-luminous supernovae to long and short GRBs.

The recently detected gravitational wave signals (GW150914 and GW151226) of the merger event of a pair of relatively massive stellar-mass black holes (BHs) calls for an investigation of the formation of such progenitor systems in general. We analyse the common envelope (CE) stage of the "traditional" formation channel in binaries where the first-formed compact object undergoes an in-spiral inside the envelope of its evolved companion star and ejects the envelope in that process. We calculate envelope binding energies of donor stars with initial masses between 4 and 115 Msun for metallicities of Z=Zsun/2 and Z=Zsun/50, and derive minimum masses of in-spiralling objects needed to eject these envelopes. We find that CE evolution, besides from producing WD-WD and NS-NS binaries, may, in principle, also produce massive BH-BH systems with individual BH component masses up to ~50-60 Msun, in particular for donor stars evolved to giants. However, the physics of envelope ejection of massive stars remains uncertain. We discuss the applicability of the energy-budget formalism, the location of the bifurcation point, the recombination energy and the accretion energy during in-spiral as possible energy sources, and also comment on the effect of inflated helium cores. Massive stars in a wide range of metallicities and with initial masses up to at least 115 Msun may possibly shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate mass and low-mass stars with degenerate cores. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 could possibly be produced from the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain. [Abridged]

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.

In an era of large spectroscopic surveys of stars and big data, sophisticated statistical methods become more and more important in order to infer fundamental stellar parameters such as mass and age. Bayesian techniques are powerful methods because they can match all available observables simultaneously to stellar models while taking prior knowledge properly into account. However, in most cases it is assumed that observables are uncorrelated which is generally not the case. Here, we include correlations in the Bayesian code BONNSAI by incorporating the covariance matrix in the likelihood function. We derive a parametrisation of the covariance matrix that, in addition to classical uncertainties, only requires the specification of a correlation parameter that describes how observables co-vary. Our correlation parameter depends purely on the method with which observables have been determined and can be analytically derived in some cases. This approach therefore has the advantage that correlations can be accounted for even if information for them are not available in specific cases but are known in general. Because the new likelihood model is a better approximation of the data, the reliability and robustness of the inferred parameters are improved. We find that neglecting correlations biases the most likely values of inferred stellar parameters and affects the precision with which these parameters can be determined. For example, we apply our technique to massive OB stars, but emphasise that it is valid for any type of stars. For effective temperatures and surface gravities determined from atmosphere modelling, we find that masses can be underestimated on average by $0.5\sigma$ and mass uncertainties overestimated by a factor of about 2 when neglecting correlations. At the same time, the age precisions are underestimated over a wide range of stellar parameters. [abridged]

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr.

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total pressure ratio inside the envelope of stellar models and the macroturbulent velocities observed in corresponding Galactic stars. To probe this connection further, we evaluated the turbulent pressure that arises in the envelope convective zones of stellar models in the mass range 1-125 Msun based on the mixing-length theory and computed for metallicities of the Large and Small Magellanic Cloud. We find that the turbulent pressure contributions in models with these metallicities located in the hot high-luminosity part of the Hertzsprung-Russel (HR) diagram is lower than in similar models with solar metallicity, whereas the turbulent pressure in low-metallicity models populating the cool part of the HR-diagram is not reduced. Based on our models, we find that the currently available observations of hot massive stars in the Magellanic Clouds appear to support a connection between macroturbulence and the turbulent pressure in stellar envelopes. Multidimensional simulations of sub-surface convection zones and a larger number of high-quality observations are necessary to test this idea more rigorously.

A significant fraction of massive main-sequence stars show strong, large-scale magnetic fields. The origin of these fields, their lifetimes, and their role in shaping the characteristics and evolution of massive stars are currently not well understood. We compile a catalogue of 389 massive main-sequence stars, 61 of which are magnetic, and derive their fundamental parameters and ages. The two samples contain stars brighter than magnitude 9 in the V band and range in mass between 5 and 100 Msun. We find that the fractional main-sequence age distribution of all considered stars follows what is expected for a magnitude limited sample, while that of magnetic stars shows a clear decrease towards the end of the main sequence. This dearth of old magnetic stars is independent of the choice of adopted stellar evolution tracks, and appears to become more prominent when considering only the most massive stars. We show that the decreasing trend in the distribution is significantly stronger than expected from magnetic flux conservation. We also find that binary rejuvenation and magnetic suppression of core convection are unlikely to be responsible for the observed lack of older magnetic massive stars, and conclude that its most probable cause is the decay of the magnetic field, over a time span longer than the stellar lifetime for the lowest considered masses, and shorter for the highest masses. We then investigate the spin-down ages of the slowly rotating magnetic massive stars and find them to exceed the stellar ages by far in many cases. The high fraction of very slowly rotating magnetic stars thus provides an independent argument for a decay of the magnetic fields.

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.

We present new binary stellar evolution models that include the effects of tidal forces, rotation, and magnetic torques with the goal of testing Planetary Nebulae (PNe) shaping via binary interaction. We explore whether tidal interaction with a companion can spin up the AGB envelope. To do so we have selected binary systems with main sequence masses of 2.5 \Mo and of 0.8 \Mo and evolve them allowing initial separations of 5, 6, 7, and 8 AU. The binary stellar evolution models have been computed all the way to the PNe formation phase or until Roche lobe overflow (RLOF) is reached, whatever happens first. We show that with initial separations of 7 and 8 AU, the binary avoids entering into RLOF, and the AGB star reaches moderate rotational velocities at the surface ($\sim 3.5 $ and $\sim 2 $ \kms respectively) during the inter-pulse phases, but after the thermal pulses it drops to a final rotational velocity of only $\sim 0.03 $ \kms. For the closest binary separations explored, 5 and 6 AU, the AGB star reaches rotational velocities of $\sim 6 $ and $\sim 4 $ \kms respectively when the RLOF is initiated. We conclude that the detached binary models that avoid entering the RLOF phase during the AGB will not shape bipolar PNe, since the acquired angular momentum is lost via the wind during the last two thermal pulses. This study rules out tidal spin-up in non-contact binaries as a sufficient condition to form bipolar PNe.

The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the Fe opacity peak leads to convective instability, and a corresponding envelope inflation may induce pulsational instability. Here, we investigate the occurrence and consequences of both instabilities in models of Wolf-Rayet stars. We determine the convective velocities in the sub-surface convective zones to estimate the amplitude of the turbulent velocity at the base of the wind that potentially leads to the formation of small-scale wind structures, as observed in several WR stars. We also investigate the effect of mass loss on the pulsations of our models. We approximated solar metallicity WR stars by models of mass-losing helium stars, and we characterized the properties of convection in the envelope adopting the standard MLT. Our results show the occurrence of sub-surface convective regions in all studied models. Small surface velocity amplitudes are predicted for models with masses below 10Msun. For models with M>10Msun, the surface velocity amplitudes are of the order of 10km/s. Moreover we find the occurrence of pulsations for stars in the mass range 9-14Msun, while mass loss appears to stabilize the more massive WR stars. We confront our results with observationally derived line variabilities of 17 WN stars. The data suggest variability to occur for stars above 10Msun, which is increasing linearly with mass above this value, in agreement with our results. We further find some of our models to be unstable to radial pulsations, and predict local magnetic fields of the order of hundreds of Gauss in WR stars more massive than 10Msun. Our study relates the surface velocity fluctuations induced by sub-surface convection to the formation of clumping in the inner part of the wind. From this mechanism, we expect a stronger variability in more massive WR stars.

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

Type Ia supernovae (SNe Ia) have been an important tool for astronomy for quite some time; however, the nature of their progenitors remains somewhat mysterious. Recent theoretical studies indicated the possibility of producing thermonuclear detonations of carbon-oxygen white dwarfs (CO WDs) at masses less than the Chandrasekhar mass through accretion of helium-rich matter, which would, depending on mass accretion rate, mass, and initial temperature of the WD, spectrally resemble either a normal SN Ia or a peculiar one. This study aims to further resolve the state of binary systems comprised of a sub-Chandrasekhar-mass CO WD and a helium star at the point where an accretion-induced detonation occurs and constrains the part of the initial parameter space where this kind of phenomenon is possible. Preexisting data obtained through simulations of single, constantly accreting CO WDs is used as an indicator for the behavior of new binary models in which the WD is treated as a point mass and which include the non-degenerate partner as a fully resolved stellar model. We parameterize the ignition of the accumulated helium layer, changes in the WD temperature, and changes in the CO core mass depending on the mass transfer rate. The initial conditions allowing for detonation do not form a single contiguous area in the parameter space, whose shape is profoundly influenced by the behavior of the donor star. Mass loss due to Nova outbursts acts in favor of detonation. According to our criteria, about 10% of the detonations in this study can be expected to show spectra consistent with ordinary SNe Ia; the rest exhibit peculiar features.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

Approximately 10% of massive OBA main-sequence (MS) and pre-MS stars harbour strong, large-scale magnetic fields. At the same time there is a dearth of magnetic stars in close binaries. A process generating strong magnetic fields only in some stars must be responsible with the merging of pre-MS and MS stars being suggested as one such channel. Stars emerging from the coalescence of two MS stars are rejuvenated, appearing younger than they are. They can therefore be identified by comparison with reference clocks. Here we predict the rejuvenation of MS merger products over a wide range of masses and binary configurations calibrated to smoothed-particle-hydrodynamical merger models. We find that the rejuvenation is of the order of the nuclear timescale and is strongest in the lowest-mass mergers and the most evolved binary progenitors with the largest mass ratios. These predictions allow us to put constraints on the binary progenitors of merger products. We show that the magnetic stars HR 2949 and $\tau$ Sco are younger than the potential binary companion HR 2948 and the Upper-Sco association, respectively, making them promising merger candidates. We find that the age discrepancies and the potential binary progenitors of both are consistent with them being rejuvenated merger products, implying that their magnetic fields may originate from this channel. Searching for age discrepancies in magnetic stars is therefore a powerful way to explore which fraction of magnetic stars may have obtained their strong magnetic field in MS mergers and to improve our understanding of magnetism in massive stars and their remnants.

We report the detection of a magnetic field in the helium-strong star CPD-57 3509 (B2 IV), a member of the Galactic open cluster NGC3293, and characterise the star's atmospheric and fundamental parameters. Spectropolarimetric observations with FORS2 and HARPSpol are analysed using two independent approaches to quantify the magnetic field strength. A high-S/N FLAMES/GIRAFFE spectrum is analysed using a hybrid non-LTE model atmosphere technique. Comparison with stellar evolution models constrains the fundamental parameters of the star. We obtain a firm detection of a surface averaged longitudinal magnetic field with a maximum amplitude of about 1 kG. Assuming a dipolar configuration of the magnetic field, this implies a dipolar field strength larger than 3.3 kG. Moreover, the large amplitude and fast variation (within about 1 day) of the longitudinal magnetic field implies that CPD-57 3509 is spinning very fast despite its apparently slow projected rotational velocity. The star should be able to support a centrifugal magnetosphere, yet the spectrum shows no sign of magnetically confined material; in particular, emission in H{\alpha} is not observed. Apparently, the wind is either not strong enough for enough material to accumulate in the magnetosphere to become observable or, alternatively, some leakage process leads to loss of material from the magnetosphere. The quantitative spectroscopic analysis of the star yields an effective temperature and a logarithmic surface gravity of 23750+-250 K and 4.05+-0.10, respectively, and a surface helium fraction of 0.28+-0.02 by number. The surface abundances of C, N, O, Ne, S, and Ar are compatible with the cosmic abundance standard, whereas Mg, Al, Si, and Fe are depleted by about a factor of 2. This abundance pattern can be understood as the consequence of a fractionated stellar wind. CPD-57 3509 is one of the most evolved He-strong stars known.

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational $\gamma$ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the $\gamma$ Dor phenomenon, calling for more detailed theoretical modelling in this direction.

The explosion of ultra-stripped stars in close binaries may explain new discoveries of weak and fast optical transients. We have demonstrated that helium star companions to neutron stars (NSs) may evolve into naked metal cores as low as ~1.5 Msun, barely above the Chandrasekhar mass limit, by the time they explode. Here we present a new systematic investigation of the progenitor evolution leading to such ultra-stripped supernovae (SNe), in some cases yielding pre-SN envelopes of less than 0.01 Msun. We discuss the nature of these SNe (electron-capture vs iron core-collapse) and their observational light-curve properties. Ultra-stripped SNe are highly relevant for binary pulsars, as well as gravitational wave detection of merging NSs by LIGO/VIRGO, since these events are expected to produce mainly low-kick NSs in the mass range 1.10-1.80 Msun.

We propose that nebular H-alpha emission as detected in the Type Ic superluminous supernova iPTF13ehe stems from matter which is stripped from a companion star when the supernova ejecta collide with it. The temporal evolution, the line broadening, and the overall blueshift of the emission are consistent with this interpretation. We scale the nebular H-alpha luminosity predicted for Type Ia supernovae in single-degenerate systems to derive the stripped mass required to explain the H-alpha luminosity of iPTF13ehe. We find a stripped mass of 0.1 - 0.9 solar masses, assuming that the supernova luminosity is powered by radioactivity or magnetar spin down. Because a central heating source is required to excite the H-alpha emission, an interaction-powered model is not favored for iPTF13ehe if the H-alpha emission is from stripped matter. We derive a companion mass of more than 20 solar masses and a binary separation of less than about 20 companion radii based on the stripping efficiency during the collision, indicating that the supernova progenitor and the companion formed a massive close binary system. If Type Ic superluminous supernovae generally occur in massive close binary systems, the early brightening observed previously in several Type Ic superluminous supernovae may also be due to the collision with a close companion. Observations of nebular hydrogen emission in future Type Ic superluminous supernovae will enable us to test this interpretation.

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

About 10$\%$ of the massive main sequence stars have recently been found to host a strong, large scale magnetic field. Both, the origin and the evolutionary consequences of these fields are largely unknown. We argue that these fields may be sufficiently strong in the deep interior of the stars to suppress convection near the outer edge of their convective core. We performed parametrised stellar evolution calculations and assumed a reduced size of the convective core for stars in the mass range 16 M$_{\odot}$ to 28 M$_{\odot}$ from the zero age main sequence until core carbon depletion. We find that such models avoid the coolest part of the main sequence band, which is usually filled by evolutionary models that include convective core overshooting. Furthermore, our `magnetic' models populate the blue supergiant region during core helium burning, i.e., the post-main sequence gap left by ordinary single star models, and some of them end their life in a position near that of the progenitor of Supernova 1987A in the HR diagram. Further effects include a strongly reduced luminosity during the red supergiant stage, and downward shift of the limiting initial mass for white dwarf and neutron star formation.

The massive red supergiant (RSG) W26 in Westerlund 1 is one of a growing number of RSGs shown to have winds that are ionized from the outside in. The fate of this dense wind material is important for models of second generation star formation in massive star clusters. Mackey et al. (2014) showed that external photoionization can stall the wind of RSGs and accumulate mass in a dense static shell. We use 1D R-HD simulations of an externally photoionized wind to predict the Halpha and [NII] emission arising from photoionized winds both with and without a dense shell. We analyse spectra of the Halpha and [NII] emission in the environment around W26 and compare them with predicted synthetic emission. Simulations of slow winds that are decelerated into a dense shell show strongly limb-brightened line emission, with line radial velocities that are independent of the wind speed. Faster winds (>22 km/s) do not form a dense shell, have less limb-brightening, and the line radial velocity is a good tracer of the wind speed. The brightness of the [NII] and Halpha lines as a function of distance from W26 agrees reasonably well with observations when only the line flux is considered. The radial velocity disagrees, however: the brightest observed emission is blueshifted by ~25 km/s relative to the radial velocity of the star, whereas a spherically symmetric wind has the brightest emission at zero radial velocity. Our results show that the bright nebula surrounding W26 must be asymmetric; we suggest it is confined by external ram pressure from the wind of the nearby supergiant W9. We obtain a lower limit on the nitrogen abundance within the nebula of 2.35 times solar. The line ratio strongly favours photoionization over shock ionization, and so even if the observed nebula is pressure confined there should still be an ionization front and a photoionization-confined shell closer to the star.

Massive, luminous stars reaching the Eddington limit in their interiors develop very dilute, extended envelopes. This effect is called envelope inflation. If the progenitors of Type Ib/c supernovae, which are believed to be Wolf-Rayet (WR) stars, have inflated envelopes then the shock breakout signals diffuse in them and can extend their rise times significantly. We show that our inflated, hydrogen-free, WR stellar models with a radius of ~Rsun can have shock breakout signals longer than ~60 s. The puzzlingly long shock breakout signal observed in the Type Ib SN 2008D can be explained by an inflated progenitor envelope, and more such events might argue in favour of existence of inflated envelopes in general.

We report the results of spectrophotometric observations of the massive star MN18 revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. Using the optical spectrum obtained with the Southern African Large Telescope, we classify this star as B1 Ia. The evolved status of MN18 is supported by the detection of nitrogen overabundance in the nebula, which implies that it is composed of processed material ejected by the star. We analysed the spectrum of MN18 by using the code CMFGEN, obtaining a stellar effective temperature of \approx 21 kK. The star is highly reddened, E(B-V)\approx 2 mag. Adopting an absolute visual magnitude of M_V=-6.8\pm0.5 (typical of B1 supergiants), MN18 has a luminosity of log L/Lsun \approx 5.42\pm0.30, a mass-loss rate of \approx (2.8-4.5)\times10^{-7} Msun/yr, and resides at a distance of \approx 5.6^{+1.5} _{-1.2} kpc. We discuss the origin of the nebula around MN18 and compare it with similar nebulae produced by other blue supergiants in the Galaxy (Sher 25, HD 168625, [SBW2007] 1) and the Large Magellanic Cloud (Sk-69 202). The nitrogen abundances in these nebulae imply that blue supergiants can produce them from the main sequence stage up to the pre-supernova stage. We also present a K-band spectrum of the candidate luminous blue variable MN56 (encircled by a ring-like nebula) and report the discovery of an OB star at \approx 17 arcsec from MN18. The possible membership of MN18 and the OB star of the star cluster Lynga 3 is discussed.

We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

After the death of a runaway massive star, its supernova shock wave interacts with the bow shocks produced by its defunct progenitor, and may lose energy, momentum, and its spherical symmetry before expanding into the local interstellar medium (ISM). We investigate whether the initial mass and space velocity of these progenitors can be associated with asymmetric supernova remnants. We run hydrodynamical models of supernovae exploding in the pre-shaped medium of moving Galactic core-collapse progenitors. We find that bow shocks that accumulate more than about 1.5 Mo generate asymmetric remnants. The shock wave first collides with these bow shocks 160-750 yr after the supernova, and the collision lasts until 830-4900 yr. The shock wave is then located 1.35-5 pc from the center of the explosion, and it expands freely into the ISM, whereas in the opposite direction it is channelled into the region of undisturbed wind material. This applies to an initially 20 Mo progenitor moving with velocity 20 km/s and to our initially 40 Mo progenitor. These remnants generate mixing of ISM gas, stellar wind and supernova ejecta that is particularly important upstream from the center of the explosion. Their lightcurves are dominated by emission from optically-thin cooling and by X-ray emission of the shocked ISM gas. We find that these remnants are likely to be observed in the [OIII] lambda 5007 spectral line emission or in the soft energy-band of X-rays. Finally, we discuss our results in the context of observed Galactic supernova remnants such as 3C391 and the Cygnus Loop.

Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and consequences for understanding Cepheid evolution.

Within the context of the collaboration "B fields in OB stars (BOB)", we used the FORS2 low-resolution spectropolarimeter to search for a magnetic field in 50 massive stars, including two reference magnetic massive stars. Because of the many controversies of magnetic field detections obtained with the FORS instruments, we derived the magnetic field values with two completely independent reduction and analysis pipelines. We compare and discuss the results obtained from the two pipelines. We obtained a general good agreement, indicating that most of the discrepancies on magnetic field detections reported in the literature are caused by the interpretation of the significance of the results (i.e., 3-4 sigma detections considered as genuine, or not), instead of by significant differences in the derived magnetic field values. By combining our results with past FORS1 measurements of HD46328, we improve the estimate of the stellar rotation period, obtaining P = 2.17950+/-0.00009 days. For HD125823, our FORS2 measurements do not fit the available magnetic field model, based on magnetic field values obtained 30 years ago. We repeatedly detect a magnetic field for the O9.7V star HD54879, the HD164492C massive binary, and the He-rich star CPD -57 3509. We obtain a magnetic field detection rate of 6+/-4%, while by considering only the apparently slow rotators we derive a detection rate of 8+/-5%, both comparable with what was previously reported by other similar surveys. We are left with the intriguing result that, although the large majority of magnetic massive stars is rotating slowly, our detection rate is not a strong function of the stellar rotational velocity.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.

The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$\alpha$ shows a variable emission. The H$\alpha$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.

The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are however found in multiple systems. By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary sub-populations with one another as well as with that of presumed single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the massive stars spin rates. We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (\vrot) for components of 114 spectroscopic binaries in 30 Doradus. The \vrot\ values are derived from the full-width at half-maximum (FWHM) of a set of spectral lines, using a FWHM vs. \vrot\ calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. The overall \vrot\ distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at $\vrot < 200$ kms) and a shoulder at intermediate velocities ($200 < \vrot < 300$ kms). The distributions of binaries and single stars however differ in two ways. First, the main peak at $\vrot \sim$100 kms is broader and slightly shifted toward higher spin rates in the binary distribution compared to that of the presumed-single stars. Second, the \vrot distribution of primaries lacks a significant population of stars spinning faster than 300 kms while such a population is clearly present in the single star sample.

Context. Oxygen sequence Wolf-Rayet (WO) stars represent a very rare stage in the evolution of massive stars. Their spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive stars, as well as their remaining lifetime and the expected properties of their supernovae. Aims. We aim to homogeneously analyse the currently known presumed-single WO stars to obtain the key stellar and outflow properties and to constrain their evolutionary state. Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code cmfgen to model the X-Shooter spectra of the WO stars and deduce the atmospheric parameters. We calculate dedicated evolutionary models to determine the evolutionary state of the stars. Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models with helium zero-age main sequence masses of $M_{\mathrm{He, ini}} = 15-25 M_{\odot}$ that exhibit a fairly strong (on the order of a few times $10^{-5} M_{\odot} \mathrm{yr}^{-1}$), homogeneous ($f_\mathrm{c} > 0.3$) stellar wind. Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of $M_{\mathrm{ini}} = 40-60 M_{\odot}$. They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years.

Massive rotating single stars with an initial metal composition appropriate for the dwarf galaxy I Zw 18 ([Fe/H]=$-$1.7) are modelled during hydrogen burning for initial masses of 9-300 M$_{\odot}$ and rotational velocities of 0-900 km s$^{-1}$. Internal mixing processes in these models were calibrated based on an observed sample of OB-type stars in the Magellanic Clouds. Even moderately fast rotators, which may be abundant at this metallicity, are found to undergo efficient mixing induced by rotation resulting in quasi chemically-homogeneous evolution. These homogeneously-evolving models reach effective temperatures of up to 90 kK during core hydrogen burning. This, together with their moderate mass-loss rates, make them Transparent Wind Ultraviolet INtense stars (TWUIN star), and their expected numbers might explain the observed HeII ionizing photon flux in I Zw 18 and other low-metallicity HeII galaxies. Our slowly rotating stars above $\sim$80 M$_{\odot}$ evolve into late B- to M-type supergiants during core hydrogen burning, with visual magnitudes up to 19$^{\mathrm{m}}$ at the distance of I Zw 18. Both types of stars, TWUIN stars and luminous late-type supergiants, are only predicted at low metallicity. Massive star evolution at low metallicity is shown to differ qualitatively from that in metal-rich environments. Our grid can be used to interpret observations of local star-forming dwarf galaxies and high-redshift galaxies, as well as the metal-poor components of our Milky Way and its globular clusters.

The evolution of massive stars even on the main sequence is not yet well understood. Due to the steep mass-luminosity relation, massive main sequence stars become very luminous. This brings their envelopes very close to the Eddington limit. We are analysing stellar evolutionary models in which the Eddington limit is reached and exceeded, and explore the rich diversity of physical phenomena which take place in their envelopes, and investigate their observational consequences. We use the grids of detailed stellar models by Brott et al. (2011) and Koehler et al. (2015), to investigate the envelope properties of core hydrogen burning massive stars. We find that at the stellar surface, the Eddington limit is almost never reached, even for stars up to 500 Msun. When an appropriate Eddington limit is defined locally in the stellar envelope, most stars more massive than 40 Msun actually exceed this limit, in particular in the partial ionization zones of iron, helium or hydrogen. While most models adjust their structure such that the local Eddington limit is exceeded at most by a few per cent, our most extreme models do so by a factor of more than seven. We find that the local violation of the Eddington limit has severe consequences for the envelope structure, as it leads to envelope inflation, convection, density inversions and possibly to pulsations. We find that all models with luminosities higher than 4*10^5 Lsun i.e. stars above 40 Msun show inflation, with a radius increase of up to a factor of about 40. We find that the hot edge of the S Dor variability region coincides with a line beyond which our models are inflated by more than a factor of two, indicating a possible connection between S Dor variability and inflation. Furthermore, our coolest models show highly inflated envelopes with masses of up to several solar masses, and appear to be candidates to produce major LBV eruptions.

The explosion of ultra-stripped stars in close binaries can lead to ejecta masses < 0.1 M_sun and may explain some of the recent discoveries of weak and fast optical transients. In Tauris et al. (2013), it was demonstrated that helium star companions to neutron stars (NSs) may experience mass transfer and evolve into naked ~1.5 M_sun metal cores, barely above the Chandrasekhar mass limit. Here we present a systematic investigation of the progenitor evolution leading to ultra-stripped supernovae (SNe). In particular, we examine the binary parameter space leading to electron-capture (EC SNe) and iron core-collapse SNe (Fe CCSNe), respectively, and determine the amount of helium ejected with applications to their observational classification as Type Ib or Type Ic. We mainly evolve systems where the SN progenitors are helium star donors of initial mass M_He = 2.5 - 3.5 M_sun in tight binaries with orbital periods of P_orb = 0.06 - 2.0 days, and hosting an accreting NS, but we also discuss the evolution of wider systems and of both more massive and lighter - as well as single - helium stars. In some cases we are able to follow the evolution until the onset of silicon burning, just a few days prior to the SN explosion. We find that ultra-stripped SNe are possible for both EC SNe and Fe CCSNe, and that the amount of helium ejected is correlated with P_orb - the tightest systems even having donors being stripped down to envelopes of less than 0.01 M_sun. We estimate the rise time of ultra-stripped SNe to be in the range 12 hr - 8 days, and light curve decay times between 1 and 50 days. Ultra-stripped SNe may produce NSs in the mass range 1.10 - 1.80 M_sun and are highly relevant for LIGO/VIRGO since most (possibly all) merging double NS systems have evolved through this phase. Finally, we discuss the low-velocity kicks which might be imparted on these resulting NSs at birth. [Abridged]

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km/s) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(obs) = 0.25 +/- 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower fractions of 0.08 +/- 0.08 and 0.10 +/- 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f_B(true) = 0.58 +/- 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

Massive B-type stars with strong magnetic fields and fast rotation are very rare and provide a mystery for theories of both star formation and magnetic field evolution. Only two such stars, called sigma Ori E analogs, were previously known. Recently, a team involved in APOGEE, one of the Sloan Digital Sky Survey III programs, announced the discovery of two additional rigidly rotating magnetosphere stars, HD23478 and HD345439. The presence of magnetic fields in these newly discovered sigma Ori E analogs was not investigated in the past. In the framework of our ESO Large Programme, and one normal ESO programme, we carried out low-resolution FORS2 spectropolarimetric observations of HD23478 and HD345439. From the measurements using hydrogen lines, we discover a rather strong longitudinal magnetic field of the order of up to 1.5kG in HD23478, and up to 1.3kG using the entire spectrum. The analysis of HD345439 using four subsequent spectropolarimetric subexposures does not reveal the presence of a magnetic field at a significance level of 3sigma. On the other hand, the inspection of individual subexposures indicates that HD345439 may host a strong magnetic field, rapidly varying over 88 minutes. A hint at the fast rotation of HD345439 is also given by the behaviour of several metallic and He I lines in the low-resolution FORS2 spectra, showing profile variations already on such a short time scale.

Accurate determinations of stellar mass functions and ages of stellar populations are crucial to much of astrophysics. We analyse the evolution of stellar mass functions of coeval main sequence stars including all relevant aspects of single- and binary-star evolution. We show that the slope of the upper part of the mass function in a stellar cluster can be quite different to the slope of the initial mass function. Wind mass loss from massive stars leads to an accumulation of stars which is visible as a peak at the high mass end of mass functions, thereby flattening the mass function slope. Mass accretion and mergers in close binary systems create a tail of rejuvenated binary products. These blue straggler stars extend the single star mass function by up to a factor of two in mass and can appear up to ten times younger than their parent stellar cluster. Cluster ages derived from their most massive stars that are close to the turn-off may thus be significantly biased. To overcome such difficulties, we propose the use of the binary tail of stellar mass functions as an unambiguous clock to derive the cluster age because the location of the onset of the binary tail identifies the cluster turn-off mass. It is indicated by a pronounced jump in the mass function of old stellar populations and by the wind mass loss peak in young stellar populations. We further characterise the binary induced blue straggler population in star clusters in terms of their frequency, binary fraction and apparent age.

The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. We find VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select HeI lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass in excess of 25Msun for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating `single' O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

Stars close to the Eddington luminosity can have large low-density inflated envelopes. We show that the rise times of shock breakout signals from supernovae can be extended significantly if supernova progenitors have an inflated stellar envelope. If the shock breakout occurs in such inflated envelopes, the shock breakout signals diffuse in them, and their rise time can be significantly extended. Then, the rise times of the shock breakout signals are dominated by the diffusion time in the inflated envelope rather than the light-crossing time of the progenitors. We show that our inflated Wolf-Rayet star models whose radii are of the order of the solar radius can have shock breakout signals which are longer than ~100 sec. The existence of inflated envelopes in Wolf-Rayet supernova progenitors may be related to the mysterious long shock breakout signal observed in Type Ib SN 2008D. Extended shock breakout signals may provide evidence for the existence of inflated stellar envelopes and can be used to constrain the physical properties of these enigmatic structures.

We present a dense model grid with tailored input chemical composition appropriate for the Large Magellanic Cloud. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses of 70-500 Msun and with initial surface rotational velocities of 0-550 km/s, covering the core-hydrogen burning phase of evolution. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 Msun, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above 150 MZAMS, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of 56000 K, at 180 Msun, and to an evolution of stars in the mass range 50-100 Msun to the regime of luminous blue variables in the HR diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above 180 Msun. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses.

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper

Classical Cepheids, like binary stars, are laboratories for stellar evolution and Cepheids in binary systems are especially powerful ones. About one-third of Galactic Cepheids are known to have companions and Cepheids in eclipsing binary systems have recently been discovered in the Large Magellanic Cloud. However, there are no known Galactic binary Cepheids with orbital periods less than one year. We compute population synthesis models of binary Cepheids to compare to the observed period and eccentricity distributions of Galactic Cepheids as well as to the number of observed eclipsing binary Cepheids in the LMC. We find that our population synthesis models are consistent with observed binary properties of Cepheids. Furthermore, we show that binary interaction on the red giant branch prevents some red giant stars from becoming classical Cepheids. Such interactions suggest that the binary fraction of Cepheids should be significantly less than that of their main-sequence progenitors, and that almost all binary Cepheids have orbital periods longer than one year. If the Galactic Cepheid spectroscopic binary fraction is about 35%, then the spectroscopic binary fraction of their intermediate mass main sequence progenitors is about 40-45%.

In stably stratified stars, numerical magneto-hydrodynamics simulations have shown that arbitrary initial magnetic fields evolve into stable equilibrium configurations, usually containing nearly axisymmetric, linked poloidal and toroidal fields that stabilize each other. In this work, we test the hypothesis that stable stratification is a requirement for the existence of such stable equilibria. For this purpose, we follow numerically the evolution of magnetic fields in barotropic (and thus neutrally stable) stars, starting from two different types of initial conditions, namely random disordered magnetic fields, as well as linked poloidal-toroidal configurations resembling the previously found equilibria. With many trials, we always find a decay of the magnetic field over a few Alfv\'en times, never a stable equilibrium. This strongly suggests that there are no stable equilibria in barotropic stars, thus clearly invalidating the assumption of barotropic equations of state often imposed on the search of magnetic equilibria. It also supports the hypothesis that, as dissipative processes erode the stable stratification, they might destabilize previously stable magnetic field configurations, leading to their decay.

Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD44743; B1 II/III) and epsilon CMa (HD52089; B1.5 II). For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field. We determined the physical parameters of both stars and characterise their X-ray spectrum. For beta CMa, our mode identification analysis led to determining a rotation period of 13.6+/-1.2 days and of an inclination angle of the rotation axis of 57.6+/-1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of ~22 degrees and a dipolar magnetic field strength (Bd) of ~100 G (60<Bd<230 G within 1 sigma), below what is typically found for other magnetic massive stars. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Both stars are expected to have a dynamical magnetosphere. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.

Interstellar bubbles around O stars are driven by a combination of the star's wind and ionizing radiation output. The wind contribution is uncertain because the boundary between the wind and interstellar medium is difficult to observe. Mid-infrared observations (e.g., of the H II region RCW 120) show arcs of dust emission around O stars, contained well within the H II region bubble. These arcs could indicate the edge of an asymmetric stellar wind bubble, distorted by density gradients and/or stellar motion. We present two-dimensional, radiation-hydrodynamics simulations investigating the evolution of wind bubbles and H II regions around massive stars moving through a dense (n=3000 cm^{-3}), uniform medium with velocities ranging from 4 to 16 km/s. The H II region morphology is strongly affected by stellar motion, as expected, but the wind bubble is also very aspherical from birth, even for the lowest space velocity considered. Wind bubbles do not fill their H II regions (we find filling factors of 10-20%), at least for a main sequence star with mass M~30 Msun. Furthermore, even for supersonic velocities the wind bow shock does not significantly trap the ionization front. X-ray emission from the wind bubble is soft, faint, and comes mainly from the turbulent mixing layer between the wind bubble and the H II region. The wind bubble radiates <1 per cent of its energy in X-rays; it loses most of its energy by turbulent mixing with cooler photoionized gas. Comparison of the simulations with the H II region RCW 120 shows that its dynamical age is <=0.4 Myr and that stellar motion <=4 km/s is allowed, implying that the ionizing source is unlikely to be a runaway star but more likely formed in situ. The region's youth, and apparent isolation from other O or B stars, makes it very interesting for studies of massive star formation and of initial mass functions.

Recent stellar evolution models show consistently that very massive metal-free stars evolve into red supergiants shortly before they explode. We argue that the envelopes of these stars, which will form pair-instability supernovae, become pulsationally unstable and that this will lead to extreme mass-loss rates despite the tiny metal content of the envelopes. We investigate the pulsational properties of such models and derive pulsationally induced mass-loss rates, which take the damping effects of the mass loss on the pulsations selfconsistently into account. We find that the pulsations may induce mass-loss rates of ~ 1e-4 - 1e-2 Msun/yr shortly before the explosions, which may create a dense circumstellar medium. Our results show that very massive stars with dense circumstellar media may stem from a wider initial mass range than pulsational-pair instability supernovae. The extreme mass loss will cease when so much of the hydrogen-rich envelope is lost that the star becomes more compact and stops pulsating. The helium core of these stars therefore remains unaffected, and their fate as pair-instability supernovae remains unaltered. The existence of dense circumstellar media around metal-free pair-instability supernovae can make them brighter and bluer, and they may be easier to detect at high redshifts than previously expected. We argue that the mass-loss enhancement in pair-instability supernova progenitors can naturally explain some observational properties of superluminous supernovae: the energetic explosions of stars within hydrogen-rich dense circumstellar media with little 56Ni production and the lack of a hydrogen-rich envelope in pair-instability supernova candidates with large 56Ni production.

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs) which have been spun up to high rotation rates via accretion of matter from a companion star in a low-mass X-ray binary (LMXB). However, many details of this recycling scenario remain to be understood. Here we investigate binary evolution in close LMXBs to study the formation of radio MSPs with low-mass helium white dwarf companions (He WDs) in tight binaries with orbital periods P_orb = 2-9 hr. In particular, we examine: i) if such observed systems can be reproduced from theoretical modelling using standard prescriptions of orbital angular momentum losses (i.e. with respect to the nature and the strength of magnetic braking), ii) if our computations of the Roche-lobe detachments can match the observed orbital periods, and iii) if the correlation between WD mass and orbital period (M_WD, P_orb) is valid for systems with P_orb < 2 days. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in ~ 400 close LMXB systems with different initial values of donor star mass, NS mass, orbital period and the so-called gamma-index of magnetic braking. Subsequently, we followed the orbital and the interior evolution of the detached low-mass (proto) He WDs, including stages with residual shell hydrogen burning. We find that a severe fine-tuning is necessary to reproduce the observed MSPs in tight binaries with He WD companions of mass < 0.20 M_sun, which suggests that something needs to be modified or is missing in the standard input physics of LMXB modelling. We demonstrate that the theoretically calculated (M_WD, P_orb)-relation is in general also valid for systems with P_orb < 2 days, although with a large scatter in He WD masses between 0.15-0.20 M_sun. The results of the thermal evolution of the (proto) He WDs are reported in a follow-up paper (Paper II).

A large number of low-mass (< 0.20 M_sun) helium white dwarfs (He WDs) have recently been discovered. The majority of these are orbiting another WD or a millisecond pulsar (MSP) in a close binary system; a few examples are found to show pulsations or to have a main-sequence star companion. There appear to be discrepancies between the current theoretical modelling of such low-mass He WDs and a number of key observed cases, indicating that their formation scenario remains to be fully understood. Here we investigate the formation of detached proto-He WDs in close-orbit low-mass X-ray binaries (LMXBs). Our prime focus is to examine the thermal evolution and the contraction phase towards the WD cooling track and investigate how this evolution depends on the WD mass. Our calculations are then compared to the most recent observational data. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in a large number of close-orbit LMXBs. Subsequently, we followed the evolution of the detached low-mass proto-He WDs, including stages with residual shell hydrogen burning and vigorous flashes caused by unstable CNO burning. We find that the time between Roche-lobe detachment until the low-mass proto-He WD reaches the WD cooling track is typically Delta_t_proto = 0.5 - 2 Gyr, depending systematically on the WD mass and therefore on its luminosity. The lowest WD mass for developing shell flashes is ~0.21 M_sun for progenitor stars of mass M2 <= 1.5 M_sun (and ~0.18 M_sun for M2 = 1.6 M_sun). The long timescale of low-mass proto-He WD evolution can explain a number of recent observations, including some MSP systems hosting He WD companions with very low surface gravities and high effective temperatures. We find no evidence for Delta_t_proto to depend on the occurrence of flashes and thus question the suggested dichotomy in thermal evolution of proto-WDs.

The distribution of stars in the Hertzsprung-Russell diagram narrates their evolutionary history and directly assesses their properties. Placing stars in this diagram however requires the knowledge of their distances and interstellar extinctions, which are often poorly known for Galactic stars. The spectroscopic Hertzsprung-Russell diagram (sHRD) tells similar evolutionary tales, but is independent of distance and extinction measurements. Based on spectroscopically derived effective temperatures and gravities of almost 600 stars, we derive for the first time the observational distribution of Galactic massive stars in the sHRD. While biases and statistical limitations in the data prevent detailed quantitative conclusions at this time, we see several clear qualitative trends. By comparing the observational sHRD with different state-of-the-art stellar evolutionary predictions, we conclude that convective core overshooting may be mass-dependent and, at high mass ($\geq 15\,M_\odot$), stronger than previously thought. Furthermore, we find evidence for an empirical upper limit in the sHRD for stars with $T_{\rm{eff}}$ between 10000 and 32000 K and, a strikingly large number of objects below this line. This over-density may be due to inflation expanding envelopes in massive main-sequence stars near the Eddington limit.

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30\,Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{km s^{-1}}$. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300\,$\rm{km s^{-1}}$, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part -- and could potentially be completely due -- to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

We report the discovery of a circular mid-infrared shell around the emission-line star Wray 16-137 using archival data of the Spitzer Space Telescope. Follow-up optical spectroscopy of Wray 16-137 with the Southern African Large Telescope revealed a rich emission spectrum typical of the classical luminous blue variables (LBVs) like P Cygni. Subsequent spectroscopic and photometric observations showed drastic changes in the spectrum and brightness during the last three years, meaning that Wray 16-137 currently undergoes an S Dor-like outburst. Namely, we found that the star has brightened by \approx 1 mag in the V and I_c bands, while its spectrum became dominated by Fe ii lines. Taken together, our observations unambiguously show that Wray 16-137 is a new member of the family of Galactic bona fide LBVs.

Powerful telescopes equipped with multi-fibre or integral field spectrographs combined with detailed models of stellar atmospheres and automated fitting techniques allow for the analysis of large number of stars. These datasets contain a wealth of information that require new analysis techniques to bridge the gap between observations and stellar evolution models. To that end, we develop BONNSAI (BONN Stellar Astrophysics Interface), a Bayesian statistical method, that is capable of comparing all available observables simultaneously to stellar models while taking observed uncertainties and prior knowledge such as initial mass functions and distributions of stellar rotational velocities into account. BONNSAI can be used to (1) determine probability distributions of fundamental stellar parameters such as initial masses and stellar ages from complex datasets, (2) predict stellar parameters that were not yet observationally determined and (3) test stellar models to further advance our understanding of stellar evolution. An important aspect of BONNSAI is that it singles out stars that cannot be reproduced by stellar models through $\chi^{2}$ hypothesis tests and posterior predictive checks. BONNSAI can be used with any set of stellar models and currently supports massive main-sequence single star models of Milky Way and Large and Small Magellanic Cloud composition. We apply our new method to mock stars to demonstrate its functionality and capabilities. In a first application, we use BONNSAI to test the stellar models of Brott et al. (2011a) by comparing the stellar ages inferred for the primary and secondary stars of eclipsing Milky Way binaries. Ages are determined from dynamical masses and radii that are known to better than 3%. We find that the stellar models reproduce the Milky Way binaries well. BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai.

At least 5 per cent of the massive stars are moving supersonically through the interstellar medium (ISM) and are expected to produce a stellar wind bow shock. We explore how the mass loss and space velocity of massive runaway stars affect the morphology of their bow shocks. We run two-dimensional axisymmetric hydrodynamical simulations following the evolution of the circumstellar medium of these stars in the Galactic plane from the main sequence to the red supergiant phase. We find that thermal conduction is an important process governing the shape, size and structure of the bow shocks around hot stars, and that they have an optical luminosity mainly produced by forbidden lines, e.g. [OIII]. The Ha emission of the bow shocks around hot stars originates from near their contact discontinuity. The H$\alpha$ emission of bow shocks around cool stars originates from their forward shock, and is too faint to be observed for the bow shocks that we simulate. The emission of optically-thin radiation mainly comes from the shocked ISM material. All bow shock models are brighter in the infrared, i.e. the infrared is the most appropriate waveband to search for bow shocks. Our study suggests that the infrared emission comes from near the contact discontinuity for bow shocks of hot stars and from the inner region of shocked wind for bow shocks around cool stars. We predict that, in the Galactic plane, the brightest, i.e. the most easily detectable bow shocks are produced by high-mass stars moving with small space velocities.

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

New spectroscopic observations of the LBV/WR multiple system HD5980 in the Small Magellanic Cloud are used to address the question of the masses and evolutionary status of the two very luminous stars in the 19.3d eclipsing binary system. Two distinct components of the N V 4944 A line are detected in emission and their radial velocity variations are used to derive masses of 61 and 66 Mo, under the assumption that binary interaction effects on this atomic transition are negligible. We propose that this binary system is the product of quasi-chemically homogeneous evolution with little or no mass transfer. Thus, both of these binary stars may be candidates for gamma-ray burst progenitors or even pair instability supernovae. Analysis of the photospheric absorption lines belonging to the third-light object in the system confirm that it consists of an O-type star in a 96.56d eccentric orbit (e=0.82) around an unseen companion. The 5:1 period ratio and high eccentricities of the two binaries suggest that they may constitute a hierarchical quadruple system.

We discuss two important effects for the astrospheres of runaway stars: the propagation of ionizing photons far beyond the astropause, and the rapid evolution of massive stars (and their winds) near the end of their lives. Hot stars emit ionizing photons with associated photoheating that has a significant dynamical effect on their surroundings. 3D simulations show that HII regions around runaway O stars drive expanding conical shells and leave underdense wakes in the medium they pass through. For late O stars this feedback to the interstellar medium is more important than that from stellar winds. Late in life, O stars evolve to cool red supergiants more rapidly than their environment can react, producing transient circumstellar structures such as double bow shocks. This provides an explanation for the bow shock and linear bar-shaped structure observed around Betelgeuse.

The enigmatic oxygen sequence Wolf-Rayet (WO) stars represent a very late stage in massive star evolution, although their exact nature is still under debate. The spectra of most of the WO stars have never been analysed through detailed modelling with a non-local thermodynamic equilibrium expanding atmosphere code. Here we present preliminary results of the first homogeneous analysis of the (apparently) single WOs.

The most massive stars on the asymptotic giant branch (AGB), so called super-AGB stars, are thought to produce supernovae (SNe) triggered by electron captures in their degenerate O+Ne+Mg cores. Super-AGB stars are expected to have slow winds with high mass-loss rates, so their wind density is high. The explosions of super-AGB stars are therefore presumed to occur in this dense wind. We provide the first synthetic light curves (LCs) for such events by exploding realistic electron-capture supernova (ecSN) progenitors within their super-AGB winds. We find that the early LC, i.e. before the recombination wave reaches the bottom of the H-rich envelope of SN ejecta (the plateau phase), is not affected by the dense wind. However, after the plateau phase, the luminosity remains higher when the super-AGB wind is taken into account. We compare our results to the historical LC of SN 1054, the progenitor of the Crab Nebula, and show that the explosion of an ecSN within an ordinary super-AGB wind can explain the LC features. We conclude that SN 1054 could have been a Type IIn SN without any extra extreme mass loss which was previously suggested to be necessary to account for its early high luminosity. We also show that our LCs match Type IIn SNe with an early plateau phase (`Type IIn-P') and suggest that they are ecSNe within super-AGB winds. Although some ecSNe can be bright in the optical spectral range due to the large progenitor radius, their X-ray luminosity from the interaction does not necessarily get as bright as other Type IIn SNe whose optical luminosities are also powered by the interaction. Thus, we suggest that optically-bright X-ray-faint Type IIn SNe can emerge from ecSNe. Optically-faint Type IIn SNe, such as SN 2008S, can also originate from ecSNe if their H-rich envelope masses are small. Some of them can be observed as `Type IIn-b' SNe due to the small H-rich envelope mass.

The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.

The circumstellar medium around massive stars is strongly impacted by stellar winds, radiation, and explosions. We use numerical simulations of these interactions to constrain the current properties and evolutionary history of various stars by comparison with observed circumstellar structures. Two- and three-dimensional simulations of bow shocks around red supergiant stars have shown that Betelgeuse has probably only recently evolved from a blue supergiant to a red supergiant, and hence its bow shock is very young and has not yet reached a steady state. We have also for the first time investigated the magnetohydrodynamics of the photoionised H II region around the nearby runaway O star Zeta Oph. Finally, we have calculated a grid of models of bow shocks around main sequence and evolved massive stars that has general application to many observed bow shocks, and which forms the basis of future work to model the explosions of these stars into their pre-shaped circumstellar medium.

Both recent observations and stellar evolution models suggest that pair-instability supernovae (PISNe) could occur in the local Universe, at metallicities below Z_Sun/3. Previous PISN models were mostly produced at very low metallicities in the context of the early Universe. We present new PISNe models at a metallicity of Z=0.001, which are relevant for the local Universe. We take the self-consistent stellar evolutionary models of pair-instability progenitors with initial masses of 150 and 250 solar masses at metallicity of Z=0.001 by Langer et al. (2007) and follow the evolution of these models through the supernova explosions, using a hydrodynamics stellar evolution code with an extensive nuclear network including 200 isotopes. Both models explode as PISNe without leaving a compact stellar remnant. Our models produce a nucleosynthetic pattern that is generally similar to that of Population III PISN models, which is mainly characterized by the production of large amounts of alpha-elements and a strong deficiency of the odd-charged elements. However, the odd-even effect in our models is significantly weaker than that found in Population III models. The comparison with the nucleosynthetic yields from core-collapse supernovae at a similar metallicity (Z=0.002) indicates that PISNe could have strongly influenced the chemical evolution below Z=0.002, assuming a standard initial mass function. The odd-even effect is predicted to be most prominent for the intermediate mass elements between silicon and calcium. With future observations of chemical abundances in Population II stars, our result can be used to constrain the number of PISNe that occurred during the past evolution of our Galaxy.

The first soft gamma-ray repeater was discovered over three decades ago, and subsequently identified as a magnetar. However there is currently no consenus on the formation channel(s) of these objects. The presence of a magnetar in the starburst cluster Westerlund 1 implies a >40Msun progenitor, favouring formation in a binary that was disrupted at supernova. To test this hypothesis we searched for the putative pre-SN companion via a radial velocity survey to identify high-velocity runaways, with subsequent atmospheric analysis of the resultant candidate, Wd1-5. Wd1-5 appears overluminous for its spectroscopic mass and we find evidence of He- and N-enrichment, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain for a single star. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~41Msun+35Msun) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field.

The blue supergiant Sher25 is surrounded by an asymmetric, hourglass-shaped circumstellar nebula, which shows similarities to the triple-ring structure seen around SN1987A. From optical spectroscopy over six consecutive nights, we detect periodic radial velocity variations in the stellar spectrum of Sher25 with a peak-to-peak amplitude of ~12 km/s on a timescale of about 6 days, confirming the tentative detec-tion of similar variations by Hendry et al. From consideration of the amplitude and timescale 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.

We report the first-ever discovery of a Wolf-Rayet (WR) star in the Large Magellanic Cloud via detection of a circular shell with the Spitzer Space Telescope. Follow-up observations with Gemini-South resolved the central star of the shell into two components separated from each other by approx 2 arcsec (or approx 0.5 pc in projection). One of these components turns out to be a WN3 star with H and He lines both in emission and absorption (we named it BAT99 3a using the numbering system based on extending the Breysacher et al. catalogue). Spectroscopy of the second component showed that it is a B0 V star. Subsequent spectroscopic observations of BAT99 3a with the du Pont 2.5-m telescope and the Southern African Large Telescope revealed that it is a close, eccentric binary system, and that the absorption lines are associated with an O companion star. We analyzed the spectrum of the binary system using the non-LTE Potsdam Wolf-Rayet (PoWR) code, confirming that the WR component is a very hot (approx 90 kK) WN star. For this star, we derived a luminosity of log L/Lsun =5.45 and a mass-loss rate of 10^{-5.8} Msun/yr, and found that the stellar wind composition is dominated by helium with 20 per cent of hydrogen. Spectroscopy of the shell revealed an He iii region centred on BAT99 3a and having the same angular radius (approx 15 arcsec) as the shell. We thereby add a new example to a rare class of high-excitation nebulae photoionized by WR stars. Analysis of the nebular spectrum showed that the shell is composed of unprocessed material, implying that the shell was swept-up from the local interstellar medium. We discuss the physical relationship between the newly identified massive stars and their possible membership of a previously unrecognized star cluster.

We investigate the formation of dust in a stellar wind during the red-supergiant (RSG) phase of a very massive Population III star with the zero-age main sequence mass of 500 M_sun. We show that, in a carbon-rich wind with a constant velocity, carbon grains can form with a lognormal-like size distribution, and that all of the carbon available for dust formation finally condense into dust for wide ranges of the mass-loss rate ((0.1-3)x10^{-3} M_sun yr^{-1}) and wind velocity (1-100 km s^{-1}). We also find that the acceleration of the wind driven by newly formed dust suppresses the grain growth but still allows more than half of gas-phase carbon to be finally locked up in dust grains. These results indicate that at most 1.7 M_sun of carbon grains can form in total during the RSG phase of 500 M_sun Population III stars. Such a high dust yield could place very massive primordial stars as important sources of dust at the very early epoch of the universe if the initial mass function of Population III stars was top-heavy. We also briefly discuss a new formation scenario of carbon-rich ultra-metal-poor stars considering the feedback from very massive Population III stars.

So called superluminous supernovae have been recently discovered in the local Universe. It appears possible that some of them originate from stellar explosions induced by the pair instability mechanism. Recent stellar evolution models also predict pair instability supernovae (PISNe) from very massive stars at fairly high metallicities (i.e. Z~0.004). We provide supernova (SN) models and synthetic light curves for two progenitor models, a 150 Msun red-supergiant and a 250 Msun yellow-supergiant at a metallicity of Z=0.001, for which the evolution from the main sequence to collapse, and the initiation of the PISN itself, has been previously computed in a realistic and self-consistent way. We are using the radiation hydrodynamics code STELLA to describe the SN evolution of both models over a time frame of about 500 days. We describe the shock-breakout phases of both SNe which are characterized by a higher luminosity, a longer duration and a lower effective temperature than those of ordinary SNeIIP. We derive the bolometric as well as the U, B, V, R and I light curves of our PISN models, which show a long-lasting plateau phase with maxima at Mbol=-19.3 mag and -21.3 mag for our lower and higher mass model, respectively. We also describe the photospheric composition and velocity as function of time. We conclude that the light curve of the explosion of our initially 150 Msun star resembles those of relatively bright SNeIIP, whereas its photospheric velocity at early times is smaller. Its 56Ni mass of 0.04 Msun also falls well into the range found in ordinary core collapse SNe. The light curve and photospheric velocity of our 250 Msun models has a striking resemblance with that of the superluminous SN2007bi, strengthening its interpretation as PISN. We conclude that PISNe may occur more frequently in the local universe than previously assumed.

The Hertzsprung-Russell diagram is an essential diagnostic diagram for stellar structure and evolution, which has now been in use for more than 100 years. Our spectroscopic Hertzsprung-Russell (sHR) diagram shows the inverse of the flux-mean gravity versus the effective temperature. Observed stars whose spectra have been quantitatively analyzed can be entered in this diagram without the knowledge of the stellar distance or absolute brightness. Observed stars can be as conveniently compared to stellar evolution calculations in the sHR diagram as in the Hertzsprung-Russell diagram. However, at the same time, our ordinate is proportional to the stellar mass-to-luminosity ratio, which can thus be directly determined. For intermediate- and low-mass star evolution at constant mass, we show that the shape of an evolutionary track in the sHR diagram is identical to that in the Hertzsprung-Russell diagram. We also demonstrate that for hot stars, their stellar Eddington factor can be directly read off the sHR diagram. For stars near their Eddington limit, we argue that a version of the sHR diagram may be useful where the gravity is exchanged by the effective gravity. We discuss the advantages and limitations of the sHR diagram, and show that it can be fruitfully applied to Galactic stars, but also to stars with known distance, e.g., in the LMC or in galaxies beyond the Local Group.

Recent magnetic field surveys in O- and B-type stars revealed that about 10% of the core-hydrogen-burning massive stars host large-scale magnetic fields. The physical origin of these fields is highly debated. To identify and model the physical processes responsible for the generation of magnetic fields in massive stars, it is important to establish whether magnetic massive stars are found in very young star-forming regions or whether they are formed in close interacting binary systems. In the framework of our ESO Large Program, we carried out low-resolution spectropolarimetric observations with FORS2 in 2013 April of the three most massive central stars in the Trifid nebula, HD164492A, HD164492C, and HD164492D. These observations indicated a strong longitudinal magnetic field of about 500-600G in the poorly studied component HD164492C. To confirm this detection, we used HARPS in spectropolarimetric mode on two consecutive nights in 2013 June. Our HARPS observations confirmed the longitudinal magnetic field in HD164492C. Furthermore, the HARPS observations revealed that HD164492C cannot be considered as a single star as it possesses one or two companions. The spectral appearance indicates that the primary is most likely of spectral type B1-B1.5V. Since in both observing nights most spectral lines appear blended, it is currently unclear which components are magnetic. Long-term monitoring using high-resolution spectropolarimetry is necessary to separate the contribution of each component to the magnetic signal. Given the location of the system HD164492C in one of the youngest star formation regions, this system can be considered as a Rosetta Stone for our understanding of the origin of magnetic fields in massive stars.

Star clusters are known as superb tools for understanding stellar evolution. In a quest for understanding the physical origin of magnetism and chemical peculiarity in about 7% of the massive main-sequence stars, we analysed two of the ten brightest members of the ~10 Myr old Galactic open cluster NGC 2264, the early B-dwarfs HD47887 and HD47777. We find accurate rotation periods of 1.95 and 2.64 days, respectively, from MOST photometry. We obtained ESPaDOnS spectropolarimetric observations, through which we determined stellar parameters, detailed chemical surface abundances, projected rotational velocities, and the inclination angles of the rotation axis. Because we found only small (<5 km/s) radial velocity variations, most likely caused by spots, we can rule out that HD47887 and HD47777 are close binaries. Finally, using the least-squares deconvolution technique, we found that both stars possess a large-scale magnetic field with an average longitudinal field strength of about 400 G. From a simultaneous fit of the stellar parameters we determine the evolutionary masses of HD47887 and HD47777 to be 9.4+/-0.7 M0 and 7.6+/-0.5 M0. Interestingly, HD47777 shows a remarkable helium underabundance, typical of helium-weak chemically peculiar stars, while the abundances of HD47887 are normal, which might imply that diffusion is operating in the lower mass star but not in the slightly more massive one. Furthermore, we argue that the rather slow rotation, as well as the lack of nitrogen enrichment in both stars, can be consistent with both the fossil and the binary hypothesis for the origin of the magnetic field. However, the presence of two magnetic and apparently single stars near the top of the cluster mass-function may speak in favour of the latter.

We have computed new stellar evolution models that include the effects of rotation and magnetic torques under different hypothesis. The goal is to test if a single star can sustain in the envelope the rotational velocities needed for the magneto hydrodynamical (MHD) simulations to shape bipolar Planetary Nebulae (PNe) when the high mass-loss rates take place. Stellar evolution models with main sequence masses of 2.5 and 5 Mo, and initial rotational velocities of 250 km/s have been followed all the way to the PNe formation phase. We find that stellar cores have to be spun down using magnetic torques in order to reproduce the rotation rates observed for white dwarfs. During the asymptotic giant branch phase and beyond, the magnetic braking of the core has a practically null effect in increasing the rotational velocity of the envelope since the stellar angular momentum is removed efficiently by the wind. We have, as well, tested best possible case scenarios in rather non-physical contexts to give enough angular momentum to the envelope. We find that we cannot get the envelope of a single star rotating at the speeds needed by the MHD simulations to form bipolar PNe. We conclude that single stellar rotators are unlikely to be the progenitors of bipolar PNe under the current MHD model paradigm.

Because the majority of massive stars are born as members of close binary systems, populations of massive main-sequence stars contain stellar mergers and products of binary mass transfer. We simulate populations of massive stars accounting for all major binary evolution effects based on the most recent binary parameter statistics and extensively evaluate the effect of model uncertainties. Assuming constant star formation, we find that $8^{+9}_{-4}\,\%$ of a sample of early type stars to be the product of a merger resulting from a close binary system. In total we find that $30^{+10}_{-15}\,\%$ of massive main-sequence stars are the product of binary interaction. We show that the commonly adapted approach to minimize the effects of binaries on an observed sample by excluding systems detected as binaries through radial velocity campaigns can be counterproductive. Systems with significant radial velocity variations are mostly pre-interaction systems. Excluding them substantially enhances the relative incidence of mergers and binary products in the non radial velocity variable sample. This poses a challenge for testing single stellar evolutionary models. It also raises the question of whether certain peculiar classes of stars, such as magnetic O-stars, are the result of binary interaction and it emphasizes the need to further study the effect of binarity on the diagnostics that are used to derive the fundamental properties (star-formation history, initial mass function, mass to light ratio) of stellar populations nearby and at high redshift.

Various types of magnetic fields occur in stars: small scale fields, large scale fields, and internal toroidal fields. While the latter may be ubiquitous in stars due to differential rotation, small scale fields (spots) may be associated with envelop convection in all low and high mass stars. The stable large scale fields found in only about 10 per cent of intermediate mass and massive stars may be understood as a consequence of dynamical binary interaction, e.g., the merging of two stars in a binary. We relate these ideas to magnetic fields in white dwarfs and neutron stars, and to their role in core-collapse and thermonuclear supernova explosions.

We present the results of a quantitative spectroscopic analysis of the oxygen-sequence Wolf- Rayet star DR1 in the low-metallicity galaxy IC 1613. Our models suggest that the strong oxygen emission lines are the result of the high temperature of this WO3 star and do not necessarily reflect a more advanced evolutionary stage than WC stars.

In this Letter, we explore the hypothesis that the smooth appearance of bow shocks around some red supergiants (RSGs) might be caused by the ionization of their winds by external sources of radiation. Our numerical simulations of the bow shock generated by IRC-10414 (the first-ever RSG with an optically detected bow shock) show that the ionization of the wind results in its acceleration by a factor of two, which reduces the difference between the wind and space velocities of the star and makes the contact discontinuity of the bow shock stable for a range of stellar space velocities and mass-loss rates. Our best fit model reproduces the overall shape and surface brightness of the observed bow shock and suggests that the space velocity and mass-loss rate of IRC-10414 are $\approx$50 ${\rm km} \, {\rm s}^{-1}$ and $\approx$$10^{-6}$ $M_\odot \, {\rm yr}^{-1}$, respectively, and that the number density of the local ISM is $\approx$3 ${\rm cm}^{-3}$. It also shows that the bow shock emission comes mainly from the shocked stellar wind. This naturally explains the enhanced nitrogen abundance in the line-emitting material, derived from the spectroscopy of the bow shock. We found that photoionized bow shocks are $\approx$15$-$50 times brighter in optical line emission than their neutral counterparts, from which we conclude that the bow shock of IRC-10414 must be photoionized.

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$.

Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show that mass transfer from a helium star to a compact companion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M_sun of material and the formation of a low-mass neutron star. We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.

Rotation is of key importance for the evolution of hot massive stars, however, the rotational velocities of these stars are difficult to determine. Based on our own data for 31 Galactic O stars and incorporating similar data for 86 OB supergiants from the literature, we aim at investigating the properties of rotational and extra line-broadening as a function of stellar parameters and at testing model predictions about the evolution of stellar rotation. Fundamental stellar parameters were determined by means of the code FASTWIND. Projected rotational and extra broadening velocities originate from a combined Ft + GOF method. Model calculations published previously were used to estimate the initial evolutionary masses. The sample O stars with Minit > 50 Msun rotate with less that 26% of their break-up velocity, and they also lack objects with v sin i < 50 km/s. For the stars with Minit > 35 Msun on the hotter side of the bi-stability jump, the observed and predicted rotational rates agree quite well; for those on the cooler side of the jump, the measured velocities are systematically higher than the predicted ones. In general, the derived extra broadening velocities decrease toward cooler Teff, whilst for later evolutionary phases they appear, at the same v sin i, higher for high-mass stars than for low-mass ones. None of the sample stars shows extra broadening velocities higher than 110 km/s. For the majority of the more massive stars, extra broadening either dominates or is in strong competition with rotation. Conclusions: For OB stars of solar metallicity, extra broadening is important and has to be accounted for in the analysis. When appearing at or close to the zero-age main sequence, most of the single and more massive stars rotate slower than previously thought. Model predictions for the evolution of rotation in hot massive stars may need to be updated.

We report the results of optical spectroscopy of a candidate evolved massive star in the Cygnus X region, TYC 3159-6-1, revealed via detection of its curious circumstellar nebula in archival data of the Spitzer Space Telecope. We classify TYC 3159-6-1 as an O9.5-O9.7 Ib star and derive its fundamental parameters by using the stellar atmosphere code FASTWIND. The He and CNO abundances in the photosphere of TYC 3159-6-1 are consistent with the solar abundances, suggesting that the star only recently evolved off the main sequence. Proper motion and radial velocity measurements for TYC 3159-6-1 show that it is a runaway star. We propose that Dolidze 7 is its parent cluster. We discuss the origin of the nebula around TYC 3159-6-1 and suggest that it might be produced in several successive episodes of enhanced mass-loss rate (outbursts) caused by rotation of the star near the critical, \Omega-limit.

Context. Oxygen sequence Wolf-Rayet (WO) stars are thought to represent the final evolutionary stage of the most massive stars. The characteristic strong O vi emission possibly originates from an enhanced oxygen abundance in the stellar wind. Alternatively, the O vi emission can be caused by the high temperature of these stars, in which case the WO stars are the high-temperature extension of the more common carbon sequence Wolf-Rayet (WC) stars. Aims. By constraining the physical properties and evolutionary status of DR1, a WO star in the low-metallicity Local Group dwarf galaxy IC 1613 and one of only two objects of its class known in a SMC-like metallicity environment, we aim to investigate the nature of WO stars and their evolutionary connection with WC stars. Methods. We use the non-LTE atmosphere code cmfgen to model the observed spectrum of DR1 and to derive its stellar and wind parameters. We compare our values with other studies of WC and WO stars, as well as with the predictions of evolutionary models. We also model the surrounding nebula using the photo-ionization code cloudy. Results. The oxygen and carbon abundances that we obtain are comparable to values found for WC stars. The temperature and luminosity are, however, higher than those of WC stars. DR1 is embedded in the hottest known H ii region in the Local Group. The nebular properties can be consistently reproduced by cloudy models adopting DR1 as central ionizing source. Conclusions. Comparison of the abundances and temperature of DR1 with core helium-burning models show that DR1 is currently well into the second half of helium burning. If the properties of DR1 are representative for the WO class, it would imply that WO stars are the high-temperature and high-luminosity extension of the WC stars, and do not necessarily represent a later evolutionary stage.

Most runaway OB stars, like the majority of massive stars residing in their parent clusters, go through the red supergiant (RSG) phase during their lifetimes. Nonetheless, although many dozens of massive runaways were found to be associated with bow shocks, only two RSG bow-shock-producing stars, Betelgeuse and \mu Cep, are known to date. In this paper, we report the discovery of an arc-like nebula around the late M-type star IRC-10414 using the SuperCOSMOS H-alpha Survey. Our spectroscopic follow-up of IRC-10414 with the Southern African Large Telescope (SALT) showed that it is a M7 supergiant, which supports previous claims on the RSG nature of this star based on observations of its maser emission. This was reinforced by our new radio- and (sub)millimeter-wavelength molecular line observations made with the Atacama Pathfinder Experiment (APEX) 12 meter telescope and the Effelsberg 100 m radio telescope, respectively. The SALT spectrum of the nebula indicates that its emission is the result of shock excitation. This finding along with the arc-like shape of the nebula and an estimate of the space velocity of IRC-10414 (\approx 70\pm20 km/s) imply the bow shock interpretation for the nebula. Thus, IRC-10414 represents the third case of a bow-shock-producing RSG and the first one with a bow shock visible at optical wavelengths. We discuss the smooth appearance of the bow shocks around IRC-10414 and Betelgeuse and propose that one of the necessary conditions for stability of bow shocks generated by RSGs is the ionization of the stellar wind. Possible ionisation sources of the wind of IRC-10414 are proposed and discussed.

It is now believed that magnetohydrodynamic equilibria can exist in stably stratified stars due to the seminal works of Braithwaite & Spruit (2004) and Braithwaite & Nordlund (2006). What is still not known is whether magnetohydrodynamic equilibria can exist in a barotropic star, in which stable stratification is not present. It has been conjectured by Reisenegger (2009) that there will likely not exist any magnetohydrodynamical equilibria in barotropic stars. We aim to test this claim by presenting preliminary MHD simulations of barotropic stars using the three dimensional stagger code of Nordlund & Galsgaard (1995).

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{\Delta RV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

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 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. 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. 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. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

At least 25 per cent of massive stars are ejected from their parent cluster, becoming runaways or exiles, travelling with often-supersonic space velocities through the interstellar medium (ISM). Their overpressurised H II regions impart kinetic energy and momentum to the ISM, compress and/or evaporate dense clouds, and can constrain properties of both the star and the ISM. Here we present one-, two-, and (the first) three-dimensional simulations of the H II region around a massive star moving supersonically through a uniform, magnetised ISM, with properties appropriate for the nearby O star Zeta Oph. The H II region leaves an expanding overdense shell behind the star and, inside this, an underdense wake that should be filled with hot gas from the shocked stellar wind. The gas column density in the shell is strongly influenced by the ISM magnetic field strength and orientation. H-alpha emission maps show the H II region remains roughly circular, although the star is displaced somewhat from the centre of emission. For our model parameters, the kinetic energy feedback from the H II region is comparable to the mechanical luminosity of the stellar wind, and the momentum feedback rate is >100X larger than that from the wind and about 10X larger than the total momentum input rate available from radiation pressure. Compared to the star's eventual supernova explosion, the kinetic energy feedback from the H II region over the star's main sequence lifetime is >100X less, but the momentum feedback is up to 4X larger. H II region dynamics are found to have only a small effect on the ISM conditions that a bow shock close to the star would encounter.

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs), formed via type Ib/c core-collapse supernovae (SNe), which have been spun up to high rotation rates via accretion from a companion star in a low-mass X-ray binary (LMXB). In an alternative formation channel, NSs are produced via the accretion-induced collapse (AIC) of a massive white dwarf (WD) in a close binary. Here we investigate binary evolution leading to AIC and examine if NSs formed in this way can subsequently be recycled to form MSPs and, if so, how they can observationally be distinguished from pulsars formed via the standard core-collapse SN channel in terms of their masses, spins, orbital periods and space velocities. Numerical calculations with a detailed stellar evolution code were used for the first time to study the combined pre- and post-AIC evolution of close binaries. We investigated the mass transfer onto a massive WD in 240 systems with three different types of non-degenerate donor stars: main-sequence stars, red giants, and helium stars. When the WD is able to accrete sufficient mass (depending on the mass-transfer rate and the duration of the accretion phase) we assumed it collapses to form a NS and we studied the dynamical effects of this implosion on the binary orbit. Subsequently, we followed the mass-transfer epoch which resumes once the donor star refills its Roche lobe and calculated the continued LMXB evolution until the end. We demonstrate that the final properties of these MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. Finally, we discuss the link between AIC and young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions. [Abridged]

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few x 0.001 Msun), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.

The nearby red supergiant (RSG) Betelgeuse has a complex circumstellar medium out to at least 0.5 parsecs from its surface, shaped by its mass-loss history within the past 0.1 Myr, its environment, and its motion through the interstellar medium (ISM). In principle its mass-loss history can be constrained by comparing hydrodynamic models with observations. Observations and numerical simulations indicate that Betelgeuse has a very young bow shock, hence the star may have only recently become a RSG. To test this possibility we calculated a stellar evolution model for a single star with properties consistent with Betelgeuse. We incorporated the resulting evolving stellar wind into 2D hydrodynamic simulations to model a runaway blue supergiant (BSG) undergoing the transition to a RSG near the end of its life. The collapsing BSG wind bubble induces a bow shock-shaped inner shell which at least superficially resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We investigate whether this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

Globular clusters contain many stars with surface abundance patterns indicating contributions from hydrogen burning products, as seen in the anti-correlated elemental abundances of e.g. sodium and oxygen, and magnesium and aluminium. Multiple generations of stars can explain this phenomenon, with the second generation forming from a mixture of pristine gas and ejecta from the first generation. We show that massive binary stars may be a source of much of the material that makes this second generation of stars. Mass transfer in binaries is often non-conservative and the ejected matter moves slowly enough that it can remain inside a globular cluster and remain available for subsequent star formation. Recent studies show that there are more short-period massive binaries than previously thought, hence also more stars that interact and eject nuclear-processed material.

Theories on the origin of magnetic fields in massive stars remain poorly developed, because the properties of their magnetic field as function of stellar parameters could not yet be investigated. To investigate whether magnetic fields in massive stars are ubiquitous or appear only in stars with a specific spectral classification, certain ages, or in a special environment, we acquired 67 new spectropolarimetric observations for 30 massive stars. Among the observed sample, roughly one third of the stars are probable members of clusters at different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. Spectropolarimetric observations were obtained during four different nights using the low-resolution spectropolarimetric mode of FORS2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of the VLT. Furthermore, we present a number of follow-up observations carried out with the high-resolution spectropolarimeters SOFIN mounted at the Nordic Optical Telescope (NOT) and HARPS mounted at the ESO 3.6m between 2008 and 2011. To assess the membership in open clusters and associations, we used astrometric catalogues with the highest quality kinematic and photometric data currently available. The presence of a magnetic field is confirmed in nine stars previously observed with FORS1/2: HD36879, HD47839, CPD-282561, CPD-472963, HD93843, HD148937, HD149757, HD328856, and HD164794. New magnetic field detections at a significance level of at least 3sigma were achieved in five stars: HD92206c, HD93521, HD93632, CPD-468221, and HD157857. Among the stars with a detected magnetic field, five stars belong to open clusters with high membership probability. According to previous kinematic studies, five magnetic O-type stars in our sample are candidate runaway stars.

Rotation is thought to be a major factor in the evolution of massive stars, especially at low metallicity, with consequences for their chemical yields, ionizing flux and final fate. Determining the natal rotation-rate distribution of stars is of high priority given its importance as a constraint on theories of massive star formation and as input for models of stellar populations in the local Universe and at high redshift. Recently, it has become clear that the majority of massive stars interact with a binary companion before they die. We investigate how this affects the distribution of rotation rates. For this purpose, we simulate a massive binary-star population typical for our Galaxy assuming continuous star formation. We find that, because of binary interaction, 20^+5_-10% of all massive main-sequence stars have projected rotational velocities in excess of 200km/s. We evaluate the effect of uncertain input distributions and physical processes and conclude that the main uncertainties are the mass transfer efficiency and the possible effect of magnetic braking, especially if magnetic fields are generated or amplified during mass accretion and stellar mergers. The fraction of rapid rotators we derive is similar to that observed. If indeed mass transfer and mergers are the main cause for rapid rotation in massive stars, little room remains for rapidly rotating stars that are born single. This implies that spin down during star formation is even more efficient than previously thought. In addition, this raises questions about the interpretation of the surface abundances of rapidly rotating stars as evidence for rotational mixing. Furthermore, our results allow for the possibility that all early-type Be stars result from binary interactions and suggest that evidence for rotation in explosions, such as long gamma-ray bursts, points to a binary origin.

Massive stars are essential to understand a variety of branches of astronomy including galaxy and star cluster evolution, nucleosynthesis and supernovae, pulsars and black holes. It has become evident that massive star evolution is very diverse, being sensitive to metallicity, binarity, rotation, and possibly magnetic fields. While the problem to obtain a good statistical observational database is alleviated by current large spectroscopic surveys, it remains a challenge to model these diverse paths of massive stars towards their violent end stage. We show that the main sequence stage offers the best opportunity to gauge the relevance of the various possible evolutionary scenarios. This also allows to sketch the post-main sequence evolution of massive stars, for which observations of Wolf-Rayet stars give essential clues. Recent supernova discoveries due to the current boost in transient searches allow tentative mappings of progenitor models with supernova types, including pair instability supernovae and gamma-ray bursts.

Aims: Projected rotational velocities (\vsini) 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 and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0\leq$\ve$\leq$100\,\kms and the high velocity component having \ve$\sim 250$\,\kms. 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.

Nearby companions alter the evolution of massive stars in binary systems. Using a sample of Galactic massive stars in nearby young clusters, we simultaneously measure all intrinsic binary characteristics relevant to quantify the frequency and nature of binary interactions. We find a large intrinsic binary fraction, a strong preference for short orbital periods and a flat distribution for the mass-ratios. Our results do not support the presence of a significant peak of equal-mass `twin' binaries. As a result of the measured distributions, we find that over seventy per cent of all massive stars exchange mass with a companion. Such a rate greatly exceeds previous estimates and implies that the majority of massive stars have their evolution strongly affected by interaction with a nearby companion.

Measurements of rates of period change of Classical Cepheids probe stellar physics and evolution. Additionally, better understanding of Cepheid structure and evolution provides greater insight into their use as standard candles and tools for measuring the Hubble constant. Our recent study of the period change of the nearest Cepheid, Polaris, suggested that it is undergoing enhanced mass loss when compared to canonical stellar evolution model predictions. In this work, we expand the analysis to rates of period change measured for about 200 Galactic Cepheids and compare them to population synthesis models of Cepheids including convective core overshooting and enhanced mass loss. Rates of period change predicted from stellar evolution models without mass loss do not agree with observed rates whereas including enhanced mass loss yields predicted rates in better agreement with observations. This is the first evidence that enhanced mass loss as suggested previously for Polaris and delta Cephei must be a ubiquitous property of Classical Cepheids.

(shortened) We investigate the properties of Galactic WR stars and their environment to identify evolutionary channels that may lead to the formation of LGRBs. To this purpose we compile available information on the spectropolarimetric properties, the presence of CS ejecta, and the CS velocities in the environment of Galactic WR stars. We use linear line-depolarization as an indicator of rotation, nebular morphology as an indicator of stellar ejecta, and velocity patterns in UV absorption features as an indicator of increased velocities in the CS environment. We find that the ~23% WR stars with "possible ejecta nebulae" dominate the population of WR stars with spectropolarimetric signatures of rotation, while WR stars without such nebulae only rarely show indications of rotation. The corresponding objects are most likely in an early stage after a preceding RSG or LBV phase, and have not yet lost their angular momenta due to the strong mass-loss in the WR phase. From their photometric periods we estimate rotation parameters in the range omega = 0.04...0.25, corresponding to moderate rotation speeds of 36...120 km/s. These values are very uncertain, but comply with the specific surface angular momentum requirement for LGRB progenitors. Our results indicate that, in the Galaxy, mainly "young" WR stars shortly after a RSG/LBV phase show spectropolarimetric signatures of rotation. Their rotation rates are thus likely enhanced with respect to the majority of Galactic WR stars. According to their estimated specific surface angular momenta, a subgroup of stars exploding in this phase may represent an evolutionary channel towards LGRBs at high metallicities, comparable to the Galaxy.

Although the first millisecond pulsars (MSPs) were discovered 30 years ago we still do not understand all details of their formation process. Here, we present new results from Tauris, Langer & Kramer (2012) on the recycling scenario leading to radio MSPs with helium or carbon-oxygen white dwarf companions via evolution of low- and intermediate mass X-ray binaries (LMXBs, IMXBs). We discuss the location of the spin-up line in the (P,Pdot)-diagram and estimate the amount of accreted mass needed to obtain a given spin period and compare with observations. Finally, we constrain the true ages of observed recycled pulsars via calculated isochrones in the (P,Pdot)-diagram.

(shortened) The first couple of stellar generations may have been massive, of order 100 Msun, and to have played a dominant role in galaxy formation and the chemical enrichment of the early Universe. Some fraction of these objects may have died as pair-instability supernovae or gamma-ray bursts. The winds if these stars may have played an important role in determining these outcomes. As the winds are driven by radiation pressure on spectral lines, their strengths are expected to vary with metallicity. Until now, most mass-loss predictions for metal-poor O-type stars have assumed a scaled-down solar-abundance pattern. However, Population III evolutionary tracks show significant surface enrichment through rotational mixing of CNO-processed material, because even metal-poor stars switch to CNO-burning early on. We address the question of whether the CNO surface enhanced self-enrichment in the first few generations of stars could impact their mass-loss properties. For this, we employ Monte Carlo simulations to establish the local line-force and solve for the momentum equation of the stellar outflow, testing whether an outflow can actually be established by assessing the net acceleration at the sonic point of the flow. Stellar evolution models of rotating metal-poor stars are used to specify the surface chemical composition, focussing on the phases of early enrichment. We find that the mass-loss rates of CNO enhanced metal-poor stars are higher than those of non-enriched stars, but they are much lower than those rates where the CNO abundance is included in the total abundance Z. We present a heuristic formula that provides mass-loss estimates for CNO-dominated winds in relation to scaled-down solar abundances.

Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^\pi\ (with 0.15 < log P < 3.5) and f_q ~ q^\kappa\ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are \pi = -0.45 +/- 0.30 and \kappa = -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

Mass-loss rate, $\dot{M}$, is one of the key parameters affecting evolution and observational manifestations of massive stars, and their impact on the ambient medium. Despite its importance, there is a factor of ~100 discrepancy between empirical and theoretical $\dot{M}$ of late-type O dwarfs, the so-called weak-wind problem. In this Letter, we propose a simple novel method to constrain $\dot{M}$ of runaway massive stars through observation of their bow shocks and Stromgren spheres, which might be of decisive importance for resolving the weak-wind problem. Using this method, we found that $\dot{M}$ of the well-known runaway O9.5 V star zeta Oph is more than an order of magnitude higher than that derived from ultraviolet (UV) line-fitting (Marcolino et al. 2009) and is by a factor of 6 to 7 lower than those based on the theoretical recipe by Vink et al. (2000) and the Halpha line (Mokiem et al. 2005). The discrepancy between $\dot{M}$ derived by our method and that based on UV lines would be even more severe if the stellar wind is clumpy. At the same time, our estimate of $\dot{M}$ agrees with that predicted by the moving reversing layer theory by Lucy (2010).

Millisecond pulsars (MSPs) are mainly characterised by their spin periods, B-fields and masses - quantities which are largely affected by previous interactions with a companion star in a binary system. In this paper, we investigate the formation mechanism of MSPs by considering the pulsar recycling process in both intermediate-mass X-ray binaries (IMXBs) and low-mass X-ray binaries (LMXBs). The IMXBs mainly lead to the formation of binary MSPs with a massive carbon-oxygen (CO) or an oxygen-neon-magnesium white dwarf (ONeMg WD) companion, whereas the LMXBs form recycled pulsars with a helium white dwarf (He WD) companion. We discuss the accretion physics leading to the spin-up line in the PPdot-diagram and demonstrate that such a line cannot be uniquely defined. We derive a simple expression for the amount of accreted mass needed for any given pulsar to achieve its equilibrium spin and apply this to explain the observed differences of the spin distributions of recycled pulsars with different types of companions. From numerical calculations we present further evidence for significant loss of rotational energy in accreting X-ray MSPs in LMXBs during the Roche-lobe decoupling phase (Tauris 2012) and demonstrate that the same effect is negligible in IMXBs. We examine the recycling of pulsars with CO WD companions via Case BB Roche-lobe overflow (RLO) of naked helium stars in post common envelope binaries. We find that such pulsars typically accrete of the order 0.002-0.007 M_sun which is just about sufficient to explain their observed spin periods. We introduce isochrones of radio MSPs in the PPdot-diagram to follow their spin evolution and discuss their true ages from comparison with observations. Finally, we apply our results of the spin-up process to the massive pulsar J1614-2230 (Paper I) and put new constraints on the birth masses of a number of recycled pulsars. [Abridged]

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, X-ray binaries and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. Over seventy per cent of all massive stars will exchange mass with a companion, leading to a binary merger in one third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

Rotational mixing in massive main sequence stars is predicted to monotonically increase their surface nitrogen abundance with time. We use this effect to design a method for constraining the age and the inclination angle of massive main sequence stars, given their observed luminosity, effective temperature, projected rotational velocity and surface nitrogen abundance. This method relies on stellar evolution models for different metallicities, masses and rotation rates. We use the population synthesis code STARMAKER to show the range of applicability of our method. We apply this method to 79 early B-type main sequence stars near the LMC clusters NGC 2004 and N 11 and the SMC clusters NGC 330 and NGC 346. From all stars within the sample, 17 were found to be suitable for an age analysis. For ten of them, which are rapidly rotating stars without a strong nitrogen enhancement, it has been previously concluded that they did not evolve as rotationally mixed single stars. This is confirmed by our analysis, which flags the age of these objects as highly discrepant with their isochrone ages. For the other seven stars, their nitrogen and isochrone ages are found to agree within error bars, what validates our method. Constraints on the inclination angle have been derived for the other 62 stars,with the implication that the nitrogen abundances of the SMC stars, for which mostly only upper limits are known, fall on average significantly below those limits. Nitrogen chronology is found to be a new useful tool for testing stellar evolution and to constrain fundamental properties of massive main sequence stars. A web version of this tool is provided.

A significant fraction of massive stars are moving supersonically through the interstellar medium (ISM), either due to disruption of a binary system or ejection from their parent star cluster. The interaction of their wind with the ISM produces a bow shock. In late evolutionary stages these stars may undergo rapid transitions from red to blue and vice versa on the Hertzsprung-Russell diagram, with accompanying rapid changes to their stellar winds and bow shocks. Recent 3D simulations of the bow shock produced by the nearby runaway red supergiant (RSG) Betelgeuse, under the assumption of a constant wind, indicate that the bow shock is very young (<30000 years old), hence Betelgeuse may have only recently become a RSG. To test this possibility, we have calculated stellar evolution models for single stars which match the observed properties of Betelgeuse in the RSG phase. The resulting evolving stellar wind is incorporated into 2D hydrodynamic simulations in which we model a runaway blue supergiant (BSG) as it undergoes the transition to a RSG near the end of its life. We find that the collapsing BSG wind bubble induces a bow shock-shaped inner shell around the RSG wind that resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We suggest that this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

[Abridged] We present a new grid of massive population III star models including the effects of rotation on the stellar structure and chemical mixing, and magnetic torques for the transport of angular momentum. Based on the grid, we also present a phase diagram for the expected final fates of rotating massive Pop III stars. Our non-rotating models become redder than the previous models in the literature, given the larger overshooting parameter adopted in this study. In particular, convective dredge-up of the helium core material into the hydrogen envelope is observed in our non-rotating very massive star models (>~200 Msun), which is potentially important for the chemical yields. On the other hand, the stars become bluer and more luminous with a higher rotational velocity. With the Spruit-Tayler dynamo, our models with a sufficiently high initial rotational velocity can reach the critical rotation earlier and lose more mass as a result, compared to the previous models without magnetic fields. The most dramatic effect of rotation is found with the so-called chemically homogeneous evolution (CHE), which is observed for a limited mass and rotational velocity range. CHE has several important consequences: 1) Both primary nitrogen and ionizing photons are abundantly produced. 2) Conditions for gamma-ray burst progenitors are fulfilled for an initial mass range of 13 - 84 Msun. 3) Pair instability supernovae of type Ibc are expected for 84 -190 Msun and 4) Both a pulsational pair instability supernova and a GRB may occur from the same progenitor of about 56 - 84 Msun, which might significantly influence the consequent GRB afterglow. We find that CHE does not occur for very massive stars (> 190 Msun), in which case the hydrogen envelope expands to the red-supergiant phase and the final angular momentum is too low to make any explosive event powered by rotation.

Betelgeuse, the bright, cool red supergiant in Orion, is moving supersonically relative to the local interstellar medium. The star emits a powerful stellar wind which collides with this medium, forming a cometary structure, a bow shock, pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. The models include realistic low temperature cooling and cover a range of plausible interstellar medium densities and stellar velocities between 0.3 - 1.9 cm-3 and 28 - 73 km/s. We show that the flow dynamics and morphology of the bow shock differ substantially due to the preferential growth of Rayleigh-Taylor or Kelvin-Helmholtz instabilities in the models. The former dominate the models with slow stellar velocities resulting in a clumpy bow shock sub-structure, whereas the latter produce a smoother, more layered sub-structure in the fast models. If the mass in the bow shock shell is low, as seems to be implied by the AKARI luminosities (~0.003 Msun), then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular nature of the bow shock shell is consistent with this conclusion. Thus, our results suggest that Betelgeuse entered the red supergiant phase only recently.

We report near-infrared radial velocity measurements of the recently identified donor star in the high mass X-ray binary system OAO 1657-415 obtained in the H band using ISAAC on the VLT. Cross-correlation methods were employed to construct a radial velocity curve with a semi-amplitude of $22.1 \pm 3.5$ km s$^{-1}$. Combined with other measured parameters of this system this provides a dynamically determined neutron star mass of $1.42 \pm 0.26 M$_{\odot}$ and a mass of $14.3 \pm 0.8 M$_{\odot}$ for the Ofpe/WN9 highly evolved donor star. OAO 1657-415 is an eclipsing High Mass X-ray binary pulsar with the largest eccentricity and orbital period of any within its class. Of the ten known eclipsing X-ray binary pulsars OAO 1657-415 becomes the ninth with a dynamically determined neutron star 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 High Mass X-ray binaries (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 OAO1657-415 remains problematic, we conclude by proposing two scenarios which may account for OAO1657-415 current orbital configuration.

We report the discovery of two new Galactic candidate luminous blue variable (cLBV) stars via detection of circular shells (typical of known confirmed and cLBVs) and follow-up spectroscopy of their central stars. The shells were detected at 22 um in the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars of the shells conducted with the renewed Southern African Large Telescope (SALT) showed that their spectra are very similar to those of the well-known LBVs P Cygni and AG Car, and the recently discovered cLBV MN112, which implies the LBV classification for these stars as well. The LBV classification of both stars is supported by detection of their significant photometric variability: one of them brightened in the R- and I-bands by 0.68\pm0.10 mag and 0.61\pm0.04 mag, respectively, during the last 13-18 years, while the second one (known as Hen 3-1383) varies its B,V,R,I and K_s brightnesses by \simeq 0.5-0.9 mag on time-scales from 10 days to decades. We also found significant changes in the spectrum of Hen 3-1383 on a timescale of \simeq 3 months, which provides additional support for the LBV classification of this star. Further spectrophotometric monitoring of both stars is required to firmly prove their LBV status. We discuss a connection between the location of massive stars in the field and their fast rotation, and suggest that the LBV activity of the newly discovered cLBVs might be directly related to their possible runaway status.

Polaris is one of the most observed stars in the night sky, with recorded observations spanning more than 200 years. From these observations, one can study the real-time evolution of Polaris via the secular rate of change of the pulsation period. However, the measurements of the rate of period change do not agree with predictions from state-of-the-art stellar evolution models. We show that this may imply that Polaris is currently losing mass at a rate of $\dot{M} \approx 10^{-6} M_\odot$ yr$^{-1}$ based on the difference between modeled and observed rates of period change, consistent with pulsation-enhanced Cepheid mass loss. A relation between the rate of period change and mass loss has important implications for understanding stellar evolution and pulsation, and provides insight into the current Cepheid mass discrepancy.

We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pcs from the X-ray pulsar PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer from the progenitor of PSR J0537-691 and the supernova explosion respectively.

Mass loss forms an important aspect of the evolution of massive stars, as well as for the enrichment of the surrounding ISM. Our goal is to predict accurate mass-loss rates and terminal wind velocities. These quantities can be compared to empirical values, thereby testing radiation-driven wind models. One specific issue is that of the "weak-wind problem", where empirically derived mass-loss rates fall orders of magnitude short of predicted values. We employ an established Monte Carlo model and a recently suggested new line acceleration formalism to solve the wind dynamics consistently. We provide a new grid of mass-loss rates and terminal wind velocities of O stars, and compare the values to empirical results. Our models fail to provide mass-loss rates for main-sequence stars below a luminosity of log(L/Lsun) = 5.2, where we run into a fundamental limit. At luminosities below this critical value there is insufficient momentum transferred in the region below the sonic point to kick-start the acceleration. This problem occurs at the location of the onset of the weak-wind problem. For O dwarfs, the boundary between being able to start a wind, and failing to do so, is at spectral type O6/O6.5. The direct cause of this failure is a combination of the lower luminosity and a lack of Fe V lines at the wind base. This might indicate that another mechanism is required to provide the necessary driving to initiate the wind. For stars more luminous than log(L/Lsun) = 5.2, our new mass-loss rates are in excellent agreement with the mass-loss prescription by Vink et al. 2000. This implies that the main assumption entering the method of the Vink et al. prescriptions - i.e. that the momentum equation is not explicitly solved for - does not compromise the reliability of the Vink et al. results for this part of parameter space (Abridged).

Some studies have claimed the existence of a stellar upper-mass limit of 150 Msun. A factor that is often overlooked concerns the issue that there might be a significant difference between the present-day and the initial mass of the most massive stars - as a result of mass loss. The upper-mass limit may be substantially higher, possibly exceeding 200 Msun. The issue of the upper mass-limit will however remain uncertain as long as there is only limited quantitative knowledge of mass loss in close proximity to the Eddington (= Gamma) limit. For this reason, we present mass-loss predictions from Monte Carlo radiative transfer models for very massive stars up to 300 Msun. Using our new dynamical approach, we find an upturn or "kink" in the mass-loss versus Gamma dependence, at the point where our model winds become optically thick. These are the first mass-loss predictions where the transition from optically thin O-star winds to optically thick Wolf-Rayet winds has been resolved.

Context. The Cepheid mass discrepancy, the difference between masses predicted from stellar evolution and stellar pulsation calculations, is a challenge for the understanding of stellar astrophysics. Recent models of the eclipsing binary Cepheid OGLE-LMC-CEP-0227 have suggested that the discrepancy may be resolved. Aims. We explore for what physical parameters do stellar evolution models agree with the measured properties of OGLE-LMC-CEP0227 and compare to canonical stellar evolution models assuming no convective core overshooting. Methods. We construct state-of-the-art stellar evolution models for varying mass, metallicity, and convective core overshooting and compare the stellar evolution predictions with the observed properties. Results. The observed mass, effective temperature, and radius of the two stars in the binary system are well fit by numerous combinations of physical parameters, suggesting a Cepheid mass discrepancy of 10-20% relative to canonical stellar evolution models. Conclusions. The properties of the observed binary Cepheid suggest that the Cepheid mass discrepancy is still a challenge and requires more specific observations, such as the rate of period change, to better constrain and understand the necessary physics for stellar evolution models to resolve the discrepancy.

NGC346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud discovered by the VLT-FLAMES survey of massive stars. Spectra obtained with VLT-FLAMES are used to construct a radial velocity curve and photometry obtained with the Faulkes Telescope South is then used to derive orbital parameters, while spectra of the secondary are compared with synthetic spectra from TLUSTY model atmospheres. The orbital period is found to be 4.20381(12) days, with masses of 19.1+/-1.0 and 11.9+/-0.6 Msun. The primary is a rapidly rotating late-O dwarf while the secondary, an early-B giant, displays near-synchronous rotation and has filled its Roche lobe, implying that it was originally the more massive component with recent mass transfer `spinning up' the primary to near-critical rotation. Comparison with synthetic spectra finds temperatures of 34.5kK and 24.5kK for the primary and secondary respectively, with the nitrogen abundance of the secondary enhanced compared to baseline values for the SMC, consistent with the predictions of models of interacting binaries. NGC346-013 likely evolved via non-conservative mass transfer in a system with initial masses ~22+15Msun, with the well-constrained orbital solution and atmospheric parameters making it an excellent candidate for tailored modelling with binary evolution codes. This system will form a cornerstone in constraining the physics of thermal timescale mass transfer, and the associated mass transfer efficiency, in massive close binary systems.

We report on the status of our spectropolarimetric observations of massive stars. During the last years, we have discovered magnetic fields in many objects of the upper main sequence, including Be stars, beta Cephei and Slowly Pulsating B stars, and a dozen O stars. Since the effects of those magnetic fields have been found to be substantial by recent models, we are looking into their impact on stellar rotation, pulsation, stellar winds, and chemical abundances. Accurate studies of the age, environment, and kinematic characteristics of the magnetic stars are also promising to give us new insight into the origin of the magnetic fields. Furthermore, longer time series of magnetic field measurements allow us to observe the temporal variability of the magnetic field and to deduce the stellar rotation period and the magnetic field geometry. Studies of the magnetic field in massive stars are indispensable to understand the conditions controlling the presence of those fields and their implications on the stellar physical parameters and evolution.

Massive stars lose a large fraction of their original mass over the course of their evolution. These stellar winds shape the surrounding medium according to parameters that are the result of the characteristics of the stars, varying over time as the stars evolve, leading to both permanent and temporary features that can be used to constrain the evolution of the progenitor star. Because long Gamma-Ray Bursts (GRBs) are thought to originate from massive stars, the characteristics of the circumstellar medium (CSM) should be observable in the signal of GRBs. This can occur directly, as the characteristics of the GRB-jet are changed by the medium it collides with, and indirectly because the GRB can only be observed through the extended circumstellar bubble that surrounds each massive star. We use computer simulations to describe the circumstellar features that can be found in the vicinity of massive stars and discuss if, and how, they may appear in GRB observations. Specifically, we make hydrodynamical models of the circumstellar environment of a rapidly rotating, chemically near-homogeneous star, which is represents a GRB progenitor candidate model. The simulations show that the star creates a large scale bubble of shocked wind material, which sweeps up the interstellar medium in an expanding shell. Within this bubble, temporary circumstellar shells, clumps and voids are created as a result of changes in the stellar wind. Most of these temporary features have disappeared by the time the star reaches the end of its life, leaving a highly turbulent circumstellar bubble behind. Placing the same star in a high density environment simplifies the evolution of the CSM as the more confined bubble prohibits the formation of some of the temporary structures.

The Cepheid mass discrepancy, the difference between mass estimates from stellar pulsation and stellar evolution models, is a long standing challenge for the understanding of stellar astrophysics. We discuss the current state of the mass discrepancy and test the role of pulsation-driven mass loss in state-of-the-art stellar evolution calculations of Large Magellanic Cepheids in resolving it. We find that Cepheid mass loss is a significant contributor to the mass discrepancy, but it is not clear if the metallicity dependence of Cepheid mass loss is consistent with the measured metallicity dependence of the mass discrepancy.

The most massive stars are thought to be hydrogen-rich Wolf-Rayet stars of late spectral subtype (WNh stars). In previous theoretical studies the enhanced mass loss of these stars has been attributed to their proximity to the Eddington limit. Here we investigate observed trends in the mass-loss properties of such young, very massive stars. We derive theoretical mass-luminosity relations for very massive stars, based on a large grid of stellar structure models. Using these relations, we estimate Eddington factors for a sample of stars, under different assumptions of their evolutionary status. We evaluate the resulting mass-loss relations, and compare them with theoretical predictions. We find observational evidence that the mass loss in the WR regime is dominated by the Eddington parameter Gamma_e, which has important consequences for the way we understand Wolf-Rayet stars and their mass loss. In addition, we derive wind masses that support the picture that the WNh stars in young stellar clusters are very massive, hydrogen-burning stars. Our findings suggest that the proximity to the Eddington limit is the physical reason for the onset of Wolf-Rayet type mass loss. This means that, e.g. in stellar evolution models, the Wolf-Rayet stage should be identified by large Eddington parameters, instead of a helium-enriched surface composition. The latter is most likely only a consequence of strong mass loss, in combination with internal mixing. For very massive stars, the enhanced Gamma-dependent mass loss is responsible for the formation of late WNh subtypes with high hydrogen surface abundances, partly close to solar. Because mass loss dominates the evolution of very massive stars, we expect a strong impact of this effect on their end products, in particular on the potential formation of black holes, and Gamma-Ray Bursts, as well as the observed upper mass limit of stars.

We present colour-magnitude diagrams (CMDs) for a sample of seven young massive clusters in the galaxies NGC 1313, NGC 1569, NGC 1705, NGC 5236 and NGC 7793. The clusters have ages in the range 5-50 million years and masses of 10^5 -10^6 Msun. Although crowding prevents us from obtaining photometry in the central regions of the clusters, we are still able to measure up to 30-100 supergiant stars in each of the richest clusters, along with the brighter main sequence stars. The resulting CMDs and luminosity functions are compared with photometry of artificially generated clusters, designed to reproduce the photometric errors and completeness as realistically as possible. In agreement with previous studies, our CMDs show no clear gap between the H-burning main sequence and the He-burning supergiant stars, contrary to predictions by common stellar isochrones. In general, the isochrones also fail to match the observed number ratios of red-to-blue supergiant stars, although the difficulty of separating blue supergiants from the main sequence complicates this comparison. In several cases we observe a large spread (1-2 mag) in the luminosities of the supergiant stars that cannot be accounted for by observational errors. This spread can be reproduced by including an age spread of 10-30 million years in the models. However, age spreads cannot fully account for the observed morphology of the CMDs and other processes, such as the evolution of interacting binary stars, may also play a role.

The recent discovery of a 2 M_sun binary millisecond pulsar (Demorest et al. 2010) has not only important consequences for the equation-of-state of nuclear matter at high densities but also raises the interesting question if the neutron star PSR J1614-2230 was born massive. The answer is vital for understanding neutron star formation in core collapse supernovae. Furthermore, this system raises interesting issues about the nature of the progenitor binary and how it evolved during its mass exchanging X-ray phase. In this paper we discuss the progenitor evolution of PSR J1614-2230. We have performed detailed stellar evolution modelling of intermediate-mass X-ray binaries undergoing Case A Roche-lobe overflow (RLO) and applied an analytic parameterization for calculating the outcome of either a common envelope evolution or the highly super-Eddington isotropic re-emission mode. We find two viable possibilities for the formation of the PSR J1614-2230 system: either it contained a 2.2-2.6 M_sun giant donor star and evolved through a common envelope and spiral-in phase or, more likely, it descended from a close binary system with a 4.0-5.0 M_sun main sequence donor star via Case A RLO. We conclude that the neutron star must have been born with a mass of ~1.95 M_sun or 1.7+-0.15 M_sun, respectively - which significantly exceeds neutron star birth masses in previously discovered radio pulsar systems. Based on the expected neutron star birth masses from considerations of stellar evolution and explosion models, we find that the progenitor star of PSR J1614-2230 is likely to have been more massive than 20 M_sun.

Some studies have claimed a universal stellar upper-mass limit of 150 Msun. A factor that is often overlooked is that there might be a difference between the current and initial masses of the most massive stars, as a result of mass loss. We present Monte Carlo mass-loss predictions for very massive stars in the range 40-300 Msun, with large luminosities and Eddington factors Gamma. Using our new dynamical approach, we find an upturn in the mass-loss vs. Gamma dependence, at the point where the winds become optically thick. This coincides with the location where wind efficiency numbers surpass the single-scattering limit of Eta = 1, reaching values up to Eta = 2.5. Our modelling suggests a transition from common O-type winds to Wolf-Rayet characteristics at the point where the winds become optically thick. This transitional behaviour is also revealed with respect to the wind acceleration parameter beta, which starts at values below 1 for the optically thin O-stars, and naturally reaches values as high as 1.5-2 for the optically thick Wolf-Rayet models. An additional finding concerns the transition in spectral morphology of the Of and WN characteristic He II line at 4686 Angstrom. When we express our mass-loss predictions as a function of the electron scattering Gamma_e (=L/M) only, we obtain a mass-loss Gamma dependence that is consistent with a previously reported power-law Mdot propto Gamma^5 (Vink 2006) that was based on our semi-empirical modelling approach. When we express Mdot in terms of both Gamma and stellar mass, we find Mdot propto M^0.8 Gamma^4.8 for our high Gamma models. Finally, we confirm that the Gamma-effect on the mass-loss predictions is much stronger than that of an increased helium abundance, calling for a fundamental revision in the way mass loss is incorporated in evolutionary models of the most massive stars.

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars - frequency, mass ratio and orbital separation - are still poorly constrained. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. We determined an orbital period of only 5.05 days - making it one of the most compact WC binaries yet identified - while we interpret a significant near-IR flare between 2004-6 as evidence for a third massive stellar component in the system in a wide eccentric orbit. The presence of a near-IR excess characteristic of hot dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short period binaries. In conjunction with previous results, the simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants and X-ray binaries. [ABRIDGED]

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems

Context. A longstanding challenge for understanding classical Cepheids is the Cepheid mass discrepancy, where theoretical mass estimates using stellar evolution and stellar pulsation calculations have been found to differ by approximately 10 - 20%. Aims. We study the role of pulsation-driven mass loss during the Cepheid stage of evolution as a possible solution to this mass discrepancy. Methods. We computed stellar evolution models with a Cepheid mass-loss prescription and various amounts of convective core overshooting. The contribution of mass loss towards the mass discrepancy is determined using these models, Results. Pulsation-driven mass loss is found to trap Cepheid evolution on the instability strip, allowing them to lose about 5 - 10% of their total mass when moderate convective core overshooting, an amount consistent with observations of other stars, is included in the stellar models. Conclusions. We find that the combination of moderate convective core overshooting and pulsation-driven mass loss can solve the Cepheid mass discrepancy.

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

When a core collapse supernova occurs in a binary system, the surviving star as well as the compact remnant emerging from the SN, may reach a substantial space velocity. With binary population synthesis modelling at solar and one fifth of solar metallicity, we predict the velocities of such runaway stars or binaries. We compile predictions for runaway OB stars, red supergiants and Wolf-Rayet stars. For those stars or binaries which undergo a second stellar explosion we compute their further evolution and the distance travelled until a Type II or Type Ibc SN or a long or short gamma-ray burst occurs. We find our predicted population of OB runaway stars broadly matches the observed population of stars but, to match the fastest observed WR runaway stars, we require that black holes receive an asymmetric kick upon formation. We find that at solar metallicity Type Ic SN progenitors travel shorter distances than the progenitors of other SN types because they are typically more massive and thus have shorter lifetimes. Those of Type IIP SN can fly farthest about 48 pc on average at solar metallicity. In considering the consequences of assuming that the progenitors of long GRBs are spun-up secondary stars that experience quasi-homogeneous evolution, we find that such evolution has a dramatic effect on the population of runaway WR stars and that some 30 per cent of GRBs could occur a hundred parsecs or more from their initial positions. We also consider mergers of double compact object binaries consisting of neutron stars and/or black holes. We find the most common type of visible mergers are neutron star--black hole mergers that are roughly ten times more common than neutron star--neutron star mergers. We also find that there may be a population of low-velocity neutron stars that are ejected from a binary rather than by their own natal kick.

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

To investigate statistically whether magnetic fields in massive stars are ubiquitous or appear in stars with specific spectral classification, certain ages, or in a special environment, we acquired 41 new spectropolarimetric observations for 36 stars. Among the observed sample roughly half of the stars are probable members of clusters at different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. Spectropolarimetric observations were obtained during three different nights using the low-resolution spectropolarimetric mode of FORS2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of the VLT. To assess the membership in open clusters and associations, we used astrometric catalogues with the best currently available kinematic and photometric data. A field at a significance level of 3sigma was detected in ten O-type stars. Importantly, the largest longitudinal magnetic fields were measured in two Of?p stars: <B_z>=-381+-122G for CPD-282561 and <B_z>=-297+-62G for HD148937, previously detected by us as magnetic. The obtained observations of HD148937 on three different nights indicate that the magnetic field is slightly variable. Our new measurements support our previous conclusion that large-scale organized magnetic fields with polar field strengths in excess of 1kG are not widespread among O-type stars. Among the stars with a detected magnetic field, only one star, HD156154, belongs to an open cluster at high membership probability. According to previous kinematic studies, four magnetic O-type stars in the sample are well-known candidate runaway stars.

We explore the role of mass loss and convective core overshoot in the evolution of Classical Cepheids. Stellar evolution models are computed with a recipe for pulsation-driven mass loss and it is found that mass loss alone is unable to account for the long-standing Cepheid mass discrepancy. However, the combination of mass loss and moderate convective core overshooting does provide a solution, bringing the amount of convective core overshooting in Cepheids closer to that found for other stars.

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

Stellar models of massive single stars are still plagued by major uncertainties. Testing and calibrating against observations is essential for their reliability. For this purpose one preferably uses observed stars that have never experienced strong binary interaction, i.e. "true single stars". However, the binary fraction among massive stars is high and identifying "true single stars" is not straight forward. Binary interaction affects systems in such a way that the initially less massive star becomes, or appears to be, single. For example, mass transfer results in a widening of the orbit and a decrease of the luminosity of the donor star, which makes it very hard to detect. After a merger or disruption of the system by the supernova explosion, no companion will be present. The only unambiguous identification of "true single stars" is possible in detached binaries, which contain two main-sequence stars. For these systems we can exclude the occurrence of mass transfer since their birth. A further advantage is that binaries can often provide us with direct measurements of the fundamental stellar parameters. Therefore, we argue these binaries are worth the effort needed to observe and analyze them. They may provide the most stringent test cases for single stellar models.

Rotation can have severe consequences for the evolution of massive stars. It is now considered as one of the main parameters, alongside mass and metallicity that determine the final fate of single stars. In massive, fast rotating stars mixing processes induced by rotation may be so efficient that helium produced in the center is mixed throughout the envelope. Such stars evolve almost chemically homogeneously. At low metallicity they remain blue and compact, while they gradually evolve into Wolf-Rayet stars and possibly into progenitors of long gamma-ray bursts. In binaries this type of evolution may occur because of (I) tides in very close binaries, as a result of (II) spin up by mass transfer, as result of (III) a merger of the two stars and (IV) when one of the components in the binary was born with a very high initial rotation rate. As these stars stay compact, the evolutionary channels are very different from what classical binary evolutionary models predict. In this contribution we discuss examples of nearly chemically homogeneous evolution in very close tidally-locked binaries. Even in such very close massive binaries, the stars may remain compact and avoid mass transfer, while Roche lobe overflow and a merger would be inevitable in the classical picture. This type of evolution may provide an alternative path to form tight Wolf-Rayet binaries and massive black hole binaries.

We investigate the evolution of Type Ib/c supernova (SN Ib/c) progenitors in close binary systems, using new evolutionary models that include the effects of rotation, with initial masses of 12 - 25 Msun for the primary components, and of single helium stars with initial masses of 2.8 - 20 Msun. We find that, despite the impact of tidal interaction on the rotation of primary stars, the amount of angular momentum retained in the core at the presupernova stage in different binary model sequences converge to a value similar to those found in previous single star models. This amount is large enough to produce millisecond pulsars, but too small to produce magnetars or long gamma-ray bursts. We employ the most up-to-date estimate for the Wolf-Rayet mass loss rate, and its implications for SN Ib/c progenitors are discussed in detail. In terms of stellar structure, SN Ib/c progenitors in binary systems are predicted to have a wide range of final masses even up to 7 Msun, with helium envelopes of 0.16 - 1.5 Msun. Our results indicate that, if the lack of helium lines in the spectra of SNe Ic were due to small amounts of helium, the distribution of both initial and final masses of SN Ic progenitors should be bimodal. Furthermore, we find that a thin hydrogen layer (0.001 - 0.01 Msun) is expected to be present in many SN Ib progenitors at the presupernova stage. We show that the presence of hydrogen, together with a rather thick helium envelope, can lead to a significant expansion of some SN Ib/c progenitors by the time of supernova explosion. This may have important consequences for the shock break-out and supernova light curve. We also argue that some SN progenitors with thin hydrogen layers produced via Case AB/B transfer might be related to Type IIb supernova progenitors with relatively small radii of about 10 Rsun.

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

Rapidly rotating stars are readily produced in binary systems. An accreting star in a binary system can be spun up by mass accretion and quickly approach the break-up limit. Mergers between two stars in a binary are expected to result in massive, fast rotating stars. These rapid rotators may appear as Be or Oe stars or at low metallicity they may be progenitors of long gamma-ray bursts. Given the high frequency of massive stars in close binaries it seems likely that a large fraction of rapidly rotating stars result from binary interaction. It is not straightforward to distinguish a a fast rotator that was born as a rapidly rotating single star from a fast rotator that resulted from some kind of binary interaction. Rapidly rotating stars resulting from binary interaction will often appear to be single because the companion tends to be a low mass, low luminosity star in a wide orbit. Alternatively, they became single stars after a merger or disruption of the binary system during the supernova explosion of the primary. The absence of evidence for a companion does not guarantee that the system did not experience binary interaction in the past. If binary interaction is one of the main causes of high stellar rotation rates, the binary fraction is expected to be smaller among fast rotators. How this prediction depend on uncertainties in the physics of the binary interactions requires further investigation.

Westerlund 1 is a young, massive Galactic starburst cluster that contains a rich coeval population of Wolf-Rayet stars, hot- and cool-phase transitional supergiants, and a magnetar. We use spectroscopic and photometric observations of the eclipsing double-lined binary W13 to derive dynamical masses for the two components, in order to determine limits for the progenitor masses of the magnetar CXOU J164710.2-455216 and the population of evolved stars in Wd1. W13 has an orbital period of 9.2709+/-0.0015 days and near-contact configuration. The shallow photometric eclipse rules out an inclination greater than 65 degrees, leading to lower limits for the masses of the emission-line optical primary and supergiant optical secondary of 21.4+/-2.6Msun and 32.8+/-4.0Msun respectively, rising to 23.2 +3.3/-3.0Msun and 35.4 +5.0/-4.6 Msun for our best-fit inclination 62 +3/-4 degrees. Comparison with theoretical models of Wolf-Rayet binary evolution suggest the emission-line object had an initial mass in excess of 35Msun, with the most likely model featuring highly non-conservative late-Case-A/Case-B mass transfer and an initial mass in excess of 40Msun. This confirms the high magnetar progenitor mass inferred from its membership in Wd1, and represents the first dynamical constraint on the progenitor mass of any magnetar. The red supergiants in Wd1 must have similar progenitor masses to W13 and are therefore amongst the most massive stars to undergo a red supergiant phase, representing a challenge for population models that suggest stars in this mass range end their redwards evolution as yellow hypergiants. [ABRIDGED]

Thermohaline mixing has recently been proposed to occur in low-mass red giants, with large consequence for the chemical yields of low-mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1Msun and 3Msun, in comparison to other mixing processes acting in these stars. We use a stellar evolution code which includes rotational mixing, internal magnetic fields and thermohaline mixing. We confirm that during the red giant stage, thermohaline mixing has the potential to decrease the abundance of ^3He which is produced earlier on the main sequence. In our models we find that this process is working on the RGB only in stars with initial mass M \simle 1.5Msun. Moreover we report that thermohaline mixing is present also during core helium burning and beyond, and has the potential to change the surface abundances of AGB stars. While we find rotational and magnetic mixing to be negligible compared to the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with the differential rotation may be essential to establish the time scale of thermohaline mixing in red giants. To explain the surface abundances observed at the bump in the luminosity function, the speed of the mixing process needs to be more than two orders of magnitude higher than in our models. However it is not clear if thermohaline mixing is the only physical process responsible for these surface abundance anomalies. Therefore, at this stage, it is not possible to calibrate the efficiency of thermohaline mixing against the observations.

The location of B supergiants in the Hertzsprung-Russell diagram (HRD) represents a long-standing problem in massive star evolution. Here we propose their nature may be revealed utilising their rotational properties, and we highlight a steep drop in massive star rotation rates at an effective temperature of 22000 K. We discuss two potential explanations for it. On the one hand, the feature might be due to the end of the main sequence, which could potentially constrain the core overshooting parameter. On the other hand, the feature might be the result of enhanced mass loss at the predicted location of the bi-stability jump. We term this effect "bi-stability breaking" and discuss its potential consequences for the evolution of massive stars.

We highlight the role of the light elements (Li, Be, B) in the evolution of massive single and binary stars, which is largely restricted to a diagnostic value, and foremost so for the element boron. However, we show that the boron surface abundance in massive early type stars contains key information about their foregoing evolution which is not obtainable otherwise. In particular, it allows to constrain internal mixing processes and potential previous mass transfer event for binary stars (even if the companion has disappeared). It may also help solving the mystery of the slowly rotating nitrogen-rich massive main sequence stars.

The formation of close massive black-hole binaries is a challenge for binary evolutionary models, especially the intriguing system M33 X-7 which harbours one of the most massive stellar-mass black holes (16 solar masses) orbiting a 70 solar mass O-star every 3.5 days. In standard binary evolution theory an episode of mass transfer or common envelope is inevitable in a binary with such a small orbital period, which complicates the formation of a black hole with such a high mass. To explain this system, we discuss a new binary evolution channel (De Mink et al. 2009), in which rotational mixing plays an important role. In very massive close binaries, tides force the rotation rate of the stars to be so high that rotationally induced mixing becomes very efficient. Helium produced in the center is mixed throughout the envelope. Instead of expanding during their mainsequence evolution (with the inevitable consequence of mass transfer), these stars stay compact, and avoid filling their Roche lobe. They gradually evolve into massive helium stars. This scenario naturally leads to the formation of very massive black holes in a very close orbit with a less evolved massive companion such as M33 X-7.

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch stars and massive stars rotating near the break-up limit. We discuss the potential of massive binaries as an interesting alternative source of processed material. We discuss observational evidence for mass shedding from interacting binaries. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material that can be ejected by massive binaries is larger than the contribution of two previously suggested sources combined.

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch (AGB) stars and massive stars rotating near the break-up limit (spin stars). We propose massive binaries as an alternative source. We compute the evolution of a 20 Msun star in a close binary and find that it sheds about 10 Msun. The ejecta are enriched in He, N, Na, and Al and depleted in C and O, similar to the abundance patterns observed in gobular cluster stars. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material ejected by massive binaries is greater than the contribution of AGB stars and spin stars combined, assuming that the majority of massive stars in a proto-globular cluster interact. If we take the possible contribution of intermediate mass stars in binaries into account and assume that the ejecta are diluted with an equal amount of unprocessed material, we find that this scenario can provide enough material to form a second generation of low-mass stars, which is as numerous as the first generation of low-mass stars. In this scenario there is no need to make commonly adopted assumptions, such as preferential loss of the first generation of stars, external pollution of the cluster, or an anomalous initial mass function. [Abridged]

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

There are large differences in the proposed progenitor models for the Cas A SNR. One of these differences is the presence or absence of a Wolf-Rayet (WR) phase of the progenitor star. The mass loss history of the progenitor star strongly affects the shape of the Supernova remnant (SNR). In this paper we investigate whether the progenitor star of Cas A had a WR phase or not and how long it may have lasted. We performed two-dimensional multi-species hydrodynamical simulations of the CSM around the progenitor star for several WR life times, each followed by the interaction of the supernova ejecta with the CSM. We then looked at the influence of the length of the WR phase and compared the results of the simulations with the observations of Cas A. The difference in the structure of the CSM, for models with different WR life times, has a strong impact on the resulting SNR. With an increasing WR life time the reverse shock velocity of the SNR decreases and the range of observed velocities in the shocked material increases. Furthermore, if a WR phase occurs, the remainders of the WR shell will be visible in the resulting SNR. Comparing our results with the observations suggests that the progenitor star of Cas A did not have a WR phase. We also find that the quasi-stationary flocculi (QSF) in Cas A are not consistent with the clumps from a WR shell that have been shocked and accelerated by the interaction with the SN ejecta. We can also conclude that for a SN explosion taking place in a CSM that is shaped by the wind during a short < 15000 yr WR phase, the clumps from the WR shell will be visible inside the SNR.

Hydrogen and/or helium accreted by a neutron star from a binary companion may undergo thermonuclear fusion. At different mass accretion rates different burning regimes are discerned. Theoretical models predict helium fusion to proceed as a thermonuclear runaway for accretion rates below the Eddington limit and as stable burning above this limit. Observations, however, place the boundary close to 10% of the Eddington limit. We study the effect of rotationally induced transport processes on the stability of helium burning. For the first time detailed calculations of thin helium shell burning on neutron stars are performed using a hydrodynamic stellar evolution code including rotation and rotationally induced magnetic fields. We find that in most cases the instabilities from the magnetic field provide the dominant contribution to the chemical mixing, while Eddington-Sweet circulations become important at high rotation rates. As helium is diffused to greater depths, the stability of the burning is increased, such that the critical accretion rate for stable helium burning is found to be lower. Combined with a higher heat flux from the crust, as suggested by recent studies, turbulent mixing could explain the observed critical accretion rate. Furthermore, close to this boundary we find oscillatory burning, which previous studies have linked to mHz QPOs. In models where we continuously lower the heat flux from the crust, the period of the oscillations increases by up to several tens of percents, similar to the observed frequency drift, suggesting that this drift could be caused by the cooling of deeper layers.

Multi wavelength observations of Cassiopeia A (Cas A) have provided us with a strong evidence for the presence of circumstellar material surrounding the progenitor star. It has been suggested that its progenitor was a massive star with a strong mass loss. But, despite the large amount of observational data from optical, IR, radio and x-ray observations, the identity of Cas A progenitor is still elusive. In this work, we compute stellar and circumstellar numerical models to look for the progenitor of Cas A. The models will be compared with the observational constraints. We have computed stellar evolution models to get time-dependent wind parameters and surface abundances. We have chosen a set of probable progenitor stars, with initial masses of 23, 28, 29, 30 and 33 \Mo, with mass loss. The derived mass loss rates and wind terminal velocities are used to simulate the evolution of the circumstellar medium. Our stellar set gives distinct SN progenitors: RSG, luminous blue super giants (LBSGs) and Wolf-Rayet (WR) stars. The 23 \Mo star explodes as a RSG in a $\rm \rho \sim r^{-2}$ dense, free streaming wind surrounded by a thin, compressed, RSG shell. The 28 $\Mo$ star explodes as a LBSG, and the SN blast wave interacts with a low density, free streaming wind surrounded by an unstable and massive ''RSG+LBSG'' shell. Finally, the 30 and 33 \Mo stars explode as WR stars surrounded by fast, WR winds that terminate in highly fragmented ''WR+RSG shell''. We have compared the surface chemical abundances of our stellar models with the observational abundances in Cas A. The abundance analysis shows that the progenitor was a star with an initial mass of the order of 30 \Mo, while the hydrodynamical analysis favors progenitors of initial masses around 23.

We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}\Lsun$, $10^{3.9}\Lsun$, and $10^{4.2}$\Lsun$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines.

Models of rotating single stars can successfully account for a wide variety of observed stellar phenomena, such as the surface enhancements of N and He. However, recent observations have questioned the idea that rotational mixing is the main process responsible for the surface enhancements, emphasizing the need for a strong and conclusive test. We investigate the consequences of rotational mixing for massive main-sequence stars in short-period binaries. In these systems the tides spin up the stars to rapid rotation. We use a state-of-the-art stellar evolution code including the effect of rotational mixing, tides, and magnetic fields. We discuss the surface abundances expected in massive close binaries (M1~20 solar masses) and we propose using such systems to test the concept of rotational mixing. As these short-period binaries often show eclipses, their parameters can be determined with high accuracy, allowing for a direct comparison with binary evolution models. In more massive close systems (M1~50 solar masses, Porb<~2 days) we find that helium is efficiently mixed throughout the envelope. The star remains blue and compact during the main-sequence phase. It stays within its Roche lobe while it gradually becomes a helium star. It is the less massive star, in which the effects of rotational mixing are less pronounced, which fills its Roche lobe first. We propose that this evolution path provides an alternative channel for the formation of tight Wolf-Rayet binaries with a main-sequence companion and might explain massive black hole binaries such as the intriguing system M33 X-7.

We have previously analysed the spectra of 135 early B-type stars in the 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 SMC metallicities with the analysis of ~50 Galactic and ~100 SMC early B-type stars with rotational velocities up to ~300km/s. The surface nitrogen abundances are utilised as a probe of the mixing process. In the SMC, we find a population of slowly rotating nitrogen-rich stars amongst the early B type core-hydrogen burning stars, similar to the LMC. In the Galactic sample we find no significant enrichment amongst the core hydrogen-burning stars, which appears to be in contrast with the expectation from both rotating single-star and close binary evolution models. However, only a small number of the rapidly rotating stars have evolved enough to produce a significant nitrogen enrichment, and these may be analogous to the non-enriched rapid rotators previously found in the LMC sample. Finally, in each metallicity regime, a population of highly enriched supergiants is observed, which cannot be the immediate descendants of core-hydrogen burning stars. Their abundances are, however, compatible with them having gone through a previous red supergiant phase. Together, these observations paint a complex picture of the nitrogen enrichment in massive main sequence and supergiant stellar atmospheres, where age and binarity cause crucial effects. Whether rotational mixing is required to understand our results remains an open question at this time, but could be answered by identifying the true binary fraction in those groups of stars that do not agree with single-star evolutionary models (abridged).

We have used VLT FLAMES data to constrain the uncertain physics of rotational mixing in stellar evolution models. We have simulated a population of single stars and find two groups of observed stars that cannot be explained: (1) a group of fast rotating stars which do not show evidence for rotational mixing and (2) a group of slow rotators with strong N enrichment. Binary effects and fossil magnetic fields may be considered to explain those two groups. We suggest that the element boron could be used to distinguish between rotational mixing and the binary scenario. Our single star population simulations quantify the expected amount of boron in fast and slow rotators and allow a comparison with measured nitrogen and boron abundances in B-stars.

One of the main uncertainties in evolutionary calculations of massive stars is the efficiency of internal mixing. It changes the chemical profile inside the star and can therefore affect the structure and further evolution. We demonstrate that eclipsing binaries, in which the tides synchronize the rotation period of the stars and the orbital period, constitute a potentially strong test for the efficiency of rotational mixing. We present detailed stellar evolutionary models of massive binaries assuming the composition of the Small Magellanic Cloud. In these models we find enhancements in the surface nitrogen abundance of up to 0.6 dex.

We present results from the first extensive study of convection zones in the envelopes of hot massive stars, which are caused by opacity peaks associated with iron and helium ionization. These convective regions can be located very close to the stellar surface. Recent observations of microturbulence in massive stars from the VLT-Flames survey are in good agreement with our predictions concerning the occurrence and the strength of sub-surface convection in hot stars. We argue further that convection close to the surface may trigger clumping at the base of the stellar wind of massive stars.

In the course of the life of a massive star, wind-wind interaction can give rise to the formation of circumstellar nebulae which are both predicted and observed in the nature. We present generic model calculations to predict the properties of such nebulae for blue supergiants. From stellar evolution calculations including rotation, we obtain the time dependence of the stellar wind properties and of the stellar radiation field. These are used as input for hydro-calculations of the circumstellar medium throughout the star's life. Here, we present the results for a rapidly rotating 12 solar masses single star. This star undergoes a blue loop during its post main sequence evolution, at the onset of which its contraction spins it up close to critical rotation. Due to the consequent anisotropic mass loss, the blue supergiant wind sweeps up the preceding slow wind into an hour glass structure. Its collision with the previously formed spherical red supergiant wind shell forms a short-lived luminous nebula consisting of two polar caps and a central inner ring. With time, the polar caps evolve into mid-latitude rings which gradually move toward the equatorial plane while the central ring is fading. These structures are reminiscent to the observed nebulae around the blue supergiant Sher 25 and the progenitor of SN 1987A. The simple model of an hour glass colliding with a spherical shell retrieves most of the intriguing nebula geometries discovered around blue supergiants, and suggests them to form an evolutionary sequence. Our results indicate that binarity is not required to obtain them.

Thermohaline mixing has recently been proposed to occur in low mass red giants, with large consequences for the chemical yields of low mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1Msun and 3Msun, in comparison to other mixing processes acting in these stars. We confirm that thermohaline mixing has the potential to destroy most of the ^3He which is produced earlier on the main sequence during the red giant stage. In our models we find that this process is working only in stars with initial mass M <~ 1.5Msun. Moreover, we report that thermohaline mixing can be present during core helium burning and beyond in stars which still have a ^3He reservoir. While rotational and magnetic mixing is negligible compared to the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with differential rotation and magnetic fields may be essential to establish the time scale of thermohaline mixing in red giants.

Rotational mixing is a very important but uncertain process in the evolution of massive stars. We propose to use close binaries to test its efficiency. Based on rotating single stellar models we predict nitrogen surface enhancements for tidally locked binaries. Furthermore we demonstrate the possibility of a new evolutionary scenario for very massive (M > 40 solar mass) close (P < 3 days) binaries: Case M, in which mixing is so efficient that the stars evolve quasi-chemically homogeneously, stay compact and avoid any Roche-lobe overflow, leading to very close (double) WR binaries.

The collapsar model requires rapidly rotating Wolf-Rayet stars as progenitors of long gamma-ray bursts. However, Galactic Wolf-Rayet stars rapidly lose angular momentum due to their intense stellar winds. We investigate whether the tidal interaction of a Wolf-Rayet star with a compact object in a binary system can spin up the Wolf-Rayet star enough to produce a collapsar. We compute the evolution of close Wolf-Rayet binaries, including tidal angular momentum exchange, differential rotation of the Wolf-Rayet star, internal magnetic fields, stellar wind mass loss, and mass transfer. The Wolf-Rayet companion is approximated as a point mass. We then employ a population synthesis code to infer the occurrence rates of the various relevant binary evolution channels. We find that the simple scenario -- i.e., the Wolf-Rayet star being tidally spun up and producing a collapsar -- does not occur at solar metallicity and may only occur with low probability at low metallicity. It is limited by the widening of the binary orbit induced by the strong Wolf-Rayet wind or by the radius evolution of the Wolf-Rayet star that most often leads to a binary merger. The tidal effects enhance the merger rate of Wolf-Rayet stars with black holes such that it becomes comparable to the occurrence rate of long gamma-ray bursts.

The VLT-FLAMES Survey of Massive Stars was an ESO Large Programme to understand rotational mixing and stellar mass-loss in different metallicity environments, in order to better constrain massive star evolution. We gathered high-quality spectra of over 800 stars in the Galaxy and in the Magellanic Clouds. A sample of this size is unprecedented, enabled by the first high-resolution, wide-field, multi-object spectrograph on an 8-m telescope. We developed spectral analysis techniques that, in combination with non-LTE, line-blanketed model atmospheres, were used to quantitatively characterise every star. The large sample, combined with the theoretical developments, has produced exciting new insights into the evolution of the most massive stars.

We review the role of rotation in massive close binary systems. Rotation has been advocated as an essential ingredient in massive single star models. However, rotation clearly is most important in massive binaries where one star accretes matter from a close companion, as the resulting spin-up drives the accretor towards critical rotation. Here, we explore our understanding of this process, and its observable consequences. When accounting for these consequences, the question remains whether rotational effects in massive single stars are still needed to explain the observations.

The prediction of the spins of the compact remnants is a fundamental goal of the theory of stellar evolution. Here, we confront the predictions for white dwarf spins from evolutionary models including rotation with observational constraints. We perform stellar evolution calculations for stars in the mass range 1... 3$\mso$, including the physics of rotation, from the zero age main sequence into the TP-AGB stage. We calculate two sets of model sequences, with and without inclusion of magnetic fields. From the final computed models of each sequence, we deduce the angular momenta and rotational velocities of the emerging white dwarfs. While models including magnetic torques predict white dwarf rotational velocities between 2 and 10 km s$^{-1}$, those from the non-magnetic sequences are found to be one to two orders of magnitude larger, well above empirical upper limits. We find the situation analogous to that in the neutron star progenitor mass range, and conclude that magnetic torques may be required in order to understand the slow rotation of compact stellar remnants in general.

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 ~300km/s 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 vsini less than 50km/s incompatible with our models ~20% of the sample). Furthermore, while we find fast rotators with enrichments in agreement with the models, the observation of evolved (log g less than 3.7dex) fast rotators that are relatively unenriched (a further ~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 are required to understand the surface properties of a population of massive main sequence stars.

We discuss how rotation and binary interactions may be related to the diversity of type Ibc supernovae and long gamma-ray bursts. After presenting recent evolutionary models of massive single and binary stars including rotation, the Tayler-Spruit dynamo and binary interactions, we argue that the nature of SNe Ibc progenitors from binary systems may not significantly differ from that of single star progenitors in terms of rotation, and that most long GRB progenitors may be produced via the quasi-chemically homogeneous evolution at sub-solar metallicity. We also briefly discuss the possible role of magnetic fields generated in the convective core of a massive star for the transport of angular momentum, which is potentially important for future stellar evolution models of supernova and GRB progenitors.

We discuss recent models on the evolution of massive stars at very low metallicity including the effects of rotation, magnetic fields and binarity. Very metal poor stars lose very little mass and angular momentum during the main sequence evolution, and rotation plays a dominant role in their evolution. In rapidly rotating massive stars, the rotationally induced mixing time scale can be even shorter than the nuclear time scale throughout the main sequence. The consequent quasi-chemically homogeneous evolution greatly differs from the standard massive star evolution that leads to formation of red giants with strong chemical stratification. Interesting outcomes of such a new mode of evolution include the formation of rapidly rotating massive Wolf-Rayet stars that emit large amounts of ionizing photons, the formation of a long gamma-ray bursts and a hypernovae, and the production of large amounts of primary nitrogen. We show that binary interactions can further enhance the effects of rotation, as mass accretion in a close binary spins up the secondary.

Rapidly rotating, chemically homogeneously evolving massive stars are considered to be progenitors of long gamma-ray bursts. We present numerical simulations of the evolution of the circumstellar medium around a rapidly rotating 20 Msol star at a metallicity of Z=0.001. Its rotation is fast enough to produce quasi-chemically homogeneous evolution. While conventionally, a star of 20 Msol would not evolve into a Wolf-Rayet stage, the considered model evolves from the main sequence directly to the helium main sequence. We use the time-dependent wind parameters, such as mass loss rate, wind velocity and rotation-induced wind anisotropy from the evolution model as input for a 2D hydrodynamical simulation. While the outer edge of the pressure-driven circumstellar bubble is spherical, the circumstellar medium close to the star shows strong non-spherical features during and after the periods of near-critical rotation. We conclude that the circumstellar medium around rapidly rotating massive stars differs considerably from the surrounding material of non-rotating stars of similar mass. Multiple blue-shifted high velocity absorption components in gamma-ray burst afterglow spectra are predicted. As a consequence of near critical rotation and short stellar evolution time scales during the last few thousand years of the star's life, we find a strong deviation of the circumstellar density profile in the polar direction from the 1/R^2 density profile normally associated with stellar winds close to the star

The collapsar model for long gamma-ray bursts requires a rapidly rotating Wolf-Rayet star as progenitor. We test the idea of producing rapidly rotating Wolf-Rayet stars in massive close binaries through mass accretion and consecutive quasi-chemically homogeneous evolution; the latter had previously been shown to provide collapsars below a certain metallicity threshold for single stars. The binary channel presented here may provide a means for massive stars to obtain the high rotation rates required to evolve quasi-chemically homogeneous and fulfill the collapsar scenario. Moreover, it suggests that a possibly large fraction of long gamma-ray bursts occurs in runaway stars.

The VLT-Flames Survey for Massive Stars (Evans05,Evans06) provides recise measurements of rotational velocities and nitrogen surface abundances of massive stars in the Magellanic Clouds. Specifically, for the first time, such abundances have been estimated for stars with significant rotational velocities. This extraordinary data set gives us the unique possibility to calibrate rotationally and magnetically induced mixing processes. Therefore, we have computed a grid of stellar evolution models varying in mass, initial rotational velocity and chemical composition. In our models we find that although magnetic fields generated by the Spruit-Taylor dynamo are essential to understand the internal angular momentum transport (and hence the rotational behavior), the corresponding chemical mixing must be neglected to reproduce the observations. Further we show that for low metallicities detailed initial abundances are of prime importance, as solar-scaled abundances may result in significant calibration errors.

Thermohaline mixing has recently been proposed to occur in low mass red giants, with large consequence for the chemical yields of low mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1 Msun and 3 Msun. We use a stellar evolution code which includes rotational mixing and internal magnetic fields. We confirm that thermohaline mixing has the potential to destroy most of the helium 3 which is produced earlier on the main sequence during the red giant stage, in stars below 1.5Msun. We find this process to continue during core helium burning and beyond. We find rotational and magnetic mixing to be negligible compared to the thermohaline mixing in the relevant layers, even if the interaction of thermohaline motions with the differential rotation may be essential to establish the time scale of thermohaline mixing in red giants.

[Abridged] We present a comprehensive study of the metallicity dependence of the mass-loss rates in stationary stellar winds of hot massive stars. Assuming a power-law dependence of mass loss on metallicity, Mdot \propto Z^{m}, and adopting a theoretical relation between the terminal velocity and metallicity, v_inf \propto Z^{0.13} (Leitherer et al.), 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 10^{5.2} Lsun. Within the errors, this result agrees with the prediction of m = 0.69 +/- 0.10 from Vink et al. For the high luminosity stars we find the mass loss rates to be greater than the predictions, implying wind clumping factors in their line-forming regions of ~4. For lower luminosity stars, the winds are so weak that their strengths cannot be reliably derived from optical lines, and one must rely on analysis of UV lines. In the low-luminosity domain the Galactic WLR is found to be much steeper than expected from theory, leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L ~ 10^{4.75} Lsun, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars are overestimated as a result of wind clumping, the degree of clumping would likely to be independent of metallicity, so the scalings derived in this study are expected to remain correct.

Pair creation supernovae (PCSN) are thought to be produced from very massive low metallicity stars. The spectacularly bright SN 2006gy does show signatures expected from PCSNe. Here, we investigate the metallicity threshold below which PCSN can form and estimate their occurrence rate. We perform stellar evolution calculations for stars of 150$\mso$ and 250$\mso$ of low metallicity (Z$_{\odot}$/5 and Z$_{\odot}$/20), and analyze their mass loss rates. We find that the bifurcation between quasi-chemically homogeneous evolution for fast rotation and conventional evolution for slower rotation, which has been found earlier for massive low metallicity stars, persists in the mass range considered here. Consequently, there are two separate PCSN progenitor types: (I) Fast rotators produce PCSNe from very massive Wolf-Rayet stars, and (II) Slower rotators that generate PCSNe in hydrogen-rich massive yellow hypergiants. We find that hydrogen-rich PCSNe could occur at metallicities as high as Z$_{\odot}$/3, which -- assuming standard IMFs are still valid to estimate their birth rates -- results in a rate of about one PCSN per 1000 supernovae in the local universe, and one PCSN per 100 supernovae at a redshift of $z=5$. PCSNe from WC-type Wolf-Rayet stars are restricted to much lower metallicity.

We study the late evolution of solar metallicity stars in the transition region between white dwarf formation and core collapse. This includes the super-asymptotic giant branch (super-AGB, SAGB) stars, which have massive enough cores to ignite carbon burning and form an oxygen-neon (ONe) core. The most massive SAGB stars have cores that may grow to the Chandrasekhar mass because of continued shell-burning. Their cores collapse, triggering a so called electron capture supernovae (ECSN). From stellar evolution models we find that the initial mass range for SAGB evolution is 7.5 ... 9.25\msun. We perform calculations with three different stellar evolution codes to investigate the sensitivity of this mass range to some of the uncertainties in current stellar models. The mass range significantly depends on the treatment of semiconvective mixing and convective overshooting. To consider the effect of a large number of thermal pulses, as expected in SAGB stars, we construct synthetic SAGB models that include a semi-analytical treatment of dredge-up, hot-bottom burning, and thermal pulse properties. This synthetic model enables us to compute the evolution of the main properties of SAGB stars from the onset of thermal pulses until the core reaches the Chandrasekhar mass or is uncovered by the stellar wind. Thereby, we determine the stellar initial mass ranges that produce ONe-white dwarfs and electron-capture supernovae. The latter is found to be 9.0 ... 9.25\msun for our fiducial model, implying that electron-capture supernovae would constitute about 4% of all supernovae in the local universe. Our synthetic approach allows us to explore the uncertainty of this number imposed by uncertainties in the third dredge-up efficiency and ABG mass loss rate. We find for ECSNe a upper limit of ~20% of all supernovae (abridged).

Van Marle et al. (2005) showed that circumstellar absorption lines in early Type Ib/c supernova and gamma-ray burst afterglow spectra may reveal the progenitor evolution of the exploding Wolf-Rayet star. While the quoted paper deals with Wolf-Rayet stars which evolved through a red supergiant stage, we investigate here the initially more massive Wolf-Rayet stars which are thought to evolve through a Luminous Blue Variable (LBV) stage. We perform hydrodynamic simulations of the evolution of the circumstellar medium around a 60 Msol star, from the main sequence through the LBV and Wolf-Rayet stages, up to core collapse. We then compute the column density of the circumstellar matter as a function of radial velocity, time and angle. This allows a comparison with the number and blue-shifts, of absorption components in the spectra of LBVs, Wolf-Rayet stars, Type Ib/c supernovae and gamma-ray burst afterglows. Our simulation for the post-LBV stage shows the formation of various absorption components, which are, however, rather short lived; they dissipate on time scales shorter than 50,000yr. As the LBV stage is thought to occur at the beginning of core helium burning, the remaining Wolf-Rayet life time is expected to be one order of magnitude larger. When interpreting the absorption components in the afterglow spectrum of GRB-021004 as circumstellar, it can be concluded that the progenitor of this source did most likely not evolve through an LBV stage. However, a close binary with late common-envelope phase (Case C) may produce a circumstellar medium that closely resembles the LBV to Wolf-Rayet evolution, but with a much shorter Wolf-Rayet period.

[Abridged] We have studied the optical spectra 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. Stellar parameters are determined using an automated fitting method, combining the stellar atmosphere code FASTWIND with the genetic-algorithm optimisation routine PIKAIA. Results for stars in the LH9 and LH10 associations of N11 are consistent with a sequential star formation scenario, in which activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2, of which the hottest is found to be ~49-54 kK (cf. ~45-46 kK for O3 stars). The masses of helium-enriched dwarfs and giants 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. For bright giants and supergiants no such mass-discrepancy is found, implying that these stars follow tracks of modestly (or non-)rotating objects. Stellar mass-loss properties were found to be intermediate to those found in massive stars in the Galaxy and the SMC, and comparisons with theoretical predictions at LMC metallicity yielded good agreement over the luminosity range of our targets, i.e. 5.0 < log L/L(sun) < 6.1.

We explore the relationship between the models for progenitor systems of Type Ia supernovae and the properties of the supernova remnants that evolve after the explosion. Most models for Type Ia progenitors in the single degenerate scenario predict substantial outflows during the presupernova evolution. Expanding on previous work, we estimate the imprint of these outflows on the structure of the circumstellar medium at the time of the supernova explosion, and the effect that this modified circumstellar medium has on the evolution of the ensuing supernova remnant. We compare our simulations with the observational properties of known Type Ia supernova remnants in the Galaxy (Kepler, Tycho, SN 1006), the Large Magellanic Cloud (0509-67.5, 0519-69.0, N103B), and M31 (SN 1885). We find that optically thick outflows from the white dwarf surface (sometimes known as accretion winds) with velocities above 200 km/s excavate large low-density cavities around the progenitors. Such large cavities are incompatible with the dynamics of the forward shock and the X-ray emission from the shocked ejecta in all the Type Ia remnants that we have examined.

We present a binary channel for the progenitors of long gamma-ray bursts. We test the idea of producing rapidly rotating Wolf-Rayet stars in massive close binaries through mass accretion and consecutive quasi-chemically homogeneous evolution. The binary channel presented here may provide a means for massive stars to obtain the required high rotation rates. Moreover, it suggests that a possibly large fraction of long gamma-ray bursts occurs in runaway stars. This can have important observational consequences for both the positions of GRBs, and their afterglow properties.

We present hydrodynamical models of circumstellar medium (CSM) of long gamma-ray burst (GRB) progenitor candidates. These are massive stars that have lost a large amount of mass in the form of stellar wind during their evolution. There are two possible ways to probe the CSM of long GRB progenitors. Firstly, the GRB afterglow consists of synchrotron radiation, emitted when the GRB jet sweeps up the surrounding medium. Therefore, the lightcurve is directly related to the density profile of the CSM. The density can either decrease with the radius squared (as is the case for a freely expanding stellar wind) or be constant (as we would expect for shocked wind or the interstellar medium). Secondly, material between the GRB and the observer will absorb part of the afterglow radiation, causing absorption lines in the afterglow spectrum. In some cases, such absorption lines are blue-shifted relative to the source indicating that the material is moving away from the progenitor star. This can be explained in terms of wind interactions in the CSM. We can use the CSM of these stars to investigate their prior evolutionary stage.

Stellar oscillations are excited in non-synchronously rotating stars in binary systems due to the tidal forces. Tangential components of the tides can drive a shear flow which behaves as a differentially forced rotating structure in a stratified outer medium. In this paper we show that our single-layer approximation for the calculation of the forced oscillations yields results that are consistent with the predictions for the synchronization timescales in circular orbits. In addition, calibrating our model results to fit Zahn's relationship between synchronization timescales and orbital separation, we are able to constrain the value of the kinematical viscosity parameter. For the 4Mo+5Mo binary modeled in this paper, the kinematical viscosity is in the range 0.0015 -- 0.0043 Ro^2/d for orbital periods in the range 2.5 -- 25 d. Furthermore, the energy dissipation rate due to the shear flow is found to decrease by ~2 orders of magnitude as synchronization is approached, implying that binary systems may approach synchronization relatively quickly but that it takes a much longer timescale to actually attain this condition. Our model can be used to calculate the energy dissipation rates throughout the orbital cycle for arbitrary values of eccentricity and stellar rotational velocity. We suggest that the asymmetric distribution of energy dissipation over the stellar surface may lead to the appearance of localized regions of enhanced surface activity.

The progenitor stars of long Gamma-Ray Bursts (GRBs) are thought to be Wolf-Rayet stars, which generate a massive and energetic wind. Nevertheless, about 25 percent of all GRB afterglows light curves indicate a constant density medium close to the exploding star. We explore various ways to produce this, by creating situations where the wind termination shock arrives very close to the star, as the shocked wind material has a nearly constant density. Typically, the distance between a Wolf-Rayet star and the wind termination shock is too large to allow afterglow formation in the shocked wind material. Here, we investigate possible causes allowing for a smaller distance: A high density or a high pressure in the surrounding interstellar medium (ISM), a weak Wolf-Rayet star wind, the presence of a binary companion, and fast motion of the Wolf-Rayet star relative to the ISM. We find that all four scenarios are possible in a limited parameter space, but that none of them is by itself likely to explain the large fraction of constant density afterglows. A low GRB progenitor metallicity, and a high GRB energy make the occurrence of a GRB afterglow in a constant density medium more likely. This may be consistent with constant densities beingpreferentially found for energetic, high redshift GRBs.

We present grids of massive star evolution models at four different metallicities (Z=0.004, 0.002, 0.001, 0.00001). The effects of rotation on the stellar structure and the transport of angular momentum and chemical elements through the Spruit-Tayler dynamo and rotationally induced instabilities are considered. After discussing uncertainties involved with the adopted physics, we elaborate the final fate of massive stars as a function of initial mass and spin rate, at each considered metallicity. In particular, we investigate for which initial conditions long gamma-ray bursts (GRBs) are expected to be produced in the frame of the collapsar model. Then, using an empirical spin distribution of young massive metal-poor stars and a specified metallicity-dependent history of star-formation, we compute the expected GRB rate as function of metallicity and redshift based on our stellar evolution models. The GRB production in our models is limited to metallicities of Z \lsim 0.004, with the consequence that about 50 % of all GRBs are predicted to be found at redshifts above z = 4, with most supernovae occurring at redshifts below z\simeq 2.2. The average GRB/SN ratio predicted by our model is about 1/200 globally, and 1/1250 at low redshift. Future strategies for testing the considered GRB progenitor scenario are briefly discussed.

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. 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 present day v_sini 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 and an effective half width of 100-150 km/s. The v_r 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, 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 v_r. The mass loss rates of the SMC objects having luminosities of log L/L_sun > 5.4 are in excellent agreement with predictions. However, for lower luminosity stars the winds are too weak to determine M_dot accurately from the optical spectrum. Two of three spectroscopically classified Vz stars from our sample are located close to the theoretical zero age main sequence, as expected.

We present 2.5D radiation-hydrodynamics simulations of the accretion-induced collapse (AIC) of white dwarfs, starting from 2D rotational equilibrium configurations of a 1.46-Msun and a 1.92-Msun model. Electron capture leads to the collapse to nuclear densities of these cores within a few tens of milliseconds. The shock generated at bounce moves slowly, but steadily, outwards. Within 50-100ms, the stalled shock breaks out of the white dwarf along the poles. The blast is followed by a neutrino-driven wind that develops within the white dwarf, in a cone of ~40deg opening angle about the poles, with a mass loss rate of 5-8 x 10^{-3} Msun/yr. The ejecta have an entropy on the order of 20-50 k_B/baryon, and an electron fraction distribution that is bimodal. By the end of the simulations, at >600ms after bounce, the explosion energy has reached 3-4 x 10^{49}erg and the total ejecta mass has reached a few times 0.001Msun. We estimate the asymptotic explosion energies to be lower than 10^{50}erg, significantly lower than those inferred for standard core collapse. The AIC of white dwarfs thus represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit weak, explosion. We document in detail the numerous effects of the fast rotation of the progenitors: The neutron stars are aspherical; the ``nu_mu'' and anti-nu_e neutrino luminosities are reduced compared to the nu_e neutrino luminosity; the deleptonized region has a butterfly shape; the neutrino flux and electron fraction depend strongly upon latitude (a la von Zeipel); and a quasi-Keplerian 0.1-0.5-Msun accretion disk is formed.

We explore the consequences of new observational and theoretical evidence that long gamma-ray bursts prefer low metallicity environments. Using recently derived mass-metallicity correlations and the mass function from SDSS studies, and adopting an average cosmic metallicity evolution from \citet{kewley2005} and \citet{savaglio2005} we derive expressions for the the relative number of massive stars formed below a given fraction of solar metallicity, $\epsilon$, as function of redshift. We demonstrate that about 1/10th of all stars form with $\epsilon < 0.1$. Therefore, a picture where the majority of GRBs form with $\epsilon < 0.1$ is not inconsistent with an empirical global SN/GRB ratio of 1/1000. It implies that (1) GRB's peak at a significantly higher redshift than supernovae; (2) massive star evolution at low metallicity may be qualitatively different and; (3) the larger the low-metallicity bias of GRBs the less likely binary evolution channels can be significant GRB producers.

Despite the growing evidence that long Gamma-Ray Bursts (GRBs) are associated with deaths of Wolf-Rayet stars, the evolutionary path of massive stars to GRBs and the exact nature of GRB progenitors remained poorly known. However, recent massive star evolutionary models indicate that -- for sufficiently low metallicity -- initially very rapidly rotating stars can satisfy the conditions for collapsar formation. Even though magnetic torques are included in these models, a strong core spin-down is avoided through quasi-chemically homogeneous evolution induced by rotational mixing. Here, we explore for which initial mass and spin-range single stars of Z=Zsun/20 are expected to produce GRBs. We further find a dichotomy in the chemical structure of GRB progenitors, where lower initial masses end their lives with a massive helium envelope which still contains some amounts of hydrogen, while higher initial masses explode with C/O-dominated hydrogen-free atmospheres.

Recent models of rotating massive stars including magnetic fields prove it difficult for the cores of single stars to retain enough angular momentum to produce a collapsar and gamma-ray burst. At low metallicity, even very massive stars may retain a massive hydrogen envelope due to the weakness of the stellar winds, posing an additional obstacle to the collapsar model. Here, we consider the evolution of massive, magnetic stars where rapid rotation induces almost chemically homogeneous evolution. We find that in this case, the requirements of the collapsar model are rather easily fulfilled if the metallicity sufficiently small: 1) Rapidly rotating helium stars are formed without the need to remove the hydrogen envelope, avoiding mass-loss induced spin-down. 2) Angular momentum transport from the helium core to hydrogen envelope by magnetic torques is insignificant. We demonstrate this by calculating evolutionary models of massive stars with various metallicities, and derive an upper metallicity limit for this scenario based on currently proposed mass loss rates. Our models also suggest the existence of a lower CO-core mass limit of about 10 Msun -- which relates to an initial mass of only about 20 Msun within our scenario -- for GRB production. We argue that the relative importance of the considered GRB progenitor channel, compared to any channel related to binary stars, may increase with decreasing metallicity, and that this channel might be the major path to GRBs from first stars.

Long gamma-ray bursts are thought to be caused by a subset of exploding Wolf-Rayet stars. We argue that the circumstellar absorption lines in early supernova and in gamma-ray burst afterglow spectra may allow us to determine the main properties of the Wolf-Rayet star progenitors which can produce those two events. To demonstrate this, we first simulate the hydrodynamic evolution of the circumstellar medium around a 40 Msun star up to the time of the supernova explosion. Knowledge of density, temperature and radial velocity of the circumstellar matter as function of space and time allows us to compute the column density in the line of sight to the centre of the nebula, as a function of radial velocity, angle, and time. Our column density profiles indicate the possible number, strengths, widths and velocities of absorption line components in supernova and gamma-ray burst afterglow spectra. Our example calculation shows four distinct line features during the Wolf-Rayet stage, at about 0, 50, 150-700 and 2200 km/s, with only those of the lowest and highest velocity present at all times. The 150-700 km/s feature decays rapidly as function of time after the onset of the Wolf-Rayet stage. It consists of a variable number of components, and, especially in its evolved stage, is depending strongly on the particular line of sight. A comparison with absorption lines detected in the afterglow of GRB 021004 suggests that the high velocity absorption component in GRB 021004 may be attributed to the free streaming Wolf-Rayet wind, which is consistent with the steep density drop indicated by the afterglow light curve. The presence of the intermediate velocity components implies that the duration of the Wolf-Rayet phase of the progenitor of GRB 021004 was much smaller than the average Wolf-Rayet life time.

This paper has been withdrawn.

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 ~ 25,000), in fields centered on the open clusters NGC 3293, NGC 4755 and NGC 6611. 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 NGC 3293 and NGC 4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared.

Since close WR+O binaries are the result of a strong interaction of both stars in massive close binary systems, they can be used to constrain the highly uncertain mass and angular momentum budget during the major mass transfer phase. We explore the progenitor evolution of the three best suited WR+O binaries HD 90657, HD 186943 and HD 211853, which are characterized by a WR/O mass ratio of $\sim$0.5 and periods of 6..10 days. We are doing so at three different levels of approximation: predicting the massive binary evolution through simple mass loss and angular momentum loss estimates, through full binary evolution models with parametrized mass transfer efficiency, and through binary evolution models including rotation of both components and a physical model which allows to compute mass and angular momentum loss from the binary system as function of time during the mass transfer process. All three methods give consistently the same answers. Our results show that, if these systems formed through stable mass transfer, their initial periods were smaller than their current ones, which implies that mass transfer has started during the core hydrogen burning phase of the initially more massive star. Furthermore, the mass transfer in all three cases must have been highly non-conservative, with on average only $\sim$10% of the transferred mass being retained by the mass receiving star. This result gives support to our system mass and angular momentum loss model, which predicts that, in the considered systems, about 90% of the overflowing matter is expelled by the rapid rotation of the mass receiver close to the $\Omega$-limit, which is reached through the accretion of the remaining 10%.

The collapsar model for gamma-ray bursts requires three essential ingredients: a massive core, removal of the hydrogen envelope, and enough angular momentum in the core. We study current massive star evolution models of solar metallicity to determine which massive star physics is capable of producing these ingredients. In particular, we investigate the role of hydrodynamic and magnetic internal angular momentum transport and binary mass and angular momentum transfer. We follow the evolution of rotating single stars and of binary systems that include rotational processes for both stars. Neglecting magnetic fields, we show that the cores of massive single stars can maintain a high specific angular momentum when evolved with the assumption that mean molecular weight gradients suppress rotational mixing processes. In binary systems that undergo mass transfer during core hydrogen burning the mass receiving star accretes large amounts of high angular momentum material, leading to a spin-up of the core. We find, however, that this merely compensates for the tidal angular momentum loss due to spin-orbit coupling, which leads to synchronous rotation before the mass transfer event. Therefore the resulting cores do not rotate faster than in single stars. We also present models that include magnetic fields generated by differential rotation and we consider the internal angular momentum transport by magnetic torques. We investigate the capability of magnetic torques to efficiently pump angular momentum into the cores of accreting stars. Despite our finding that this mechanism works, the magnetic coupling of core and envelope after the accreting star ends core hydrogen burning leads to slower rotation than in the non-magnetic case.

A recent study by Yoon & Langer (2004a) indicated that the inner cores of rapidly accreting (Mdot > 10^{-7} M_sun/yr) CO white dwarfs may rotate differentially, with a shear rate near the threshold value for the onset of the dynamical shear instability. Such differentially rotating white dwarfs obtain critical masses for thermonuclear explosion or electron-capture induced collapse which significantly exceed the canonical Chandrasekhar limit. Here, we construct two-dimensional differentially rotating white dwarf models with rotation laws resembling those of the one-dimensional models of Yoon & Langer (2004a). We derive analytic relations between the white dwarf mass, its angular momentum, and its rotational-, gravitational- and binding energy. We show that these relations are applicable for a wide range of angular velocity profiles, including solid body rotation. We demonstrate that pre-explosion and pre-collapse conditions of both, rigidly and differentially rotating white dwarfs are well established by the present work, which may facilitate future multi-dimensional simulations of Type Ia supernova explosions and studies of the formation of millisecond pulsars and gamma-ray bursts from collapsing white dwarfs.Our results lead us to suggest various possible evolutionary scenarios for progenitors of Type Ia supernovae, leading to a new paradigm of a variable mass of exploding white dwarfs, at values well above the classical Chandrasekhar mass. Based on our 2D-models, we argue for the supernova peak brightness being proportional to the white dwarf mass, which could explain various aspects of the diversity of Type Ia supernovae, such as their variation in brightness, the dependence of their mean luminosity on the host galaxy type, and the weak correlation between ejecta velocity and peak brightness.

The currently favored scenario for the progenitor evolution of Type Ia supernovae (SNe Ia) presumes that white dwarfs in close binary systems grow to the Chandrasekhar limit via mass accretion from their non-degenerate companions. However, the accreted hydrogen and/or helium usually participate thermally unstable or even violent nuclear reactions in a geometrically confined region, due to the compactness of the white dwarf. Since shell flashes induced by the thermal instability might induce significant loss of mass, efficient mass increase of white dwarfs by hydrogen and/or helium accretion has been seriously questioned. A good understanding of the stability of thermonuclear shell sources is therefore crucial in order to investigate the evolution of accreting white dwarfs as SNe Ia progenitors. Here, we present a quantitative criterion for the thermal stability of thermonuclear shell sources, and discuss the effects of rotation on the stability of helium shell burning in helium accreting CO white dwarfs with Mdot \approx 10^{-7} ... 10^{-6} Msun/yr. In particular, we show that, if the effects of rotation are properly considered, helium shell sources are significantly stabilized, which might increase the likelihood for accreting white dwarfs to grow to the Chandrasekhar limit.

We investigate the effects of rotation on the behavior of the helium burning shell source in accreting carbon-oxygen white dwarfs, in the context of the single degenerate Chandrasekhar mass progenitor scenario for Type Ia supernovae (SNe Ia). We model the evolution of helium accreting white dwarfs of initially 1 Msun, assuming four different constant accretion rates (2, 3, 5 and 10 times10^{-7} Msun/yr). In a one-dimensional approximation, we compute the mass accretion and subsequent nuclear fusion of helium into carbon and oxygen, as well as angular momentum accretion, angular momentum transport inside the white dwarf, and rotationally induced chemical mixing. Our models show two major effects of rotation: a) The helium burning nuclear shell source in the rotating models is much more stable than in corresponding non-rotating models -- which increases the likelihood of accreting white dwarfs to reach the stage of central carbon ignition. This effect is mainly due to rotationally induced mixing at the CO/He interface which widens the shell source, and due to the centrifugal force lowering the density and degeneracy at the shell source location. b) The C/O-ratio in the layers which experience helium shell burning -- which may affect the energy of a SN Ia explosion -- is strongly decreased by the rotationally induced mixing of alpha-particles into the carbon-rich layers. We discuss implications of our results for the evolution of SNe Ia progenitors.

We present a quantitative criterion for the thermal stability of thermonuclear shell sources. We find the thermal stability of shell sources to depend on exactly three factors: they are more stable when they are geometrically thicker, less degenerate and hotter. This confirms and unifies previously obtained results in terms of the geometry, temperature and density of the shell source, by a simplified but quantitative approach to the physics of shell nuclear burning. We present instability diagrams in the temperature-density plane for hydrogen and helium shell burning, which allow a simple evaluation of the stability conditions of such shell sources in stellar models. The performance of our stability criterion is demonstrated in various numerical models: in a 3 Msun AGB star, in helium accreting CO white dwarfs, in a helium white dwarf which is covered by a thin hydrogen envelope, and in a 1.0 Msun giant.

We systematically examine how the presence in a binary affects the final core structure of a massive star and its consequences for the subsequent supernova explosion. Interactions with a companion star may change the final rate of rotation, the size of the helium core, the strength of carbon burning and the final iron core mass. Stars with initial masses larger than \sim 11\Ms that experiece core collapse will generally have smaller iron cores at the time of the explosion if they lost their envelopes due to a previous binary interaction. Stars below \sim 11\Ms, on the other hand, can end up with larger helium and metal cores if they have a close companion, since the second dredge-up phase which reduces the helium core mass dramatically in single stars does not occur once the hydrogen envelope is lost. We find that the initially more massive stars in binary systems with masses in the range 8 - 11\Ms are likely to undergo an electron-capture supernova, while single stars in the same mass range would end as ONeMg white dwarfs. We suggest that the core collapse in an electron-capture supernova (and possibly in the case of relatively small iron cores) leads to a prompt explosion rather than a delayed neutrino-driven explosion and that this naturally produces neutron stars with low-velocity kicks. This leads to a dichotomous distribution of neutron star kicks, as inferred previously, where neutron stars in relatively close binaries attain low kick velocities. We illustrate the consequences of such a dichotomous kick scenario using binary population synthesis simulations and discuss its implications. This scenario has also important consequences for the minimum initial mass of a massive star that becomes a neutron star. (Abbreviated.)

Today, Type Ia supernovae are essential tools for cosmology, and recognized as major contributors to the chemical evolution of galaxies. The construction of detailed supernova progenitor models, however, was so far prevented by various physical and numerical difficulties in simulating binary systems with an accreting white dwarf component, e.g., unstable helium shell burning which may cause significant expansion and mass loss. Here, we present the first binary evolution calculation which models both stellar components and the binary interaction simultaneously, and where the white dwarf mass grows up to the Chandrasekhar limit by mass accretion. Our model starts with a 1.6 Msun helium star and a 1.0 Msun CO white dwarf in a 0.124 day orbit. Thermally unstable mass transfer starts when the CO core of the helium star reaches 0.53 Msun, with mass transfer rates of 1...8 times 10^{-6} Msun/yr. The white dwarf burns the accreted helium steadily until the white dwarf mass has reached ~ 1.3 Msun and weak thermal pulses follow until carbon ignites in the center when the white dwarf reaches 1.37 Msun. Although the supernova production rate through this channel is not well known, and this channel can not be the only one as its progenitor life time is rather short (~ 10^7 - 10^8 yr), our results indicate that helium star plus white dwarf systems form a reliable route for producing Type Ia supernovae.

We present evolutionary models of helium accreting carbon-oxygen white dwarfs in which we include the effects of the spin-up of the accreting star induced by angular momentum accretion, rotationally induced chemical mixing and rotational energy dissipation. Initial masses of 0.6 Msun and 0.8 Msun and constant accretion rates of a few times 10^{-8} Msun/yr of helium rich matter have been considered, which is typical for the sub-Chandrasekhar mass progenitor scenario for Type Ia supernovae. It is found that the helium envelope in an accreting white dwarf is heated efficiently by friction in the differentially rotating spun-up layers. As a result, helium ignites much earlier and under much less degenerate conditions compared to the corresponding non-rotating case. Consequently, a helium detonation may be avoided, which questions the sub-Chandrasekhar mass progenitor scenario for Type Ia supernovae. We discuss implications of our results for the evolution of helium star plus white dwarf binary systems as possible progenitors of recurrent helium novae.

We discuss the effects of rotation on the evolution of accreting carbon-oxygen white dwarfs, with the emphasis on possible consequences in Type Ia supernova (SN Ia) progenitors. Starting with a slowly rotating white dwarf, we simulate the accretion of matter and angular momentum from a quasi-Keplerian accretion disk. The role of the various rotationally induced hydrodynamic instabilities for the transport of angular momentum inside the white dwarf is investigated. We find that the dynamical shear instability is the most important one in the highly degenerate core. Our results imply that accreting white dwarfs rotate differentially throughout,with a shear rate close to the threshold value for the onset of the dynamical shear instability. As the latter depends on the temperature of the white dwarf, the thermal evolution of the white dwarf core is found to be relevant for the angular momentum redistribution. As found previously, significant rotation is shown to lead to carbon ignition masses well above 1.4 Msun. Our models suggest a wide range of white dwarf explosion masses, which could be responsible for some aspects of the diversity observed in SNe Ia. We analyze the potential role of the bar-mode and the r-mode instability in rapidly rotating white dwarfs, which may impose angular momentum loss by gravitational wave radiation. We discuss the consequences of the resulting spin-down for the fate of the white dwarf, and the possibility to detect the emitted gravitational waves at frequencies of 0.1 >... 1.0 Hz in nearby galaxies with LISA. Possible implications of fast and differentially rotating white dwarf cores for the flame propagation in exploding white dwarfs are also briefly discussed.

(abridged) We model the nucleosynthesis during the thermal pulse phase of a rotating, solar metallicity AGB star of 3M_sun. Rotationally induced mixing during the thermal pulses produces a layer (~2E-5M_sun) on top of the CO-core where large amounts of protons and C12 co-exist. We follow the abundance evolution in this layer, in particular that of the neutron source C13 and of the neutron poison N14. In our AGB model mixing persists during the entire interpulse phase due to the steep angular velocity gradient at the core-envelope interface. We follow the neutron production during the interpulse phase, and find a resulting maximum neutron exposure of tau_max =0.04 mbarn^-1, which is too small to produce any significant s-process. In parametric models, we then investigate the combined effects of diffusive overshooting from the convective envelope and rotationally induced mixing. Models with overshoot and weaker interpulse mixing - as perhaps expected from more slowly rotating stars - yield larger neutron exposures. We conclude that the incorporation of rotationally induce mixing processes has important consequences for the production of heavy elements in AGB stars. Through a distribution of initial rotation rates it may lead to a natural spread in the neutron exposures obtained in AGB stars of a given mass - as appears to be required by observations. Our results suggest that both processes, diffusive overshoot and rotational mixing, may be required to obtain a consistent description of the s-process in AGB stars which fulfils all observational constraints. Finally, we find that mixing due to rotation within our current framework does increase the production of N15 in the partial mixing zone, however still falling short of what seems required by observations.

We model the nucleosynthesis during a radiative interpulse phase of a rotating 3Msun Asymptotic Giant Branch (AGB) star. We find an enhanced production of the neutron source species C13 compared to non-rotating models due to shear mixing of protons and C12 at the core-envelope interface. We estimate that the resulting total production of heavy elements by slow neutron capture s-process is too low to account for most observations. This due to the fact that rotationally induced mixing during the interpulse phase causes a pollution of the C13 pocket layer with the neutron poison N14. As a result we find a maximum neutron exposure of tau_max=0.04mbarn^{-1} in the s-process layer of our solar metallicity model with rotation. This is about a factor of 5 to 10 less than required to reproduce the observed stellar s-process abundance patterns. We compare our results with models that include hydrodynamic overshooting mixing, and with simple parametric models including the combined effects of overshooting and mixing in the interpulse. Within the parametric model a range of mixing efficiencies during the interpulse phase correlates with a spread in the s-process efficiency. Such a spread is observed in AGB and post-AGB stars as well as in pre-solar SiC grains.

We discuss the first available binary evolution models which include up-to-date rotational physics for both components, as well as angular momentum accretion and spin-orbit coupling. These models allow a self-consistent computation of the mass transfer efficiency during Roche-lobe overflow phases, and a determination of the transition from quasi-conservative to non-conservative evolution. Applications to massive binary systems lead to predictions for the spin rates of compact objects in binaries, and for the occurrence of gamma-ray bursts from collapsars in binaries. Rotational effects in accreting white dwarfs are found to stabilise the shell burning and decrease the carbon abundance in progenitor models for Chandrasekhar-mass Type Ia supernovae, and to potentially avoid a detonation of the white dwarf within the sub-Chandrasekhar mass scenario.

Rotation in massive stars has been studied on the main sequence and during helium burning for decades, but only recently have realistic numerical simulations followed the transport of angular momentum that occurs during more advanced stages of evolution. The results affect such interesting issues as whether rotation is important to the explosion mechanism, whether supernovae are strong sources of gravitational radiation, the star's nucleosynthesis, and the initial rotation rate of neutron stars and black holes. We find that when only hydrodynamic instabilities (shear, Eddington-Sweet, etc.) are included in the calculation, one obtains neutron stars spinning at close to critical rotation at their surface -- or even formally in excess of critical. When recent estimates of magnetic torques (Spruit 2002) are added, however, the evolved cores spin about an order of magnitude slower. This is still more angular momentum than observed in young pulsars, but too slow for the collapsar model for gamma-ray bursts.

How massive stars die -- what sort of explosion and remnant each produces -- depends chiefly on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar winds are the only means of mass loss, and these are chiefly a function of the metallicity of the star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects the evolution and final fate of massive stars. We map, as a function of mass and metallicity, where black holes and neutron stars are likely to form and where different types of supernovae are produced. Integrating over an initial mass function, we derive the relative populations as a function of metallicity. Provided single stars rotate rapidly enough at death, we speculate upon stellar populations that might produce gamma-ray bursts and jet-driven supernovae.

We review the current basic picture of the evolution of massive stars and how their evolution and structure changes as a function of initial mass. We give an overview of the fate of modern (Pop I) and primordial (Pop III) stars with emphasis on massive and very massive stars. For single stars we show how the type of explosions, the type of remnant and their frequencies changes for different initial metallicities.

We present the results of a detailed non-LTE analysis of the UV and optical spectrum of the O6.5Iaf+ star HD153919 - the mass donor in the high-mass X-ray binary 4U1700-37. Given the eclipsing nature of the system these results allow us to determine the most likely masses of both components of the binary via Monte Carlo simulations. These suggest a mass for HD153919 of 58+/-11M_sun - implying the initial mass of the companion was rather high (>60 M_sun). The most likely mass for the compact companion is found to be 2.44+/-0.27M_sun, with only 3.5 per cent of the trials resulting in a mass less than 2.0M_sun and none less than 1.65M_sun. Our observational data is inconsistent with the canonical neutron star mass and the lowest black hole mass observed (>4.4M_sun; Nova Vel). Significantly changing observational parameters can force the compact object mass into either of these regimes but this results in the O-star mass changing by factors of greater than 2, well beyond the limits determined from its evolutionary state and surface gravity. The low mass of the compact object implies that it is difficult to form high mass black holes through both the Case A & B mass transfer channels and, if the compact object is a neutron star, would significantly constrain the high density nuclear EoS

We discuss the occurrence of the s-process during the radiative interpulse phase of rotating AGB stars. Due to differential rotation, protons are mixed into C12-rich layers after thermal pulses, in the course of the so called third dredge up episode. We follow the time evolution of key isotope abundances in the relevant layers with a post-processing code which includes time dependant mixing and nucleosynthesis. In rotating AGB models, the mixing persists during the entire interpulse phase due to the steep gradient of angular velocity at the envelope-core interface. As the layers containing protons and C12, which are formed this way, become hotter, a C13-pocket is formed in a natural way. However, in this situation also N14 is formed and spread over the entire C13-pocket. We include the neutron consuming N14(n,p) reaction in our network and determine to what extent it reduces the production of trans-iron elements. We propose that rotation may be responsible for the spread of efficiencies of the C13 neutron source as required by observations.

We present chemical abundances for four main sequence B stars in the young cluster NGC 2004 in the Large Magellanic Cloud (LMC). Apart from H II regions, unevolved OB-type stars are currently the only accessible source of present-day CNO abundances for the MCs not altered by stellar evolution. Using UVES on the VLT, we obtained spectra of sufficient resolution (R = 20 000) and signal-to-noise (S/N > 100) to derive abundances for a variety of elements (He, C, N, O, Mg and Si) with NLTE line formation. This study doubles the number of main sequence B stars in the LMC with detailed chemical abundances. More importantly and in contrast to previous studies, we find no CNO abundance anomalies brought on by e.g. binary interaction or rotational mixing. Thus, this is the first time that abundances from H II regions in the LMC can sensibly be cross-checked against those from B stars by excluding evolutionary effects. We confirm the H II-region CNO abundances to within the errors, in particular the extraordinarily low nitrogen abundance of epsilon(N) around 7.0. Taken at face value, the nebular carbon abundance is 0.16 dex below the B-star value which could be interpreted in terms of interstellar dust depletion. Oxygen abundances from the two sources agree to within 0.03 dex. In comparison with the Galactic thin disk at MC metallicities, the Magellanic Clouds are clearly nitrogen-poor environments.

We present models for the complete life and death of a 60 solar mass star evolving in a close binary system, from the main sequence phase to the formation of a compact remnant and fallback of supernova debris. After core hydrogen exhaustion, the star expands, loses most of its envelope by Roche lobe overflow, and becomes a Wolf-Rayet star. We study its post-mass transfer evolution as a function of the Wolf-Rayet wind mass loss rate (which is currently not well constrained and will probably vary with initial metallicity of the star). Varying this mass loss rate by a factor 6 leads to stellar masses at collapse that range from 3.1 to 10.7 solar masses. Although the iron core masses at collapse are generally larger for stars with larger final masses, they do not depend monotonically on the final stellar mass or even the C/O-core mass. We then compute the evolution of all models through collapse and bounce. The results range from strong supernova explosions for the lower final masses to the direct collapse of the star into a black hole for the largest final mass. Correspondingly, the final remnant masses, which were computed by following the supernova evolution and fallback of material for a time scale of about one year, are between 1.2 and 10 solar masses. We discuss the remaining uncertainties of this result and outline the consequences of our results for the understanding of the progenitor evolution of X-ray binaries and gamma-ray burst models.

Boron abundances have been derived for seven main sequence B-type stars from HST STIS spectra around the B III 2066 A line. In two stars, boron appears to be undepleted with respect to the presumed initial abundance. In one star, boron is detectable but it is clearly depleted. In the other four stars, boron is undetectable implying depletions of 1 to 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. Also, 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 analysed the halo B-type star, PG0832+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-AGB star.

The discrepancy in the past years of many more black-hole soft X-ray transients (SXTs), of which a dozen have now been identified, had challenged accepted wisdom in black hole evolution. The sheer number of inferred SXTs requires that many, if not most, stars of ZAMS masses $20-35\msun$ end up as black holes (Portegies Zwart et al. 1997; Ergma & van den Heuvel 1998). In this paper we show that this can be understood by challenging the accepted wisdom that the result of helium core burning in a massive star is independent of whether the core is covered by a hydrogen envelope, or "naked" while it burns. The latter case occurs in binaries when the envelope of the more massive star is transferred to the companion by Roche Lobe overflow while in either main sequence or red giant stage. For solar metallicity, whereas the helium cores which burn while naked essentially never go into high-mass black holes, those that burn while clothed do so, beginning at ZAMS mass $\sim 20\msun$, the precise mass depending on the $^{12}C(\alpha,\gamma)^{16}$O rate as we outline. In this way the SXTs can be evolved, provided that the H envelope of the massive star is removed only following the He core burning. We also show that in order to evolve a black hole of mass $\gsim 10\msun$ such as observed in Cyg X-1, even employing extremely massive progenitors of ZAMS mass $\gsim 60\msun$ for the black hole, the core must be covered by hydrogen during a substantial fraction of the core burning. In other words, the progenitor must be a WNL star. We evolve Cyg X-1 in an analogous way to which the SXTs are evolved, the difference being that the companion in Cyg X-1 is more massive than those in the SXTs, so that Cyg X-1 shines continuously.

We use the analysis of the heavy element abundances (C, N, O, S) in circumstellar nebulae around Luminous Blue Variables to infer the evolutionary phase in which the material has been ejected. (1) We discuss the different effects that may have changed the gas composition of the nebula since it was ejected (2) We calculate the expected abundance changes at the stellar surface due to envelope convection in the red supergiant phase. If the observed LBV nebulae are ejected during the RSG phase, the abundances of the LBV nebulae require a significantly smaller amount of mass to be lost than assumed in evolutionary models. (3) We calculate the changes in the surface composition during the main sequence phase by rotation induced mixing. If the nebulae are ejected at the end of the MS-phase, the abundances in LBV nebulae are compatible with mixing times between 5 x 10^6 and 1 x 10^7 years. The existence of ON stars supports this scenario. (4) The predicted He/H ratio in the nebulae are significantly smaller than the current observed photospheric values of their central stars. Combining various arguments we show that the LBV nebulae are ejected during the blue SG phase and that the stars have not gone through a RSG phase. The chemical enhancements are due to rotation induced mixing, and the ejection is possibly triggered by near-critical rotation. During the ejection, the outflow was optically thick, which resulted in a large effective radius and a low effective temperature. This also explains the observed properties of LBV dust.

We review recent studies on the morphology and distribution with respect to the galactic plane of planetary nebulae (PNs), as well as recent advances in MHD modeling of PNs. We discuss a tentative explanation for the connection between morphological classes and their galactic distribution

We present evolutionary calculations for 74 close binaries systems with initial primary masses in the range 12...25 M_sun, and initial secondary masses between 6 and 24 M_sun. The initial periods were chosen such that mass overflow starts during the core hydrogen burning phase of the primary (Case A), or shortly thereafter (Case B). We assume conservative evolution for contact-free systems, i.e., no mass or angular momentum loss from those system except due to stellar winds. We investigate the borderline between contact-free evolution and contact, as a function of the initial system parameters. We also investigate the effect of the treatment of convection, and found it relevant for contact and supernova order in Case A systems, particularly for the highest considered masses. For Case B systems we find contact for initial periods above approximate 10 days and below. However, in that case (and for not too large periods) contact occurs only after the mass ratio has been reversed, due to the increased fraction of the donor's convective envelope. As most In all Cases we find contact for mass ratios below approximate 0.65. We derive the observable properties of our systems after the major mass transfer event, where the mass gainer is a main sequence or supergiant O or early B type star, and the mass loser is a helium star. We point out that the assumption of conservative evolution for contact-free systems could be tested by finding helium star companions to O stars.

A thermal pulse during the post-AGB phase of stellar evolution may lead to a unique mode of light element nucleosynthesis. The stage is set by the ingestion of the unprocessed envelope material into the hot He-flash convection zone below. If the temperature is sufficiently large and the C12 abundance high enough (e.g. T_8 > 0.8, X(C12) ~ 0.4 and X(H) ~ 1E-3) protons react faster with C12 and form C13 than destroying Be7 The latter forms by alpha-capture of He3 after an initial reduction of the He3 abundance to about 3E-5 X(He4) by the ppI reaction He3(He3,2p)He4 (for T_8 ~ 1). All He3 is burned within minutes to weeks depending on the temperature. Be7 is now present at about the previously mentioned level of He3. Its further fate is determined by the reactions Be7(e-,nu)Li7 and the alpha-capture reactions of Be7 and Li7. These captures lead to the production of Be11 which in turn is finally destroyed by Be11(alpha,n)N14$. The details of this mechanism of light element production in real stars is expected to be fairly dependent on the description of mixing.

Close binaries consisting of a main sequence star and a white dwarf are considered as candidates for Type~Ia supernova progenitors. We present selfconsistent calculations of the time dependence of the structure of the main sequence star, the mass transfer rate, and the orbit by means of a binary stellar evolution program. In contrast to results based on simple estimates of the mass transfer rate in systems of the considered type, our results allow for the possibility that even systems with rather small initial white dwarf masses (~ 0.7 M_sun) may produce Type Ia supernovae. We present results for two different metallicities, Z=0.02 and Z=0.001. We find that for systems with the lower metallicity, the mass transfer rates are on average five times larger than in comparable system at solar metallicity. This leads to a systematic shift of the supernova Ia progenitor population. Firstly, while for Z=0.02 donor star initial masses in supernova progenitor systems are restricted to the range 1.6...2.3 M_sun, they are in the interval 1.4...1.8 M_sun at low Z. Secondly, the initial white dwarf masses need, on average, to be larger by 0.2 M_sun at low Z in order to obtain a Chandrasekhar mass white dwarf. This metallicity dependences may be responsible for a drop of the Type Ia supernova rate for low metallicity, and may introduce a Z-dependence in the properties of supernovae. We also estimate the X-ray luminosities of the computed systems, and investigate their donor star and orbital properties.

(abridged) Observations clearly show that low-mass AGB stars can provide a nucleosynthesis site of the s-process. Recent stellar evolution models indicate that radiative burning of C13 between thermal pulses in low-mass AGB stars may indeed provide the needed neutrons. Some mixing between the proton-rich envelope and the carbon-rich core may lead to the production of C13. However, the responsible physical mechanism is not yet unambiguously identified. We present stellar model calculations with overshoot and rotation. Overshoot, with a time-dependent and exponentially decaying efficiency, leads to a partial mixture of protons and C12 during the third dredge-up. According to the depth-dependent ratio of protons and C12, a small C13-pocket forms underneath a N14-rich layer. Overshoot does not allow for any mixing during the interpulse phase. Rotation introduces mixing driven by large angular velocity gradients which form at the envelope-core interface in AGB stars, in particular after a thermal pulse. This leads to partial mixing after a pulse, as in the case of overshoot. However, rotation continues to mix the region of the C13-pocket with a diffusion coefficient of log D ~ 2...3 [cm**2 / s]. This does not only spread the C13-pocket, but also the more massive N14-rich layer, and finally leads to mixture of both layers. By the time when the temperature there has risen to about 9*10E7K and neutron production sets in, the N14 abundance exceeds the C13 abundance by a factor of 5...10. We analyze the role of N14 as a neutron poison by considering the recycling of neutrons via N14(n,p)C14 and C12(p,gamma)N13(beta+)C13 qualitatively. We find that as long as X(N14) << X(C12), the s-process will still be possible to occur under radiative conditions.

We investigate the evolution of the surface properties of models for rotating massive stars, i.e., their luminosities, effective temperatures, surface rotational velocities, and surface abundances of all isotopes, from the zero age main sequence to the supernova stage. Our results are based on the grid of stellar models by Heger, Langer, & Woosley, which covers solar metallicity stars in the initial mass range 8 - 25 solar masses. Results are parameterized by initial mass, initial rotational velocity and major uncertainties in the treatment of the rotational mixing inside massive stars. Rotationally induced mixing processes widen the main sequence and increase the core hydrogen burning lifetime, similar to the effects of convective overshooting. It can also significantly increase the luminosity during and after core hydrogen burning, and strongly affects the evolution of the effective temperature. Our models predict surface rotational velocities for various evolutionary stages, in particular for blue supergiants, red supergiants, and for the immediate presupernova stage. We discuss the changes of the surface abundances due to rotationally induced mixing for main sequence and post main sequence stars. We single out two characteristics by which the effect of rotational mixing can be distinguished from that of massive close binary mass transfer, the only alternative process leading to non-standard chemical surface abundances in massive stars. A comparison with observed abundance anomalies in various types of massive stars supports the concept of rotational mixing in massive stars and indicates that it is responsible for most of the observed abundance anomalies.

We present new, NTT longslit spectroscopy of the B[e] supergiant in the binary system R4 in the Small Magellanic Cloud. The data show extended, forbidden N and S emissions which are typical signatures of circumstellar matter. Their extension along the space axis of the slit defines an angular size of 8.6'' which translates into a linear size of 2.4 pc. The N emission lines also show the velocity structure of a bipolar outflow expanding at 100 km/s on average. This implies that, for a measured radius of 1.2 pc, the outflow originated about 12000 yrs ago. The line flux ratio [NII]6584/[SII]6717 indicates that the nebula is nitrogen enriched and therefore it has been ejected from the central star. This is the first bipolar, ejection nebula detected around a well-established B[e] supergiant. The bipolar morphology and the chemical enrichment shown by the nebula associated with R4 are consistent with the picture of a binary merger (Langer & Heger 1998), in which R4 was originally a system composed by a close pair and a third star (the observed A companion). The close pair merged into a single star and the merging process produced a circumstellar nebula that was later shaped by the ensueing B star wind.

The progenitor evolution of the massive X-ray binary Wray 977 is investigated using new models of massive close binary evolution. These models yield constraints on the mass limit for neutron star/black hole formation in single stars, M_BH. We argue for quasi-conservative evolution in this system, and we find M_BH > 13...21 M_sun from the existence of a neutron star in Wray 977, with the uncertainty being due to uncertainties in the treatment of convection. Our results revise earlier published much larger values of M_BH derived from the parameters of Wray 977. Then, on the basis of a grid of 37 evolutionary models for massive close binaries with various initial masses, mass ratios and periods, we derive primary initial-final mass, initial mass-final helium core mass, and initial mass-final CO-core mass relations for the various mass transfer cases of close binary evolution. From these models we derive for single stars that M_BH =< 25 M_sun, independent of whether most black hole binaries formed through the Case A/B or the Case C binary channel. Using our grid of binary models, we obtain a consistent scenario for the formation of black holes in binary systems.

We present an evolutionary sequence of a low mass star from the Asymptotic Giant Branch (AGB) through its post-AGB stage, during which its surface chemical composition changes from hydrogen-rich to strongly hydrogen-deficient as consequence of a very late thermal pulse, following the so-called born-again scenario. The internal structure and abundance changes during this pulse are computed with a numerical method which allows the physically consistent calculation of stellar layers where thermonuclear and mixing time scale are comparable - a situation which occurs when the helium flash driven convection zone extends to the hydrogen-rich surface layers during the pulse peak. The final surface mass fractions are [He/C/O]=[0.38/0.36/0.22], where the high oxygen abundance is due to diffusive overshoot employed during the AGB evolution. These models are the first to achieve general agreement with the surface abundance pattern observed in hydrogen-deficient post-AGB stars - e.g. the PG 1159 stars or the WR-type central stars of planetary nebulae -, confirming the born-again scenario with a physically consistent calculation and supporting the occurrence of convective overshooting in thermally pulsing AGB stars.

The evolution of rotating stars with zero-age main sequence (ZAMS) masses in the range 8 to 25 M_sun is followed through all stages of stable evolution. The initial angular momentum is chosen such that the star's equatorial rotational velocity on the ZAMS ranges from zero to ~ 70 % of break-up. Redistribution of angular momentum and chemical species are then followed as a consequence of rotationally induced circulation and instablities. The effects of the centrifugal force on the stellar structure are included. Uncertain mixing efficiencies are gauged by observations. We find, as noted in previous work, that rotation increases the helium core masses and enriches the stellar envelopes with products of hydrogen burning. We determine, for the first time, the angular momentum distribution in typical presupernova stars along with their detailed chemical structure. Angular momentum loss due to (non-magnetic) stellar winds and the redistribution of angular momentum during core hydrogen burning are of crucial importance for the specific angular momentum of the core. Neglecting magnetic fields, we find angular momentum transport from the core to the envelope to be unimportant after core helium burning. We obtain specific angular momenta for the iron core and overlaying material of 1E16...1E17 erg s. These values are insensitive to the initial angular momentum. They are small enough to avoid triaxial deformations of the iron core before it collapses, but could lead to neutron stars which rotate close to break-up. They are also in the range required for the collapsar model of gamma-ray bursts. The apparent discrepancy with the measured rotation rates of young pulsars is discussed.

We present the first evolutionary models of intermediate mass stars up to their thermal pulses which include effects of rotation on the stellar structure as well as rotationally induced mixing of chemical species and angular momentum. We find a significant angular momentum transport from the core to the hydrogen-rich envelope and obtain a white dwarf rotation rate comparable to current observational upper limits of 50 km/s. Large angular momentum gradients at the bottom of the convective envelope and the tip of the pulse driven convective shell are shown to produce marked chemical mixing between the proton-rich and the 12C-rich layers during the so called third dredge-up. This leads to a subsequent production of 13C which is followed by neutron production through 13C(alpha,n) in radiative layers in between thermal pulses. Although uncertainties in the efficiency of rotational mixing processes persist, we conclude that rotation is capable of producing a 13C-rich layer as required for the occurrence of the s-process in TP-AGB stars.

The observed giant outbursts of Luminous Blue Variables (LBVs) may occur when these massive stars approach their Eddington limits. When this happens, they must reach a point where the centrifugal force and the radiative acceleration cancel out gravity at the equator. We call this the Omega-limit. When stars are close to the Omega-limit, strong non-spherical mass loss should occur. This suggests a scenario where a slow and very dense wind, strongly confined to the equatorial plane, is followed by a fast and almost spherical wind. We compute two-dimensional hydrodynamic models of the evolution of the nebula formed from such interacting winds, using parameters consistent with the outburst of eta Carinae in the last century. This outburst gave birth to the Homunculus, the hourglass-shaped inner part of a highly structured circumstellar nebula. Assuming the star was very close to the Omega-limit during outburst, our models produce gas distributions that strongly resemble the Homunculus on large and small scale. This supports the general conjecture that giant outbursts in LBVs occur when they approach the Eddington limit. Our models constrains the average mass loss rate since the outburst to values smaller than the present-day mass loss rate and suggest that eta Car is approaching another outburst. Our models imply that the occurrence of giant LBV outbursts depends on the initial stellar rotation rate, and that the initial angular momentum is as important to the evolution of very massive stars as their initial mass or metallicity.

Detections of extragalactic $^{15}$N are reported from observations of the rare hydrogen cyanide isotope HC$^{15}$N toward the Large Magellanic Cloud (LMC) and the core of the (post-) starburst galaxy NGC 4945. Accounting for optical depth effects, the LMC data from the massive star-forming region N113 infer a $^{14}N/^{15}$N ratio of 111 $\pm$ 17, about twice the $^{12}C/^{13}$C value. For the LMC star-forming region N159HW and for the central region of NGC 4945, $^{14}N/^{15}$N ratios are also $\approx$ 100. The $^{14}N/^{15}$N ratios are smaller than all interstellar nitrogen isotope ratios measured in the disk and center of the Milky Way, strongly supporting the idea that $^{15}$N is predominantly of `primary' nature, with massive stars being its dominant source. Although this appears to be in contradiction with standard stellar evolution and nucleosynthesis calculations, it supports recent findings of abundant $^{15}$N production due to rotationally induced mixing of protons into the helium-burning shells of massive stars.

We present the first diffraction-limited images of the mass-loss envelope of the red supergiant star VY CMa. The two-dimensional optical and NIR images were reconstructed from 3.6 m telescope speckle data using bispectrum speckle interferometry. At the wavelengths ~0.8 \mum (RG 780 filter), 1.28 \mum, and 2.17 \mum the diffraction-limited resolutions of 46 mas, 73 mas, and 124 mas were achieved. All images clearly show that the circumstellar envelope of VY CMa is non-spherical. The RG 780, 1.28 \mum, and 2.17 mum FWHM Gauss fit diameters are 67 mas * 83 mas, 80 mas * 116 mas and 138 mas * 205 mas, respectively, or 100 AU * 125 AU, 120 AU * 174 AU and 207 AU * 308 AU (for a distance of 1500 pc). We discuss several interpretations for the asymmetric morphology. Combining recent results about the angular momentum evolution of red supergiants and their pulsational properties, we suggest that VY CMa is an immediate progenitor of IRC +10420, a post red supergiant during its transformation into a Wolf-Rayet star.

We report on a mechanism which may lead to a spin-up of the surface of a rotating single star leaving the Hayashi line, which is much stronger than the spin-up expected from the mere contraction of the star. By analyzing rigidly rotating, convective stellar envelopes, we qualitatively work out the mechanism through which these envelopes may be spun up or down by mass loss through their lower or upper boundary, respectively. We find that the first case describes the situation in retreating convective envelopes, which tend to retain most of the angular momentum while becoming less massive, thereby increasing the specific angular momentum in the convection zone and thus in the layers close to the stellar surface. We explore the spin-up mechanism quantitatively in a stellar evolution calculation of a rotating 12 M_sun star, which is found to be spun up to critical rotation after leaving the red supergiant branch. We discuss implications of this spin-up for the circumstellar matter around several types of stars, i.e., post-AGB stars, B[e] stars, pre-main sequence stars, and, in particular, the progenitor of Supernova 1987A.

We summarize the present status of the predictions of massive star models for the evolution of their surface properties. After discussing luminosity, temperature and chemical composition, we focus on the question whether massive stars may arrive at critical rotation during their evolution, either on the main sequence or in later stages. We find both cases to be possible and briefly discuss observable consequences.

In this paper, we investigate the evolutionary status of B[e]~stars from the point of view of stellar evolution theory. We try to answer to the question of how massive hot supergiants --- i.e. evolved stars --- can be capable of producing a circumstellar disk. We find and discuss three possibilities: very massive evolved main sequence stars close to critical rotation due to their proximity to their Eddington-limit, blue supergiants which have just left the red supergiant branch, and single star merger remnants of a close binary system. While the latter process seems to be required to understand the properties of the spectroscopic binary R4 in the LMC, the other two scenarios may be capable of explaining the distribution of the B[e] stars in the HR~diagram. The three scenarios make different predictions about the duration of the B[e]~phase, the time integrated disk mass and the stellar properties during the B[e]~phase, which may ultimately allow to distinguish them observationally.

Ten years later, astronomers are still puzzled by the stellar evolution that produced SN 1987A --- a blue supergiant. In single star models, the new OPAL opacities make blue solutions more difficult to achieve, though still possible for certain choices of convection physics. We also consider rotation, which has the desirable effect of producing large surface enhancements of nitrogen and helium, but the undesirable effect of increasing the helium-core mass at the expense of the envelope. The latter makes blue solutions more difficult. Still, we seek a model that occurs with high probability in the LMC and for which the time-scale for making the last transition from red to blue, $\sim$ 20,000 years, has a physical interpretation --- the Kelvin-Helmholtz time of the helium core. Single star models satisfy both criteria and might yet prove to be the correct explanation for Sk -69 202, provided new rotational or convection physics can simultaneously give a blue star and explain the ring structure. Some speculations on how this might be achieved are presented and some aspects of binary models briefly discussed.

We investigate the pulsational properties of RSG models --- which we evolve from ZAMS masses in the range 10 to $20 \Msun$ --- by means of linear and non-linear calculations. We find period and growth rate of the dominant fundamental mode to increase with increasing luminosity-to-mass ratio $L/M$. Our models obtain relatively large $L/M$ values due to the inclusion of rotation in the evolutionary calculations; however, the largest values are obtained at and beyond central He-exhaustion due to major internal rearrangements of the nuclear burning regions. Our non-linear calculations as well as the behavior of the linear period and growth rate of the pulsations for periods approaching the Kelvin-Helmholtz time scale of the H-rich stellar envelope point towards the possibility of large amplitude pulsations. Such properties are similar to that found in AGB stars and suggest the possibility of a ``superwind'' to occur before the RSGs explode as supernovae. We conclude that changes in global stellar properties during the last few $10^4 \yr$ before core collapse may lead to drastic changes in the pulsational and wind properties of pre-supernova stars, with marked consequences for the immediate pre-supernova structure of the star and the circumstellar medium. We compare our results with observations of long-period OH/IR variables and discuss observational evidence for our scenario from observed supernova light curves, spectra and remnants.

The first detection of a $^{36}$S-bearing molecule in interstellar space is reported. The $J$=$2-1$ and $3-2$ transitions of C$^{36}$S have been observed toward eight Galactic molecular hot cores. From a comparison with other optically thin isotopic species of CS, the abundance ratio of $^{34}$S/$^{36}$S is 115$(\pm 17)$. This is smaller than the solar system ratio of 200 and supports the idea that $^{36}$S is, unlike the other stable sulfur isotopes, a purely secondary nucleus that is produced by s-process nucleosynthesis in massive stars.

A 12C32S, 13C32S, 12C34S, and 12C33S J = 2 - 1 line survey has been made to study interstellar 32S/34S and 34S/33S ratios from the galactic disk. The four CS isotopomers were detected in 20 star forming regions with galactocentric distances between 3 and 9 kpc. From a comparison of line velocities, the C33S J = 2 - 1 rest frequency is about 250 kHz below the value given in the Lovas (1992) catalog. Taking 12C/13C ratios from Wilson & Rood (1994) and assuming equal 12C32S and 13C32S excitation temperatures and beam filling factors, 12C32S opacities are in the range 3 to 15; average 32S/34S and 34S/33S isotope ratios are 24.4 +/- 5.0 and 6.27 +/- 1.01, respectively. While no systematic variation in the 34S/33S isotope ratio is found, the 32S/34S ratio increases with galactocentric distance when accounting for the 12C/13C gradient of the galactic disk. A fit to the unweighted data yields 32S/34S = 3.3 +/- 0.5 (dGC/kpc) + 4.1 +/- 3.1 with a correlation coefficient of 0.84. Since the interstellar sulfur (S) isotopes are synthesized by oxygen burning in massive stars, consequences for nucleosynthesis and models of chemical evolution are briefly discussed.