Abstract We present the successful recovery of common-envelope ejection efficiency assumed in a simulated population of double white dwarf (DWD) binaries like those which may be observed by the future Laser Interferometer Space Antenna (LISA) mission. We simulate the formation of DWD binaries by using the COSMIC population synthesis code to sample binary formation conditions such as initial mass function, metallicity of star formation, initial orbital period, and initial eccentricity. These binaries are placed in the m12i synthetic Milky Way–like galaxy, and their signal-to-noise ratio (SNR) for the LISA instrument is estimated, considering a Galactic gravitational-wave foreground informed by the population. Through the use of Fisher estimates, we construct a likelihood function for the measurement error of the LISA-bright DWD binaries (≥20 SNR,fGW≥ 5 mHz), in their gravitational-wave frequency (fGW) and chirp mass. By repeating this process for different assumptions of the common-envelope ejection efficiency, we apply Bayesian hierarchical inference to find the best match to an injected astrophysical assumption for a fiducial population model. We conclude that the impact of common-envelope ejection efficiency on the mass-transfer processes involved in DWD formation may be statistically relevant in the future observed LISA population, and that constraints on binary formation may be found by comparing simulated populations to a future observed population.
more »
« less
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
- PAR ID:
- 10299603
- Date Published:
- Journal Name:
- Astronomical Journal
- ISSN:
- 2027-5943
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Common envelope (CE) evolution is an outstanding open problem in stellar evolution, critical to the formation of compact binaries including gravitational-wave sources. In the “classical” isolated binary evolution scenario for double compact objects, the CE is usually the second mass transfer phase. Thus, the donor star of the CE is the product of a previous binary interaction, often stable Roche lobe overflow (RLOF). Because of the accretion of mass during the first RLOF, the main-sequence core of the accretor star grows and is “rejuvenated.” This modifies the core-envelope boundary region and decreases significantly the envelope binding energy for the remaining evolution. Comparing accretor stars from self-consistent binary models to stars evolved as single, we demonstrate that the rejuvenation can lower the energy required to eject a CE by ∼42%–96% for both black hole and neutron star progenitors, depending on the evolutionary stage and final orbital separation. Therefore, binaries experiencing first stable mass transfer may more easily survive subsequent CE events and result in possibly wider final separations compared to current predictions. Despite their high mass, our accretors also experience extended “blue loops,” which may have observational consequences for low-metallicity stellar populations and asteroseismology.more » « less
-
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
-
Common-envelope evolution is a stage in binary system evolution in which a giant star engulfs a companion. The standard energy formalism is an analytical framework to estimate the amount of energy transferred from the companion's shrinking orbit into the envelope of the star that engulfed it. We show analytically that this energy transfer is larger than predicted by the standard formalism. As the orbit of the companion shrinks, the mass it encloses becomes smaller, and the companion is less bound than if the enclosed mass had remained constant. Therefore, more energy must be transferred to the envelope for the orbit to shrink further. We derive a revised energy formalism that accounts for this effect, and discuss its consequences in two contexts: the formation of neutron star binaries, and the engulfment of planets and brown dwarfs by their host stars. The companion mass required to eject the stellar envelope is smaller by up to 50% , leading to differences in common-envelope evolution outcomes. The energy deposition in the outer envelope of the star, which is related to the transient luminosity and duration, is up to a factor of ≈7 higher. Common-envelope efficiency values above unity, as defined in the literature, are thus not necessarily unphysical, and result at least partly from an incomplete description of the energy deposition. The revised energy formalism presented here can improve our understanding of stellar merger and common-envelope observations and simulations.more » « less
An official website of the United States government

