In a continuing effort to investigate the role of magnetic fields in evolved low and intermediate mass stars (principally regarding the shaping of their envelopes), we present new ALMA high resolution polarization data obtained for the nebula OH 231.8+4.2. We found that the polarized emission likely arises from aligned grains in the presence of magnetic fields rather than radiative alignment and self scattering. The ALMA data show well organized electric field orientations in most of the nebula and the inferred magnetic field vectors (rotated by 90 degrees) trace an hourglass morphology centred on the central system of the nebula. One region in the southern part of OH 231.8+4.2 shows a less organized distribution probably due to the shocked environment. These findings, in conjunction with earlier investigations (maser studies and dust emission analysis at other scales and wavelengths) suggest an overall magnetic hourglass located inside a toroidal field. We propose the idea that the magnetic field structure is closely related to the architecture of a magnetic tower and that the outflows were therefore magnetically launched. While the current dynamical effect of the fields might be weak in the equatorial plane principally due to the evolution of the envelope, it would still be affecting the outflows. In that regard, the measurement of the magnetic field at the stellar surface, which is still missing, combined with a full MHD treatment are required to better understand and constrain the events occurring in OH 231.8+4.2.

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

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.

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

We study the formation, evolution and physical properties of accretion disks formed via wind capture in binary systems. Using the AMR code AstroBEAR, we have carried out high resolution 3D simulations that follow a stellar mass secondary in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in order to properly resolve the formation of accretion disks around the secondary. We then compare simulations of binaries with three different orbital radii (10, 15, 20 AU). Disks are formed in all three cases, however the size of the disk and, most importantly, its accretion rate decreases with orbital radii. In addition, the shape of the orbital motions of material within the disk becomes increasingly elliptical with increasing binary separation. The flow is mildly unsteady with "fluttering" around the bow shock observed. The disks are generally well aligned with the orbital plane after a few binary orbits. We do not observe the presence of any large scale, violent instabilities (such as the flip-flop mode). For the first time, moreover, it is observed that the wind component that is accreted towards the secondary has a vortex tube-like structure, rather than a column-like one as it was previously thought. In the context of AGB binary systems that might be precursors to Pre-Planetary and Planetary Nebula, we find that the wind accretion rates at the chosen orbital separations are generally too small to produce the most powerful outflows observed in these systems if the companions are main sequence stars but marginally capable if the companions are white dwarfs. It is likely that many of the more powerful PPN and PN involve closer binaries than the ones considered here. The results also demonstrate principles of broad relevance to all wind-capture binary systems.

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.

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.

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.

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

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

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

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

We present 3-D MHD and synthetic numerical simulations to follow the evolution of randomly magnetized intra-cluster medium plasma under the effects of powerful, light, hypersonic and bipolar jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with power law energy spectrum tuned to the expectation for Kolmogorov turbulence. We investigate the power of jets and the viewing angle used for the synthetic Rotation Measure (RM) observations. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties; the average RM and the RM standard deviation are increased by the action of the jets. This may lead to overestimations of the CMFs' strength up to 70%. The effect correlates with the jet power. Jets distort and amplify CMFs especially near the edges of the lobes and the jets' heads. Thus the RM structure functions are flattened at scales comparable to the source size. Jet-produced RM enhancements depend on the orientation of the jet axis to the line of sight and are thus more apparent in quasars than in radio galaxies.

