sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Nordhaus, Jason T.

Normalized to: Nordhaus, J.

53 article(s) in total. 174 co-authors, from 1 to 28 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2009557
Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution
Comments: 15 pages, 13 figures
Submitted: 2019-12-03
The morphology of bipolar planetary nebulae can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope evolution. Here we use the 3-D hydrodynamic adaptive mesh refinement code AsrtroBEAR to study the possibility that bipolar preplanetary nebula outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a common envelope event. The output of a single common envelope simulation via the smoothed particle hydrodynamics code PHANTOM serves as the initial conditions for our simulations. 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 common envelope binary remnant into the ejecta. We compare cases with two different values of 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 common envelope ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases, reflecting the influence of asymmetries in the turbulent dynamics of the common envelope interaction. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post common envelope winds may not be spherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers.
[2]  [pdf] - 1971328
Hydrodynamic Simulations of Disrupted Planetary Accretion Discs Inside the Core of an AGB Star
Comments: Submitted to MNRAS
Submitted: 2019-07-31
Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc.Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.
[3]  [pdf] - 1914736
Astro2020 APC White Paper: Accessible Astronomy: Policies, Practices, and Strategies to Increase Participation of Astronomers with Disabilities
Comments: 10 pages, 1 figure, Astro2020 APC White paper, State of the Profession Consideration
Submitted: 2019-07-10
(Abridged) In this white paper, we outline the major barriers to access within the educational and professional practice of astronomy. We present current best practices for inclusivity and accessibility, including classroom practices, institutional culture, support for infrastructure creation, hiring processes, and outreach initiatives. We present specific ways--beyond simple compliance with the ADA--that funding agencies, astronomers, and institutions can work together to make astronomy as a field more accessible, inclusive, and equitable. In particular, funding agencies should include the accessibility of institutions during proposal evaluation, hold institutions accountable for inaccessibility, and support efforts to gather data on the status and progress of astronomers and astronomy students with disabilities.
[4]  [pdf] - 1863939
Energy Budget and Core-Envelope Motion in Common Envelope Evolution
Comments: 17 pages, 10 figures, 5 tables, accepted for publication in MNRAS
Submitted: 2018-12-28, last modified: 2019-03-25
We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total envelope energy is nearly constant during this time because positive energy transferred to the gas from the core particles is counterbalanced by the negative binding energy from the closer proximity of the inner layers to the plunged-in secondary. During the subsequent slow spiral-in phase, energy continues to transfer to the envelope from the red giant core and secondary core particles. We also propose that relative motion between the centre of mass of the envelope and the centre of mass of the particles could account for the offsets of planetary nebula central stars from the nebula's geometric centre.
[5]  [pdf] - 1846833
The Role of Convection in Determining the Ejection Efficiency of Common Envelope Interactions
Comments: Accepted to MNRAS. 10 pages, 7 figures
Submitted: 2018-11-07, last modified: 2019-02-26
A widely used method for parameterizing the outcomes of common envelopes (CEs) involves defining an ejection efficiency, $\bar{\alpha}_{\mathrm{eff}}$, that represents the fraction of orbital energy used to unbind the envelope as the orbit decays. Given $\bar{\alpha}_{\mathrm{eff}}$, a prediction for the post-CE orbital separation is possible with knowledge of the energy required to unbind the primary's envelope from its core. Unfortunately, placing observational constraints on $\bar{\alpha}_{\mathrm{eff}}$ is challenging as it requires knowledge of the primary's structure at the onset of the common envelope phase. Numerical simulations have also had difficulties reproducing post-CE orbital configurations as they leave extended, but still bound, envelopes. Using detailed stellar interior profiles, we calculate $\bar{\alpha}_{\mathrm{eff}}$ values for a matrix of primary-companion mass pairs when the primary is at maximal extent in its evolution. We find that the ejection efficiency is most sensitive to the properties of the surface-contact convective region (SCCR). In this region, the convective transport timescales are often short compared to orbital decay timescales, thereby allowing the star to effectively radiate orbital energy and thus lower $\bar{\alpha}_{\mathrm{eff}}$. The inclusion of convection in numerical simulations of CEs may aid ejection without the need for additional energy sources as the orbit must shrink substantially further before the requisite energy can be tapped to drive ejection. Additionally, convection leads to predicted post-CE orbital periods of less than a day in many cases, an observational result that has been difficult to reproduce in population studies where $\bar{\alpha}_{\mathrm{eff}}$ is taken to be constant. Finally, we provide a simple method to calculate $\bar{\alpha}_{\mathrm{eff}}$ if the properties of the SCCR are known.
[6]  [pdf] - 2050118
Accretion in Common Envelope Evolution
Comments: 4 pages, 3 figures, to appear in the proceedings of IAU Symposium 343: Why Galaxies Care About AGB Stars, (eds.) F. Kerschbaum, M. Groenewegen and H. Olofsson
Submitted: 2018-10-10
Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.
[7]  [pdf] - 1716290
Planetary Nebulae Shaped By Common Envelope Evolution
Comments: To be published in proceedings of Aspherical Planetary Nebula 7
Submitted: 2018-07-16
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.
[8]  [pdf] - 1724963
Accretion in common envelope evolution
Comments: 14 pages, submitted to MNRAS
Submitted: 2018-05-09
Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the envelope. Additional power sources might be necessary and accretion onto the inspiraling companion is one such source. Accretion is likely common in post-asymptotic giant branch (AGB) binary interactions but how it operates and how its consequences depend on binary separation remain open questions. Here we use high resolution global 3D hydrodynamic simulations of CEE with the adaptive mesh refinement (AMR) code AstroBEAR, to bracket the range of CEE companion accretion rates by comparing runs that remove mass and pressure via a subgrid accretion model with those that do not. The results show that if a pressure release valve is available, super-Eddington accretion may be common. Jets are a plausible release valve in these environments, and they could also help unbind and shape the envelopes.
[9]  [pdf] - 1583131
Wind-accelerated orbital evolution in binary systems with giant stars
Comments: 10 pages, 10 figures and 6 tables submitted to MNRAS
Submitted: 2017-05-04, last modified: 2017-08-31
Using 3D radiation-hydrodynamic simulations and analytic theory, we study the orbital evolution of asymptotic-giant-branch (AGB) binary systems for various initial orbital separations and mass ratios, and thus different initial accretion modes. The time evolution of binary separations and orbital periods are calculated directly from the averaged mass loss rate, accretion rate and angular momentum loss rate. We separately consider spin-orbit synchronized and zero spin AGB cases. We find that the the angular momentum carried away by the mass loss together with the mass transfer can effectively shrink the orbit when accretion occurs via wind-Roche-lobe overflow. In contrast, the larger fraction of mass lost in Bondi-Hoyle-Lyttleton accreting systems acts to enlarge the orbit. Synchronized binaries tend to experience stronger orbital period decay in close binaries. We also find that orbital period decay is faster when we account for the nonlinear evolution of the accretion mode as the binary starts to tighten. This can increase the fraction of binaries that result in common envelope, luminous red novae, Type Ia supernovae and planetary nebulae with tight central binaries. The results also imply that planets in the the habitable zone around white dwarfs are unlikely to be found.
[10]  [pdf] - 1562773
Mass Transfer and Disc Formation in AGB Binary Systems
Comments: 14 pages, 31 figures and 4 tables. accepted by MNRAS
Submitted: 2017-02-20, last modified: 2017-04-13
We investigate mass transfer and the formation of disc in binary systems using a combination of numerical simulations and theory. We consider six models distinguished by binary separation, secondary mass and outflow mechanisms. Each system consists of an asymptotic-giant-branch (AGB) star and an accreting secondary. The AGB star loses its mass via a wind. In one of our six models, the AGB star incurs a short period of outburst. In all cases, the secondary accretes part of the ejected mass and also influences the mass-loss rate of the AGB star. The ejected mass may remain gravitationally bound to the binary system and form a circumbinary disk, or contribute to an accretion disk around the secondary. In other cases, the ejecta will escape the binary system. The accretion rate on to the secondary changes non-linearly with binary separation. In our closest binary simulations, our models exemplify the wind Roche lobe overflow while in our wide binary cases, the mass transfer exhibits Bondi-Hoyle accretion. The morphologies of the outflows in the binary systems are varied. The variety may provide clues to how the late AGB phase influences planetary nebulae shaping. We employ the adaptive-mesh-refinement code ASTROBEAR for our simulations and include ray-tracing, radiation transfer, cooling and dust formation. To attain the highest computational efficiency and the most stable results, all simulations are run in the corotating frame.
[11]  [pdf] - 1581315
Mass transfer in asymptotic-giant-branch binary systems
Comments: 10 figures, 2 pages, revised for IAU323 proceedings
Submitted: 2017-01-30
Binary stars can interact via mass transfer when one member (the primary) ascends onto a giant branch. The amount of gas ejected by the binary and the amount of gas accreted by the secondary over the lifetime of the primary influence the subsequent binary phenomenology. Some of the gas ejected by the binary will remain gravitationally bound and its distribution will be closely related to the formation of planetary nebulae. We investigate the nature of mass transfer in binary systems containing an AGB star by adding radiative transfer to the AstroBEAR AMR Hydro/MHD code.
[12]  [pdf] - 1580864
Post-common envelope PN, fundamental or irrelevant?
Comments: 4 pages, 3 figures, proceedings of the IAU Symposium 323, Planetary Nebulae: Multiwavelength Probes of Stellar and Galactic Evolution
Submitted: 2016-12-11
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.
[13]  [pdf] - 1490661
Three-dimensional hydrodynamic simulations of L2 Puppis
Comments: 6 pages, 12 figures, submitted to MNRAS
Submitted: 2016-02-19, last modified: 2016-10-04
Recent observations of the L2 Puppis system suggest that this Mira-like variable may be in the early stages of forming a bipolar planetary nebula (PN). As one of nearest and brightest AGB stars, thought be a binary, L2 Puppis serves as a benchmark object for studying the late-stages of stellar evolution. We perform global, three-dimensional, adaptive-mesh-refinement hydrodynamic simulations of the L2 Puppis system with AstroBEAR. We use the radiative transfer code RADMC-3D to construct the broad-band spectral-energy-distribution (SED) and synthetic observational images from our simulations. Given the reported binary parameters, we are able to reproduce the current observational data if a short pulse of dense material is released from the AGB star with a velocity sufficient to escape the primary but not the binary. Such a situation could result from a thermal pulse, be induced by a periastron passage of the secondary, or could be launched if the primary ingests a planet.
[14]  [pdf] - 1391718
Radio Emission from Red-Giant Hot Jupiters
Comments: published in ApJ, corrected typo in abstract
Submitted: 2016-01-20, last modified: 2016-04-18
When planet-hosting stars evolve off the main sequence and go through the red-giant branch, the stars become orders of magnitudes more luminous and, at the same time, lose mass at much higher rates than their main-sequence counterparts. Accordingly, if planetary companions exist around these stars at orbital distances of a few AU, they will be heated up to the level of canonical hot Jupiters and also be subjected to a dense stellar wind. Given that magnetized planets interacting with stellar winds emit radio waves, such "Red-Giant Hot Jupiters" (RGHJs) may also be candidate radio emitters. We estimate the spectral auroral radio intensity of RGHJs based on the empirical relation with the stellar wind as well as a proposed scaling for planetary magnetic fields. RGHJs might be intrinsically as bright as or brighter than canonical hot Jupiters and about 100 times brighter than equivalent objects around main-sequence stars. We examine the capabilities of low-frequency radio observatories to detect this emission and find that the signal from an RGHJ may be detectable at distances up to a few hundred parsecs with the Square Kilometer Array.
[15]  [pdf] - 1416749
The Creation of AGB Fallback Shells
Comments: 7 pages, 13 figures, accepted by MNRAS on Jan. 4th 2016
Submitted: 2015-05-12, last modified: 2016-01-22
The possibility that mass ejected during Asymptotic Giant Branch (AGB) stellar evolution phases falls back towards the star has been suggested in applications ranging from the formation of accretion disks to the powering of late-thermal pulses. In this paper, we seek to explicate the properties of fallback flow trajectories from mass-loss events. We focus on a transient phase of mass ejection with sub-escape speeds, followed by a phase of a typical AGB wind. We solve the problem using both hydrodynamic simulations and a simplified one-dimensional analytic model that matches the simulations. For a given set of initial wind characteristics, we find a critical shell velocity that distinguishes between "shell fallback" and "shell escape". We discuss the relevance of our results for both single and binary AGB stars. In particular, we discuss how our results help to frame further studies of fallback as a mechanism for forming the substantial population of observed post-AGB stars with dusty disks.
[16]  [pdf] - 972092
Proposed Nomenclature for Surface Features on Pluto and Its Satellites and Names for Newly Discovered Satellites
Comments: 6 pages, white paper
Submitted: 2015-03-26
In anticipation of the July 2015 flyby of the Pluto system by NASA's New Horizons mission, we propose naming conventions and example names for surface features on Pluto and its satellites (Charon, Nix, Hydra, Kerberos, Styx) and names for newly discovered satellites.
[17]  [pdf] - 1203409
A call for a paradigm shift from neutrino-driven to jet-driven core-collapse supernova mechanisms
Comments: final version; accepted for publication in mnras
Submitted: 2014-02-18, last modified: 2015-01-28
Three-dimensional (3D) simulations in recent years have shown severe difficulties producing 10^51 erg explosions of massive stars with neutrino based mechanisms while on the other hand demonstrated the large potential of mechanical effects, such as winds and jets in driving explosions. In this paper we study the typical time-scale and energy for accelerating gas by neutrinos in core-collapse supernovae (CCSNe) and find that under the most extremely favorable (and probably unrealistic) conditions, the energy of the ejected mass can reach at most 5X10^50 erg. More typical conditions yield explosion energies an order-of-magnitude below the observed 10^51 erg explosions. On the other hand, non-spherical effects with directional outflows hold promise to reach the desired explosion energy and beyond. Such directional outflows, which in some simulations are produced by numerical effects of 2D grids, can be attained by angular momentum and jet launching. Our results therefore call for a paradigm shift from neutrino-based explosions to jet-driven explosions for CCSNe.
[18]  [pdf] - 1223548
The Chandra Planetary Nebulae Survey (ChanPlaNS): III. X-ray Emission from the Central Stars of Planetary Nebulae
Comments: A-c-c-e-p-t-e-d by ApJ; ChanPlaNS website:; remember to tip your servers!
Submitted: 2014-12-08
We present X-ray spectral analysis of 20 point-like X-ray sources detected in Chandra Planetary Nebula Survey (ChanPlaNS) observations of 59 planetary nebulae (PNe) in the solar neighborhood. Most of these 20 detections are associated with luminous central stars within relatively young, compact nebulae. The vast majority of these point-like X-ray-emitting sources at PN cores display relatively "hard" ($\geq0.5$~keV) X-ray emission components that are unlikely to be due to photospheric emission from the hot central stars (CSPN). Instead, we demonstrate that these sources are well modeled by optically-thin thermal plasmas. From the plasma properties, we identify two classes of CSPN X-ray emission: (1) high-temperature plasmas with X-ray luminosities, $L_{\rm X}$, that appear uncorrelated with the CSPN bolometric luminosity, $L_{\rm bol}$; and (2) lower-temperature plasmas with $L_{\rm X}/L_{\rm bol}\sim10^{-7}$. We suggest these two classes correspond to the physical processes of magnetically active binary companions and self-shocking stellar winds, respectively. In many cases this conclusion is supported by corroborative multiwavelength evidence for the wind and binary properties of the PN central stars. By thus honing in on the origins of X-ray emission from PN central stars, we enhance the ability of CSPN X-ray sources to constrain models of PN shaping that invoke wind interactions and binarity.
[19]  [pdf] - 1157523
On the Orbits of Low-mass Companions to White Dwarfs and the Fates of the Known Exoplanets
Comments: Replaced with version in Journal
Submitted: 2012-11-05, last modified: 2014-09-02
The ultimate fates of binary companions to stars (including whether the companion survives and the final orbit of the binary) are of interest in light of an increasing number of recently discovered, low-mass companions to white dwarfs (WDs). In this Letter, we study the evolution of a two-body system wherein the orbit adjusts due to structural changes in the primary, dissipation of orbital energy via tides, and mass loss during the giant phases; previous studies have not incorporated changes in the primary's spin. For companions ranging from Jupiter's mass to ~0.3 Msun and primaries ranging from 1-3 Msun, we determine the minimum initial semimajor axis required for the companion to avoid engulfment by the primary during post-main-sequence evolution, and highlight the implications for the ultimate survival of the known exoplanets. We present regions in secondary mass and orbital period space where an engulfed companion might be expected to survive the common envelope phase (CEP), and compare with known M dwarf+WD short-period binaries. Finally, we note that engulfed Earth-like planets cannot survive a CEP. Detection of a first-generation terrestrial planet in the white dwarf habitable zone requires scattering from a several-AU orbit to a high-eccentricity orbit (with a periastron of ~Rsun) from which it is damped into a circular orbit via tidal friction, possibly rendering it an uninhabitable, charred ember.
[20]  [pdf] - 873260
The Chandra Planetary Nebula Survey (ChanPlaNS). II. X-ray Emission from Compact Planetary Nebulae
Comments: 28 pages, 5 figures, accepted to ApJ
Submitted: 2014-07-15, last modified: 2014-08-20
We present results from the most recent set of observations obtained as part of the Chandra X-ray observatory Planetary Nebula Survey (ChanPlaNS), the first comprehensive X-ray survey of planetary nebulae (PNe) in the solar neighborhood (i.e., within ~1.5 kpc of the Sun). The survey is designed to place constraints on the frequency of appearance and range of X-ray spectral characteristics of X-ray-emitting PN central stars and the evolutionary timescales of wind-shock-heated bubbles within PNe. ChanPlaNS began with a combined Cycle 12 and archive Chandra survey of 35 PNe. ChanPlaNS continued via a Chandra Cycle 14 Large Program which targeted all (24) remaining known compact (R_neb <~ 0.4 pc), young PNe that lie within ~1.5 kpc. Results from these Cycle 14 observations include first-time X-ray detections of hot bubbles within NGC 1501, 3918, 6153, and 6369, and point sources in HbDs 1, NGC 6337, and Sp 1. The addition of the Cycle 14 results brings the overall ChanPlaNS diffuse X-ray detection rate to ~27% and the point source detection rate to ~36%. It has become clearer that diffuse X-ray emission is associated with young (<~5x10^3 yr), and likewise compact (R_neb<~0.15 pc), PNe with closed structures and high central electron densities (n_e>~1000 cm^-3), and rarely associated with PNe that show H_2 emission and/or pronounced butterfly structures. Hb 5 is one such exception of a PN with a butterfly structure that hosts diffuse X-ray emission. Additionally, of the five new diffuse X-ray detections, two host [WR]-type CSPNe, NGC 1501 and NGC 6369, supporting the hypothesis that PNe with central stars of [WR]-type are likely to display diffuse X-ray emission.
[21]  [pdf] - 1208543
Dusty Disks Around Central Stars of Planetary Nebulae
Comments: 11 pages, 5 figures, AJ, in press
Submitted: 2014-03-23
Only a few percent of cool, old white dwarfs (WDs) have infrared excesses interpreted as originating in small hot disks due to the infall and destruction of single asteroids that come within the star's Roche limit. Infrared excesses at 24 micron were also found to derive from the immediate vicinity of younger, hot WDs, most of which are still central stars of planetary nebulae (CSPN). The incidence of CSPN with this excess is 18%. The Helix CSPN, with a 24 micron excess, has been suggested to have a disk formed from collisions of Kuiper belt-like objects (KBOs). In this paper, we have analyzed an additional sample of CSPN to look for similar infrared excesses. These CSPN are all members of the PG 1159 class and were chosen because their immediate progenitors are known to often have dusty environments consistent with large dusty disks. We find that, overall, PG 1159 stars do not present such disks more often than other CSPN, although the statistics (5 objects) are poor. We then consider the entire sample of CSPN with infrared excesses, and compare it to the infrared properties of old WDs, as well as cooler post-AGB stars. We conclude with the suggestion that the infrared properties of CSPN more plausibly derive from AGB-formed disks rather than disks formed via the collision of KBOs, although the latter scenario cannot be ruled out. We finally remark that there seems to be an association between CSPN with a 24 micron excess and confirmed or possible binarity of the central star.
[22]  [pdf] - 1943187
The Herschel Planetary Nebula Survey (HerPlaNS) I. Data Overview and Analysis Demonstration with NGC 6781
Comments: 27 pages, 16 figures, accepted for publication in A&A
Submitted: 2014-03-11
This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.
[23]  [pdf] - 1173500
Size of discs formed by wind accretion in binaries can be underestimated if the role of wind-driving force is ignored
Comments: 8 pages, 2 figures, accepted by MNRAS; in press
Submitted: 2013-08-19, last modified: 2013-09-11
Binary systems consisting of a secondary accreting form a wind-emitting primary are ubiquitous in astrophysics. The phenomenology of such Bondi-Hoyle-Lyttleton (BHL) accretors is particularly rich when an accretion disc forms around the secondary. The outer radius of such discs is commonly estimated from the net angular momentum produced by a density variation of material across the BHL or Bondi accretion cylinder, as the latter is tilted with respect to the direction to the primary due to orbital motion. But this approach has ignored the fact that the wind experiences an outward driving force that the secondary does not. In actuality, the accretion stream falls toward a retarded point in the secondary's orbit as the secondary is pulled toward the primary relative to the stream. The result is a finite separation or "accretion stream impact parameter" (ASIP) separating the secondary and stream. When the orbital radius a_o exceeds the BHL radius r_b, the ratio of outer disc radius estimated as the ASIP to the conventional estimate a_o^{1/2}/r_b^{1/2}>1. We therefore predict that discs will form at larger radii from the secondary than traditional estimates. This agrees with the importance of the ASIP emphasized by Huarte-Espinosa et al. (2013) and the practical consequence that resolving the initial outer radius of such an accretion disc in numerical simulations can be less demanding than what earlier estimates would suggest.
[24]  [pdf] - 1157608
The Formation and Evolution of Wind-Capture Disks In Binary Systems
Comments: 13 pages, 9 figures, in press, MNRAS
Submitted: 2012-11-07, last modified: 2013-05-03
We study the formation, evolution and physical properties of accretion disks formed via wind capture in binary systems. Using the AMR code AstroBEAR, we have carried out high resolution 3D simulations that follow a stellar mass secondary in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in order to properly resolve the formation of accretion disks around the secondary. We then compare simulations of binaries with three different orbital radii (10, 15, 20 AU). Disks are formed in all three cases, however the size of the disk and, most importantly, its accretion rate decreases with orbital radii. In addition, the shape of the orbital motions of material within the disk becomes increasingly elliptical with increasing binary separation. The flow is mildly unsteady with "fluttering" around the bow shock observed. The disks are generally well aligned with the orbital plane after a few binary orbits. We do not observe the presence of any large scale, violent instabilities (such as the flip-flop mode). For the first time, moreover, it is observed that the wind component that is accreted towards the secondary has a vortex tube-like structure, rather than a column-like one as it was previously thought. In the context of AGB binary systems that might be precursors to Pre-Planetary and Planetary Nebula, we find that the wind accretion rates at the chosen orbital separations are generally too small to produce the most powerful outflows observed in these systems if the companions are main sequence stars but marginally capable if the companions are white dwarfs. It is likely that many of the more powerful PPN and PN involve closer binaries than the ones considered here. The results also demonstrate principles of broad relevance to all wind-capture binary systems.
[25]  [pdf] - 1152268
Dimensional Dependence of the Hydrodynamics of Core-Collapse Supernovae
Comments: Submitted to the Astrophysical Journal
Submitted: 2012-10-18
The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of explosions. Unfortunately, much of this work has necessarily been performed assuming axisymmetry and it remains unclear whether or not this compromises those results. In this work, we present analyses of simplified two- and three-dimensional CCSN models with the goal of comparing the multidimensional hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial to understanding the explosion mechanism, others are quite dramatic and make interpreting 2D CCSN models problematic. In particular, we find that imposing axisymmetry artificially produces excess power at the largest spatial scales, power that has been deemed critical in the success of previous explosion models and has been attributed solely to the standing accretion shock instability. Nevertheless, our 3D models, which have an order of magnitude less power on large scales compared to 2D models, explode earlier. Since we see explosions earlier in 3D than in 2D, the vigorous sloshing associated with the large scale power in 2D models is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. Possibly related to the earlier explosions in 3D, we find that about 25% of the accreted material spends more time in the gain region in 3D than in 2D, being exposed to more integrated heating and reaching higher peak entropies, an effect we associate with the differing characters of turbulence in 2D and 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.
[26]  [pdf] - 1118315
The Chandra X-ray Survey of Planetary Nebulae (ChanPlaNS): Probing Binarity, Magnetic Fields, and Wind Collisions
Comments: 41 pages, 6 figures; submitted to the Astronomical Journal
Submitted: 2012-04-25, last modified: 2012-05-31
We present an overview of the initial results from the Chandra Planetary Nebula Survey (ChanPlaNS), the first systematic (volume-limited) Chandra X-ray Observatory survey of planetary nebulae (PNe) in the solar neighborhood. The first phase of ChanPlaNS targeted 21 mostly high-excitation PNe within ~1.5 kpc of Earth, yielding 4 detections of diffuse X-ray emission and 9 detections of X-ray-luminous point sources at the central stars (CSPNe) of these objects. Combining these results with those obtained from Chandra archival data for all (14) other PNe within ~1.5 kpc that have been observed to date, we find an overall X-ray detection rate of ~70%. Roughly 50% of the PNe observed by Chandra harbor X-ray-luminous CSPNe, while soft, diffuse X-ray emission tracing shocks formed by energetic wind collisions is detected in ~30%; five objects display both diffuse and point-like emission components. The presence of X-ray sources appears correlated with PN density structure, in that molecule-poor, elliptical nebulae are more likely to display X-ray emission (either point-like or diffuse) than molecule-rich, bipolar or Ring-like nebulae. All but one of the X-ray point sources detected at CSPNe display X-ray spectra that are harder than expected from hot (~100 kK) central star photospheres, possibly indicating a high frequency of binary companions to CSPNe. Other potential explanations include self-shocking winds or PN mass fallback. Most PNe detected as diffuse X-ray sources are elliptical nebulae that display a nested shell/halo structure and bright ansae; the diffuse X-ray emission regions are confined within inner, sharp-rimmed shells. All sample PNe that display diffuse X-ray emission have inner shell dynamical ages <~5x10^3 yr, placing firm constraints on the timescale for strong shocks due to wind interactions in PNe.
[27]  [pdf] - 524909
Precision Astrometry of the Exoplanet Host Candidate GD 66
Comments: 6 pages, 3 figures, 1 table, accepted to MNRAS
Submitted: 2012-05-03
The potential existence of a giant planet orbiting within a few AU of a stellar remnant has profound implications for both the survival and possible regeneration of planets during post-main sequence stellar evolution. This paper reports Hubble Space Telescope Fine Guidance Sensor and U.S. Naval Observatory relative astrometry of GD 66, a white dwarf thought to harbor a giant planet between 2 and 3 AU based on stellar pulsation arrival times. Combined with existing infrared data, the precision measurements here rule out all stellar-mass and brown dwarf companions, implying that only a planet remains plausible, if orbital motion is indeed the cause of the variations in pulsation timing.
[28]  [pdf] - 1092390
The Hydrodynamic Origin of Neutron Star Kicks
Comments: Accepted to MNRAS, replaced with accepted manuscript
Submitted: 2011-12-14, last modified: 2012-03-28
We present results from a suite of axisymmetric, core-collapse supernova simulations in which hydrodynamic recoil from an asymmetric explosion produces large proto-neutron star (PNS) velocities. We use the adaptive-mesh refinement code CASTRO to self-consistently follow core-collapse, the formation of the PNS and its subsequent acceleration. We obtain recoil velocities of up to 620 km/s at ~1 s after bounce. These velocities are consistent with the observed distribution of pulsar kicks and with PNS velocities obtained in other theoretical calculations. Our PNSs are still accelerating at several hundred km/s at the end of our calculations, suggesting that even the highest velocity pulsars may be explained by hydrodynamic recoil in generic, core-collapse supernovae.
[29]  [pdf] - 1083605
A Binary Scenario for the Formation of Strongly Magnetized White Dwarfs
Comments: To appear in the Proceedings of the 2nd International Symposium on Strong Electromagnetic Fields and Neutron Stars, Varadero, Cuba
Submitted: 2011-08-24
Since their initial discovery, the origin of isolated white dwarfs (WDs) with magnetic fields in excess of $\sim$1 MG has remained a mystery. Recently, the formation of these high-field magnetic WDs has been observationally linked to strong binary interactions incurred during post-main-sequence evolution. Planetary, brown dwarf or stellar companions located within a few AU of main-sequence stars may become engulfed during the primary's expansion off the main sequence. Sufficiently low-mass companions in-spiral inside a common envelope until they are tidally shredded near the natal white dwarf. Formation of an accretion disk from the disrupted companion provides a source of turbulence and shear which act to amplify magnetic fields and transport them to the WD surface. We show that these disk-generated fields explain the observed range of magnetic field strengths for isolated, high-field magnetic WDs. Additionally, we discuss a high-mass binary analogue which generates a strongly-magnetized WD core inside a pre-collapse, massive star. Subsequent core-collapse to a neutron star may produce a magnetar.
[30]  [pdf] - 1077050
The Dust Properties of Two Hot R Coronae Borealis Stars and a Wolf-Rayet Central Star of a Planetary Nebula: in Search of a Possible Link
Comments: 20 pages, 5 figures, accepted for publication in The Astronomical Journal
Submitted: 2011-06-03
We present new Spitzer/IRS spectra of two hot R Coronae Borealis (RCB) stars, one in the Galaxy,V348 Sgr, and one lying in the LMC, HV 2671. These two objects may constitute a link between the RCB stars and the late Wolf-Rayet ([WCL]) class of central stars of planetary nebula (CSPNe) such as CPD -56 8032 that has little or no hydrogen in their atmospheres. HV 2671 and V348 Sgr are members of a rare subclass that has significantly higher effective temperatures than most RCB stars, but sharing the traits of hydrogen deficiency and dust formation that define the cooler RCB stars. The [WC] CSPNe star, CPD -56 8032, displays evidence for dual-dust chemistry showing both PAHs and crystalline silicates in its mid-IR spectrum. HV 2671 shows strong PAH emission but shows no sign of having crystalline silicates. The spectrum of V348 Sgr is very different from those of CPD -56 8032 and HV 2671. The PAH emission seen strongly in the other two stars is not present. Instead, the spectrum is dominated by a broad emission centered at about 8.2 micron. The mid-IR spectrum of CPD -56 8032 shows emission features that may be associated with C60. The other two stars do not show evidence for C60. HV 2671 has also been detected by Herschel/PACS and SPIRE. V348 Sgr and CPD -56 8032 have been detected by AKARI/FIS. These data were combined with Spitzer, IRAS, 2MASS and other photometry to produce their spectral energy distributions from the visible to the far-IR. Monte Carlo radiative transfer modeling was used to study the circumstellar dust around these stars. HV 2671 and CPD -56 8032 require both a flared inner disk with warm dust and an extended diffuse envelope with cold dust to to fit their SEDs. The SED of V348 Sgr can be fit with a much smaller disk and envelope.
[31]  [pdf] - 1041516
Induced Rotation in 3D Simulations of Core Collapse Supernovae: Implications for Pulsar Spins
Comments: Accepted to the Astrophysical Journal
Submitted: 2010-10-25, last modified: 2011-03-09
It has been suggested that the observed rotation periods of radio pulsars might be induced by a non-axisymmetric spiral-mode instability in the turbulent region behind the stalled supernova bounce shock, even if the progenitor core was not initially rotating. In this paper, using the three-dimensional AMR code CASTRO with a realistic progenitor and equation of state and a simple neutrino heating and cooling scheme, we present a numerical study of the evolution in 3D of the rotational profile of a supernova core from collapse, through bounce and shock stagnation, to delayed explosion. By the end of our simulation ($\sim$420 ms after core bounce), we do not witness significant spin up of the proto-neutron star core left behind. However, we do see the development before explosion of strong differential rotation in the turbulent gain region between the core and stalled shock. Shells in this region acquire high spin rates that reach $\sim$$150\,$ Hz, but this region contains too little mass and angular momentum to translate, even if left behind, into rapid rotation for the full neutron star. We find also that much of the induced angular momentum is likely to be ejected in the explosion, and moreover that even if the optimal amount of induced angular momentum is retained in the core, the resulting spin period is likely to be quite modest. Nevertheless, induced periods of seconds are possible.
[32]  [pdf] - 320864
VLTI observations of the dust geometry around R Coronae Borealis stars
Comments: Accepted for publication in MNRAS; 14 pages, 10 figures, 6 tables
Submitted: 2011-02-21
We are investigating the formation and evolution of dust around the hydrogen-deficient supergiants known as R Coronae Borealis (RCB) stars. We aim to determine the connection between the probable merger past of these stars and their current dust-production activities. We carried out high-angular resolution interferometric observations of three RCB stars, namely RY Sgr, V CrA, and V854 Cen with the mid-IR interferometer, MIDI on the VLTI, using two telescope pairs. The baselines ranged from 30 to 60 m, allowing us to probe the dusty environment at very small spatial scales (~ 50 mas or 400 stellar radii). The observations of the RCB star dust environments were interpreted using both geometrical models and one-dimensional radiative transfer codes. From our analysis we find that asymmetric circumstellar material is apparent in RY Sgr, may also exist in V CrA, and is possible for V854 Cen. Overall, we find that our observations are consistent with dust forming in clumps ejected randomly around the RCB star so that over time they create a spherically symmetric distribution of dust. However, we conclude that the determination of whether there is a preferred plane of dust ejection must wait until a time series of observations are obtained.
[33]  [pdf] - 339431
Towards a New Prescription for the Tidal Capture of Planets, Brown Dwarfs and Stellar Companions
Comments: Part of PlanetsbeyondMS/2010 proceedings
Submitted: 2011-02-02
We study the evolution of the orbit of a planetary companion as the primary star evolves during the red giant branch phase. Our aim is to determine a prescription for the capture radius and capture timescales for companions in a wide range of masses, from planet to star, including the effect of a non constant f parameter. Our initial results reproduce past findings. According to our study the mass-loss during the red giant branch does not significantly affect the semi major axis of the planet and for a 1.5 M star the capture radius is approximately 300 R or 1.4 AU.
[34]  [pdf] - 1041144
The formation of high-field magnetic white dwarfs from common envelopes
Comments: Accepted to Proceedings of the National Academy of Sciences. Under PNAS embargo until time of publication
Submitted: 2010-10-07, last modified: 2011-01-12
The origin of highly-magnetized white dwarfs has remained a mystery since their initial discovery. Recent observations indicate that the formation of high-field magnetic white dwarfs is intimately related to strong binary interactions during post-main-sequence phases of stellar evolution. If a low-mass companion, such as a planet, brown dwarf, or low-mass star is engulfed by a post-main-sequence giant, the hydrodynamic drag in the envelope of the giant leads to a reduction of the companion's orbit. Sufficiently low-mass companions in-spiral until they are shredded by the strong gravitational tides near the white dwarf core. Subsequent formation of a super-Eddington accretion disk from the disrupted companion inside a common envelope can dramatically amplify magnetic fields via a dynamo. Here, we show that these disk-generated fields are sufficiently strong to explain the observed range of magnetic field strengths for isolated, high-field magnetic white dwarfs. A higher-mass binary analogue may also contribute to the origin of magnetar fields.
[35]  [pdf] - 955439
Theoretical Support for the Hydrodynamic Mechanism of Pulsar Kicks
Comments: Replaced with Phys. Rev. D accepted version
Submitted: 2010-10-04, last modified: 2010-11-05
The collapse of a massive star's core, followed by a neutrino-driven, asymmetric supernova explosion, can naturally lead to pulsar recoils and neutron star kicks. Here, we present a two-dimensional, radiation-hydrodynamic simulation in which core collapse leads to significant acceleration of a fully-formed, nascent neutron star (NS) via an induced, neutrino-driven explosion. During the explosion, a ~10% anisotropy in the low-mass, high-velocity ejecta lead to recoil of the high-mass neutron star. At the end of our simulation, the NS has achieved a velocity of ~150 km s$^{-1}$ and is accelerating at ~350 km s$^{-2}$, but has yet to reach the ballistic regime. The recoil is due almost entirely to hydrodynamical processes, with anisotropic neutrino emission contributing less than 2% to the overall kick magnitude. Since the observed distribution of neutron star kick velocities peaks at ~300-400 km s$^{-1}$, recoil due to anisotropic core-collapse supernovae provides a natural, non-exotic mechanism with which to obtain neutron star kicks.
[36]  [pdf] - 1025124
Tides and Tidal Engulfment in Post-Main Sequence Binaries: Period Gaps for Planets and Brown Dwarfs Around White Dwarfs
Comments: Replaced with accepted version and correct Fig. 8
Submitted: 2010-02-10, last modified: 2010-08-09
The presence of a close, low-mass companion is thought to play a substantial and perhaps necessary role in shaping post-Asymptotic Giant Branch and Planetary Nebula outflows. During post-main-sequence evolution, radial expansion of the primary star, accompanied by intense winds, can significantly alter the binary orbit via tidal dissipation and mass loss. To investigate this, we couple stellar evolution models (from the zero-age main-sequence through the end of the post-main sequence) to a tidal evolution code. The binary's fate is determined by the initial masses of the primary and the companion, the initial orbit (taken to be circular), and the Reimers mass-loss parameter. For a range of these parameters, we determine whether the orbit expands due to mass loss or decays due to tidal torques. Where a common envelope (CE) phase ensues, we estimate the final orbital separation based on the energy required to unbind the envelope. These calculations predict period gaps for planetary and brown dwarf companions to white dwarfs. The upper end of the gap is the shortest period at which a CE phase is avoided. The lower end is the longest period at which companions survive their CE phase. For binary systems with 1 $M_\odot$ progenitors, we predict no Jupiter-mass companions with periods $\lesssim$270 days. Once engulfed, Jupiter-mass companions do not survive a CE phase. For binary systems consisting of a 1 $M_\odot$ progenitor with a companion 10 times the mass of Jupiter, we predict a period gap between $\sim$0.1 and $\sim$380 days. These results are consistent with both the detection of a $\sim$50 $M_{\rm J}$ brown dwarf in a $\sim$0.003 AU ($\sim$0.08 day) orbit around the white dwarf WD 0137-349 and the tentative detection of a $\sim$2 $M_{\rm J}$ planet in a $\gtrsim$2.7 AU ($\gtrsim$4 year) orbit around the white dwarf GD66.
[37]  [pdf] - 207128
Two Populations of Companions around White Dwarfs: The Effect of Tides and Tidal Engulfment
Comments: To appear in proceedings of "Asymmetric Planetary Nebulae V"; A. A. Zijlstra, I. McDonald and E. Lagadec eds
Submitted: 2010-07-27
During post-main-sequence evolution, radial expansion of the primary star, accompanied by intense winds, can significantly alter the binary orbit via tidal dissipation and mass loss. The fate of a given binary system is determined by the initial masses of the primary and companion, the initial orbit (taken to be circular), the Reimers mass-loss parameter, and the tidal prescription employed. For a range of these parameters, we determine whether the orbit expands due to mass loss or decays due to tidal torques. Where a common envelope (CE) phase ensues, we estimate the final orbital separation based on the energy required to unbind the envelope. These calculations predict period gaps for planetary and brown dwarf companions to white dwarfs. In particular, the lower end of the gap is the longest period at which companions survive their CE phase while the upper end of the gap is the shortest period at which a CE phase is avoided. For binary systems with 1 $M_\odot$ progenitors, we predict no Jupiter-mass companions with periods $\lesssim$270 days. For binary systems consisting of a 1 $M_\odot$ progenitor with a 10 Jupiter-mass companion, we predict a close, post-CE population with periods $\lesssim$0.1 days and a far population with periods $\gtrsim$380 days. These results are consistent with the detection of a $\sim$50 $M_{\rm J}$ brown dwarf in a $\sim$0.08 day orbit around the white dwarf WD 0137-349 and the tentative detection of a $\sim$2 $M_{\rm J}$ planet in a $\sim$4 year orbit around the white dwarf GD66.
[38]  [pdf] - 903232
Dimension as a Key to the Neutrino Mechanism of Core-Collapse Supernova Explosions
Comments: Accepted to ApJ July 7th, Replaced with accepted version
Submitted: 2010-06-18, last modified: 2010-07-07
We explore the dependence on spatial dimension of the viability of the neutrino heating mechanism of core-collapse supernova explosions. We find that the tendency to explode is a monotonically increasing function of dimension, with 3D requiring $\sim$40$-$50\% lower driving neutrino luminosity than 1D and $\sim$15$-$25\% lower driving neutrino luminosity than 2D. Moreover, we find that the delay to explosion for a given neutrino luminosity is always shorter in 3D than 2D, sometimes by many hundreds of milliseconds. The magnitude of this dimensional effect is much larger than the purported magnitude of a variety of other effects, such as nuclear burning, inelastic scattering, or general relativity, which are sometimes invoked to bridge the gap between the current ambiguous and uncertain theoretical situation and the fact of robust supernova explosions. Since real supernovae occur in three dimensions, our finding may be an important step towards unraveling one of the most problematic puzzles in stellar astrophysics. In addition, even though in 3D we do see pre-explosion instabilities and blast asymmetries, unlike the situation in 2D, we do not see an obvious axially-symmetric dipolar shock oscillation. Rather, the free energy available to power instabilites seems to be shared by more and more degrees of freedom as the dimension increases. Hence, the strong dipolar axisymmetry seen in 2D and previously identified as a fundamental characteristic of the shock hydrodynamics may not survive in 3D as a prominent feature.
[39]  [pdf] - 33646
A New Method for Determining the Milky Way Bar Pattern Speed
Comments: 6 pages, 2 figures. To appear in "Tumbling, twisting, and winding galaxies: Pattern speeds along the Hubble sequence", E. M. Corsini and V. P. Debattista (eds.), Memorie della Societa` Astronomica Italiana
Submitted: 2010-02-08
Previous work has related the Galactic bar to structure in the local stellar velocity distribution. Here we show that the bar also influences the spatial gradients of the velocity vector via the Oort constants. By numerical integration of test-particles we simulate measurements of the Oort C value in a gravitational potential including the Galactic bar. We account for the observed trend that C is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C with our simulations we improve on previous models of the bar, estimating that the bar pattern speed is Omega_b/Omega_0=1.87\pm0.04, where Omega_0 is the local circular frequency, and the bar angle lies within 20[deg]<phi_0<45[deg]. We find that the Galactic bar affects measurements of the Oort constants A and B less than ~2 km/s/kpc for the hot stars.
[40]  [pdf] - 21177
Is the Milky Way ringing? The hunt for high velocity streams
Comments: 5 pages, 4 figures, Accepted to MNRAS Letters, added references
Submitted: 2009-02-10, last modified: 2009-03-24
We perform numerical simulations of a stellar galactic disk with initial conditions chosen to represent an unrelaxed population which might have been left following a merger. Stars are unevenly distributed in radial action angle, though the disk is axisymmetric. The velocity distribution in the simulated Solar neighborhood exhibits waves traveling in the direction of positive v, where u,v are the radial and tangential velocity components. As the system relaxes and structure wraps in phase space, the features seen in the uv-plane move closer together. We show that these results can be obtained also by a semi-analytical method. We propose that this model could provide an explanation for the high velocity streams seen in the Solar neighborhood at approximate v in km/s, of -60 (HR 1614), -80 (Arifyanto and Fuchs 2006), -100 (Arcturus), and -160 (Klement et al. 2008). In addition, we predict four new features at v ~ -140, -120, 40 and 60 km/s. By matching the number and positions of the observed streams, we estimate that the Milky Way disk was strongly perturbed ~1.9 Gyr ago. This event could have been associated with Galactic bar formation.
[41]  [pdf] - 13842
Magnetic Mixing in Red Giant and Asymptotic Giant Branch Stars
Comments: Replaced to correct small error in published version: In \S 2.1, paragraphs 2 and 3 incorrectly refer to the poloidal field when qualitatively discussing magnetic diffusion in the shear zone. The correct physical interpretation is that the toroidal field diffuses through the shear zone consistent with the value of $\beta_\phi$
Submitted: 2008-06-24, last modified: 2008-11-25
The available information on isotopic abundances in the atmospheres of low-mass Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars requires that episodes of extensive mixing occur below the convective envelope, reaching down to layers close to the hydrogen burning shell (Cool Bottom Processing). Recently \cite{Busso:2007jw} suggested that dynamo-produced buoyant magnetic flux tubes could provide the necessary physical mechanisms and also supply sufficient transport rates. Here, we present an $\alpha-\Omega$ dynamo in the envelope of an RGB/AGB star in which shear and rotation drain via turbulent dissipation and Poynting flux. In this context, if the dynamo is to sustain throughout either phase, convection must resupply shear. Under this condition, volume-averaged, peak toroidal field strengths of $<B_\phi>\simeq3\times10^3$ G (RGB) and $<B_\phi>\simeq5\times10^3$ G (AGB) are possible at the base of the convection zone. If the magnetic fields are concentrated in flux tubes, the corresponding field strengths are comparable to those required by Cool Bottom Processing.
[42]  [pdf] - 18728
A dense disk of dust around the born-again Sakurai's object
Comments: A&A Letter, accepted
Submitted: 2008-11-20
In 1996, Sakurai's object (V4334 Sgr) suddenly brightened in the centre of a faint Planetary Nebula (PN). This very rare event was interpreted as the reignition of a hot white dwarf that caused a rapid evolution back to the cool giant phase. From 1998 on, a copious amount of dust has formed continuously, screening out the star which has remained embedded in this expanding high optical depth envelope. The new observations, reported here, are used to study the morphology of the circumstellar dust in order to investigate the hypothesis that Sakurai's Object is surrounded by a thick spherical envelope of dust. We have obtained unprecedented, high-angular resolution spectro-interferometric observations, taken with the mid-IR interferometer MIDI/VLTI, which resolve the dust envelope of Sakurai's object. We report the discovery of a unexpectedly compact (30 x 40 milliarcsec, 105 x 140 AU assuming a distance of 3.5 kpc), highly inclined, dust disk. We used Monte Carlo radiative-transfer simulations of a stratified disk to constrain its geometric and physical parameters, although such a model is only a rough approximation of the rapidly evolving dust structure. Even though the fits are not fully satisfactory, some useful and robust constraints can be inferred. The disk inclination is estimated to be 75+/-3 degree with a large scale height of 47+/-7 AU. The dust mass of the disk is estimated to be 6 10^{-5} solar mass. The major axis of the disk (132+/-3 degree) is aligned with an asymmetry seen in the old PN that was re-investigated as part of this study. This implies that the mechanism responsible for shaping the dust envelope surrounding Sakurai's object was already at work when the old PN formed.
[43]  [pdf] - 334
Crystalline silicates and dust processing in the protoplanetary disks of the Taurus young cluster
Comments: 116 pages, including 11 figures and 4 tables
Submitted: 2007-04-11, last modified: 2008-08-20
We characterize the crystalline silicate content and spatial distribution of small dust grains in a large sample of protoplanetary disks in the Taurus-Auriga young cluster, using Spitzer Space Telescope mid-infrared spectra. In turn we use the results to analyze the evolution of structure and composition of these 1-2 Myr-old disks around Solar- and later-type young stars, and test the standard models of dust processing which result in the conversion of originally amorphous dust into minerals. We find strong evidence of evolution of the dust crystalline mass fraction in parallel with that of the structure of the disks, in the sense that increasing crystalline mass fraction is strongly linked to dust settling to the disk midplane. We also confirm that the crystalline silicates are confined to small radii, r < 10 AU. However, we see no significant correlation of crystalline mass fraction with stellar mass or luminosity, stellar accretion rate, disk mass, or disk/star mass ratio, as would be expected in the standard models of dust processing based upon photo-evaporation and condensation close to the central star, accretion-heating-driven annealing at r < 1 AU, or spiral-shock heating at r < 10 AU, with or without effective radial mixing mechanisms. Either another grain-crystallizing mechanism dominates over these, or another process must be at work within the disks to erase the correlations they produce. We propose one of each sort that seem to be worth further investigation, namely X-ray heating and annealing of dust grains, and modulation of disk structure by giant-planetary formation and migration.
[44]  [pdf] - 9188
Towards a Spectral Technique for Determining Material Geometry Around Evolved Stars: Application to HD 179821
Comments: 18 pages, 4 figures, 3 tables; accepted to MNRAS
Submitted: 2008-01-18, last modified: 2008-05-12
HD 179821 is an evolved star of unknown progenitor mass range (either post-Asymptotic Giant Branch or post-Red Supergiant) exhibiting a double peaked spectral energy distribution (SED) with a sharp rise from $\sim8-20$ $\mu$m. Such features have been associated with ejected dust shells or inwardly truncated circumstellar discs. In order to compare SEDs from both systems, we employ a spherically symmetric radiative transfer code and compare it to a radiative, inwardly truncated disc code. As a case study, we model the broad-band SED of HD 179821 using both codes. Shortward of 40 $\mu$m, we find that both models produce equivalent fits to the data. However, longward of 40 $\mu$m, the radial density distribution and corresponding broad range of disc temperatures produce excess emission above our spherically symmetric solutions and the observations. For HD 179821, our best fit consists of a $T_{eff}=7000$ K central source characterized by $\tau_V\sim1.95$ and surrounded by a radiatively driven, spherically symmetric dust shell. The extinction of the central source reddens the broad-band colours so that they resemble a $T_{eff}=5750$ K photosphere. We believe that HD 179821 contains a hotter central star than previously thought. Our results provide an initial step towards a technique to distinguish geometric differences from spectral modeling.
[45]  [pdf] - 436169
The Bipolar Engines of post-AGB stars: Transient Dynamos and Common Envelopes
Comments: 8 pages, 4 figures. To appear in Asymmetrical Planetary Nebulae IV. R. L. M. Corradi, A. Manchado, N. Soker eds - replaced to properly link to NASA ADS
Submitted: 2007-07-25, last modified: 2007-12-14
The physical mechanism(s) responsible for transitioning from a spherical Asymptotic Giant Branch (AGB) star to an asymmetric post-AGB (pAGB) object is poorly understood. In particular, excess momenta in the outflows of pAGB objects suggest that a binary may be required to supply an additional source of energy and angular momentum. The extraction of rotational energy from the engine is likely fundamental and may be facilitated if a dynamo is operating in the interior. In this regard, single star magnetic outflow models have been proposed as mechanisms for producing and shaping PNe, however these models neglect the back-reaction of the large-scale magnetic field on the flow. Here we present a transient $\alpha-\Omega$ dynamo operating in the envelope of an AGB star in (1) an isolated setting and (2) a common envelope in which the secondary is a low-mass companion in-spiraling in the AGB interior. The back reaction of the fields on the shear is included and differential rotation and rotation deplete via turbulent dissipation and Poynting flux. For an isolated star, the shear must be resupplied in order to sufficiently sustain the dynamo. We comment on the energy requirements that convection must satisfy to accomplish this. For the common envelope case, a robust dynamo can result as the companion provides an additional source of energy and angular momentum.
[46]  [pdf] - 8139
Dynamos and Chemical Mixing in Evolved Stars
Comments: 7 pages, 3 figures. To appear in AIP Proceedings of the IXth Torino Workshop on AGB Nucleosynthesis
Submitted: 2007-12-14
In low-mass Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars, anomalous mixing must transport material near the hydrogen-burning shell to the convective envelope. Recently, it was suggested that buoyant magnetic flux tubes could supply the necessary transport rate (Busso et al. 2007). The fields are assumed to originate from a dynamo operating in the stellar interior. Here, we show what is required of an $\alpha-\Omega$ dynamo in the envelope of an AGB star to maintain these fields. Differential rotation and rotation drain via turbulent dissipation and Poynting flux, so if shear can be resupplied by convection, then large-scale toroidal field strengths of $\left<B_\phi\right>\simeq3\times10^4$ G can be sustained at the base of the convection zone.
[47]  [pdf] - 90614
New Constraints on the Galactic Bar
Comments: 4 pages, 2 figures, Accepted to ApJ Letters. Replaced with accepted version
Submitted: 2007-03-30, last modified: 2007-12-11
Previous work has related the Galactic Bar to structure in the local stellar velocity distribution. Here we show that the Bar also influences the spatial gradients of the velocity vector via the Oort constants. By numerical integration of test-particles we simulate measurements of the Oort C value in a gravitational potential including the Galactic Bar. We account for the observed trend that C is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C with our simulations we improve on previous models of the Bar, estimating that the Bar pattern speed is Omega_b/Omega_0=1.87\pm0.04, where Omega_0 is the local circular frequency, and the Bar angle lies within 20<phi_0<45 deg. We find that the Galactic Bar affects measurements of the Oort constants A and B less than ~2 km/s/kpc for the hot stars.
[48]  [pdf] - 4406
Planetary Nebulae Principles & Paradigms: Binaries, Accretion, Magnetic Fields
Comments: 8 pages, 1 figure; for the proceedings of "Asymmetric Planetary Nebuale IV," R. L. M. Corradi, A. Manchado, N. Soker eds
Submitted: 2007-08-30, last modified: 2007-09-30
Observations suggest that many, if not all, post AGB systems evolve through an aspherical outflow phase. Such outflows require a sufficient engine rotational energy which binaries can provide. Via common envelope evolution, binaries can directly eject equatorial outflows or produce poloidal outflows from magnetized accretion disks around the primary or secondary. We discuss how accretion driven magnetohydrodynamic outflow models all make similar predictions for the outflow power and speed and we distinguish between the launch vs. propagation regimes of such outflows. We suggest that the high velocity bipolar outflows observed in planetary nebulae (PNe) and the lower velocity but higher power bipolar outflows observed in pre-PNe (pPNe) are kinematically consistent with time dependent accretion onto a white dwarf (WD) within a depleting envelope. Since the WD primary core is always present in all post-AGB systems, accretion onto this core is potentially common. Previous work has focused on core accretion from sub-stellar companions, but low mass stellar companions may be more important, and further work is needed.
[49]  [pdf] - 4958
The Formation of Crystalline Dust in AGB Winds from Binary Induced Spiral Shocks
Comments: Submitted to ApJL
Submitted: 2007-09-14
As stars evolve along the Asymptotic Giant Branch, strong winds are driven from the outer envelope. These winds form a shell, which may ultimately become a planetary nebula. Many planetary nebulae are highly asymmetric, hinting at the presence of a binary companion. Some post-Asymptotic Giant Branch objects are surrounded by torii of crystalline dust, but there is no generally accepted mechanism for annealing the amorphous grains in the wind to crystals. In this Letter, we show that the shaping of the wind by a binary companion is likely to lead to the formation of crystalline dust in the orbital plane of the binary.
[50]  [pdf] - 85318
Isolated vs. Common Envelope Dynamos in Planetary Nebula Progenitors
Comments: 25 pages, 7 figures, submitted to MNRAS; revised to fix Fig. 5
Submitted: 2006-09-26, last modified: 2006-10-02
The origin, evolution and role of magnetic fields in the production and shaping of proto-planetary and planetary nebulae (PPNe, PNe) is a subject of active research. Most PNe and PPNe are axisymmetric with many exhibiting highly collimated outflows, however, it is important to understand whether such structures can be generated by isolated stars or require the presence of a binary companion. Toward this end we study a dynamical, large-scale alpha-Omega interface dynamo operating in a 3.0 M_sun Asymptotic Giant Branch star (AGB) in both an isolated setting and one in which a low-mass companion is embedded inside the envelope. The back reaction of the fields on the shear is included and differential rotation and rotation deplete via turbulent dissipation and Poynting flux. For the isolated star, the shear must be resupplied in order to sufficiently sustain the dynamo. Furthermore, we investigate the energy requirements that convection must satisfy to accomplish this by analogy to the sun. For the common envelope case, a robust dynamo results, unbinding the envelope under a broad range of conditions. Two qualitatively different types of explosion may arise: (i) magnetically dominated, possibly resulting in collimated bipolar outflows and (ii) thermally induced from turbulent dissipation, possibly resulting in quasi-spherical outflows. A range of models is presented for a variety of companion masses.
[51]  [pdf] - 81514
Low-Mass Binary Induced Outflows from Asymptotic Giant Branch Stars
Comments: 24 pages, 6 figures, submitted to MNRAS
Submitted: 2006-04-20
A significant fraction of planetary nebulae (PNe) and proto-planetary nebulae (PPNe) exhibit aspherical, axisymmetric structures, many of which are highly collimated. The origin of these structures is not entirely understood, however recent evidence suggests that many observed PNe harbor binary systems, which may play a role in their shaping. In an effort to understand how binaries may produce such asymmetries, we study the effect of low-mass (< 0.3 M_sun) companions (planets, brown dwarfs and low-mass main sequence stars) embedded into the envelope of a 3.0 M_sun star during three epochs of its evolution (Red Giant Branch, Asymptotic Giant Branch (AGB), interpulse AGB). We find that common envelope evolution can lead to three qualitatively different consequences: (i) direct ejection of envelope material resulting in a predominately equatorial outflow, (ii) spin-up of the envelope resulting in the possibility of powering an explosive dynamo driven jet and (iii) tidal shredding of the companion into a disc which facilitates a disc driven jet. We study how these features depend on the secondary's mass and discuss observational consequences.
[52]  [pdf] - 1233432
Extracting Rotational Energy in Supernova Progenitors: Transient Poynting Flux Growth vs. Turbulent Dissipation
Comments: 26 pages, submitted to New Astronomy (revised to incorporate the depletion of shear energy via both turbulent diffusion and magnetic stress)
Submitted: 2004-10-28, last modified: 2005-09-03
Observational evidence for anisotropy in supernovae (SN) may signal the importance of angular momentum and differential rotation in the progenitors. Free energy in differential rotation and rotation can be extracted magnetically or via turbulent dissipation. The importance that magnetohydrodyamic jets and coronae may play in driving SN motivates understanding large scale dynamos in SN progenitors. We develop a dynamical large scale interface dynamo model in which the differential rotation and rotation deplete both through Poynting flux and turbulent diffusion. We apply the model to a differentially rotating core surrounded by a convection zone of a SN progenitor from a initial 15$M_\odot$ star. Unlike the Sun, the dynamo is transient because the differential rotation is primarily due to the initial collapse. Up to $\sim 10^{51}$erg can be drained into time-integrated Poynting flux and heat, the relative fraction of which depends on the relative amount of turbulence in the shear layer vs. convection zone and the fraction of the shear layer into which the magnetic field penetrates. Both sinks can help facilitate explosions and could lead to different levels of anisotropy and pulsar kicks. In all cases, the poloidal magnetic field is much weaker than the toroidal field, and the Poynting flux is lower than previous estimates which invoke the magnitude of the total magnetic energy. A signature of a large scale dynamo is that the oscillation of the associated Poynting flux on $\sim 1$ sec time scales, implying the same for the energy delivery to a SN.
[53]  [pdf] - 48844
Magnetorotational Instability in Liquid Metal Couette Flow
Comments: 12 pages, 22 figures, 1 table. To appear in the Astrophysical Journal
Submitted: 2002-04-17
Despite the importance of the magnetorotational instability (MRI) as a fundamental mechanism for angular momentum transport in magnetized accretion disks, it has yet to be demonstrated in the laboratory. A liquid sodium alpha-omega dynamo experiment at the New Mexico Institute of Mining and Technology provides an ideal environment to study the MRI in a rotating metal annulus (Couette flow). A local stability analysis is performed as a function of shear, magnetic field strength, magnetic Reynolds number, and turbulent Prandtl number. The later takes into account the minimum turbulence induced by the formation of an Ekman layer against the rigidly rotating end walls of a cylindrical vessel. Stability conditions are presented and unstable conditions for the sodium experiment are compared with another proposed MRI experiment with liquid gallium. Due to the relatively large magnetic Reynolds number achievable in the sodium experiment, it should be possible to observe the excitation of the MRI for a wide range of wavenumbers and further to observe the transition to the turbulent state.