The central radio source in M87 provides the best opportunity to study jet formation because it has a large angular size for the gravitational radius of the black hole and has a bright jet that is well resolved by VLBI observations. We present intensive monitoring observations from 2007 and 2008, plus roughly annual observations that span 17 years, all made with the the Very Long Baseline Array at 43 GHz with a resolution of about 30 by 60 Rs. Our high-dynamic-range images clearly show the wide-opening-angle structure and the counter-jet. The jet and counter-jet are nearly symmetric in the inner 1.5 milli-arcseconds (mas; 0.12 pc in projection) with both being edge brightened. Both show deviations from parabolic shape in the form of an initial rapid expansion and subsequent contraction followed by further rapid expansion and, beyond the visible counter-jet, subsequent collimation. Proper motions and counter-jet/jet intensity ratios both indicate acceleration from apparent speeds of $\lesssim 0.5c$ to $\gtrsim 2c$ in the inner about 2 mas and suggest a helical flow. The jet displays a sideways shift with an approximately 8 to 10 year quasi-periodicity. The shift propagates outwards non-ballistically and significantly more slowly than the flow speed revealed by the fastest moving components. Polarization data show a systematic structure with magnetic field vectors that suggest a toroidal field close to the core.

In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global jets containing helical magnetic fields.

We study the interaction of relativistic jets with their environment, using 3-dimensional relativistic particle-in-cell simulations for two cases of jet composition: (i) electron-proton ($e^{-}-p^{+}$) and (ii) electron-positron ($e^{\pm}$) plasmas containing helical magnetic fields. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability and the Mushroom instability. We have found that these kinetic instabilities are suppressed and new types of instabilities can grow. For the $e^{-}-p^{+}$ jet, a recollimation-like instability occurs and jet electrons are strongly perturbed, whereas for the $e^{\pm}$ jet, a recollimation-like instability occurs at early times followed by kinetic instability and the general structure is similar to a simulation without a helical magnetic field. We plan to perform further simulations using much larger systems to confirm these new findings.

M87 is the best source in which to study a jet at high resolution in gravitational units because it has a very high mass black hole and is nearby. The angular size of the black hole is second only to Sgr A*, which does not have a strong jet. The jet structure is edge brightened with a wide opening angle base and a weak counterjet. We have roughly annual observations for 17 years plus intensive monitoring at three week intervals for a year and five day intervals for 2.5 months made with the Very Long Baseline Array (VLBA) at 43 GHz. The inner jet shows very complex dynamics, with apparent motions both along and across the jet. Speeds from zero to over 2c are seen, with acceleration observed over the first 3 milli-arcseconds. The counterjet decreases in brightness much more rapidly than the main jet, as is expected from relativistic beaming in an accelerating jet oriented near the line-of-sight. Details of the structure and dynamics are discussed. The roughly annual observations show side-to-side motion of the whole jet with a characteristic time scale of about 9 years.

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

Very long baseline interferometry (VLBI) imaging of radio emission from extragalactic jets provides a unique probe of physical mechanisms governing the launching, acceleration, and collimation of relativistic outflows. The two-dimensional structure and kinematics of the jet in M\,87 (NGC\,4486) have been studied by applying the Wavelet-based Image Segmentation and Evaluation (WISE) method to 11 images obtained from multi-epoch Very Long Baseline Array (VLBA) observations made in January-August 2007 at 43 GHz ($\lambda = 7$ mm). The WISE analysis recovers a detailed two-dimensional velocity field in the jet in M\,87 at sub-parsec scales. The observed evolution of the flow velocity with distance from the jet base can be explained in the framework of MHD jet acceleration and Poynting flux conversion. A linear acceleration regime is observed up to $z_{obs} \sim 2$\,mas. The acceleration is reduced at larger scales, which is consistent with saturation of Poynting flux conversion. Stacked cross correlation analysis of the images reveals a pronounced stratification of the flow. The flow consists of a slow, mildly relativistic layer (moving at $\beta \sim 0.5\,c$), associated either with instability pattern speed or an outer wind, and a fast, accelerating stream line (with $\beta \sim 0.92$, corresponding to a bulk Lorentz factor $\gamma \sim 2.5$). A systematic difference of the apparent speeds in the northern and southern limbs of the jet is detected, providing evidence for jet rotation. The angular velocity of the magnetic field line associated with this rotation suggests that the jet in M87 is launched in the inner part of the disk, at a distance $r_0 \sim 5\, R_\mathrm{s}$ from the central engine. The combined results of the analysis imply that MHD acceleration and conversion of Poynting flux to kinetic energy play the dominant roles in collimation and acceleration of the flow in M\,87.

One of the key open questions in the study of relativistic jets is their interaction with the environment. Here, we study the initial evolution of both electron-proton and electron-positron relativistic jets, focusing on their lateral interaction with the ambient plasma. We trace the generation and evolution of the toroidal magnetic field generated by both kinetic Kelvin-Helmholtz (kKH) and Mushroom instabilities (MI). This magnetic field collimates the jet. We show that in electron-proton jet, electrons are perpendicularly accelerated with jet collimation. The magnetic polarity switches from the clockwise to anti-clockwise in the middle of jet, as the instabilities weaken. For the electron-positron jet, we find strong mixture of electron-positron with the ambient plasma, that results in the creation of a bow shock. Merger of magnetic field current filaments generate density bumps which initiate a forward shock. The strong mixing between jet and ambient particles prevents full development of the jet on the studied scale. Our results therefore provide a direct evidence for both jet collimation and particle acceleration in the created bow shock. Differences in the magnetic field structures generated by electron-proton and electron-positron jets may contribute to observable differences in the polarized properties of emission by electrons.

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants.

