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.

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.

We present the results of hydrodynamic simulations of the interaction between a 10 Jupiter mass planet and a red or asymptotic giant branch stars, both with a zero-age main sequence mass of 3.5 $M_\odot$. Dynamic in-spiral timescales are of the order of few years and a few decades for the red and asymptotic giant branch stars, respectively. The planets will eventually be destroyed at a separation from the core of the giants smaller than the resolution of our simulations, either through evaporation or tidal disruption. As the planets in-spiral, the giant stars' envelopes are somewhat puffed up. Based on relatively long timescales and even considering the fact that further in-spiral should take place before the planets are destroyed, we predict that the merger would be difficult to observe, with only a relatively small, slow brightening. Very little mass is unbound in the process. These conclusions may change if the planet's orbit enhances the star's main pulsation modes. Based on the angular momentum transfer, we also suspect that this star-planet interaction may be unable to lead to large scale outflows via the rotation-mediated dynamo effect of Nordhaus and Blackman. Detectable pollution from the destroyed planets would only result for the lightest, lowest metallicity stars. We furthermore find that in both simulations the planets move through the outer stellar envelopes at Mach-3 to Mach-5, reaching Mach-1 towards the end of the simulations. The gravitational drag force decreases and the in-spiral slows down at the sonic transition, as predicted analytically.

We assess the systematic uncertainties in stellar evolutionary calculations for low- to intermediate-mass, main-sequence stars. We compare published stellar tracks from several different evolution codes with our own tracks computed using the stellar codes STARS and MESA. In particular, we focus on tracks of 1 and 3 solar masses at solar metallicity. We find that the spread in the available 1 solar mass tracks (computed before the recent solar composition revision by Asplund et al.) can be covered by tracks between 0.97-1.01 solar masses computed with the STARS code. We assess some possible causes of the origin of this uncertainty, including how the choice of input physics and the solar constraints used to perform the solar calibration affect the tracks. We find that for a 1 solar mass track, uncertainties of around 10% in the initial hydrogen abundance and initial metallicity produce around a 2% error in mass. For the 3 solar mass tracks, there is very little difference between the tracks from the various different stellar codes. The main difference comes in the extent of the main sequence, which we believe results from the different choices of the implementation of convective overshooting in the core. Uncertainties in the initial abundances lead to a 1-2% error in the mass determination. These uncertainties cover only part of the total error budget, which should also include uncertainties in the input physics (e.g., reaction rates, opacities, convective models) and any missing physics (e.g., radiative levitation, rotation, magnetic fields). Uncertainties in stellar surface properties such as luminosity and effective temperature will further reduce the accuracy of any potential mass determinations.

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.

As part of a larger program aimed at better quantifying the uncertainties in stellar computations, we attempt to calibrate the extent of convective overshooting in low to intermediate mass stars by means of eclipsing binary systems. We model 12 such systems, with component masses between 1.3 and 6.2 solar masses, using the detailed binary stellar evolution code STARS, producing grids of models in both metallicity and overshooting parameter. From these, we determine the best fit parameters for each of our systems. For three systems, none of our models produce a satisfactory fit. For the remaining systems, no single value for the convective overshooting parameter fits all the systems, but most of our systems can be well described with an overshooting parameter between 0.09 and 0.15, corresponding to an extension of the mixed region above the core of about 0.1-0.3 pressure scale heights. Of the nine systems where we are able to obtain a good fit, seven can be reasonably well fit with a single parameter of 0.15. We find no evidence for a trend of the extent of overshooting with either mass or metallicity, though the data set is of limited size. We repeat our calculations with a second evolution code, MESA, and we find general agreement between the two codes. For the extension of the mixed region above the convective core required by the MESA models is about 0.15-0.4 pressure scale heights. For the system EI Cep, we find that MESA gives an overshooting region that is larger than the STARS one by about 0.1 pressure scale heights for the primary, while for the secondary the difference is only 0.05 pressure scale heights.

We describe and implement an adaptive particle-mesh algorithm to solve the Poisson equation for grid-based hydrodynamics codes with nested grids. The algorithm is implemented and extensively tested within the astrophysical code Enzo against the multigrid solver available by default. We find that while both algorithms show similar accuracy for smooth mass distributions, the adaptive particle-mesh algorithm is more accurate for the case of point masses, and is generally less noisy. We also demonstrate that the two-body problem can be solved accurately in a configuration with nested grids. In addition, we discuss the effect of subcycling, and demonstrate that evolving all the levels with the same timestep yields even greater precision.

In an attempt to determine how many planetary nebulae derive from binary interactions, we have started a project to measure their unbiased binary fraction. This number, when compared to the binary fraction of the presumed parent population can give a first handle on the origin of planetary nebulae. By detecting 27 bona fide central stars in the I band we have found that 30% of our sample have an I band excess between one and a few sigmas, possibly denoting companions brighter than M3-4V and with separations smaller than approximately 1000 AU. By accounting for the undetectable companions, we determine a de-biased binary fraction of 67-78% for all companions at all separations. We compare this number to a main sequence binary fraction of (50+/-4)% determined for spectral types F6V-G2V, appropriate if the progenitors of today's PN central star population is indeed the F6V-G2V stars. The error on our estimate could be between 10 and 30%. We conclude that the central star binary fraction may be larger than expected from the putative parent population. Using the more sensitive J band of a subset of 11 central stars, the binary fraction is 54% for companions brighter than approximately M5-6V and with separations smaller than about 900 AU. De-biassing this number we obtain a binary fraction of 100-107%. The two numbers should be the same and the discrepancy is likely due to small number statistics. We also present an accurately vetted compilation of observed main sequence star magnitudes, colours and masses, which can serve as a reference for future studies. We also present synthetic colours of hot stars as a function of temperature (20-170kK) and gravity (log g= 6-8) for Solar and PG1159 compositions.

