Development of convective envelopes in massive stars: Implications for gravitational wave sources

TL;DR

This study refines the understanding of convective envelope development in massive stars and its impact on binary evolution, significantly reducing predicted black hole merger rates and influencing gravitational wave source predictions.

Contribution

It introduces a self-consistent criterion for convective envelope presence in massive stars, improving binary evolution models and gravitational wave source predictions.

Findings

Revised CE model reduces BH-BH merger rate by a factor of 20.

Convective envelopes are highly sensitive to superadiabacity and mixing length treatments.

Predicted TZO/QS formation rate in the Galaxy is very low, making detection unlikely.

Abstract

The structure of stellar envelopes strongly influences the course and outcome of binary mass transfer, in particular of common envelope (CE) evolution. Convective envelopes can most easily be ejected during CE events, leading to short-period binaries and potentially gravitational wave (GW) sources. Conversely, radiative envelope are thought to lead to CE mergers and Thorne-Zytkow objects (TZOs) or quasi-stars (QS). Rapid binary models based on Hurley et al. (2000) often assume that any CE event with a Hertzsprung gap donor results in a CE merger, in tension with literature. We improve this with a more self-consistent criterion based on the presence of a convective envelope. Using 1D stellar models (MESA), we systematically investigate the development of convective envelopes in massive stars. We provide fitting formulae for rapid binary codes and implement them into the StarTrack population synthesis code to refine the CE treatment and examine the impact on GW sources, TZOs, and QSs. We show that convective envelopes in massive stars are highly sensitive to the treatment of superadiabacity and the mixing length. Our revised CE model significantly reduces (factor 20) the predicted merger rate of binary black hole (BH-BH) mergers with total masses between roughly 20 and 50 Msun. This leads to a bimodal mass distribution with a strong metallicity dependence. We also predict that the current TZO/QS formation rate in the Galaxy (up to roughly 10-4 yr-1), combined with their predicted lifetimes, makes their detection unlikely. Our study strongly suggests that the role of CE evolution in the formation of BH-BH mergers has been considerably overestimated for BH-BH mergers with Mtot > 20 Msun. We highlight that any prediction from the CE channel for massive BH-BH mergers (>50 Msun) heavily hinges on our limited understanding of stellar structure and mass loss close to the Eddington limit.