A combination of three-dimensional (3D) magnetohydrodynamics (MHD) and synthetic numerical simulations are presented to follow the evolution of a randomly magnetised plasma that models the intra-cluster medium (ICM), under the isolated effects of powerful, light, hypersonic and bipolar Fanaroff-Riley class II (FR II) jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with a Kolmogorov-like energy spectrum. Both the power of the jets and the viewing angle that is used for the synthetic Rotation Measure (RM) observations are investigated. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties which we analyse as a function of time as well as the viewing angle. The average RM and the RM standard deviation are increased by the action of the jets. Energetics, RM statistics and magnetic power spectral analysis consistently show that the effects also correlate with the jets' power, and that the lightest, fastest jets produce the strongest changes in their environment. We see jets distort and amplify the CMFs especially near the edges of the lobes and the jets' heads. This process leads to a flattening of the RM structure functions at scales comparable to the source size. The edge features we find are similar to ones observed in Hydra A. The results show that jet-produced RM enhancements are more apparent in quasars than in radio galaxies. Globally, jets tend to enhance the RM standard deviation which may lead to overestimations of the CMFs' strength by about 70%. This study means to serve as a pathfinder for the SKA, EVLA and LOFAR to follow the evolution of cosmic magnetic fields.

Radio observations of Fanaroff-Riley class II sources often show correlations between the synchrotron emission and the linear-polarimetric distributions. Magnetic position vectors seem to align with the projected emission of both the radio jets and the sources' edges. Using statistics we study such relation as well as its unknown time evolution via synthetic polarisation maps of model FR II sources formed in 3D-MHD numerical simulations of bipolar, hypersonic and weakly magnetised jets. The magnetic field is initially random with a Kolmogorov power spectrum, everywhere. We investigate the structure and evolution of magnetic fields in the sources as a function of the power of jets and the observational viewing angle. Our synthetic polarisation maps agree with observations, showing B-field vectors which are predominantly aligned with the jet axis, and show that magnetic fields inside sources are shaped by the jets' backflow. Polarimetry is found to correlate with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than inside more spherical ones. We see jets increase the magnetic energy in cocoons in proportion to the jet velocity and the cocoon width. Filaments in the synthetic emission maps suggest turbulence develops in evolved sources.

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.

3D-MHD numerical simulations of bipolar, hypersonic, weakly magnetized jets and synthetic synchrotron observations are presented to study the structure and evolution of magnetic fields in FR II radio sources. The magnetic field setup in the jet is initially random. The power of the jets as well as the observational viewing angle are investigated. We find that synthetic polarization maps agree with observations and show that magnetic fields inside the sources are shaped by the jets' backflow. Polarimetry statistics correlates with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than for ones with fatter cocoons. Jets increase the magnetic energy in cocoons, in proportion to the jet velocity. Both, filaments in synthetic emission maps and 3D magnetic power spectra suggest that turbulence develops in evolved sources.

We investigate the effects of intermittent strong jets from an Active Galactic Nuclei (AGN) of a massive galaxy in the core of a cool core galaxy cluster, on the dynamics and metal distribution of the intra-cluster medium (ICM). We use a simple model for the metal distribution within the host galaxy which includes metal injection via star formation. We carry out 2D axisymmetric hydrodynamic simulations of these systems. After having established a cooling flow, two light jets are injected in opposite directions with a range of (intermittent) active phases. We follow the time evolution of the system from the jets' active phases up to 3 Gyr. The general metallicity evolution for all our simulations is very similar on large-scales. The convective flows driven by the jets advect gas and metals from the central galaxy to distances beyond 1.5 Mpc within the cluster. Intermittent jets are able to distribute metals to grater radii. The metal injection has effects on the general metal abundances present in the ICM, the faster the metal replenishment time the higher the metal abundance. Although metallicity gradients in the very central regions of cool core clusters are likely to be shaped by less-energetic phenomena, we find evidence in our simulations for metallicity gradients similar to those observed out to ~400 kpc radius. The metal distribution details of the central galaxy have no effects on the dynamical evolution of the ICM metals.

The idea of a unified model for all astrophysical jets has been considered for quite some time. We present here a few scaling laws relevant to all type of astrophysical jets, analogous to those of \citet{sams96} which are widely used for astrophysical black holes. We use Buckingham's $ \Pi $ theorem of dimensional analysis to obtain a family of dimensional relations among the physical quantities associated to astrophysical jets.