The buildup of small-scale magnetic helicity which accompanies the oppositely signed growth on large scales is central to conventional dynamical quenching theories of mean-field dynamos. However, the conventional formalism presumes isotropic turbulence and thereby excludes part of the magnetic Lorentz back-reaction. This renders it insufficient to predict the full quenching for general anisotropic flows. To overcome this deficiency, we derive a new generalized quenching formalism that includes the full back-reacting Lorentz force, and a new "selective-damping-$\tau$" closure which conserves magnetic helicity. We apply the formalism to examples of $\bm\alpha^2$ dynamos and show its predicted quenching for different cases of turbulence---isotropic helical, anisotropic helical, and anisotropic non-helical. It predicts stronger-than-conventional quenching in general, but reduces to the conventional case in the helical isotropic limit.

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.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Molecular cloud complexes exhibit both (i) an unfettered Larson-type spectrum over much of their dynamic range, whilst (ii) still producing a much lower star-formation rate than were this cascade to remain unfettered all the way down to star-forming scales. Here we explain the compatibility of these attributes with minimalist considerations of a mass-conserving fragmentation cascade, combined with estimates of stellar feedback. Of importance is that the amount of feedback needed to abate fragmentation and truncate the complex decreases with decreasing scale. The scale at which the feedback momentum matches the free-fall momentum marks a transition scale below most of the cascade is truncated and the molecular cloud complex dissipated. For a $10^6M_\odot$ GMC complex starting with radius of $\sim 50$pc, the combined feedback from young stellar objects, supernovae, radiation, and stellar winds for a GMC cloud complex can truncate the cascade within an outer free-fall time but only after the cascade reaches parsec scales.

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 kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent velocity and the size scale. We derive a dispersion relation for the growth of entropic instabilities in a spherical self-gravitating cloud and find that there is a critical maximum dissipation time scale, about equal to the crossing time, that allows for fragmentation and subsequent star formation. However, the time scale for the loss of turbulent energy may be shorter or longer, for example with rapid thermal cooling or the injection of mechanical energy. Differences in the time scale for energy loss in different star-forming regions may result in differences in the outcome, for example, in the initial mass function.

Constraining dynamo theories of magnetic field origin by observation is indispensable but challenging, in part because the basic quantities measured by observers and predicted by modelers are different. We clarify these differences and sketch out ways to bridge the divide. Based on archival and previously unpublished data, we then compile various important properties of galactic magnetic fields for nearby spiral galaxies. We consistently compute strengths of total, ordered, and regular fields, pitch angles of ordered and regular fields, and we summarize the present knowledge on azimuthal modes, field parities, and the properties of non-axisymmetric spiral features called magnetic arms. We review related aspects of dynamo theory, with a focus on mean-field models and their predictions for large-scale magnetic fields in galactic discs and halos. Further, we measure the velocity dispersion of HI gas in arm and inter-arm regions in three galaxies, M 51, M 74, and NGC 6496, since spiral modulation of the root-mean-square turbulent speed has been proposed as a driver of non-axisymmetry in large-scale dynamos. We find no evidence for such a modulation and place upper limits on its strength, helping to narrow down the list of mechanisms to explain magnetic arms. Successes and remaining challenges of dynamo models with respect to explaining observations are briefly summarized, and possible strategies are suggested. With new instruments like the Square Kilometre Array (SKA), large data sets of magnetic and non-magnetic properties from thousands of galaxies will become available, to be compared with theory.

The association of starspots with magnetic fields leads to an expectation that quantities which correlate with magnetic field strength may also correlate with {starspot} coverage. Since younger stars spin faster and are more magnetically active, assessing whether {starspot} coverage correlates with shorter rotation periods and stellar youth tests these principles. Here we analyze the {starspot} covering fraction versus stellar age for M{-}, G{-}, K{-}, and F{-}type stars based on previously determined variability and rotation periods of over 30,000 {\textit{Kepler}} main-sequence stars. We determine the correlation between age and variability using single and dual power law best fits. We find that {starspot} coverage does indeed decrease with age. Only when the data {are} binned in an effort to remove the effects of activity cycles of individual stars, do statistically significant power law fits emerge for each stellar type. {Using bin averages,} we then find that the {starspot} covering fraction scales with the {X}-ray to bolometric ratio to the power $\lambda$ with {$0.22\pm 0.03 < \lambda < 0.32\pm 0.09$} for {G-type} stars of rotation period below 15 days and for the full range of F{-} and M{-type} stars. For K{-}type stars, we find two branches of $\lambda$ separated by variability bins, with the lower branch showing nearly constant starspot coverage and the upper branch {$\lambda \sim 0.35\pm 0.04$.} G{-}type stars with periods longer than $15$ days exhibit a transition to steeper power law of {$ \lambda \sim 2.4 \pm 1.0$.} The potential connection to previous rotation-age measurements suggesting a magnetic breaking transition at the solar age, corresponding to period of $24.5$ is also of interest.

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 show that non-axisymmetric, non-helical perturbations in an unstratified shear flow produce a shear-plane averaged electromotive force (EMF) proportional to a spatially dependent kinetic helicity. This new "shear-driven $\alpha$-effect" can amplify shear-plane averaged magnetic fields, but is undiscoverable when invoking the homogeneity assumption for dynamo coefficients. It is a candidate mechanism underlying large scale field generation in astrophysical, laboratory, and simulated shear flows, and highlights the essentiality of correctly averaging dynamo coefficients and fields.

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.

Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two questions: (I) what are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) how precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel.Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise.The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work.

The Magnetorotational Instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sources in the universe. But the standard MRI under a purely vertical magnetic field has heretofore never been directly measured despite numerous efforts over more than a decade. Here we report an unambiguous laboratory demonstration of the spring-mass analogue to the standard MRI by comparing motion of a spring-tethered ball within different rotating flows. The experiment corroborates the theory: efficient outward angular momentum transport manifests only for cases with a weak spring in quasi-Keperian flow. Our experimental method accomplishes this in a new way, thereby connecting solid and fluid mechanics to plasma astrophysics.

Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean field dynamo theories used to model the growth of large scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coefficients with successively higher order spatial derivatives of the mean magnetic field. One ingredient of standard models is the kinematic coefficient of the zeroth order term, namely the averaged kinetic pseudotensor $\bm{\alpha}$, that converts toroidal to poloidal fields. Here we demonstrate an efficient way to calculate this quantity for rotating stratified turbulence, whereby the pre-averaged quantity is calculated for the motion of a single plume, and the average is then taken over an ensemble of plumes of different orientations. We calculate the plume dynamics in the most convenient frame, before transforming back to the lab frame and averaging. Our concise configuration space calculation gives essentially identical results to previous lengthier approaches. The present application exemplifies what is a broadly applicable method.

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 extent to which a magnetosphere protects its planetary atmosphere from stellar wind ablation depends upon how well it prevents energy and momentum exchange with the atmosphere and how well it traps otherwise escaping plasma. We focus on the former, and provide a formalism for estimating approximate upper limits on mass, energy, and momentum capture, and use them to constrain loss rates. Our approach quantifies a competition between the local deflection of incoming plasma by a planetary magnetic field and the increase in area for solar wind mass capture provided by this magnetosphere when the wind and planetary fields incur magnetic reconnection. Solar wind capture can be larger for a magnetized planet versus an unmagnetized planet with small ionosphere, even if the rate of energy transfer is less. We find this to be the case throughout Earth's history since the lunar forming impact-- the likely start of the geodynamo. Focusing of incoming charged particles near the magnetic poles can increase the local energy flux for a strong enough planetary field and weak enough wind, and we find that this may be the case for Earth's future. The future protective role of Earth's magnetic field would then depend on its ability to trap outgoing plasma. This competition between increased collection area and reduced inflow speed from a magnetosphere is likely essential in determining the net protective role of planetary fields and its importance in sustaining habitability, and can help explain the current ion loss rates from Mars, Venus and Earth.

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.

The initial magnetic field of previous magnetorotational instability (MRI) simulations has always included a significant system-scale component, even if stochastic. However, it is of conceptual and practical interest to assess whether the MRI can grow when the initial field is turbulent. The ubiquitous presence of turbulent or random flows in astrophysical plasmas generically leads to a small-scale dynamo (SSD), which would provide initial seed turbulent velocity and magnetic fields in the plasma that becomes an accretion disc. Can the MRI grow from these more realistic initial conditions? To address this we supply a standard shearing box with isotropically forced SSD generated magnetic and velocity fields as initial conditions, and remove the forcing. We find that if the initially supplied fields are too weak or too incoherent, they decay from the initial turbulent cascade faster than they can grow via the MRI. When the initially supplied fields are sufficient to allow MRI growth and sustenance, the saturated stresses, large-scale fields, and power spectra match those of the standard zero net flux MRI simulation with an initial large scale vertical field.

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.

Accretion disc theory is less developed than stellar evolution theory although a similarly mature phenomenological picture is ultimately desired. While the interplay of theory and numerical simulations has amplified community awareness of the role of magnetic fields in angular momentum transport, there remains a long term challenge to incorporate insight gained from simulations back into improving practical models for comparison with observations. Here we emphasize the need to incorporate the role of non-local transport more precisely. To show where large scale transport would fit into the theoretical framework and how it is currently missing, we review why the wonderfully practical approach of Shakura-Sunyaev (1973,SS73) is necessarily a mean field theory, and one which does not include large scale transport. Observations of coronae and jets combined with the interpretation of results even from shearing box simulations of the magnetorotational instability (MRI) suggest that a significant fraction of disc transport is indeed non-local. We show that the Maxwell stresses in saturation are dominated by large scale contributions and the physics of MRI transport is not fully captured by a viscosity. We also clarify the standard physical interpretation of the MRi as it applies to shearing boxes. Computational limitations have so far focused most attention toward local simulations but the next generation of global simulations should help to inform improved mean field theories. Mean field accretion theory and mean field dynamo theory should in fact be unified into a single theory that predicts the time evolution of spectra and luminosity from separate disc, corona, and outflow contributions. Finally, we note that any mean field theory has a finite predictive precision that needs to be quantified when comparing the predictions to observations.

Late-type main sequence stars exhibit an x-ray to bolometric flux ratio that depends on ${\tilde Ro}$, the ratio of rotation period to convective turnover time, as ${\tilde Ro}^{-\zeta}$ with $2\le \zeta \le 3$ for ${\tilde Ro} >0.13$, but saturates with $|\zeta| <0.2$ for ${\tilde Ro} < 0.13$. Saturated stars are younger than unsaturated stars and show a broader spread of rotation rates and x-ray activity. The unsaturated stars have magnetic fields and rotation speeds that scale roughly with the square root of their age, though possibly flattening for stars older than the sun. The connection between faster rotators, stronger fields, and higher activity has been established observationally, but a theory for the unified time-evolution of x-ray luminosity, rotation, magnetic field and mass loss that captures the above trends has been lacking. Here we derive a minimalist holistic framework for the time evolution of these quantities built from combining a Parker wind with new ingredients: (1) explicit sourcing of both the thermal energy launching the wind and the x-ray luminosity via dynamo produced magnetic fields; (2) explicit coupling of x-ray activity and mass loss saturation to dynamo saturation (via magnetic helicity build-up and convection eddy shredding); (3) use of coronal equilibrium to determine how magnetic energy is divided into wind and x-ray contributions. For solar-type stars younger than the sun, we infer conduction to be a subdominant power loss compared to x-rays and wind. For older stars, conduction is more important, possibly quenching the wind and reducing angular momentum loss. We focus on the time evolution for stars younger than the sun, highlighting what is possible for further generalizations. Overall, the approach shows promise toward a unified explanation of all of the aforementioned observational trends.

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.

Galactic dynamo models sustained by supernova (SN) driven turbulence and differential rotation have revealed that the sustenance of large scale fields requires a flux of small scale magnetic helicity to be viable. Here we generalize a minimalist analytic version of such galactic dynamos to explore some heretofore unincluded contributions from shear on the total turbulent energy and turbulent correlation time, with the helicity fluxes maintained by either winds, diffusion, or magnetic buoyancy. We construct an analytic framework for modeling the turbulent energy and correlation time as functions of SN rate and shear. We compare our prescription with previous approaches that only include rotation. The solutions depend separately on the rotation period and the eddy turnover time and not just on their ratio (the Rossby number). We consider models in which these two time scales are allowed to be independent and also a case in which they are mutually dependent on radius when a radial dependent SN rate model is invoked. For the case of a fixed rotation period (or fixed radius) we show that the influence of shear is dramatic for low Rossby numbers, reducing the correlation time of the turbulence, which in turn, strongly reduces the saturation value of the dynamo compared to the case when the shear is ignored. We also show that even in the absence of winds or diffusive fluxes, magnetic buoyancy may be able to sustain sufficient helicity fluxes to avoid quenching.

From numerical simulations, we show that non-rotating magnetohydrodynamic shear flows are unstable to finite amplitude velocity perturbations and become turbulent, leading to the growth and sustenance of magnetic energy, including large scale fields. This supports the concept that sustained magnetic energy from turbulence is independent of the driving mechanism for large enough magnetic Reynolds numbers.

We revisit the nature of large dips in flux from extinction by dusty circumstellar material that is observed by Kepler for many young stars in the Upper Sco and $\rho$ Oph star formation regions. These young, low-mass "dipper" stars are known to have low accretion rates and primarily host moderately evolved dusty circumstellar disks. Young low mass stars often exhibit rotating star spots that cause quasi-periodic photometric variations. We found no evidence for periods associated with the dips that are different from the star spot rotation period in spectrograms constructed from the light curves. The material causing the dips in most of these light curves must be approximately corotating with the star. We find that disk temperatures computed at the disk corotation radius are cool enough that dust should not sublime. Crude estimates for stellar magnetic field strengths and accretion rates are consistent with magnetospheric truncation near the corotation radius. Magnetospheric truncation models can explain why the dips are associated with material near corotation and how dusty material is lifted out of the midplane to obscure the star which would account for the large fraction of young low mass stars that are dippers. We propose that variations in disk orientation angle, stellar magnetic field dipole tilt axis, and disk accretion rate are underlying parameters accounting for differences in the dipper light curves.

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.

We study the stability properties of Rayleigh unstable flows both in the purely hydrodynamic and magnetohydrodynamic (MHD) regimes for two different values of the shear $q=2.1, 4.2$ ($q = - d\ln\Omega / d\ln r$) and compare it with the Keplerian case $q=1.5$. We find that the $q>2$ regime is unstable both in the hydrodynamic and in the MHD limit (with an initially weak magnetic field). In this regime, the velocity fluctuations dominate the magnetic fluctuations. In contrast, in the $q<2$ (magnetorotational instability (MRI)) regime the magnetic fluctuations dominate. This highlights two different paths to MHD turbulence implied by the two regimes, suggesting that in the $q>2$ regime the instability produces primarily velocity fluctuations that cause magnetic fluctuations, with the causality reversed for the $q<2$ MRI unstable regime. We also find that the magnetic field correlation is increasingly localized as the shear is increased in the Rayleigh unstable regime. In calculating the time evolution of spatial averages of different terms in the MHD equations, we find that the $q>2$ regime is dominated by terms which are nonlinear in the fluctuations, whereas for $q<2$, the linear terms play a more significant role.

