We study the formation of dust in the expanding gas ejected as a result of a common envelope binary interaction. In our novel approach, we apply the dust formation model of Nozawa et al. to the outputs of the 3D hydrodynamic SPH simulation performed by Iaconi et al., that involves a giant of 0.88~\ms \ and 83~\rs, with a companion of 0.6~\ms \ placed on the surface of the giant in circular orbit. After simulating the dynamic in-spiral phase we follow the expansion of the ejecta for $\simeq 18\,000$~days. During this period the gas is able to cool down enough to reach dust formation temperatures. Our results show that dust forms efficiently in the window between $\simeq 300$~days (the end of the dynamic in-spiral) and $\simeq 5000$~days. The dust forms in two separate populations; an outer one in the material ejected during the first few orbits of the companion inside the primary's envelope and an inner one in the rest of the ejected material. The inner dust population dominates the grain size distribution at the end of the simulation. We are able to fit the grain size distribution at the end of the simulation with a double power law. The slope of the power law for smaller grains is flatter than that for larger grains, creating a knee-shaped distribution. The power law indexes are however different from the classical values determined for the interstellar medium. We also estimate that the contribution to cosmic dust by common envelope events is not negligible and comparable to that of novae and supernovae.

During the common envelope binary interaction, the expanding layers of the gaseous common envelope recombine and the resulting recombination energy has been suggested as a contributing factor to the ejection of the envelope. In this paper we perform a comparative study between simulations with and without the inclusion of recombination energy. We use two distinct setups, comprising 0.88-M$_{\odot}$ and 1.8-M$_{\odot}$ giants, that have been studied before and can serve as benchmarks. In so doing we conclude that (i) the final orbital separation is not affected by the choice of equation of state. In other words, simulations that unbind but a small fraction of the envelope result in similar final separations to those that, thanks to recombination energy, unbind a far larger fractions. (ii) The adoption of a tabulated equation of state results in a much greater fraction of unbound envelope and we demonstrate the cause of this to be the release of recombination energy. (iii) The fraction of hydrogen recombination energy that is allowed to do work should be about half of that which our adiabatic simulations use. (iv) However, for the heavier star simulation we conclude that it is helium and not hydrogen recombination energy that unbinds the gas and we determine that all helium recombination energy is thermalised in the envelope and does work. (v)~The outer regions of the expanding common envelope are likely to see the formation of dust. This dust would promote additional unbinding and shaping of the ejected envelope into axisymmetric morphologies.

We present a comparative study between the results of most hydrodynamic simulations of the common envelope binary interaction to date and observations of post common envelope binaries. The goal is to evaluate whether this dataset indicates the existence of a formula that may predict final separations of post-common envelope systems as a function of pre-common envelope parameters. Some of our conclusions are not surprising while others are more subtle. We find that: (i) Values of the final orbital separation derived from common envelope simulations must at this time be considered upper limits. Simulations that include recombination energy do not seem to have systematically different final separations; these and other simulations imply $\alpha_{\rm CE} < 0.6-1.0$. At least one simulation, {applicable to double-degenerate systems}, implies $\alpha_{\rm CE} < 0.2$. (ii) Despite large reconstruction errors, the post-RGB observations reconstructed parameters are in agreement with some of the simulations. The post-AGB observations behave instead as if they had a systematically lower value of $\alpha_{\rm CE}$. The lack of common envelope simulations with low mass AGB stars leaves us with no insight as to why this is the case. (iii) The smallest mass companion that survives the common envelope with intermediate mass giants is 0.05-0.1~\ms. (iv) Observations of binaries with separations larger than $\sim$10~\rs, tend to have high $M_2/M_1$ mass ratios and may go through a relatively long phase of unstable Roche lobe mass transfer followed by a weakened common envelope (or with no common envelope at all). (v) The effect of the spatial resolution and of the softening length on simulation results remains poorly quantified.

