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.

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.

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

Stripped-envelope stars (SESs) form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict them to be common, but only a handful of subdwarfs (i.e., SESs) with B-type companions are known. This could be the result of observational biases hindering detection, or an incorrect understanding of binary evolution. We reanalyze the well-studied post-interaction binary phi Persei. Recently, new data improved the orbital solution of the system, which contains a ~1.2 Msun SES and a rapidly rotating ~9.6 Msun Be star. We compare with an extensive grid of evolutionary models using a Bayesian approach and find initial masses of the progenitor of 7.2+/-0.4 Msun for the SES and 3.8+/-0.4 Msun for the Be star. The system must have evolved through near-conservative mass transfer. These findings are consistent with earlier studies. The age we obtain, 57+/-9 Myr, is in excellent agreement with the age of the alpha Persei cluster. We note that neither star was initially massive enough to produce a core-collapse supernova, but mass exchange pushed the Be star above the mass threshold. We find that the subdwarf is overluminous for its mass by almost an order of magnitude, compared to the expectations for a helium core burning star. We can only reconcile this if the subdwarf is in a late phase of helium shell burning, which lasts only 2-3% of the total lifetime as a subdwarf. This could imply that up to ~50 less evolved, dimmer subdwarfs exist for each system similar to phi Persei. Our findings can be interpreted as a strong indication that a substantial population of SESs indeed exists, but has so far evaded detection because of observational biases and lack of large-scale systematic searches.

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.

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}$.

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, $15^{+9}_{-8}$%, of core-collapse supernovae are `late', that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass ($4-8M_{\odot}$). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by $14^{+15}_{-14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that $\phi$ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.