One in 5 planetary nebulae are ejected from common envelope binary interactions but Kepler Space Telescope results are already showing this proportion to be larger. Their properties, such as abundances can be starkly different from those of the general population, so they should be considered separately when using PN as chemical or population probes. Unfortunately post-common envelope PN cannot be discerned using only their morphologies, but this will change once we couple our new common envelope simulations with PN formation models.

Approximate empirical models of thermal convection can allow us to identify the essential properties of the flow in simplified form, and to produce empirical estimates using only a few parameters. Such "low-dimensional" empirical models can be constructed systematically by writing numerical or experimental measurements as superpositions of a set of appropriate basis modes, a process known as Galerkin projection. For three-dimensional convection in a cylinder, those basis modes should be vector-valued, mutually orthogonal, and defined in cylindrical coordinates. Here we construct such a basis set and demonstrate that it has these desired properties and boundary conditions when the exact constraint of incompressibility is relaxed. We show its use for representing sample simulation data and point out its potential for low-dimensional convection models.

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.

We study the dynamo generation (exponential growth) of large scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal ($x$-$y$) averaging, we demonstrate the presence of large scale fields when either horizontal or vertical ($y$-$z$) averaging is employed. By computing planar averaged fields and power spectra, we find large scale dynamo action in the early MRI growth phase---a previously unidentified feature. Fast growing horizontal low modes and fiducial vertical modes over a narrow range of wave numbers amplify these planar averaged fields in the MRI growth phase, before turbulence sets in. The large scale field growth requires linear fluctuations but not nonlinear turbulence (as defined by mode-mode coupling) and grows as a direct global mode of the MRI. Only by vertical averaging, can it be shown that the growth of horizontal low wavenumber MRI modes directly feed-back to the initial vertical field providing a clue as to why the large scale vertical field sustains against turbulent diffusion in the saturation regime. We compute the terms in the planar averaged mean field equations to identify the individual contributions to large scale field growth for both vertical and horizontal averaging. The large scale fields obtained from such vertical averaging are found to compare well with global cylindrical simulations and quasilinear analytical analysis from a previous study by Ebrahimi \& Blackman. We discuss the potential implications of these new results for understanding large scale MRI dynamo saturation and turbulence.

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.

For cylindrical differentially rotating plasmas, we study large-scale magnetic field generation from finite amplitude non-axisymmetric perturbations by comparing numerical simulations with quasi-linear analytic theory. When initiated with a vertical magnetic field of either zero or finite net flux, our global cylindrical simulations exhibit the magnetorotational instability (MRI) and large scale dynamo growth of radially alternating mean fields, averaged over height and azimuth. This dynamo growth is explained by our analytic calculations of a non-axisymmetric fluctuation-induced EMF that is sustained by azimuthal shear of the fluctuating fields. The standard "Omega effect" (shear of the mean field by differential rotation) is unimportant. For the MRI case, we express the large-scale dynamo field as a function of differential rotation. The resulting radially alternating large-scale fields may have implications for angular momentum transport in disks and corona. To connect with previous work on large scale dynamos with local linear shear and identify the minimum conditions needed for large scale field growth, we also solve our equations in local Cartesian coordinates. We find that large scale dynamo growth in a linear shear flow without rotation can be sustained by shear plus non-axisymmetric fluctuations -- even if not helical, a seemingly previously unidentified distinction. The linear shear flow dynamo emerges as a more restricted version of our more general new global cylindrical calculations.

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.

The extent to which large scale magnetic fields are susceptible to turbulent diffusion is important for interpreting the need for in situ large scale dynamos in astrophysics and for observationally inferring field strengths compared to kinetic energy. By solving coupled equations for magnetic energy and magnetic helicity in a system initiated with isotropic turbulence and an arbitrarily helical large scale field, we quantify the decay rate of the latter for a bounded or periodic system. The energy associated with the non-helical magnetic field rapidly decays by turbulent diffusion, but the decay rate of the helical component depends on whether the ratio of its magnetic energy to the turbulent kinetic energy exceeds a critical value given by M_{1,c} =(k_1/k_2)^2, where k_1 and k_2 are the wave numbers of the large and forcing scales. Turbulently diffusing helical fields to small scales while conserving magnetic helicity requires a rapid increase in total magnetic energy. As such, only when the helical fields are sub-critical can they so diffuse. When super-critical, the large scale helical field decays slowly, at a rate determined by microphysical dissipation even when macroscopic turbulence is present. Amplification of small scale magnetic helicity abates the turbulent diffusion. Two implications are that: (1) Standard arguments supporting the need for in situ large scale dynamos based on the otherwise rapid turbulent diffusion of large scale fields require re-thinking since only the non-helical field is so diffused in a closed system. Boundary terms could however provide potential pathways for rapid change of the large scale helical field. (2) Since M_{1,c} <<1 for k_1 << k_2, the presence of long-lived ordered large scale helical fields, as in extragalactic jets, does not guarantee that the magnetic field dominates the kinetic energy.

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.

The onset and nature of the earliest geomagnetic field is important for understanding the evolution of the core, atmosphere and life on Earth. A record of the early geodynamo is preserved in ancient silicate crystals containing minute magnetic inclusions. These data indicate the presence of a geodynamo during the Paleoarchean, between 3.4 and 3.45 billion years ago. While the magnetic field sheltered Earth's atmosphere from erosion at this time, standoff of the solar wind was greatly reduced, and similar to that during modern extreme solar storms. These conditions suggest that intense radiation from the young Sun may have modified the atmosphere of the young Earth by promoting loss of volatiles, including water. Such effects would have been more pronounced if the field were absent or very weak prior to 3.45 billion years ago, as suggested by some models of lower mantle evolution. The frontier is thus trying to obtain geomagnetic field records that are >>3.45 billion-years-old, as well as constraining solar wind pressure for these times. In this review we suggest pathways for constraining these parameters and the attendant history of Earth's deep interior, hydrosphere and atmosphere. In particular, we discuss new estimates for solar wind pressure for the first 700 million years of Earth history, the competing effects of magnetic shielding versus solar ion collection, and bounds on the detection level of a geodynamo imposed by the presence of external fields. We also discuss the prospects for constraining Hadean-Paleoarchean magnetic field strength using paleointensity analyses of zircons.

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.

The magnetorotational instability (MRI) is a shear instability and thus its sensitivity to the shear parameter $q = - d\ln\Omega/d\ln r $ is of interest to investigate. Motivated by astrophysical disks, most (but not all) previous MRI studies have focused on the Keplerian value of $ q=1.5$. Using simulation with 8 vertical density scale heights, we contribute to the subset of studies addressing the the effect of varying $q$ in stratified numerical simulations. We discuss why shearing boxes cannot easily be used to study $q>2$ and thus focus on $q<2$. As per previous simulations, which were either unstratified or stratified with a smaller vertical domain, we find that the $q$ dependence of stress for the stratified case is not linear, contrary to the Shakura-Sunyaev model. We find that the scaling agrees with \cite{1996MNRAS.281L..21A} who found it to be proportional to the shear to vorticity ratio $q/(2-q)$. We also find however, that the shape of the magnetic and kinetic energy spectra are relatively insensitive to $q$ and that the ratio of Maxwell stress to magnetic energy ratio also remains nearly independent of $q$. This is consistent with a theoretical argument in which the rate of amplification of the azimuthal field depends linearly on $q$ and the turbulent correlation time $\tau$ depends inversely on $q$. As such, we measure the correlation time of the turbulence and find that indeed it is inversely proportional to $q$.

Observations of late-type main-sequence stars have revealed empirical scalings of coronal activity versus rotation period or Rossby number $Ro$ (a ratio of rotation period to convective turnover time) which has hitherto lacked explanation. For $Ro >> 1$, the activity observed as X-ray to bolometric flux varies as $Ro^{-q}$ with $2\le q \le 3$, whilst $|q| < 0.12$ for $Ro << 1$. Here we explain the transition between these two regimes and the power law in the $Ro >> 1$ regime by constructing an expression for the coronal luminosity based on dynamo magnetic field generation and magnetic buoyancy. We explain the $Ro<<1$ behavior from the inference that observed rotation is correlated with internal differential rotation and argue that once the shear time scale is shorter than the convective turnover time, eddies will be shredded on the shear time scale and so the eddy correlation time actually becomes the shear time and the convection time drops out of the equations. We explain the $Ro >> 1$ behavior using a dynamo saturation theory based on magnetic helicity buildup and buoyant loss.

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.

Blackman & Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here we present such a schematic ribbon diagram for the $\alpha^2$ dynamo that tracks magnetic helicity and provides distinct scales of large scale magnetic helicity, small scale magnetic helicity, and kinetic helicity involved in the process. This also motivates our construction of a new "2.5 scale" minimalist generalization of the helicity-evolving equations for the \alpha^2 dynamo that separately allows for these three distinct length scales while keeping only two dynamical equations. We solve these equations and, as in previous studies, find that the large scale field first grows at a rate independent of the magnetic Reynolds number R_M before quenching to an R_M dependent regime. But we also show that the larger the ratio of the wavenumber where the small scale current helicity resides to that of the forcing scale, the earlier the non-linear dynamo quenching occurs, and the weaker the large scale field is at the turnoff from linear growth. The harmony between the theory and the schematic diagram exemplifies a general lesson that magnetic fields in MHD are better visualized as two-dimensional ribbons (or pairs of lines) rather than single lines.

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.

The extent to which angular momentum transport in accretion discs is primarily local or non-local and what determines this is an important avenue of study for understanding accretion engines. Taking a step along this path, we analyze simulations of the magnetorotational instability (MRI) by calculating energy and stress power spectra in stratified isothermal shearing box simulations in several new ways. We divide our boxes in two regions, disc and corona where the disc is the MRI unstable region and corona is the magnetically dominated region. We calculate the fractional power in different quantities, including magnetic energy and Maxwell stresses and find that they are dominated by contributions from the lowest wave numbers. This is even more dramatic for the corona than the disc, suggesting that transport in the corona region is dominated by larger structures than the disc. By calculating averaged power spectra in one direction of $k$ space at a time, we also show that the MRI turbulence is strongly anisotropic on large scales when analyzed by this method, but isotropic on small scales. Although the shearing box itself is meant to represent a local section of an accretion disc, the fact that the stress and energy are dominated by the largest scales highlights that the locality is not captured within the box. This helps to quantify the intuitive importance of global simulations for addressing the question of locality of transport, for which similar analyses can be performed.

Interaction between a central outflow and a surrounding wind is common in astrophysical sources powered by accretion. Understanding how the interaction might help to collimate the inner central outflow is of interest for assessing astrophysical jet formation paradigms. In this context, we studied the interaction between two nested supersonic plasma flows generated by focusing a long pulse high-energy laser beam onto a solid target. A nested geometry was created by shaping the energy distribution at the focal spot with a dedicated phase plate. Optical and X-ray diagnostics were used to study the interacting flows. Experimental results and numerical hydrodynamic simulations indeed show the formation of strongly collimated jets. Our work experimentally confirms the "shock-focused inertial confinement" mechanism proposed in previous theoretical astrophysics investigations.

Magnetic fields of laboratory, planetary, stellar, and galactic plasmas commonly exhibit significant order on large temporal or spatial scales compared to the otherwise random motions within the hosting system. Such ordered fields can be measured in the case of planets, stars, and galaxies, or inferred indirectly by the action of their dynamical influence, such as jets. Whether large scale fields are amplified in situ or a remnant from previous stages of an object's history is often debated for objects without a definitive magnetic activity cycle. Magnetic helicity, a measure of twist and linkage of magnetic field lines, is a unifying tool for understanding large scale field evolution for both mechanisms of origin. Its importance stems from its two basic properties: (1) magnetic helicity is typically better conserved than magnetic energy; and (2) the magnetic energy associated with a fixed amount of magnetic helicity is minimized when the system relaxes this helical structure to the largest scale available. Here I discuss how magnetic helicity has come to help us understand the saturation of and sustenance of large scale dynamos, the need for either local or global helicity fluxes to avoid dynamo quenching, and the associated observational consequences. I also discuss how magnetic helicity acts as a hindrance to turbulent diffusion of large scale fields, and thus a helper for fossil remnant large scale field origin models in some contexts. I briefly discuss the connection between large scale fields and accretion disk theory as well. The goal here is to provide a conceptual primer to help the reader efficiently penetrate the literature.

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

Some combination of binary interactions and accretion plausibly conspire to produce the ubiquitous collimated outflows from planetary nebulae (PN) and their presumed pre-planetary nebulae (PPN) precursors. But which accretion engines are viable? The difficulty in observationally resolving the engines warrants the pursuit of indirect constraints. We show how kinematic outflow data for 19 PPN can be used to determine the minimum required accretion rates. We consider main sequence (MS) and white dwarf (WD) accretors and five example accretion rates inferred from published models to compare with the minima derived from outflow momentum conservation. While our primary goal is to show the method in anticipation of more data and better theoretical constraints, taking the present results at face value already rule out modes of accretion: Bondi-Hoyle Lyttleton (BHL) wind accretion and wind Roche lobe overflow (M-WRLOF, based on Mira parameters) are too feeble for all 19/19 objects for a MS accretor. For a WD accretor, BHL is ruled out for 18/19 objects and M-WRLOF for 15/19 objects. Roche lobe overflow (RLOF) from the primary at the Red Rectangle level can accommodate 7/19 objects, though RLOF modes with higher accretion rates are not yet ruled out. Accretion modes operating from within common envelope evolution can accommodate all 19 objects, if jet collimation can be maintained. Overall, sub-Eddington rates for a MS accretor are acceptable but 8/19 would require super-Eddington rates for a WD.

The recent study of Blackman and Subramanian (Paper I) indicates that large scale helical magnetic fields are resilient to turbulent diffusion in the sense that helical fields stronger than a critical value, decay on slow (~resistively mediated), rather than fast ($\sim$ turbulent) time scales. This gives more credence to potential fossil field origin models of the magnetic fields in stars, galaxies and compact objects. We analyze a suite of direct numerical simulations (DNS) of decaying large scale helical magnetic fields in the presence of non-helical turbulence to further study the physics of helical field decay. We study two separate cases: (1) the initial field is large enough to decay resistively, is tracked until it transitions to decay fast, and the critical large scale helical field at that transition is sought; (2) the case of Paper I, wherein there is a critical initial helical field strength below which the field undergoes fast decay right from the beginning. In case (1), both DNS and solutions of the two scale model (from Paper 1), reveal that the transition energy, $E_{c1}$, is independent of the scale of the turbulent forcing, within a small range of $\Rm$. We also find that the kinetic alpha, $\alpha_K$, is subdominant to magnetic alpha, $\alpha_M$, in the DNS, justifying an assumption in the two scale model. For case (2), we show exact solutions of two scale model in the limit of $\eta \rightarrow 0$ in fully helical case, leading to the transition energy, $E_{c2} = (k_1/\kf)^2 M_{eq}$, where $k_1$ and $\kf$ are the large scale and small turbulent forcing scale respectively and $M_{eq}$ is the equipartition magnetic energy. The DNS in this case agree qualitatively with the two scale model but the $R_M$ currently achievable, is too small to satisfy a condition $3/R_M << (k_1/k_f)^2$, necessary to robustly reveal the transition, $E_{c2}$ (Abridged).

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.

