The impact of envelope binding energies on the merger rate density of binary compact objects

TL;DR

This study assesses how different stellar models and definitions of core boundaries influence envelope binding energies, significantly affecting the predicted merger rates of binary compact objects.

Contribution

It provides a comprehensive set of envelope binding energies from extensive stellar models and demonstrates their impact on binary merger rate predictions.

Findings

Internal energy sources can change binding energies by over an order of magnitude.

Different core boundary definitions have secondary effects on binding energies.

Using consistent stellar models alters merger rate density predictions by over an order of magnitude.

Abstract

The common envelope (CE) phase plays a key role in the formation of binary compact object systems. Its final outcome strongly depends on the envelope binding energy, but this quantity is often estimated using fitting formulas that are not fully consistent with the underlying stellar evolution models adopted in population-synthesis codes. Here, we investigate envelope binding energies across the most extensive stellar grid considered to date. Our stellar tracks, evolved with PARSEC v2.0, include hydrogen (H) -rich stars with metallicities ranging from $Z = 10^{-11}$ (Population III stars) to $Z = 0.03$, and initial masses between 2 and 2000 M$_\odot$, as well as pure-helium stars with masses from 0.36 to 350 M$_\odot$. We examine the sensitivity of the envelope binding energies to the selected core-envelope boundary definition and to different internal energy source contributions. For H-rich stars, we find that internal energy sources can alter the envelope binding energy by more than an order of magnitude, whereas the core boundary criteria play a secondary role. In contrast, for pure helium stars, the core-boundary criterion becomes the dominant factor. The envelope binding energies derived from different stellar tracks can show deviations of several orders of magnitude, with larger differences for more massive stars and higher metallicities.Finally, by implementing our new envelope binding energy prescriptions into the binary population synthesis code SEVN, we show that the predicted merger rate densities of compact binaries can differ by more than an order of magnitude compared to previous models. Our results highlight the importance of using envelope binding energies that are consistent with the underlying stellar evolution models and caution against extrapolating empirical fits beyond the considered parameter space.