The peculiar molecular envelope around the post-AGB star IRAS 08544--4431

Circumbinary disks have been hypothesized to exist around a number of binary post-AGB stars. Although most of the circumbinary disks have been inferred through the near IR excess, a few of them are strong emitters of molecular emission. Here we present high angular resolution observations of the emission of $^{12}$CO and its isotopomer $^{13}$CO J=2--1 line from the circumstellar envelope around the binary post-AGB star IRAS 08544$-$4431, which is one of the most prominent members of this class of objects. We find that the envelope is resolved in our observations and two separate components can be identified: (a) a central extended and strong component with very narrow linewidth between 2 - 6 \kms; (b) a weak bipolar outflow with expansion velocity up to 8 \kms. The central compact component possesses low and variable $^{12}$CO/$^{13}$CO J=2--1 line ratio, indicating optically thick emission of the main isotope. We estimate a molecular gas mass of 0.0047 M$_\odot$ for this component based on the optically thinner $^{13}$CO J=2--1 line. We discuss the relation of the molecular envelope and the circumbinary disk inferred from near IR excess and compare with other known cases where the distribution of molecular gas has been imaged at high angular resolution.Comment: 14 pages, 4 figures. Accepted for publication in Astrophysical Journa

The Mass-Loss Induced Eccentric Kozai Mechanism: A New Channel for the Production of Close Compact Object-Stellar Binaries

Over a broad range of initial inclinations and eccentricities an appreciable fraction of hierarchical triple star systems with similar masses are essentially unaffected by the Kozai-Lidov mechanism (KM) until the primary in the central binary evolves into a compact object. Once it does, it may be much less massive than the other components in the ternary, enabling the "eccentric Kozai mechanism (EKM):" the mutual inclination between the inner and outer binary can flip signs driving the inner binary to very high eccentricity, leading to a close binary or collision. We demonstrate this "Mass-loss Induced Eccentric Kozai" (MIEK) mechanism by considering an example system and defining an ad-hoc minimal separation between the inner two members at which tidal affects become important. For fixed initial masses and semi-major axes, but uniform distributions of eccentricity and cosine of the mutual inclination, ~10% of systems interact tidally or collide while the primary is on the MS due to the KM or EKM. Those affected by the EKM are not captured by earlier quadrupole-order secular calculations. We show that fully ~30% of systems interact tidally or collide for the first time as the primary swells to AU scales, mostly as a result of the KM. Finally, ~2% of systems interact tidally or collide for the first time after the primary sheds most of its mass and becomes a WD, mostly as a result of the MIEK mechanism. These findings motivate a more detailed study of mass-loss in triple systems and the formation of close NS/WD-MS and NS/WD-NS/WD binaries without an initial common envelope phase.Comment: 12 pages, 6 figures, 1 table. Accepted for publication in ApJ. For a brief video explaining this paper, see http://youtu.be/4CdTOF17q5

The early blast wave of the 2010 explosion of U Scorpii

Three-dimensional hydrodynamic simulations exploring the first 18 hours of the 2010 January 28 outburst of the recurrent nova U Scorpii have been performed. Special emphasis was placed on capturing the enormous range in spatial scales in the blast. The pre-explosion system conditions included the secondary star and a flared accretion disk. These conditions can have a profound influence on the evolving blast wave. The blast itself is shadowed by the secondary star, which itself gives rise to a low-temperature bow-shock. The accretion disk is completely destroyed in the explosion. A model with a disk gas density of 10^{15} cm^{-3} produced a blast wave that is collimated and with clear bipolar structures, including a bipolar X-ray emitting shell. The degree of collimation depends on the initial mass of ejecta, energy of explosion, and circumstellar gas density distribution. It is most pronounced for a model with the lowest explosion energy (10^{43} erg) and mass of ejecta (10^{-8} M_{\odot}). The X-ray luminosities of three of six models computed are close to, but consistent with, an upper limit to the early blast X-ray emission obtained by the Swift satellite, the X-ray luminosity being larger for higher circumstellar gas density and higher ejecta mass. The latter consideration, together with estimates of the blast energy from previous outbursts, suggests that the mass of ejecta in the 2010 outburst was not larger than 10^{-7} M_{\odot}.Comment: 6 pages, 4 Figures; accepted for publication on ApJL. Version with full resolution images can be found at http://www.astropa.unipa.it/~orlando/PREPRINTS/u-sco.pd

