Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics. (arXiv:2011.06630v1 [astro-ph.HE])
<a href="http://arxiv.org/find/astro-ph/1/au:+Law_Smith_J/0/1/0/all/0/1">Jamie A. P. Law-Smith</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Everson_R/0/1/0/all/0/1">Rosa Wallace Everson</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Ramirez_Ruiz_E/0/1/0/all/0/1">Enrico Ramirez-Ruiz</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Mink_S/0/1/0/all/0/1">Selma E. de Mink</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Son_L/0/1/0/all/0/1">Lieke A. C. van Son</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Gotberg_Y/0/1/0/all/0/1">Ylva G&#xf6;tberg</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Zellmann_S/0/1/0/all/0/1">Stefan Zellmann</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Vigna_Gomez_A/0/1/0/all/0/1">Alejandro Vigna-G&#xf3;mez</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Renzo_M/0/1/0/all/0/1">Mathieu Renzo</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Wu_S/0/1/0/all/0/1">Samantha Wu</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Schroder_S/0/1/0/all/0/1">Sophie L. Schr&#xf8;der</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Foley_R/0/1/0/all/0/1">Ryan J. Foley</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Hutchinson_Smith_T/0/1/0/all/0/1">Tenley Hutchinson-Smith</a>

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_odot$
red supergiant and a 1.4$M_odot$ 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 $approx 1.1$-$2.8
R_odot$, leading to a binary neutron star that will merge within 0.01-1 Gyr.
We find an $alpha_{rm CE}$-equivalent efficiency of $approx 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 $lesssim 0.005 R_odot$ 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.

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_odot$
red supergiant and a 1.4$M_odot$ 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 $approx 1.1$-$2.8
R_odot$, leading to a binary neutron star that will merge within 0.01-1 Gyr.
We find an $alpha_{rm CE}$-equivalent efficiency of $approx 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 $lesssim 0.005 R_odot$ 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.

http://arxiv.org/icons/sfx.gif