We present radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C303. From this data we derive the magnetoplasma and electrodynamic parameters of this 50 kpc long jet. For a {$\sim 2$} kpc segment of this jet we obtain for the first time a direct determination of a {\it galactic}-scale electric current ($\sim 3 \times 10^{18}$ A), and its direction $-$ {\it positive} away from the AGN. Our analysis strongly supports a model where the jet energy flow is mainly electromagnetic.

The $\Omega$-phase of the liquid sodium $\alpha$-$\Omega$ dynamo experiment at NMIMT in cooperation with LANL has successfully demonstrated the production of a high toroidal field, $B_{\phi} \simeq 8\times B_r$ from the radial component of an applied poloidal magnetic field, $B_r$. This enhanced toroidal field is produced by rotational shear in stable Couette flow within liquid sodium at $Rm \simeq 120$. The small turbulence in stable Taylor-Couette flow is caused by Ekman flow where $ (\delta v/v)^2 \sim 10^{-3} $. This high $\Omega$-gain in low turbulence flow contrasts with a smaller $\Omega$-gain in higher turbulence, Helmholtz-unstable shear flows. This result supports the ansatz that large scale astrophysical magnetic fields are created within semi-coherent large scale motions in which turbulence plays only a smaller diffusive role that enables magnetic flux linkage.

In this White Paper, we emphasize the need for and the important role of plasma astrophysics in the studies of formation, evolution of, and feedback by Active Galaxies. We make three specific recommendations: 1) We need to significantly increase the resolution of VLA, perhaps by building an EVLA-II at a modest cost. This will provide the angular resolution to study jets at kpc scales, where, for example, detailed Faraday rotation diagnosis can be done at 1GHz transverse to jets; 2) We need to build coordinated programs among NSF, NASA, and DOE to support laboratory plasma experiments (including liquid metal) that are designed to study key astrophysical processes, such as magneto-rotational instability (origin of angular momentum transport), dynamo (origin of magnetic fields), jet launching and stability. Experiments allowing access to relativistic plasma regime (perhaps by intense lasers and magnetic fields) will be very helpful for understanding the stability and dissipation physics of jets from Supermassive Black Holes; 3) Again through the coordinated support among the three Agencies, we need to invest in developing comprehensive theory and advanced simulation tools to study the accretion disks and jets in relativistic plasma physics regime, especially in connecting large scale fluid scale phenomena with relativistic kinetic dissipation physics through which multi-wavelength radiation is produced.

It is shown that a dynamo can operate in an Active Galactic Nuclei accretion disk due to the Keplerian shear and due to the helical motions of expanding and twisting plumes of plasma heated by many star passages through the disk. Each plume rotates a fraction of the toroidal flux into poloidal flux, always in the same direction, through a finite angle, and proportional to its diameter. The predicted growth rate of poloidal magnetic flux, based upon two analytic approaches and numerical simulations, leads to a rapid exponentiation of a seed field, \sim 0.1 to \sim 0.01 per Keplerian period of the inner part of the disk. The initial value of the seed field may therefore be arbitrarily small yet reach, through dynamo gain, saturation very early in the disk history. Because of tidal disruption of stars close to the black hole, the maximum growth rate occurs at a radius of about 100 gravitational radii from the central object. The generated mean magnetic field, a quadrupole field, has predominantly even parity so that the radial component does not reverse sign across the midplane. The linear growth is predicted to be the same by each of the following three theoretical analyses: the flux conversion model, the mean field approach, and numerical modeling. The common feature is the conducting fluid flow, considered in companion Paper I (Pariev and Colgate 2006) where two coherent large scale flows occur naturally: the differential winding of Keplerian motion and differential rotation of expanding plumes.

