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

Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primaries. But CE evolution with AGB progenitors is desirable to simulate because such events are the likely progenitors of most bi-polar planetary nebulae (PNe), and prominent observational testing grounds for CE physics. Here we present a high resolution global simulation of CE evolution involving an AGB primary and $1\,\mathrm{M}_\odot$ secondary, evolved for $20$ orbital revolutions. During the last $16$ of these orbits, the envelope unbinds at an almost constant rate of about $0.1$-$0.2\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$. If this rate were maintained, the envelope would be unbound in less than $10\,\mathrm{yr}$. The dominant source of this unbinding is consistent with inspiral; we assess the influence of the ambient medium to be subdominant. We compare this run with a previous run that used an RGB phase primary evolved from the same $2\,\mathrm{M}_\odot$ main sequence star to assess the influence of the evolutionary state of the primary. When scaled appropriately, the two runs are quite similar, but with some important differences.

We model the settlement of the galaxy by space-faring civilizations in order to address issues related to the Fermi Paradox. We explore the problem in a way that avoids assumptions about the intent and motivation of any exo-civilization seeking to settle other planetary systems. We first consider the speed of an advancing settlement via probes of finite velocity and range to determine if the galaxy can become inhabited with space-faring civilizations on timescales shorter than its age. We also include the effect of stellar motions on the long term behavior of the settlement front which adds a diffusive component to its advance. The results of these models demonstrate that the Milky Way can be readily 'filled-in' with settled stellar systems under conservative assumptions about interstellar spacecraft velocities and launch rates. We then consider the question of the galactic steady-state achieved in terms of the fraction of settled planets. We do this by considering the effect of finite settlement civilization lifetimes on the steady states. We find a range of parameters for which the galaxy supports a population of interstellar space-faring civilizations even though some settleable systems are uninhabited. Both results point to ways in which Earth might remain unvisited in the midst of an inhabited galaxy. Finally we consider how our results can be combined with the finite horizon for evidence of previous settlements in Earth's geologic record. Our steady-state model can constrain the probabilities for an Earth visit by a settling civilization before a given time horizon. These results break the link between Hart's famous "Fact A" (no interstellar visitors on Earth now) and the conclusion that humans must, therefore, be the only technological civilization in the galaxy.

The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outward from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-$\alpha$ flux required to drive the photo-evaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-$\alpha$ fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.

The influences and consequences of toroidal magnetic fields in shaping the visible lobes of pre planetary nebulae ("prePNe") are explored in this, the last of a series of papers of parameter studies of prePN evolution. To probe these influences we start with the steady, diverging, and field-free wind model of our previous papers and add weak to moderate toroidal fields to the winds in order to generate arrays of outcomes after 500 y, after which the structures grow almost homologously. As expected, toroidal fields in the stellar winds invariably form very thin and dense axial features whose structure is best described as a thin cold jet with an ultra-dense and neutral leading knot, or "bullet", at its tip. The speed of the leading knot depends only on the density contrast (the ratio of injected to ambient gas densities at the nozzle) and wind injection speed, but not on the field strength or opening angle. The lobes formed by the ram pressure of the winds take a variety of forms and sizes that depend primarily on the geometric structure of the injected gas and the density contrast. About 20% of the HST images of prePNe show the unique signatures of shaping by toroidal fields. Pairs of low-ionization knots seen along the major axis of fully ionized PNe, often called "FLIERs" are easily explained as the very dense, cold , and neutral remnants of magnetically formed knots.

We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium.

We compute the forces, torque and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably approximated by conventional analytic theory and peaks at values proportional to the companion mass. Good agreement between global and local 3D "wind tunnel" simulations, including similar net drag force and flow pattern, is obtained for comparable regions of parameter space. However, subsequent to the first periastron passage, the drag force is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and depends only weakly on companion mass. The discrepancy is exacerbated for larger companion mass and when the inter-particle separation reduces to the Bondi-Hoyle-Lyttleton accretion radius, creating a turbulent thermalized region. Greater flow symmetry during this phase leads to near balance of opposing gravitational forces in front of and behind the companion, hence a small net drag. The reduced drag force at late times helps explain why companion-core separations necessary for envelope ejection are not reached by the end of limited duration CE simulations.

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.

We report on simulations of laboratory experiments in which magnetized supersonic flows are driven through a wire mesh. The goal of the study was to investigate the ability of such a configuration to generate supersonic, MHD turbulence. We first report on the morphological structures that develop in both magnetized and non-magnetized cases. We then analyze the flow using a variety of statistical measures, including power spectra and probability distribution functions of the density. Using these results we estimate the sonic mach number in the flows downstream of the wire mesh. We find the initially hypersonic (M=20) planar shock through the wire mesh does lead to downstream turbulent conditions. However, in both magnetized and non-magnetized cases, the resultant turbulence was marginally supersonic to transonic (M~1), and highly anisotropic in structure.

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

We have used the Hubble Space Telescope to acquire new broad-band and narrow-band images of the optical line emission and red continuum associated with the HH 7-11 stellar jet in the NGC 1333 star formation region. Combining the new narrow-band images of H$\alpha$, [O~I] $\lambda$6300 and [S II] $\lambda$6716 allows us to measure electron densities and excitations at each point in the outflow with the spatial resolution of HST, while the I-band image traces out the boundary of the cavity evacuated by the outflow. Comparing these images with those taken $\sim$ 20 years ago yields high precision proper motions for all the HH objects in the outflow. HH 11 is a bullet-like clump, and emerges from the exciting source SVS 13A towards the Earth at 24 degrees to line of sight. In contrast, HH 8 and HH 10 consist of two rings of shocked gas that show no bulk proper motions even though the emitting gas is blueshifted. The HH 8 rings are expanding with time. These shocks mark places where ambient material located along the path of the jet redirects the outflow. HH 7 consists of multiple shells, and emits strongly in H$_2$ in what appears to be a terminal bow shock for the outflow, implying that the jet has yet to fully break out of its nascent cloud core. The jet largely fragments into clumps by the time it reaches HH 7. As in the case of HH 110, deflection from ambient material plays a key role in producing observable shock waves in the HH 7-11 outflow.

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.

We develop a physical framework for interpreting high-resolution images of pre planetary nebule ("prePNe") with pairs of candle shaped lobes. We use hydrodynamical models to infer the historical properties of the flows injected from the nucleus that shape the lobes into standard forms. First, we find a suitable set of parameters of a fast, collimated, tapered flow that is actively reshaped by an exterior slow AGB wind and that nicely fits the basic shape, kinematics, mass, and momenta of this class of prePNe. Next we vary the most influential parameters of this "baseline" model-such as density, speed, and geometry-to see how changes in the flow parameters affect the nebular observables after 900y. Several generic conlusions emerge, such as the injected flows that create the hollow candle-shaped lobes must be light, "tapered", and injected considerably faster than the lobe expansion speed. Multi-polar and starfish prePNe probably evolve from wide angle flows in which thin-shell instabilites corrugate their leading edges. We show how the common linear relationship of Doppler shift and position along the lobe is a robust outcome the interaction of tapered diverging streamlines with the lobes' curved walls. Finally we probe how magnetic fields affect the basline model by adding a toroidal field to the injected baseline flow. Examples of prePNe and PNe that may have been magnetically shaped are listed. We conclude that the light, field-free, tapered baseline flow model is an successful and universal pardigm for unravelling the histories of lobe formation in prePNe.

The photoionization-driven evaporation of planetary atmospheres has emerged as a potentially fundamental process for planets on short period orbits. While 1-D studies have proven the effectiveness of stellar fluxes at altering the atmospheric mass and composition for sub-Jupiter mass planets, there remains much that is uncertain with regard to the larger-scale, multidimensional nature of such "planetary wind" flows. In this paper we use a new radiation-hydrodynamic platform to simulate atmospheric evaporative flows. Using the AstroBEAR AMR multiphysics code in a co-rotating frame centered on the planet, we model the transfer of ionizing photons into the atmosphere, the subsequent launch of the wind and the wind's large scale evolution subject to tidal and non-inertial forces. We run simulations for planets of 0.263 and 0.07 Jupiter masses and stellar fluxes of $2 \times 10^{13}$ and $2 \times 10^{14}$ photons/cm^2/s. Our results reveal new, potentially observable planetary wind flow patterns, including the development, in some cases, of an extended neutral tail lagging behind the planet in its orbit.

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.

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.

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.

Observations of transiting extra-solar planets provide rich sources of data for probing the in-system environment. In the WASP-12 system, a broad depression in the usually-bright MgII h&k lines has been observed, in addition to atmospheric escape from the extremely hot Jupiter WASP-12b. It has been hypothesized that a translucent circumstellar cloud is formed by the outflow from the planet, causing the observed signatures. We perform 3D hydrodynamic simulations of the full system environment of WASP-12, injecting a planetary wind and stellar wind from their respective surfaces. We find that a torus of density high enough to account for the lack of MgII h&k line core emission in WASP-12 can be formed in approximately 13 years. We also perform synthetic observations of the Lyman-alpha spectrum at different points in the planet's orbit, which demonstrate that significant absorption occurs at all points in the orbit, not just during transits, as suggested by the observations.

If an industrial civilization had existed on Earth many millions of years prior to our own era, what traces would it have left and would they be detectable today? We summarize the likely geological fingerprint of the Anthropocene, and demonstrate that while clear, it will not differ greatly in many respects from other known events in the geological record. We then propose tests that could plausibly distinguish an industrial cause from an otherwise naturally occurring climate event.

M2-9, or the "Butterfly Nebula" is one of the most iconic outflow sources from an evolved star. In this paper we present a hydrodynamic model of M2-9 in which the nebula is formed and shaped by a steady, low-density ("light"), mildly collimated "spray" of gas injected at 200 km s^-1 that interacts with a far denser, intrinsically simple pre-existing AGB wind has slowly formed all of the complex features within M2-9's lobes (including the knot pairs N3/S3 and N4/S4 at their respective leading edges, and the radial gradient of Doppler shifts within 20" of the nucleus). We emphasize that the knot pairs are not ejected from the star but formed in situ. In addition, the observed radial speed of the knots is only indirectly related to the speed of the gas injected by the star. The model allows us to probe the early history of the wind geometry and lobe formation. We also formulate a new estimate of the nebular distance D = 1.3 kpc. The physical mechanism that accounts for the linear radial speed gradient in M2-9 applies generally to many other pre planetary nebulae whose hollow lobes exhibit similar gradients along their edges.

Human influence on the biosphere has been evident at least since the development of widespread agriculture, and some stratigraphers have suggested that the activities of modern civilization indicate a geological epoch transition. The study of the anthropocene as a geological epoch, and its implication for the future of energy-intensive civilizations, is an emerging transdisciplinary field in which astrobiology can play a leading role. Habitability research of Earth, Mars, and exoplanets examines extreme cases relevant for understanding climate change as a planetary process. Energy-intensive civilizations will also face thermodynamic limits to growth, which provides an important constraint for estimating the longevity of human civilization and guiding the search for extraterrestrial intelligence. We recommend that missions concepts such as LUVOIR, HabEx, and OST be pursued in order to make significant progress toward understanding the future evolution of life on our planet and the possible evolution of technological, energy-intensive life elsewhere in the universe.

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.

We develop a classification scheme for the evolutionary state of planets based on the non-equilibrium thermodynamics of their coupled systems, including the presence of a biosphere and the possibility of what we call an agency-dominated biosphere (i.e. an energy-intensive technological species). The premise is that Earths entry into the Anthropocene represents what might be from an astrobiological perspective a predictable planetary transition. We explore this problem from the perspective of the solar system and exoplanet studies. Our classification discriminates planets by the forms of free energy generation driven from stellar forcing. We then explore how timescales for global evolutionary processes on Earth might be synchronized with ecological transformations driven by increases in energy harvesting and its consequences (which might have reached a turning point with global urbanization). Finally, we describe quantitatively the classification scheme based on the maintenance of chemical disequilibrium in the past and current Earth systems and on other worlds in the solar system. In this perspective, the beginning of the Anthropocene can be seen as the onset of the hybridization of the planet - a transitional stage from one class of planetary systems interaction to another. For Earth, this stage occurs as the effects of human civilization yield not just new evolutionary pressures, but new selected directions for novel planetary ecosystem functions and their capacity to generate disequilibrium and enhance planetary dissipation.

We present numerical simulations of reorienting oblique shocks that form in the collision layer between magnetized colliding flows. Reorientation aligns parsec-scale post-shock filaments normal to the background magnetic field. We find that reorientation begins with pressure gradients between the collision region and the ambient medium. This drives a lateral expansion of post-shock gas, which reorients the growing filament from the outside-in (i.e. from the flow-ambient boundary, toward the colliding flows axis). The final structures of our simulations resemble polarization observations of filaments in Taurus and Serpens South, as well as the integral-shaped filament in Orion A. Given the ubiquity of colliding flows in the interstellar medium, shock reorientation may be relevant to the formation of filaments normal to magnetic fields.

We present a detailed hydrodynamic model that matches the present structure of the well-observed preplanetary nebula OH231.8+04.2. The purpose of the model is to present a physically justified and coherent picture of its evolutionary history from about 100 years of the start of the formation of its complex outer structures to the present. We have adopted a set of initial conditions that are heavily constrained by high-quality observations of its present structure and kinematics. The shaping of the nebula occurs while the densities of the flows are "light": less than the surrounding AGB-wind environment. The simulations show that pairs of essentially coeval clumps and sprays of the same extent and density but different outflow speeds sculpted both the pair of thin axial flow "or spine" and the bulbs. The total ejected mass and momentum in the best fit model are surprisingly large -- 3 M_sun and 2.2 x 10^41 gm cm s-1, respectively -- however, these values are reduced by up to a factor of ten in other models that fit the data almost as well. Our ultimate goal is to combine the present model results of masses, momenta, flow speeds, and flow geometries for OH231 with those of other models to be published in the future in order to find common attributes of their ejection histories.

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.

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.

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.

The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this shear creates inhospitable environments for star formation. As the degree of shear increases (i.e. the obliquity of the shock increases), we find that it takes longer for sink particles to form, they form in lower numbers, and they tend to be less massive. With regard to magnetic fields, we find that even a weak field stalls gravitational collapse within forming clouds. Additionally, an initially oblique collision interface tends to reorient over time in the presence of a magnetic field, so that it becomes normal to the oncoming flows. This was demonstrated by our most oblique shock interface, which became fully normal by the end of the simulation.

Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-{\alpha} absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption.

The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and second, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here.

Supersonic outflows from objects as varied as stellar jets, massive stars and novae often exhibit multiple shock waves that overlap one another. When the intersection angle between two shock waves exceeds a critical value, the system reconfigures its geometry to create a normal shock known as a Mach stem where the shocks meet. Mach stems are important for interpreting emission-line images of shocked gas because a normal shock produces higher postshock temperatures and therefore a higher-excitation spectrum than an oblique one does. In this paper we summarize the results of a series of numerical simulations and laboratory experiments designed to quantify how Mach stems behave in supersonic plasmas that are the norm in astrophysical flows. The experiments test analytical predictions for critical angles where Mach stems should form, and quantify how Mach stems grow and decay as intersection angles between the incident shock and a surface change. While small Mach stems are destroyed by surface irregularities and subcritical angles, larger ones persist in these situations, and can regrow if the intersection angle changes to become more favorable. The experimental and numerical results show that although Mach stems occur only over a limited range of intersection angles and size scales, within these ranges they are relatively robust, and hence are a viable explanation for variable bright knots observed in HST images at the intersections of some bow shocks in stellar jets.

