Runaway Coalescence of Pre-Common-Envelope Stellar Binaries

TL;DR

This paper combines hydrodynamic simulations and semi-analytic modeling to understand runaway mass transfer in binary star systems, successfully reproducing observed behaviors and providing a tool for analyzing merging binaries.

Contribution

It introduces a calibrated semi-analytic model of runaway binary coalescence based on hydrodynamic simulations, and releases it as an accessible Python package.

Findings

Semi-analytic model closely matches hydrodynamic results.

Model accurately reproduces orbit evolution variations with binary properties.

Application to V1309 Sco explains its orbital decay and merger timing.

Abstract

We study the process of runaway, unstable Roche lobe overflow in coalescing binary systems and its dependence on the properties of the binary involved. We create three-dimensional hydrodynamic models of binary coalescences, and follow them through a phase of increasing Roche lobe overflow until the accretor is engulfed by the donor at the onset of a common envelope phase. In these models, we vary binary properties of mass ratio, donor structure and spin, and equation of state through the gas adiabatic index. We compare the numerical results to semi-analytic models of binary orbit evolution based on mass and angular momentum exchange between two point masses. Using our hydrodynamic simulations, we measure the key parameters: the donor mass loss rate and the angular momentum exchanged per unit mass loss from the donor. Using these calibrations, the semi-analytic model closely reproduces the escalating mass loss and runaway orbital decay observed in the hydrodynamic models. The semi-analytic model accurately reproduces the major differences in orbit evolution that arise with varying mass ratio and donor structure. We encapsulate the semi-analytic model in a publicly-released python package RLOF. We apply this model to the observed period decay and subsequent merger of the binary V1309 Sco, and find that it can simultaneously reproduce the observed orbital decay and time of outburst. We further demonstrate that there is a relationship between period derivative and second derivative that can be a useful metric for evaluating candidate merging binaries.