We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas we find generation of strong large-scale DC currents and magnetic fields which extend over the entire shear-surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates transverse structure similar to that produced by the Weibel instability.

Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code Ratpenat. We study the evolution of Fanaroff-Riley Class I (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FRI jets within the first kiloparsec. However, powerful FRI jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow-shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kpc scales is initiated by mass loading in the weaker FRI jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes.

We have investigated the influence of velocity shear and a radial density profile on the spatial development of the current driven kink instability along helically magnetized relativistic jets via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we use a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers growth of the kink instability. If the velocity shear radius is located inside the characteristic radius of the helical magnetic field, a static non-propagating current driven kink is excited as the perturbation propagates down the jet. Temporal growth disrupts the initial flow across the computational grid not too far from the inlet. On the other hand, if the velocity shear radius is outside the characteristic radius of the helical magnetic field, the kink is advected with the flow and grows spatially down the jet. In this case flow is maintained to much larger distances from the inlet. The effect of different radial density profiles is more subtle. When the density increases with radius, the kink appears to saturate by the end of the simulation without apparent disruption of the helical twist. This behavior suggests that relativistic jets consisting of a tenuous spine surrounded by a denser medium with a velocity shear radius outside the radius of maximum toroidal magnetic field have a relatively stable configuration.

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($\gamma_{\rm j} = 1.5$) is larger than for a relativistic jet case ($\gamma_{\rm j} = 15$).

We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants.

The remarkable stability of extragalactic jets is surprising, given the reasonable possibility of the growth of instabilities. In addition, much work in the literature has invoked this possibility in order to explain observed jet structures and obtain information from these structures. For example, it was recently shown that the observed helical structures in the jet in S5 0836+710 could be associated with helical pressure waves generated by Kelvin-Helmholtz instability. Our aim is to resolve the arc-second structure of the jet in the quasar S5 0836+710 and confirm the lack of a hot-spot (reverse jet-shock) found by present observing arrays, as this lack implies a loss of jet collimation before interaction with the intergalactic medium. In this work, we use an observation performed in 2008 using EVN and MERLIN. The combined data reduction has provided a complete image of the object at arc-second scales. The lack of a hot-spot in the arc-second radio structure is taken as evidence that the jet losses its collimation between the VLBI region and the region of interaction with the ambient medium. This result, together with the previous identification of the helical structures in the jet with helical pressure waves that grow in amplitude with distance, allow us to conclude that the jet is probably disrupted by the growth of Kelvin-Helmholtz instability. This observational evidence confirms that the physical parameters of jets can be extracted using the assumption that instability is present in jets and can be the reason for many observed structures. Interestingly, the observed jet is classified as a FRII object in terms of its luminosity, but its large-scale morphology does not correspond to this classification. The implications of this fact are discussed.

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. We find that, in accord with linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outwards from the axis, the initial small perturbations grow strongly, and if multiple wavelengths are excited non-linear interaction eventually disrupts the initial cylindrical configuration. When the profile of the poloidal field is shallow, the instability develops slowly and eventually saturates. We briefly discuss implications of our findings for Poynting dominated jets.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

Helical structures are common in extragalactic jets. They are usually attributed in the literature to periodical phenomena in the source (e.g., precession). In this work, we use VLBI data of the radio-jet in the quasar S5 0836+710 and hypothesize that the ridge-line of helical jets like this corresponds to a pressure maximum in the jet and assume that the helically twisted pressure maximum is the result of a helical wave pattern. For our study, we use observations of the jet in S5 0836+710 at different frequencies and epochs. The results show that the structures observed are physical and not generated artificially by the observing arrays. Our hypothesis that the observed intensity ridge-line can correspond to a helically twisted pressure maximum is confirmed by our observational tests. This interpretation allows us to explain jet misalignment between parsec and kiloparsec scales when the viewing angle is small, and also brings us to the conclusion that high-frequency observations may show only a small region of the jet flow concentrated around the maximum pressure ridge-line observed at low frequencies. Our work provides a potential explanation for the apparent transversal superluminal speeds observed in several extragalactic jets by means of transversal shift of an apparent core position with time.

lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti-parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Radio and optical images of the M87 jet show bright filaments, twisted into an apparent double helix, extending from HST-1 to knot A. Proper motions within the jet suggest a decelerating jet flow passing through a slower, accelerating wave pattern. We use these observations to develop a mass and energy flux conserving model describing the jet flow and conditions along the jet. We determine the cocoon conditions required if the twisted filaments are the result of the Kelvin-Helmholtz (KH) unstable elliptical mode. We find that the cocoon must be cooler than the jet at HST-1 but must be about as hot as the jet at knot A. Under these conditions we find that the observed filament wavelength is near the elliptical mode maximum growth rate and growth is rapid enough for the filaments to develop and saturate well before HST-1. We generate a pseudo-synchrotron image of a model jet carrying a combination of normal modes of the KH instability. The pseudo-synchrotron image of the jet reveals: (1) that a slow decline in the model jet's surface brightness is still about five times faster than the real jet; (2) that KH produced dual helically twisted filaments can appear qualitatively similar to those on the real jet if any helical perturbation to the jet is very small or nonexistent inside knot A; (3) that the knots in the real jet cannot be associated with the twisted filamentary features and are unlikely to be the result of a KH instability. The existence of the knots in the real jet, the limb brightening of the real jet in the radio, and the slower decline of the surface brightness of the real jet indicate that additional processes --- such as unsteady jet flow and internal particle acceleration --- are occurring within the jet. Disruption of the real jet beyond knot A by KH instability is consistent with the jet and cocoon conditions we find at knot A.

We have investigated the influence of a velocity shear surface on the linear and non-linear development of the CD kink instability of force-free helical magnetic equilibria in 3D. In this study we follow the temporal development within a periodic computational box and concentrate on flows that are sub-Alfvenic on the cylindrical jet's axis. Displacement of the initial force-free helical magnetic field leads to the growth of CD kink instability. We find that helically distorted density structure propagates along the jet with speed and flow structure dependent on the radius of the velocity shear surface relative to the characteristic radius of the helically twisted force-free magnetic field. At small velocity shear surface radius the plasma flows through the kink with minimal kink propagation speed. The kink propagation speed increases as the velocity shear radius increases and the kink becomes more embedded in the plasma flow. A decreasing magnetic pitch profile and faster flow enhance the influence of velocity shear. Simulations show continuous transverse growth in the nonlinear phase of the instability. The growth rate of the CD kink instability and the nonlinear behavior also depend on the velocity shear surface radius and flow speed, and the magnetic pitch radial profile. Larger velocity shear radius leads to slower linear growth, makes a later transition to the nonlinear stage, and with larger maximum amplitude than occur for a static plasma column. However, when the velocity shear radius is much greater than the characteristic radius of the helical magnetic field, linear and non-linear development can be similar to the development of a static plasma column.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system. We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system generates a jitter/synchrotron spectrum.

We have investigated the relaxation of a hydrostatic hot plasma column containing toroidal magnetic field by the Current-Driven (CD) kink instability as a model of pulsar wind nebulae. In our simulations the CD kink instability is excited by a small initial velocity perturbation and develops turbulent structure inside the hot plasma column. We demonstrate that, as envisioned by Begelman, the hoop stress declines and the initial gas pressure excess near the axis decreases. The magnetization parameter \sigma, the ratio of the Poynting to the kinetic energy flux, declines from an initial value of 0.3 to about 0.01 when the CD kink instability saturates. Our simulations demonstrate that axisymmetric models strongly overestimate the elongation of the pulsar wind nebulae. Therefore, the previous requirement for an extremely low pulsar wind magnetization can be abandoned. The observed structure of the pulsar wind nebulae do not contradict the natural assumption that the magnetic energy flux still remains a good fraction of the total energy flux after dissipation of alternating fields.

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

The blazars provide a considerable opportunity to peer into the workings within a few tens of parsecs of the central engine in AGN. This considerable opportunity involves significant challenges as different macroscopic dynamical processes and microscopic physical processes operating at different locations can be responsible for the observed emission. In this proceedings article I review recent theoretical and numerical results relevant to dynamics inside the blazar zone, review the particle acceleration processes capable of producing the high energy particles required by the observed emission, discuss some of the progress made at the microphysical level, and consider what recent TeV and radio observations of M87 can tell us about the blazar zone.

The TeV radio galaxy M87 is the first radio galaxy detected in the TeV regime. The structure of its jet, which is not pointing towards our line of sight, is spatially resolved in X-ray (by Chandra), optical and radio observations. In 2008, the three main Atmospheric Cherenkov Telescope observatories VERITAS, MAGIC and H.E.S.S. coordinated their observations in a joint campaign from January to May with a total observation time of approx. 120 hours. In February, strong and rapid day-scale TeV flares were detected. VLBA monitoring observations during the same period showed that the 43 GHz radio flux density of the unresolved core began to rise at the time of the TeV flares and eventually reached levels above any previously seen with VLBI. New jet components appeared during the flare. The localization accuracy of the TeV instruments of many arcseconds, even for strong sources, is inadequate to constrain the origin of the emission in the inner jets of AGNs. For M87, with a 6 billion solar mass black hole and a distance of 16.7 Mpc, the VLBA resolution instead corresponds to 30 by 60 Schwarzschild radii. This is starting to resolve the jet collimation region. The temporal coincidence of the TeV and radio flares indicates that they are related and provides the first direct evidence that the TeV radiation from this source is produced within a few tens of RS of the radio core, thought to be coincident to within the VLBA resolution with the black hole.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. We have also performed simulations with electron-ion jets. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability for electron-positron jets and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks for both cases. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value for pair plasmas. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system with two different cases for Lorentz factors (15 and 100). We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability.

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context.

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.

We solve the Riemann problem for the deceleration of an arbitrarily magnetized relativistic flow injected into a static unmagnetized medium in one dimension. We find that for the same initial Lorentz factor, the reverse shock becomes progressively weaker with increasing magnetization \sigma (the Poynting-to kinetic energy flux ratio), and the shock becomes a rarefaction wave when \sigma exceeds a critical value, \sigma_c, defined by the balance between the magnetic pressure in the flow and the thermal pressure in the forward shock. In the rarefaction wave regime, we find that the rarefied region is accelerated to a Lorentz factor that is significantly larger than the initial value. This acceleration mechanism is due to the strong magnetic pressure in the flow. We discuss the implications of these results for models of gamma-ray bursts and active galactic nuclei.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless, relativistic shocks, particle (electron, positron, and ion) acceleration is due to plasma waves and their associated instabilities (e.g., the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform, small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting "jitter" radiation from deflected electrons has different properties compared to synchrotron radiation, which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general, and supernova remnants.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the presence of relativistic jets, instabilities such as the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability create collisionless shocks, which are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons in small-scale magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation, a case of diffusive synchrotron radiation, may be important to understand the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

M87 provides the best opportunity to study the base of a jet where it is collimated and accelerated. The size of that region scales with the mass of the black hole, and M87 has the best combination of high mass, proximity to the Earth, and presence of a bright jet. VLBI observations of M87 can probe regions under 100 gravitational radii where theoretical studies suggest that the jet formation and acceleration occurs. A one-year sequence of 43 GHz observations every 3 weeks on the VLBA is being used to study the structure and dynamics in this region. Initial results from that effort are reported here, including the observation of rapid motions - sufficiently rapid that more frequent observations are planned for early 2008. The contribution ends with a discussion of prospects for future VLBI observations of M87 with VSOP2. For VSOP2, a strong recommendation is made that a series of daily observations at 43 GHz be planned.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks e.g. gamma-ray bursts (GRBs) active galactic nuclei (AGNs) and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless relativistic shocks particle (electron, positron and ion) acceleration is due to plasma waves and their associated instabilities (e.g. the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting ``jitter'' radiation from deflected electrons has different properties compared to synchrotron radiation which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general and supernova remnants.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons has different properties than synchrotron radiation which assumes a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

We performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field changes the properties of the shock interface between the tenuous, overpressured jet ($V^z_j$) flowing tangentially to a dense external medium. Magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A ``poloidal'' magnetic field ($B^z$), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and a stronger inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to larger Lorentz factors than those obtained in the pure-hydrodynamic case. Likewise, a strong ``toroidal'' magnetic field ($B^y$), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case. Overall, the acceleration efficiency in the ``poloidal'' case is less than that of the ``toroidal'' case but both geometries still result in higher Lorentz factors than the pure-hydrodynamic case. Thus, the presence and relative orientation of a magnetic field in relativistic jets can significant modify the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

With its high black hole mass, proximity, and bright jet, M87 provides the best prospect for a direct imaging study of the acceleration and collimation region of a jet. Previous VLBI observations have shown an edge brightened structure with a wide opening angle at small scales. An effort to measure component speeds in this region using existing VLBA data at 43 GHz gave tentative results of 0.25c to 0.4c but also indicated that faster sampling is needed. Here we provide a progress report on a project to make a properly sampled movie of motions in the inner jet using the VLBA at 43 GHz. A pilot project during 2006 measured speeds of about 0.6c and was used to set a frame interval of 3 weeks for the movie. The movie observations began in January 2007. Results from the pilot and from the first frame of the movie are presented. The goal of the project is to provide observations of the structure and dynamics of the jet on scales from under 100 to a few hundred Schwarzschild radii that can be compared with expectations from theoretical studies and numerical modeling.

A new general relativistic magnetohydrodynamics (GRMHD) code ``RAISHIN'' used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, c/2, on the KH instability associated with a relativistic, Lorentz factor equal 2.5, jet spine-sheath interaction. In the simulations sound speeds up to c/1.7 and Alfven wave speeds up to 0.56 c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from c/2 sheath speeds is found.

The linearized relativistic magnetohydrodynamic (RMHD) equations describing a uniform axially magnetized cylindrical relativistic jet spine embedded in a uniform axially magnetized relativistically moving sheath are derived. The displacement current is retained in the equations so that effects associated with Alfven wave propagation near light speed can be studied. A dispersion relation for the normal modes is obtained. Analytical solutions for the normal modes in the low and high frequency limits are found and a general stability condition is determined. A trans-Alfvenic and even a super-Alfvenic relativistic jet spine can be stable to velocity shear driven Kelvin-Helmholtz modes. The resonance condition for maximum growth of the normal modes is obtained in the kinetically and magnetically dominated regimes. Numerical solution of the dispersion relation verifies the analytical solutions and is used to study the regime of high sound and Alfven speeds.

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor $\gamma = 2.5$ jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are $\lesssim c/\sqrt 3$ in the weakly magnetized simulations and $\lesssim 0.3c$ in the strongly magnetized simulations. The Alfven wave speed is $\lesssim 0.07c$ in the weakly magnetized simulations and $\lesssim 0.56c$ in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

We investigate jet formation in black-hole systems using 3-D General Relativistic Particle-In-Cell (GRPIC) and 3-D GRMHD simulations. GRPIC simulations, which allow charge separations in a collisionless plasma, do not need to invoke the frozen condition as in GRMHD simulations. 3-D GRPIC simulations show that jets are launched from Kerr black holes as in 3-D GRMHD simulations, but jet formation in the two cases may not be identical. Comparative study of black hole systems with GRPIC and GRMHD simulations with the inclusion of radiate transfer will further clarify the mechanisms that drive the evolution of disk-jet systems.

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using the RAISHIN code. In the numerical simulations a jet with Lorentz factor \gamma=2.5 is precessed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized moving sheath. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

We have performed several simulations of black hole systems (non-rotating, black hole spin parameter a=0.0 and rapidly rotating, a=0.95) with a geometrically thin Keplerian disk using the newly developed RAISHIN code. The simulation results show the formation of jets driven by the Lorentz force and the gas pressure gradient. The jets have mildly relativistic speed (\sim 0.4c). The matter is continuously supplied from the accretion disk and the jet propagates outward until each applicable terminal simulation time (non-rotating: t/\tau_S = 275 and rotating: t/tau_S = 200, \tau_S \equiv r_S/c). It appears that a rotating black hole creates an additional, faster, and more collimated matter-dominated inner outflow (\sim 0.5c) formed and accelerated by the twisted magnetic field resulting from frame-dragging in the black hole ergosphere. This is the first known simulation confirming the formation of an inner magnetically-driven, matter-dominated jet by the frame-dragging effect from a black hole co-rotating with a thin Keplerian disk threaded by a vertical magnetic field. This result indicates that jet kinematic structure depends on black hole rotation and on the initial magnetic field configuration and strength.

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well-known quasar 3C 273. We present here the results from a 3D relativistic hydrodynamics numerical simulation based on the parameters given for the jet by Lobanov & Zensus (2001), and one in which the effects of jet precession and the injection of discrete components have been taken into account. We compare the model with the structures observed in 3C 273 using very long baseline interferometry and constrain the basic properties of the flow. We find growing perturbation modes in the simulation with similar wavelengths to those observed, but with a different set of wave speeds and mode identification. If the observed longest helical structure is produced by the precession of the flow, longer precession periods should be expected. Our results show that some of the observed structures could be explained by growing Kelvin-Helmholtz instabilities in a slow moving region of the jet. However, we point towards possible errors in the mode identification that show the need of more complete linear analysis in order to interpret the observations. We conclude that, with the given viewing angle, superluminal components and jet precession cannot explain the observed structures.

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks may be responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated long-term particle acceleration associated with relativistic electron-ion or electron-positron jet fronts propagating into an unmagnetized ambient electron-ion or electron-positron plasma. These simulations have been performed with a longer simulation system than our previous simulations in order to investigate the nonlinear stage of the Weibel instability and its particle acceleration mechanism. The current channels generated by the Weibel instability are surrounded by toroidal magnetic fields and radial electric fields. This radial electric field is quasi stationary and accelerates particles which are then deflected by the magnetic field.

Weibel instability created in collisionless shocks is responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-ion jet fronts propagating into an ambient plasma without initial magnetic fields with a longer simulation system in order to investigate nonlinear stage of the Weibel instability and its acceleration mechanism. The current channels generated by the Weibel instability induce the radial electric fields. The z component of the Poynting vector (E x B) become positive in the large region along the jet propagation direction. This leads to the acceleration of jet electrons along the jet. In particular the E x B drift with the large scale current channel generated by the ion Weibel instability accelerate electrons effectively in both parallel and perpendicular directions.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity nearly 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating into an ambient plasma. We find that the growth times of the Weibel instability in electron-positron jets are not affected by the (electron-positron or electron-ion) ambient plasmas. However, the amplitudes of generated local magnetic fields in the electron-ion ambient plasma are significantly larger than those in the electron-positron ambient plasma.

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating into an ambient plasma. We find that the growth times of Weibel instability are proportional to the Lorentz factors of jets. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity $\sim 0.3c$) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units ($\tau_{\rm S} = r_{\rm S}/c$ where $r_{\rm S} \equiv 2GM/c^2$) which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around $r = 3 r_{\rm S}$. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfv\'{e}n waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating into an ambient plasma. We find small differences in the results for no ambient and modest ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The small scale magnetic field structure generated by the Weibel instability is appropriate to the generation of ``jitter'' radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. New simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. Furthermore, the non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at the comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. Additionally, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by the Weibel instability scale proportional to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields which contribute to the electron's (positron's) transverse deflection behind the jet head.

The predicted spatial development of helical structures along an expanding jet are used to model observed structures and motions in the 3C 120 jet. Proper motions of helical components associated with the pattern and of other components that move through the pattern provide estimates of flow and helical pattern speeds. Theoretical modeling allows determination of sound speeds as a function of jet viewing angle. Jet morphology argues against viewing angles larger than 12 deg. At a 12 deg viewing angle the most likely scenario involves a decline in jet sound speed from c/3 < a_j < c/1.7 at ~0.5 mas from the core to 0.1c < a_j < 0.25c at ~25 mas from the core accompanied by some acceleration in the jet flow from Lorentz factor < 5 to Lorentz factor > 7. At smaller viewing angles the jet sound speed is lower and at a 6 deg viewing angle the jet sound speed is about a factor 2 less but the flow Lorentz factor is comparable. The decline in radio intensity is on the order of what would be associated with isothermal jet expansion. Knot interknot intensity variations are greater than would be expected from adiabatic compressions associated with the helical twist and we infer the presence of a shock along the leading edge of the helical twist. Our results imply that the macroscopic heating of the expanding jet fluid is less than the microscopic energization of the synchrotron radiating relativistic electrons.

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g.,Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating into an ambient plasma with and without initial magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock at the same simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale (mainly transverse) magnetic fields which contribute to the electron's (positron's) transverse deflection behind the jet head. This small scale magnetic field structure is appropriate to the generation of ``jitter'' radiation from deflected electrons (positrons) as opposed to synchrotron radiation.

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Collapsars may be a source for the long Gamma Ray Bursts (GRBs) in the BATSE catalog. Collapsars may radiate gamma rays anisotropically by beamed jet emission close to the observer's line of sight. These jets must penetrate the radiation-dominated medium of their collapsar progenitor and break through its atmosphere in order to produce a GRB. We present a study of previously published, axisymmetric, collapsar jet simulations. Here we use the linearized, relativistic fluid equations to find the Kelvin-Helmholtz modes that are triggered by recollimation shocks within the jet. The modes will grow as they propagate with the jet. These are of interest since the light curves of GRBs show evidence of a variable flow from the GRB engine. We also evaluate effects of grid scaling in the numerical simulation and show that short wavelength modes are suppressed by grid scaling before the jet breaks out of the Helium shell.

Many jets exhibit twisted helical structures. Where superluminal motions are detected, jet orientation and pattern/flow speed are considerably constrained. In this case modeling efforts can place strong limits on conditions in the jet and in the external environment. This can be done by modeling the spatial development of helical structures which are sensitively dependent on these conditions. Along an expanding jet this sensitivity manifests itself in predictable changes in pattern speed and observed wavelength. In general, twists of low frequency relative to the local resonant frequency are advected along the expanding jet into a region in which the twist frequency is high relative to the local resonant frequency. The wave speed can be very different in these two frequency regimes. Potential effects include helical twists with a nearly constant apparent wavelength, an apparent wavelength scaling approximately with the jet radius for up to two orders of magnitude of jet expansion, or multiple twist wavelengths with vastly different intrinsic scale and vastly different wave speeds that give rise to similar observed twist wavelengths but with very different observed motion. In this paper I illustrate the basic intrinsic and observed behavior of these structures and show how to place constraints on jet conditions in superluminal jets using the apparent structures and motions in the inner 3C 120 jet.

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation (Medvedev 2000) from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

The wealth of high quality data now available on the M87 jet inspired us to carry out a detailed analysis of the plasma physical conditions in the jet. In a companion paper (Lobanov, Hardee & Eilek, this proceedings) we identify a double-helix structure within the jet, and apply Kelvin-Helmholtz stability analysis to determine the physical state of the jet plasma. In this paper we treat the jet as a test case for in situ particle acceleration. We find that plasma turbulence is likely to exist at levels which can maintain the energy of electrons radiating in the radio to optical range, consistent with the broadband spectrum of the jet.

Relativistic jets in Galactic superluminals and extragalactic AGN may be surrounded by a wind near to the central engine. Theoretical analysis and numerical simulation reveal considerable stabilization of relativistic jet flow by a wind to helical and higher order asymmetric modes of jet distortion. When velocities are measured in the source (inlet) frame, reduction in the absolute velocity difference between jet and wind, $\Delta v = v_{jet} - v_{wnd}$, provides stabilization in addition to stabilization provided by a high jet Lorentz factor, but a high Lorentz factor wind is not needed to stabilize a high Lorentz factor jet. However, the fundamental pinch mode is not similarly affected and knots with spacing a few times the jet radius are anticipated to develop in such flows. Thus, we identify a mechanism that can suppress large scale asymmetric structures while allowing axisymmetric structures to develop. Relativistic jets surrounded by outflowing winds will be more stable than if surrounded by a stationary or backflowing external medium. Knotty structures along a straight jet like that in 3C 175 could be triggered by pinching of an initially low Mach number jet surrounded by a suitable wind. As the jet enters the radio lobe, suppression of any surrounding outflow or backflow associated with the high pressure lobe triggers exponential growth of helical twisting.

Two three-dimensional magnetohydrodynamical simulations of strongly magnetized conical jets, one with a poloidal and one with a helical magnetic field, have been performed. In the poloidal simulation a significant sheath (wind) of magnetized moving material developed and partially stabilized the jet to helical twisting. The fundamental pinch mode was not similarly affected and emission knots developed in the poloidal simulation. Thus, astrophysical jets surrounded by outflowing winds could develop knotty structures along a straight jet triggered by pinching. Where helical twisting dominated the dynamics, magnetic field orientation along the line-of-sight could be organized by the toroidal flow field accompanying helical twisting. On astrophysical jets such structure could lead to a reversal of the direction of Faraday rotation in adjacent zones along a jet. Theoretical analysis showed that the different dynamical behavior of the two simulations could be entirely understood as a result of dependence on the velocity shear between jet and wind which must exceed a surface Alfven speed before the jet becomes unstable to helical and higher order modes of jet distortion.

Results are presented from long term VLBI monitoring of the parsec-scale radio jet in 3C120, primarily at 5 and 1.7 GHz. Numerous superluminal features are followed, some over large fractional distances from the core. The speeds of the features are not all the same, but vary by less than a factor of 2. The superluminal features extend to core distances of over 0.2 arcseconds and the region of slowing has not been found. The 1.7 GHz images show evidence for stationary features in addition to the superluminal features. The structure of the jet in the vicinity of the most likely stationary feature is suggestive of a helical pattern seen in projection. The deprojected wavelength of the pattern is very long relative to the jet radius, unlike the situation in sources such as M87. If the 3C120 jet does contain a slowly-moving, helical structure, then theory suggests that the jet resides in a relatively cool medium, not in a relativistically hot cocoon or lobe.

Three dimensional numerical simulations of the response of a Lorentz factor 2.5 relativistic jet to precession at three different frequencies have been performed. Low, moderate and high precession frequencies have been chosen relative to the maximally unstable frequency predicted by a Kelvin-Helmholtz stability analysis. Transverse motion and velocity decreases as the precession frequency increases. Although helical displacement of the jet decreases in amplitude as the precession frequency increases, a helical shock is generated in the medium external to the jet at all precession frequencies. Complex pressure and velocity structure inside the jet is shown to be produced by a combination of the helical surface and first body modes predicted by a normal mode analysis of the relativistic hydrodynamic equations. The surface and first body mode have different wave speed and wavelength, are launched in phase by the periodic precession, and exhibit beat patterns in synthetic emission images. Wave (pattern) speeds range from $0.41c$ to $0.86c$ but beat patterns remain stationary. Thus, we find a mechanism that can produce differentially moving and stationary features in the jet.

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component, instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness, and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component.

The tailed radio galaxies that have been called ``Type I'' are not a uniform set. To study their dynamics, we have used the Ledlow-Owen data set, which provides a new sample of 250 radio galaxies in nearby Abell clusters. These sources divide into two clear categories based on their radio morphology. Type A sources (``straight'') contain nearly straight jets which are embedded in outer radio lobe. Type B sources (``tailed'') have a well-collimated jet flow which undergoes a sudden transition, at an inner hot spot, to a less collimated flow which continues on and forms a radio tail. We have not found any separation of these classes in terms of radio power, radio flux size, galaxy power or external gas density. We propose the difference is due to the development, or not, of a disruptive flow instability, such as Kelvin-Helmholtz, and the saturation of the instability when it develops.

The effect of optically thin radiative cooling on the Kelvin-Helmholtz instability of three dimensional jets is investigated via linear stability theory and nonlinear hydrodynamical simulation. Two different cooling functions are considered: radiative cooling is found to have a significant effect on the stability of the jet in each case. The wavelengths and growth rates of unstable modes in the numerical simulations are found to be in good agreement with theoretical predictions. Disruption of the jet is found to be sensitive to the precessional frequency at the origin with lower frequencies leading to more rapid disruption. Strong nonlinear effects are observed as the result of the large number of normal modes in three dimensions which provide rich mode-mode interactions. These mode-mode interactions provide new mechanisms for the formation of knots in the flows. Significant structural features found in the numerical simulations appear similar to structures observed on protostellar jets.

We extend investigations of mass entrainment by jets, which previously have focused on cylindrical supermagnetosonic jets and expanding trans-Alfv\'enic jets, to a set of expanding supermagnetosonic jets. We precess these jets at the origin to excite the helical mode of the Kelvin-Helmholtz (or KH) instability, in order to compare the results with predictions from linear stability analysis. We analyze this simulation set for the spatial development of magnetized mass, which we interpret as jet plus entrained, initially unmagnetized external mass. As with the previous simulation sets, we find that the growth of magnetized mass is associated with the growth of the KH instability through linear, nonlinear, and saturated stages and with the expansion of magnetized material in simulated observations of the jet. From comparison of measured wavelengths and wave speeds with the predictions from linear stability analysis, we see evidence that the KH instability is the primary cause for mass entrainment in these simulations, and that the expansion reduces the rate of mass entrainment. This reduced rate can be observed as a somewhat greater distance between the two transition points separating the three stages of expansion.

We present the results of a series of axisymmetric time-dependent magnetohydrodynamic (MHD) simulations of the propagation of cooling, overdense jets, motivated by the properties of outflows associated with young stellar objects. A variety of initial field strengths and configurations are explored for both steady and time-variable (pulsed) jets. Even apparently weak magnetic fields with strengths B < 60 micro-G in the pre-shocked jet beam can have a significant effect on the dynamics, for example by altering the density, width, and fragmentation of thin shells formed by cooling gas. A linear analysis predicts that axisymmetric pinch modes of the MHD Kelvin-Helmholtz instability should grow only slowly for the highly supermagnetosonic jets studied here; we find no evidence for them in our simulations. Some of our models appear unstable to current-driven pinch modes, however the resulting pressure and density variations induced in the jet beam are not large, making this mechanism an unlikely source of emission knots in the jet beam. In the case of pulsed jets, radial hoop stresses confine shocked jet material in the pulses to the axis, resulting in a higher density in the pulses in comparison to purely hydrodynamic models.

The properties of wave-like helically twisted normal mode structures on steady relativistic jets are summarized. Wave speeds are a function of the wavelength and less than the jet speed. However, normal mode interference can lead to both stationary and superluminal phase effects. A maximum pressure fluctuation criterion suggested by numerical simulations of axisymmetric relativistic jets is used to find the maximum asymmetric jet distortions and velocity fluctuations. Cyclic transverse velocity fluctuation can lead to variation in the flow direction on the order of the relativistic beaming angle. Resulting variation in the Doppler boost factor can lead to significant brightness asymmetries as helical structures twist around the jet beam. Growth of these structures is reduced as the jet density, Lorentz factor or Mach number are increased. Maximum jet distortion is reduced as the Lorentz factor increases and this suggests a reduction in mass entrainment or other non-linear disruptive processes that influence the morphological development of radio sources.

We investigate the growth of jet plus entrained mass in simulations of supermagnetosonic cylindrical and expanding jets. The entrained mass spatially grows in three stages: from an initially slow spatial rate to a faster rate and finally at a flatter rate. These stages roughly coincide with the similar rates of expansion in simulated radio intensity maps, and also appear related to the growth of the Kelvin-Helmholtz instability through linear, nonlinear, and saturated regimes. In the supermagnetosonic cylindrical jets, we found that a jet with an embedded primarily toroidal magnetic field is more stable than a jet with a primarily axial magnetic field. Also, pressure-matched expanding jets are more stable and entrain less mass than cylindrical jets with equivalent inlet conditions.

Three-dimensional magnetohydrodynamical simulations of strongly magnetized ``light'' conical jets have been performed. An investigation of the transition from sub-Alfv\'enic to super-Alfv\'enic flow has been made for nearly poloidal and for helical magnetic fields. The jets are stable to asymmetric modes of jet distortion provided they are sub-Alfv\'enic over most of their interior but destabilize rapidly when they become on average super-Alfv\'enic. The jets are precessed at the origin and the resulting small amplitude azimuthal motion is communicated down the jet to the Alfv\'en point where it couples to a slowly moving and rapidly growing helical twist. Significant jet rotation can contribute to destabilization via increase in the velocity shear between the jet and the external medium. Destabilization is accompanied by significant mass entrainment and the jets slow down significantly as denser external material is entrained. Synchrotron intensity images satisfactorily reveal large scale helical structures but have trouble distinguishing a large amplitude elliptical jet distortion that appears as an apparent pinching in an intensity image. Smaller scale jet distortions are not clearly revealed in intensity images, largely as a result of the relatively small total pressure variations that accompany destabilization and growing distortions. Fractional polarization is high as a result of the strong ordered magnetic fields except where the intensity image suggests cancellation of polarization vectors by integration through twisted structures.

The amplitude of jet distortions and accompanying pressure and velocity fluctuations resulting from Kelvin-Helmholtz instability of three dimensional relativistic jets are explored. The effect of instability on jets as they accelerate from sub- to super-Alfv\'enic speeds is explored and a relativistic stabilization mechanism for trans-Alfv\'enic jets is proposed. The level to which asymmetric instabilities on supermagnetosonic relativistic jets will grow is predicted theoretically and a Doppler boosted ``apparent'' emissivity is computed. Effects due to helically twisted filamentary structure produced by asymmetric modes of instability should be readily observable on relativistic jets.

We compare results from a relativistic and a nonrelativistic set of 2D axisymmetric jet simulations. For a set of five relativistic simulations that either increase the Lorentz factor or decrease the adiabatic index we compute nonrelativistic simulations with equal useful power or thrust. We examine these simulations for morphological and dynamical differences, focusing on the velocity field, the width of the cocoon, the age of the jets, and the internal structure of the jet itself. The primary result of these comparisons is that the velocity field of nonrelativistic jet simulations cannot be scaled up to give the spatial distribution of Lorentz factors seen in relativistic simulations. Since the local Lorentz factor plays a major role in determining the total intensity for parsec scale extragalactic jets, this suggests that a nonrelativistic simulation cannot yield the proper intensity distribution for a relativistic jet. Another general result is that each relativistic jet and its nonrelativistic equivalents have similar ages (in dynamical time units, = R/a_a, where R is the initial radius of a cylindrical jet and a_a is the sound speed in the ambient medium). In addition to these comparisons, we have completed four new relativistic simulations to investigate the effect of varying thermal pressure on relativistic jets. The simulations generally confirm that faster (larger Lorentz factor) and colder jets are more stable, with smaller amplitude and longer wavelength internal variations. The apparent stability of these jets does not follow from linear normal mode analysis, which suggests that there are available growing Kelvin-Helmholtz modes. (Abridged.)

A set of 3D MHD simulations of magnetized jets has been performed. The jets contain an equipartition primarily poloidal magnetic field and the effect of jet density on jet dynamics and structure is evaluated. The jet is precessed at the origin to excite Kelvin-Helmholtz unstable helical modes. We extensively compare the structure in these simulations with linear stability theory. The jet that is dense with respect to the external medium develops a high speed core surrounded by a less dense sheath consisting of slower moving jet fluid. These simulations suggest that extended extragalactic jets propagate to such large distances because they are surrounded by a lobe or cocoon whose density is less than the jet density. (Abridged abstract.)