High resolution observations of Young Stellar Object (YSO) jets show them to be composed of many small-scale knots or clumps. In this paper we report results of 3-D numerical simulations designed to study how such clumps interact and create morphologies and kinematic patterns seen in emission line observations. Our simulations focus on clump scale dynamics by imposing velocity differences between spherical, over-dense regions which then lead to the formation of bow shocks as faster clumps overtake slower material. We show that much of the spatial structure apparent in emission line images of jets arises from the dynamics and interactions of these bow shocks. Our simulations show a variety of time-dependent features, including bright knots associated with Mach stems where the shocks intersect, a "frothy" emission structure that arises from the presence of the Non-linear Thin Shell Instability (NTSI) along the surfaces of the bow shocks, and the merging and fragmentation of clumps. Our simulations use a new non-equilibrium cooling method to produce synthetic emission maps in H$\alpha$ and [S II]. These are directly compared with multi-epoch Hubble Space Telescope (HST) observations of Herbig-Haro (HH) jets. We find excellent agreement between features seen in the simulations and the observations in terms of both proper motion and morphologies. Thus we conclude that YSO jets may be dominated by heterogeneous structures and that interactions between these structures and the shocks they produce can account for many details of YSO jet evolution.

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.

In this paper we address the cosmic frequency of technological species. Recent advances in exoplanet studies provide strong constraints on all astrophysical terms in the Drake Equation. Using these and modifying the form and intent of the Drake equation we show that we can set a firm lower bound on the probability that one or more additional technological species have evolved anywhere and at any time in the history of the observable Universe. We find that as long as the probability that a habitable zone planet develops a technological species is larger than ~$10^{-24}$, then humanity is not the only time technological intelligence has evolved. This constraint has important scientific and philosophical consequences.

We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.

Hubble Space Telescope observations show bright knots of H$\alpha$ emission within outflowing young stellar jets. Velocity variations in the flow create secondary bow shocks that may intersect and lead to enhanced emission. When the bow shocks intersect at or above a certain critical angle, a planar shock called a Mach stem is formed. These shocks could produce brighter H$\alpha$ emission since the incoming flow to the Mach stem is parallel to the shock normal. In this paper we report first results of a study using 2-D numerical simulations designed to explore Mach stem formation at the intersection of bow shocks formed by hypersonic "bullets" or "clumps". Our 2-D simulations show how the bow shock shapes and intersection angles change as the adiabatic index $\gamma$ changes. We show that the formation or lack of a Mach stem in our simulations is consistent with the steady-state Mach stem formation theory. Our ultimate goal, which is part of an ongoing research effort, is to characterize the physical and observational consequences of bow shock intersections including the formation of Mach stems.