A magnetic field dynamo in the inner regions of the accretion disk surrounding the supermassive black holes in AGNs may be the mechanism for the generation of magnetic fields in galaxies and in extragalactic space. We argue that the two coherent motions produced by 1) the Keplerian motion and 2) star-disk collisions, numerous in the inner region of AGN accretion disks, are both basic to the formation of a robust, coherent dynamo and consequently the generation of large scale magnetic fields. They are frequent enough to account for an integrated dynamo gain, e^{10^{9}} at 100 gravitational radii of a central black hole, many orders of magnitude greater than required to amplify any seed field no matter how small. The existence of extra-galactic, coherent, large scale magnetic fields whose energies greatly exceed all but massive black hole energies is recognized. In paper II (Pariev, Colgate, and Finn 2006) we argue that in order to produce a dynamo that can access the free energy of black hole formation and produce all the magnetic flux in a coherent fashion the existence of these two coherent motions in a conducting fluid is required. The differential winding of Keplerian motion is obvious, but the disk structure depends upon the model of "alpha", the transport coefficient of angular momentum chosen. The counter rotation of driven plumes in a rotating frame is less well known, but fortunately the magnetic effect is independent of the disk model. Both motions are discussed in this paper, paper I. The description of the two motions are preliminary to two theoretical derivations and one numerical simulation of the alpha-omega dynamo in paper II. (Abridged)

We propose that nearly every accelerated CR was part of the parallel current that maintains all force-free (f-f) magnetic fields. Charged particles are accelerated by the E-parallel (to the magnetic filed B) produced by reconnection. The inferred total energy in extra-galactic cosmic rays is 10^(60) ergs per galaxy spacing volume, provided that acceleration mechanisms assumed do not preferentially only accelerate ultra high energy cosmic rays (UHECRs). This total energy is about 10^5 times the parent galactic CR or magnetic energy. The formation energy of supermassive black holes (SMBHs) at galaxy centers, 10^(62) ergs, becomes the only feasible source. An efficient dynamo process converts gravitational free energy into magnetic energy in an accretion disk around a SMBH. Aided by Keplerian winding, this dynamo converts a poloidal seed field into f-f fields, which are transported into the general inter-galactic medium (IGM). This magnetic energy is also efficiently converted into particle energies, as evidenced by the radiation from energetic particles. CRs of the IGM are then the result of the continuing dissipation, in a Hubble time, of this magnetic energy, by acceleration within the f-f fields confined within the super-galactic walls and filaments of large scale structures. In addition, most UHECRs are diffusively lost to the galactic voids before the GZK attenuation time, 10^8 years and hence are lost from view. Similarly within the galaxy we expect that the winding by accretion are the sources of galactic CR acceleration. [edited]

A scenario is put forth for the formation of supermassive black holes at the centers of galaxies. It depends upon the formation of a Mestel disk with a flat rotation curve, where the mass is proportional to r and thickness to 1/r. Such disks should form from the collapse of uniformly rotating, isolated, gaseous clouds, either proto-galactic, galaxy-mass damped Lyman alpha clouds or the gas that survives galaxy mergers. We propose that in any case the disk will be unstable to the Rossby vortex instability (RVI). This instability grows from any large, steep pressure gradient in an optically thick disk. Such pressure gradients either occur adjacent to compact objects or could be triggered by individual supernovae in and around the disk. Upon excitation, the RVI transports angular momentum far more efficiently than the turbulence of a viscosity-based Shakura-Sunyaev disk. The critical column density condition in a L* Mestel disk predicts that the black hole mass is ~ 10^8 Msolar and is also proportional to the 1/4 power of the velocity dispersion, in accordance with observations.

Despite the importance of the magnetorotational instability (MRI) as a fundamental mechanism for angular momentum transport in magnetized accretion disks, it has yet to be demonstrated in the laboratory. A liquid sodium alpha-omega dynamo experiment at the New Mexico Institute of Mining and Technology provides an ideal environment to study the MRI in a rotating metal annulus (Couette flow). A local stability analysis is performed as a function of shear, magnetic field strength, magnetic Reynolds number, and turbulent Prandtl number. The later takes into account the minimum turbulence induced by the formation of an Ekman layer against the rigidly rotating end walls of a cylindrical vessel. Stability conditions are presented and unstable conditions for the sodium experiment are compared with another proposed MRI experiment with liquid gallium. Due to the relatively large magnetic Reynolds number achievable in the sodium experiment, it should be possible to observe the excitation of the MRI for a wide range of wavenumbers and further to observe the transition to the turbulent state.

