ABSTRACT Tidal dissipation due to turbulent viscosity in the convective regions of giant stars plays an important role in shaping the orbits of pre-common-envelope systems. Such systems are possible sources of transients and close compact binary systems that will eventually merge and produce detectable gravitational wave signals. Most previous studies of the onset of common envelope episodes have focused on circular orbits and synchronously rotating donor stars under the assumption that tidal dissipation can quickly spin-up the primary and circularize the orbit before the binary reaches Roche lobe overflow (RLO). We test this assumption by coupling numerical models of the post-main-sequence stellar evolution of massive stars with the model for tidal dissipation in convective envelopes developed in Vick & Lai – a tidal model that is accurate even for highly eccentric orbits with small pericentre distances. We find that, in many cases, tidal dissipation does not circularize the orbit before RLO. For a $10\, {\rm M}_{\odot }$ ($15\, {\rm M}_{\odot }$) primary star interacting with a $1.4\, {\rm M}_{\odot }$ companion, systems with pericentre distances within 3 au (6 au) when the primary leaves the main sequence will retain the initial orbital eccentricity when the primary grows to the Roche radius. Even inmore »
Dynamical tides in eccentric binaries containing massive main-sequence stars: analytical expressions
Tidal evolution of eccentric binary systems containing at least one massive main-sequence (MS) star plays an important role in the formation scenarios of merging compact-object binaries. The dominant dissipation mechanism in such systems involves tidal excitation of outgoing internal gravity waves at the convective-radiative boundary and dissipation of the waves at the stellar envelope/surface. We have derived analytical expressions for the tidal torque and tidal energy transfer rate in such binaries for arbitrary orbital eccentricities and stellar rotation rates. These expressions can be used to study the spin and orbital evolution of eccentric binaries containing massive MS stars, such as the progenitors of merging neutron star binaries. Applying our results to the PSR J0045-7319 system, which has a massive B-star companion and an observed, rapidly decaying orbit, we find that for the standard radius of convective core based on non-rotating stellar models, the B-star must have a significant retrograde and differential rotation in order to explain the observed orbital decay rate. Alternatively, we suggest that the convective core may be larger as a result of rapid stellar rotation and/or mass transfer to the B-star in the recent past during the post-MS evolution of the pulsar progenitor.
- Publication Date:
- NSF-PAR ID:
- 10361782
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 510
- Issue:
- 4
- Page Range or eLocation-ID:
- p. 4943-4951
- ISSN:
- 0035-8711
- Publisher:
- Oxford University Press
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT At least $70\, {\rm per\, cent}$ of massive OBA-type stars reside in binary or higher order systems. The dynamical evolution of these systems can lend insight into the origins of extreme phenomena such as X-ray binaries and gravitational wave sources. In one such dynamical process, the Eccentric Kozai–Lidov (EKL) mechanism, a third companion star alters the secular evolution of a binary system. For dynamical stability, these triple systems must have a hierarchical configuration. We explore the effects of a distant third companion’s gravitational perturbations on a massive binary’s orbital configuration before significant stellar evolution has taken place (≤10 Myr). We include tidal dissipation and general relativistic precession. With large (38 000 total) Monte Carlo realizations of massive hierarchical triples, we characterize imprints of the birth conditions on the final orbital distributions. Specifically, we find that the final eccentricity distribution over the range of 0.1–0.7 is an excellent indicator of its birth distribution. Furthermore, we find that the period distributions have a similar mapping for wide orbits. Finally, we demonstrate that the observed period distribution for approximately 10-Myr-old massive stars is consistent with EKL evolution.
-
ABSTRACT Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary mass ratios, eccentricities and background disc properties. We use the γ-law equation of state and adopt a robust post-processing treatment to evaluate the mass accretion rate, torque and energy transfer rate on the binary to determine its long-term orbital evolution. We find that circular comparable-mass binaries contract, with an orbital decay rate of a few times the mass doubling rate. Eccentric binaries always experience eccentricity damping. Prograde binaries with higher eccentricities or smaller mass ratios generally have slower orbital decay rates, with some extreme cases exhibiting orbital expansion. The averaged binary mass accretion rate depends on the physical size of the accretor. The accretion flows are highly variable, and the dominant variability frequency is the apparent binary orbital frequency (in the rotating frame around the central massive BH) for circular binaries but gradually shifts to the radial epicyclic frequency as the binary eccentricity increases.more »
-
Abstract The evolution of many close binary and multiple star systems is defined by phases of mass exchange and interaction. As these systems evolve into contact, tidal dissipation is not always sufficient to bring them into circular, synchronous orbits. In these cases, encounters of increasing strength occur while the orbit remains eccentric. This paper focuses on the outcomes of close tidal passages in eccentric orbits. Close eccentric passages excite dynamical oscillations about the stars’ equilibrium configurations. These tidal oscillations arise from the transfer of orbital energy into oscillation mode energy. When these oscillations reach sufficient amplitude, they break near the stellar surface. The surface wave-breaking layer forms a shock-heated atmosphere that surrounds the object. The continuing oscillations in the star’s interior launch shocks that dissipate into the atmosphere, damping the tidal oscillations. We show that the rapid, nonlinear dissipation associated with the wave breaking of fundamental oscillation modes therefore comes with coupled mass loss to the wave-breaking atmosphere. The mass ratio is an important characteristic that defines the relative importance of mass loss and energy dissipation and therefore determines the fate of systems evolving under the influence of nonlinear dissipation. The outcome can be rapid tidal circularization (
q ≪ 1) ormore » -
Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed subpopulation of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ≈10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ≈85% and ≈60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction f B ∼ 0.05 our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift z ∈ [0, 2.5], a tentative comparison with the metallicitymore »