Understanding large scale magnetic field growth in turbulent plasmas in the magnetohydrodynamic limit is a goal of magnetic dynamo theory. In particular, assessing how well large scale helical field growth and saturation in simulations matches that predicted by existing theories is important for progress. Using numerical simulations of isotropically forced turbulence without large scale shear with the implications, we focus on several aspects of this comparison that have not been previously tested: (1) Leading mean field dynamo theories which break the field into large and small scales predict that large scale helical field growth rates are determined by the difference between kinetic helicity and current helicity with no dependence on the non-helical energy in small scale magnetic fields. Our simulations show that the growth rate of the large scale field from fully helical forcing is indeed unaffected by the presence or absence of small scale magnetic fields amplified in a precursor non-helical dynamo. However, because the precursor non helical dynamo in our simulations produced fields that were strongly sub-equipartition with respect to the kinetic energy, we cannot yet rule out the potential influence of stronger non- helical small scale fields. (2) We have identified two features in our simulations which cannot be explained by the most minimalist versions of two-scale mean field theory: (i) fully helical small scale forcing produces significant non-helical large scale magnetic energy and (ii) the saturation of the large scale field growth is time-delayed with respect to what minimalist theory predicts. We comment on desirable generalizations to the theory in this context and future desired work.

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.

The magnetic flux freezing theorem or Alfv\'en's theorem, is a basic principle of ideal magnetohydrodynamics (MHD), a commonly used approximation to describe the aspects of astrophysical and laboratory plasmas. The theorem states that the magnetic flux---the integral of magnetic field penetrating a surface---is conserved in time as that surface is distorted in in time by fluid motions. Pedagogues of MHD commonly derive flux freezing without showing how to take the material derivative of a general flux integral and/or assuming a vanishing field divergence from the outset. Here I avoid these shortcomings and derive flux freezing by direct differentiation, explicitly using a Jacobian to transform between the evolving field-penetrating surface at different times. The approach is instructive for its generality and helps elucidate the role of magnetic monopoles in breaking flux freezing. The paucity of appearances of this derivation in standard MHD texts suggests that its pedagogic value is underappreciated.

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.

Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new 2D particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simulations show that electrons above a certain energy threshold experience shock-drift-acceleration. The transition energy between the thermal and non-thermal spectrum and the spectral index from the simulations are consistent with some of the X-ray spectra from RHESSI in the energy regime of $E\lesssim 40\sim 100$ keV. Plasma instabilities associated with the shock structure such as the modified-two-stream and the electron whistler instabilities are identified using numerical solutions of the kinetic dispersion relations. We also show that the results from PIC simulations with reduced ion/electron mass ratio can be scaled to those with the realistic mass ratio.

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.

Conversion of gravitational energy into radiation in accretion discs and the origin of large scale magnetic fields in astrophysical rotators have often been distinct topics of research. In semi-analytic work on both problems it has been useful to presume large scale symmetries, necessarily resulting in mean field theories. MHD turbulence makes the underlying systems locally asymmetric and nonlinear. Synergy between theory and simulations should aim for the development of practical mean field models that capture essential physics and can be used for observational modeling. Mean field dynamo (MFD) theory and alpha-viscosity accretion theory exemplify such ongoing pursuits. 21st century MFD theory has more nonlinear predictive power compared to 20th century MFD theory, whereas accretion theory is still in a 20th century state. In fact, insights from MFD theory are applicable to accretion theory and the two are artificially separated pieces of what should be a single theory. I discuss pieces of progress that provide clues toward a unified theory. A key concept is that large scale magnetic fields can be sustained via local or global magnetic helicity fluxes or via relaxation of small scale magnetic fluctuations, without the kinetic helicity driver of 20th century textbooks. These concepts may help explain the formation of large scale fields that supply non-local angular momentum transport via coronae and jets in a unified theory of accretion and dynamos. In diagnosing the role of helicities and helicity fluxes in disk simulations, each disk hemisphere should be studied separately to avoid being misled by cancelation that occurs as a result of reflection asymmetry. The fraction of helical field energy in disks is expected to be small compared to the total field in each hemisphere as a result of shear, but can still be essential for large scale dynamo action.

Progress toward understanding principles of nonlinear growth and saturation of large scale magnetic fields has emerged from comparison of theoretical models that incorporate the evolution of magnetic helicity with numerical simulations for problems that are more idealized than expected in astrophysical circumstances, but still fully non-linear. We carry out a new comparison of this sort for the magnetic field growth from forced isotropic helical turbulence in a periodic box. Previous comparisons be- tween analytic theory and simulations of this problem have shown that a two-scale model compares well with the simulations in agreeing that the driver of large scale field growth is the difference between kinetic and current helicities associated with the small scale field, and that the backreaction that slows the growth of the large scale field as the small scale current helicity grows. However, the two-scale model artificially re- stricts the small scale current helicity to reside at the same scale as the driving kinetic helicity. In addition, previous comparisons have focused on the late time saturation regime, and less on the early-time growth regime. Here we address these issues by comparing a three scale model to new simulations for both early and late time growth regimes. We find that the minimalist extension to three scales provides a better model for the field evolution than the two-scale model because the simulations show that the small scale current helicity does not reside at the same scale as that of the driving kinetic helicity. The simulations also show that the peak of the small scale current helicity migrates toward lower wave numbers as the growth evolves from the fast to saturated growth regimes.

Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f_h,C has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f_h,C decreases as the ratio of forcing to large scale wave numbers k_F/k_min increases. From the condition that a large scale helical dynamo must overcome the backreaction from any non-helical field on the large scales, we develop a theory that can explain the simulations. For k_F/k_min>8 we find f_h,C< 3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate scale separation.

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.

Astrophysical shocks are often studied in the high Mach number limit but weakly compressive fast shocks can occur in magnetic reconnection outflows and are considered to be a site of particle energization in solar flares. Here we study the microphysics of such perpendicular, low Mach number collisionless shocks using two-dimensional particle-in-cell (PIC) simulations with a reduced ion/electron mass ratio and employ a moving wall boundary method for initial generation of the shock. This moving wall method allows for more control of the shock speed, smaller simulation box sizes, and longer simulation times than the commonly used fixed wall, reflection method of shock formation. Our results, which are independent of the shock formation method, reveal the prevalence shock drift acceleration (SDA) of both electron and ions in a purely perpendicular shock with Alfv\'en Mach number $M_A=6.8$ and ratio of thermal to magnetic pressure $\beta=8$. We determine the respective minimum energies required for electrons and ions to incur SDA. We derive a theoretical electron distribution via SDA that compares to the simulation results. We also show that a modified two-stream instability due to the incoming and reflecting ions in the shock transition region acts as the mechanism to generate collisionless plasma turbulence that sustains the shock.

Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called "{\alpha}2" LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations. When the system is helically kinetically forced (KF), the growth of the large scale helical field is accompanied by growth of small scale magnetic (and current) helicity which ultimately quench the LSD. Here, using both simulations and theory, we study the complementary magnetically forced(MF) case in which the system is forced with an electric field that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the back-reactor that saturates the LSD. Simulations of both MF and KF cases can be approximately modeled with the same equations of magnetic helicity evolution, but with complementary initial conditions. A key difference between KF and MF cases is that the helical large scale field in the MF case grows with the same sign of injected magnetic helicity, whereas the large and small scale magnetic helicities grow with opposite sign for the KF case. The MF case can arise even when the thermal pressure is approximately smaller than the magnetic pressure, and requires only that helical small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are actually MF LSDs.

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.

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.

Magnetic reconnection is commonly perceived to drive flow and particle acceleration in flares of solar, stellar, and astrophysical disk coronae but the relative roles of different acceleration mecha- nisms in a given reconnection environment are not well understood. We show via direct numerical simulations that reconnection outflows produce weak fast shocks, when conditions for fast recon- nection are met and the outflows encounter an obstacle. The associated compression ratios lead to a Fermi acceleration particle spectrum that is significantly steeper than the strong fast shocks commonly studied, but consistent with the demands of solar flares. While this is not the only acceleration mechanism operating in a reconnection environment, it is plausibly a ubiquitous one.