A magnetic dynamo experiment is under construction at the New Mexico Institute of Mining and Technology. The experiment is designed to demonstrate in the laboratory the alpha-omega magnetic dynamo, which is believed to operate in many rotating and conducting astrophysical objects. The experiment uses the Couette flow of liquid sodium between two cylinders rotating with different angular velocities to model the omega-effect. The alpha-effect is created by the rising and expanding jets of liquid sodium driven through a pair of orifices in the end plates of the cylindrical vessel, presumably simulating plumes driven by buoyancy in astrophysical objects. The water analog of the dynamo device has been constructed and the flow necessary for the dynamo has been demonstrated. Results of the numerical simulations of the kinematic dynamo are presented. The toroidal field produced by the omega-effect is predicted to be B_{\phi} \simeq (R_m/2\pi) B_{poloidal}\simeq 20 \times B_{poloidal} for the expected magnetic Reynolds number of R_m \sim 120. The critical rate of jets necessary for the dynamo self-excitation is predicted from the calculations to be a pair of plumes every 4 revolutions of the outer cylinder. For reasonable technical limitations on the strength of materials and the power of the drive, the self-excitation of the dynamo appears to be feasible.

We study the evolution of a magnetic arcade that is anchored to an accretion disk and is sheared by the differential rotation of a Keplerian disk. By including an extremely low external plasma pressure at large distances, we obtain a sequence of axisymmetric magnetostatic equilibria and show that there is a fundamental difference between field lines that are affected by the plasma pressure and those are not (i.e., force-free). Force-free fields, while being twisted by the differential rotation of the disk, expand outward at an angle of $\sim 60^\circ$ away from the rotation axis, consistent with the previous studies. These force-free field lines, however, are enclosed by the outer field lines which originate from small disk radii and come back to the disk at large radii. These outer fields experience most of the twist, and they are also affected most by the external plasma pressure. At large cylindrical radial distances, magnetic pressure and plasma pressure are comparable so that any further radial expansion of magnetic fields is prevented or slowed down greatly by this pressure. This hindrance to cylindrical radial expansion causes most of the added twist to be distributed on the ascending portion of the field lines, close to the rotation axis. Since these field lines are twisted most, the increasing ratio of the toroidal $B_{\phi}$ component to the poloidal component $B_{R,z}$ eventually results in the collimation of magnetic energy and flux around the rotation axis. We discuss the difficulty with adding a large number of twists within the limitations of the magnetostatic approximation.

A quantitative analysis of two radio source samples having opposite extremes of ambient gas density leads to important new conclusions about the magnetic energy in the IGM. We conclude that giant sources in rarefied IGM environments, which contain magnetic energies E_B ~ 10^60-61 ergs, can be viewed as important "calorimeters" of the minimum energy a black hole (BH) accretion disk system injects into the IGM. In contrast to the radiation energy released by BH accretion, most of the magnetic energy is "trapped" initially in a volume, up to ~10^73 cm^3, around the host galaxy. But since these large, Mpc scale radio lobes are still overpressured after the AGN phase, their subsequent expansion and diffusion will magnetize a large fraction of the entire IGM. This suggests that the energy stored in intergalactic magnetic fields will have a major, as yet underestimated effect on the evolution of subsequently forming galaxies. Comparison with the second sample, consisting of sources within 150 kpc of rich cluster cores, shows that the minimum magnetic energy E_B can be a strongly variable fraction of the inferred accretion energy E_acc, and that it depends on the ambient IGM environment. AGNs inject significant energy as PdV work on the thermal ICM gas, and their magnetic energy, even ignoring the contribution from stellar and starburst outflows, is sufficient to account for that recently found beyond the inner cores of galaxy clusters. Other loss processes in the course of the lobe expansion are considered. We conclude that the aggregate IGM magnetic energy derived purely from galactic black holes since the first epoch of significant galaxy BH formation is sufficiently large that it will have an important influence on the process of both galaxy and visible structure formation on scales up to ~ 1Mpc.

