Post-dynamical inspiral phase of common envelope evolution

Abstract

After the companion dynamically plunges through the primary’s envelope, the two cores remain surrounded by a common envelope and the decrease of the orbital period Porb stalls. The subsequent evolution has never been systematically explored with multidimensional simulations. For this study, we performed 3D hydrodynamical simulations of an envelope evolving under the influence of a central binary star using an adaptively refined spherical grid. We followed the evolution over hundreds of orbits of the central binary to characterize the transport of angular momentum by advection, gravitational torques, turbulence, and viscosity. We find that local advective torques from the mean flow and Reynolds stresses associated with the turbulent flow dominate the angular momentum transport, which occurs outward in a disk-like structure about the orbital plane and inward along the polar axis. Turbulent transport is less efficient, but can locally significantly damp or enhance the net angular momentum radial transport and may even reverse its direction. Short-term variability in the envelope is remarkably similar to circumbinary disks, including the formation and destruction of lump-like overdensities, which enhance mass accretion and contribute to the outward transport of eccentricity generated in the vicinity of the binary. If the accretion onto the binary is allowed, the orbital decay timescale settles to a nearly constant value τb ∼ 103 to 104 Porb, while preventing accretion leads to a slowly increasing τb ∼ 105 Porb at the end of our simulations. Our results suggest that the post-dynamical orbital contraction and envelope ejection will slowly continue while the binary is surrounded by gas and that τb is often much shorter than the thermal timescale of the envelope.