We study the response of giant stars to mass loss. One-dimensional simulations of red and asymptotic giant branch stars with mass loss rates from $10^{-3}$ up to a few \msun/yr show in no case any significant radius increase. The largest radius increase of 0.2% was found in the case with the lowest mass loss rate. For dynamical-timescale mass loss rates, that may be encountered during a common envelope phase, the evolution is not adiabatic. The superadiabatic outer layer of the giant's envelope has a local thermal timescale comparable to the dynamical timescale. Therefore, this layer has enough time to readjust thermally. Moreover, the giant star is driven out of hydrostatic equilibrium and evolves dynamically. In these cases no increase of the stellar radius with respect to its initial value is found. If the mass loss rate is high enough, the superadiabaticity of the outer layer is lost progressively and a radiative zone forms due to a combination of thermal and dynamical readjustment. Conditions for unstable mass transfer based on adiabatic mass loss models that predict a significant radius increase, may need to be re-evaluated.

The recent discovery of two Earth-mass planets in close orbits around an evolved star has raised questions as to whether substellar companions can survive encounters with their host stars. We consider whether these companions could have been stripped of significant amounts of mass during the phase when they orbited through the dense inner envelopes of the giant. We apply the criterion derived by Murray et al. for disruption of gravitationally bound objects by ram pressure, to determine whether mass loss may have played a role in the histories of these and other recently discovered low-mass companions to evolved stars. We find that the brown dwarf and Jovian mass objects circling WD 0137-349, SDSS J08205+0008, and HIP 13044 are most unlikely to have lost significant mass during the common envelope phase. However, the Earth-mass planets found around KIC 05807616 could well be the remnant of one or two Jovian mass planets that lost extensive mass during the common envelope phase.

Recently, there has been a significant level of discussion of the correct treatment of Kelvin-Helmholtz instability in the astrophysical community. This discussion relies largely on how the KHI test is posed and analyzed. We pose a stringent test of the initial growth of the instability. The goal is to provide a rigorous methodology for verifying a code on two dimensional Kelvin-Helmholtz instability. We ran the problem in the Pencil Code, Athena, Enzo, NDSPHMHD, and Phurbas. A strict comparison, judgment, or ranking, between codes is beyond the scope of this work, though this work provides the mathematical framework needed for such a study. Nonetheless, how the test is posed circumvents the issues raised by tests starting from a sharp contact discontinuity yet it still shows the poor performance of Smoothed Particle Hydrodynamics. We then comment on the connection between this behavior to the underlying lack of zeroth-order consistency in Smoothed Particle Hydrodynamics interpolation. We comment on the tendency of some methods, particularly those with very low numerical diffusion, to produce secondary Kelvin-Helmholtz billows on similar tests. Though the lack of a fixed, physical diffusive scale in the Euler equations lies at the root of the issue, we suggest that in some methods an extra diffusion operator should be used to damp the growth of instabilities arising from grid noise. This statement applies particularly to moving-mesh tessellation codes, but also to fixed-grid Godunov schemes.

During the past 20 years, the idea that non-spherical planetary nebulae (PN) may need a binary or planetary interaction to be shaped was discussed by various authors. It is now generally agreed that the varied morphologies of PN cannot be fully explained solely by single star evolution. Observationally, more binary central stars of planetary nebulae (CSPN) have been discovered, opening new possibilities to understand the connections between binarity and morphology. So far, \simeq 45 binary CSPN have been detected, most being close systems detected via flux variability. To determine the PN binary fraction, one needs a method to detect wider binaries. We present here recent results obtained with the various techniques described, concentrating on binary infrared excess observations aimed at detecting binaries of any separation.

We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

The {\alpha}-formalism is a common way to parametrize the common envelope interaction between a giant star and a more compact companion. The {\alpha} parameter describes the fraction of orbital energy released by the companion that is available to eject the giant star's envelope. By using new, detailed stellar evolutionary calculations we derive a user-friendly prescription for the {\lambda} parameter and an improved approximation for the envelope binding energy, thus revising the {\alpha} equation. We then determine {\alpha} both from simulations and observations in a self consistent manner. By using our own stellar structure models as well as population considerations to reconstruct the primary's parameters at the time of the common envelope interaction, we gain a deeper understanding of the uncertainties. We find that systems with very low values of q (the ratio of the companion's mass to the mass of the primary at the time of the common envelope interaction) have higher values of {\alpha}. A fit to the data suggests that lower mass companions are left at comparable or larger orbital separations to more massive companions. We conjecture that lower mass companions take longer than a stellar dynamical time to spiral in to the giant's core, and that this is key to allowing the giant to use its own thermal energy to help unbind its envelope. As a result, although systems with light companions might not have enough orbital energy to unbind the common envelope, they might stimulate a stellar reaction that results in the common envelope ejection.