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. More than 70% 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.

Context. With growing evidence for the existence of very massive stars at subsolar metallicity, there is an increased need for corresponding stellar evolution models. Aims. We present a dense model grid with a tailored input chemical composition appropriate for the Large Magellanic Cloud (LMC). Methods. 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 from 70 to 500 M and with initial surface rotational velocities from 0 to 550 km s −1 , covering the core-hydrogen burning phase of evolution. Results. 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 M , 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 M at zero age, 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 ZAMS surface temperature of ∼56 000 K, at ∼180 M , and to an evolution of stars in the mass range 50 M . . . 100 M to the regime of luminous blue variables in the Hertzsprung-Russell diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above ∼180 M . Conclusions. 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.

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.

We introduce a Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS spectroscopy of R136 using 17 contiguous 52 arcsec × 0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90 per cent 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 C IVλλ1548-51, including 16 O2-3 stars. For the first time, we spectroscopically classify all Weigelt and 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 He II λ1640 emission line. This emission is totally dominated by very massive stars with initial masses above ∼100 M . The presence of strong He II λ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 M . We include montages of ultraviolet spectroscopy for Large Magellanic Cloud O stars in the appendix. Future studies in this series will focus on optical STIS medium resolution observations.