In this first paper from forthcoming series of works devoted to radio image of relativistic jets from active galactic nuclei the role of internal structure of a flow is discussed. We determine the radial profiles of all physical values for reasonable Michel magnetization parameter $\sigma_{\rm M}$ and ambient pressure $P_{\rm ext}$. Maps of Doppler boosting factor $\delta$ and observed directions of linear polarization of synchrotron emission are also constructed.

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.

We report experimental observation of an instability in a Couette-Taylor flow of a polymer fluid in a thin gap between two coaxially rotating cylinders in a regime where their angular velocity decreases with the radius while the specific angular momentum increases with the radius. In the considered regime, neither the inertial Rayleigh instability nor the purely elastic instability are possible. We propose that the observed "elasto-rotational" instability is an analog of the magnetorotational instability which plays a fundamental role in astrophysical Keplerian accretion disks.

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)

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

Most of energy carried by relativistic AGN jets remains undetected until hundreds of kiloparsecs where interaction with intergalactic medium produces hot spots. The jet's hidden energy is only partially dissipated at smaller scales, from parsecs to kiloparsecs. Several media may play the role of the ``prime mover'': ions, pairs or large scale magnetic fields. Analyzing VLBI polarization structures of relativistic parsec scale jets we conclude that large-scale magnetic fields can explain the salient polarization properties of parsec-scale AGN jets. This implies that large-scale magnetic fields carry a non-negligible fraction of jet luminosity. We also discuss the possibility that relativistic AGN jets may be electromagnetically (Poynting flux) dominated. In this case, dissipation of the toroidal magnetic field (and not fluid shocks) may be responsible for particle acceleration.

(Abridged) We consider the polarization properties of optically thin synchrotron radiation emitted by relativistically moving electron--positron jets carrying large-scale helical magnetic fields. In our model, the jet is cylindrical, and the emitting plasma moves parallel to the jet axis with a characteristic Lorentz factor $\Gamma$. We draw attention to the strong influence that the bulk relativistic motion of the emitting relativistic particles has on the observed polarization. We conclude that large-scale magnetic fields can explain the salient polarization properties of parsec-scale AGN jets. Since the typical degrees of polarization are $\leq 15%$, the emitting parts of the jets must have comparable rest-frame toroidal and poloidal fields. In this case, most relativistic jets are strongly dominated by the toroidal magnetic field component in the observer's frame, $B_\phi/B_z \sim \Gamma$. We also discuss the possibility that relativistic AGN jets may be electromagnetically (Poynting flux) dominated. In this case, dissipation of the toroidal magnetic field (and not fluid shocks) may be responsible for particle acceleration.

Observations by the {\RHESSI} satellite of large polarization of the prompt $\gamma$-ray emission from the Gamma Ray Burst GRB021206 \citep{coburn03} imply that the magnetic field coherence scale is larger than the size of the visible emitting region, $\sim R/\Gamma$, where $R$ is the radius of the flow, $\Gamma$ is the associated Lorentz factor. Such fields cannot be generated in a causally disconnected, hydrodynamically dominated outflow. Electromagnetic models of GRBs \citep{lyutikov02}, in which large scale, dynamically dominant, magnetic fields are present in the outflow from the very beginning, provide a natural explanation of this large reported linear polarization. We derive Stokes parameters of synchrotron emission of a relativistically moving plasma with a given magnetic field configuration and calculate the pulse averaged polarization fraction of the emission from a relativistically expanding shell carrying global toroidal magnetic field. For viewing angles larger than $1/\Gamma$ the observed patch of the emitting shell has almost homogeneous magnetic field, producing a large fractional polarization ($56% $ for a power-law energy distribution of relativistic particles $dn/d\epsilon \propto \epsilon^{-3}$). The maximum polarization is smaller than the theoretical upper limit for a stationary plasma in uniform magnetic field due to relativistic kinematic effects.

We consider the propagation and stability of the helical electromagnetic perturbations in a relativistic force-free Pointing flux dominated jet (magnetic helix). Perturbations are found to dissipate energy in the vicinity of the resonant surfaces inside the jet, where the resonance of the perturbation with the Alfven waves occurs. The energy of the perturbations is converted into the energy of accelerated particles, which emit synchrotron radiation. The intensity and polarization of this radiation is modulated with the perturbation. Perturbations with vanishing group velocity are the most important and could be associated with the bright knots in parsec scale jets.

