skip to main content


Title: Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics
The coalescence of two neutron stars was recently observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12 M⊙ red supergiant and a 1.4 M⊙ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and a final orbital separation of ≈1.1 - 2.8R⊙ , leading to a binary neutron star that will merge within 0.01-1 Gyr. We find an αCE -equivalent efficiency of ≈0.1 - 0.4 for the models we study, but this may be specific for these extended progenitors. We fully resolve the core of the star to ≲0.005R⊙ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources.  more » « less
Award ID(s):
1911206
NSF-PAR ID:
10299603
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Astronomical Journal
ISSN:
2027-5943
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    In their most recent observing run, the LIGO-Virgo-KAGRA Collaboration observed gravitational waves from compact binary mergers with highly asymmetric mass ratios, including both binary black holes (BBHs) and neutron star-black holes (NSBHs). It appears that NSBHs with mass ratiosq≃ 0.2 are more common than equally asymmetric BBHs, but the reason for this remains unclear. We use the binary population synthesis codecosmicto investigate the evolutionary pathways leading to the formation and merger of asymmetric compact binaries. We find that within the context of isolated binary stellar evolution, most asymmetric mergers start off as asymmetric stellar binaries. Because of the initial asymmetry, these systems tend to first undergo a dynamically unstable mass transfer phase. However, after the first star collapses into a compact object, the mass ratio is close to unity and the second phase of mass transfer is usually stable. According to our simulations, this stable mass transfer fails to shrink the orbit enough on its own for the system to merge. Instead, the natal kick received by the second-born compact object during its collapse is key in determining how many of these systems can merge. For the most asymmetric systems with mass ratios ofq≤ 0.1, the merging systems in our models receive an average kick magnitude of 255 km s−1during the second collapse, while the average kick for non-merging systems is 59 km s−1. Because lower mass compact objects, like neutron stars, are expected to receive larger natal kicks than higher mass BHs, this may explain why asymmetric NSBH systems merge more frequently than asymmetric BBH systems.

     
    more » « less
  2. 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 in systems that tidally circularize before RLO, the donor star may be rotating subsynchronously at the onset of mass transfer. Our results demonstrate that some possible precursors to double neutron star systems are likely eccentric at the Roche radius. The effects of pre-common-envelope eccentricity on the resulting compact binary merit further study. 
    more » « less
  3. Abstract

    One of the most mysterious astrophysical states is the common envelope (CE) phase of binary evolution, in which two stars are enshrouded by the envelope shed by one of them. Interactions between the stars and the envelope shrinks the orbit. The CE can lead to mergers or to a subsequent phase of interactions. Mergers may involve any combination of two compact objects and/or stars. Some involving white dwarfs may produce Type Ia supernovae, while merging neutron stars may yield gamma-ray bursts, and merging compact objects of all kinds produce gravitational radiation. Since CEs can arise from a variety of different initial conditions, and due to the complexity of the processes involved, it is difficult to predict their end states. When many systems are being considered, as in population synthesis calculations, conservation principles are generally employed. Here we use angular momentum in a new way to derive a simple expression for the final orbital separation. This method provides advantages for the study of binaries and is particularly well suited to higher-order multiples, now considered to be important in the genesis of potential mergers. Here we focus on CEs in binaries, and the follow-up paper extends our formalism to multiple-star systems within which a CE occurs.

     
    more » « less
  4. Abstract

    The fate of stars in the zero-age main-sequence (ZAMS) range ≈8–12Mis unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSN progenitors, few have been observationally confirmed, likely due to the faintness and rapid evolution of some of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. These sources have a [Caii]λλ7291, 7324/[Oi]λλ6300, 6364 flux ratio of ≳2 in their nebular spectra. Comparing the measured [Oi] luminosity (≲1039erg s−1) and derived oxygen mass (≈0.01M) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12M. The ejecta properties (Mej≲ 1MandEkin∼ 1050erg) are also consistent. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main-sequence companion and result in the formation of a neutron star−main sequence binary. Such binaries have been suggested to be progenitors of neutron star−white dwarf systems that could merge within a Hubble time and be detectable with LISA.

     
    more » « less
  5. ABSTRACT

    Post-common envelope binaries (PCEBs) containing a white dwarf (WD) and a main-sequence (MS) star can constrain the physics of common envelope evolution and calibrate binary evolution models. Most PCEBs studied to date have short orbital periods (Porb ≲ 1 d), implying relatively inefficient harnessing of binaries’ orbital energy for envelope expulsion. Here, we present follow-up observations of five binaries from 3rd data release of Gaia mission containing solar-type MS stars and probable ultramassive WDs ($M\gtrsim 1.2\ {\rm M}_{\odot}$) with significantly wider orbits than previously known PCEBs, Porb = 18–49 d. The WD masses are much higher than expected for systems formed via stable mass transfer at these periods, and their near-circular orbits suggest partial tidal circularization when the WD progenitors were giants. These properties strongly suggest that the binaries are PCEBs. Forming PCEBs at such wide separations requires highly efficient envelope ejection, and we find that the observed periods can only be explained if a significant fraction of the energy released when the envelope recombines goes into ejecting it. Our one-dimensional stellar models including recombination energy confirm prior predictions that a wide range of PCEB orbital periods, extending up to months or years, can potentially result from Roche lobe overflow of a luminous asymptotic giant branch (AGB) star. This evolutionary scenario may also explain the formation of several wide WD + MS binaries discovered via self-lensing, as well as a significant fraction of post-AGB binaries and barium stars.

     
    more » « less