Search

Post-asymptotic giant branch (post-AGB) binaries are surrounded by dusty circumbinary disks, and exhibit unexpected orbital properties resulting from poorly understood binary interaction processes. Re-accreted gas from the circumbinary disk alters the photospheric chemistry of the post-AGB star, producing a characteristic underabundance of refractory elements that correlates with condensation temperature $\unicode{x2013}$a phenomenon known as chemical depletion. This work investigates how re-accretion from a disk drives chemical depletion, and the impact accreted matter has on post-AGB evolution. We used the MESA code to evolve 0.55 and 0.60 M$_{\odot}$ post-AGB stars with the accretion of refractory element-depleted gas from a circumbinary disk. Our study adopts observationally-constrained initial accretion rates and disk masses to reproduce the chemical depletion patterns of six well-studied post-AGB binary stars: EP Lyr, HP Lyr, IRAS 17038-4815, IRAS 09144-4933, HD 131356, and SX Cen. We find high accretion rates ($>\,$10$^{-7}$ M$_{\odot}$yr$^{-1}$) and large disk masses ($\geq\,$10$^{-2}$ M$_{\odot}$) necessary to reproduce observed depletion, particularly in higher-mass, hotter post-AGB stars (T$_{\textrm{eff}}\geq$ 6000 K). A slower evolution (lower core mass) is required to reproduce cooler (T$_{\textrm{eff}}\leq$ 5000 K) depleted post-AGB stars. Rapid accretion significantly impacts post-AGB evolution, stalling stars at cooler effective temperatures and extending post-AGB lifetimes by factors of around 3 to 10. Despite this, extended post-AGB timescales remain within or below the planetary nebula (PN) visibility timescale, suggesting accretion cannot account for the observed lack of ionised PNe in post-AGB binaries. Our findings constrain accretion-flow parameters and advance our understanding of disk-binary interactions in post-AGB systems., Comment: 14 pages, 8 figures, 5 tables

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 consider what properties of magnetic fields and jets may be supported by our disk and discuss prospects for future work., Comment: 17 pages, 20 figures

We present the first interferometric imaging of molecular line emission from the Ring Nebula, NGC~6720, in the form of Submillimeter Array (SMA) observations of CO $J=2\rightarrow 1$ emission. The SMA $^{12}$CO(2--1) mapping data, with $\sim$3$''$ spatial resolution and 2 km s$^{-1}$ velocity resolution, provide an unprecedentedly detailed, 3D view of the Ring's clumpy molecular envelope. The morphology of the velocity-integrated SMA $^{12}$CO(2--1) image closely resembles those of near-IR H$_2$ and PAH emission in JWST/NIRCam imaging of NGC~6720, with the molecular gas forming a geometrically thin layer surrounding the ionized gas imaged by HST and JWST. A simple, geometrical model of the $^{12}$CO(2--1) data shows that the intrinsic structure of NGC~6720's molecular envelope closely resembles a truncated, triaxial ellipsoid that is viewed close to pole-on, and that the dynamical age of the molecular envelope is $\sim$6000 yr. The SMA $^{12}$CO(2--1) data furthermore reveal that filamentary features seen projected in the Ring's interior in JWST imaging are in fact fast-moving polar knots or bullets with radial velocities of $\pm$45--50 km s$^{-1}$ relative to systemic, and that the hot progenitor star remnant is positioned at the precise geometric center of the clumpy, ellipsoidal molecular shell. We assert that the Ring's molecular envelope was formed via a relatively sudden, AGB-terminating mass ejection event $\sim$6000 yr ago, and that this ellipsoidal envelope was then punctured by fast, collimated polar outflows resulting from interactions between the progenitor and one or more companion stars. Such an evolutionary scenario may describe most molecule-rich, ``Ring-like'' planetary nebulae., Comment: 17 pages, 10 figures; to appear in The Astrophysical Journal (accepted 20 Jan. 2025)

Planetary nebulae are formed by the matter ejected by low-to-intermediate mass stars (~0.8-8 times the mass of the Sun) towards the end of their lives. As hydrogen and then helium fuel sources run out, stars expand. During these giant phases stars also lose sizable amounts of mass. During the second giant phase, after the exhaustion of core helium, the mass loss is so great that stars lose a large fraction of their mass (50 - 90%), leaving behind a small, hot core, known as a white dwarf, surrounded by a nebula. Planetary nebulae are the result of many processes that shape and alter their ionization structure and chemical composition. The resulting nebula, illuminated by the ultraviolet-rich spectrum of the remnant very hot stellar core, is a spectacle of beauty and science. In this chapter, we show that these objects are invaluable laboratories for astrophysics, astrochemistry, and astromineralogy studies, with impact in many areas of Astronomy., Comment: 37 pages, 25 figures; this is a pre-print of a chapter for the Encyclopedia of Astrophysics (edited by I. Mandel, section editor F.R.N. Schneider) to be published by Elsevier as a Reference Module