We investigate the common envelope binary interaction, that leads to the formation of compact binaries, such as the progenitor of Type Ia supernovae or of mergers that emit detectable gravitational waves. In this work we diverge from the classic numerical approach that models the dynamic in-spiral. We focus instead on the asymptotic behaviour of the common envelope expansion after the dynamic in-spiral terminates. We use the SPH code {\sc phantom} to simulate one of the setups from Passy et al., with a 0.88~\ms, 83~\rs \ RGB primary and a 0.6~\ms \ companion, then we follow the ejecta expansion for $\simeq 50$~yr. Additionally, we utilise a tabulated equation of state including the envelope recombination energy in the simulation (Reichardt et al.), achieving a full unbinding. We show that, as time passes, the envelope's radial velocities dominate over the tangential ones, hence allowing us to apply an homologous expansion kinematic model to the ejecta. The external layers of the envelope become homologous as soon as they are ejected, but it takes $\simeq 5000$~days ($\simeq 14$~yr) for the bulk of the unbound gas to achieve an homologous expanding regime. We observe that the complex distribution generated by the dynamic in-spiral evolves into a more ordered, ring-like shaped one in the asymptotic regime. We show that the thermodynamics of the expanding envelope are in very good agreement with those expected for an adiabatically expanding sphere under the homologous condition and give a prediction for the location and temperature of the photosphere assuming dust to be the main source of opacity.

We conduct three-dimensional hydrodynamic simulations of the common envelope binary interaction and show that if the companion were to launch jets while interacting with the giant primary star's envelope, the jets would remove a substantial fraction of the envelope's gas. We use the setup and numerical code of an earlier common envelope study that did not include jets, with a 0.88Mo, 83Ro red giant star and a 0.3Mo companion. The assumption is that the companion star accretes mass via an accretion disk that is responsible for launching the jets which, in the simulations, are injected numerically. For the first time we conduct simulations that include jets as well as the gravitational energy released by the inspiraling core-companion system. We find that simulations with jets unbind approximately three times as much envelope mass than identical simulations that do not include jets, though the total fraction of unbound gas remains below 50 per cent for these particular simulations. The jets generate high velocity outflows in the polar directions. The jets also increase the final core-companion orbital separation and lead to a kick velocity of the core-companion binary system. Our results show that, if able to form, jets could play a crucial role in ejecting the envelope and in shaping the outflow.

We have simulated a common envelope interaction of a 0.88-M$_{\odot}$, 90-R$_{\odot}$, red giant branch star and a 0.6-M$_{\odot}$, compact companion with the smoothed particle hydrodynamics code, Phantom, from the beginning of the Roche lobe overflow phase to the beginning of the self-regulated inspiral, using three different resolutions. The duration of the Roche lobe overflow phase is resolution dependent and would lengthen with increased resolution beyond the $\sim$20 years observed, while the inspiral phase and the post-common envelope separation are largely independent of resolution. Mass transfer rates through the Lagrangian points drive the orbital evolution during the Roche lobe overflow phase, as predicted analytically. The absolute mass transfer rate is resolution dependent, but always within an order of magnitude of the analytical value. Similarly, the gravitational drag in the simulations is close to the analytical approximation. This gives us confidence that simulations approximate reality. The $L_2$ and $L_3$ outflow observed during Roche lobe overflow remains bound, forming a circumbinary disk that is largely disrupted by the common envelope ejection. However, a longer phase of Roche lobe overflow and weaker common envelope ejection typical of a more stable binary may result in a surviving circumbinary disk. Finally, we examine the density distribution resulting from the interaction for simulations that include or omit the phase of Roche lobe overflow. We conclude that the degree of stability of the Roche lobe phase may modulate the shape of the subsequent planetary nebula, explaining the wide range of post-common envelope planetary nebula shapes observed.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