Strong Variable Ultraviolet Emission from Y Gem: Accretion Activity in an AGB Star with a Binary Companion?

Binarity is believed to dramatically affect the history and geometry of mass loss in AGB and post-AGB stars, but observational evidence of binarity is sorely lacking. As part of a project to look for hot binary companions to cool AGB stars using the GALEX archive, we have discovered a late-M star, Y Gem, to be a source of strong and variable UV emission. Y Gem is a prime example of the success of our technique of UV imaging of AGB stars in order to search for binary companions. Y Gem's large and variable UV flux makes it one of the most prominent examples of a late AGB star with a mass accreting binary companion. The UV emission is most likely due to emission associated with accretion activity and a disk around a main-sequence companion star. The physical mechanism generating the UV emission is extremely energetic, with an integrated luminosity of a few L(sun) at its peak. We also find weak CO J=2-1 emission from Y Gem with a very narrow line profile (FWHM of 3.4 km/s). Such a narrow line is unlikely to arise in an outflow, and is consistent with emission from an orbiting, molecular reservoir of radius 300 AU. Y Gem may be the progenitor of the class of post-AGB stars which are binaries and possess disks but no outflows.Comment: 2 figures (Fig. 1 in color

The Red Rectangle: Its Shaping Mechanism and its Source of Ultraviolet Photons

The proto-planetary Red Rectangle nebula is powered by HD 44179, a spectroscopic binary (P = 318 d), in which a luminous post-AGB component is the primary source of both luminosity and current mass loss. Here, we present the results of a seven-year, eight-orbit spectroscopic monitoring program of HD 44179, designed to uncover new information about the source of the Lyman/far-ultraviolet continuum in the system as well as the driving mechanism for the bipolar outflow producing the current nebula. Our observations of the H-alpha line profile around the orbital phase of superior conjunction reveal the secondary component to be the origin of the fast (max. v~560$km s$^{-1}$) bipolar outflow in the Red Rectangle. The variation of total H-alpha flux from the central H II region with orbital phase also identifies the secondary or its surroundings as the source of the far-ultraviolet ionizing radiation in the system. The estimated mass of the secondary (~0.94 M$\sun$) and the speed of the outflow suggest that this component is a main sequence star and not a white dwarf, as previously suggested. We identify the source of the Lyman/far-ultraviolet continuum in the system as the hot, inner region (T$_{max} \ge 17,000$K) of an accretion disk surrounding the secondary, fed by Roche lobe overflow from the post-AGB primary at a rate of about$2 - 5\times10^{-5}$M$\sun$yr$^{-1}$. The total luminosity of the accretion disk around the secondary is estimated to be at least 300 L$\sun$, about 5% of the luminosity of the entire system. (abridged)Comment: Accepted for publication in Ap

Continuous and Burst-like Accretion onto Substellar Companions in Mira Winds

We present numerical hydrodynamical modeling of the effects of a giant planet or brown dwarf companion orbiting within the extended atmosphere and wind formation zone of an approximately solar-mass Mira variable star. The large-scale, time-dependent accretion flows within the radially oscillating and outflowing circumstellar gas around Miras are related to Bondi-Hoyle-Lyttleton flows, but have not, to our knowledge, been previously modelled. The new models presented in this paper illustrate the changes in accretion and wake dynamics as the companion mass is varied over a range from 10 to 50 Jupiter masses ($M_J$). The character of the accretion onto the companion changes greatly as the companion mass is increased. At the lowest companion masses considered here, a low continuous rate of mass accretion is punctuated by large, nearly periodic bursts of accretion. When the companion mass is large, the mass accretion has both a continuous part, and a rapidly varying, nearly stochastic part. These trends can be understood as the result of the interplay between the shocks and radial oscillations in the circumstellar gas, and the wake flow behind the companion. Models with accretion bursts may produce observable optical brightenings, and may affect SiO maser emission. (Abridged)Comment: 16 pgs., 10 figures with low resolution versions of Figs. 1, 9. Accepted for publication in MNRA

Surface Modeling to Support Small-Body Spacecraft Exploration and Proximity Operations

