A synthetic population of Wolf-Rayet stars in the LMC based on detailed single and binary star evolution models

TL;DR

This study uses detailed binary and single star evolution models to create a synthetic population of Wolf-Rayet stars in the Large Magellanic Cloud, revealing the significant role of binary interactions in WR star formation.

Contribution

It provides a comprehensive binary evolution model-based synthetic WR population for the LMC, highlighting the importance of binary interactions in WR star formation and population characteristics.

Findings

Models reproduce the LMC WR population in number and luminosity.

Binary interactions are crucial for WR formation below certain luminosities.

Discrepancies suggest the need for considering observational biases and improved physics.

Abstract

Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate for the Large Magellanic Cloud (LMC). The binary models are calculated through Roche-lobe overflow and mass transfer, until the initially more massive star exhausts helium in its core. We distinguish models of WR and helium stars based on the estimated stellar wind optical depth. We use these models to build a synthetic WR population, assuming constant star formation. Our models can reproduce the WR population of the LMC to significant detail, including the number and luminosity functions of the main WR subtypes. We find that for binary fractions of 100% (50%), all LMC WR stars below $10^6\,L_{\odot}$ ($10^{5.7}\,L_{\odot}$) are stripped binary mass donors. We also identify several insightful mismatches. With a single star fraction of 50\%, our models produce too many yellow supergiants, calling either for a larger initial binary fraction, or for enhanced mass-loss near the Humphreys-Davidson limit. Our models predict more long-period WR binaries than observed, arguably due to an observational bias towards short periods. Our models also underpredict the shortest-period WR binaries, which may have implications for understanding the progenitors of double black hole mergers. The fraction of binary produced WR stars may be larger than often assumed, and outline the risk to mis-calibrate stellar physics when only single star models are used to reproduce the observed WR stars.