The powerful narrow jets observed to emanate from many compact accreting objects may arise from the twisting of a magnetic field threading a differentially rotating accretion disk which acts to magnetically extract angular momentum and energy from the disk. Two main regimes have been discussed, {\it hydromagnetic outflows}, which have a significant mass flux and have energy and angular momentum carried by both the matter and the electromagnetic field and, Poynting outflows, where the mass flux is negligible and energy and angular momentum are carried predominantly by the electromagnetic field. Here we consider a Keplerian disk initially threaded by a dipole-like magnetic field and we present solutions of the force-free Grad-Shafranov equation for the coronal plasma. We find solutions with Poynting jets where there is a continuous outflow of energy and toroidal magnetic flux from the disk into the external space. This behavior contradicts the commonly accepted ``theorem'' of Solar plasma physics that the motion of the footpoints of a magnetic loop structure leads to a stationary magnetic field configuration with zero power and flux outflows. In addition we discuss recent magnetohydrodynamic (MHD) simulations which establish that quasi-stationary collimated Poynting jets similar to our Grad-Shafranov solutions arise from the inner part of a disk threaded by a dipole-like magnetic field. At the same time we find that there is a steady uncollimated hydromagnetic outflow from the outer part of the disk. The Poynting jets represent a likely model for the jets from active galactic nuclei, microquasars, and gamma ray burst sources.

(abridged) The interpretation of Faraday rotation measure maps of AGNs within galaxy clusters has revealed regions, $\sim 50-100$ kpc, that are populated with large, $\sim 30 \mu$ G magnetic fields. The magnetic energy of these coherent regions is $\sim 10^{59-60}$ ergs, and the total magnetic energy over the whole cluster ($\sim 1$ Mpc across) is expected to be even larger. A sequence of physical processes that are responsible for the production, redistribution and dissipation of these magnetic fields is proposed. These fields are associated with single AGNs within the cluster and therefore with all galaxies during their AGN phase, simply because only the central supermassive black holes ($\sim 10^8 M_{\odot}$) have an accessible energy, $\sim 10^{61}$ ergs. We propose an $\alpha-\Omega$ dynamo process in an accretion disk. The disk rotation naturally provides a large winding number, $\sim 10^{11}$ turns, sufficient to make both large gain and large flux. The helicity of the dynamo can be generated by the differential plume rotation derived from star-disk collisions. This helicity generation process has been demonstrated in the laboratory and the dynamo gain was simulated numerically. A liquid sodium analog of the dynamo is being built. Speculations are that the back reaction of the saturated dynamo will lead to the formation of a force-free magnetic helix, which will carry the energy and flux of the dynamo away from the accretion disk and redistribute the field within the clusters and galaxy walls. The magnetic reconnection of a small fraction of this energy logically is the source of the AGN luminosity, and the remainder of the field energy should then dominate the free energy of the present-day universe.

(abridged) We study the nonlinear evolution of the Rossby wave instability in thin disks using global 2D hydrodynamic simulations. The key questions we are addressing in this paper are: (1) What happens when the instability becomes nonlinear? Specifically, does it lead to vortex formation? (2) What is the detailed behavior of a vortex? (3) Can the instability sustain itself and can the vortex last a long time? Among various initial equilibria that we have examined, we generally find that there are three stages of the disk evolution: (1) The exponential growth of the initial small amplitude perturbations. This is in excellent agreement with the linear theory; (2) The production of large scale vortices and their interactions with the background flow, including shocks. Significant accretion is observed due to these vortices. (3) The coupling of Rossby waves/vortices with global spiral waves, which facilitates further accretion throughout the whole disk. Even after more than 20 revolutions at the radius of vortices, we find that the disk maintains a state that is populated with vortices, shocks, spiral waves/shocks, all of which transport angular momentum outwards. We elucidate the physics at each stage and show that there is an efficient outward angular momentum transport in stages (2) and (3) over most parts of the disk, with an $\alpha$ in the range from $10^{-4}$ to $10^{-2}$. We show why such vortices prove to be almost ideal ``units'' in transporting angular momentum outwards, namely by positively correlating the radial and azimuthal velocity components. We find some special cases in which entropy can remain the same while angular momentum is transported. We conclude that Rossby wave/vortex instability is an efficient, purely hydrodynamic mechanism for angular momentum transport in thin disks.