In order to simulate physically plausible surfaces that represent geologically evolved surfaces, demonstrating demanding surface-relative guidance navigation and control (GN&C) actions, such surfaces must be made to mimic the geological processes themselves. A report describes how, using software and algorithms to model body surfaces as a series of digital terrain maps, a series of processes was put in place that evolve the surface from some assumed nominal starting condition. The physical processes modeled in this algorithmic technique include fractal regolith substrate texturing, fractally textured rocks (of empirically derived size and distribution power laws), cratering, and regolith migration under potential energy gradient. Starting with a global model that may be determined observationally or created ad hoc, the surface evolution is begun. First, material of some assumed strength is layered on the global model in a fractally random pattern. Then, rocks are distributed according to power laws measured on the Moon. Cratering then takes place in a temporal fashion, including modeling of ejecta blankets and taking into account the gravity of the object (which determines how much of the ejecta blanket falls back to the surface), and causing the observed phenomena of older craters being progressively buried by the ejecta of earlier impacts. Finally, regolith migration occurs which stratifies finer materials from coarser, as the fine material progressively migrates to regions of lower potential energy

Wide Binary Effects on Asymmetries in Asymptotic Giant Branch Circumstellar Envelopes

Observations of increasingly higher spatial resolution reveal the existence of asymmetries in the circumstellar envelopes of a small fraction of asymptotic giant branch (AGB) stars. Although there is no general consensus for their origin, a binary companion star may be responsible. Within this framework, we investigate the gravitational effects associated with a sufficiently wide binary system, where Roche lobe overflow is unimportant, on the outflowing envelopes of AGB stars using three dimensional hydrodynamic simulations. The effects due to individual binary components are separately studied, enabling investigation of the stellar and circumstellar characteristics in detail. The reflex motion of the AGB star alters the wind velocity distribution, thereby, determining the overall shape of the outflowing envelope. On the other hand, the interaction of the companion with the envelope produces a gravitational wake, which exhibits a vertically thinner shape. The two patterns overlap and form clumpy structures. To illustrate the diversity of shapes, we present the numerical results as a function of inclination angle. Not only is spiral structure produced by the binary interaction, but arc patterns are also found that represent the former structure when viewed at different inclinations. The arcs reveal a systematic shift of their centers of curvature for cases when the orbital speed of the AGB star is comparable to its wind speed. They take on the shape of a peanut for inclinations nearly edge-on. In the limit of slow orbital motion of the AGB star relative to the wind speed, the arc pattern becomes nearly spherically symmetric. We find that the aspect ratio of the overall oblate shape of the pattern is an important diagnostic probe of the binary as it can be used to constrain the orbital velocity of the AGB star, and moreover the binary mass ratio.Comment: 33 pages, 11 figures, 4 tables, accepted for publication in the Astrophysical Journa

Radiation pressure and pulsation effects on the Roche lobe

Several observational pieces of evidence indicate that specific evolutionary channels which involve Roche lobe overflow are not correctly accounted for by the classical Roche model. We generalize the concept of Roche lobe in the presence of extra forces (caused by radiation pressure or pulsations). By computing the distortion of the equipotential surfaces, we are able to evaluate the impact of these perturbing forces on the stability of Roche-lobe overflow (RLOF). Radiative forces are parametrized through the constant reduction factor that they impose on the gravitational force from the radiating star (neglecting any shielding in case of large optical thickness). Forces imparted by pulsations are derived from the velocity profile of the wind that they trigger. We provide analytical expressions to compute the generalized Roche radius. Depending on the extra force, the Roche-lobe radius may either stay unchanged, become smaller, or even become meaningless (in the presence of a radiatively- or pulsation-driven wind). There is little impact on the RLOF stability.Comment: 11 pages, 13 Postscript figure

A Phenomenological Model for the Extended Zone Above AGB Stars

I suggest the existence of an extended zone above the surface of asymptotic giant branch (AGB), as well as similar stars experiencing high mass loss rates. In addition to the escaping wind, in this zone there are parcels of gas that do not reach the escape velocity. These parcels of dense gas rise slowly and then fall back. The wind and bound gas exist simultaneously to distances of ~100AU. I term this region the effervescent zone. In this phenomenological study I find that the density of the bound material in the effervescent zone falls as ~r^{-5/2}, not much faster than the wind density. The main motivation to propose the effervescent model is to allow wide binary companions to influence the morphology of the descendant planetary nebulae (PN) by accreting mass from the effervescent zone. Accretion from the effervescent zone is more efficient than accretion from the wind in forming an accretion disk around the companion. The companion might then blow two jets that will shape the descendant PN.Comment: New Astronomy, in pres