Search

Very low-mass main-sequence stars reveal some curious trends in observed rotation period distributions that require abating the spin-down that standard rotational evolution models would otherwise imply. By dynamically coupling magnetically mediated spin-down to tidally induced spin-up from close orbiting substellar companions, we show that tides from sub-stellar companions may explain these trends. In particular, brown dwarf companions can delay the spin-down and explain the observed bimodality in rotation period distribution of old, late-type M stars. We find that tidal forces also strongly influence stellar X-ray activity evolution, so that methods of gyrochronological aging must be generalized for stars with even sub-stellar companions. We also discuss how the theoretical predictions of the spin evolution model can be used with future data to constrain the population distribution of companion orbital separations., Comment: 23 pages, 7 figures, 1 table. Submitted to Nature Astronomy

Magnetic reconnection is a ubiquitous fundamental process in space and astrophysical plasmas that rapidly converts magnetic energy into some combination of flow energy, thermal energy, and non-thermal energetic particles. Over the past decade, a new experimental platform has been developed to study magnetic reconnection using strong coil currents powered by high power lasers at low plasma beta, typical conditions under which reconnection is energetically important in astrophysics. KJ-class lasers were used to drive parallel currents to reconnect MG-level magnetic fields in a quasi-axisymmetric geometry, similar to the Magnetic Reconnection Experiment or MRX, and thus this platform is named micro-MRX. This presentation summarizes two major findings from micro-MRX: direct measurement of accelerated electrons and observation of ion acoustic waves during anti-parallel reconnection. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that direct acceleration by the out-of-plane reconnection electric field is at work. Furthermore, a sudden onset of ion acoustic bursts has been measured by collective Thomson scattering in the exhaust of magnetic reconnection, followed by electron acoustic bursts with electron heating and bulk acceleration. These results demonstrate that the micro-MRX platform offers a novel and unique approach to study magnetic reconnection in the laboratory in addition to the capabilities provided by traditional magnetized plasma experiments such as MRX and the upcoming FLARE (Facility for Laboratory Reconnection experiments). Future approaches to study other particle acceleration mechanisms and ion acoustic waves from magnetic reconnection are also discussed., Comment: 16 pages, 13 figures, 89 references, accepted for publication in Physics of Plasmas

We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum, and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at $\approx 100$ km/s. The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$. As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value $\approx 20$ km/s. We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements., Comment: Final Production Files accepted by the Journal of Plasma Physics on 18th May 2024. 9 figures

Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we discuss several features of the plasma structure in these experiments through order-of-magnitude models. First, we show that the observed rotation velocity would produce a centrifugal force strong enough to support the hollow density profile. Second, we show that the axial jet should diverge much faster than what was observed, were it not for a magnetized halo with 3T which surrounds the jet and exerts pressure on the interface.Finally, we discuss the temperature structure in the axial jet and plasma halo.We show that a 3T magnetic field would also suppress electron heat conduction,leading to the flat profile observed experimentally. We also find that the axial jet is efficiently radiatively cooled,whereas the halo is not, which would explain the thermal decoupling between the two regions., Comment: Second version with clarifications on text and figures, and explicit calculations regarding the omega-effect. Manuscript re-submitted to IEEE Transactions On Plasma Science. 4 Figures

We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calculation. We systematically vary plasma beta ($\beta_m$) in the flows, the strength of the cooling ($\Lambda_0$) and the exponent ($\alpha$) of temperature-dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a "spine" is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depends on initial flow parameters and provide a simple analytic model that captures the basic features of the flow., Comment: 14 pages, 16 figures. Submitted to MNRAS

Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding flows %radius in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows (SLOs) being deflected from the horizontal when the flows are of differing cross-section.

We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gradients launch a rotating plasma jet upwards, which is confined by a combination of ram, thermal, and magnetic pressure of a surrounding plasma halo. The jet has subsonic rotation, with a maximum rotation velocity $23 \pm 3$ km/s. The rotational velocity profile is quasi-Keplerian with a positive Rayleigh discriminant $\kappa^2 \propto r^{-2.8\pm0.8}$ rad$^2$/s$^2$. The plasma completes $0.5 - 2$ full rotations in the experimental time frame ($\sim 150$ ns)., Comment: 4 figures

Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently rapidly the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one and three dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period., Comment: 12 pages, 21 figures. Accepted for publication in MNRAS

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events., Comment: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process., Comment: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

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., Comment: 4 pages

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., Comment: 10 figures, 2 pages, revised for IAU323 proceedings

One in 5 planetary nebulae are ejected from common envelope binary interactions but Kepler Space Telescope results are already showing this proportion to be larger. Their properties, such as abundances can be starkly different from those of the general population, so they should be considered separately when using PN as chemical or population probes. Unfortunately post-common envelope PN cannot be discerned using only their morphologies, but this will change once we couple our new common envelope simulations with PN formation models., Comment: 4 pages, 3 figures, proceedings of the IAU Symposium 323, Planetary Nebulae: Multiwavelength Probes of Stellar and Galactic Evolution

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., Comment: 10 pages, 6 figures (submitted to MNRAS)

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., Comment: To appear in the proceedings of IAU Symposium 314 "Young Stars and Planets Near the Sun"

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., Comment: 6 pages, 5 figures, 1 table; accepted by Physical Review Letters

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., Comment: 27 pages, 16 figures, accepted for publication in A&A