In this paper we explore the dynamics of radiative axisymmetric MHD jets at high resolution using AMR methods. The goal of the study is to determine both the dynamics and emission properties of such jets. To that end we have implemented microphysics enabling us to produce synthetic maps of H$\alpha$ and [S II]. The jets are pulsed either sinusoidally or randomly via a time-dependent ejection velocity which leads to a complicated structure of internal shocks and rarefactions as has been seen in previous simulations. The high resolution of our simulations allows us to explore in great detail the effect of pinch forces (due to the jet's toroidal magnetic field) within the "working surfaces" where pulses interact. We map the strong H$\alpha$ emission marking shock fronts and the strong [S II] emission inside cooling regions behind shocks as observed with high-resolution images of jets. We find that pinch forces in the stronger field cases produce additional emission regions along the axis as compared with purely hydrodynamic runs. These simulations are a first step to understanding the full 3-D emission properties of radiative MHD jets.

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.

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.

We consider the interaction between a marginally stable Bonnor-Ebert (BE) sphere and the surrounding ambient medium. In particular, we explore how the infall from an evolving ambient medium can trigger collapse of the sphere using 3-D adaptive mesh refinement simulations. We find the resulting collapse dynamics to vary considerably with ambient density. In the highest ambient density cases, infalling material drives a strong compression wave into the cloud. It is the propagation of this wave through the cloud interior that triggers the subsequent collapse. For lower ambient densities, we find the main trigger of collapse to be a quasistatic adjustment of the BE sphere to gravitational settling of the ambient gas. In all cases, we find that the classic "outside-in" collapse mode for super-critical BE spheres is recovered before a protostar (i.e. sink particle) forms. Our work supports scenarios in which BE dynamics naturally begins with either a compression-wave or infall dominated phase, and only later assumes the usual outside-in collapse behavior.

Star formation can be triggered by compression from wind or supernova driven shock waves that sweep over molecular clouds. Because these shocks will likely contain processed elements, triggered star formation has been proposed as an explanation for short lived radioactive isotopes (SLRI) in the Solar System. Previous studies have tracked the triggering event to the earliest phases of collapse and have focused on the shock properties required for both successful star formation and mixing of SLRI's. In this paper, we use Adaptive Mesh Refinement (AMR) simulation methods, including sink particles, to simulate the full collapse and subsequent evolution of a stable Bonnor-Ebert sphere subjected to a shock and post-shock wind. We track the flow of the cloud material after a star (a sink particle) has formed. For non-rotating clouds we find robust triggered collapse and little bound circumstellar material remaining around the post-shock collapsed core. When we add initial cloud rotation we observe the formation of disks around the collapsed core which then interact with the post-shock flow. Our results indicate that these circumstellar disks are massive enough to form planets and are long-lived, in spite of the ablation driven by post-shock flow ram pressure. As a function of the initial conditions, we also track the time evolution of the accretion rates and particle mixing between between the ambient wind and cloud material. The latter is maximized for cases of highest mach number.

In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observations and their key points of contact with theory. Our review of jet/outflow phenomena is organized into three length-scale domains: Source and Disk Scales ($0.1-10^2$ au) where the connection with protostellar and disk evolution theories is paramount; Envelope Scales ($10^2-10^5$ au) where the chemistry and propagation shed further light on the jet launching process, its variability and its impact on the infalling envelope; Parent Cloud Scales ($10^5-10^6$ au) where global momentum injection into cluster/cloud environments become relevant. Issues of feedback are of particular importance on the smallest scales where planet formation regions in a disk may be impacted by the presence of disk winds, irradiation by jet shocks or shielding by the winds. Feedback on envelope scales may determine the final stellar mass (core-to-star efficiency) and envelope dissipation. Feedback also plays an important role on the larger scales with outflows contributing to turbulent support within clusters including alteration of cluster star formation efficiencies (feedback on larger scales currently appears unlikely). A particularly novel dimension of our review is that we consider results on jet dynamics from the emerging field of High Energy Density Laboratory Astrophysics (HEDLA). HEDLA is now providing direct insights into the 3-D dynamics of fully magnetized, hypersonic, radiative outflows.

We explore how questions related to developing a sustainable human civilization can be cast in terms of astrobiology. In particular we show how ongoing astrobiological studies of the coupled relationship between life, planets and their co-evolution can inform new perspectives and direct new studies in sustainability science. Using the Drake Equation as a vehicle to explore the gamut of astrobiology, we focus on its most import factor for sustainability: the mean lifetime <L> of an ensemble of Species with Energy-Intensive Technology (SWEIT). We then cast the problem into the language of dynamical system theory and introduce the concept of a trajectory bundle for SWEIT evolution and discuss how astrobiological results usefully inform the creation of dynamical equations, their constraints and initial conditions. Three specific examples of how astrobiological considerations can be folded into discussions of sustainability are discussed: (1) concepts of planetary habitability, (2) mass extinctions and their possible relation to the current, so-called Anthropocene epoch, and (3) today's changes in atmospheric chemisty (and the climate change it entails) in the context of pervious epochs of biosphere-driven atmospheric and climate alteration (i.e. the Great Oxidation Event).

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.

We study the interaction of strong shock waves with magnetized clumps. Previous numerical work focused on the simplified scenario in which shocked clumps are immersed in a globally uniform magnetic field that extends through both the clump and the ambient medium. Here we consider the complementary circumstance in which the field is completely self-contained within the clumps. This could arise naturally during clump formation via dynamical or thermal instabilities for example as magnetic field pinches off from the ambient medium. Using our AMR MHD code AstroBEAR, we carry out a series of simulations with magnetized clumps that have different self-contained magnetic field configurations. We find that the clump and magnetic evolution are sensitive to the fraction of magnetic field aligned with versus perpendicular to the shock normal. The relative strength of magnetic pressure and tension in the different field configurations allows us to analytically understand the different cases of post-shock evolution. We also show how turbulence and the mixing it implies depends of the initial field configuration and suggest ways in which observed shock-clump morphology may be used as a proxy for identifying internal field topologies a posteriori.

Our ultimate goal is to probe the nature of the collimator of the outflows in the pre PN CRL618. CRL618 is uniquely suited for this purpose owing to its multiple, bright, and carefully studied finger-shaped outflows east and west of its nucleus. We compare new HST images to images in the same filters observed as much as 11 y previously to uncover large proper motions and surface brightness changes in its multiple finger-shaped outflows. The expansion age of the ensemble of fingers is close to 100y. We find strong brightness variations at the fingertips during the past decade. Deep IR images reveal a multiple ring- like structure of the surrounding medium into which the outflows propagate and interact. Tightly constrained three-dimensional ("3D") hydrodynamic models link the properties of the fingers to their possible formation histories. We incorporate previously published complementary information to discern whether each of the fingers of CRL618 are the results of steady, collimated outflows or a brief ejection event that launched a set of bullets about a century ago. Finally, we argue on various physical grounds that fingers of CRL618 are likely to be the result of a spray of clumps ejected at the nucleus of CRL618 since any mechanism that form a sustained set of unaligned jets is unprecedented.

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.

Observational evidence from local star-forming regions mandates that star formation occurs shortly after, or even during, molecular cloud formation. Models of the formation of molecular clouds in large-scale converging flows have identified the physical mechanisms driving the necessary rapid fragmentation. They also point to global gravitational collapse driving supersonic turbulence in molecular clouds. Previous cloud formation models have focused on turbulence generation, gravitational collapse, magnetic fields, and feedback. Here, we explore the effect of structure in the flow on the resulting clouds and the ensuing gravitational collapse. We compare two extreme cases, one with a collision between two smooth streams, and one with streams containing small clumps. We find that structured converging flows lead to a delay of local gravitational collapse ("star formation"). Thus, more gas has time to accumulate, eventually leading to a strong global collapse, and thus to a high star formation rate. Uniform converging flows fragment hydrodynamically early on, leading to the rapid onset of local gravitational collapse and an overall low sink formation rate.

Current Adaptive Mesh Refinement (AMR) simulations require algorithms that are highly parallelized and manage memory efficiently. As compute engines grow larger, AMR simulations will require algorithms that achieve new levels of efficient parallelization and memory management. We have attempted to employ new techniques to achieve both of these goals. Patch or grid based AMR often employs ghost cells to decouple the hyperbolic advances of each grid on a given refinement level. This decoupling allows each grid to be advanced independently. In AstroBEAR we utilize this independence by threading the grid advances on each level with preference going to the finer level grids. This allows for global load balancing instead of level by level load balancing and allows for greater parallelization across both physical space and AMR level. Threading of level advances can also improve performance by interleaving communication with computation, especially in deep simulations with many levels of refinement. While we see improvements of up to 30% on deep simulations run on a few cores, the speedup is typically more modest (5-20%) for larger scale simulations. To improve memory management we have employed a distributed tree algorithm that requires processors to only store and communicate local sections of the AMR tree structure with neighboring processors. Using this distributed approach we are able to get reasonable scaling efficiency (> 80%) out to 12288 cores and up to 8 levels of AMR - independent of the use of threading.

Jets from active galaxies propagate from the central black hole out to the radio lobes on scales of hundreds of kiloparsecs. The jets may encounter giant stars with strong stellar winds and produce observable signatures. For strong winds and weak jets, the interaction may truncate the jet flow during its transit via the mass loading. For weaker jets, the interaction can produce knots in the jet. We present recent 3DMHD numerical simulations to model the evolution of this jet-wind interaction and its observational consequences. We explore (i) the relative mechanical luminosity of the radio jets and the stellar winds (ii) the impact parameter between the jets' axis and the stellar orbital path (iii) the relative magnetic field strength of the jets and the stellar winds.

Modern theoretical models of astrophysical jets combine accretion, rotation, and magnetic fields to launch and collimate supersonic flows from a central source. Near the source, magnetic field strengths must be large enough to collimate the jet requiring that the Poynting flux exceeds the kinetic-energy flux. The extent to which the Poynting flux dominates kinetic energy flux at large distances from the engine distinguishes two classes of models. In magneto-centrifugal launch (MCL) models, magnetic fields dominate only at scales $\lesssim 100$ engine radii, after which the jets become hydrodynamically dominated (HD). By contrast, in Poynting flux dominated (PFD) magnetic tower models, the field dominates even out to much larger scales. To compare the large distance propagation differences of these two paradigms, we perform 3-D ideal MHD AMR simulations of both HD and PFD stellar jets formed via the same energy flux. We also compare how thermal energy losses and rotation of the jet base affects the stability in these jets. For the conditions described, we show that PFD and HD exhibit observationally distinguishable features: PFD jets are lighter, slower, and less stable than HD jets. Unlike HD jets, PFD jets develop current-driven instabilities that are exacerbated as cooling and rotation increase, resulting in jets that are clumpier than those in the HD limit. Our PFD jet simulations also resemble the magnetic towers that have been recently created in laboratory astrophysical jet experiments.

Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both tower flows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars.

The majority of Proto-planetary nebulae (PPN) are observed to have bipolar morphologies. The majority of mature PN are observed to have elliptical shapes. In this paper we address the evolution of PPN/PN morphologies attempting to understand if a transition from strongly bipolar to elliptical shape can be driven by changes in the parameters of the mass loss process. To this end we present 2.5D hydrodynamical simulations of mass loss at the end stages of stellar evolution for intermediate mass stars. We track changes in wind velocity, mass loss rate and mass loss geometry. In particular we focus on the transition from mass loss dominated by a short duration jet flow (driven during the PPN phase) to mass loss driven by a spherical fast wind (produced by the central star of the PN). We address how such changes in outflow characteristics can change the nebula from a bipolar to an elliptical morphology. Our results show that including a period of jet formation in the temporal sequence of PPN to PN produces realistic nebular synthetic emission geometries. More importantly such a sequence provides insight, in principle, into the apparent difference in morphology statistics characterizing PPN and PN systems. In particular we find that while jet driven PPN can be expected to be dominated by bipolar morphologies, systems that begin with a jet but are followed by a spherical fast wind will evolve into elliptical nebulae. Furthermore, we find that spherical nebulae are highly unlikely to ever derive from either bipolar PPN or elliptical PN.

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.

We present new results of 3-D AMR MHD simulations focusing on two distinct aspects of PPN evolution. We first report new simulations of collimated outflows driven entirely by magnetic fields. These Poynting flux dominated "magnetic towers" hold promise for explaining key properties of PPN flows. Our simulations address magnetic tower evolution and stability. We also present results of a campaign of simulations to explore the development of accretion disks formed via wind capture. Our result focus on the limits of disk formation and the range of disk properties.

An emerging theme in modern astrophysics is the connection between astronomical observations and the underlying physical phenomena that drive our cosmos. Both the mechanisms responsible for the observed astrophysical phenomena and the tools used to probe such phenomena - the radiation and particle spectra we observe - have their roots in atomic, molecular, condensed matter, plasma, nuclear and particle physics. Chemistry is implicitly included in both molecular and condensed matter physics. This connection is the theme of the present report, which provides a broad, though non-exhaustive, overview of progress in our understanding of the cosmos resulting from recent theoretical and experimental advances in what is commonly called laboratory astrophysics. This work, carried out by a diverse community of laboratory astrophysicists, is increasingly important as astrophysics transitions into an era of precise measurement and high fidelity modeling.

Astrophysical jets have been studied with observations, theoretical models and numerical simulations for decades. Recently, supersonic magnetized jets have been formed in laboratory experiments of high-energy density plasmas. I will review these studies and discuss the experimental setup that has been used to form millimeter-scale jets driven by strong toroidal magnetic fields in a MAGPIE generator. The physical conditions of these experiments are such that they can be scaled to astrophysical scenarios. These laboratory jets provide insights on the underlying physics of magnetic tower jets and help constrain some models of astrophysical jets. In this context, we also discuss the connection between the laboratory jets and recent 3D-MHD numerical simulations of Poynting flux dominated jets. The simulations allow us to investigate the effects of thermal energy losses and base rotation on the growth rate of kink mode perturbations, and to compare the evolution of PFD jets with a hydrodynamic counterpart of the same energy flux.

We present 3D-MHD AMR simulations of Poynting flux dominated (PFD) jets formed by injection of magnetic energy. We compare their evolution with a hydrodynamic jet which is formed by injecting kinetic energy with the same energy flux than the PFD jets. We predict characteristic emission distributions for each of these jets. Current-driven perturbations in PFD jets are amplified by both cooling and rotation for the regimes studied: Shocks and thermal pressure support are weakened by cooling, making the jets more susceptible to kinking. Rotation amplifies the toroidal magnetic field which also exacerbates the kink instability.

We present four-color images of CRL 2688 obtained in 2009 using the Wide-Field Camera 3 on HST. The F606W image is compared with archival images in very similar filters to monitor the proper motions of nebular structure. We find that the bright N-S lobes have expanded uniformly by 2.5% and that the ensemble of rings has translated radially by 0.07 in 6.65 y. The rings were ejected every 100y for ~4 millennia until the lobes formed 250y ago. Starlight scattered from the edges of the dark E-W dust lane is coincident with extant H2 images and leading tips of eight pairs of CO outflows. We interpret this as evidence that fingers lie within geometrically opposite cones of opening angles {\approx} 30{\circ} like those in CRL618. By combining our results of the rings with 12CO absorption from the extended AGB wind we ascertain that the rings were ejected at ~18 km s-1 with very little variation and that the distance to CRL2688, v_{exp}$ / ${\dot\theta}_exp$, is 300 - 350 pc. Our 2009 imaging program included filters that span 0.6 to 1.6{\mu}m. We constructed a two-dimensional dust scattering model of stellar radiation through CRL2688 that successfully reproduces the details of the nebular geometry, its integrated spectral energy distribution, and nearly all of its color variations. The model implies that the optical opacity of the lobes >~ 1, the dust particle density in the rings decreases as radius^{-3} and that the mass and momentum of the AGB winds and their rings have increased over time.

Current AMR simulations require algorithms that are highly parallelized and manage memory efficiently. As compute engines grow larger, AMR simulations will require algorithms that achieve new levels of efficient parallelization and memory management. We have attempted to employ new techniques to achieve both of these goals. Patch or grid based AMR often employs ghost cells to decouple the hyperbolic advances of each grid on a given refinement level. This decoupling allows each grid to be advanced independently. In AstroBEAR we utilize this independence by threading the grid advances on each level with preference going to the finer level grids. This allows for global load balancing instead of level by level load balancing and allows for greater parallelization across both physical space and AMR level. Threading of level advances can also improve performance by interleaving communication with computation, especially in deep simulations with many levels of refinement. To improve memory management we have employed a distributed tree algorithm that requires processors to only store and communicate local sections of the AMR tree structure with neighboring processors.

Interfaces between hot and cold magnetized plasmas exist in various astrophysical contexts, for example where hot outflows impinge on an ambient interstellar medium (ISM). It is of interest to understand how the structure of the magnetic field spanning the interface affects the temporal evolution of the temperature gradient. Here we explore the relation between the magnetic field topology and the heat transfer rate by adding various fractions of tangled vs. ordered field across a hot-cold interface allow the system to evolve to a steady state. We find a simple mathematical relation for the rate of heat conduction as a function of the initial ratio of ordered to tangled field across the interface. We discuss potential implications for the astrophysical context of magnetized wind blown bubbles (WBB) around evolved stars.

We present new, third-epoch HST H-alpha and [S II] images of three HH jets (HH 1 and 2, HH 34, and HH 47) and compare these images with those from the previous epochs. The high-spatial resolution, coupled with a time-series whose cadence is of order both the hydrodynamical and radiative cooling timescales of the flow allows us to follow the hydrodynamical/magnetohydrodynamical evolution of an astrophysical plasma system in which ionization and radiative cooling play significant roles. Cooling zones behind the shocks are resolved, so it is possible to identify which way material flows through a given shock wave. The images show that heterogeneity is paramount in these jets, with clumps dominating the morphologies of both bow shocks and their Mach disks. This clumpiness exists on scales smaller than the jet widths and determines the behavior of many of the features in the jets. Evidence also exists for considerable shear as jets interact with their surrounding molecular clouds, and in several cases we observe shock waves as they form and fade where material emerges from the source and as it proceeds along the beam of the jet. Fine-structure within two extended bow shocks may result from Mach stems that form at the intersection points of oblique shocks within these clumpy objects. Taken together, these observations represent the most significant foray thus far into the time domain for stellar jets, and comprise one of the richest data sets in existence for comparing the behavior of a complex astrophysical plasma flows with numerical simulations and laboratory experiments.

We investigate the role discrete clumps embedded in an astrophysical jet play on the jet's morphology and line emission characteristics. By varying clumps' size, density, position, and velocity, we cover a range of parameter space motivated by observations of objects such as the Herbig Haro object HH~34. We here extend the results presented in Yirak et al. 2009, including how analysis of individual observations may lead to spurious sinusoidal variation whose parameters vary widely over time, owing chiefly to interacts between clumps. The goodness of the fits, while poor in all simulations, are best when clump-clump collisions are minimal. Our results indicate that a large velocity dispersion leads to a clump-clump collision-dominated flow which disrupts the jet beam. Finally, we present synthetic emission images of H-$\alpha$ and [SII] and note an excess of [SII] emission along the jet length as compared to observations. This suggests that observed beams undergo earlier processing, if they are present at all.

We present experimental results on the formation of supersonic, radiatively cooled jets driven by pressure due to the toroidal magnetic field generated by the 1.5 MA, 250 ns current from the MAGPIE generator. The morphology of the jet produced in the experiments is relevant to astrophysical jet scenarios in which a jet on the axis of a magnetic cavity is collimated by a toroidal magnetic field as it expands into the ambient medium. The jets in the experiments have similar Mach number, plasma beta and cooling parameter to those in protostellar jets. Additionally the Reynolds, magnetic Reynolds and Peclet numbers are much larger than unity, allowing the experiments to be scaled to astrophysical flows. The experimental configuration allows for the generation of episodic magnetic cavities, suggesting that periodic fluctuations near the source may be responsible for some of the variability observed in astrophysical jets. Preliminary measurements of kinetic, magnetic and Poynting energy of the jets in our experiments are presented and discussed, together with estimates of their temperature and trapped toroidal magnetic field.

The explosive BN/KL outflow emerging from OMC1 behind the Orion Nebula may have been powered by the dynamical decay of a non-hierarchical multiple system $\sim$500 years ago that ejected the massive stars I, BN, and source n, with velocities of about 10 to 30 km s$^{-1}$. New proper motion measurements of H$_2$ features show that within the errors of measurement, the outflow originated from the site of stellar ejection. Combined with published data, these measurements indicate an outflow age of $\sim$500 years, similar to the time since stellar ejection. The total kinetic energy of the ejected stars and the outflow is about 2 to $6 \times 10^{47}$ ergs. It is proposed that the gravitational potential energy released by the formation of a short-period binary, most likely source I, resulted in stellar ejection and powered the outflow. A scenario is presented for the formation of a compact, non-hierarchical multiple star system, its decay into an ejected binary and two high-velocity stars, and launch of the outflow. Three mechanisms may have contributed to the explosion in the gas: (i) Unbinding of the circum-cluster envelope following stellar ejection, (ii) disruption of circumstellar disks and high-speed expulsion of the resulting debris during the final stellar encounter, and (iii) the release of stored magnetic energy. Plausible proto-stellar disk end envelope properties can produce the observed outflow mass, velocity, and kinetic energy distributions. The ejected stars may have acquired new disks by fall-back or Bondi-Hoyle accretion with axes roughly orthogonal to their velocities. The expulsion of gas and stars from OMC1 may have been driven by stellar interactions.

Using 3D-MHD Eulerian-grid numerical simulations, we study the formation and evolution of rising magnetic towers propagating into an ambient medium. The towers are generated from a localized injection of pure magnetic energy. No rotation is imposed on the plasma. We compare the evolution of a radiatively cooling tower with an adiabatic one, and find that both bend due to pinch instabilities. Collimation is stronger in the radiative cooling case; the adiabatic tower tends to expand radially. Structural similarities are found between these towers and the millimeter scale magnetic towers produced in laboratory experiments.

Proto-planetary nebulae (pPN) and planetary nebulae (PN) seem to be formed by interacting winds from asymptotic giant branch (AGB) stars. The observational issue that most pPN are bipolar but most older PN are elliptical is addressed. We present 2.5D hydrodynamical numerical simulations of episodic cooling interacting winds to investigate the long term evolution of PN morphologies. We track wind acceleration, decrease in mass-loss and episodic change in wind geometry from spherical (AGB) to collimated (pPN) and back to spherical again (PN). This outflow sequence is found to produce realistic PN dynamics and morphological histories. Effects from different AGB distributions and jet duty cycles are also investigated.

Feedback from protostellar outflows can influence the nature of turbulence in star forming regions even if they are not the primary source of velocity dispersion for all scales of molecular clouds. For the rate and power expected in star forming regions, we previously (Carroll et al. 2009) demonstrated that outflows could drive supersonic turbulence at levels consistent with the scaling relations from Matzner 2007 although with a steeper velocity power spectrum than expected for an isotropically driven supersonic turbulent cascade. Here we perform higher resolution simulations and combine simulations of outflow driven turbulence with those of isotropically forced turbulence. We find that the presence of outflows within an ambient isotropically driven turbulent environment produces a knee in the velocity power spectrum at the outflow scale and a steeper slope at sub-outflow scales than for a purely isotropically forced case. We also find that the presence of outflows flattens the density spectrum at large scales effectively reducing the formation of large scale turbulent density structures. These effects are qualitatively independent of resolution. We have also carried out Principal Component Analysis (PCA) for synthetic data from our simulations. We find that PCA as a tool for identifying the driving scale of turbulence has a misleading bias toward low amplitude large scale velocity structures even when they are not necessarily the dominant energy containing scales. This bias is absent for isotropically forced turbulence but manifests strongly for collimated outflow driven turbulence.

Numeric convergence studies demonstrate that the evolution of an adiabatic clump is well-captured by roughly 100 cells per clump radius. The presence of radiative cooling, however, imposes limits on the problem due to the removal of thermal energy. Numerical studies which include radiative cooling typically adopt the 100--200 cells per clump radius resolution. In this paper we present the results of a convergence study for radiatively cooling clumps undertaken over a broad range of resolutions, from 12 to 1,536 cells per clump radius, employing adaptive mesh refinement (AMR) in a 2D axisymmetric geometry ("2.5D"). We also provide a fully 3D simulation, at 192 cells per clump radius, which supports our 2.5D results. We find no appreciable self-convergence at ~100 cells per clump radius as small-scale differences owing to increasingly resolving the "cooling length" have global effects. We therefore conclude that self-convergence is an insufficient criterion to apply on its own when addressing the question of sufficient resolution for radiatively cooled shocked clump simulations. We suggest the adoption of alternate criteria to support a statement of sufficient resolution, such as the demonstration of adequate resolution of the cooling layers behind shocks. We discuss an associated refinement criteria for AMR codes.

The WGLA of the AAS (http://www.aas.org/labastro/) promotes collaboration and exchange of knowledge between astronomy and planetary sciences and the laboratory sciences (physics, chemistry, and biology). Laboratory data needs of ongoing and next generation planetary science missions are carefully evaluated and recommended in this white paper submitted by the WGLA to Planetary Decadal Survey.

Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets, and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. Laboratory jets can be studied for a variety of initial conditions, arbitrary viewing angles, and different times, attributes especially helpful for interpreting astronomical images where the viewing angle and initial conditions are fixed and the time domain is limited. Experiments are also a powerful way to test numerical fluid codes in a parameter range where the codes must perform well. In this paper we combine images from a series of laboratory experiments of deflected supersonic jets with numerical simulations and new spectral observations of an astrophysical example, the young stellar jet HH 110. The experiments provide key insights into how deflected jets evolve in 3-D, particularly within working surfaces where multiple subsonic shells and filaments form, and along the interface where shocked jet material penetrates into and destroys the obstacle along its path. The experiments also underscore the importance of the viewing angle in determining what an observer will see. The simulations match the experiments so well that we can use the simulated velocity maps to compare the dynamics in the experiment with those implied by the astronomical spectra. The experiments support a model where the observed shock structures in HH 110 form as a result of a pulsed driving source rather than from weak shocks that may arise in the supersonic shear layer between the Mach disk and bow shock of the jet's working surface.

We present results of a series of magnetohydrodynamic (MHD) and hydro- dynamic (HD) 2.5D simulations of the morphology of outflows driven by nested wide-angle winds - i.e. winds which eminate from a central star as well as from an orbiting accretion disk. While our results are broadly relevent to nested wind systems we have tuned the parameters of the simulations to touch on issues in both Young Stellar Objects and Planetary Nebula studies. In particular our studies connect to open issues in the early evolution of Planetary Nebulae. We find that nested MHD winds exhibit marked morphological differences from the single MHD wind case along both dimensions of the flow. Nested HD winds on the other hand give rise mainly to geometric distortions of an outflow that is topologically similar to the flow arising from a single stellar HD wind. Our MHD results are insensitive to changes in ambient temperature between ionized and un-ionized circumstellar environments. The results are sensitive to the relative mass-loss rates, and to the relative speeds of the stellar and disk winds. We also present synthetic emission maps of both nested MHD and HD simulations. We find that nested MHD winds show knots of emission appearing on-axis that do not appear in the HD case.

We report on experiments in which magnetically driven radiatively cooled plasma jets were produced by a 1 MA, 250 ns current pulse on the MAGPIE pulsed power facility. The jets were driven by the pressure of a toroidal magnetic field in a ''magnetic tower'' jet configuration. This scenario is characterized by the formation of a magnetically collimated plasma jet on the axis of a magnetic ''bubble'', confined by the ambient medium. The use of a radial metallic foil instead of the radial wire arrays employed in our previous work allows for the generation of episodic magnetic tower outflows which emerge periodically on timescales of ~30 ns. The subsequent magnetic bubbles propagate with velocities reaching ~300 km/s and interact with previous eruptions leading to the formation of shocks.

Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The mission enabling impact of laboratory astrophysics ranges from the scientific conception stage for airborne and space-based observatories, all the way through to the scientific return of these missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA. These efforts are necessary for the success of astronomical research being funded by NASA. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics (WGLA), have prepared this White Paper on the laboratory astrophysics infrastructure needed to maximize the scientific return from NASA's space and Earth sciences program.

A description is given of the algorithms implemented in the AstroBEAR adaptive mesh refinement code for ideal magnetohydrodynamics. The code provides several high resolution, shock capturing schemes which are constructed to maintain conserved quantities of the flow in a finite volume sense. Divergence free magnetic field topologies are maintained to machine precision by collating the components of the magnetic field on a cell-interface staggered grid and utilizing the constrained transport approach for integrating the induction equations. The maintenance of magnetic field topologies on adaptive grids is achieved using prolongation and restriction operators which preserve the divergence and curl of the magnetic field across co-located grids of different resolution. The robustness and correctness of the code is demonstrated by comparing the numerical solution of various tests with analytical solutions or previously published numerical solutions obtained by other codes.

Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The impact of laboratory astrophysics ranges from the scientific conception stage for ground-based, airborne, and space-based observatories, all the way through to the scientific return of these projects and missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA, NSF, and DOE. These efforts are necessary for the success of astronomical research being funded by the agencies. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics, have prepared this State of the Profession Position Paper on the laboratory astrophysics infrastructure needed to ensure the advancement of astronomy and astrophysics in the next decade.

Almost all stars in the 1-8 Msun range evolve through the Asymptotic Giant Branch (AGB), preplanetary nebula (PPN) and planetary nebula (PN) evolutionary phases. Most stars that leave the main sequence in a Hubble time will end their lives in this way. The heavy mass loss which occurs during the AGB phase is important across astrophysics, and the particulate matter crucial for the birth of new solar systems is made and ejected by AGB stars. Yet stellar evolution from the beginning of the AGB phase to the PN phase remains poorly understood. We do not understand how the mass-loss (rate, geometry, temporal history) depends on fundamental stellar parameters or the presence of a binary companion. While the study of evolved non-massive stars has maintained a relatively modest profile in recent decades, we are nonetheless in the midst of a quiet but exciting revolution in this area, driven by new observational results, such as the discovery of jets and disks in stellar environments where these were never expected, and by the recognition of new symmetries such as multipolarity and point-symmetry occuring frequently in the nebulae resulting from the outflows. In this paper we summarise the major unsolved problems in this field, and specify the areas where allocation of effort and resources is most likely to help make significant progress.

As the panel on Planetary Systems and Star Formation (PSF) is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of solar system bodies (other than the Sun) and extrasolar planets, debris disks, exobiology, the formation of individual stars, protostellar and protoplanetary disks, molecular clouds and the cold ISM, dust, and astrochemistry. Central to the progress in these areas are the corresponding advances in laboratory astro- physics which are required for fully realizing the PSF scientific opportunities in the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics and chemistry which produce the observed spectra and describe the astrophysical processes. We discuss four areas of laboratory astrophysics relevant to the PSF panel: atomic, molecular, solid matter, and plasma physics. Section 2 describes some of the new opportunities and compelling themes which will be enabled by advances in laboratory astrophysics. Section 3 provides the scientific context for these opportunities. Section 4 discusses some experimental and theoretical advances in laboratory astrophysics required to realize the PSF scientific opportunities of the next decade. As requested in the Call for White Papers, we present in Section 5 four central questions and one area with unusual discovery potential. We give a short postlude in Section 6.

As the Galactic Neighborhood (GAN) panel is fully aware, the next decade will see major advances in our understanding of this area of research. To quote from their charge, these advances will occur in studies of the galactic neighborhood, including the structure and properties of the Milky Way and nearby galaxies, and their stellar populations and evolution, as well as interstellar media and star clusters. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GAN scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics and chemistry that produces the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the GAN panel are atomic, molecular, solid matter, plasma, and nuclear physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GAN scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

As the Cosmology and Fundamental Physics (CFP) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the early universe, the microwave background, the reionization and galaxy formation up to virialization of protogalaxies, large scale structure, the intergalactic medium, the determination of cosmological parameters, dark matter, dark energy, tests of gravity, astronomically determined physical constants, and high energy physics using astronomical messengers. Central to the progress in these areas are the corresponding advances in laboratory astrophysics which are required for fully realizing the CFP scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics which produce the observed astrophysical processes. The 5 areas of laboratory astrophysics which we have identified as relevant to the CFP panel are atomic, molecular, plasma, nuclear, and particle physics. Here, Section 2 describes some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the CFP scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

As the Galaxies across Cosmic Time (GCT) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the formation, evolution, and global properties of galaxies and galaxy clusters, as well as active galactic nuclei and QSOs, mergers, star formation rate, gas accretion, and supermassive black holes. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GCT scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produce the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the CFP panel are atomic, molecular, solid matter, plasma, nuclear, and particle physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GCT scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

As the Stars and Stellar Evolution (SSE) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the Sun as a star, stellar astrophysics, the structure and evolution of single and multiple stars, compact objects, SNe, gamma-ray bursts, solar neutrinos, and extreme physics on stellar scales. Central to the progress in these areas are the corresponding advances in laboratory astrophysics, required to fully realize the SSE scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produces the observed astrophysical processes. The 6 areas of laboratory astrophysics, which we have identified as relevant to the CFP panel, are atomic, molecular, solid matter, plasma, nuclear physics, and particle physics. In this white paper, we describe in Section 2 the scientific context and some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we discuss some of the experimental and theoretical advances in laboratory astrophysics required to realize the SSE scientific opportunities of the next decade. As requested in the Call for White Papers, Section 4 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 5.

We report the serendipitous detection by the XMM-Newton X-ray Observatory of an X-ray source at the position of the Type I (He- and N-rich) bipolar planetary nebula Hb 5. The Hb 5 X-ray source appears marginally resolved. While the small number of total counts (~170) and significant off-axis angle of the X-ray source (~7.8') precludes a definitive spatial analysis, the morphology of the X-ray emission appears to trace the brightest features seen in optical images of Hb 5. The X-ray spectrum is indicative of a thermal plasma at a temperature between 2.4 and 3.7 MK and appears to display strong Neon emission. The inferred X-ray luminosity is L_X = 1.5 x 10^32 ergs/s. These results suggest that the detected X-ray emission is dominated by shock-heated gas in the bipolar nebula, although we cannot rule out the presence of a point-like component at the position of the central star. The implications for and correspondence with current models of shock-heated gas in planetary nebulae is discussed.

Collimated outflows (jets) are ubiquitous in the universe appearing around sources as diverse as protostars and extragalactic supermassive blackholes. Jets are thought to be magnetically collimated, and launched from a magnetized accretion disk surrounding a compact gravitating object. We have developed the first laboratory experiments to address time-dependent, episodic phenomena relevant to the poorly understood jet acceleration and collimation region. The experimental results show the periodic ejections of magnetic bubbles naturally evolving into a heterogeneous jet propagating inside a channel made of self-collimated magnetic cavities. The results provide a unique view of the possible transition from a relatively steady-state jet launching to the observed highly structured outflows.

Herbig-Haro (HH) jets are commonly thought of as homogeneous beams of plasma traveling at hypersonic velocities. Structure within jet beams is often attributed to periodic or ``pulsed'' variations of conditions at the jet source. Simulations based on this scenario result in knots extending across the jet diameter. Observations and recent high energy density laboratory experiments shed new light on structures below this scale and indicate they may be important for understanding the fundamentals of jet dynamics. In this paper we offer an alternative to ``pulsed'' models of protostellar jets. Using direct numerical simulations we explore the possibility that jets are chains of sub-radial clumps propagating through a moving inter-clump medium. Our models explore an idealization of this scenario by injecting small ($r<r_{jet}$), dense ($\rho>\rho_{jet}$) spheres embedded in an otherwise smooth inter-clump jet flow. The spheres are initialized with velocities differing from the jet velocity by $\sim15$%. We find the consequences of shifting from homogeneous to heterogeneous flows are significant as clumps interact with each other and with the inter-clump medium in a variety of ways. Structures which mimic what is expected from pulsed-jet models can form, as can previously unseen ``sub-radial'' behaviors including backward facing bow shocks and off-axis working surfaces. While these small-scale structures have not been seen before in simulation studies, they are found in high resolution jet observations. We discuss implications of our simulations for the interpretation of protostellar jets with regard to characterization of knots by a ``lifetime'' or ``velocity history'' approach as well as linking observed structures with central engines which produce the jets.

The link between turbulence in star formatting environments and protostellar jets remains controversial. To explore issues of turbulence and fossil cavities driven by young stellar outflows we present a series of numerical simulations tracking the evolution of transient protostellar jets driven into a turbulent medium. Our simulations show both the effect of turbulence on outflow structures and, conversely, the effect of outflows on the ambient turbulence. We demonstrate how turbulence will lead to strong modifications in jet morphology. More importantly, we demonstrate that individual transient outflows have the capacity to re-energize decaying turbulence. Our simulations support a scenario in which the directed energy/momentum associated with cavities is randomized as the cavities are disrupted by dynamical instabilities seeded by the ambient turbulence. Consideration of the energy power spectra of the simulations reveals that the disruption of the cavities powers an energy cascade consistent with Burgers'-type turbulence and produces a driving scale-length associated with the cavity propagation length. We conclude that fossil cavities interacting either with a turbulent medium or with other cavities have the capacity to sustain or create turbulent flows in star forming environments. In the last section we contrast our work and its conclusions with previous studies which claim that jets can not be the source of turbulence.

In this paper we explore the relationship between protostellar outflows and turbulence in molecular clouds. Using 3-D numerical simulations we focus on the hydrodynamics of multiple outflows interacting within a parsec scale volume. We explore the extent to which transient outflows injecting directed energy and momentum into a sub-volume of a molecular cloud can be converted into random turbulent motions. We show that turbulence can readily be sustained by these interactions and show that it is possible to broadly characterize an effective driving scale of the outflows. We compare the velocity spectrum obtained in our studies to that of isotropically forced hydrodynamic turbulence finding that in outflow driven turbulence a power law is indeed achieved. However we find a steeper spectrum (beta ~ 3) is obtained in outflow driven turbulence models than in isotropically forced simulations (beta ~ 2). We discuss possible physical mechanisms responsible for these results as well and their implications for turbulence in molecular clouds where outflows will act in concert with other processes such as gravitational collapse.

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.

The first laboratory astrophysics experiments to produce a radiatively cooled plasma jet with dynamically significant angular momentum are discussed. A new configuration of wire array z-pinch, the twisted conical wire array, is used to produce convergent plasma flows each rotating about the central axis. Collision of the flows produces a standing shock and jet that each have supersonic azimuthal velocities. By varying the twist angle of the array, the rotation velocity of the system can be controlled, with jet rotation velocities reaching ~20% of the propagation velocity.

The study of PN has been confronting a growing list of dilemmas which have yet to find coherent resolution. These issues are both observational and theoretical and can be stated as a series of "facts" which can not, as of yet, be accounted for via a single framework. We review these facts and propose a skeleton framework for developing a new understanding post-AGB stars, PPN and PN. Our framework represents an attempt to articulate a a global perspective on the late stages of stellar evolution that can embrace both the nature of the central engine and the outflows they produce. Our framework focuses on interacting binary central stars which drive collimated outflows through MHD processes. We propose that the field of AGB/PN studies now faces a "Grand Challenge" in articulating the observational systematics of these objects in a way that can address issues related to binarity and magnetic shaping. A theoretical Grand Challenge is also faced in the form of integrated studies which can explicate the highly non-linear processes associated with MHD outflows driven by interacting binaries. These issues include the generation of magnetic fields via dynamo processes, the creation of accretion disks, the dynamics of Common Envelope ejection and the creation of magnetized jets.

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.

Many proto-planetary nebulae (PPN) appear as narrow collimated structures sometimes showing multiple, roughly aligned lobes. In addition, many PPN flows have been shown to have short acceleration times. In this paper we explore whether jet or ``bullet'' (a massive clump) models fit the observations of individual collimated lobes adequately by comparing simulations of both radiatively cooled (stable) jets and bullets. We find that the clump model is somewhat favored over jets because (1) it leads to greater collimation of outflows (2) it accounts better and more naturally for ring-like structures observed in the PPN CRL 618, and (3) it is more successful in reproducing the Hubble-flow character of observed kinematics in some PPN. In addition, bullets naturally account for observed multipolar flows, since the likely MHD launch mechanisms required to drive outflows make multiple non-aligned jets unlikely. Thus we argue that PPN outflows may be driven by explosive MHD launch mechanisms such as those discussed in the context of supernovae (SNe) and gamma-ray bursts(GRB).

We present computational results showing eccentricity growth in the inner portions of a protoplanetary disc. We attribute this to the evolving surface density of the disc. The planet creates a gap, which adjusts the balance between the 3:1 (eccentricity exciting) and 2:1 (eccentricity damping) resonances. The eccentricity of the inner disc can rise as high as 0.3, which is sufficient to cause it to be accreted onto the star. This offers an alternative mechanism for producing the large holes observed in the discs of CoKu Tau/4, GM Aur and DM Tau.

We have performed 2D hydrodynamic simulations of a pulsed astrophysical jet propagating through a medium that is populated with spherical inhomogeneities, or "clumps," which are smaller than the jet radius. The clumps are seen to affect the jet in several ways, such as impeding jet propagation and deflecting the jet off-axis. While there has been some debate as to the prevalence of these types of condensations in the ISM or in molecular clouds, the exploration of this region of parameter space nonetheless both shows the potential for these clumps to disrupt astrophysical jets and yields results which recover aspects of recent observations of Herbig-Haro objects. We find that the propagation of the jet and the vorticity induced in the clump/ambient medium correlate well with a "dynamic filling function" $f_d$ across all the simulations.

Although several lines of evidence suggest that jets from young stars are driven magnetically from accretion disks, existing observations of field strengths in the bow shocks of these flows imply that magnetic fields play only a minor role in the dynamics at these locations. To investigate this apparent discrepancy we performed numerical simulations of expanding magnetized jets with stochastically variable input velocities with the AstroBEAR MHD code. Because the magnetic field B is proportional to the density n within compression and rarefaction regions, the magnetic signal speed drops in rarefactions and increases in the compressed areas of velocity-variable flows. In contrast, B ~ n^0.5 for a steady-state conical flow with a toroidal field, so the Alfven speed in that case is constant along the entire jet. The simulations show that the combined effects of shocks, rarefactions, and divergent flow cause magnetic fields to scale with density as an intermediate power 1 > p > 0.5. Because p > 0.5, the Alfven speed in rarefactions decreases on average as the jet propagates away from the star. This behavior is extremely important to the flow dynamics because it means that a typical Alfven velocity in the jet close to the star is significantly larger than it is in the rarefactions ahead of bow shocks at larger distances, the one place where the field is a measurable quantity. We find that the observed values of weak fields at large distances are consistent with strong fields required to drive the observed mass loss close to the star. For a typical stellar jet the crossover point inside which velocity perturbations of 30 - 40 km/s no longer produce shocks is ~ 300 AU from the source.

The evolution of laboratory produced magnetic jets is followed numerically through three-dimensional, non-ideal magnetohydrodynamic simulations. The experiments are designed to study the interaction of a purely toroidal field with an extended plasma background medium. The system is observed to evolve into a structure consisting of an approximately cylindrical magnetic cavity with an embedded magnetically confined jet on its axis. The supersonic expansion produces a shell of swept-up shocked plasma which surrounds and partially confines the magnetic tower. Currents initially flow along the walls of the cavity and in the jet but the development of current-driven instabilities leads to the disruption of the jet and a re-arrangement of the field and currents. The top of the cavity breaks-up and a well collimated, radiatively cooled, 'clumpy' jet emerges from the system.

We investigate the evolution of fossil cavities produced by extinct YSO jets and wide angle outflows. Fossil cavities are ellipsoidal or cylindrical shells of swept-up of ambient (molecular cloud) material moving at low velocities. The cavities form when the momentum in a YSO jet or wide angle outflow decays in time allowing the bowshock or swept-up shell to decelerate to velocities near the turbulent speed in the cloud. It has been suggested in previous studies that cavities provide efficient coupling between the jets/outflows and the cloud and, as such, are the agents by which cloud turbulence is can be re-energized. In this paper we carry forward a series of numerical simulations of jets and outflows whose momentum flux decrease in time. We compare simulations with decaying momentum fluxes to those with constant flux. We show that decaying flux models exhibit deceleration of the outflow head and backfilling via expansion off of the cavity walls. They also have lower density contrast, are longer lived and wider than their continuously driven counterparts. The simulations recover the basic properties of observed fossil cavities. In addition, we provide synthetic observations in terms of P-V diagrams which demonstrate that fossil cavities form both jets and wide angle outflows are characterized by linear "Hubble-law" expansions patterns superimposed on "spur" patterns indicative of the head of a bow shock.

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.

We present Chandra X-ray Observatory monitoring observations of the recent accretion outburst displayed by the pre-main sequence (pre-MS) star V1647 Ori. The X-ray observations were obtained over a period beginning prior to outburst onset in late 2003 and continuing through its apparent cessation in late 2005, and demonstrate that the mean flux of the spatially coincident X-ray source closely tracked the near-infrared luminosity of V1647 Ori throughout its eruption. We find negligible likelihood that the correspondence between X-ray and infrared light curves over this period was the result of multiple X-ray flares unrelated to the accretion burst. The recent Chandra data confirm that the X-ray spectrum of V1647 Ori hardened during outburst, relative both to its pre-outburst state and to the X-ray spectra of nearby pre-MS stars in the L1630 cloud. We conclude that the observed changes in the X-ray emission from V1647 Ori over the course of its 2003-2005 eruption were generated by a sudden increase and subsequent decline in its accretion rate. These results for V1647 Ori indicate that the flux of hard X-ray emission from erupting low-mass, pre-MS stars, and the duration and intensity of such eruptions, reflect the degree to which star-disk magnetic fields are reorganized before and during major accretion events.

We present an atlas of Spitzer Space Telescope Infrared Spectrograph (IRS) spectra of highly luminous, compact mid-infrared sources in the Large Magellanic Cloud. Sources were selected on the basis of infrared colors and 8 micron (MSX) fluxes indicative of highly evolved, intermediate- to high-mass stars with current or recent mass loss at large rates. We determine the chemistry of the circumstellar envelope from the mid-IR continuum and spectral features and classify the spectral types of the stars. In the sample of 60 sources, we find 21 Red Supergiants (RSGs), 16 C-rich Asymptotic Giant Branch (AGB) stars, 11 HII regions, 4 likely O-rich AGB stars, 4 Galactic O-rich AGB stars, 2 OH/IR stars, and 2 B[e] supergiants with peculiar IR spectra. We find that the overwhelming majority of the sample AGB stars (with typical IR luminosities ~1.0E4 L_sun) have C-rich envelopes, while the O-rich objects are predominantly luminous RSGs with L_IR ~ 1.0E5 L_sun. We determine mean bolometric corrections to the stellar K-band flux densities and find that for carbon stars, the bolometric corrections depend on the infrared color, whereas for RSGs, the bolometric correction is independent of IR color. Our results reveal that objects previously classified as PNe on the basis of IR colors are in fact compact HII regions with very red IRS spectra that include strong atomic recombination lines and PAH emission features. We demonstrate that the IRS spectral classes in our sample separate clearly in infrared color-color diagrams that use combinations of 2MASS data and synthetic IRAC/MIPS fluxes derived from the IRS spectra. On this basis, we suggest diagnostics to identify and classify, with high confidence levels, IR-luminous evolved stars and HII regions in nearby galaxies using Spitzer and near-infrared photometry.

Jets and outflows are thought to be an integral part of accretion phenomena and are associated with a large variety of objects. In these systems, the interaction of magnetic fields with an accretion disk and/or a magnetized central object is thought to be responsible for the acceleration and collimation of plasma into jets and wider angle flows. In this paper we present three-dimensional MHD simulations of magnetically driven, radiatively cooled laboratory jets that are produced on the MAGPIE experimental facility. The general outflow structure comprises an expanding magnetic cavity which is collimated by the pressure of an extended plasma background medium, and a magnetically confined jet which develops within the magnetic cavity. Although this structure is intrinsically transient and instabilities in the jet and disruption of the magnetic cavity ultimately lead to its break-up, a well collimated, knotty jet still emerges from the system; such clumpy morphology is reminiscent of that observed in many astrophysical jets. The possible introduction in the experiments of angular momentum and axial magnetic field will also be discussed.

We present the results of a numerical magnetohydrodynamic simulation that demonstrates a mechanism by which magnetic fields tap rotational energy of a stellar core and expel the envelope. Our numerical setup, designed to focus on the basic physics of the outflow mechanism, consists of a solid, gravitating sphere, which may represent the compact core of a star, surrounded by an initially hydrostatic envelope of ionized gas. The core is threaded by a dipolar magnetic field that also permeates the envelope. At the start of the simulation, the core begins to rotate at 10% of the escape speed. The magnetic field is sufficiently strong to drive a magneto-rotational explosion, whereby the entire envelope is expelled, confirming the expectation of analytical models. Furthermore, the dipolar nature of the field results in an explosion that is enhanced simultaneously along the rotation axis (a jet) and along the magnetic equator. While the initial condition is simplified, the simulation approximates circumstances that may arise in astrophysical objects such as Type II supernovae, gamma ray bursts, and proto-planetary nebulae.

We review work on the evolution of planetary nebulae and proto-planetaries via magneto-rotational mechanisms showing that a dynamo generated magnetic field can produce the energy and momentum needed to drive pPN and PNe outflows. Angular momentum considerations lead to the conclusion that single stars may not be capable of supporting strong fields for long times. Thus we take the working hypothesis that most PN may form via binary stars. We propose that the grand challenge for PN studies is fully understanding the diverse physical processes at work in binary late stage evolution including the development of disks, fields and outflows.

We discuss the formation of globules in planetary nebulae, typified by those observed in the Helix Nebula. We show that the properties of the globules, their number, mass, separation, and overall geometry strongly support a scenario in which globules are formed by the fragmentation of a swept-up shell as opposed to models in which the knots form in the AGB wind. We show that the RT or other instabilities which lead to the break-up of shells formed in the nebulae by fast winds or ionization fronts can produce arrays of globules with the overall geometry and within the mass range observed. We also show that the presence of a magnetic field in the circumstellar gas may play an important role in controlling the fragmentation process. Using field strengths measured in the precursor AGB envelopes, we find that close to the central star where the fields are relatively strong, the wavelengths of unstable MRT modes are larger than the shell dimensions, and the fragmentation of the shell is suppressed. The wavelength of the most unstable MRT mode decreases with increasing distance from the star, and when it becomes comparable to the shell thickness, it can lead to the sudden, rapid break-up of an accelerating shell. For typical nebula parameters, the model results in numerous fragments with a mass scale and a separation scale similar to those observed. Our results provide a link between global models of PN shaping in which shells form via winds and ionization fronts, and the formation of small scale structures in the nebulae.

We have produced an atlas of Spitzer Infrared Spectrograph (IRS) spectra of mass-losing, evolved stars in the Large Magellanic Cloud. These stars were selected to have high mass-loss rates and so contribute significantly to the return of processed materials to the ISM. Our high-quality spectra enable the determination of the chemistry of the circumstellar envelope from the mid-IR spectral features and continuum. We have classified the spectral types of the stars and show that the spectral types separate clearly in infrared color-color diagrams constructed from 2MASS data and synthetic IRAC/MIPS fluxes derived from our IRS spectra. We present diagnostics to identify and classify evolved stars in nearby galaxies with high confidence levels using Spitzer and 2MASS photometry. Comparison of the spectral classes determined using IRS data with the IR types assigned based on NIR colors also revealed a significant number of misclassifications and enabled us to refine the NIR color criteria resulting in more accurate NIR color classifications of dust-enshrouded objects.

We present a series of numerical studies of the interaction of colliding radiative, hydrodynamic young stellar outflows. We study the effect of the collision impact parameter on the acceleration of ambient material and the degree to which the flow is isotropized by the collision as a mechanism for driving turbulence in the parent molecular cloud. Our results indicate that the high degrees of compression of outflow material, achieved through radiative shocks near the vertex of the interaction, prevents the redirected outflow from spraying over a large spatial region. Furthermore, the collision reduces the redirected outflow's ability to entrain and impart momentum into the ambient cloud. Consideration of the probabilities of outflow collisions leads us to conclude that individual low velocity fossil outflows are the principle coupling between outflows and the cloud.

Recent spectral observations by the Spitzer Space Telescope (SST) reveal that some discs around young ($\sim {\rm few} \times 10^6$ yr old) stars have remarkably sharp transitions to a low density inner region in which much of the material has been cleared away. It has been recognized that the most plausible mechanism for the sharp transition at a specific radius is the gravitational influence of a massive planet. This raises the question of whether the planet can also account for the hole extending all the way to the star. Using high resolution numerical simulations, we show that Jupiter-mass planets drive spiral waves which create holes on time scales $\sim 10$ times shorter than viscous or planet migration times. We find that the theory of spiral-wave driven accretion in viscous flows by Takeuchi et al. (1996) can be used to provide a consistent interpretation of the simulations. In addition, although the hole surface densities are low, they are finite, allowing mass accretion toward the star. Our results therefore imply that massive planets can form extended, sharply bounded spectral holes which can still accommodate substantial mass accretion rates. The results also imply that holes are more likely than gaps for Jupiter mass planets around solar mass stars.

With 2D hydrodynamical simulations of disks perturbed externally by stars, brown dwarfs or planets we investigate possible scenarios that can account for the spiral structure in circumstellar disks. We consider two scenarios, spiral structure driven by an external bound planet or low mass star and that excited by a previous stellar close encounter or flyby. We find that both scenarios produce morphology similar to that observed in the outer disks of HD 141569A and HD 100546; moderately open 2-armed outer spiral structure. The outer two-armed spiral structure observed in the disk of HD 141569A is qualitatively reproduced with tidal perturbations from its companion binary HD 141569B,C on a prograde orbit near periapse. Our simulation accounts for the outer spiral arms, but is less successful than the secular model of Augereau and Papaloizou at matching the lopsidedness or asymmetry of the disk edge at 300AU. The disk has been previously truncated by the tidal force from the binary. A bound object (stellar or planetary) is unlikely to explain the spiral structure in HD 100546. A co-eval planet or brown dwarf in the disk of sufficient mass to account for the amplitude of the spiral structure would be detectable in NICMOS and STIS images, however existing images reveal no such object. A previous encounter could explain the observed structure, provided that the encounter occurred less than a few thousand year ago. The object responsible for causing the spiral structure must then be within a few arcminutes of the star. However, the USNO-B proper motion survey reveals no candidate object. Moreover, the probability that a field star encountered HD 100546 in the past few thousand years is very low.

We investigate gap formation in gaseous protostellar disks by a planet in a circular orbit in the limit of low disk viscosity. This regime may be appropriate to an aging disk after the epoch of planet formation. We find that the distance of planet to the gap outer boundary can be between the location of the $m=2$ and $m=1$ outer Lindblad resonances. This distance is weakly dependent upon both the planet's mass and disk viscosity. We find that the evolution of the disk edge takes place on two timescales. The first timescale is set by the spiral density waves driven by the planet. The second timescale depends on the viscosity of the disk. The disk approaches a state where the outward angular momentum flux caused by the disk viscosity is balanced by the dissipation of spiral density waves which are driven at the Lindblad resonances. This occurs inefficiently however because of the extremely low gas density near the planet. We find that the distance between the planet and the peak density at the disk outer edge is only weakly dependent on the viscosity and planet mass, however the ratio of the gas density near the planet to that in the disk (or the slope of density along the disk edge) is strongly dependent upon both quantities. We find that the disk density profile along the edge scales approximately with disk viscosity divided by the square of the planet mass. We account for this behavior with a simple scenario in which the dissipation of angular momentum from the spiral density waves is balanced against diffusion in the steep edge of the disk.

Outflows from young stellar objects have been identified as a possible source of turbulence in molecular clouds. To investigate the relationship between outflows, cloud dynamics and turbulence, we compare the kinematics of the molecular gas associated with NGC 1333, traced in 13CO(1-0), with the distribution of young stellar objects (YSOs) within. We find a velocity dispersion of ~ 1-1.6 km/s in 13CO that does not significantly vary across the cloud, and is uncorrelated with the number of nearby young stellar outflows identified from optical and submillimeter observations. However, from velocity channel maps we identify about 20 cavities or depressions in the 13CO intensity of scales > 0.1-0.2 pc and velocity widths 1-3 km/s. The cavities exhibit limb brightened rims in both individual velocity channel maps and position velocity diagrams, suggesting that they are slowly expanding. We interpret these cavities to be remnants of past YSO outflow activity: If these cavities are presently empty, they would fill in on time scales of a million years. This can exceed the lifetime of a YSO outflow phase, or the transit time of the central star through the cavity, explaining the the absence of any clear correlation between the cavities and YSO outflows. We find that the momentum and energy deposition associated with the expansion of the cavities is sufficient to power the turbulence in the cloud. In this way we conclude that the cavities are an important intermediary step between the conversion of YSO outflow energy and momentum into cloud turbulent motions.

We carry out high resolution simulations of the inner regions of a wide angle wind driven bipolar outflow using an Adaptive Mesh Refinement code. Our code follows H-He gas with molecular, atomic and ionic components and the associated time dependent molecular chemistry and ionization dynamics with radiative cooling. Our simulations explore the nature of the outflow when a spherical wind expands into a rotating, collapsing envelope. We compare with key observational properties of the outflow system of Source I in the BN/KL region. Our calculations show that the wind evacuates a bipolar outflow cavity in the infalling envelope. We find the head of the outflow to be unstable and that it rapidly fragments into clumps. We resolve the dynamics of the strong shear layer which defines the side walls of the cavity. We conjecture that this layer is the likely site of maser emission and examine its morphology and rotational properties. The shell of swept up ambient gas that delineates the cavity edge retains its angular momentum. This rotation is roughly consistent with that observed in the Source I SiO maser spots. The observed proper motions and line-of-sight velocity are approximately reproduced by the model. The cavity shell at the base of the flow assumes an X-shaped morphology which is also consistent with Source I. We conclude that the wide opening angle of the outflow is evidence that a wide-angle wind drives the Source I outflow and not a collimated jet.

We present the first results of high energy density laboratory astrophysics experiments which explore the evolution of collimated outflows and jets driven by a toroidal magnetic field. The experiments are scalable to astrophysical flows in that critical dimensionless numbers such as the Mach number, the plasma beta and the magnetic Reynolds number are all in the astrophysically appropriate ranges. Our experiments use the MAGPIE pulsed power machine and allow us to explore the role of magnetic pressure in creating and collimating the outflow as well as showing the creation of a central jet within the broader outflow cavity. We show that currents flow along this jet and we observe its collimation to be enhanced by the additional hoop stresses associated with the generated toroidal field. Although at later times the jet column is observed to go unstable, the jet retains its collimation. We also present simulations of the magnetic jet evolution using our two-dimensional resistive magneto-hydrodynamic (MHD) laboratory code. We conclude with a discussion of the astrophysical relevance of the experiments and of the stability properties of the jet.

We report a ~38 ks X-ray observation of McNeil's Nebula obtained with XMM on 2004 April 4. V1647 Ori, the young star in outburst illuminating McNeil's Nebula, is detected with XMM and appears variable in X-rays. We investigate the hardness ratio variability and time variations of the event energy distribution with quantile analysis, and show that the large increase of the count rate from V1647 Ori observed during the second half of the observation is not associated with any large plasma temperature variations as for typical X-ray flares from young low-mass stars. X-ray spectral fitting shows that the bulk (~75%) of the intrinsic X-ray emission in the 0.5-8 keV energy band comes from a soft plasma component (0.9 keV) reminiscent of the X-ray spectrum of the classical T Tauri star TW Hya, for which X-ray emission is believed to be generated by an accretion shock onto the photosphere of a low-mass star. The hard plasma component (4.2 keV) contributes ~25% of the total X-ray emission, and can be understood only in the framework of plasma heating sustained by magnetic reconnection events. We find a hydrogen column density of NH=4.1E22 cm-2, which points out a significant excess of hydrogen column density compared to the value derived from optical/IR observations, consistent with the picture of the rise of a wind/jet unveiled from ground optical spectroscopy. The X-ray flux observed with XMM ranges from roughly the flux observed by Chandra on 2004 March 22 (~10 times greater than the pre-outburst X-ray flux) to a value two times greater than that caught by Chandra on 2004 March 7 (~200 times greater than the pre-outburst X-ray flux). We have investigated the possibility that V1647 Ori displays a periodic variation in X-ray brightness as suggested by the combined Chandra+XMM data set (abridged).

Young, low-mass stars are luminous X-ray sources whose powerful X-ray flares may exert a profound influence over the process of planet formation. The origin of such emission is uncertain. Although many or perhaps most recently formed, low-mass stars emit X-rays as a consequence of solar-like coronal activity, it has also been suggested that X-ray emission may be a direct result of mass accretion onto the forming star. Here we report X-ray imaging spectroscopy observations which reveal a factor ~50 increase in the X-ray flux from a young star that is presently undergoing a spectacular optical/IR outburst. The outburst is thought to be due to the sudden onset of a phase of rapid accretion. The coincidence of a surge in X-ray brightness with the optical/IR eruption demonstrates that strongly enhanced high-energy emission from young stars can occur as a consequence of high accretion rates. We suggest that such accretion- enhanced X-ray emission from erupting young stars may be short-lived, because intense star-disk magnetospheric interactions are quenched rapidly by the subsequent accretion flood.

Spitzer observations of the young star CoKuTau/4 reveal a disk with a 10 AU hole that is most likely caused by a newly formed planet. Assuming that the planet opened a gap in the viscous disk, we estimate that the planet mass is greater than 0.1 Jupiter masses. This estimate depends on a lower limit to the disk viscosity derived from the time scale needed to accrete the inner disk, creating the now detectable hole. The planet migration time scale must at least modestly exceed the time for the spectrally inferred hole to clear. The proximity of the planet to the disk edge implied by our limits suggests that the latter is perturbed by the nearby planet and may exhibit a spiral pattern rotating with the planet. This pattern might be resolved with current ground based mid-infrared cameras and optical cameras on the Hubble Space Telescope. The required sub-Myr planet formation may challenge core accretion formation models. However, we find that only if the planet mass is larger than about 10 Jupiter masses, allowing for a high enough surface density without inducing migration, would formation by direct gravitational instability be possible.

We present two dimensional cylindrically symmetric hydrodynamic simulations and synthetic emission maps of a stellar wind propagating into an infalling, rotating environment. The resulting outflow morphology, collimation and stability observed in these simulations have relevance to the study of young stellar objects, Herbig-Haro jets and molecular outflows. Our code follows hydrogen gas with molecular, atomic and ionic components tracking the associated time dependent molecular chemistry and ionization dynamics with radiative cooling appropriate for a dense molecular gas. We present tests of the code as well as new simulations which indicate the presence of instabilities in the wind-blown bubble's swept-up shell.

We present new data from High Energy Density (HED) laboratory experiments designed to explore the interaction of a heavy hypersonic radiative jet with a cross wind. The jets are generated with the MAGPIE pulsed power machine where converging conical plasma flows are produced from a cylindrically symmetric array of inclined wires. Radiative hypersonic jets emerge from the convergence point. The cross wind is generated by ablation of a plastic foil via soft-X-rays from the plasma convergence region. Our experiments show that the jets are deflected by the action of the cross wind with the angle of deflection dependent on the proximity of the foil. Shocks within the jet beam are apparent in the data. Analysis of the data shows that the interaction of the jet and cross wind is collisional and therefore in the hydro-dynamic regime. MHD plasma code simulations of the experiments are able to recover the deflection behaviour seen in the experiments. We consider the astrophysical relevance of these experiments applying published models of jet deflection developed for AGN and YSOs. Fitting the observed jet deflections to quadratic trajectories predicted by these models allows us to recover a set of plasma parameters consistent with the data. We also present results of 3-D numerical simulations of jet deflection using a new astrophysical Adaptive Mesh Refinement code. These simulations show highly structured shocks occurring within the beam similar to what was observed in the experiments

Winds from accretion disks have been proposed as the driving source for precessing jets and extreme bipolar morphologies in Planetary Nebulae (PNe) and proto-PNe (pPNe). In this paper we address the applicability MHD disk wind models to PNe and pPNe. We first review Pelletier & Pudritz's (1992) non self-similar approach to magneto-centrifugal disk winds. We then discuss conditions which may occur in PNe and pPNe accretion disks that form via binary interactions. Finally, examining conditions on dimensionless parameters needed for the magneto-centrifugal disk wind model, we show that such winds can recover the observed momentum and energy input rates for PNe and pPNe. High accretion rates ($M_a \approx 10^{-4}$ \mdot) may be required in the latter case. We find that the observed {\it total} energy and momentum in pPNe can be recovered with disk wind models using existing disk formation scenarios. When combined with existing scenarios for accretion disk formation from disrupted stellar companions, our models may provide an explanation for the existence of high speed polar knots (FLIERS) observed in some PNe.

We report the first detection of X-ray emission in a pre-planetary nebula, Hen 3-1475. Pre-planetary nebulae are rare objects in the short transition stage between the Asymptotic Giant Branch and planetary nebula evolutionary phases, and Hen 3-1475, characterised by a remarkable S-shaped chain of optical knots, is one of the most noteworthy members of this class. Observations with the Advanced CCD Imaging Spectrometer (ACIS) onboard the Chandra X-Ray observatory show the presence of compact emission coincident with the brightest optical knot in this bipolar object, which is displaced from the central star by 2.7 arcsec along the polar axis. Model fits to the X-ray spectrum indicate an X-ray temperature and luminosity, respectively, of (4.3-5.7) 10^6 K and (4+/-1.4) 10^{31} (D/5 kpc)^2 erg s^{-1}, respectively. Our 3-sigma upper limit on the luminosity of compact X-ray emission from the central star in Hen 3-1475 is ~5 10^{31} (D/5 kpc)^2 erg s^{-1}. The detection of X-rays in Hen 3-1475 is consistent with models in which fast collimated post-AGB outflows are crucial to the shaping of planetary nebulae; we discuss such models in the context of our observations.

Winds from accretion disks have been proposed as the driving source for precessing jets and extreme bipolar morphologies in Planetary Nebulae (PNe) and proto-PNe (pPNe). Here we apply MHD disk wind models to PNe and pPNe by estimating separately the asymptotic MHD wind velocities and mass loss rates. We show that the resulting winds can recover the observed momentum and energy input rates for PNe and pPNe.

Many astrophysical flows occur in inhomogeneous media. We briefly discuss some general properties of the adiabatic and radiative inhomogeneous systems and discuss the relevance of those properties to the planetary nebulae systems. We then focus on radiative hypersonic bullets and the applicability of this model to planetary and protoplanetary systems such as CRL 618, NGC 6543, Hen 3-1475.

The Accretion-Ejection Instability (AEI) has been proposed to explain the low frequency Quasi-Periodic Oscillation (QPO) observed in low-mass X-Ray Binaries. Its frequency, typically a fraction of the Keplerian frequency at the disk inner radius, is in the right range indicated by observations. With numerical simulation we will show how this instability is able to produce a modulation of the X-ray flux and what are the characteristic of this modulation. More simulations are required, especially 3D MHD simulations. We will briefly present a new code in development: AstroBear which will allow us to create synthetic spectra.

We present the results of a numerical study focusing on the propagation of a hypersonic bullet subject to radiative cooling. Our goal is to explore the feasibility of such a model for the formation of ``strings'' observed in the Eta Carinae Homunculus nebula. Our simulations were performed in cylindrical symmetry with the adaptive mesh refinement code AstroBEAR. The radiative cooling of the system was followed using the cooling curve by Dalgarno & McCray (1972). In this letter we discuss the evolution and overall morphology of the system as well as key kinematic properties. We find that radiative bullets can produce structures with properties similar to those of the Eta Carinae strings, i.e. high length-to-width ratios and Hubble-type flows in the form of a linear velocity increase from the base of the wake to the bullet head. These features, along with the appearance of periodic ``ring-like'' structures, may also make this model applicable to other astrophysical systems such as planetary nebulae, e.g. CRL 618 and NGC 6543, young stellar objects, etc.

We discuss a mechanism by which a giant star can expel its envelope in an outburst, leaving its core exposed. The outburst is powered by rotational kinetic energy of the core, transferred to the envelope via the twisting of magnetic fields. We show that, if the core is magnetized, and if it has sufficient angular momentum, this mechanism may be triggered at the end of the asymptotic giant branch phase, and drive a proto-planetary nebula (pPN) outflow. This explosion of magnetic energy self-consistently explains some of the asymmetries and dynamics of pPNe.

We present three-dimensional hydrodynamic simulations of the interaction of a slow wind from an asymptotic giant branch(AGB) star and a jet blown by an orbiting companion. The jet or "Collimated Fast Wind" is assumed to originate from an accretion disk which forms via Bondi accretion of the AGB wind or Roche lobe overflow. We present two distinct regimes in the wind-jet interaction determined by the ratio of the AGB wind to jet momentum flux. Our results show that when the wind momentum flux overwhelms the flux in the jet a more dis-ordered outflow outflow results with the jet assuming a corkscrew pattern and multiple shock structures driven into the AGB wind. In the opposite regime the jet dominates and will drive a highly collimated narrow waisted outflow. We compare our results with scenarios described by Soker & Rappaport (2000) and extrapolate the structures observed in PNe and Symbiotic stars.

We report the discovery, by the Chandra X-ray Observatory, of X-ray emission from the bipolar planetary nebula Menzel 3. In Chandra CCD imaging, Mz 3 displays hot (3-6x10^6 K) gas within its twin, coaxial bubbles of optical nebulosity, as well as a compact X-ray source at the position of its central star(s). The brightest diffuse X-ray emission lies along the polar axis of the optical nebula, suggesting a jet-like configuration. The observed combination of an X-ray-emitting point source and possible X-ray jet(s) is consistent with models in which accretion disks and, potentially, magnetic fields shape bipolar planetary nebulae via the generation of fast, collimated outflows.

We present a design for high energy density laboratory experiments studying the interaction of hypersonic shocks with a large number of inhomogeneities. These ``clumpy'' flows are relevant to a wide variety of astrophysical environments including the evolution of molecular clouds, outflows from young stars, Planetary Nebulae and Active Galactic Nuclei. The experiment consists of a strong shock (driven by a pulsed power machine or a high intensity laser) impinging on a region of randomly placed plastic rods. We discuss the goals of the specific design and how they are met by specific choices of target components. An adaptive mesh refinement hydrodynamic code is used to analyze the design and establish a predictive baseline for the experiments. The simulations confirm the effectiveness of the design in terms of articulating the differences between shocks propagating through smooth and clumpy environments. In particular, we find significant differences between the shock propagation speeds in a clumpy medium compared to a smooth one with the same average density. The simulation results are of general interest for foams in both inertial confinement fusion and laboratory astrophysics studies. Our results highlight the danger of using average properties of inhomogeneous astrophysical environments when comparing timescales for critical processes such as shock crossing and gravitational collapse times.

We describe a new mechanism for pulsations in evolved stars: relaxation oscillations driven by a coupling between the luminosity-dependent mass-loss rate and the H fuel abundance in a nuclear-burning shell. When mass loss is included, the outward flow of matter can modulate the flow of fuel into the shell when the stellar luminosity is close to the Eddington luminosity $L_{\rm Edd}$. When the luminosity drops below $L_{\rm Edd}$, the mass outflow declines and the shell is re-supplied with fuel. This process can be repetitive. We demonstrate the existence of such oscillations and discuss the dependence of the results on the stellar parameters. In particular, we show that the oscillation period scales specifically with the mass of the H-burning relaxation shell (HBRS), defined as the part of the H-burning shell above the minimum radius at which the luminosity from below first exceeds the Eddington threshold at the onset of the mass loss phase. For a stellar mass $M_*\sim 0.7\Msun$, luminosity $L_*\sim 10^4\Lsun$, and mass loss rate $|\dot M|\sim 10^{-5}\Msun$ yr$^{-1}$, the oscillations have a recurrence time $\sim 1400$ years $\sim 57\tau_{\rm fsm}$, where $\tau_{\rm fsm}$ is the timescale for modulation of the fuel supply in the HBRS by the varying mass-loss rate. This period agrees with the $\sim$ 1400-year period inferred for the spacings between the shells surrounding some planetary nebulae, and the the predictied shell thickness, of order 0.4 times the spacing, also agrees reasonably well.

Powerful, highly collimated jets, surrounded by bipolar molecular outflows, are commonly observed near Young Stellar Objects (YSOs). In the usual theoretical picture of star formation, a jet is ejected from a magnetized accretion disk, with a molecular outflow being driven either by the jet or by a wider wind coming from the disk. Here, we propose an alternative global model for the flows surrounding YSOs. In addition to a central accretion-ejection engine driving the jet, the molecular outflow is powered by the infalling matter and follows a circulation pattern around the central object without necessarily being entrained by a jet. It is shown that the model produces a heated pressure-driven outflow with magneto-centrifugal acceleration and collimation. We report solutions for the three different parts of this self-similar model, i.e. the jet, the infalling envelope and the circulating matter that eventually forms the molecular outflow. This new picture of the accretion/outflow phase provides a possible explanation for several observed properties of YSO outflows. The most relevant ones are the presence of high mass molecular outflows around massive protostars, and a realistic fraction (typically 0.1) of the accretion flow that goes into the jet.

We present two-dimensional, cylindrically symmetric simulations of hydrodynamic and magnetohydrodynamic (MHD) wide-angle winds interacting with a collapsing environment. These simulations have direct relevance to young stellar objects (YSOs). The results may also be of use in the study of collimated outflows from proto-planetary and planetary nebulae. We study a range of wind configurations consistent with asymptotic MHD wind collimation. The degree of collimation is parameterized by the ratio of the wind density at the pole to that of the equator. We find that a toroidal magnetic field can have a significant influence on the resulting outflow, giving rise to a very dense, jet-like flow in the post-shock region. The properties of the flow in this region are similar to the asymptotic state of a collimated MHD wind. We conclude that wide-angle MHD winds are quite likely capable of driving molecular outflows. Due to difficulty in treating MHD winds ab-initio in simulations we choose magnetic field strengths in the wind consistent slow magnetic rotators. While MHD launched winds will be in the fast rotator regime we discuss how our results, which rely on toroidal pinch effects, will hold for stronger field strengths.

We consider the stability of an accretion disk wind to cloud formation when subject to a central radiation force. For a vertical launch velocity profile that is Keplerian or flatter and the presence of a significant radiation pressure, the wind flow streamlines cross in a conical layer. We argue that such regions are highly unstable, and are natural sites for supersonic turbulence and, consequently, density compressions. We suggest that combined with thermal instability these will all conspire to produce clouds. Such clouds can exist in dynamical equilibrium, constantly dissipating and reforming. As long as there is an inner truncation radius to the wind, our model emerges with a biconical structure similar to that inferred by Elvis (2000) for the broad line region (BLR) of active galactic nuclei (AGN). Our results may also apply to other disk-wind systems.

Many astrophysical flows occur in inhomogeneous (clumpy) media. We present results of a numerical study of steady, planar shocks interacting with a system of embedded cylindrical clouds. Our study uses a two-dimensional geometry. Our numerical code uses an adaptive mesh refinement allowing us to achieve sufficiently high resolution both at the largest and the smallest scales. We neglect any radiative losses, heat conduction, and gravitational forces. Detailed analysis of the simulations shows that interaction of embedded inhomogeneities with the shock/postshock wind depends primarily on the thickness of the cloud layer and arrangement of the clouds in the layer. The total cloud mass and the total number of individual clouds is not a significant factor. We define two classes of cloud distributions: thin and thick layers. We define the critical cloud separation along the direction of the flow and perpendicular to it distinguishing between the interacting and noninteracting regimes of cloud evolution. Finally we discuss mass-loading and mixing in such systems.

We present first results of astrophysically relevant experiments where highly supersonic plasma jets are generated via conically convergent flows. The convergent flows are created by electrodynamic acceleration of plasma in a conical array of fine metallic wires (a modification of the wire array Z-pinch). Stagnation of plasma flow on the axis of symmetry forms a standing conical shock effectively collimating the flow in the axial direction. This scenario is essentially similar to that discussed by Canto\' ~and collaborators as a purely hydrodynamic mechanism for jet formation in astrophysical systems. Experiments using different materials (Al, Fe and W) show that a highly supersonic ($M\sim 20$), well-collimated jet is generated when the radiative cooling rate of the plasma is significant. We discuss scaling issues for the experiments and their potential use for numerical code verification. The experiments also may allow direct exploration of astrophysically relevant issues such as collimation, stability and jet-cloud interactions.

Many astrophysical flows occur in inhomogeneous media. The broad-line regions (BLR) of active galactic nuclei (AGNs) are one of the important examples where emission-line clouds interact with the outflow. We present results of a numerical study of the interaction of a steady, planar shock / supersonic postshock flow with a system of embedded cylindrical clouds in a two-dimensional geometry. Detailed analysis shows that the interaction of embedded inhomogeneities with the shock / postshock wind depends primarily on the thickness of the cloud layer and the arrangement of the clouds in the layer, as opposed to the total cloud mass and the total number of individual clouds. This allows us to define two classes of cloud distributions: thin and thick layers. We define the critical cloud separation along the direction of the flow and perpendicular to it. This definition allows us to distinguish between the interacting and noninteracting regimes of cloud evolution. Finally we discuss mass-loading in such systems.

We discuss the propagation of radiative MHD jets and outflows focusing on outflows driven by magneto-centrifugal rotators. Our goal is to link the properties of the jets with the physics of the sources which produce them. We find that density and magnetic field stratification (with radius) in jets from magnetized rotators leads to new behavior including the development of a dense inner jet core and a low density collar. We also report on more general studies of ambipolar diffusion and field geometry in pulsed jets. Finally we describe a new work designed to study the effects of magnetized winds on circumstellar environments appropriate to YSOs and PNe.

Recent studies have focused on the the role of initially weak toroidal magnetic fields embedded in a stellar wind as the agent for collimation in planetary nebulae. In these models the wind is assumed to be permeated by a helical magnetic field in which the poloidal component falls off faster than the toroidal component. The collimation only occurs after the wind is shocked at large distances from the stellar source. In this paper we re-examine assumptions built into this ``Magnetized Wind Blown Bubble'' (MWBB) model. We show that a self-consistent study of the model leads to a large parameter regime where the wind is self-collimated before the shock wave is encountered. We also explore the relation between winds in the MWBB model and those which are produced via magneto-centrifugal processes. We conclude that a more detailed examination of the role of self-collimation is needed in the context of PNe studies.

Planetary nebulae are thought to be formed when a slow wind from the progenitor giant star is overtaken by a subsequent fast wind generated as the star enters its white dwarf stage$^{1}$. A shock forms near the boundary between the winds, which creates a relatively dense shell that provides the characteristic appearance of a planetary nebula. A spherically symmetric wind will produce a spherically symmetric shell, yet over half of known planetary nebulae are not spherical; rather, they are elliptical or bipolar in shape$^{2}$. While a magnetic field could launch and collimate a bipolar outflow, the origin of such a field has hitherto been unclear, as previous work suggested that a field could not be generated${^3}$. Here we show that an asymptotic-giant-branch (AGB) star can indeed generate a strong magnetic field, in a dynamo at the interface between a rapidly rotating core and the more slowly rotating envelope of the star. The field is strong enough to shape the bipolar outflows that produce the observed bipolar planetary nebulae. Magnetic braking of the stellar core during this process may also explain the puzzlingly$^{4}$ slow rotation of most white dwarf stars.

MHD winds can emanate from both stars and surrounding accretion disks. It is of interest to know how much wind power is available and which (if either) of the two rotators dominates that power. We investigate this in the context of multi-polar planetary nebulae (PNe) and proto-planetary nebulae (PPNe), for which recent observations have revealed the need for a wind power source in excess of that available from radiation driving, and a possible need for magnetic shaping. We calculate the MHD wind power from a coupled disk and star, where the former results from binary disruption. The resulting wind powers depend only on the accretion rate and stellar properties. We find that if the stellar envelope were initially slowly rotating, the disk wind would dominate throughout the evolution. If the envelope of the star were rapidly rotating, the stellar wind could initially be of comparable power to the disk wind until the stellar wind carries away the star's angular momentum. Since an initially rapidly rotating star can have its spin and magnetic axes misaligned to the disk, multi-polar outflows can result from this disk wind system. For times greater than a spin-down time, the post-AGB stellar wind is slaved to the disk for both slow and rapid initial spin cases and the disk wind luminosity dominates. We find a reasonably large parameter space where a hybrid star+disk MHD driven wind is plausible and where both or either can account for PPNe and PNe powers. We also speculate on the morphologies which may emerge from the coupled system. The coupled winds might help explain the shapes of a number of remarkable multi-shell or multi-polar nebulae. Magnetic activity such as X-ray flares may be associated with the both central star and the disk and would be a valuable diagnostic for the dynamical role of MHD processes in PNe.

We investigate through high resolution 3D simulations the nonlinear evolution of compressible magnetohydrodynamic flows subject to the Kelvin-Helmholtz instability. We confirm in 3D flows the conclusion from our 2D work that even apparently weak magnetic fields embedded in Kelvin-Helmholtz unstable plasma flows can be fundamentally important to nonlinear evolution of the instability. In fact, that statement is strengthened in 3D by this work, because it shows how field line bundles can be stretched and twisted in 3D as the quasi-2D Cat's Eye vortex forms out of the hydrodynamical motions. In our simulations twisting of the field may increase the maximum field strength by more than a factor of two over the 2D effect. If, by these developments, the Alfv\'en Mach number of flows around the Cat's Eye drops to unity or less, our simulations suggest magnetic stresses will eventually destroy the Cat's Eye and cause the plasma flow to self-organize into a relatively smooth and apparently stable flow that retains memory of the original shear. For our flow configurations the regime in 3D for such reorganization is $4\lesssim M_{Ax} \lesssim 50$, expressed in terms of the Alfv\'en Mach number of the original velocity transition and the initial Alfv\'en speed projected to the flow plan. For weaker fields the instability remains essentially hydrodynamic in early stages, and the Cat's Eye is destroyed by the hydrodynamic secondary instabilities of a 3D nature. Then, the flows evolve into chaotic structures that approach decaying isotropic turbulence. In this stage, there is considerable enhancement to the magnetic energy due to stretching, twisting, and turbulent amplification, which is retained long afterwards. The magnetic energy eventually catches up to the kinetic energy, and the nature of flows become magnetohydrodynamic.

Hypersonic, highly collimated, mass outflows ({\it jets}) are a ubiquitous phenomena in astrophysics. While the character of the jets differ, many exhibit some form of quasi-periodic clumping indicating the jet source is episodic or {\it pulsed}. The presence of pulsed jets in so many astrophysical contexts suggests a common formation mechanism. Such a process seems to have been found in {\it Magneto-centrifugal launching}, the combination of magnetic and centrifugal forces that occurs when a magnetized gaseous accretion disk orbits a central gravitating source. Observations of strong magnetic fields in jets are, however, rare or indirect. Thus the presence and effects of magnetic fields in YSO jets remains an unresolved issue of the highest importance. In this letter we focus on what should be expected of the structure of the fields in pulsed YSO jets. We show that combining velocity variability with an initial field configuration consistent with collimated, Magneto-centrifugally launched jets leads to a clear set of predictions concerning the geometry and relative strength of the magnetic field components in evolving YSO (and perhaps other) jets.

We present simulations of the propagation of magnetized jets. This work differs from previous studies in that the cross-sectional distributions of the jets's state variables are derived from analytical models for magneto-centrifugal launching. The source is a magnetized rotator whose properties are specfied as boundary conditions. The jets in these simulations are considerably more complex than the ``top-hat''constant density etc. profiles used in previous work. We find that density and magnetic field stratification (with radius) in the jet leads to new behavior including the separation of an inner jet core from a low density collar. We find this {\it jet within a jet} structure, along with the magnetic stresses, leads to propagation behaviors not observed in previous simulation studies. Our methodology allows us to compare MHD jets from different types of sources whose properties could ultimately be derived from the behavior of the propagating jets.

MHD jet equilibria that depend on source properties are obtained using a simplified model for stationary, axisymmetric and rotating magnetized outflows. The present rotation laws are more complex than previously considered and include a Keplerian disc. The ensuing jets have a dense, current-carrying central core surrounded by an outer collar with a return current. The intermediate part of the jet is almost current-free and is magnetically dominated. Most of the momentum is located around the axis in the dense core and this region is likely to dominate the dynamics of the jet. We address the linear stability and the non-linear development of instabilities for our models using both analytical and 2.5-D numerical simulation's. The instabilities seen in the simulations develop with a wavelength and growth time that are well matched by the stability analysis. The modes explored in this work may provide a natural explanation for knots observed in astrophysical jets.

Theoretical Models for the shaping of PNe are reviewed in light of new high resolution images. The new data indicate the purely hydrodynamic interacting stellar winds model can not recover the full variety of shapes and kinematics. New models, some speculative, others more firmly grounded are discussed. In particular, accretion disks and magnetic fields are identified as two of the most promising avenues of future research. Outstanding issues such as jet formation by PNe disks and dynamo activity in P-AGB stars remain to studied. Finally, new simulations of the Egg Nebula are presented as an example of a ``paleontological'' study designed to recover the history of an individual object.

The interaction of a stellar or disk wind with a collapsing environment holds promise for explaining a variety of outflow phenomena observed around young stars. In this paper we present the first simulations of these interactions. The focus here is on exploring how ram pressure balance between wind and ambient gas and post-shock cooling affects the shape of the resulting outflows. In our models we explore the role of ram pressure and cooling by holding the wind speed constant and adjusting the ratio of the inflow mass flux to the wind mass flux (Mdot_a/Mdot_w) Assuming non-spherical cloud collapse, we find that relatively strong winds can carve out wide, conical outflow cavities and that relatively weak winds can be strongly collimated into jet-like structures. If the winds become weak enough, they can be cut off entirely by the infalling environment. We identify discrepancies between results from standard snowplow models and those presented here that have important implications for molecular outflows. We also present mass vs. velocity curves for comparison with observations.

We have carried out simulations of the nonlinear evolution of the magnetohydrodynamic (MHD) Kelvin-Helmholtz (KH) instability for compressible fluids in $2\frac{1}{2}$-dimensions, extending our previous work by Frank et al (1996) and Jones \etal (1997). In the present work we have simulated flows in the x-y plane in which a ``sheared'' magnetic field of uniform strength ``smoothly'' rotates across a thin velocity shear layer from the z direction to the x direction, aligned with the flow field. We focus on dynamical evolution of fluid features, kinetic energy dissipation, and mixing of the fluid between the two layers, considering their dependence on magnetic field strength for this geometry. The introduction of magnetic shear can allow a Cat's Eye-like vortex to form, even when the field is stronger than the nominal linear instability limit given above. For strong fields that vortex is asymmetric with respect to the preliminary shear layer, however, so the subsequent dissipation is enhanced over the uniform field cases of comparable field strength. In fact, so long as the magnetic field achieves some level of dynamical importance during an eddy turnover time, the asymmetries introduced through the magnetic shear will increase flow complexity, and, with that, dissipation and mixing. The degree of the fluid mixing between the two layers is strongly influenced by the magnetic field strength. Mixing of the fluid is most effective when the vortex is disrupted by magnetic tension during transient reconnection, through local chaotic behavior that follows.

We address the issue of ambipolar diffusion in Herbig-Haro jets. The current consensus holds that these jets are launched and collimated via MHD forces. Observations have, however, shown that the jets can be mildly to weakly ionized. Beginning with a simple model for cylindrical equilibrium between neutral, plasma and magnetic pressures we calculate the characteristic time-scale for ambipolar diffusion. Our results show that a significant fraction of HH jets will have ambipolar diffusion time-scales equivalent to, or less than the dynamical time-scales. This implies that MHD equilibria established at the base of a HH jet may not be maintained as the jet propagates far from its source. For typical jet parameters one finds that the length scale where ambipolar diffusion should become significant corresponds to the typical size of large (parsec) scale jets. We discuss the significance of these results for the issue of magnetic fields in parsec-scale jets.

In this paper we provide a possible link between the structure of the bipolar nebula surrounding SN1987A and the properties of its progenitor star. A Wind Blwon Bubble (WBB) scenario is emplyed, in which a fast, tenuous wind from a Blue Supergiant expands into a slow, dense wind, expelled during an earlier Red Supergiant phase. The bipolar shapre develops due to a pole-to-equator density contrast in the slow wind (ie, the slow wind forms a slow torus). We use the Wind Compressed Disk (WCD) model of Bjorkman & Cassinelli (1992) to determine the shape of the slow torus. In the WCD scenario, the shape of the torus is determined by the rotation of the progenitor star. We then use a self-similar semi-analytical method for wind blown bubble evolution to determine the shape of the resulting bipolar nebula. We find that the union of the wind-compressed-disk and bipolar-wind-blown- bubble models allows us to recover the salient properties of SN1987A's circumstellar nebula. In particular, the size, speed and density of SN1987A's inner ring are easily reproduced in our calculations. An exploration of parameter space shows the the red supergiant progenitor must be been rotating at > 0.3 of its breakup speed. We conclude that the progenitor was most likely spun up by a merger with a binary companion. Using a simple model for the binary merger we find that the companion is likely to have had a mass > 0.5 M_sun.

A description is given for preserving ${\bmsy\nabla}\cdot{\vec B}=0$ in a magnetohydrodynamic (MHD) code that employs the upwind, Total Variation Diminishing (TVD) scheme and the Strang-type operator splitting for multi-dimensionality. The method is based on the staggered mesh technique to constrain the transport of magnetic field: the magnetic field components are defined at grid interfaces with their advective fluxes on grid edges, while other quantities are defined at grid centers. The magnetic field at grid centers for the upwind step is calculated by interpolating the values from grid interfaces. The advective fluxes on grid edges for the magnetic field evolution are calculated from the upwind fluxes at grid interfaces. Then, the magnetic field can be maintained with ${\bmsy\nabla}\cdot{\vec B}=0$ exactly, if this is so initially, while the upwind scheme is used for the update of fluid quantities. The correctness of the code is demonstrated through tests comparing numerical solutions either with analytic solutions or with numerical solutions from the code using an explicit divergence-cleaning method. Also the robustness is shown through tests involving realistic astrophysical problems.

Hypersonic bipolar outflows are a ubiquitous phenomena associated with both young and highly evolved stars. Observations of Planetary Nebulae, the nebulae surrounding Luminous Blue Variables such as $\eta$ Carinae, Wolf Rayet bubbles, the circumstellar environment of SN 1987A and Young Stellar Objects all revealed high velocity outflows with a wide range of shapes. In this paper I review the current state of our theoretical understanding of these outflows. Beginning with Planetary Nebulae considerable progress has been made in understanding bipolar outflows as the result of stellar winds interacting with the circumstellar environment. In what has been called the "Generalized Wind Blown Bubble" (GWBB) scenario, a fast tenuous wind from the central star expands into a ambient medium with an aspherical (toroidal) density distribution. Inertial gradients due to the gaseous torus quickly lead to an expanding prolate or bipolar shell of swept-up gas bounded by strong shock waves. Numerical simulations of the GWBB scenario show a surprisingly rich variety of gasdynamical behavior, allowing models to recover many of the observed properties of stellar bipolar outflows including the development of collimated supersonic jets. In this paper we review the physics behind the GWBB scenario in detail and consider its strengths and weakness. Alternative models involving MHD processes are also examined. Applications of these models to each of the principle classes of stellar bipolar outflow (YSO, PNe, LBV, SN87A) are then reviewed. Outstanding issues in the study of bipolar outflows are considered as are those questions which arise when the outflows are viewed as a single class of phenomena occuring across the HR diagram.

In this paper we address the issue of the origin of LBV bipolar bubbles. Previous studies have explained the shapes of LBV nebulae, such as $\eta$ Car, by invoking the interaction of an isotropic fast wind with a previously deposited, slow aspherical wind (a ``slow torus''). In this paper we focus on the opposite scenario where an aspherical fast wind expands into a previously deposited isotropic slow wind. Using high resolution hydrodynamic simulations, which include the effects of radiative cooling, we have completed a series of numerical experiments to test if and how aspherical fast winds effect wind blown bubble morphologies. Our experiments explore a variety of models for the latitudinal variations of fast wind flow parameters. The simulations demonstrate that aspherical fast winds can produce strongly bipolar outflows. In addition the properties of outflows recover some important aspects of LBV bubbles which the previous "slow torus" models can not.

Collimated outflows from Young Stellar Objects (YSOs) can be seen as tracers of the accretion powered systems which drive them. In this paper I review some theoretical and observational aspects of YSO outflows through the prism of questions relating to the protostellar source. The issue I address is: can collimated outflows be used as ``fossils'' allowing the history of protostellar evolution to be recovered? Answering this question relies on accurately identifying where theoretical tools and observational diagnostics converge to provide unique solutions of the protostellar physics. I discuss potential links between outflow and source including the time and direction variability of jets, the jet/molecular outflow connection, and the the effect of magnetic fields. I also discuss models of the jet/outflow collimation mechanism.

We present multi-dimensional simulations of magnetized radiative jets appropriate to Young Stellar Objects. Magnetized jets subject to collisionally excited radiative losses have not, as yet, received extensive scrutiny. The purpose of this letter is to articulate the propagation dynamics of radiative MHD jets in the context of the extensive jet literature. Most importantly, we look for morphological and kinematic diagnostics that may distinguish hydrodynamic protostellar jets from their magnetically dominated cousins. Our simulations are axisymmetric (2.5-D). A toroidal B_phi field geometry is used. Our models have high sonic Mach numbers (Ms approx 10), but lower fast mode Mach number (Mf approx 5). This is approximately the case for jets formed via disk-wind or X-wind models - currently the consensus choice for launching and collimating YSO jets. Time-dependent radiative losses are included via a coronal cooling curve. Our results demonstrate that the morphology and propagation characteristics of strongly magnetized radiative jets can differ significantly from jets with weak fields. In particular the formation of nose-cones via post-shock hoop stresses leads to narrow bow shocks and enhanced bow shock speeds. In addition, the hoop stresses produce strong shocks in the jet beam which constrasts with the relatively unperturbed beam in radiative hydrodynamic jets. Our simulations show that pinch modes produced by magnetic tension can strongly effect magnetized protostellar jets.

In this paper we explore the effect of radiative losses on purely hydrodynamic jet collimation models applicable to Young Stellar Objects (YSOs). In our models aspherical bubbles form from the interaction of a central YSO wind with an aspherical circum-protostellar density distribution. Building on a previous non-radiative study (Frank & Mellema 1996) we demonstrate that supersonic jets are a natural and robust consequence of aspherical wind-blown bubble evolution. The simulations show that the addition of radiative cooling makes the hydrodynamic collimation mechanisms studied by Frank & Mellema (1996) more effective. An unexpected result is the production of cool, dense jets forming through conical converging flows at the poles of the bubbles. For steady winds the formation of these jets occurs early in the bubble evolution. At later times we find that the dynamical and cooling time scales for the jet material become similar. Since observations of YSO jets show considerable variability in the jet beam we present a simple one-dimensional (1-D) model for the evolution of a variable wind interacting with an accreting environment. These models indicate that, in the presence of a varying protostellar wind, the hydrodynamic collimation processes studied in our simulations can produce cool jets with sizes and time scales consistent with observations.

Previous studies have explained the shapes of LBV nebulae, such as eta Car, by invoking the interaction of an isotropic fast wind with a previously deposited, slow aspherical wind (a "slow torus"). In this work I focus on the opposite scenario where an aspherical fast wind expands into a previously deposited isotropic slow wind. Using high resolution hydrodynamic simulations which include the effects of radiative cooling I have completed a series of numerical experiments to test if and how aspherical fast winds effects wind blown bubble morphologies. The simulations demonstrate that aspherical fast winds can produce strongly bipolar outflows and recover some important aspects

We present the results of numerical hydrodynamic models for the collimation of outflows from young stellar objects. We show that the presence of a toroidal environment can lead to efficient formation of jets and bipolar outflows from initially uncollimated central winds. The interaction between the wind and the environment leads to two types of collimation, one which is dominated by radiative cooling effects, and one which works when cooling is less efficient. We describe the two types of jets as they appear in the simulations and we suggest a description for the long term evolution of these structures in more realistic time-dependent wind sources.

We have carried out high resolution MHD simulations of the nonlinear evolution of Kelvin-Helmholtz unstable flows in 2 1/2 dimensions. The modeled flows and fields were initially uniform except for a thin shear layer with a hyperbolic tangent velocity profile and a small, normal mode perturbation. The calculations consider periodic sections of flows containing magnetic fields parallel to the shear layer, but projecting over a full range of angles with respect to the flow vectors. They are intended as preparation for fully 3D calculations and to address two specific questions raised in earlier work: 1) What role, if any, does the orientation of the field play in nonlinear evolution of the MHD Kelvin-Helmholtz instability in 2 1/2 D. 2) Given that the field is too weak to stabilize against a linear perturbation of the flow, how does the nonlinear evolution of the instability depend on strength of the field. The magnetic field component in the third direction contributes only through minor pressure contributions, so the flows are essentially 2D. Even a very weak field can significantly enhance the rate of energy dissipation. In all of the cases we studied magnetic field amplification by stretching in the vortex is limited by tearing mode, ``fast'' reconnection events that isolate and then destroy magnetic flux islands within the vortex and relax the fields outside the vortex. If the magnetic tension developed prior to reconnection is comparable to Reynolds stresses in the flow, that flow is reorganized during reconnection. Otherwise, the primary influence on the plasma is generation of entropy. The effective expulsion of flux from the vortex is very similar to that shown by Weiss for passive fields in idealized vortices with large magnetic Reynolds numbers. We demonstrated that this

We present 3-D numerical hydrodynamical simulations of precessing supersonic heavy jets to explore how well they serve as a model for generating molecular outflows from Young Stellar Objects. The dynamics are studied with a number of high resolution simulations on a Cartesian grid (128x128x128 zones) using a high order finite difference method. A range of cone angles and precession rates were included in the study. Two higher resolution runs (256x256x256 zones) were made for comparison in order to confirm numerical convergence of global flow characteristics. Morphological, kinematical and dynamical characteristics of precessing jets are described and compared to important properties of straight jets and also to observations of YSOs. In order to examine the robustness of precessing jets as a mean to produce molecular outflows around Young Stellar Objects, ``synthetic observations'' of the momentum distributions of the simulated precessing jets are compared to observations of molecular outflows. It is found that precessing jets match better the morphology, highly forward driven momentum and momentum distributions along the long axis of molecular outflows than do wind-driven or straight jet-driven flow models.

We explore the physics of time-dependent hydrodynamic collimation of jets from Young Stellar Objects (YSOs). Using parameters appropriate to YSOs we have carried out high resolution hydrodynamic simulations modeling the interaction of a central wind with an environment characterized by a moderate opening angle toroidal density distribution. The results show that the the wind/environment interaction produces strongly collimated supersonic jets. The jet is composed of shocked wind gas. Using analytical models of wind blown bubble evolution we show that the scenario studied here should be applicable to YSOs and can, in principle, initiate collimation on the correct scales (R ~ 100 AU). The simulations reveal a number of time-dependent non-linear features not anticipated in previous analytical studies including: a prolate wind shock; a chimney of cold swept-up ambient material dragged into the bubble cavity; a plug of dense material between the jet and bow shocks. We find that the collimation of the jet occurs through both de Laval nozzles and focusing of the wind via the prolate wind shock. Using an analytical model for shock focusing we demonstrate that a prolate wind shock can, by itself, produce highly collimated supersonic jets.

Using a new numerical code we have carried out two-dimensional simulations of the nonlinear evolution of unstable sheared magnetohydrodynamic flows. We considered two cases: a strong magnetic field (Alfven Mach number, M_a = 2.5) and a weak field (M_a =5). Each flow rapidly evolves until it reaches a nearly steady condition, which is fundamentally different from the analogous gasdynamic state. Both MHD flows relax to a stable, laminar flow on timescales less than or of the order of 15 linear growth times, measured from saturation of the instability. That timescale is several orders of magnitude less than the nominal dissipation time for these simulated flows, so this condition represents an quasi-steady relaxed state. The strong magnetic field case reaches saturation as magnetic tension in the displaced flow boundary becomes sufficient to stabilize it. That flow then relaxes in a straightforward way to the steady, laminar flow condition. The weak magnetic field case, on the other hand, begins development of the vortex expected for gasdynamics, but that vortex is destroyed by magnetic stresses that locally become strong. Magnetic topologies lead to reconnection and dynamical alignment between magnetic and velocity fields. Together these processes produce a sequence of intermittent vortices and subsequent relaxation to a nearly laminar flow condition in which the magnetic cross helicity is nearly maximized. Remaining irregularities consist of a pair of flux tubes straddling the shear layer. Fluctuations within those features are closely aligned, representing Alfv\'en waves propagating locally downstream.

We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. The constraint of a divergence-free field is enforced exactly by calculating a correction via a gauge transformation in each time step. Results from two-dimensional shock tube tests show that the code captures correctly discontinuities in all three MHD waves families as well as contact discontinuities. The numerical viscosities and resistivity in the code, which are useful in order to understand simulations involving turbulent flows, are estimated through the decay of two-dimensional linear waves. Finally, the robustness of the code in two-dimensions is demonstrated through calculations of the Kelvin-Helmholtz instability and the Orszag-Tang vortex.

This paper looks into various aspects brought to light by numerical work on the generalized interacting winds model for planetary nebulae. First, a detailed comparison between radiative and non-radiative models is made, showing that one's naive expectations of the effects of radiative heating and cooling are not always true. Secondly, we consider the evolution of the slow wind after it has gotten ionized. It is found that the initial aspherical density distribution is smoothed out after ionization, thus requiring the aspherical nebula to be formed within the first few thousand years of PN evolution. Thirdly, the nature of the flow in the hot bubble is investigated. Both one and two-dimensional numerical models show strong signs of instabilities or turbulent flow in the hot bubble. Although observationally hard to prove or disprove this turbulent structure is critically examined. It is found that although the turbulence is not a numerical artifact, the full three-dimensional picture will most definitely differ from what is found in two dimensions. The implications for the interpretation of the models are discussed. Finally, the issue of the soft X-ray emission from PNe is considered. It is found that soft X-rays originate mainly from the thin interface between the hot bubble and the actual nebula.

The status of numerical hydrodynamical models for Planetary Nebulae is reviewed. Since all of the numerical work is based on the interacting winds model, we start with a description of this model and give an overview of the early analytical and numerical models. Subsequently we address the numerical models which include radiation effects, first of all the ones which neglect any effects of stellar evolution. These `constant environment' models are shown to closely match typical observed nebulae, both in images and kinematic data. This shows that the basic generalized interacting winds model gives a good description of the situation in aspherical PNe. Next we discuss models that do include the effects of stellar and fast wind evolution. This introduces several new effects, the most important of which are the formation of a surrounding attached envelope, and the modification of the expansion of the nebula, which helps in creating aspherical PNe very early on in their evolution. The ionization of the slow wind also leads to a gradual smoothing out of its aspherical character, working against aspherical PNe forming in later stages. Finally we discuss some applications of the model to nebular problems.

We present the results of two-dimensional radiation-gasdynamic simulations of aspherical Planetary Nebulae (PNe) evolution. These simulations were constructed using the Generalized Interacting Stellar Winds (GISW) scenario of Balick (1987) where a fast, tenuous wind from the central star expands into a toroidal, slow, dense wind. We demonstrate that the GISW model can produce a wide range of aspherical flow patterns. We have constructed self-consistent synthetic observations of the models from forbidden line emissivities used in the energy loss term. We present integrated intensity and long-slit spectrum, (Position-Velocity) maps of the models projected at different angles on the sky. These synthetic observations are compared with real intensity and Position-Velocity maps of PNe. We find that there is a very good match between the synthetic and real observations in terms of morphologies, kinematics, and physical conditions. From the results of these simulations we conclude that the GISW scenario can account for most, if not all, PNe morphologies, thus confirming Balick's (1987) conjecture.

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic ray (CR) modified plane shock evolution using the two-fluid model for diffusive particle acceleration. The numerical shocks evolve to published analytical steady state properties. In order to probe the dynamical role of magnetic fields we have explored for these time asymptotic states the parameter space of upstream fast mode Mach number, $M_f$, and plasma $\beta$, compiling the results into maps of dynamical steady state CR acceleration efficiency, $\epsilon_c$. These maps, along with additional numerical experiments, show that $\epsilon_c$ is reduced through the action of compressive work on tangential magnetic fields in CR-MHD shocks. Thus $\epsilon_c$ in low $\beta$, moderate $M_f$ shocks tends to be smaller in quasi perpendicular shocks than it would be high $\beta$ shocks of the same $M_f$. This result supports earlier conclusions that strong, oblique magnetic fields inhibit diffusive shock acceleration. For quasi parallel shocks with $\beta < 1$, on the other hand, $\epsilon_c$ seems to be increased at a given $M_f$ when compared to high $\beta$ shocks. The apparent contradiction to the first conclusion results, however, from the fact that for small $\beta$ quasi parallel shocks, the fast mode Mach number is not a good measure of compression through the shock. That is better reflected in the sonic Mach number, which is greater. Acceleration efficiencies for high and low $\beta$ having comparable sonic Mach numbers are more similar. Time evolution of CR-MHD shocks is qualitatively similar to CR-gasdynamical shocks. However, several potentially interesting differences are apparent.