The Helix is a visually striking and the nearest planetary nebula, yet any companions responsible for its asymmetric morphology have yet to be identified. In 2020, low-amplitude photometric variations with a periodicity of 2.8 d were reported based on Cycle 1 TESS observations. In this work, with the inclusion of two additional sectors, these periodic light curves are compared with lcurve simulations of irradiated companions in such an orbit. Based on the light curve modelling, there are two representative solutions: i) a Jupiter-sized body with 0.102 Rsol and an arbitrarily small orbital inclination i=1 deg, and ii) a 0.021 Rsol exoplanet with i approx. 25 deg, essentially aligned with the Helix Nebular inclination. Irradiated substellar companion models with equilibrium temperature 4970 K are constructed and compared with existing optical spectra and infrared photometry, where Jupiter-sized bodies can be ruled out, but companions modestly larger than Neptune are still allowed. Additionally, any spatially-unresolved companions are constrained based on the multi-wavelength, photometric spectral energy distribution of the central star. No ultracool dwarf companion earlier than around L5 is permitted within roughly 1200 au, leaving only faint white dwarfs and cold brown dwarfs as possible surviving architects of the nebular asymmetries. While a planetary survivor is a tantalizing possibility, it cannot be ruled out that the light curve modulation is stellar in nature, where any substellar companion requires confirmation and may be possible with JWST observations., Comment: 9 pages, 7 figures, 1 table, accepted to MNRAS

As mass-losing asymptotic giant branch (AGB) stars evolve to planetary nebulae (PNe), the mass outflow geometries transform from nearly spherical to extreme aspherical. The physical mechanisms governing this transformation are widely believed to be linked to binarity and the associated production of disks and fast jets during transitional (post-AGB) evolutionary stages. We are carrying out a systematic ALMA survey ($P$re-planet$A$ry $N$ebulae high-angular-res$O$lution su$R$vey with $A$L$MA$ or PANORAMA) of a representative sample of bipolar and multipolar post-AGB objects. We have obtained high angular-resolution (0".1-0".4) observations of the CO(3--2) and/or 6--5 emission in order to probe the spatio-kinematic structure of the collimated outflows and the central disk/torii. The results are remarkable, generally showing the presence of bipolar or multipolar high-velocity outflows, dense toroidal waists, and in one case, a geometrically-thin circular ring around the central bipolar nebula. A high degree of point-symmetry characterizes the morphology of the mass ejecta. In this contribution, we present these and other highlights from our survey. We aim to use 2D/3D radiative transfer modeling in order to derive accurate outflow momenta, masses and mass-loss rates for our sample, and build hydrodynamical models that can explain the observed spatio-kinematic structures. These results will then be used to distinguish between different classes of PN-shaping binary interaction models.

We performed numerical simulations of the common envelope (CE) interaction between two intermediate-mass asymptotic giant branch (AGB) stars and their low-mass companions. For the first time, formation and growth of dust in the envelope is calculated explicitly. We find that the first dust grains appear as early as $\sim$1-3 yrs after the onset of the CE, and are smaller than grains formed later. As the simulations progress, a high-opacity dusty shell forms, resulting in the CE photosphere being up to an order of magnitude larger than it would be without the inclusion of dust. At the end of the simulations, the total dust yield is $0.0082~M_{\odot}$ ($0.022~M_{\odot}$) for a CE with a $1.7~M_{\odot}$ ($3.7~M_{\odot}$) AGB star. Dust formation does not substantially lead to more mass unbinding or substantially alter the orbital evolution., Comment: 6 pages, 6 figures, contributed talk to the Proceedings IAU Symposium No. 384

We study the properties of double white dwarf (DWD) mergers by performing hydrodynamic simulations using the new and improved adaptive mesh refinement code Octo-Tiger. We follow the orbital evolution of DWD systems of mass ratio q=0.7 for tens of orbits until and after the merger to investigate them as a possible origin for R Coronae Borealis (RCB) type stars. We reproduce previous results, finding that during the merger, the Helium WD donor star is tidally disrupted within 20-80 minutes since the beginning of the simulation onto the accretor Carbon-Oxygen WD, creating a high temperature shell around the accretor. We investigate the possible Helium burning in this shell and the merged object's general structure. Specifically, we are interested in the amount of Oxygen-16 dredged-up to the hot shell and the amount of Oxygen-18 produced. This is critical as the discovery of very low Oxygen-16 to Oxygen-18 ratios in RCB stars pointed out the merger scenario as a favorable explanation for their origin. A small amount of hydrogen in the donor may help keep the Oxygen-16 to Oxygen-18 ratios within observational bounds, even if moderate dredge-up from the accretor occurs. In addition, we perform a resolution study to reconcile the difference found in the amount of Oxygen-16 dredge-up between smoothed-particle hydrodynamics and grid-based simulations., Comment: 32 pages, Accepted to MNRAS

We present Submillimeter Array (SMA) mapping of $^{12}$CO $J=2\rightarrow 1$, $^{13}$CO $J=2\rightarrow 1$, and CN $N=2\rightarrow 1$ emission from the Ring-like planetary nebula (PN) NGC 3132, one of the subjects of JWST Early Release Observation (ERO) near-infrared imaging. The $\sim$5$''$ resolution SMA data demonstrate that the Southern Ring's main, bright, molecule-rich ring is indeed an expanding ring, as opposed to a limb-brightened shell, in terms of its intrinsic (physical) structure. This suggests that NGC 3132 is a bipolar nebula viewed more or less pole-on (inclination $\sim$15--30$^\circ$). The SMA data furthermore reveal that the nebula harbors a second expanding molecular ring that is aligned almost orthogonally to the main, bright molecular ring. We propose that this two-ring structure is the remnant of an ellipsoidal molecular envelope of ejecta that terminated the progenitor star's asymptotic giant branch evolution and was subsequently disrupted by a series of misaligned fast, collimated outflows or jets resulting from interactions between the progenitor and one or more companions., Comment: 20 pages, 12 figures; accepted by The Astrophysical Journal