The common envelope binary interaction remains one of the least understood phases in the evolution of compact binaries, including those that result in Type Ia supernovae and in mergers that emit detectable gravitational waves. In this work we continue the detailed and systematic analysis of 3D hydrodynamic simulations of the common envelope interaction aimed at understanding the reliability of the results. Our first set of simulations replicate the 5 simulations of Passy et al. (a 0.88Msun, 90Rsun RGB primary with companions in the range 0.1 to 0.9Msun) using a new AMR gravity solver implemented on our modified version of the hydrodynamic code Enzo. Despite smaller final separations obtained, these more resolved simulations do not alter the nature of the conclusions that are drawn. We also carry out 5 identical simulations but with a 2.0Msun primary RGB star with the same core mass as the Passy et al. simulations, isolating the effect of the envelope binding energy. With a more bound envelope all the companions in-spiral faster and deeper though relatively less gas is unbound. Even at the highest resolution, the final separation attained by simulations with a heavier primary is similar to the size of the smoothed potential even if we account for the loss of some angular momentum by the simulation. As a result we suggest that a ~2.0$Msun RGB primary may possibly end in a merger with companions as massive as 0.6Msun, something that would not be deduced using analytical arguments based on energy conservation.

The common envelope binary interaction occurs when a star transfers mass onto a companion that cannot fully accrete it. The interaction can lead to a merger of the two objects or to a close binary. The common envelope interaction is the gateway of all evolved compact binaries, all stellar mergers and likely many of the stellar transients witnessed to date. Common envelope simulations are needed to understand this interaction and to interpret stars and binaries thought to be the byproduct of this stage. At this time, simulations are unable to reproduce the few observational data available and several ideas have been put forward to address their shortcomings. The need for more definitive simulation validation is pressing, and is already being fulfilled by observations from time-domain surveys. In this article, we present an initial method and its implementation for post-processing grid-based common envelope simulations to produce the light-curve so as to compare simulations with upcoming observations. Here we implemented a zeroth order method to calculate the light emitted from common envelope hydrodynamic simulations carried out with the 3D hydrodynamic code Enzo used in uni-grid mode. The code implements an approach for the computation of luminosity in both optically thick and optically thin regimes and is tested using the first 135 days of the common envelope simulation of Passy et al. (2012), where a 0.8 solar masses red giant branch star interacts with a 0.6 solar masses companion. This code is used to highlight two large obstacles that need to be overcome before realistic light curves can be calculated. We explain the nature of these problems and the attempted solutions and approximations in full detail to enable the next step to be identified and implemented. We also discuss our simulation in relation to recent data of transients identified as common envelope interactions.

One in 5 planetary nebulae are ejected from common envelope binary interactions but Kepler Space Telescope results are already showing this proportion to be larger. Their properties, such as abundances can be starkly different from those of the general population, so they should be considered separately when using PN as chemical or population probes. Unfortunately post-common envelope PN cannot be discerned using only their morphologies, but this will change once we couple our new common envelope simulations with PN formation models.

We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the simulation outcome. We conclude that for both codes the post-common envelope separation is somewhat larger and the amount of unbound mass slightly greater when the initial separation is wide enough that the giant does not yet overflow or just overflows its Roche lobe. PHANTOM has been adapted to the common envelope problem here for the first time and a full comparison with ENZO is presented, including an investigation of convergence as well as energy and angular momentum conservation. We also set our simulations in the context of past simulations. This comparison reveals that it is the expansion of the giant before rapid in-spiral and not spinning up of the star that causes a larger final separation. We also suggest that the large range in unbound mass for different simulations is difficult to explain and may have something to do with simulations that are not fully converged.

The Rotten Egg Nebula has at its core a binary composed of a Mira star and an A-type companion at a separation >10 au. It has been hypothesized to have formed by strong binary interactions between the Mira and a companion in an eccentric orbit during periastron passage ~800 years ago. We have performed hydrodynamic simulations of an asymptotic giant branch star interacting with companions with a range of masses in orbits with a range of initial eccentricities and periastron separations. For reasonable values of the eccentricity, we find that Roche lobe overflow can take place only if the periods are <<100 years. Moreover, mass transfer causes the system to enter a common envelope phase within several orbits. Since the central star of the Rotten Egg nebula is an AGB star, we conclude that such a common envelope phase must have lead to a merger, so the observed companion must have been a tertiary companion of a binary that merged at the time of nebula ejection. Based on the mass and timescale of the simulated disc formed around the companion before the common envelope phase, we analytically estimate the properties of jets that could be launched. Allowing for super-Eddington accretion rates, we find that jets similar to those observed are plausible, provided that the putative lost companion was relatively massive.