The recent paper by Heyer et al (2009) indicates that observations of size, linewidth and column density of interstellar clouds do not agree with simple virial equilibrium (VE) as a balance between gravitational and kinetic energies in the sense that the clouds either have too much kinetic energy or too little mass to be bound. This may be explained by violation of VE as suggested by Dobbs et al 2011, by observational underestimation of the masses as suggested by Heyer et al 2009, or by an external pressure acting as an additional confining force as suggested earlier by Heyer et al 2004. The data of Heyer et al. 2009 cannot be explained with a single value for the external pressure, but if different clouds in the sample have different external pressures in the range of Pe/k = E4 to E7 cm-3 K, then most of the clouds could be in pressure virial equilibrium (PVE). In this paper we discuss two consequences of the external pressure. First, we show that the observational data are consistent with the hypothesis (Chie\'ze 1987) that most clouds are at a critical mass for dynamical stability determined solely by the pressure. Above this mass a cloud is unstable to gravitational collapse or fragmentation. Second, we show that the external pressure modifies the well-known size-linewidth relationship first proposed by Larson (1981) so that the proportionality is no longer constant but depends on the external pressure.

In their engaging recountals of Chandrasekhar's extraordinary career (Physics Today, vol 63, Issue 12, Dec 2010), neither Dyson nor Wali mention that Chandrasekhar was the third person not the first, to publish a white dwarf mass limit incorporating a relativistic treatment of degenerate electrons. As it has become a common misconception that Chandrasekhar was the first, a clarifying reminder on this historical point is warranted. In short, the white dwarf mass limit widely attributed to Chandrasekhar (1931) should be the specific white dwarf mass limit calculated for a polytrope. The insight that a relativistic treatment of degeneracy leads to the existence of a white dwarf mass limit first appeared in papers of W. Anderson (1929) and E.C. Stoner (1930) for a uniform density star. Accordingly, Chandrasekhar (1931) cites Stoner (1930) and points out that the polytrope white dwarf mass limit is less than Stoner's uniform density white dwarf mass limit by about 20%.

The origin of highly-magnetized white dwarfs has remained a mystery since their initial discovery. Recent observations indicate that the formation of high-field magnetic white dwarfs is intimately related to strong binary interactions during post-main-sequence phases of stellar evolution. If a low-mass companion, such as a planet, brown dwarf, or low-mass star is engulfed by a post-main-sequence giant, the hydrodynamic drag in the envelope of the giant leads to a reduction of the companion's orbit. Sufficiently low-mass companions in-spiral until they are shredded by the strong gravitational tides near the white dwarf core. Subsequent formation of a super-Eddington accretion disk from the disrupted companion inside a common envelope can dramatically amplify magnetic fields via a dynamo. Here, we show that these disk-generated fields are sufficiently strong to explain the observed range of magnetic field strengths for isolated, high-field magnetic white dwarfs. A higher-mass binary analogue may also contribute to the origin of magnetar fields.

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.

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.

If accretion disc emission results from turbulent dissipation, then axisymmetric accretion theory must be used as a mean field theory: turbulent flows are at most axisymmetric only when suitably averaged. Spectral predictions therefore have an intrinsic imprecision that must be quantified to interpret the variability exhibited by a source observed at different epochs. We quantify contributions to the stochastic imprecision that come from azimuthal and radial averaging and show that the imprecision is minimized for a particular choice of radial averaging, which in turn, corresponds to an optimal spectral resolution of a telescope for a spatially unresolved source. If the optimal spectral resolution is less than that of the telescope then the data can be binned to compare to the theoretical prediction of minimum imprecision. Little stochastic variability is predicted at radii much larger than that at which the dominant eddy turnover time ($\sim$ orbit time) exceeds the time interval between observations; the epochs would then be sampling the same member of the stochastic ensemble. We discuss the application of these principles to protoplanetary discs for which there is presently a paucity of multi-epoch data but for which such data acquisition projects are underway.

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.

We further develop the model of molecular cloud fragmentation introduced in Field, Blackman and Keto (2007; FBK). We show that external pressure acting on fragments establishes a scale-dependent critical mass. Fragments with masses less than the critical value are confined largely by pressure, while those with masses greater than or equal to the critical value collapse under self gravitation. Both types of fragments are commonly observed. Without specifying the source of the external pressure, and without assuming any other scaling relations, we predict the power - law index in the relation between the rms velocity of supersonic motions and the size of fragments . We then investigate the possibility that the external pressure is due to the kinetic energy of H atoms released by photodissociation of hydrogen molecules in the fragment. This can account approximately for the observed values of external pressure and two additional observations: the value of the scaling coefficient in the power law mentioned above, and the observation of outflowing atomic hydrogen around molecular clouds. A further prediction is HI at fragment edges with column densities of order 1E20 per sq. cm and velocities of a few km/s that should be detectable with high resolution 21 cm observations. Finally, we predict the magnitude of the coefficient of dissipation in the observed supersonic flows.

The origin of large scale magnetic fields in astrophysical rotators, and the conversion of gravitational energy into radiation near stars and compact objects via accretion have been subjects of active research for a half century. Magnetohydrodynamic turbulence makes both problems highly nonlinear, so both subjects have benefitted from numerical simulations.However, understanding the key principles and practical modeling of observations warrants testable semi-analytic mean field theories that distill the essential physics. Mean field dynamo (MFD) theory and alpha-viscosity accretion disc theory exemplify this pursuit. That the latter is a mean field theory is not always made explicit but the combination of turbulence and global symmetry imply such. The more commonly explicit presentation of assumptions in 20th century textbook MFDT has exposed it to arguably more widespread criticism than incurred by 20th century alpha-accretion theory despite complementary weaknesses. In the 21st century however, MFDT has experienced a breakthrough with a dynamical saturation theory that consistently agrees with simulations. Such has not yet occurred in accretion disc theory, though progress is emerging. Ironically however, for accretion engines, MFDT and accretion theory are presently two artificially uncoupled pieces of what should be a single coupled theory. Large scale fields and accretion flows are dynamically intertwined because large scale fields likely play a key role in angular momentum transport. I discuss and synthesize aspects of recent progress in MFDT and accretion disc theory to suggest why the two likely conspire in a unified theory.

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.

Non-thermal X-ray emission in compact accretion engines can be interpreted to result from magnetic dissipation in an optically thin magnetized corona above an optically thick accretion disk. If coronal magnetic field originates in the disk and the disk is turbulent, then only magnetic structures large enough for their turbulent shredding time to exceed their buoyant rise time survive the journey to the corona. We use this concept and a physical model to constrain the minimum fraction of magnetic energy above the critical scale for buoyancy as a function of the observed coronal to bolometric emission. Our results suggest that a significant fraction of the magnetic energy in accretion disks resides in large scale fields, which in turn provides circumstantial evidence for significant non-local transport phenomena and the need for large scale magnetic field generation. For the example of Seyfert AGN, for which of order 30 per cent of the bolometric flux is in the X-ray band, we find that more than 20 per cent of the magnetic energy must be of large enough scale to rise and dissipate in the corona.

Faraday rotation (FR) is widely used to infer the orientation and strength of magnetic fields in astrophysical plasmas. Although the absence of electron-positron pairs is a plausible assumption in many astrophysical environments, the magnetospheres of pulsars and black holes and their associated jets may involve a significant pair plasma fraction. This motivates being mindful of the effect of positrons on FR. Here we derive and interpret exact expressions of FR for a neutral plasma of arbitrary composition. We focus on electron-ion-positron plasmas in which charge neutrality is maintained by an arbitrary combination of ions and positrons. Because a pure electron-positron plasma has zero FR, the greater the fraction of positrons the higher the field strength required to account for the same FR. We first obtain general formulae and then specifically consider parameters relevant to active galctic nuclei (AGN) jets to illustrate the significant differences in field strengths that FR measurements from radio frequency measurements. Complementarily, using galaxy cluster core plasmas as examples, we discuss how plasma composition can be constrained if independent measurements of the field strength and number density are available and combined with FR.

Models of accretion discs and their associated outflows often incorporate assumptions of axisymmetry and symmetry across the disc plane. However, for turbulent discs these symmetries only apply to averaged quantities and do not apply locally. The local asymmetries can induce local imbalances in outflow power across the disc mid-plane, which can in turn induce local tilting torques. Here we calculate the effect of the resulting stochastic torques on disc annuli that are a consequence of standard mean field accretion disc models. The torques induce a random walk of the vector perpendicular to the plane of each averaged annulus. This random walk is characterized by a radially dependent diffusion coefficient which we calculate for small angle tilt. We use the coefficient to calculate a radially dependent time scale for annular tilt and associated jet wobble. The wobble time depends on the square of the wander angle so the age of a given system determines the maximum wobble angle. We apply this to examples of blazars, young stellar objects and binary engines of pre-planetary nebulae and microquasars. It is noteworthy that for an averaging time $t_w\sim 3$ days, we estimate a wobble angle for jets in SS433 of $\theta\sim 0.8$ degrees, not inconsistent with observational data. In general the non-periodic nature of the stochastic wobble could distinguish it from faster periodic jet precession.

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.

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.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's Astronomy and Astrophysics 2010 Decadal Survey. The white paper is endorsed by the American Physical Society's (APS) Topical Group on Plasma Astrophysics (GPAP).

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.

Many, if not all, post AGB stellar systems swiftly transition from a spherical to a powerful aspherical pre-planetary nebula (pPNE) outflow phase before waning into a PNe. The pPNe outflows require engine rotational energy and a mechanism to extract this energy into collimated outflows. Just radiation and rotation are insufficient but a symbiosis between rotation, differential rotation and large scale magnetic fields remains promising. Present observational evidence for magnetic fields in evolved stars is suggestive of dynamically important magnetic fields, but both theory and observation are rife with research opportunity. I discuss how magnetohydrodynamic outflows might arise in pPNe and PNe and distinguish different between approaches that address shaping vs. those that address both launch and shaping. Scenarios involving dynamos in single stars, binary driven dynamos, or accretion engines cannot be ruled out. One appealing paradigm involves accretion onto the primary post-AGB white dwarf core from a low mass companion whose decaying accretion supply rate owers first the pPNe and then the lower luminosity PNe. Determining observational signatures of different MHD engines is a work in progress. Accretion disk theory and large scale dynamos pose many of their own fundamental challenges, some of which I discuss in a broader context.

The available information on isotopic abundances in the atmospheres of low-mass Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars requires that episodes of extensive mixing occur below the convective envelope, reaching down to layers close to the hydrogen burning shell (Cool Bottom Processing). Recently \cite{Busso:2007jw} suggested that dynamo-produced buoyant magnetic flux tubes could provide the necessary physical mechanisms and also supply sufficient transport rates. Here, we present an $\alpha-\Omega$ dynamo in the envelope of an RGB/AGB star in which shear and rotation drain via turbulent dissipation and Poynting flux. In this context, if the dynamo is to sustain throughout either phase, convection must resupply shear. Under this condition, volume-averaged, peak toroidal field strengths of $<B_\phi>\simeq3\times10^3$ G (RGB) and $<B_\phi>\simeq5\times10^3$ G (AGB) are possible at the base of the convection zone. If the magnetic fields are concentrated in flux tubes, the corresponding field strengths are comparable to those required by Cool Bottom Processing.

Supermassive black holes in giant elliptical galaxies are remarkably faint given their expected accretion rates. This motivates models of radiatively inefficient accretion, due to either ion-electron thermal decoupling, generation of outflows that inhibit accretion, or settling of gas to a gravitationally unstable disk that forms stars in preference to feeding the black hole. The latter model predicts the presence of cold molecular gas in a thin disk around the black hole. Here we report Submillimeter Array observations of the nucleus of the giant elliptical galaxy M87 that probe 230 GHz continuum and CO(J=2--1) line emission. Continuum emission is detected from the nucleus and several knots in the jet, including one that has been undergoing flaring behavior. We estimate a conservative upper limit on the mass of molecular gas within ~100pc and +-400km/s line of sight velocity of the central black hole of ~8x10^6Msun, which includes an allowance for possible systematic errors associated with subtraction of the continuum. Ignoring such errors, we have a 3 sigma sensitivity to about 3x10^6Msun. In fact, the continuum-subtracted spectrum shows weak emission features extending up to 4 sigma above the RMS dispersion of the line-free channels. These may be artifacts of the continuum subtraction process. Alternatively, if they are interpreted as CO emission, then the implied molecular gas mass is ~5x10^6Msun spread out over a velocity range of 700km/s. These constraints on molecular gas mass are close to the predictions of the model of self-gravitating, star-forming accretion disks fed by Bondi accretion (Tan & Blackman 2005).

Turbulent viscosity is frequently used in accretion disk theory to replace the microphysical viscosity in order to accomodate the observational need for in- stabilities in disks that lead to enhanced transport. However, simply replacing the microphysical transport coefficient by a single turbulent transport coeffi- cient hides the fact that the procedure should formally arise as part of a closure in which the hydrodynamic or magnetohydrodynamic equations are averaged, and correlations of turbulent fluctuations are replaced by transport coefficients. Here we show how a mean field approach leads quite naturally two transport coefficients, not one, that govern mass and angular momentum transport. In particular, we highlight that the conventional approach suffers from a seemingly inconsistent neglect of turbulent diffusion in the surface density equation. We constrain these new transport coefficients for specific cases of inward, outward, and zero net mass transport. In addition, we find that one of the new transport terms can lead to oscillations in the mean surface density which then requires a constant or small inverse Rossby number for disks to maintain a monotonic power-law surface density.

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 magnetic Reynolds number R_M, is defined as the product of a characteristic scale and associated flow speed divided by the microphysical magnetic diffusivity. For laminar flows, R_M also approximates the ratio of advective to dissipative terms in the total magnetic energy equation, but for turbulent flows this latter ratio depends on the energy spectra and approaches unity in a steady state. To generalize for flows of arbitrary spectra we define an effective magnetic dissipation number, R_{M,e}, as the ratio of the advection to microphysical dissipation terms in the total magnetic energy equation, incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, and distinct pairs of inner and outer scales for magnetic and kinetic spectra. As expected, for a substantial parameter range R_{M,e}\sim {O}(1) << R_M. We also distinguish R_{M,e} from {\tilde R}_{M,e} where the latter is an effective magnetic Reynolds number for the mean magnetic field equation when a turbulent diffusivity is explicitly imposed as a closure. That R_{M,e} and {\tilde R}_{M,e} approach unity even if R_M>>1 highlights that, just as in hydrodynamic turbulence,energy dissipation of large scale structures in turbulent flows via a cascade can be much faster than the dissipation of large scale structures in laminar flows. This illustrates that the rate of energy dissipation by magnetic reconnection is much faster in turbulent flows, and much less sensitive to microphysical reconnection rates compared to laminar flows.

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.

In low-mass Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars, anomalous mixing must transport material near the hydrogen-burning shell to the convective envelope. Recently, it was suggested that buoyant magnetic flux tubes could supply the necessary transport rate (Busso et al. 2007). The fields are assumed to originate from a dynamo operating in the stellar interior. Here, we show what is required of an $\alpha-\Omega$ dynamo in the envelope of an AGB star to maintain these fields. Differential rotation and rotation drain via turbulent dissipation and Poynting flux, so if shear can be resupplied by convection, then large-scale toroidal field strengths of $\left<B_\phi\right>\simeq3\times10^4$ G can be sustained at the base of the convection zone.

The physical mechanism(s) responsible for transitioning from a spherical Asymptotic Giant Branch (AGB) star to an asymmetric post-AGB (pAGB) object is poorly understood. In particular, excess momenta in the outflows of pAGB objects suggest that a binary may be required to supply an additional source of energy and angular momentum. The extraction of rotational energy from the engine is likely fundamental and may be facilitated if a dynamo is operating in the interior. In this regard, single star magnetic outflow models have been proposed as mechanisms for producing and shaping PNe, however these models neglect the back-reaction of the large-scale magnetic field on the flow. Here we present a transient $\alpha-\Omega$ dynamo operating in the envelope of an AGB star in (1) an isolated setting and (2) a common envelope in which the secondary is a low-mass companion in-spiraling in the AGB interior. The back reaction of the fields on the shear is included and differential rotation and rotation deplete via turbulent dissipation and Poynting flux. For an isolated star, the shear must be resupplied in order to sufficiently sustain the dynamo. We comment on the energy requirements that convection must satisfy to accomplish this. For the common envelope case, a robust dynamo can result as the companion provides an additional source of energy and angular momentum.

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.

By combining data from different published 3-D simulations of Keplerian shearing boxes unstable to the magnetorotational instability (MRI), we highlight tight anti-correlations between the total effective inferred angular momentum transport parameter, $\alpha_{tot}$, its separate Maxwell and Reynolds contributions $\alpha_{mag}$ and $\alpha_{kin}$, and the kinetic to magnetic pressure ratio $\beta$, defined with the initial or saturated (when available) thermal pressure. Plots of $Log (\alpha_{kin}), Log (\alpha_{mag})$, and $Log (\alpha_{tot})$ vs $Log (\beta)$ are well fit by straight lines even as $\alpha_{kin}$, $\alpha_{mag}$,and $\alpha_{tot}$ vary by four orders of magnitude over the simulations included. The ratio $\alpha_{kin}/\alpha_{mag}$ and the product $\alpha_{tot}\beta$ are quite constant and largely independent of the presence or absence of weak mean fields, the choice of initial and boundary conditions, and the resolution. In short, simulations have more strongly constrained the product $\alpha_{tot}\beta$ than $\alpha_{tot}$ itself.

We present a model in which the supersonic motions observed in molecular clouds are driven by gravitational energy released as large structures fragment into smaller ones. The fragmentation process begins in large molecular clouds, and continues down to fragments of a critical mass, at which gravitational confinement may be replaced by pressure confinement. The power laws that describe the scaling of density and mass, and number spectra and mass spectra of the fragments are given in terms of that of the observed velocity dispersion of the fragments. The results agree reasonably well with observations over the range from several to about a third of a million solar masses.

Observations suggest that many, if not all, post AGB systems evolve through an aspherical outflow phase. Such outflows require a sufficient engine rotational energy which binaries can provide. Via common envelope evolution, binaries can directly eject equatorial outflows or produce poloidal outflows from magnetized accretion disks around the primary or secondary. We discuss how accretion driven magnetohydrodynamic outflow models all make similar predictions for the outflow power and speed and we distinguish between the launch vs. propagation regimes of such outflows. We suggest that the high velocity bipolar outflows observed in planetary nebulae (PNe) and the lower velocity but higher power bipolar outflows observed in pre-PNe (pPNe) are kinematically consistent with time dependent accretion onto a white dwarf (WD) within a depleting envelope. Since the WD primary core is always present in all post-AGB systems, accretion onto this core is potentially common. Previous work has focused on core accretion from sub-stellar companions, but low mass stellar companions may be more important, and further work is needed.

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.

We consider the dissipation by Fermi acceleration of magnetosonic turbulence in the Reynolds Layer of the interstellar medium. The scale in the cascade at which electron acceleration via stochastic Fermi acceleration (STFA) becomes comparable to further cascade of the turbulence defines the inner scale. For any magnetic turbulent spectra equal to or shallower than Goldreich-Sridhar this turns out to be $\ge 10^{12}$cm, which is much larger than the shortest length scales observed in radio scintillation measurements. While STFA for such spectra then contradict models of scintillation which appeal directly to an extended, continuous turbulent cascade, such a separation of scales is consistent with the recent work of \citet{Boldyrev2} and \citet{Boldyrev3} suggesting that interstellar scintillation may result from the passage of radio waves through the galactic distribution of thin ionized boundary surfaces of HII regions, rather than density variations from cascading turbulence. The presence of STFA dissipation also provides a mechanism for the non-ionizing heat source observed in the Reynolds Layer of the interstellar medium \citep{Reynolds}. STFA accommodates the proper heating power, and the input energy is rapidly thermalized within the low density Reynolds layer plasma.

Most large scale dynamo research for astrophysical rotators focuses on interior flow driven helical dynamos (FDHDs), but larger scale coronal fields most directly influence observations. It is thus important to understand the relationship between coronal and interior fields. Coronal field relaxation is actually a type of magnetically dominated helical dynamo (MDHD). MDHDs also occur in fusion plasma devices where they drive a system toward its relaxed state in response to magnetic helicity injection that otherwise drives the system away from this state. Global scale fields of astrophysical rotators and jets are thus plausibly produced by a direct coupling between an interior FDHD and a coronal MDHD, interfaced by magnetic helicity transport through their mutual boundary. Tracking the magnetic helicity also elucidates how both FDHD and MDHDs evolve and saturate. The utility of magnetic helicity is unhampered by its non-gauge invariance since physical fields can always be recovered.

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 present preliminary results of particle-in-cell simulations of collisionless magnetic reconnection for conditions that apply to magnetically dominated pair plasma jets such as those in gamma-ray bursts. We find similar behaviour to previous authors but with reduced reconnection efficiency. This results because we include the full electromagnetic field dynamically and allow electric field fluctuations to grow. Although weak, these fluctuations impede early x-point formation by periodically decelerating and even reversing the inflow.

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 linear growth rate and saturation level of magnetic fields for Weibel instabilities driven by ion temperature anisotropy, defined as $\alpha=(T_\perp/T_\parallel)-1$ where $T_\perp$ and $T_\parallel$ are ion temperatures perpendicular and parallel to the wave vector, are derived in the small $\alpha$-limit. It is shown that the ratio of the saturated magnetic energy to the initial ion energy scales as the fourth power of the electron to ion mass ratio, $m/M$, for an initially unmagnetized plasma with $\alpha\leq M/m$. Particle-in-cell simulations confirm the mass scaling and also show that the electron energy gain is of the same order of magnitude as the magnetic field energy. This implies that the Weibel instabilities cannot provide a faster-than-Coulomb collisionless mechanism to equilibrate ion-electron plasmas with ions initially much hotter than electrons, a key component in low-luminosity astrophysical accretion flows. The results here also show that the large $\alpha$-limit formulae used in the study of magnetic field generation in collisionless shocks are only valid if $\alpha\geq M/m$.

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.

It is difficult to imagine a planet formation model that does not at some stage include a gravitationally unstable disc. Initially unstable gas-dust discs may form planets directly, but the high surface density required has motivated the alternative that gravitational instability occurs in a dust sub-layer only after grains have grown large enough by electrostatic sticking. Although such growth up to the instability stage is efficient for laminar discs, concern has mounted as to whether realistic disc turbulence catastrophically increases the settling time, thereby requiring additional processes to facilitate planet formation on the needed time scales. To evaluate this concern, we develop a model for grain growth that accounts for the influence of turbulence on the collisional velocity of grains and on the scale height of the dust layer. The relative effect on these quantities depends on the grain size. The model produces a disc-radius dependent time scale to reach a gravitationally unstable phase of planet formation. For a range of dust sticking and disc parameters, we find that for viscosity parameters $\alpha \le 10^{-3}$, this time scale is short enough over a significant range in radii $R$ that turbulence does not catastrophically slow the early phases of planet formation, even in the absence of agglomeration enhancement agents like vortices.

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.

The absence of other viable momentum sources for collimated flows leads to the likelihood that magnetic fields play a fundamental role in jet launch and/or collimation in astrophysical jets. To best understand the physics of jets, it is useful to distinguish between the launch region where the jet is accelerated and the larger scales where the jet propagates as a collimated structure. Observations presently resolve jet propagation, but not the launch region. Simulations typically probe the launch and propagation regions separately, but not both together. Here, I identify some of the physics of jet launch vs. propagation and what laboratory jet experiments to date have probed. Reproducing an astrophysical jet in the lab is unrealistic, so maximizing the benefit of the experiments requires clarifying the astrophysical connection.

Active Galactic Nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centers implies that there will be many jet/star interactions, which can mass-load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass-loading by stellar winds for a broader spectrum of wind mass-loss rates than has been previously considered. Given the observed stellar mass distributions in galactic centers, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses $M_{BH} \la 10^4M_{\odot}$ are rapidly mass loaded and quenched by stellar winds. For $10^4 M_{\odot}<M_{BH}<10^8 M_{\odot}$, the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass-loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon dominated composition on scales $\ga 2$ kpc therein even if the jet is initially pair plasma dominated.

A significant fraction of planetary nebulae (PNe) and proto-planetary nebulae (PPNe) exhibit aspherical, axisymmetric structures, many of which are highly collimated. The origin of these structures is not entirely understood, however recent evidence suggests that many observed PNe harbor binary systems, which may play a role in their shaping. In an effort to understand how binaries may produce such asymmetries, we study the effect of low-mass (< 0.3 M_sun) companions (planets, brown dwarfs and low-mass main sequence stars) embedded into the envelope of a 3.0 M_sun star during three epochs of its evolution (Red Giant Branch, Asymptotic Giant Branch (AGB), interpulse AGB). We find that common envelope evolution can lead to three qualitatively different consequences: (i) direct ejection of envelope material resulting in a predominately equatorial outflow, (ii) spin-up of the envelope resulting in the possibility of powering an explosive dynamo driven jet and (iii) tidal shredding of the companion into a disc which facilitates a disc driven jet. We study how these features depend on the secondary's mass and discuss observational consequences.

The term ``dynamo'' means different things to the laboratory fusion plasma and astrophysical plasma communities. To alleviate the resulting confusion and to facilitate interdisciplinary progress, we pinpoint conceptual differences and similarities between laboratory plasma dynamos and astrophysical dynamos. We can divide dynamos into three types: 1. magnetically dominated helical dynamos which sustain a large scale magnetic field against resistive decay and drive the magnetic geometry toward the lowest energy state, 2. flow-driven helical dynamos which amplify or sustain large scale magnetic fields in an otherwise turbulent flow, and 3. flow-driven nonhelical dynamos which amplify fields on scales at or below the driving turbulence. We discuss how all three types occur in astrophysics whereas plasma confinement device dynamos are of the first type. Type 3 dynamos requires no magnetic or kinetic helicity of any kind. Focusing on type 1 and 2 dynamos, we show how different limits of a unified set of equations for magnetic helicity evolution reveal both types. We explicitly describe a steady-state example of a type 1 dynamo, and three examples of type 2 dynamos: (i) closed volume and time dependent; (ii) steady-state with open boundaries; (iii) time dependent with open boundaries.

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.

We propose a new two-stage model for acceleration of electrons in solar flares. In the first stage, electrons are accelerated stochastically in a post-reconnection turbulent downflow. The second stage is the reprocessing of a subset of these electrons as they pass through a weakly compressive fast shock above the apex of the closed flare loop on their way to the chromosphere. We call this the "shock reprocessing" model. The model reproduces the energy dependent arrival time delays observed for both the pulsed and smooth components of impulsive solar flare x-rays with physically reasonable parameters for the downflow region. The model also predicts an emission site above the loop-top, as seen in the Masuda flare. The loop-top source distinguishes the shock reprocessing model from previous models. The model makes testable predictions for the energy dependence of footpoint pulse strengths and the location and spectrum of the loop-top emission, and can account for the observed soft-hard-soft trend in the spectral evolution of footpoint emission. Our model highlights the concept that reconnection is an acceleration environment rather than a single process. Which combination of processes operate may depend on the initial conditions that determine, for example, whether the reconnection downflow is turbulent. The shock reprocessing model comprises one such combination.

Observational evidence for anisotropy in supernovae (SN) may signal the importance of angular momentum and differential rotation in the progenitors. Free energy in differential rotation and rotation can be extracted magnetically or via turbulent dissipation. The importance that magnetohydrodyamic jets and coronae may play in driving SN motivates understanding large scale dynamos in SN progenitors. We develop a dynamical large scale interface dynamo model in which the differential rotation and rotation deplete both through Poynting flux and turbulent diffusion. We apply the model to a differentially rotating core surrounded by a convection zone of a SN progenitor from a initial 15$M_\odot$ star. Unlike the Sun, the dynamo is transient because the differential rotation is primarily due to the initial collapse. Up to $\sim 10^{51}$erg can be drained into time-integrated Poynting flux and heat, the relative fraction of which depends on the relative amount of turbulence in the shear layer vs. convection zone and the fraction of the shear layer into which the magnetic field penetrates. Both sinks can help facilitate explosions and could lead to different levels of anisotropy and pulsar kicks. In all cases, the poloidal magnetic field is much weaker than the toroidal field, and the Poynting flux is lower than previous estimates which invoke the magnitude of the total magnetic energy. A signature of a large scale dynamo is that the oscillation of the associated Poynting flux on $\sim 1$ sec time scales, implying the same for the energy delivery to a SN.

Using particle-in-cell (PIC) simulations, we study the evolution of a strongly magnetized plasma slab propagating into a finite density ambient medium. Like previous work, we find that the slab breaks into discrete magnetic pulses. The subsequent evolution is consistent with diamagnetic relativistic pulse acceleration of \cite{liangetal2003}. Unlike previous work, we use the actual electron to proton mass ratio and focus on understanding trapping vs. transmission of the ambient plasma by the pulses and on the particle acceleration spectra. We find that the accelerated electron distribution internal to the slab develops a double-power law. We predict that emission from reflected/trapped external electrons will peak after that of the internal electrons. We also find that the thin discrete pulses trap ambient electrons but allow protons to pass through, resulting in less drag on the pulse than in the case of trapping of both species. Poynting flux dominated scenarios have been proposed as the driver of relativistic outflows and particle acceleration in the most powerful astrophysical jets.

The luminosities of the centers of nearby elliptical galaxies are very low compared to models of thin disc accretion to their black holes at the Bondi rate, typically a few hundredths to a few tenths of a solar mass per year. This has motivated models of inefficiently-radiated accretion that invoke weak electron-ion thermal coupling, and/or inhibited accretion rates due to convection or outflows. Here we point out that even if such processes are operating, a significant fraction of the accreting gas is prevented from reaching the central black hole because it condenses into stars in a gravitationally unstable disc. Star formation occurs inside the Bondi radius (typically ~100pc in giant ellipticals), but still relatively far from the black hole in terms of Schwarzschild radii. Star formation depletes and heats the gas disc, eventually leading to a marginally stable, but much reduced, accretion flow to the black hole. We predict the presence of cold (~100K), dusty gas discs, containing clustered H-alpha emission and occasional type II supernovae, both resulting from the presence of massive stars. Star formation accounts for several features of the M87 system: a thin disc, traced by H-alpha emission, is observed on scales of about 100pc, with features reminiscent of spiral arms and dust lanes; the star formation rate inferred from the intensity of H-alpha emission is consistent with the Bondi accretion rate of the system. Star formation may therefore help suppress accretion onto the central engines of massive ellipticals. We also discuss some implications for the fueling of the Galactic center and quasars.

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.

Non-axisymmetric accretion discs can lead to flux variability. Here we provide an analytic framework for modeling non-axisymmetric structures like hotspots and spiral waves and their influence on observed timing measurements. The presently unexplained Low-frequency Quasi-Periodic Oscillations (LFQPO), observed in X-ray binaries and cataclysmic variables, could be the result of such discs. Our framework serves as a guide to quantify the properties that non-axisymmetric structures produced by nonlinear accretion disc models must have in order to explain observed features such as LFQPOs. The results from our microquasar appliations also provide analogous predictions for X-ray modulation in active galactic nuclei. The formalism and physical interpration is of practical use for generic non-axisymmetric accretion disc systems.

In the Hamiltonian treatment of purely mechanical systems, the canonical and actual momentum of a particle are the same. In contrast, for a plasma of charged particles and electromagnetic fields, those two momenta are different. We show how this distinction is fundamental in identifying the limitations of a recent attempt by Binney (2003) to rule out two-temperature collisionless astrophysical accretion flows from Hamiltonian theory. This illustrates the Hamiltonian method for astrophysical plasmas, its relation to the equations of motion, and its role in practical calculations. We also discuss how the complete Hamiltonian treatment of a plasma should couple the particle motion to a fully dynamical treatment of the electromagnetic fields. Our results stand independent from the discussion of Quataert (2003) who argued that time scale calculated in Binney (2003) is not the equipartition time as claimed.

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.

Considering a plasma with an initially weak large scale field subject to nonhelical turbulent stirring, Zeldovich (1957), for two-dimensions, followed by others for three dimensions, and Zeldovich et al. (1983) have presented formulae of the form $<b^2>=f(R_M){Bbar}^2$. Such ``Zeldovich relations'' have sometimes been interpreted to provide steady-state relations between the energy associated with the fluctuating magnetic field and that associated with a large scale or mean field multiplied by a function $f$ that depends on spatial dimension and a magnetic Reynolds number $R_M$. Here we dissect the origin of these relations and pinpoint pitfalls that show why they are inapplicable to realistic, dynamical MHD turbulence and that they disagree with many numerical simulations. For 2-D, we show that when the total magnetic field is determined by a vector potential, the standard Zeldovich relation applies only transiently, characterizing a maximum possible value that the field energy can reach before necessarily decaying. in relation to a seed value $Bbar$. In 3-D, we show that the standard Zeldovich relations are derived by balancing subdominant terms. In contrast, balancing the dominant terms shows that the fluctuating field can grow to a value independent of $R_M$ and the initially imposed $Bbar$, as seen in numerical simulations. We also emphasize that these Zeldovich relations of nonhelical turbulence imply nothing about the amount mean field growth in a helical dynamo. In short, by re-analyzing the origin of the Zeldovich relations, we highlight that they are inapplicable to realistic steady-states of large $R_M$ MHD turbulence.

A unified, three-scale system of equations accommodating nonlinear velocity driven helical dynamos, as well as time-dependent relaxation of magnetically dominated unihelical or bihelical systems is derived and solved herein. When opposite magnetic helicities of equal magnitude are injected on the intermediate and small scales, the large scale magnetic helicity grows kinematically (independent of the magnetic Reynolds number) to equal that on the intermediate scale. For both free and driven relaxation large scale fields are rapidly produced. Subsequently, a dissipation-limited dynamo, driven by growth of small scale kinetic helicity, further amplifies the large scale field. The results are important for astrophysical coronae fed with bihelical structures by dynamos in their host rotators. The large scale for the rotator corresponds to the intermediate scale for the corona. That bihelical magnetic relaxation can produce global scale fields may help to explain the formation of astrophysical coronal holes and magnetohydrodynamic outflows.

We reexamine stochastic Fermi acceleration (STFA) in the low energy (Newtonian) regime in the context of solar flares. The particle energization rate depends a dispersive term and a coherent gain term. The energy dependence of pitch angle scattering is important for determining the electron energy spectrum. For scattering by whistler wave turbulence, STFA produces a quasi-thermal spectrum. A second well-constrained scattering mechanism is needed for STFA to match the observed 10-100keV non-thermal spectrum. We suggest that STFA most plausibly acts as phase one of a two phase particle acceleration engine in impulsive flares: STFA can match the thermal spectrum below 10kev, and possibly the power law spectrum between 10 and 100keV, given the proper pitch angle scattering. However, a second phase, such as shock acceleration at loop tops, is likely required to match the spectrum above the observed knee at 100keV. Understanding this knee, if it survives further observations, is tricky.

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.

A non-linear, time-dependent, magnetically driven dynamo theory which shows how magnetically dominated configurations can relax to become helical on the largest scale available is presented. Coupled time-dependent differential equations for large scale magnetic helicity, small scale magnetic helicity, velocity, and the electromotive force are solved. The magnetic helicity on small scales relaxes to drive significant large scale helical field growth on dynamical (Alfv\'en crossing) time scales, independent of the magnitude of finite microphysical transport coefficients, after which the growing kinetic helicity slows the growth to a viscously limited pace. This magnetically driven dynamo complements the nonlinear kinetic helicity driven dynamo; for the latter, the growing magnetic helicity fluctuations suppress, rather than drive, large scale magnetic helicity growth. A unified set of equations accommodates both types of dynamos.

We directly determine the sensitivity and time delay of Earth's surface temperature response to annual solar irradiance variations from 60 years of data. A two-layer energy balance model is developed to interpret the results. Explaining both the resulting low sensitivity and time delay of 1-2 months requires negative feedback.

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 numerically investigate the possibility that a close stellar encounter could account for the high inclinations of the Kuiper belt, as originally proposed by Ida, Larwood and Burkert, however we consider encounters with pericenters within those explored by previous works. A star of mass 0.2 Msol entering the solar system on a nearly parabolic, low inclination, retrograde orbit with perihelion of 50 AU can scatter ~30% of the Kuiper belt into moderate inclination (i>10 degrees) and eccentricity orbits, while leaving the rest of the belt at low eccentricities and inclinations. This scenario provides a possible explanation for the dichotomy of the Kuiper belt, accounting for the scattered and classical populations, including objects with high eccentricities and inclinations. High eccentricity objects that were not placed into resonance with Neptune are more likely to be removed from the system by subsequent long timescale evolution. Consequently this scenario suggests that some Plutinos could have originally been objects in circular orbits which were scattered into resonant high eccentricity and high inclination orbits by the stellar encounter. Objects could have been placed into stable resonant regions with Neptune by the encounter. However long lived regions of dynamical stability should differ from those selected by resonance capture. This scenario predicts the presence of objects in resonances with Neptune that are unlikely to have been captured via Neptune's migration.

Magnetically mediated disk outflows are a leading paradigm to explain winds and jets in a variety of astrophysical sources, but where do the fields come from? Since accretion of mean magnetic flux may be disfavored in a thin turbulent disk, and only fields generated with sufficiently large scale can escape before being shredded by turbulence, in situ field production is desirable. Nonlinear helical inverse dynamo theory can provide the desired fields for coronae and outflows. We discuss the implications for contemporary protostellar disks, where the MRI (magneto-rotational instability) can drive turbulence in the inner regions, and primordial protostellar disks, where gravitational instability drives the turbulence. We emphasize that helical dynamos are compatible with the magneto-rotational instability, and clarify the relationship between the two.

It is shown that in magnetohydrodynamics, the shapes of the magnetic energy and magnetic helicity spectra can influence whether the magnetic helicity resistive dissipation time scale is longer, equal to, or shorter than that of magnetic energy. The calculations highlight that magnetic helicity need not always dissipate significantly more slowly than magnetic energy. While this may be implicitly understood in detailed studies of magnetic helicity in plasma devices, it is not sufficiently emphasized in elementary discussions comparing magnetic helicity vs. magnetic energy evolution. Measuring magnetic energy and helicity decay times may provide a posteriori insight into the initial spectra.

It is argued that much of the observed magnetic helicity losses at the solar surface may represent a reduction of an otherwise more dominant nonlinearity of solar and stellar dynamos. This nonlinearity is proportional to the internal twist (as opposed to writhe) of helical and sigmoidal surface structures.

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.

For many astrophysical sources with jets, there is evidence for the contemporaneous presence of disks. In contrast, pulsars such as the Crab and Vela show jets but have not yet revealed direct evidence for accretion disks. Here we show that for such pulsars, an accretion disk radiating below detectable thresholds may simultaneously account for (1) observed deviations in the braking indices from that of the simple dipole, (2) observed pulsar timing ages, and (3) possibly even the jet morphology via a disk outflow that interacts with the pulsar wind within, collimating and/or redirecting it.

Are magnetic fields important in primordial star formation? Assuming that star formation occurs via an accretion disk that is turbulent, initially because of local gravitational instability, we calculate the disk structure for realistic accretion rates. We predict that local gravitational viscosity is able to drive accretion, without the disk fragmenting. We then estimate the rate of dynamo amplification of seed magnetic field. Turbulence in a stratified disk can be helical, with different signs of the helicity in each hemisphere. This provides a key ingredient for production of global scale magnetic fields whose sign of flux is sustained over many orbit times. The resulting fields can drive collimated protostellar outflows that reduce the star formation efficiency from the initial gas cloud, especially once the protostar has contracted to the main sequence, at ~100 Msun. We estimate that the outflows are powerful enough to eject some material from the host dark matter halo and to initiate relatively strong magnetization of the local intergalactic medium. Close to the protostar, the outflow acts to shield the disk and equatorial regions from radiative feedback, such as ionizing photons, and this may enable accretion up to relatively large stellar masses. We conclude that magnetic fields cannot be ignored from models of primordial star formation.

We investigate field-line separation in strong MHD turbulence using direct numerical simulations. We find that in the static-magnetic-field approximation the thermal conductivity in galaxy clusters is reduced by a factor of about 50 relative to the Spitzer thermal conductivity of a non-magnetized plasma. This value is too small for heat conduction to balance radiative cooling in clusters.

The plethora of asymmetric planetary nebulae and the curiously high momenta of pre-planetary nebula outflows suggest that rotational energy is extracted from the engines. Magneto-rotational outflows driven by dynamos might be operating therein. I summarize scenarios involving (a) an isolated AGB star and (b) a binary. The efficacy of (a) requires re-establishing differential rotation over the AGB star's lifetime, whereas (b) delivers the angular momentum at the end of the AGB phase when needed, and may involve turbulent disks. Both can produce fields that drive outflows of the required mechanical luminosity and momentum, though weak points and open questions require further study.

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 develop a simple mean field approach to the transport of a passive scalar for which the fundamental equation is a second order differential equation in the transported quantity, not a first order equation. Triple correlations are included, as they must be for any realistic description of turbulence. No correlation time enters the theory, only an eddy turnover time. The approach is simpler than standard approaches which incorporate triple correlations, but more realistic than Gaussian or short correlation time closures which do not. A similar approach has proven useful in magnetohydrodynamics.

Recent MHD dynamo simulations for magnetic Prandtl number $>1$ demonstrate that when MHD turbulence is forced with sufficient kinetic helicity, the saturated magnetic energy spectrum evolves from having a single peak below the forcing scale to become doubly peaked with one peak at the system (=largest) scale and one at the forcing scale. The system scale field growth is well modeled by a recent nonlinear two-scale nonlinear helical dynamo theory in which the system and forcing scales carry magnetic helicity of opposite sign. But a two-scale theory cannot model the shift of the small-scale peak toward the forcing scale. Here I develop a four-scale helical dynamo theory which shows that the small-scale helical magnetic energy first saturates at very small scales, but then successively saturates at larger values at larger scales, eventually becoming dominated by the forcing scale. The transfer of the small scale peak to the forcing scale is completed by the end of the kinematic growth regime of the large scale field, and does not depend on magnetic Reynolds number $R_M$ for large $R_M$. The four-scale and two-scale theories subsequently evolve almost identically, and both show significant field growth on the system and forcing scales that is independent of $R_M$. In the present approach, the helical and nonhelical parts of the spectrum are largely decoupled. Implications for fractionally helical turbulence are discussed.

We develop a model of thin turbulent accretion discs supported by magnetic pressure of turbulent magnetic fields. This applies when the turbulent kinetic and magnetic energy densities are greater than the thermal energy density in the disc. Whether such discs survive in nature or not remains to be determined, but here we simply demonstrate that self-consistent solutions exist when the alpha-prescription for the viscous stress, similar to that of the original Shakura-Sunyaev model, is used. We show that \alpha \sim 1 for the strongly magnetized case and we calculate the radial structure and emission spectra from the disc in the regime when it is optically thick. Strongly magnetized optically thick discs can apply to the full range of disc radii for objects < 10^{-2} of the Eddington luminosity or for the outer parts of discs in higher luminosity sources. In the limit that the magnetic pressure is equal to the thermal or radiation pressure, our strongly magnetized disc model transforms into the Shakura-Sunyaev model with \alpha=1. Our model produces spectra quite similar to those of standard Shakura-Sunyaev models. In our comparative study, we also discovered a small discrepancy in the spectral calculations of Shakura and Sunyaev (1973).

There is mounting evidence that the ejection of magnetic helicity from the solar surface is important for the solar dynamo. Observations suggest that in the northern hemisphere the magnetic helicity flux is negative. We propose that this magnetic helicity flux is mostly due to small scale magnetic fields; in contrast to the more systematic large scale field of the 11 year cycle, whose helicity flux may be of opposite sign, and may be excluded from the observational interpretation. Using idealized simulations of MHD turbulence as well as a simple two-scale model, we show that shedding small scale (helical) field has two important effects. (i) The strength of the large scale field reaches the observed levels. (ii) The evolution of the large scale field proceeds on time scales shorter than the resistive time scale, as would otherwise be enforced by magnetic helicity conservation. In other words, the losses ensure that the solar dynamo is always in the near-kinematic regime. This requires, however, that the ratio of small scale to large scale losses cannot be too small, for otherwise the large scale field in the near-kinematic regime will not reach the observed values.

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.

Over the past few years there has been growing interest in helical magnetic field structures seen at the solar surface, in coronal mass ejections, as well as in the solar wind. Although there is a great deal of randomness in the data, on average the extended structures are mostly left-handed on the northern hemisphere and right-handed on the southern. Surface field structures are also classified as dextral (= right bearing) and sinistral (= left bearing) occurring preferentially in the northern and southern hemispheres respectively. Of particular interest here is a quantitative measurement of the associated emergence rates of helical structures, which translate to magnetic helicity fluxes. In this review, we give a brief survey of what has been found so far and what is expected based on models. Particular emphasis is put on the scale dependence of the associated fields and an attempt is made to estimate the helicity flux of the mean field vs. fluctuating field.

Two questions about the solar magnetic field might be answered together once their connection is identified. The first is important for large scale dynamo theory: what prevents the magnetic backreaction forces from shutting down the dynamo cycle? The second question is: what determines the handedness of twist and writhe in magnetized coronal ejecta? Magnetic helicity conservation is important for answering both questions. Conservation implies that dynamo generation of large scale writhed structures is accompanied by the oppositely signed twist along these structures. The latter is associated with the backreaction force. We suggest that coronal mass ejections (CME's) simultaneously liberate small scale twist and large scale writhe of opposite sign, helping to prevent the cycle from quenching and enabling a net magnetic flux change in each hemisphere. Observations and helicity spectrum measurements from a simulation of a rising flux ribbon support this idea. We show a new pictorial of dynamo flux generation that includes the backreaction and magnetic helicity conservation by characterizing the field as a 2-D ribbon rather than a 1-D line.

Line profiles represent an important probe of the geometry and dynamics of the central accretion engines in AGN. Although a flat disk illuminated from above can explain the relativistic Fe K alpha line profile in some cases, current data, plausible disk physics, and future comparison to the iron and other AGN lines invite some consideration of other geometries. General features of line profiles calculated fom concave discs, warped discs, thick clumpy (inflow/outflow) discs, and disks with spiral structure, are briefly summarized.

Magnetic fields play a dual role in gamma-ray bursts (GRBs). First, GRB and afterglow spectra (the latter interpreted as emission from external shocks) imply synchrotron radiation in a magnetic field that is a significant fraction of equipartition with the particle energy density. Second, magnetized rotators with $\sim 10^{15}$ Gauss field may power GRB by transporting Poynting flux to large distances where it dissipates and also drives an external shock. The field amplification at external shocks and in the engine involve separate processes. External shock fields are likely either seeded by a pre-GRB wind, or are amplified by two-stream plasma instabilities with MHD turbulence playing a subsequent role. In the engine, the large scale fields are likely produced by MHD helical dynamos, since flux accretion cannot easily compete with turbulent diffusion, and because structures must be large enough to rise to coronae before diffusing. Why helical dynamos are feasible, and their relation to the magnetorotational instability are among the points discussed.

We develop a new nonlinear mean field dynamo theory that couples field growth to the time evolution of the magnetic helicity and the turbulent electromotive force, $\emfb$. We show that the difference between kinetic and current helicities emerges naturally as the growth driver when the time derivative of $\emfb$ is coupled into the theory. The solutions predict significant field growth in a kinematic phase and a saturation rate/strength that is magnetic Reynolds number dependent/independent in agreement with numerical simulations. The amplitude of early time oscillations provides a diagnostic for the closure.

The Chandra AO1 HETGS observation of the micro-quasar GRS 1915+105 in the low hard state reveals (1) neutral K absorption edges from Fe, Si, Mg, and S in cold gas, and (2) highly ionized (Fe XXV and Fe XXVI) absorption attributed to a hot disk, disk wind, or corona. The neutral edges reveal anomalous Si and Fe abundances which we attribute to surrounding cold material in/near the environment of GRS 1915+105. We also point out the exciting possibility for the first astrophysical detection of XAFS attributed to material in interstellar grains. We place constraints on the ionization parameter, temperature, and hydrogen equivalent number density of the absorber near the accretion disk based on the detection of the H- and He-like Fe absorption. Observed spectral changes in the ionized lines which track the light curve point to changes in both the ionizing flux and density of the absorber, supporting the presence of a flow. Details can be found in Lee et al., 2002, ApJ., 567, 1102

We supplement the mean field dynamo growth equation with the total magnetic helicity evolution equation. This provides an explicitly time dependent model for alpha quenching in dynamo theory. For dynamos without shear, this approach accounts for the observed large scale field growth and saturation in numerical simulations. After a significant kinematic phase, the dynamo is resistively quenched, i.e. the saturation time depends on the microscopic resistivity. This is independent of whether or not the turbulent diffusivity is resistively quenched. We find that the approach is also successful for dynamos that include shear and exhibit migratory waves (cycles). In this case however, whether or not the cycle period remains of the order of the dynamical time scale at large magnetic Reynolds numbers does depend how on how the turbulent magnetic diffusivity quenches. Since this is unconstrained by magnetic helicity conservation, the diffusivity is presently an input parameter. Comparison to current numerical experiments suggests a turbulent diffusivity that depends only weakly on the magnetic Reynolds number, but higher resolution simulations are needed.

Some recent results and open issues in magnetic dynamo theory are addressed. The distinction between small-scale and mean-field dynamo (MFD) action in forced turbulent flows is emphasized. Though useful, the MFD has been controversial. This is partly due to concerns about its need for helically forced turbulence, but mainly because simple "textbook" treatments are kinematic and linear. The non-linear backreaction of the growing magnetic field has been suspected to prematurely quench MFD action. To resolve the controversy, we must first understand those non-linear MFDs which can be numerically simulated. Recently, there has been progress on this front. For simple MFDs in closed systems, dynamical quenching models that incorporate a transfer of magnetic helicity between small and large scales agree reasonably well with fully 3-D numerical periodic box simulations. Unresolved issues such as the quenching of turbulent diffusion and the additional physics needed to apply these results to real systems are also discussed herein. The following summarizes conceptual progress in describing mean-field magnetic energy growth in the simplest MFDs: For a closed turbulent flow, the non-linear mean-field dynamo, is first fast and kinematic, then slow and dynamic, and magnetic helicity transfer makes it so.

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.

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.

We present a two-scale approximation for the dynamics of a nonlinear $\alpha^2$ dynamo. Solutions of the resulting nonlinear equations agree with the numerical simulations of Brandenburg (2001), and show that $\alpha$ is quenched by the buildup of magnetic helicity at the forcing scale $1/k_2$ as the $\alpha$ effect transfers it from the large scale $1/k_1$. For times $t > (k_1/k_2)R_{M,2}$ in eddy turnover units (where $R_{M,2}$ is the magnetic Reynolds number of the forcing scale), $\alpha$ is resistively limited in the form predicted for the steady-state case. However, for $t << R_{M,2}$, $\alpha$ takes on its kinematic value, independent of $R_{M,2}$, allowing the production of large-scale magnetic energy equal to $k_1/k_2$ times equipartition. Thus the dynamic theory of $\alpha$ predicts substantial "fast" growth of large-scale field despite being "slow" at large times.

Magnetic field amplification in astrophysics ultimately requires an understanding of magnetohydrodynamic turbulence. Kinetic helicity has long been known to be important for large scale field growth in forced MHD turbulence, and has been recently demonstrated numerically to be asymptotically consistent with slow mean field dynamo action in a periodic box. Here we show numerically that the magnetic spectrum at and below the forcing scale is also strongly influenced by kinetic helicity. We identify a critical value, $f_{h,crit}$ above which the magnetic spectrum develops maxima at wavenumber $= 1$ scale {\it and} at the forcing scale, For $f< f_{h,crit}$ the field peaks only at the resistive scale. Kinetic helicity may thus be important not only for generating a large scale field, but also for establishing observed peaks in magnetic spectra at the forcing scale. The turbulent Galactic disk provides an example where both large scale ($>$ supernova forcing scale) fields and small scale ($\le$ forcing scale, with peak at forcing scale) fields are observed. We discuss this, and the potential application to the protogalaxy, but also emphasize the limitations in applying our results to these systems.

The time-averaged 30 ks Chandra HETGS observation of the micro-quasar GRS 1915+105 in the low hard state reveals for the first time in this source neutral K absorption edges from Fe, Si, Mg, & S. Ionized resonance absorption from H-, and He-like Fe (XXV, XXVI), Ca XX and possibly emission from neutral Fe Kalpha and ionized Fe XXV (forbidden, or the resonance emission component of a P-Cygni profile) are also seen. We report the tentative detection of the first astrophysical signature of XAFS in the photoelectric edge of Si (and possibly Fe and Mg), attributed to material in grains. The large column densities measured from the neutral edges reveal anomalous Si and Fe abundances. Scenarios for which the anomalous abundances can be attributed to surrounding cold material associated with GRS 1915+105 and/or that the enrichment may signify either a highly unusual supernova/hypernova, or external supernova activity local to the binary are discussed. We attribute the ionized features to a hot disk, disk-wind, or corona environment. These features allow for constraints on the ionization parameter (log xi > 4.15), temperature (T > 2.4 x 10^6 K), and hydrogen equivalent number density (n > 10^{12} cm^{-3}) for this region. Variability studies with simultaneous RXTE data show that the light curve count rate tracks changes in the disk blackbody and the power-law flux. Spectral changes in the Chandra data also track the behavior of the light curve, and may point to changes in both the ionizing flux and density of the absorber. A 3.69 Hz QPO and weak first harmonic is seen in the RXTE data.

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 calculate the iron line profiles from accretion discs with spiral velocity structures around Schwarzschild black holes. We find that quasi-periodic bumps appear in the the profiles, thereby providing a test for spiral wave patterns. This study is motivated by recent work showing that spiral density waves can result from MHD instabilities even in non-self-gravitating discs, and by improved spectral resolution of forthcoming X-ray missions.

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.

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.

A Chandra AO1 30ks HETGS observation of the X-ray transient micro-quasar GRS 1915+105 reveals absorption edges and faint line emission over the HETG energy range. We find from a preliminary analysis evidence for prominent neutral K edges associated with iron, silicon, magnesium, and tentatively sulphur. The column densities assuming solar abundances are consistent with ~ few x $10^{22} cm^{-2}$ in excess of the Galactic value, and may point to surrounding cold material associated with GRS 1915+105. Neutral Fe K$\alpha$ emission, and ionized absorption from Fe XXV and Fe XXVI are resolved. We limit our discussion to the Chandra results.

Mean field dynamo theory is a leading candidate to explain the observed large scale magnetic fields of galaxies and stars. However, controversy arises over the extent of premature quenching by the backreaction of the growing field. We distinguish between rapid mean field dynamo action, which is required by astrophysical systems, and resistively limited action. We show how the flow of magnetic helicity is important for rapid action. Existing numerical and analytic work suggesting that mean field dynamos are prematurely quenched and resistively limited include approximations or boundary conditions which suppress the magnetic helicity flow from the outset. Thus they do not unambiguously reveal whether real astrophysical mean field dynamos are dynamically suppressed when the helicity flow is allowed. An outflow of helicity also implies an outflow of magnetic energy and so active coronae or winds should accompany mean field dynamos. Open boundaries alone may not be sufficient for rapid dynamo action and the additional physics of buoyancy and outflows may be required. Possible simulation approaches to test some of the principles are briefly discussed. Some limitations of the ``Zeldovich relation'' are also addressed.

In addition to the scalar Shakura-Sunyaev $\alpha_{ss}$ turbulent viscosity transport term used in simple analytic accretion disc modeling, a pseudoscalar transport term also arises. The essence of this term can be captured even in simple models for which vertical averaging is interpreted as integration over a half-thickness and one separately studies each hemisphere. The additional term highlights a complementarity between mean field magnetic dynamo theory and accretion disc theory treated as a mean field theory. Such pseudoscalar terms have been studied, and can lead to large scale magnetic field and vorticity growth. Here it is shown that vorticity can grow even in the simplest azimuthal and half-height integrated disc model, for which mean quantities depend only on radius. The simplest vorticity growth solutions seem to have scales and vortex survival times consistent those required for facilitating planet formation. Also it is shown that when the magnetic back-reaction is included to lowest order, the pseudoscalar driving the magnetic field growth and that driving the vorticity growth will behave differently with respect to shearing and non-shearing flows: the former can reverse sign in the two cases, while the latter will have the same sign.

Mean field dynamos may explain the origin of large scale magnetic fields of galaxies, but controversy arises over the extent of dynamo quenching by the growing field. Here we explain how apparently conflicting results may be mutually consistent, by showing the role of magnetic helicity conservation and boundary terms usually neglected. We estimate the associated magnetic energy flowing out of the Galaxy but emphasize that the mechanism of field escape needs to be addressed.

The central engines of active galactic nuclei (AGN) contain cold, dense material as well as hot X-ray emitting gas. The standard paradigm for the engine geometry is a cold thin disc sandwiched between hot X-ray coronae. Strong support for this geometry in Seyferts comes from the study of fluorescent iron line profiles, although the evidence is not ubiquitously air tight. The thin disc model of line profiles in AGN and in X-ray binaries should be bench marked against other plausible possibilities. One proposed alternative is an engine consisting of dense clouds embedded in an optically thin, geometrically thick X-ray emitting engine. This model is further motivated by studies of geometrically thick engines such as advection dominated accretion flows (ADAFs). Here we compute the reprocessed iron line profiles from dense clouds embedded in geometrically thick, optically thin X-ray emitting discs near a Schwarzchild black hole. We consider a range of cloud distributions and disc solutions, including ADAFs, pure radial infall, and bipolar outflows. We find that such models can reproduce line profiles similar to those from geometrically thin, optically thick discs and might help alleviate some of the problems encountered from the latter.

Magnetic fields likely play a fundamental intermediary role between gravity and radiation in many astrophysical rotators. They can, among other things, 1) induce and be amplified by turbulence, 2) energize coronae, 3) launch and collimate outflows in ``spring'' or ``fling'' mechanisms. The first is widely recognized to be important for angular momentum transport, but can also produce intrinsic variability and possibly vorticity growth. The second leads to the production of high energy flares, and also facilitates a test of general relativity from AGN observations. The third can operate from rotators with a large scale fields, though the origin of the requisite large scale fields is somewhat unresolved. I discuss these three points in more detail below, emphasizing some open questions.

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.