All experiments, which have been proposed so far to model the magnetorotational instability (MRI) in the laboratory, involve a Couette flow of liquid metals in a rotating annulus. All liquid metals have small magnetic Prandtl numbers, Pm, of about 10^{-6} (the ratio of kinematic viscosity to magnetic diffusivity). With plasmas both large and small Pm are achievable by varying the temperature and the density of plasma. Compressibility and fast rotation of the plasma result in radial stratification of the equilibrium plasma density. Evolution of perturbations in radially stratified viscous and resistive plasma permeated by an axial uniform magnetic field is considered. The differential rotation of the plasma is induced by the ExB drift in applied radial electric field. Global unstable eigenmodes are calculated by our newly developed matrix code. The plasma is shown to be MRI unstable for parameters easily achievable in experimental setup.

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).

We calculate the optically thin synchrotron emission of fast electrons and positrons in a spiral stationary magnetic field and a radial electric field of a rotating relativistic strongly magnetized force-free jet consisting of electron-positron pair plasma. The magnetic field has a helical structure with a uniform axial component and a toroidal component that is maximal inside the jet and decreasing to zero towards the boundary of the jet. Doppler boosting and swing of the polarization angle of synchrotron emission due to the relativistic motion of the emitting volume are calculated. The distribution of the plasma velocity in the jet is consistent with the electromagnetic field structure. Two spatial distributions of fast particles are considered: uniform, and concentrated in the vicinity of the Alfven resonance surface. The latter distribution corresponds to the regular acceleration by an electromagnetic wave in the vicinity of its Alfven resonance surface inside the jet. The polarization properties of the radiation have been obtained and compared with the existing VLBI polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our results give a natural explanation of the observed bimodality in the alignment between the electric field vector of the polarized radiation and the projection of the jet axis on the plane of the sky. We interpret the motion of bright knots as a phase velocity of standing spiral eigenmodes of electromagnetic perturbations in a cylindrical jet. The degree of polarization and the velocity of the observed proper motion of bright knots depend upon the angular rotational velocity of the jet. The observed polarizations and velocities of knots indicate that the magnetic field lines are bent in the direction opposite to the direction of the jet rotation.

We develop a theory of particle acceleration inside relativistic rotating electron-positron force-free jets with spiral magnetic fields. We considered perturbation of the stationary magnetic field structure and found that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is fulfilled, i.e. where the local Alfven speed is equal to the phase speed of an eigenmode. The acceleration mechanism is regular acceleration by the electric field of the electromagnetic wave, which is the eigenmode of the force-free cylindrical jet: particles drift out of the region of the large wave amplitude near the Alfven resonant surface and gain energy. Acceleration in the strong electric field near the Alfven resonance and synchrotron losses combine to form a power-law energy spectrum of ultra-relativistic electrons and positrons with index between 2 and 3 depending upon the initial energy of the injected particles. The power law distribution ranges from 10 MeV to 1000 MeV.

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 new kind of dynamo utilizing flowing laboratory plasmas has been identified. Conversion of plasma kinetic energy to magnetic energy is verified numerically by kinematic dynamo simulations for magnetic Reynolds numbers above 210. As opposed to intrinsically-turbulent liquid-sodium dynamos, the proposed plasma dynamos correspond to laminar flow topology. Modest plasma parameters, 1-20 eV temperatures, 10^{19}-10^{20} m^{-3} densities in 0.3-1.0 m scale-lengths driven by velocities on the order of the Alfven Critical Ionization Velocity (CIV), self-consistently satisfy the conditions needed for the magnetic field amplication. Growth rates for the plasma dynamos are obtained numerically with different geometry and magnetic Reynolds numbers. Magnetic-field-free coaxial plasma guns can be used to sustain the plasma flow and the dynamo.

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.

