Three-dimensional simulations of accretion disks in pre-CE systems

TL;DR

This study uses 3D hydrodynamic simulations to analyze the formation and properties of accretion disks in pre-common envelope binary systems, providing new insights into their structure and evolution.

Contribution

It presents the first detailed 3D simulation of RLOF-induced accretion disks in a pre-CE binary, with specific focus on a red supergiant and neutron star system.

Findings

Pre-CE disk mass estimated at ~5×10⁻³ M☉

Disk radius approximately 40 R☉

Temperature profile shallower than Shakura-Sunyaev predictions

Abstract

Before a binary system enters into a common envelope (CE) phase, accretion from the primary star onto the companion star through Roche Lobe overflow (RLOF) will lead to the formation of an accretion disk, which may generate jets. Accretion before and during the CE may alter the outcome of the interaction. Previous studies have considered different aspects of this physical mechanism. Here we study the properties of an accretion disk formed via 3D hydrodynamic simulations of the RLOF mass transfer between a 7 M$_\odot$, red supergiant star and a 1.4 M$_\odot$, neutron star companion. We simulate only the volume around the companion for improved resolution. We use a 1D implicit MESA simulation of the evolution of the system during 30,000 years between the on-set of the RLOF and the CE to guide the binary parameters and the mass-transfer rate, while we simulate only 21 years of the last part of the RLOF in 3D using an ideal gas isothermal equation of state. We expect that a pre-CE disk under these parameters will have a mass of $\sim 5\times 10^{-3}$ M$_\odot$ and a radius of $\sim$40 R$_\odot$ with a scale height of $\sim$5 R$_\odot$. The temperature profile of the disk is shallower than that predicted by the formalism of Shakura and Sunyaev, but more reasonable cooling physics would need to be included. We stress test these results with respect to a number of physical and numerical parameters, as well as simulation choices, and we expect them to be reasonable within a factor of a few for the mass and 15% for the radius. We also contextualize our results within those presented in the literature, in particular with respect to the dimensionality of simulations and the adiabatic index. We discuss the measured accretion rate in the context of the Shakura and Sunyaev formalism and debate the viscous mechanisms at play, finishing with a list of prospects for future work.