We performed numerical simulations of the common envelope (CE) interaction between thermally-pulsing asymptotic giant branch (AGB) stars of 1.7~\Msun and 3.7~\Msun, respectively, and a 0.6~\Msun compact companion. We use tabulated equations of state to take into account recombination energy. For the first time, formation and growth of dust is calculated explicitly, using a carbon dust nucleation network with a C/O abundance ratio of 2.5 (by number). The first dust grains appear within $\sim$1--3~yrs after the onset of the CE, forming an optically thick shell at $\sim$10--20~au, growing in thickness and radius to values of $\sim$400--500~au over $\sim$40~yrs, with temperatures around 400~K. Most dust is formed in unbound material, having little effect on mass ejection or orbital evolution. By the end of the simulations, the total dust yield is $\sim8.4\times10^{-3}$~\Msun and $\sim2.2\times10^{-2}$~\Msun for the CE with a 1.7~\Msun and a 3.7~\Msun AGB star, respectively, corresponding to a nucleation efficiency close to 100\%, if no dust destruction mechanism is considered. Despite comparable dust yields to single AGB stars, \textit{in CE ejections the dust forms a thousand times faster, over tens of years as opposed to tens of thousands of years}. This rapid dust formation may account for the shift in the infrared of the spectral energy distribution of some optical transients known as luminous red novae. Simulated dusty CEs support the idea that extreme carbon stars and "water fountains" may be objects observed after a CE event., Comment: 18 pages, 17 figures, 2 tables

We carried out 3D smoothed particle hydrodynamics simulations of the common envelope binary interaction using the approximation of Bowen to calculate the dust opacity in order to investigate the resulting dust-driven accelerations. We have simulated two types of binary star: a 1.7 and a 3.7 $M_{\odot}$ thermally-pulsating, asymptotic giant branch stars with a 0.6 $M_{\odot}$ companion. We carried out simulations using both an ideal gas and a tabulated equations of state, with the latter considering the recombination energy of the envelope. We found that the dust-driven wind leads to a relatively small increase in the unbound gas, with the effect being smaller for the tabulated equation of state simulations and for the more massive primary. Dust acceleration does contribute to envelope expansion with only a slightly elongated morphology, if we believe the results from the tabulated equation of state as more reliable. The Bowen opacities in the outer envelopes of the two models, at late times, are large enough that the photosphere of the post-inspiral object is about ten times larger compared to the same without accounting for the dust opacities. As such, the prediction of the appearance of the transient would change substantially if dust is included., Comment: 14 Pages, 13 Figures. Accepted for publication in MNRAS

Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems., Comment: 32 pages, 5 figures for the main article. 12 pages 8 figures for the supplementary material

When the primary star in a close binary system evolves into a giant and engulfs its companion, its core and the companion temporarily orbit each other inside a common envelope. Drag forces transfer orbital energy and angular momentum to the envelope material. Depending on the efficiency of this process, the envelope may be ejected leaving behind a tight remnant binary system of two stellar cores, or the cores merge retaining part of the envelope material. The exact outcome of common-envelope evolution is critical for in the formation of X-ray binaries, supernova progenitors, the progenitors of compact-object mergers that emit detectable gravitational waves, and many other objects of fundamental astrophysical relevance. The wide ranges of spatial and temporal timescales that characterize common-envelope interactions and the lack of spatial symmetries present a substantial challenge to generating consistent models. Therefore, these critical phases are one of the largest sources for uncertainty in classical treatments of binary stellar evolution. Three-dimensional hydrodynamic simulations of at least part of the common envelope interaction are the key to gain predictive power in modeling common-envelope evolution. We review the development of theoretical concepts and numerical approaches for such three-dimensional hydrodynamic simulations. The inherent multi-physics, multi-scale challenges have resulted in a wide variety of approximations and numerical techniques to be exercised on the problem. We summarize the simulations published to date and their main results. Given the recent rapid progress, a sound understanding of the physics of common-envelope interactions is within reach and thus there is hope that one of the remaining fundamental problems of stellar astrophysics may be solved before long., Comment: submitted to Living Reviews in Computational Astrophysics, comments welcome (can be taken into account in the resubmission if received within 3 days after this posting)

We focus on post-asymptotic giant branch (post-AGB) binaries and study the interaction between the different components of these complex systems. These components comprise the post-AGB primary, a main sequence secondary, a circumbinary disk, as well as a fast bipolar outflow (jet) launched by the companion. We obtained well-sampled time series of high resolution optical spectra over the last decade and these spectra provide the basis of our study. The jet is detected in absorption, at superior conjunction, when the line of sight towards the primary goes through the bipolar cone. Our spectral time series scan the jets during orbital motion. Our spatio-kinematic model is constrained by these dynamical spectra. We complement this with a radiative-transfer model in which the Balmer series are used to derive total mass-loss rates in the jets. The jets are found to be wide and display an angle-dependent density structure with a dense and slower outer region near the jet cone and a fast inner part along the jet symmetry axes. The deprojected outflow velocities confirm that the companions are main sequence companions. The total mass-loss rates are large (10^{-8} and 10^{-5}\,solar mass per year), from which we can infer that the mass-accretion rates onto the companion star must be high as well. The circumbinary disk is likely the main source for the accretion disk around the companion. All systems with full disks that start near the sublimation radius show jets, whereas for systems with evolved transition disks, this lowers to a detection rate of 50%. Objects without an infrared excess do not show jets. We conclude that jet creation in post-AGB binaries is a mainstream process. The interaction between the circumbinary disks and the central binary provide the needed accretion flow, but the presence of a circumbinary disk does not seem to be the only prerequisite to launch a jet., Comment: 17 pages, 12 figures, accepted by Astronomy and Astrophysics