(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.

The observed iron K-alpha fluorescence lines in Seyfert-1 galaxies provide strong evidence for an accretion disk near a supermassive black hole as a source of the emission. Here we present an analysis of the geometrical and kinematic properties of the disk based on the extreme frequency shifts of a line profile as determined by measurable flux in both the red and blue wings. The edges of the line are insensitive to the distribution of the X-ray flux over the disk, and hence provide a robust alternative to profile fitting of disk parameters. Our approach yields new, strong bounds on the inclination angle of the disk and the location of the emitting region. We apply our method to interpret observational data from MCG-6-30-15 and find that the commonly assumed inclination 30 deg for the accretion disk in MCG-6-30-15 is inconsistent with the position of the blue edge of the line at a 3 sigma level. A thick turbulent disk model or the presence of highly ionized iron may reconcile the bounds on inclination from the line edges with the full line profile fits based on simple, geometrically thin disk models. The bounds on the innermost radius of disk emission indicate that the black hole in MCG-6-30-15 is rotating faster than 30 % of theoretical maximum. When applied to data from NGC 4151, our method gives bounds on the inclination angle of the X-ray emitting inner disk of 50\pm 10 deg, consistent with the presence of an ionization cone grazing the disk as proposed by Pedlar et al. (1993). The frequency extrema analysis also provides limits to the innermost disk radius in another Seyfert-1 galaxy, NGC 3516.

The observed iron K-alpha fluorescence lines in Seyfert-1 galaxies provide strong evidence for an accretion disk near a supermassive black hole as a source of the line emission. These lines serve as powerful probes for examining the structure of inner regions of accretion disks. Previous studies of line emission have considered geometrically thin disks only, where the gas moves along geodesics in the equatorial plane of a black hole. Here we extend this work to consider effects on line profiles from finite disk thickness, radial accretion flow and turbulence. We adopt the Novikov and Thorne (1973) solution, and find that within this framework, turbulent broadening is the dominant new effect. The most prominent change in the skewed, double-horned line profiles is a substantial reduction in the maximum flux at both red and blue peaks. The effect is most pronounced when the inclination angle is large, and when the accretion rate is high. Thus, the effects discussed here may be important for future detailed modeling of high quality observational data.

The observed iron K-alpha fluorescence lines in Seyfert galaxies provide strong evidence for an accretion disk near a supermassive black hole as a source of the line emission. Previous studies of line emission have considered only geometrically thin disks, where the gas moves along geodesics in the equatorial plane of a black hole. Here we extend this work to include effects on line profiles from finite disk thickness, radial accretion flow and turbulence. We adopt the Novikov-Thorne solution, and find that within this framework, turbulent broadening is the most significant effect. The most prominent changes in the skewed, double-horned line profiles is a substantial reduction in the maximum flux at both red and blue peaks. We show that at the present level of signal-to-noise in X-ray spectra, proper treatment of the actual structure of the accretion disk can change estimates of the inclination angle of the disk. Thus these effects will be important for future detailed modeling of high quality observational data.

We developed theory of particle acceleration inside the relativistic rotating electron-positron force-free jet with spiral magnetic field. We considered perturbation of stationary magnetic field structure and found that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is fulfilled. Acceleration process and synchrotron losses combined together form power law energy spectrum of ultrarelativistic electrons and positrons with index between 2 and 3 depending upon initial energy of injected particles. Synchrotron emission of these electrons and positrons in spiral magnetic field of rotating force-free jet has been calculated. Polarization properties of the radiation has been obtained and compared with existing VLBI polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our results give explanation of observed bimodality in alignment between electric field vector and jet axis. Degree of polarization and velocity of observed proper motion of bright knots depend upon angular rotational velocity of the jet. Thus, comparing them to each other, we can estimate angular rotational velocity in jets. We obtain that the fact, that generally in BL Lac objects electric field vector is oriented parallel to the jet axis while in quasars perpendicular to the jet axis, may be due to intrinsically larger angular rotational velocity and large winding of magnetic field in BL Lac jets than in quasar jets.

We investigate the linear stability of a hydrodynamic relativistic flow of magnetized plasma in the force--free approximation. We considered the case of light cylindrical jet in cold and dense environment, so the jet boundary remains at rest. Continuous and discrete spectra of frequencies are investigated analytically. An infinite sequence of eigenfrequencies is found near the edge of Alfv\'en continuum. Numerical calculations showed that modes are stable and have attenuation increment $\gamma$ small. The dispersion curves $\omega =\omega (k_\parallel )$ have a minimum for $k_{{\parallel}_0}\simeq 1/R$ ($R$ is the jet radius ). This results in accumulation of perturbations inside the jet with wavelength of the order of the jet radius. The wave crests of the perturbation pattern formed in such a way move along the jet with the velocity exceeding light speed. If one has relativistic electrons emitting synchrotron radiation inside the jet, than this pattern will be visible. This provide us with the new type of superluminal source. If the jet is oriented close to the line of sight, than the observer will see knots moving backward to the core.

Bondi type hydrodynamic accretion of the surrounding matter onto Kerr black hole with an arbitrary rotational parameter is considered. The effects of viscosity, thermal conductivity and interaction with radiation field are neglected. The black hole is supposed to be at rest with respect to matter at infinity. The flow is adiabatic and has no angular momentum. The fact that usually in astrophysics substance far from the black hole has nonrelativistic temperature introduces small parameter to the problem and allows to search for the solution as a perturbation to the accretion of a cold, that is dust--like, matter. However, far from the black hole on the scales of order of the radius of the sonic surface the expansion must be performed with respect to Bondi spherically symmetrical solution for the accretion on a Newtonian gravitating centre. The equations thus obtained are solved analytically. The conditions of the regularity of the solution at the sonic surface and at infinity allow to specify unique solution, to find the shape of the sonic surface and to determine the corrections to Bondi accretion rate.