Poynting flux driven outflows from magnetized rotators are a plausible explanation for gamma-ray burst engines. We suggest a new possibility for how such outflows might transfer energy into radiating particles. We argue that the Poynting flux drives non-linearly unstable large amplitude electromagnetic waves (LAEMW) which ``break'' at radii $r_t \sim 10^{14} $ cm where the MHD approximation becomes inapplicable. In the ``foaming'' (relativisticly reconnecting) regions formed during the wave breaks the random electric fields stochastically accelerate particles to ultrarelativistic energies which then radiate in turbulent electromagnetic fields. The typical energy of the emitted photons is a fraction of the fundamental Compton energy $ \epsilon \sim f \hbar c/r_e $ with $f \sim 10^{-3}$ plus additional boosting due to the bulk motion of the medium. The emission properties are similar to synchrotron radiation, with a typical cooling time $\sim 10^{-4}$ sec. During the wave break, the plasma is also bulk accelerated in the outward radial direction and at larger radii can produce afterglows due to the interactions with external medium. The near equipartition fields required by afterglow models maybe due to magnetic field regeneration in the outflowing plasma (similarly to the field generation by LAEMW of laser-plasma interactions) and mixing with the upstream plasma.

We consider the backreaction of the magnetic field on the magnetic dynamo coefficients and the role of boundary conditions in interpreting whether numerical evidence for suppression is dynamical. If a uniform field in a periodic box serves as the initial condition for modeling the backreaction on the turbulent EMF, then the magnitude of the turbulent EMF and thus the dynamo coefficient $\a$, have a stringent upper limit that depends on the magnetic Reynolds number $R_M$ to a power of order -1. This is not a dynamic suppression but results just because of the imposed boundary conditions. In contrast, when mean field gradients are allowed within the simulation region, or non-periodic boundary are used, the upper limit is independent of $R_M$ and takes its kinematic value. Thus only for simulations of the latter types could a measured suppression be the result of a dynamic backreaction. This is fundamental for understanding a long-standing controversy surrounding $\alpha$ suppression. Numerical simulations which do not allow any field gradients and invoke periodic boundary conditions appear to show a strong $\alpha$ suppression (e.g. Cattaneo & Hughes 1996). Simulations of accretion discs which allow field gradients and allow free boundary conditions (Brandenburg & Donner 1997) suggest a dynamo $\alpha$ which is not suppressed by a power of $R_M$. Our results are consistent with both types of simulations.