At least one in five of all planetary nebulae are the product of a common envelope (CE) interaction, where the companion in-spirals into the envelope of an asymptotic giant branch (AGB) star ejecting the nebula and leaving behind a compact binary. In this work we carry out 3D smoothed particle hydrodynamics simulations of the CE interaction between a 1.7 $M_{\odot}$ AGB star and a 0.6 $M_{\odot}$ companion. We model the AGB structure using a 1D stellar model taken at the seventh thermal pulse. The interaction takes place when the giant is on the expanding phase of the seventh thermal pulse and has a radius of 250 $R_{\odot}$. The post-CE orbital separations varies between 20 and 31 $R_{\odot}$, with the inclusion of recombination energy resulting in wider separations. Based on the observed short in-spiral time-scales, we suggest that thermal pulses can trigger CEs, extending the ability of AGB stars to capture companions into CEs, that would lead to the prediction of a larger population of post-AGB, post-CE binaries. Simulations that include a tabulated equation of state unbind a great deal more gas, likely unbinding the entire envelope on short time-scales. The shape of the CE after the in-spiral is more spherical for AGB than red giant branch stars, and even more so if recombination energy is included. We expect that the planetary nebula formed from this CE will have different features from those predicted by Zou et al. 2020., Comment: 17 pages excluding appendices, 17 figures, Accepted on MNRAS

We perform 3D hydrodynamical simulations of a common-envelope event involving a 12 solar mass red supergiant donor. Massive stars are expected to be qualitatively different from low-mass stars as their envelopes have significant support from radiation pressure, which increases both the final separation and amount of mass ejected through the common-envelope interaction. We perform adiabatic simulations that include radiation energy through the equation of state, which results in ejecting 60 per cent more mass (up to two thirds of the total envelope mass becoming unbound, or more) and yield a 10 per cent larger final separation compared to simulations that assume an ideal gas. When also including recombination energy, we find that at least three quarters of the envelope, and possibly the entire envelope, may be unbound. The final separation further increases by almost 20 per cent. The additional amount of ejected material is mainly due to energy injected from helium recombination. Hydrogen recombination plays a comparatively small role, as it mainly occurs in gas that has already become unbound. We conclude that the internal energy of the envelope can be a significant energy source for ejecting the common envelope, but ultimately radiation transport and convection need to be included., Comment: Accepted for publication in MNRAS, 17 pages excluding appendices, 15 figures