Recent rotation measure observations of a dozen or so galaxy clusters have revealed a surprisingly large amount of magnetic fields, whose estimated energy and flux are, on average, $\sim 10^{58}$ ergs and $\sim 10^{41}$ G cm$^2$, respectively. These quantities are so much larger than any coherent sums of individual galaxies within the cluster that an efficient galactic dynamo is required. We associate these fields with single AGNs within the cluster and therefore with all galaxies during their AGN phase. Only the central, massive black hole (BH) has the necessary binding energy, $\sim 10^{61}$ ergs. Only the accretion disk during the BH formation has the winding number, $\sim 10^{11}$ turns, necessary to make the gain and magnetic flux. We present a model of the BH accretion disk dynamo that might create these magnetic fields, where the helicity of the $\alpha - \Omega$ dynamo is driven by star-disk collisions. The back reaction of the saturated dynamo forms a force-free field helix that carries the energy and flux of the dynamo and redistributes them within the clusters.

The mechanism for efficient and coherent angular momentum transport remains one of the unsolved puzzles in astrophysics despite the enormous efforts that have been made. We suggest that important new insight could be gained in this problem through an experimental and theoretical study of a laboratory device (Ranque-Hilsch tube) that displays a similar enhanced angular momentum transfer which cannot be explained by a simple turbulent model. There is already good experimental evidence to suggest that the cause of this enhancement is the formation of aligned vortices that swirl around the symmetry axes very much like virtual paddle blades.

In earlier work we identified a global, non-axisymmetric instability associated with the presence of an extreme in the radial profile of the key function ${\cal L}(r) \equiv (\Sigma \Omega/\kappa^2) S^{2/\Gamma}$ in a thin, inviscid, nonmagnetized accretion disk. Here, $\Sigma(r)$ is the surface mass density of the disk, $\Omega(r)$ the angular rotation rate, $S(r)$ the specific entropy, $\Gamma$ the adiabatic index, and $\kappa(r)$ the radial epicyclic frequency. The dispersion relation of the instability was shown to be similar to that of Rossby waves in planetary atmospheres. In this paper, we present the detailed linear theory of this Rossby wave instability and show that it exists for a wider range of conditions, specifically, for the case where there is a ``jump'' over some range of $r$ in $\Sigma(r)$ or in the pressure $P(r)$. We elucidate the physical mechanism of this instability and its dependence on various parameters, including the magnitude of the ``bump'' or ``jump,'' the azimuthal mode number, and the sound speed in the disk. We find large parameter range where the disk is stable to axisymmetric perturbations, but unstable to the non-axisymmetric Rossby waves. We find that growth rates of the Rossby wave instability can be high, $\sim 0.2 \Omega_{\rm K}$ for relative small ``jumps'' or ``bumps''. We discuss possible conditions which can lead to this instability and the consequences of the instability.

The lack of our knowledge on how angular momentum is transported in accretion disks around black holes has prevented us from fully understanding their high energy emissions. We briefly highlight some theoretical models, emphasizing the energy flow and electron energization processes. More questions and uncertainties are raised from a plasma physics point of view.