We show that the presence of a steady $\ao$ dynamo in astrophysical rotators likely leads to an outflow of relative magnetic helicity and thus magnetic energy available for particle acceleration in a corona. The connection between energy and magnetic helicity arises because the generation of large scale field is linked to a generation of large scale magnetic helicity. In a steady state, the relative magnetic helicity associated with the large scale field can escape, accompanied by an equal and opposite small scale contribution from the field, since the total magnetic helicity in large magnetic Reynolds number flows is conserved. From the helicity flow, a lower limit on the magnetic energy deposited in the corona can be estimated. Steady coronal activity and the dissipation of magnetic energy is therefore a signature of an internal dynamo. Our theoretical estimate of the power delivered by a mean field dynamo is consistent with that inferred from observations to be delivered to the solar corona, the Galactic corona, and Seyfert I AGN coronae.

A generalized mean magnetic field induction equation for differential rotators is derived, including a compressibility, and the anisotropy induced on the turbulent quantities from the mean magnetic field itself and a mean velocity shear. Derivations of the mean field equations often do not emphasize that there must be anisotropy and inhomogeneity in the turbulence for mean field growth. The anisotropy from shear is the source of a term involving the product of the mean velocity gradient and the cross-helicity correlation of the isotropic parts of the fluctuating velocity and magnetic field, $\lb{\bfv}\cdot{\bfb}\rb^{(0)}$. The full mean field equations are derived to linear order in mean fields, but it is also shown that the cross-helicity term survives to all orders in the velocity shear. This cross-helicity term can obviate the need for a pre-existing seed mean magnetic field for mean field growth: though a fluctuating seed field is necessary for a non-vanishing cross-helicity, the term can produce linear (in time) mean field growth of the toroidal field from zero mean field. After one vertical diffusion time, the cross-helicity term becomes sub-dominant and dynamo exponential amplification/sustenance of the mean field can subsequently ensue. The cross-helicity term should produce odd symmetry in the mean magnetic field, in contrast to the usually favored even modes of the dynamo amplification in sheared discs. This may be important for the observed mean field geometries of spiral galaxies. The strength of the mean seed field provided by the cross- helicity depends linearly on the magnitude of the cross-helicity.