We present the identification of 34 likely binary central stars (CSs) of planetary nebulae (PNe) from {\it Kepler/K2} data, seven of which show eclipses. Of these, 29 are new discoveries. Two additional CSs with more complicated variability are also presented. We examined the light curves of all `possible', `likely' and `true' PNe in every {\it Kepler/K2} campaign (0 through 19) to identify CS variability that may indicate a binary CS. For Campaigns 0, 2, 7, 15, and 16 we find 6 likely or confirmed variables among 21 PNe. Our primary effort, though, was focused on Campaign 11 which targeted a Galactic bulge field containing approximately 183 PNe, in which we identified 30 candidate variable CSs. The periods of these variables range from 2.3~h to 30~d, and based on our analysis, most are likely to be close binary star systems. We present periods and preliminary classifications (eclipsing, double degenerate, or irradiated systems) for the likely binaries based on light curve shape. From our total sample of 204 target PNe, with a correction for incompleteness due to magnitude limits, we calculate a binary fraction of PN central stars to be 20.7 percent for all the observed PNe, or 23.5 percent if we limit our sample only to `true' PNe. However these fractions are almost certainly lower limits due to the large angular size of the \emph{Kepler} pixels, which leads to reduced sensitivity in detecting variability, primarily as a result of dilution and noise from the nebula and neighbouring stars. We discuss the binary population of CSs based on these results as part of the total known sample of close binary CSs., Comment: 21 pages, 11 figures, 4 tables. Accepted for publication in MNRAS

OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs or main sequence stars. The gravity solver conserves angular momentum to machine precision, thanks to a correction algorithm. This code uses HPX parallelization, allowing the overlap of work and communication and leading to excellent scaling properties, allowing for the computation of large problems in reasonable wall-clock times. In this paper, we investigate the code performance and precision by running benchmarking tests. These include simple problems, such as the Sod shock tube, as well as sophisticated, full, white-dwarf binary simulations. Results are compared to analytic solutions, when known, and to other grid based codes such as FLASH. We also compute the interaction between two white dwarfs from the early mass transfer through to the merger and compare with past simulations of similar systems. We measure OCTO-TIGERs scaling properties up to a core count of 80,000, showing excellent performance for large problems. Finally, we outline the current and planned areas of development aimed at tackling a number of physical phenomena connected to observations of transients., Comment: 38 pages, 24 figures, Co-Lead Authors: Dominic C. Marcello and Sagiv Shiber

Jets are a commonly observed phenomenon in post-asymptotic giant branch (post-AGB) binaries. Due to the orbital motion of the binary, the jet causes variable absorption in the Balmer profiles. In previous work, we have developed spatio-kinematic and radiative transfer models to reproduce the observed Balmer line variability and derive the spatio-kinematic structure of the jet and its mass-loss rate. Here, we apply our jet model to five post-AGB binaries with distinct H{\alpha} line variability and diverse orbital properties. Our models fit the H{\alpha} line variations very well. We estimate jet mass-loss rates between 10-8 Mdot yr-1 and 10-4 Mdot yr-1 , from which we deduce accretion rates onto the companion between 10-7 Mdot yr-1 and 10-3 Mdot yr-1 . These accretion rates are somewhat higher than can be comfortably explained with reasonable sources of accretion, but we argue that the circumbinary disc in these systems is most-likely the source feeding the accretion, although accretion from the post-AGB star cannot be ruled out. The diversity of the variability in the five objects is due to their wide ejection cones combined with a range of viewing angles, rather than inherent differences between the objects. The nature of the observations does not let us easily distinguish which jet launching model (stellar jet, disc wind, or X-wind) should be favoured. In conclusion, we show that our jet model includes the physical parameters to successfully reproduce the H{\alpha} line variations and retrieve the structure and mass-loss rates of the jet for all five objects that are representative of the diverse sample of Galactic post-AGB binaries., Comment: 20 pages, 9 figures, accepted for publication in MNRAS

Aims. In this study, we determine the morphology and mass-loss rate of jets emanating from the companion in post-asymptotic giant branch (post-AGB) binary stars with a circumbinary disk. In doing so, we also determine the mass accretion rates on to the companion and investigate the source feeding the circum-companion accretion disk. Methods. We perform a spatio-kinematic modelling of the jet of two well-sampled post-AGB binaries, BD+46442 and IRAS19135+3937, by fitting the orbital phased time series of H-alpha spectra. Once the jet geometry, velocity and scaled density structure are computed, we carry out radiative transfer modelling of the jet for the first four Balmer lines to determine the jet densities, thus allowing us to compute the jet mass-loss rates and mass accretion rates. Results. The spatio-kinematic model of the jet reproduces the observed absorption feature in the H-alpha lines. In both objects, the jets have an inner region with extremely low density. Using our radiative transfer model, we find the full three-dimensional density structure of both jets. From these results, we can compute mass-loss rates of the jets, which are of the order of 10^-7 - 10^-5 M_sol/yr. We estimate mass accretion rates onto the companion of 10^-6 - 10^-4 M_sol/yr. Conclusions. Based on the mass accretion rates found for these two objects, we conclude that the circumbinary disk is most likely the source feeding the circum-companion accretion disk. This is in agreement with the observed depletion patterns in post-AGB binaries. The high accretion rates from the circumbinary disk imply that the lifetime of the disk will be short. Mass-transfer from the post-AGB star cannot be excluded in these systems, but it is unlikely to provide a sufficient mass-transfer rate to sustain the observed jet mass-loss rates., Comment: 16 pages, 13 figures, accepted for publication in A&A

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. 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., Comment: 15 pages, 8 figures; accepted by MNRAS

The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes., Comment: 17 pages, 16 figures

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., Comment: 18 pages, 11 figures. We welcome any comments or questions

Aims. We aim to determine the geometry, density gradient, and velocity structure of jets in post-asymptotic giant branch (post-AGB) binaries. Methods. Our high cadence time series of high-resolution optical spectra of jet-creating post-AGB binary systems provide us with a unique tomography of the jet. We determine the spatio-kinematic structure of the jets based on these data by fitting the synthetic spectral line profiles created by our model to the observed, orbital phase-resolved, H{\alpha}-line profiles of these systems. The fitting routine is provided with an initial spectrum and is allowed to test three configurations, derived from three specific jet launching models: a stellar jet launched by the star, an X-wind, and a disk wind configuration. We apply a Markov-chain Monte Carlo routine in order to fit our model to the observations. Our fitting code is tested on the post-AGB binary IRAS19135+3937. Results. We find that a model using the stellar jet configuration gives a marginally better fit to our observations. The jet has a wide half-opening angle of about 76 degrees and reaches velocities up to 870 km/s. Conclusions. Our methodology is successful in determining some parameters for jets in post-AGB binaries. The model for IRAS19135+3937 includes a transparent, low density inner region (for a half-opening angle < 40 degrees). The source feeding the accretion disk around the companion is most likely the circumbinary disk. We will apply this jet fitting routine to other jet-creating post-AGB stars in order to provide a more complete description of these objects., Comment: 14 pages, 12 figures, accepted for publication in A&A

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., Comment: Accepted by MNRAS

Carbon stars (with C/O> 1) were long assumed to all be giants, because only AGB stars dredge up significant carbon into their atmospheres. The case is nearly iron-clad now that the formerly mysterious dwarf carbon (dC) stars are actually far more common than C giants, and have accreted carbon-rich material from a former AGB companion, yielding a white dwarf and a dC star that has gained both significant mass and angular momentum. Some such dC systems have undergone a planetary nebula phase, and some may evolve to become CH, CEMP, or Ba giants. Recent studies indicate that most dCs are likely from older, metal-poor kinematic populations. Given the well-known anti-correlation of age and activity, dCs would not be expected to show significant X-ray emission related to coronal activity. However, accretion spin-up might be expected to rejuvenate magnetic dynamos in these post mass-transfer binary systems. We describe our Chandra pilot study of six dCs selected from the SDSS for Halpha emission and/or a hot white dwarf companion, to test whether their X-ray emission strength and spectral properties are consistent with a rejuvenated dynamo. We detect all 6 dCs in the sample, which have X-ray luminosities ranging from logLx= 28.5 - 29.7, preliminary evidence that dCs may be active at a level consistent with stars that have short rotation periods of several days or less. More definitive results require a sample of typical dCs with deeper X-ray observations to better constrain their plasma temperatures., Comment: 13 pages, 5 figures. Revised and resubmitted June 20, accepted June 21, 2019 to ApJ

Cool, evolved stars are the main source of chemical enrichment of the interstellar medium (ISM), and understanding their mass loss and structure offers a unique opportunity to study the cycle of matter in the Universe. Pulsation, convection, and other dynamic processes in cool evolved stars create an atmosphere where molecules and dust can form, including those necessary to the formation of life (e.g.~Carbon-bearing molecules). Understanding the structure and composition of these stars is thus vital to several aspects of stellar astrophysics, ranging from ISM studies to modeling young galaxies and exoplanet research. Recent modeling efforts and increasingly precise observations now reveal that our understanding of cool stars photospheric, chromospheric, and atmospheric structures is limited by inadequate knowledge of the dynamic and chemical processes at work. Here we outline promising scientific opportunities for the next decade. We identify and discuss the following main opportunities: (1) identify and model the physical processes that must be included in current 1D and 3D atmosphere models of cool, evolved stars; (2) refine our understanding of photospheric, chromospheric, and outer atmospheric regions of cool evolved stars, their properties and parameters, through high-resolution spectroscopic observations, and interferometric observations at high angular resolution; (3) include the neglected role of chromospheric activity in the mass loss process of red giant branch and red super giant stars, and understand the role played by their magnetic fields; (4) identify the important shaping mechanisms for planetary nebulae and their relation with the parent asymptotic giant branch stars., Comment: 10 pages, 2 figures, White Paper submitted to the Astronomy and Astrophysics Decadal Survey (Astro2020)

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., Comment: Accepted for publication in MNRAS

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., Comment: Accepted by MNRAS

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., Comment: Submitted to MNRAS, comments and feedback are welcomed

The morphologies of planetary nebula have long been believed to be due to wind shaping processes in which a fast wind from the central star impacts a previously ejected envelope. Asymmetries assumed to exist in the slow wind envelope lead to inertial confinement shaping the resulting interacting winds flow. We present new results demonstrating the effectiveness of Common Envelope Evolution (CEE) at producing aspherical envelopes which, when impinged upon by a spherical fast stellar wind, produce highly bipolar, jet-like outflows. We have run two simple cases using the output of a single PHANTOM SPH CEE simulation. Our work uses the Adaptive Mesh Refinement code AstroBEAR to track the interaction of the fast wind and CEE ejecta allowing us to follow the morphological evolution of the outflow lobes at high resolution in 3-D. Our two models bracket low and high momentum output fast winds. We find the interaction leads to highly collimated bipolar outflows. In addition, the bipolar morphology depends on the fast wind momentum injection rate. With this dependence comes the initiation of significant symmetry breaking between the top and bottom bipolar lobes. Our simulations, though simplified, confirm the long-standing belief that CEE can plan a major role in PPN and PN shaping., Comment: To be published in proceedings of Aspherical Planetary Nebula 7

We have analysed the full suite of Australia Telescope Compact Array data for the Stingray planetary nebula. Data were taken in the 4- to 23-GHz range of radio frequencies between 1991 and 2016. The radio flux density of the nebula generally declined during that period, but between 2013 and 2016 it shows signs of halting that decline. We produced the first spatially resolved radio images of the Stingray nebula from data taken in 2005. A ring structure, which appears to be associated with the ring seen in HST images, was visible. In addition, we found a narrow extension to the radio emission towards the eastern and western edges of the nebula. We derived the emission measure of the nebula - this decreased between 1992 and 2011, suggesting that the nebula is undergoing recombination. The radio spectral index is broadly consistent with a free-free emission mechanism, however a single data point hints that a steeper spectral index has possibly emerged since 2013, which could indicate the presence of synchrotron emission. If a non-thermal component component has emerged, such as one associated with a region that is launching a jet or outflow, we predict that it would intensify in the years to come., Comment: Accepted to MNRAS. 9 pages, 6 figures

The majority of planetary nebulae (PNe) are not spherical, and current single-star models cannot adequately explain all the morphologies we observe. This has led to the Binary Hypothesis, which states that PNe are preferentially formed by binary systems. This hypothesis can be corroborated or disproved by comparing the estimated binary fraction of all PNe central stars (CS) to that of the supposed progenitor population. One way to quantify the rate of CS binarity is to detect near infra-red (IR) excess indicative of a low-mass main sequence companion. In this paper, a sample of known PNe within data release 2 of the ongoing VPHAS+ are investigated. We give details of the method used to calibrate VPHAS+ photometry, and present the expected colours of CS and main sequence stars within the survey. Objects were scrutinized to remove PN mimics from our sample and identify true CS. Within our final sample of 7 CS, 6 had previously either not been identified or confirmed. We detected an $i$ band excess indicative of a low-mass companion star in 3 CS, including one known binary, leading us to to conclude that VPHAS+ provides the precise photometry required for the IR excess method presented here, and will likely improve as the survey completes and the calibration process finalised. Given the promising results from this trial sample, the entire VPHAS+ catalogue should be used to study PNe and extend the IR excess-tested CS sample., Comment: 21 pages, 9 figures, accepted for publication in MNRAS

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., Comment: Accepted by MNRAS

We present GADRAD, a Python module that adopts heuristic search techniques in the form of genetic algorithms, to efficiently model post-asymptotic giant branch (post- AGB) disc environments. GADRAD systematically constructs the multi-dimensional pa- parameter probability density functions that arise from the fitting of radiative transfer and geometric models to optical interferometric data products. The result provides unbiased descriptions of the object's potential morphology, component luminosities and temperatures, dust composition, disc density profiles and mass. Correlation in the estimated parameters as well as potential degeneracies are revealed. Estimated probability distributions of the post-AGB environment parameters provide insight into the shaping processes that may occur in the transition from the post-AGB to the planetary nebula phase. We test parameter recovery on simulated artificial data products of a typical post-AGB environment. We then use GADRAD to model the mid-infrared spectrum and visibilities of the Ant nebula (Mz3), taken with the Very Large Telescope Interferometer's instrument MIDI. Our result is consistent with a large dusty disc with similar parameter values to those previously found by Chesneau et al., except for a larger dust mass of $3.5^{+7.5}_{-2.2}\times10^{-5}$ M$_{\odot}$. The parameter confidence intervals determined by GADRAD, can however be relied upon to impose additional constraints on all disc and system parameters. Based on our analysis and other considerations, we tentatively suggest that Mz3 is a pre-PN ejected during a magnetic (polar) common envelope interaction, where the binary may or may not have survived at the core of the nebula., Comment: Could not address reviewer's comments in timely manner

The importance of long-period binaries on the formation and evolution of planetary nebulae is still rather poorly understood, in part due to the lack of central star systems known to comprise such long-period binaries. Here, we report on the latest results from the on-going Mercator-HERMES survey for variability in the central stars of planetary nebulae. We present a study of the central stars of NGC 1514, BD+30$^\circ$623, the spectrum of which shows features associated with a hot nebular progenitor as well as a possible A-type companion. Cross-correlation of high-resolution HERMES spectra against synthetic spectra shows the system to be a highly eccentric ($e\sim0.5$), double-lined binary with a period of $\sim$3300 days. Previous studies indicated that the cool component might be a Horizontal Branch star of mass $\sim$0.55 M$_\odot$ but the observed radial velocity amplitudes rule out such a low mass. Assuming the nebular symmetry axis and binary orbital plane are perpendicular, the data are more consistent with a post-main-sequence star ascending towards the Giant Branch. We also present the continued monitoring of the central star of LoTr 5, HD 112313, which has now completed one full cycle, allowing the orbital period (P$\sim$2700 days) and eccentricity ($e\sim0.3$) to be derived. To date, the orbital periods of BD+30$^\circ$623 and HD 112313 are the longest to have been measured spectroscopically in the central stars of planetary nebulae. Furthermore, these systems, along with BD+33$^\circ$2642, comprise the only spectroscopic wide-binary central stars currently known., Comment: 4 pages, 4 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics Letters

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., Comment: Accepted by The Astrophysical Journal Supplement. 22 pages, 16 figures, 1 appendix

The globular cluster (GC) system of our Galaxy contains four planetary nebulae (PNe): K 648 (or Ps 1) in M15, IRAS 18333-2357 in M22, JaFu 1 in Pal 6, and JaFu 2 in NGC 6441. Because single-star evolution at the low stellar mass of present-epoch GCs was considered incapable of producing visible PNe, their origin presented a puzzle. We imaged the PN JaFu 1 with the Hubble Space Telescope (HST) to obtain photometry of its central star (CS) and high-resolution morphological information. We imaged IRAS 18333-2357 with better depth and resolution, and we analyzed its archival HST spectra to constrain its CS temperature and luminosity. All PNe in Galactic GCs now have high-quality HST data, allowing us to improve CS mass estimates. We find reasonably consistent masses between 0.53 and 0.58 Msun for all four objects, though estimates vary when adopting different stellar evolutionary calculations. The CS mass of IRAS 18333-2357, though, depends strongly on its temperature, which remains elusive due to reddening uncertainties. For all four objects, we consider their CS and nebular masses, their morphologies, and other incongruities to assess the likelihood that these objects formed from binary stars. Although generally limited by uncertainties (~0.02 Msun) in post-AGB tracks and core mass vs. luminosity relations, the high-mass CS in K 648 indicates a binary origin. The CS of JaFu 1 exhibits compact bright [O III] and Halpha emission, like EGB 6, suggesting a binary companion or disk. Evidence is weaker for a binary origin of JaFu 2., Comment: Accepted by Astrophysical Journal

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., Comment: 4 pages, 3 figures, proceedings of the IAU Symposium 323, Planetary Nebulae: Multiwavelength Probes of Stellar and Galactic Evolution

Stellar astrophysicists are increasingly taking into account the effects of orbiting companions on stellar evolution. New discoveries, many thanks to systematic time-domain surveys, have underlined the role of binary star interactions in a range of astrophysical events, including some that were previously interpreted as due uniquely to single stellar evolution. Here, we review classical binary phenomena such as type Ia supernovae, and discuss new phenomena such as intermediate luminosity transients, gravitational wave-producing double black holes, or the interaction between stars and their planets. Finally, we examine the reassessment of well-known phenomena in light of interpretations that include both single and binary stars, for example supernovae of type Ib and Ic or luminous blue variables. At the same time we contextualise the new discoveries within the framework and nomenclature of the corpus of knowledge on binary stellar evolution. The last decade has heralded an era of revival in stellar astrophysics as the complexity of stellar observations is increasingly interpreted with an interplay of single and binary scenarios. The next decade, with the advent of massive projects such as the Large Synoptic Survey Telescope, the Square Kilometre Array, the James Webb Space Telescope and increasingly sophisticated computational methods, will see the birth of an expanded framework of stellar evolution that will have repercussions in many other areas of astrophysics such as galactic evolution and nucleosynthesis., Comment: 40 pages, 18 figures accepted for publication as a Dawes Review of the Publications of the Astronomical Society of Australia

A current issue in the study of planetary nebulae with close binary central stars is the extent to which the binaries affect the shaping of the nebulae. Recent studies have begun to show a high coincidence rate between nebulae with large-scale axial or point symmetries and close binary stars. In addition, combined binary-star and spatio-kinematic modeling of the nebulae have demonstrated that all of the systems studied to date appear to have their central binary axis aligned with the primary axis of the nebula. Here we add two more systems to the list, the central stars and nebulae of NGC 6337 and Sp 1. We show both systems to be low inclination, with their binary axis nearly aligned with our line-of-sight. Their inclinations match published values for the inclinations of their surrounding nebulae. Including these two systems with the existing sample statistically demonstrates a direct link between the central binary and the nebular morphology. In addition to the systems' inclinations we give ranges for other orbital parameters from binary modeling, including updated orbital periods for the binary central stars of NGC 6337 and Sp 1., Comment: Accepted for publication in the Astrophysical Journal. In the previously uploaded version of this paper, Figure 7 was displayed rotated 180 degrees. That and minor additional editorial items have been corrected

Six planetary nebulae (PN) are known in the Kepler space telescope field of view, three newly identified. Of the 5 central stars of PN with useful Kepler data, one, J193110888+4324577, is a short-period, post common envelope binary exhibiting relativistic beaming effects. A second, the central star of the newly identified PN Pa5, has a rare O(He) spectral type and a periodic variability consistent with an evolved companion, where the orbital axis is almost aligned with the line of sight. The third PN, NGC~6826 has a fast rotating central star, something that can only be achieved in a merger. Fourth, the central star of the newly identified PN Kn61, has a PG1159 spectral type and a mysterious semi-periodic light variability which we conjecture to be related to the interplay of binarity with a stellar wind. Finally, the central star of the circular PN A61 does not appear to have a photometric variability above 2 mmag. With the possible exception of the variability of Kn61, all other variability behaviour, whether due to binarity or not, would not easily have been detected from the ground. We conclude, based on very low numbers, that there may be many more close binary or close binary products to be discovered with ultra-high precision photometry. With a larger number of high precision photometric observations we will be able to determine how much higher than the currently known 15 per cent, the short period binary fraction for central stars of PN is likely to be., Comment: 17 pages, 12 figures

We have used the tidal equations of Zahn to determine the maximum orbital distance at which companions are brought into Roche lobe contact with their giant primary, when the primary expands during the giant phases. This is a key step when determining the rates of interaction between giants and their companions. Our stellar structure calculations are presented as maximum radii reached during the red and asymptotic giant branch (RGB and AGB, respectively) stages of evolution for masses between 0.8 and 4.0 Mo (Z=0.001 - 0.04) and compared with other models to gauge the uncertainty on radii deriving from details of these calculations. We find overall tidal capture distances that are typically 1-4 times the maximum radial extent of the giant star, where companions are in the mass range from 1 Jupiter mass to a mass slightly smaller than the mass of the primary. We find that only companions at initial orbital separations between ~320 and ~630 Ro will be typically captured into a Roche lobe-filling interaction or a common envelope on the AGB. Comparing these limits with the period distribution for binaries that will make PN, we deduce that in the standard scenario where all ~1-8 Mo stars make a PN, at most 2.5 per cent of all PN should have a post-common envelope central star binary, at odds with the observational lower limit of 15-20 per cent. The observed over-abundance of post-interaction central stars of PN cannot be easily explained considering the uncertainties. We examine a range of explanations for this discrepancy., Comment: 19 pages, 16 figures, accepted by Monthly Notices of the Royal Astronomical Society

The common envelope interaction is thought to be the gateway to all evolved compact binaries and mergers. Hydrodynamic simulations of the common envelope interaction between giant stars and their companions are restricted to the dynamical, fast, in-spiral phase. They find that the giant envelope is lifted during this phase, but remains mostly bound to the system. At the same time, the orbital separation is greatly reduced, but in most simulations it levels off? at values larger than measured from observations. We conjectured that during the post-in-spiral phase the bound envelope gas will return to the system. Using hydrodynamic simulations, we generate initial conditions for our simulation that result in a fall-back disk with total mass and angular momentum in line with quantities from the simulations of Passy et al. We find that the simulated fall-back event reduces the orbital separation efficiently, but fails to unbind the gas before the separation levels off once again. We also find that more massive fall-back disks reduce the orbital separation more efficiently, but the efficiency of unbinding remains invariably very low. From these results we deduce that unless a further energy source contributes to unbinding the envelope (such as was recently tested by Nandez et al.), all common envelope interactions would result in mergers. On the other hand, additional energy sources are unlikely to help, on their own, to reduce the orbital separation. We conclude by discussing our dynamical fall-back event in the context of a thermally-regulated post-common envelope phase., Comment: 12 pages, 12 pages, Accepted to MNRAS

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., Comment: Accepted by MNRAS

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., Comment: Accepted for publication in MNRAS

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., Comment: accepted for publication in MNRAS

A growing number of close binary stars are being discovered among central stars of planetary nebulae. Recent and ongoing surveys are finding new systems and contributing to our knowledge of the evolution of close binary systems. The push to find more systems was largely based on early discoveries which suggested that 10 to 15% of all central stars are close binaries. One goal of this series of papers is confirmation and classification of these systems as close binaries and determination of binary system parameters. Here we provide time-resolved multi-wavelength photometry of the central star of Abell 65 as well as further analysis of the nebula and discussion of possible binary--nebula connections. Our results for Abell 65 confirm recent work showing that it has a close, cool binary companion, though several of our model parameters disagree with the recently published values. With our longer time baseline of photometric observations from 1989--2009 we also provide a more precise orbital period of 1.0037577 days., Comment: Accepted for publication in the Astronomical Journal