The accretion induced collapse (AIC) of a white dwarf into a neutron star has been invoked to explain gamma-ray bursts, Type Ia supernovae, and a number of problematic neutron star populations and specific binary systems. The ejecta from this collapse has also been claimed as a source of r-process nucleosynthesis. So far, most AIC studies have focussed on determining the event rates from binary evolution models and less attention has been directed toward understanding the collapse itself. However, the collapse of a white dwarf into a neutron star is followed by the ejection of rare neutron-rich isotopes. The observed abundance of these chemical elements may set a more reliable limit on the rate at which AICs have taken place over the history of the galaxy. In this paper, we present a thorough study of the collapse of a massive white dwarf in 1- and 2-dimensions and determine the amount and composition of the ejected material. We discuss the importance of the input physics (equation of state, neutrino transport, rotation) in determining these quantities. These simulations affirm that AICs are too baryon rich to produce gamm-ray bursts and do not eject enough nickel to explain Type Ia supernovae (with the possible exception of a small subclass of extremely low-luminosity Type Ias). Although nucleosynthesis constraints limit the number of neutron stars formed via AICs to <0.1% of the total galactic neutron star population, AICs remain a viable scenario for forming systems of neutron stars which are difficult to explain with Type II core-collapse supernovae.

We find a linear instability of non-axisymmetric Rossby waves in a thin non-magnetized Keplerian disk when there is a local maximum in the radial profile of a key function ${\cal L}(r) \equiv {\cal F}(r) S^{2/\Gamma}(r)$, where ${\cal F}^{-1} = \hat {\bf z}\cdot ({\bf \nabla}\times {\bf v}) /\Sigma$ is the potential vorticity, $S = P/\Sigma^\Gamma$ is the entropy, $\Sigma$ is the surface mass density, $P$ is the vertically integrated pressure, and $\Gamma$ is the adiabatic index. We consider in detail the special case where there is a local maximum in the disk entropy profile $S(r)$. This maximum acts to trap the waves in its vicinity if its height to width ratio ${\rm max}(S)/\Delta r$ is larger than a threshold value. The pressure gradient derived from this entropy variation provides the restoring force for the wave growth. We show that the trapped waves act to transport angular momentum outward. A plausible way to produce an entropy variation is when an accretion disk is starting from negligible mass and temperature, therefore negligible entropy. As mass accumulates by either tidal torquing, magnetic torquing, or Roche-lobe overflow, confinement of heat will lead to an entropy maximum at the outer boundary of the disk. Possible nonlinear developments from this instability include the formation of Rossby vortices and the formation of spiral shocks. What remains to be determined from hydrodynamic simulations is whether or not Rossby wave packets (or vortices) ``hold together'' as they propagate radially inward.

Neutron stars formed in Type II supernovae are likely to be initially obscured by late-time fallback. Although much of the late-time fallback is quickly accreted via neutrino cooling, some material remains on the neutron star, forming an atmosphere which slowly accretes through photon emission. In this paper, we derive structure equations of the fallback atmosphere and present results of one-dimensional simulations of that fallback. The atmosphere remaining after neutrino cooling becomes unimportant (less than the Compton Eddington limit) is only a fraction of the total mass accreted (10^-8 of the accreted mass or 10^-9 solar masses.) Recombined iron dominates the opacity in the outer regions leading to an opacity 1000-10,000 times higher than that of electron scattering alone. The resultant photon emission of the remnant atmosphere is limited to 1/1000th the Compton Eddington Luminosity. The late-time evolution of this system leads to the formation of a photon-driven wind from the accretion of the inner portion of the atmosphere, leaving, for most cases, a bare neutron star on timescales shorter than a year. The degenerate remnant of 1987a may not be a black hole. Instead, the fallback material may have already accreted or blown off in the accretion-driven wind. If the neutron star has either a low magnetic field or a low rotational spin frequency, we would not expect to see the neutron star remnant of 1987a.

Condensed Abstract: We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 Msun and 25 Msun) in two dimensions. Our calculations begin at the onset of core collapse and stop several 100 ms after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. (...) Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. (...) In essence, a Carnot cycle is established in which convection allows out-of-equilibrium heat transfer mediated by neutrinos to drive low entropy matter to higher entropy and therefore extracts mechanical energy from the heat generated by gravitational collapse. We argue that supernova explosions are nearly guaranteed and self-regulated by the high efficiency of the thermodynamic engine. (...) Convection continues to accumulate energy exterior to the neutron star until a successful explosion has occurred. At this time, the envelope is expelled and therefore uncoupled from the heat source (the neutron star) and the energy input ceases. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae which is effective under a wide range of physical parameters.