Flourescent iron line profiles currently provide the best diagnostic for active galactic nuclei (AGN) engine geometries. Here we construct a method for calculating the relativistic iron line profile from an arbitrarily warped accretion disc, illuminated from above and below by hard X-ray sources. This substantially generalises previous calculations of reprocessing by accretion discs by including non-axisymmetric effects. We include a relativistic treatment of shadowing by ray-tracing photon paths along Schwarzchild geodesics. We apply this method to two classes of warped discs, and generate a selection of resulting line profiles. New profile features include the possibility of sharper red, and softer blue fall-offs, a time varying line profile if the warp precesses about the disc, and some differences between `twisted' and `twist-free' warps. We discuss some qualitative implications of the line profiles in the context of Type I and II Seyfert AGN.

Spectra of Seyfert Is are commonly modelled as emission from an X-ray illuminated flat accretion disc orbiting a central black hole. This provides both a reprocessed and direct component of the X-ray emission as required by observations of individual objects and possibly a fraction of the cosmological X-ray background. There is some observational motivation to at least consider the role that an effectively concave disc surface might play: (1) a reprocessed fraction $\gsim 1/2$ in some Seyferts and possibly in the X-ray background, and (2) the commonality of a sharp iron line peak for Seyferts at 6.4KeV despite a dependence of peak location on inclination angle for flat disc models. Here it is shown that a concave disc may not only provide a larger total fraction of reprocessed photons, but can also reprocess a much larger fraction of photons in its outer regions when compared to a flat disc. This reduces the sensitivity of the 6.4KeV peak location to the inner disc inclination angle because the outer regions are less affected by Doppler and gravitational effects. If the X-ray source is isotropic, the reprocessed fraction is directly determined by the concavity. If the X-ray source is anisotropic, the location of iron line peak can still be determined by concavity but the total reflected fraction need not be as large as for the isotropic emitter case. The geometric calculations herein are applicable to general accretion disc systems illuminated from the center.

An unsolved problem in turbulent dynamo theory is the ``back reaction'' problem: to what degree does the mean magnetic field suppress the turbulent dynamo coefficients which are needed to drive its growth? The answer will ultimately derive from a combination of numerical and analytical studies. Here we show that analytic approaches to the dynamo and back reaction problems require one to separate turbulent quantities into two components: those influenced by the mean field (which are therefore anisotropic) and those independent of the mean field (and are therefore isotropic), no matter how weak the mean field is. Upon revising the standard formalism to meet this requirement, we find that: (1) The two types of components often appear in the same equation, so that standard treatments, which do not distinguish between them, are ambiguous. (2) The usual first-order smoothing approximation that is necessary to make progress in the standard treatment is unnecessary when the distinction is made. (3) In contrast to previous suggestions, the correction to the dynamo $\a$ coefficient found by Pouquet et al (1976) is actually independent of the mean field, and therefore cannot be interpreted as a quenching.