We use computational simulations to compare the impact response of different football and U.S. Army helmet pad materials. We conduct experiments to characterize the material response of different helmet pads. We simulate experimental helmet impact tests performed by the U.S. Army to validate our methods. We then simulate a cylindrical impactor striking different pads. The acceleration history of the impactor is used to calculate the Head Injury Criterion for each pad. We conduct sensitivity studies exploring the effects of pad composition, geometry, and material stiffness. We find that: (1) The football pad materials do not outperform the currently used military pad material in militarily-relevant impact scenarios; (2) Optimal material properties for a pad depend on impact energy; and (3) Thicker pads perform better at all velocities. Our analysis suggests that by using larger helmet shells with correspondingly thicker pads, impact-induced traumatic brain injury may be significantly reduced. Keywords: helmet, pad, head injury, traumatic brain injury, head injury criterion, impact, Comment: 24 pages, 12 figures, 1 table

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., Comment: 12 pages, 34 figures

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., Comment: This paper has been withdrawn by the author due to some typos

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., Comment: Accepted to Proceedings of the National Academy of Sciences. Under PNAS embargo until time of publication

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., Comment: Submitted to MNRAS

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., Comment: 6 pages, 5 figures. Accepted for publication in Astrophysics & Space Science

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 $\simeq3\times10^3$ G (RGB) and $\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., Comment: Replaced to correct small error in published version: In \S 2.1, paragraphs 2 and 3 incorrectly refer to the poloidal field when qualitatively discussing magnetic diffusion in the shear zone. The correct physical interpretation is that the toroidal field diffuses through the shear zone consistent with the value of $\beta_\phi$

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., Comment: 18 pages, 4 figures, 3 tables; accepted to MNRAS

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\simeq3\times10^4$ G can be sustained at the base of the convection zone., Comment: 7 pages, 3 figures. To appear in AIP Proceedings of the IXth Torino Workshop on AGB Nucleosynthesis

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., Comment: 8 pages, 4 figures. To appear in Asymmetrical Planetary Nebulae IV. R. L. M. Corradi, A. Manchado, N. Soker eds - replaced to properly link to NASA ADS

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., Comment: 34 Pages, 10 Figure, Accepted for publication in Physics of Plasmas

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., Comment: 25 pages, 7 figures, submitted to MNRAS; revised to fix Fig. 5

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., Comment: 15 pages, 4 figures, accepted by Astrophysics and Space Science for Special Issue High Energy Density Laboratory Astrophysics Conference

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., Comment: 24 pages, 6 figures, submitted to MNRAS

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., Comment: accepted by ApJ

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., Comment: Accepted by MNRAS. 17 pages without figures. Full version with figures can be found at http://www.pas.rochester.edu/~afrank/labastro/MF230rv.pdf

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., Comment: 13 pages, accepted to MNRAS

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, Comment: Submitted to ApJ. For a version with figures go to http://web.pas.rochester.edu/~afrank/labastro/CW/Jet-Wind-Frank.pdf

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., Comment: 4 pages, 2 figures, to appear in Magnetic field and Helicity in the Sun and the Heliosphere, eds: D. Rust & B. Schmieder, ASP Conference Series, XXV JD3 (2003)

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., Comment: 10 pages, 5 figures, to appear in Adv. Space Sci. (Cospar 2002, ed. Buchner)

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)., Comment: 27 pages, 11 figures, Astron. Astroph. in press; shortened version accepted by A&A, all calculations and conclusions are unchanged

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, Comment: Proceedings of the 4th Microquasar Workshop, 2002, eds. Durouchoux, Fuchs & Rodriguez, published by the Center for Space Physics : Kolkata ; The main findings as presented in Lee et al., 2002, ApJ., 567, 1102 are summarized

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., Comment: 15 pages TeX, 1 ps figure submitted to ApJ

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., Comment: Submitted to the Astrophysical Journal Letters

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., Comment: Submitted to the Astrophysical Journal; for additional information, higher quality images, and simulation animations see http://www.pas.rochester.edu/~wma

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., Comment: 13 Pages, (inc 4 figs), LaTex, Submitted to ApJ Lett

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., Comment: 6 pages, 2 figures; ASP Conference Series, 2001, "Mass Outflow in Active Galactic Nuclei: New Perspectives", eds. D. M. Crenshaw, S. B. Kraemer, I. M. George

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., Comment: 6 pages, 4 figures, Contributed talk at "X-ray Astronomy 2000" meeting, Palermo, Italy, Sept. 2000. Proceedings to be published in ASP Conference series, eds. R. Giaconni, L. Stella, & S. Serio, 2000

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., Comment: 14 pages, 3 figures. To appear in ApJ, Part I

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., Comment: 5 pages LaTex, and 2 .ps figures, submitted to MNRAS, 4/97