We extend the standard two-scale theory of the turbulent dynamo coefficient $\alpha$ to include the nonlinear back reaction of the mean field $\bar B$ on the turbulence. We calculate the turbulent emf as a power series in $\bar B$, assuming that the base state of the turbulence ($\bar B=0$) is isotropic, and, for simplicity, that the magnetic diffusivity equals the kinematic viscosity. The power series converges for all $\bar B$, and for the special case that the spectrum of the turbulence is sharply peaked in $k$, our result is proportional to a tabulated function of the magnetic Reynolds number $R_M$ and the ratio $\beta$ of $\bar B$ (in velocity units) to the rms turbulent velocity $v_0$. For $\beta\to 0$ (linear regime) we recover the results of Steenbeck et al. (1966) as modified by Pouquet et al. (1976). For $R_M\gg 1$, the usual astrophysical case, $\alpha$ starts to decrease at $\beta \sim 1$, dropping like $\beta^{-2}$ as $\beta \to \infty$. Hence for large $R_M$, $\alpha$ saturates at $\bar B\sim v_0$, as estimated by Kraichnan (1979), rather than at $\bar B\sim R^{-1/2}_Mv_0$, as inferred by Cattaneo and Hughes (1996) from their numerical simulations at $R_M$=100. We plan to carry out simulations with various values of $R_M$ to investigate the discrepency.

Two-temperature advection dominated accretion flow (ADAF) or hot ion tori (HIT) models help explain low luminosity stellar and galactic accreting sources and may complement observational support for black holes in nature. But low radiative efficiencies demand that ions receive a fraction $\eta \gsim 99%$ of energy dissipated in the turbulent accretion. The $\eta$ depends on the ratio of particle to magnetic pressure. If compressive modes of dissipation, like magnetic mirroring, dominate incompressible modes, then even when the pressure ratio is O(1), the required large $\eta$ can be attained. However the relative importance of compressive vs. incompressible modes is hard to estimate. The larger up in the turbulent cascade, the more compressible the turbulence. The relevant length scale for particle energization can be determined by equating the dominant eddy turnover time to the time for which an energy equal to that in the turbulence can be drained. Based on the large scales resulting from this estimate, it is suggested that compressive mirroring may be important. Also, regardless of the precise $\eta$ or dissipation mechanism, non-thermal protons seem natural in two- temperature discs because all dissipation mechanisms, and the use of an isotopic pressure, require wave-particle resonances that operate only on a subset of the particles. Finally, it is briefly mentioned how mirroring may help to generate an ADAF or HIT in the first place.

We report VLA 1 sigma upper limits of 1.5% and 3% on the intrinsic circular and linear fractional polarizations, respectively, of the water vapor maser emission 0.2 pc from the central engine of NGC 4258. A corresponding 0.5% upper limit on any Zeeman-splitting-induced circular polarization translates to a 1 sigma upper limit on the parallel, or toroidal, component of the magnetic field of 300 mG. Assuming magnetic and thermal pressure balance in the disk, this magnetic field upper limit corresponds to a model-dependent estimate of the accretion rate through the molecular disk of 10^-1.9 alpha solar masses per year for the case where the magnetic field lies along the line of sight.

The $\alpha$ turbulent viscosity formalism for accretion discs must be interpreted as a mean field theory. The extent to which the disc scale exceeds that of the turbulence determines the precision of the predicted luminosity $L_\nu$. The assumption of turbulence and use of $\a$ implies: (1) Field line stretching generates a magnetic pressure $\gsim \a^2/6$ of the total pressure generally, and a 1 to 1 relation between $\a$ and the pressure ratio when shearing instabilities dominate the viscosity. (2) Large eddy sizes and speeds in typical advection dominated accretion flows (ADAFs) lead to a lower precision in $L_\nu$ than for thin discs of a given total observation duration and central mass. The allowed variability (relative precision) at a particular frequency increases (decreases) with the size of the contributing region. For X-ray binary type ADAFs, the allowed variability is $\sim 5$% at $R \le 1000$ Schwarzchild radii for averages over $\gsim 1000$sec. But for large galactic nuclei like NGC 4258 and M87, the relative precision error can approach $50-100%$ even at $R \le 100 R_S$ for currently available observation durations. More data are then required to compare with ADAF predictions.

Cosmological gamma-ray bursts (GRBs) and their afterglows seem to result from dissipation of bulk energy in relativistic outflows, but their engine has not been unambiguously identified. The engine could be a young pulsar formed from accretion induced collapse with a dynamo amplified field. Elsewhere, we suggest that such a ``Usov type'' strong field pulsar may help explain the bimodal distribution in GRB durations. Here we discuss possible roles of a pulsar for the afterglow. We derive the expected bolometric luminosity decay. The extracted rotational energy could dissipate by shocks or by large amplitude electromagnetic waves (LAEMW). The simplest LAEMW approach predicts a slower decay in observed afterglow peak frequency and faster decay in flux than the simplest blast-wave model, though more complicated models of both can provide different dependences. LAEMW do not require the rapid magnetic field amplification demanded of the blast-wave approach because the emission originates from a nearly fixed radius. Different time dependent behavior of GRB and post-GRB emission is also predicted. Observational evidence for a pulsar in a GRB would make some GRB engine models, such as neutron star mergers and black holes unlikely. Therefore, the question of whether a pulsar is present is an important one even if it could drive a canonical fireball.

The origin and sustenance of large scale galactic magnetic fields has been a long standing and controversial astrophysical problem. Here an alternative to the ``standard'' $\a-\Omega$ mean field dynamo and primordial theories is pursued. The steady supply of supernovae induced turbulence exponentiates the total field energy, providing a significant seed mean field that can be linearly stretched by shear. The observed micro-Gauss fields would be produced primarily within one vertical diffusion time since it is only during this time that linear stretching can compete with diffusion. This approach does not invoke exponential mean field dynamo growth from the helicity $\a$-effect but does employ turbulent diffusion, which limits the number of large scale reversals. The approach could be of interest if the helicity effect is suppressed independently of the turbulent diffusion. This is an important but presently unresolved issue.

The importance of a teaching a clear definition of the ``observer'' in special relativity is highlighted using a simple astrophysical example from the exciting current research area of ``Gamma-Ray Burst'' astrophysics. The example shows that a source moving relativistically toward a single observer at rest exhibits a time ``contraction'' rather than a ``dilation'' because the light travel time between the source and observer decreases with time. Astrophysical applications of special relativity complement idealized examples with real applications and very effectively exemplify the role of a finite light travel time.

Cosmological gamma-ray bursts (GRBs) could be driven by dissipation of pure electromagnetic energy (Poynting flux) extracted from rapidly rotating compact objects with strong magnetic fields. One such possibility is a young millisecond pulsar (MSP) formed from accretion-induced collapse (AIC) of a white dwarf. The combination of an efficient magnetic dynamo, likely operating during the first seconds of the initially hot and turbulent MSP interior, and the subsequent modest beaming of gamma-ray emitting outflows, would easily account for energy constraints. But the remarkable feature of such models is that they may naturally explain the hitherto unexplained bimodal distribution in GRB time durations. The two burst classes could correspond to MSPs that form spinning above and below a gravitationally unstable limit respectively. In the former case, the spin-down time scale is due to gravitational radiation emission ($<1s$) while the spin-down time scale of the latter is due to electromagnetic dipole emission ($\gg 1s$). These two time scales account for the short and long GRB durations, i.e. the observed bimodal GRB duration distribution. A natural prediction is that the short duration GRBs would be accompanied by strong gravitational radiation emission which is absent from the longer class. Both would show millisecond variabilities.

We generalize the mean field magnetic dynamo to include local evolution of the mean vorticity in addition to the mean magnetic field. The coupled equations exhibit a general mean field dynamo instability that enables the transfer of turbulent energy to the magnetic field and vorticity on larger scales. The growth of the vorticity and magnetic field both require helical turbulence which can be supplied by an underlying global rotation. The dynamo coefficients are derived including the backreaction from the mean magnetic field to lowest order. We find that a mean vorticity field can actually seed exponential growth of mean magnetic field from only a small scale seed magnetic field, without the need for a seed mean magnetic field. The equations decouple when the fluctuating velocity and magnetic field cross-correlations vanish, resulting in the separate mean field magnetic and mean field vorticity dynamo equations.

We show that coronal magnetic dissipation in thin active sheets that sandwich standard thin accretion disks in active galactic nuclei may account for canonical electron temperatures of a few $\times 10^9$K if protons acquire most of the dissipated energy. Coulomb collisions transfer energy from the ions to the electrons, which subsequently cool rapidly by inverse-Compton scattering. In equilibrium, the proton energy density likely exceeds that of the magnetic field and both well exceed the electron and photon energy densities. The Coulomb energy transfer from protons to electrons is slow enough to maintain a high proton temperature, but fast enough to explain observed rapid X-ray variabilities in Seyferts. The $\sim 10^9$K electron temperature is insensitive to the proton temperature when the latter is $\ge 10^{12}$K.

We study accretion induced collapse of magnetized white dwarfs as an origin of millisecond pulsars. We apply magnetized accretion disk models to the pre-collapse accreting magnetic white dwarfs and calculate the white dwarf spin evolution. If the pulsar magnetic field results solely from the flux-frozen fossil white dwarf field, a typical millisecond pulsar is born with a field strength $\sim 10^{11}-10^{12}G$. The uncertainty in the field strength is mainly due to the uncertain physical parameters of the magnetized accretion disk models. A simple correlation between the pulsar spin $\Omega_*$ and the magnetic field $B_*$, $(\Omega_*/10^4s^{-1})\sim (B_{*}/10^{11}G)^{-4/5}$, is derived for a typical accretion rate $\sim 5\times 10^{-8}M_{\sun}/yr$. This correlation remains valid for a wide pre-collapse physical conditions unless the white dwarf spin and the binary orbit are synchronized prior to accretion induced collapse. We critically examine the possibility of spin-orbit synchronization in close binary systems. Using idealized homogeneous ellipsoid models, we compute the electromagnetic and gravitational wave emission from the millisecond pulsars and find that electromagnetic dipole emission remains nearly constant while millisecond pulsars may spin up rather than spin down as a result of gravitational wave emission. We also derive the physical conditions under which electromagnetic emission from millisecond pulsars formed by accretion induced collapse can be a source of cosmological gamma-ray bursts. We find that relativistic beaming of gamma-ray emission and precession of gamma-ray emitting jets are required unless the dipole magnetic field strengths are $>10^{15}$G; such strong dipole fields are in excess of those allowed from the accretion induced collapse formation process except in spin-orbit synchronization.

Observations show that magnetic reconnection and its slow shocks occur in solar flares. The basic magnetic structures are similar for long duration event (LDE) flares and faster compact impulsive (CI) flares, but the former require less non-thermal electrons than the latter. Slow shocks can produce the required non-thermal electron spectrum for CI flares by Fermi acceleration if electrons are injected with large enough energies to resonate with scattering waves. The dissipation region may provide the injection electrons, so the overall number of non-thermal electrons reaching the footpoints would depend on the size of the dissipation region and its distance from the chromosphere. In this picture, the LDE flares have converging inflows toward a dissipation region that spans a smaller overall length fraction than for CI flares. Bright loop-top X-ray spots in some CI flares can be attributed to particle trapping at fast shocks in the downstream flow, the presence of which is determined by the angle of the inflow field and velocity to the slow shocks.

We discuss the possibility that gamma-ray bursts may result from cosmological relativistic blob emitting neutron star jets that precess past the line of sight. Beaming reduces the energy requirements, so that the jet emission can last longer than the observed burst duration. One precession mode maintains a short duration time scale, while a second keeps the beam from returning to the line of sight, consistent with the paucity of repeaters. The long life of these objects reduces the number required for production as compared to short lived jets. Blobs can account for the time structure of the bursts. Here we focus largely on kinematic and time scale considerations of beaming, precession, and blobs--issues which are reasonably independent of the acceleration and jet collimation mechanisms. We do suggest that large amplitude electro-magnetic waves could be a source of blob acceleration.

The range of microphysical and global dynamical timescales in the central regions of Active Galactic Nuclei (AGN) is sufficiently wide to permit the existence of multiphase structure. In particular, very dense, cool clouds can coexist with a hot, magnetically-dominated medium and can thereby efficiently reprocess the continuum radiation generated in this primary source region. The strong dynamical forces in this central magnetosphere can give rise to extremely small clouds. Microphysical processes then determine whether such clouds can indeed survive, in spite of their extremely contrasting properties relative to the surrounding environment, for long enough to produce potentially observable thermal reprocessing signatures. We examine specific physical constraints on the thicknesses of such reprocessing clouds. Our results are plotted to show the range of conditions that is representative of the central regions of AGN. We find a parameter subspace in the extreme high density regime for which the effects of microphysical diffusion processes can be overcome and for which cool gas can maintain pressure equilibrium with the ambient magnetosphere.

Standard magnetohydrodynamic theories, such as the mean field dynamo theory, have been criticized when the back reaction of the magnetic field on turbulent motions is neglected. For the dynamo, this back reaction has been argued to suppress the turbulent motions required for optimal mean field production. Here it is suggested that if the magnetic field is spatially intermittent, for example residing in flux tubes, the back reaction on turbulent flows may be significantly reduced.

A torus of reduced differential rotation can form in the inner $\siml 10$pc core of active galactic nuclei incurring a density enhancement that can account for obscuration of X-rays in Seyferts when the initial inner core to black hole mass ratio $\gsim 50$. The same density enhancement and reduction in differential rotation can also lead to dynamo growth of poloidal fields which attain a magnitude $\sim 10^4$G when accreted onto the central engine. As radio jet models often employ poloidal fields as agents in extracting power for the jet luminosity, we suggest that jetted AGN might require this poloidal field production. Although the poloidal field would be originally produced in the obscuring torus, jetted objects are less likely to have obscuring tori: The poloidal field would only aid in powering jet emission after it accretes with the torus matter to the central engine. Thus, only during the relatively short torus accretion time scale could there be both a jet and and torus.

Many extended extragalactic radio sources require a local {\it in situ\/} acceleration mechanism for electrons, in part because the synchrotron lifetimes are shorter than the bulk travel time across the emitting regions. If the magnetic field in these sources is localized in flux tubes, reconnection may occur between regions of plasma $\be$ (ratio of particle to magnetic pressure) $<<1$, even though $\beta$ averaged over the plasma volume may be $\gsim 1$. Reconnection in low $\beta$ regions is most favorable to acceleration from reconnection shocks. The reconnection X-point regions may provide the injection electrons for their subsequent non-thermal shock acceleration to distributions reasonably consistent with observed spectra. Flux tube reconnection might therefore be able to provide $in\ situ$ acceleration required by large scale jets and lobes.

Reconnection shocks in a magnetically dominated plasma must be compressive. Non-thermal ion acceleration can occur across built-in slow shocks, and across outflow fast shocks when the outflow is supermagnetosonic and the field is line-tied. Electron acceleration may be initiated by injection from the dissipation region. Reconnection and shock acceleration thus cooperate and non-thermal acceleration should be a characteristic feature.

We derive the fully relativistic Ohm's law for an electron-positron plasma. The absence of non-resistive terms in Ohm's law and the natural substitution of the 4-velocity for the velocity flux in the relativistic bulk plasma equations do not require the field gradient length scale to be much larger than the lepton inertial lengths, or the existence of a frame in which the distribution functions are isotropic.