We investigate the responsiveness of the 2798AA Mg II broad emission line in AGN on timescales of several years. Our study is based on a sample of extremely variable AGN as well as a broad population sample. The observed response of the line in previous studies has been mixed. By focussing on extreme variability ($|\Delta g|>$ 1) we find that Mg II clearly does respond to the continuum. However, the degree of responsiveness varies strikingly from one object to another; we see cases of Mg II changing by as much as the continuum, more than the continuum, or very little at all. In 74% of the highly variable sample the behaviour of Mg II corresponds with that of H$\beta$, with 30% of the objects showing large variations in both lines. We do not detect any change in the line width that would correspond to Broad Line Region `breathing', in accordance with results from literature. Some of the objects in our highly variable sample show a clear asymmetry in the Mg II profile. This skewness can be both to the blue and the red of the line centre. Results from our broad population sample show that highly variable quasars have lower Eddington ratios. This result holds for the variability of the continuum, but the correlation is significantly reduced for the variability of the Mg II line. For the first time, we present an overview of the value of the intrinsic Baldwin Effect for Mg II in a large sample.

We have carried out a Chandra X-ray and multi-frequency radio VLBA study of the AGN TXS 0128+554, which is associated with the Fermi gamma-ray source 4FGL J0131.2+5547. The AGN is unresolved in a target 19.3 ks Chandra image, and its spectrum is well fit by a simple absorbed power law model, with no distinguishable spectral features. Its relatively soft X-ray spectrum compared to other CSOs may be indicative of a thermal emission component, for which we were able to obtain an upper temperature limit of kT = 0.08 keV. The compact radio morphology and measured advance speed of 0.32c +- 0.07c indicate a kinematic age of only 82 y +- 17 y, placing TXS 0128+554 among the youngest members of the compact symmetric object (CSO) class. The lack of compact, inverted spectrum hotspots and an emission gap between the bright inner jet and outer radio lobe structure indicate that the jets have undergone episodic activity, and were re-launched a decade ago. The predicted gamma-ray emission from the lobes, based on an inverse Compton-emitting cocoon model, is three orders of magnitude below the observed Fermi LAT flux. A comparison to other Fermi-detected and non-Fermi detected CSOs with redshift z<0.1 indicates that the gamma-ray emission likely originates in the inner jet/core region, and that nearby, recently launched AGN jets are primary candidates for detection by the Fermi LAT instrument.

We present results from a parsec-scale jet kinematics study of 409 bright radio-loud AGNs based on 15 GHz VLBA data obtained between 1994 August 31 and 2016 December 26 as part of the 2cm VLBA survey and MOJAVE programs. We tracked 1744 individual bright features in 382 jets over at least five epochs. A majority (59%) of the best-sampled jet features showed evidence of accelerated motion at the >3sigma level. Although most features within a jet typically have speeds within ~40% of a characteristic median value, we identified 55 features in 42 jets that had unusually slow pattern speeds, nearly all of which lie within 4 pc (100 pc de-projected) of the core feature. Our results combined with other speeds from the literature indicate a strong correlation between apparent jet speed and synchrotron peak frequency, with the highest jet speeds being found only in low-peaked AGNs. Using Monte Carlo simulations, we find best fit parent population parameters for a complete sample of 174 quasars above 1.5 Jy at 15 GHz. Acceptable fits are found with a jet population that has a simple unbeamed power law luminosity function incorporating pure luminosity evolution, and a power law Lorentz factor distribution ranging from 1.25 to 50 with slope -1.4 +- 0.2. The parent jets of the brightest radio quasars have a space density of 261 +- 19 Gpc$^{-3}$ and unbeamed 15 GHz luminosities above ~$10^{24.5}$ W/Hz, consistent with FR II class radio galaxies.

Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full SDSS footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude (|\Delta g|>1 mag, |\Delta r|>0.5 mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted S/N threshold on variability of broad H\beta emission, we find 17 new CLQs, yielding a confirmation rate of >~ 20%. These candidates are at lower Eddington ratio relative to the overall quasar population which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 Angstroms increases from 1.5 to 6. We release a catalog of over 200 highly variable candidates to facilitate future CLQ searches.

We present 5321 milliarcsecond-resolution total intensity and linear polarization maps of 437 active galactic nuclei (AGNs) obtained with the VLBA at 15 GHz as part of the MOJAVE survey, and also from the NRAO data archive. The former is a long-term program to study the structure and evolution of powerful parsec-scale outflows associated with AGNs. The targeted AGNs are drawn from several flux-limited radio and gamma-ray samples, and all have correlated VLBA flux densities greater than about 50 mJy at 15 GHz. Approximately 80% of these AGNs are associated with gamma-ray sources detected by the Fermi LAT instrument. The vast majority were observed with the VLBA on 5 to 15 occasions between 1996 January 19 and 2016 December 26, at intervals ranging from a month to several years, with the most typical sampling interval being six months. A detailed analysis of the linear and circular polarization evolution of these AGN jets are presented in other papers in this series.

We aim to better understand the dynamics within relativistic magneto-hydrodynamical flows in the extreme environment and close vicinity of supermassive black holes. To do so, we analyze the peculiar radio galaxy 3C 111, for which long-term polarimetric observations are available. We make use of the high spatial resolution of the VLBA network and the MOJAVE monitoring program, which provides high data quality also for single sources and allows us to study jet dynamics on parsec scales in full polarization with an evenly sampled time-domain. We additionally consider data from the IRAM 30-m Telescope as well as the SMA. Jet properties such as the electric vectors, the (polarized) flux density, feature size, and brightness temperature, describe a complex evolution of the polarized jet. The electric vector position angles (EVPAs) of features traveling down the jet perform a large and smooth rotation of $\gtrsim 180^{\circ}$ across a distance of about 20 pc. In contrast, the EVPAs are strongly variable within the first parsecs of the jet. We find a tendency towards transverse EVPAs across the jet with a local anomaly of aligned vectors in between. The transverse extent of the flow decreases coincident with a jump in brightness temperature around where we observe the EVPAs to turn into alignment with the jet flow. Also the gradients of the feature size and particle density with distance steepen in that region. We interpret the propagating polarized features with shocks and the observed local anomalies with the interaction of these shocks with a recollimation shock of the underlying flow. Together with a sheared magnetic field, this shock-shock interaction can explain the large rotation of the EVPA. The superimposed variability of the EVPAs close to the core is likely related to a clumpy Faraday screen, which also contributes significantly to the observed EVPA rotation in that region.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present 1625 new 15 GHz (2 cm) VLBA images of 295 jets associated with active galactic nuclei (AGNs) from the MOJAVE and 2 cm VLBA surveys, spanning observations between 1994 Aug 31 and 2013 Aug 20. For 274 AGNs with at least 5 VLBA epochs, we have analyzed the kinematics of 961 individual bright features in their parsec-scale jets. A total of 122 of these jets have not been previously analyzed by the MOJAVE program. In the case of 451 jet features that had at least 10 epochs, we also examined their kinematics for possible accelerations. At least half of the well-sampled features have non-radial and/or accelerating trajectories, indicating that non-ballistic motion is common in AGN jets. Since it is impossible to extrapolate any accelerations that occurred before our monitoring period, we could only determine reliable ejection dates for about 24% of those features that had significant proper motions. The distribution of maximum apparent jet speeds in all 295 AGNs measured by our program to date is peaked below 5c, with very few jets with apparent speeds above 30c. The fastest speed in our survey is about 50c, measured in the jet of the quasar PKS 0805-07, and is indicative of a maximum jet Lorentz factor of about 50 in the parent population. The Fermi LAT-detected gamma-ray AGNs in our sample have, on average, higher jet speeds than non LAT-detected AGNs, indicating a strong correlation between pc-scale jet speed and gamma-ray Doppler boosting factor. We have identified 11 moderate-redshift (z<0.35) AGNs with fast apparent speeds (>10c) that are strong candidates for future TeV gamma-ray detection. Of the five gamma-ray loud narrow-lined Seyfert I AGNs in our sample, three show highly superluminal jet motions, while the others have sub-luminal speeds. (abridged)

We study the kinematics of ridge lines on the pc-scale jet of the active galactic nucleus BL Lac. We show that the ridge lines display transverse patterns that move superluminally downstream, and that the moving patterns are analogous to waves on a whip. Their apparent speeds $\beta_\mathrm{app}$ (units of $c$) range from 3.9 to 13.5, corresponding to $\beta_\mathrm{wave}^\mathrm{gal}= 0.981 - 0.998$ in the galaxy frame. We show that the magnetic field in the jet is well-ordered with a strong transverse component, and assume that it is helical and that the transverse patterns are Alfv\'en waves propagating downstream on the longitudinal component of the magnetic field. The wave-induced transverse speed of the jet is non-relativistic ($\beta_\mathrm{tr}^\mathrm{gal} \lesssim 0.09$). In 2010 the wave activity subsided and the jet then displayed a mild wiggle that had a complex oscillatory behaviour. The Alfv\'en waves appear to be excited by changes in the position angle of the recollimation shock, in analogy to exciting a wave on a whip by shaking the handle. A simple model of the system with plasma sound speed $\beta_\mathrm{s}=0.3$ and apparent speed of a slow MHD wave $\beta_\mathrm{app,S}=4$ yields Lorentz factor of the beam $\Gamma_\mathrm{beam} \sim 4.5$, pitch angle of the helix (in the beam frame) $\alpha\sim 67^\circ$, Alfv\'en speed $\beta_\mathrm{A}\sim 0.64$, and magnetosonic Mach number $M_\mathrm{ms}\sim 4.7$. This describes a plasma in which the magnetic field is dominant and in a rather tight helix, and Alfv\'en waves are responsible for the moving transverse patterns.

Several narrow-line Seyfert 1 galaxies (NLS1s) have now been detected in gamma rays, providing firm evidence that at least some of this class of active galactic nuclei (AGN) produce relativistic jets. The presence of jets in NLS1s is surprising, as these sources are typified by comparatively small black hole masses and near- or super-Eddington accretion rates. This challenges the current understanding of the conditions necessary for jet production. Comparing the properties of the jets in NLS1s with those in more familiar jetted systems is thus essential to improve jet production models. We present early results from our campaign to monitor the kinematics and polarization of the parsec-scale jets in a sample of 15 NLS1s through multifrequency observations with the Very Long Baseline Array. These observations are complemented by fast-cadence 15 GHz monitoring with the Owens Valley Radio Observatory 40m telescope and optical spectroscopic monitoring with with the 2m class telescope at the Guillermo Haro Astrophysics Observatory in Cananea, Mexico.

We report on the acceleration properties of 329 features in 95 blazar jets from the MOJAVE VLBA program. Nearly half the features and three-quarters of the jets show significant changes in speed and/or direction. In general, apparent speed changes are distinctly larger than changes in direction, indicating that changes in the Lorentz factors of jet features dominate the observed speed changes rather than bends along the line of sight. Observed accelerations tend to increase the speed of features near the jet base, $\lesssim 10-20$ parsecs projected, and decrease their speed at longer distances. The range of apparent speeds at fixed distance in an individual jet can span a factor of a few, indicating that shock properties and geometry may influence the apparent motions; however, we suggest that the broad trend of jet features increasing their speed near the origin is due to an overall acceleration of the jet flow out to de-projected distances of order $10^2$ parsecs, beyond which the flow begins to decelerate or remains nearly constant in speed. We estimate intrinsic rates of change of the Lorentz factors in the galaxy frame of order $\dot{\Gamma}/\Gamma \simeq 10^{-3}$ to $10^{-2}$ per year which can lead to total Lorentz factor changes of a factor of a few on the length scales observed here. Finally, we also find evidence for jet collimation at projected distances of $\lesssim 10$ parsecs in the form of the non-radial motion and bending accelerations that tend to better align features with the inner jet.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

Parsec-scale VLBA images of BL Lac at 15 GHz show that the jet contains a permanent quasi-stationary emission feature 0.26 mas (0.34 pc projected) from the core, along with numerous moving features. In projection, the tracks of the moving features cluster around an axis at position angle -166.6 deg that connects the core with the standing feature. The moving features appear to emanate from the standing feature in a manner strikingly similar to the results of numerical 2-D relativistic magneto-hydrodynamic (RMHD) simulations in which moving shocks are generated at a recollimation shock. Because of this, and the close analogy to the jet feature HST-1 in M87, we identify the standing feature in BL Lac as a recollimation shock. We assume that the magnetic field dominates the dynamics in the jet, and that the field is predominantly toroidal. From this we suggest that the moving features are compressions established by slow and fast mode magneto-acoustic MHD waves. We illustrate the situation with a simple model in which the slowest moving feature is a slow-mode wave, and the fastest feature is a fast-mode wave. In the model the beam has Lorentz factor about 3.5 in the frame of the host galaxy, and the fast mode wave has Lorentz factor about 1.6 in the frame of the beam. This gives a maximum apparent speed for the moving features 10c. In this model the Lorentz factor of the pattern in the galaxy frame is approximately 3 times larger than that of the beam itself.

We have obtained milliarcsecond-scale spectral index distributions for a sample of 190 extragalactic radio jets through the MOJAVE project. The sources were observed in 2006 at 8.1, 8.4, 12.1, and 15.4 GHz, and we study the four-frequency spectrum in individual jet features. We performed detailed simulations to study the effects of image alignment and (u,v)-plane coverage on the spectral index maps to verify our results. The core spectral indices are on average flat, while the jet spectrum is in general steep with a mean index of -1.04. A simple power-law fit is often inadequate for the core regions, as expected if the cores are partially self-absorbed. The overall jet spectrum steepens at a rate of about -0.001 to -0.004 per deprojected parsec when moving further out from the core with flat spectrum radio quasars having significantly steeper spectra than the BL Lac objects. The spectrum in both types of objects flattens at the locations of the jet components indicating particle acceleration or density enhancements along the jet. The mean spectral index at the component locations of -0.81 corresponds to a power-law index of ~2.6 for the electron energy distribution. We find a significant trend that jet components with linear polarization parallel to the jet (magnetic field perpendicular to the jet) have flatter spectra, as expected for transverse shocks. Compared to quasars, BL Lacs have more jet components with perpendicular magnetic field alignment, which may explain their generally flatter spectra. The overall steepening of the spectra with distance can be explained with radiative losses if the jets are collimating or with the evolution of the high-energy cutoff in the electron spectrum if the jets are conical. This interpretation is supported by a significant correlation with the age of the component and the spectral index, with older components having steeper spectra. (abridged)

We have developed a new Bayesian Markov Chain Monte Carlo method to deconvolve light curves of blazars into individual flares, including proper estimation of the fit errors. We use the method to fit 15\,GHz light curves obtained within the OVRO 40-m blazar monitoring program where a large number of AGN have been monitored since 2008 in support of the Fermi Gamma-Ray Space Telescope mission. The time scales obtained from the fitted models are used to calculate the variability brightness temperature of the sources. Additionally, we have calculated brightness temperatures of a sample of these objects using Very Long Baseline Array data from the MOJAVE survey. Combining these two data sets enables us to study the intrinsic brightness temperature distribution in these blazars at 15 GHz. Our preliminary results indicate that the mean intrinsic brightness temperature in a sample of 14 sources is near the equipartition brightness temperature of $\sim 10^{11}$K.

We describe the parsec-scale kinematics of 200 AGN jets based on 15 GHz VLBA data obtained between 1994 Aug 31 and 2011 May 1. We present new VLBA 15 GHz images of these and 59 additional AGN from the MOJAVE and 2 cm Survey programs. Nearly all of the 60 most heavily observed jets show significant changes in their innermost position angle over a 12 to 16 year interval, ranging from 10 deg to 150 deg on the sky, corresponding to intrinsic variations of ~0.5 deg to ~2 deg. The BL Lac jets show smaller variations than quasars. Roughly half of the heavily observed jets show systematic position angle trends with time, and 20 show indications of oscillatory behavior. The time spans of the data sets are too short compared to the fitted periods (5 to 12 y), however, to reliably establish periodicity. The rapid changes and large jumps in position angle seen in many cases suggest that the superluminal AGN jet features occupy only a portion of the entire jet cross section, and may be energized portions of thin instability structures within the jet. We have derived vector proper motions for 887 moving features in 200 jets having at least five VLBA epochs. For 557 well-sampled features, there are sufficient data to additionally study possible accelerations. We find that the moving features are generally non-ballistic, with 70% of the well-sampled features showing either significant accelerations or non-radial motions. Inward motions are rare (2% of all features), are slow (<0.1 mas per y), are more prevalent in BL Lac jets, and are typically found within 1 mas of the unresolved core feature. There is a general trend of increasing apparent speed with distance down the jet for both radio galaxies and BL Lac objects. In most jets, the speeds of the features cluster around a characteristic value, yet there is a considerable dispersion in the distribution. (abridged)

We report observations of Faraday rotation measures (RMs) for a sample of 191 extragalactic radio jets observed within the MOJAVE program. Multifrequency VLBA observations were carried out over twelve epochs in 2006 at four frequencies between 8 and 15 GHz. We detect parsec-scale Faraday RMs in 149 sources and find the quasars to have larger RMs on average than BL Lac objects. The median core RMs are significantly higher than in the jet components. This is especially true for quasars where we detect a significant negative correlation between the magnitude of the RM and the de-projected distance from the core. We perform detailed simulations of the observational errors of total intensity, polarization and Faraday rotation, and concentrate on the errors of transverse Faraday RM gradients in unresolved jets. Our simulations show that the finite image restoring beam size has a significant effect on the observed RM gradients, and spurious gradients can occur due to noise in the data if the jet is less than two beams wide in polarization. We detect significant transverse RM gradients in four sources (0923+392, 1226+023, 2230+114 and 2251+158). In 1226+023 the RM is for the first time seen to change sign from positive to negative over the transverse cuts, which supports the presence of a helical magnetic field in the jet. In this source we also detect variations in the jet RM over a time scale of three months, which are difficult to explain with external Faraday screens and suggest internal Faraday rotation. By comparing fractional polarization changes in jet components between the four frequency bands to depolarization models we find that an external purely random Faraday screen viewed through only a few lines of sight can explain most of our polarization observations but in some sources, such as 1226+023 and 2251+158, internal Faraday rotation is needed.

We have conducted a survey of Faraday rotation in a sample of 191 compact radio-loud AGNs as part of the MOJAVE (Monitoring of Jets in Active galactic nuclei with VLBA Experiments) project. The observations were carried out with the VLBA at 8.1, 8.4, 12.1 and 15.3 GHz over 12 epochs in 2006. We detect sufficiently strong linear polarization in 159 out of 211 observations to calculate the rotation measure values, resulting in a large enough sample for statistical analysis of the Faraday rotation in blazars. These Faraday rotation measures can be used to study the intrinsic magnetic field order and orientation in parsec-scale blazar jets. Our sample includes 119 sources listed in the 1FGL or 2FGL catalogs and we detect rotation measure values in 111 out of 131 maps. Of the 72 sources that are not in the gamma-ray catalogs we detect RM in 48 out of 80 maps. The median RM values of the LAT-detected sources do not differ significantly from the non-LAT-detected sources. Nine of the sources in our sample have resolved enough jets to study the transverse Faraday rotation structure, and we detect significant transverse rotation measure gradients in four sources. In two of these (3C~273 and 3C~454.3) there is additional evidence to support helical magnetic field in the parsec-scale jets. The two others (0923+392 and 2230+114) require further observations to identify the nature of the gradient. It is interesting that three of the four sources with significant rotation measure gradients are sources that have shown large gamma-ray flares.

Motivated by recent observations that show increasing fractional linear polarization with increasing wavelength in a small number of optically thin jet features, i.e. "inverse depolarization", we present a physical model that can explain this effect and may provide a new and complementary probe of the low energy particle population and possible helical magnetic fields in extragalactic radio jets. In our model, structural inhomogeneities in the jet magnetic field create cancellation of polarization along the line of sight. Internal Faraday rotation, which increases like wavelength squared, acts to align the polarization from the far and near sides of the jet, leading to increased polarization at longer wavelengths. Structural inhomogeneities of the right type are naturally produced in helical magnetic fields and will also appear in randomly tangled magnetic fields. We explore both alternatives and find that, for random fields, the length scale for tangling cannot be too small a fraction of the jet diameter and still be consistent with the relatively high levels of fractional polarization observed in these features. We also find that helical magnetic fields naturally produce transverse structure for inverse depolarization which may be observable even in partially resolved jets.

I review constraints on the physical properties of AGN jets revealed through Very Long Baseline Interferometry (VLBI) studies of the structure and time-evolution of parsec-scale jets, including recent results from the MOJAVE program. In particular I focus on constraints available from very long time baseline studies which probe a wide range of jet behavior over many outbursts. Kinematic studies of propagating jet features find an apparent speed distribution that peaks around 10c for blazars, with speeds up to 50c observed. These observed speeds require Lorentz factors at least as large, implying that parsec-scale Lorentz factors up to 10-20 are common for blazars with a tail up to ~ 50. Jet flows are still becoming organized on these scales as evidenced by the high incidence of non-radial motions and/or accelerations of jet features (including increases and decreases in apparent speed and direction). Changes in Lorentz factors of propagating jet features appear to play a significant role in the observed accelerations, and while the connection between acceleration of jet features and the underlying flow is not clear, the pattern of observed accelerations suggest the flow may increase in speed near the base of the jet and decrease further out. In some jets, ejections of new features span a range of ejection angles over many epochs, tracing out wider opening angles on parsec-scales than are apparent in single epoch observations.

Radio polarimetric observations of Active Galactic Nuclei can reveal the magnetic field structure in the parsec-scale jets of these sources. We have observed the gamma-ray blazar 3C 273 as part of our multi-frequency survey with the Very Long Baseline Array to study Faraday rotation in a large sample of jets. Our observations re-confirm the transverse rotation measure gradient in 3C 273. For the first time the gradient is seen to cross zero which is further indication for a helical magnetic field and spine-sheath structure in the jet. We believe the difference to previous epochs is due to a different part of the jet being illuminated in our observations.

We investigate the dependence of gamma-ray brightness of blazars on intrinsic properties of their parsec-scale radio jets and the implication for relativistic beaming. By combining apparent jet speeds derived from high-resolution VLBA images from the MOJAVE program with millimetre-wavelength flux density monitoring data from Metsahovi Radio Observatory, we estimate the jet Doppler factors, Lorentz factors, and viewing angles for a sample of 62 blazars. We study the trends in these quantities between the sources which were detected in gamma-rays by the Fermi Large Area Telescope (LAT) during its first three months of science operations and those which were not detected. The LAT-detected blazars have on average higher Doppler factors than non-LAT-detected blazars, as has been implied indirectly in several earlier studies. We find statistically significant differences in the viewing angle distributions between gamma-ray bright and weak sources. Most interestingly, gamma-ray bright blazars have a distribution of comoving frame viewing angles that is significantly narrower than that of gamma-ray weak blazars and centred roughly perpendicular to the jet axis. The lack of gamma-ray bright blazars at large comoving frame viewing angles can be explained by relativistic beaming of gamma-rays, while the apparent lack of gamma-ray bright blazars at small comoving frame viewing angles, if confirmed with larger samples, may suggest an intrinsic anisotropy or Lorentz factor dependence of the gamma-ray emission.

We discuss acceleration measurements for a large sample of extragalactic radio jets from the MOJAVE program which studies the parsec-scale jet structure and kinematics of a complete, flux-density-limited sample of Active Galactic Nuclei (AGN). Accelerations are measured from the apparent motion of individual jet features or "components" which may represent patterns in the jet flow. We find that significant accelerations are common both parallel and perpendicular to the observed component velocities. Parallel accelerations, representing changes in apparent speed, are generally larger than perpendicular acceleration that represent changes in apparent direction. The trend for larger parallel accelerations indicates that a significant fraction of these changes in apparent speed are due to changes in intrinsic speed of the component rather than changes in direction to the line of sight. We find an overall tendency for components with increasing apparent speed to be closer to the base of their jets than components with decreasing apparent speed. This suggests a link between the observed pattern motions and the underlying flow which, in some cases, may increase in speed close to the base and decrease in speed further out; however, common hydro-dynamical processes for propagating shocks may also play a role. About half of the components show "non-radial" motion, or a misalignment between the component's structural position angle and its velocity direction, and these misalignments generally better align the component motion with the downstream emission. Perpendicular accelerations are closely linked with non-radial motion. When observed together, perpendicular accelerations are usually in the correct direction to have caused the observed misalignment.

We discuss the jet kinematics of a complete flux-density-limited sample of 135 radio-loud active galactic nuclei (AGN) resulting from a 13 year program to investigate the structure and evolution of parsec-scale jet phenomena. Our analysis is based on new 2 cm Very Long Baseline Array (VLBA) images obtained between 2002 and 2007, but includes our previously published observations made at the same wavelength, and is supplemented by VLBA archive data. In all, we have used 2424 images spanning the years 1994-2007 to study and determine the motions of 526 separate jet features in 127 jets. The data quality and temporal coverage (a median of 15 epochs per source) of this complete AGN jet sample represents a significant advance over previous kinematics surveys. In all but five AGNs, the jets appear one-sided, most likely the result of differential Doppler boosting. In general the observed motions are directed along the jet ridge line, outward from the optically thick core feature. We directly observe changes in speed and/or direction in one third of the well-sampled jet components in our survey. While there is some spread in the apparent speeds of separate features within an individual jet, the dispersion is about three times smaller than the overall dispersion of speeds among all jets. This supports the idea that there is a characteristic flow that describes each jet, which we have characterized by the fastest observed component speed. The observed maximum speed distribution is peaked at ~10c, with a tail that extends out to ~50c. This requires a distribution of intrinsic Lorentz factors in the parent population that range up to ~50. We also note the presence of some rare low-pattern speeds or even stationary features in otherwise rapidly flowing jets... (abridged)

In its first three months of operations, the Fermi Gamma-Ray Observatory has detected approximately one quarter of the radio-flux-limited MOJAVE sample of bright flat-spectrum active galactic nuclei (AGNs) at energies above 100 MeV. We have investigated the apparent parsec-scale jet speeds of 26 MOJAVE AGNs measured by the Very Long Baseline Array (VLBA) that are in the LAT bright AGN sample (LBAS). We find that the gamma-ray bright quasars have faster jets on average than the non-LBAS quasars, with a median of 15 c, and values ranging up to 34 c. The LBAS AGNs in which the LAT has detected significant gamma-ray flux variability generally have faster jets than the nonvariable ones. These findings are in overall agreement with earlier results based on nonuniform EGRET data which suggested that gamma-ray bright AGNs have preferentially higher Doppler boosting factors than other blazar jets. However, the relatively low LAT detection rates for the full MOJAVE sample (24%) and previously known MOJAVE EGRET-detected blazars (43%) imply that Doppler boosting is not the sole factor that determines whether a particular AGN is bright at gamma-ray energies. The slower apparent jet speeds of LBAS BL Lac objects and their higher overall LAT detection rate as compared to quasars suggest that the former are being detected by Fermi because of their higher intrinsic (unbeamed) gamma-ray to radio luminosity ratios.

We have compared the radio emission from a sample of parsec-scale AGN jets as measured by the VLBA at 15 GHz, with their associated gamma-ray properties that are reported in the Fermi LAT 3-month bright source list. We find in our radio-selected sample that the gamma-ray photon flux correlates well with the quasi-simultaneously measured compact radio flux density. The LAT-detected jets in our radio-selected complete sample generally have higher compact radio flux densities, and their parsec-scale cores are brighter (i.e., have higher brightness temperature) than the jets in the LAT non-detected objects. This suggests that the jets of bright gamma-ray AGN have preferentially higher Doppler-boosting factors. In addition, AGN jets tend to be found in a more active radio state within several months from LAT-detection of their strong gamma-ray emission. This result becomes more pronounced for confirmed gamma-ray flaring sources. We identify the parsec-scale radio core as a likely location for both the gamma-ray and radio flares, which appear within typical timescales of up to a few months of each other.

We report the results of parsec-scale, multi-frequency VLBA observations of the core region of 3C 279 in Stokes I, linear polarization, and circular polarization. These full polarization spectra are modeled by radiative transfer simulations to constrain the magnetic field and particle properties of the parsec-scale jet in 3C 279. The polarization properties of the core region, including the amount of linear polarization, the amount and sign of Faraday rotation, and the amount and sign of circular polarization can be explained by a consistent physical picture. The base of the jet is modeled as an inhomogeneous Blandford-Konigl style conical jet dominated by a vector-ordered poloidal magnetic field along the jet axis, and we estimate its net magnetic flux. This poloidal field is responsible for the linear and circular polarization from this inhomogeneous component. Farther down the jet the magnetic field in two homogeneous features is dominated by local shocks and a smaller fraction of vector-ordered poloidal field remains along the jet axis. In this picture, we find the jet to be kinetically dominated by protons with the radiating particles being dominated by electrons at an approximate fraction of >~ 75%. Based on the amounts of Faraday conversion deduced for the homogeneous components, we find a plausible range for the lower cutoff in the relativistic particle energy spectrum to be 5 <~ gamma_l <~ 35. The physical picture described here is not unique if the observed Faraday rotation and depolarization occur in screens external to the jet; however, we find the joint explanation of linear and circular polarization observations from a single set of magnetic fields and particle properties internal to the jet to be compelling evidence for this picture. (Abridged)

We present images from a long term program (MOJAVE: Monitoring of Jets in AGN with VLBA Experiments) to survey the structure and evolution of parsec-scale jet phenomena associated with bright radio-loud active galaxies in the northern sky. The observations consist of 2424 15 GHz VLBA images of a complete flux-density limited sample of 135 AGN above declination -20 degrees, spanning the period 1994 August to 2007 September. These data were acquired as part of the MOJAVE and 2 cm Survey programs, and from the VLBA archive. The sample selection criteria are based on multi-epoch parsec-scale (VLBA) flux density, and heavily favor highly variable and compact blazars. The sample includes nearly all the most prominent blazars in the northern sky, and is well-suited for statistical analysis and comparison with studies at other wavelengths. Our multi-epoch and stacked-epoch images show 94% of the sample to have apparent one-sided jet morphologies, most likely due to the effects of relativistic beaming. Of the remaining sources, five have two-sided parsec-scale jets, and three are effectively unresolved by the VLBA at 15 GHz, with essentially all of the flux density contained within a few tenths of a milliarcsecond.

The major multi-epoch VLBA programs are described and discussed in terms of relativistic beaming models. Broadly speaking the observed kinematics are consistent with models having a parent population which is only mildly relativistic but with Lorentz factors extending up to about 30. While the collimation and acceleration appears to mainly occur close to the central engine, there is evidence of accelerations up to 1 kpc downstream. Generally the motion appears to be linear, but in some sources the motion follows a curved trajectory. In other sources, successive features appear to be ejected in different directions possibly the result of a precessing nozzle. The launch of GLAST in 2008 will offer new opportunities to study the relation between radio and gamma-ray activity, and possibly to locate the source of the gamma-ray emission. VSOP-2 will give enhanced resolution and will facilitate the study of the two-dimensional structure of relativistic jets, while RadioAstron will provide unprecedented resolution to study the fine scale structure of the jet base.

In 1996, a major radio flux-density outburst occured in the broad-line radio galaxy 3C111. It was followed by a particularly bright plasma ejection associated with a superluminal jet component, which has shaped the parsec-scale structure of 3C111 for almost a decade. Here, we present results from 18 epochs of Very Long Baseline Array (VLBA) observations conducted since 1995 as part of the VLBA 2 cm Survey and MOJAVE monitoring programs. This major event allows us to study a variety of processes associated with outbursts of radio-loud AGN in much greater detail than has been possible in other cases: the primary perturbation gives rise to the formation of a leading and a following component, which are interpreted as a forward and a backward-shock. Both components evolve in characteristically different ways and allow us to draw conclusions about the work flow of jet-production events; the expansion, acceleration and recollimation of the ejected jet plasma in an environment with steep pressure and density gradients are revealed; trailing components are formed in the wake of the primary perturbation possibly as a result of coupling to Kelvin-Helmholtz instability pinching modes from the interaction of the jet with the external medium. The interaction of the jet with its ambient medium is further described by the linear-polarization signature of jet components traveling along the jet and passing a region of steep pressure/density gradients.

We discuss results from a decade long program to study the fine-scale structure and the kinematics of relativistic AGN jets with the aim of better understanding the acceleration and collimation of the relativistic plasma forming AGN jets. From the observed distribution of brightness temperature, apparent velocity, flux density, time variability, and apparent luminosity, the intrinsic properties of the jets including Lorentz factor, luminosity, orientation, and brightness temperature are discussed. Special attention is given to the jet in M87, which has been studied over a wide range of wavelengths and which, due to its proximity, is observed with excellent spatial resolution. Most radio jets appear quite linear, but we also observe curved non-linear jets and non-radial motions. Sometimes, different features in a given jet appear to follow the same curved path but there is evidence for ballistic trajectories as well. The data are best fit with a distribution of Lorentz factors extending up to gamma ~30 and intrinsic luminosity up to ~10^26 W/Hz. In general, gamma-ray quasars may have somewhat larger Lorentz factors than non gamma-ray quasars. Initially the observed brightness temperature near the base of the jet extend up to ~5x10^13 K which is well in excess of the inverse Compton limit and corresponds to a large excess of particle energy over magnetic energy. However, more typically, the observed brightness temperatures are ~2x10^11 K, i.e., closer to equipartition.

We report new 2 cm VLBA images of the inner radio jet of M87 showing a limb brightened structure and unambiguous evidence for a faint 3 mas long counter-feature which also appears limb brightened. Multi-epoch observations of seven separate jet features show typical speeds of less than a few percent of the speed of light, despite the highly asymmetric jet structure and the implications of the canonical relativistic beaming scenario. The observed morphology is consistent with a two stream spine-sheath velocity gradient across the jet, as might be expected from the recently discovered strong and variable TeV emission as well as from numerical modeling of relativistic jets. Considering the large jet to counter-jet flux density ratio and lack of observed fast motion in the jet, we conclude that either the inner part of the M87 jet is intrinsically asymmetric or that the bulk plasma flow speed is much greater than any propagation of shocks or other pattern motions.

Relations between the observed quantities for a beamed radio jet, apparent transverse speed and apparent luminosity (beta_app,L), and the intrinsic quantities, Lorentz factor and intrinsic luminosity (gamma,L_o), are investigated. The inversion from measured to intrinsic values is not unique, but approximate limits to gamma and L_o can be found using probability arguments. Roughly half the sources in a flux density--limited, beamed sample have a value of gamma close to the measured beta_app. The methods are applied to observations of 119 AGN jets made with the VLBA at 15 GHz during 1994-2002. The results strongly support the common relativistic beam model for an extragalactic radio jet. The (beta_app,L) data are closely bounded by a theoretical envelope, an aspect curve for gamma=32, L_o= 10^25 W/Hz. This gives limits to the maximum values of gamma and L_o in the sample: gamma_max about 32, and L_o,max ~ 10^26 W/Hz. No sources with both high beta_app and low L are observed. This is not the result of selection effects due to the observing limits, which are flux density S>0.5 Jy, and angular velocity mu<4 mas/yr. Many of the fastest quasars have a pattern Lorentz factor gamma_p close to that of the beam, gamma_b, but some of the slow quasars must have gamma_p<<gamma_b. Three of the 10 galaxies in the sample have a superluminal feature, with speeds up to beta_app about 6. The others are at most mildly relativistic. The galaxies are not off-axis versions of the powerful quasars, but Cygnus A might be an exception.

We present a new method for studying the intrinsic brightness temperatures of the parsec-scale jet cores of Active Galactic Nuclei (AGN). Our method uses observed superluminal motions and observed brightness temperatures for a large sample of AGN to constrain the characteristic intrinsic brightness temperature of the sample as a whole. To study changes in intrinsic brightness temperature, we assume that the Doppler factors of individual jets are constant in time as justified by their relatively small changes in observed flux density. We find that in their median-low brightness temperature state, the sources in our sample have a narrow range of intrinsic brightness temperatures centered on a characteristic temperature, T_int = 3 x 10^10 K, which is close to the value expected for equipartition, when the energy in the radiating particles equals the energy stored in the magnetic fields. However, in their maximum brightness state, we find that sources in our sample have a characteristic intrinsic brightness temperature greater than 2 x 10^11 K, which is well in excess of the equipartition temperature. In this state, we estimate the energy in radiating particles exceeds the energy in the magnetic field by a factor of ~ 10^5. We suggest that the excess of particle energy when sources are in their maximum brightness state is due to injection or acceleration of particles at the base of the jet. Our results suggest that the common method of estimating jet Doppler factors by using a single measurement of observed brightness temperature and/or the assumption of equipartition may lead to large scatter or systematic errors in the derived values.

We report first-epoch circular polarization results for 133 active galactic nuclei in the MOJAVE program to monitor the structure and polarization of a flux limited sample of extra-galactic radio jets with the VLBA at 15 GHz. We found strong circular polarization ($\geq 0.3$%) in approximately 15% of our sample. The circular polarization was usually associated with jet cores; however, we did find a few strong jet components to be circularly polarized. The levels of circular polarization were typically in the range of $0.3-0.5$% of the local Stokes-$I$. We found no strong correlations between fractional circular polarization of jet cores and source type, redshift, EGRET detections, linear polarization, or other observed parsec-scale jet properties. There were differences between the circular-to-linear polarization ratios of two nearby galaxies versus more distant quasars and BL Lac objects. We suggest this is because the more distant sources either have (1) less depolarization of their linear polarization, and/or (2) poorer effective linear resolution and therefore their VLBA cores apparently contain a larger amount of linearly polarized jet emission. The jet of 3C 84 shows a complex circular polarization structure, similar to observations by Homan & Wardle five years earlier; however, much of the circular polarization seems to have moved, consistent with a proper motion of 0.06$c$. The jet of 3C 273 also has several circularly polarized components, and we find their fractional circular polarization decreases with distance from the core.

We have used VLBA fringe visibility data obtained at 15 GHz to examine the compact structure in 250 extragalactic radio sources. For 171 sources in our sample, more than half of the total flux density seen by the VLBA remains unresolved on the longest baselines. There are 163 sources in our list with a median correlated flux density at 15 GHz in excess of 0.5 Jy on the longest baselines. For about 60% of the sources, we have at least one observation in which the core component appears unresolved (generally smaller than 0.05 mas) in one direction, usually transverse to the direction into which the jet extends. BL Lacs are on average more compact than quasars, while active galaxies are on average less compact. Also, in an active galaxy the sub-milliarcsecond core component tends to be less dominant. IDV sources typically have a more compact, more core-dominated structure on sub-milliarcsecond scales than non-IDV sources, and sources with a greater amplitude of intra-day variations tend to have a greater unresolved VLBA flux density. The objects known to be GeV gamma-ray loud appear to have a more compact VLBA structure than the other sources in our sample. This suggests that the mechanisms for the production of gamma-ray emission and for the generation of compact radio synchrotron emitting features are related. The brightness temperature estimates and lower limits for the cores in our sample typically range between 10^11 and 10^13 K, but they extend up to 5x10^13 K, apparently in excess of the equipartition brightness temperature, or the inverse Compton limit for stationary synchrotron sources. The largest component speeds are observed in radio sources with high observed brightness temperatures, as would be expected from relativistic beaming (abridged).

We present sub-milliarcsecond resolution total intensity and linear polarization VLBI images of 3C273, using concurrent 43 and 86 GHz data taken with the Very Long Baseline Array in May 2002. The structure seen in the innermost jet suggest that we have fortuitously caught the jet in the act of changing direction. The polarization images confirm that the core is unpolarized (fractional polarization m < 1 %) at 86 GHz, but also show well ordered magnetic fields (m ~ 15 %) in the inner jet, at a projected distance of 2.3 pc from the core. In this strongly polarized region, the rotation measure changes across the jet by 4.2 x 10^{4} rad m^{-2} over an angular width of about 0.3 milliarcseconds. If the lack of polarization in the core is also attributed to a Faraday screen, then a rotation measure dispersion > 5.2 x 10^{4} rad m^{-2} must be present in or in front of that region. These are among the highest rotation measures reported so far in the nucleus of any active galaxy or quasar, and must occur outside (but probably close to) the radio emitting region. The transverse rotation measure gradient is in the same sense as that observed by Asada et al and by Zavala and Taylor at greater core distances. The magnitude of the transverse gradient decreases rapidly with distance down the jet, and appears to be variable.

We present first-epoch, milliarcsecond-scale linear polarization images at 15 GHz of 133 jets associated with AGN in the MOJAVE survey (Monitoring of Jets in Active Galactic Nuclei with VLBA Experiments). The sample consists of all known AGN with galactic latitude |b| > 2.5 deg., J2000 declination > -20 deg., and 15 GHz VLBA flux density exceeding 1.5 Jy (2 Jy for sources with declination < 0) at any epoch during the period 1994-2003. Because of strong selection effects, the sample primarily consists of blazars with parsec-scale morphologies consisting of a bright core component at the extreme end of a one-sided jet. At least one third of the compact cores are completely unresolved on the longest VLBA baselines, indicating brightness temperatures above 10^{11} K. The unresolved cores tend to have electric vectors that are aligned with the inner jet direction, suggesting the presence of a stationary transverse shock near the base of the jet. Many of the extended jet regions display exceedingly high fractional polarizations (> 50%) and electric vectors aligned with the jet ridge line, consistent with optically thin emission from transverse shocks. Both weak- and strong-lined blazars show a general increase in fractional polarization with distance down the jet, but BL Lac jets are generally more polarized and have electric vectors preferentially aligned with the local jet direction. We show that these differences are intrinsic to the jets, and not due to sample bias. Distinct features in the jets of gamma-ray loud (EGRET) blazars are typically twice as luminous as those in non-EGRET blazars, and are more polarized. These differences can be adequately explained if EGRET blazars have slightly higher Doppler boosting factors than those not yet detected in gamma-rays.

We present a status report of dedicated circular polarization (CP) global VLBI observations of SS433. The total intensity image of the source is presented. We do not detect linear polarization in the core region of the source. We estimate the errors in D-terms for the EVN, and its effect on the final right-left circular polarization gain calibration. The data from this and earlier observations are analyzed with the "zero-V" self-calibration technique, in order to obtain indications for detectable CP in different regions of the source. Our data show no evidence for strong gyro-synchrotron (GS) emission on the size scale probed by our beam.

We report the results of a 15 GHz (2 cm) multi-epoch VLBA program, begun in 1994 to study the outflow in radio jets ejected from quasars and active galaxies. The observed flow of 208 distinct features measured in 110 quasars, active galaxies, and BL Lac objects shows highly collimated relativistic motion with apparent transverse velocities typically between zero and about 15c, with a tail extending up to about 34c. Within individual jets, different features appear to move with a similar characteristic velocity which may represent an underlying continuous jet flow, but we also see some stationary and even apparently inward moving features which co-exist with the main features. Comparison of our observations with published data at other wavelengths suggests that there is a systematic decrease in apparent velocity with increasing wavelength, probably because the observations at different wavelengths sample different parts of the jet structure. The observed distribution of linear velocities is not consistent with any simple ballistic model. Either there is a rather broad range of Lorentz factors, a significant difference between the velocity of the bulk relativistic flow and the pattern speed of underlying shocks, or a combination of these options. Assuming a ballistic flow, comparison of observed apparent velocities and Doppler factors computed from the time scale of flux density variations is consistent with a steep power law distribution of intrinsic Lorentz factors, an isotropic distribution of orientations of the parent population, and intrinsic brightness temperatures about an order of magnitude below the canonical inverse Compton limit. It appears that the parent population of radio jets is not dominated by highly relativistic flows (abridged).

The sensitivity, stability, and uniformity of calibration of the VLBA has revolutionized parsec-scale polarization studies of AGN jets. Not only does polarization probe the magnetic field structures of jets, serving as a hydrodynamic tracer of shocks, bends, and shear, but polarization also probes the medium through which it propagates by encoding the signature of Faraday effects along the line of sight. I review advances made by the VLBA in studying the polarization of jets to probe their magnetic field structures, properties of the jet plasma, and properties of the external environment. This review covers both linear and circular polarization and is set in the context of outstanding questions in the field.

We present multi-frequency, high resolution VLBA circular polarization images of the radio source 3C 84 in the center of NGC 1275. Our images reveal a complex distribution of circular polarization in the inner parsec of the radio jet, with local levels exceeding 3% polarization, the highest yet detected with VLBI techniques. The circular polarization changes sign along the jet, making 3C 84 also the first radio jet to show both signs of circular polarization simultaneously. The spectrum and changing sign of the circular polarization indicate that it is unlikely to be purely intrinsic to the emitted synchrotron radiation. The Faraday conversion process makes a significant and perhaps dominant contribution to the circular polarization, and the observed spectrum suggests the conversion process is near saturation. The sign change in the circular polarization along the jet may result from this saturation or may be due to a change in magnetic field order after an apparent bend in the jet. From the small spatial scales probed here, ~ 0.15 pc, and the comparably high levels of circular polarization inferred for the intra-day variable source PKS 1519-273, we suggest a connection between small spatial scales and efficient production of circular polarization.

This contribution presents preliminary results from coordinated polarization sensitive VSOP and VLBA imaging of the blazars 3C~279 and 3C~345 at multiple frequencies.

We discuss the production of circular polarization in compact radio sources both by the intrinsic mechanism and by Faraday conversion. We pay particular attention to the magnetic field structure, considering partially ordered fiel ds and Laing sheets, and distinguishing between uniform and unidirectional fields. (The latter can be constrained b y flux conservation arguments.) In most cases, Faraday conversion is the more important mechanism. Conversion opera tes on Stokes U, which can be generated by internal Faraday rotation, or by magnetic field fluctuations, which can therefore produce circular polarization even in a pure pair plasma. We also show that the spectrum of circular pola rization in an inhomogeneous jet can be quite different from that in a uniform source, being flat or even inverted.

We investigate the relation between the Doppler factor determined from variations in total flux at 22 and 37 GHz, and the apparent transverse velocity determined from VLBA observations at 2 cm. The data are consistent with the relativistic beaming theory for compact radio sources, in that the distribution of beta_{app}/delta_{var}, for 30 quasars, is roughly consistent with a Monte Carlo simulation. The intrinsic temperature appears to be ~2x10^{10} K, close to the "equipartition value" calculated by Readhead (1994). We deduce the distribution of Lorentz factors for a group of 48 sources; the values range up to about gamma=40.

High resolution studies of circular polarization allow us see where it arises in a jet, study its local fractional level and spectrum, and compare these results to local measures of linear polarization and Faraday rotation. Here we not only review past results from Very Long Baseline Array (VLBA) circular polarization studies, but we also present preliminary new results on two quasars. In the core of PKS 0607-157, we find strong circular polarization at 8 GHz and much weaker levels at 15 GHz. Combined with the linear polarization data, we favor a simple model where the circular is produced by Faraday conversion driven by a small amount of Faraday rotation. In the core of 3C345, we find strong circular polarization at 15 GHz in a component with distinct linear polarization. This core component is optically thick at 8 GHz, where we detect no circular polarization. With opposite trends in frequency for PKS 0607-157 and 3C345, it seems clear that local conditions in a jet can have a strong effect on circular polarization and need to be taken into account when studying inhomogeneous objects with multi-frequency observations.

We report a change in the trajectory of a well studied jet component of the quasar 3C279. The component changes in apparent projected speed and direction, and we find it to be moving with a Lorentz factor gamma >~ 15 at an initial angle of <~ 1 degree to the line of sight. The new trajectory of the component has approximately the same speed and direction as an earlier superluminal feature, originally seen in the early 1970s. The new direction for the component is also much better aligned with larger scale VLBA and VLA structure out to 0.1 arcseconds. We suggest that the trajectory change is a collimation event occurring at >~ 1 kiloparsec (deprojected) along the jet. While the change in trajectory on the sky appears to be 26 degrees the intrinsic change is <~ 1 degree. We estimate the Doppler factor prior to the change in direction to be D >~ 28 and after the change to be D >~ 23. Comparison to independent constraints on the Doppler factor suggest that the energy in the radiating particles cannot greatly exceed the energy in the magnetic field unless the volume filling factor is very much less than one.

High resolution radio observations remain the most direct way to study the formation and evolution of radio jets associated with the accretion onto massive black holes. We report preliminary results of our seven year VLBA observational program to understand the nature of relativistic beaming in blazars and the surrounding environment of massive black holes. Most blazars show an apparent outward flow away from an active core. However, in a few sources the motion appears inward, most likely the result of projection of a curved trajectory which bends back toward along the line of sight. The apparent motion of jet features is not always oriented along the direction separating the feature from the core, and in a few cases we have observed a clear change in the direction and velocity of a feature as it flows along the jet. In other sources, the motion appears to follow a simple ballistic trajectory. We find no simple relation between the time scales of flux density changes and apparent component velocities.

The gamma-ray blazar PKS 1510-089 has a highly superluminal milli-arcsecond jet at a position angle (PA) of -28 degrees and an arcsecond jet with an initial PA of 155 degrees. With a PA difference of 177 degrees between the arcsecond and milli-arcsecond jets, PKS 1510-089 is perhaps the most highly misaligned radio jet ever observed and serves as a graphic example of projection effects in a highly beamed relativistic jet. Here we present the results of observations designed to bridge the gap between the milli-arcsecond and arcsecond scales. We find that a previously detected ``counter-feature'' to the arcsecond jet is directly fed by the milli-arcsecond jet. This feature is located 0.3" from the core, corresponding to a de-projected distance of 30 kiloparsecs. The feature appears to be dominated by shocked emission and has an almost perfectly ordered magnetic field along its outside edge. We conclude that it is most likely a shocked bend, viewed end-on, where the jet crosses our line of sight to form the southern arcsecond jet. While the bend appears to be nearly 180 degrees when viewed in projection, we estimate the intrinsic bending angle to be between 12 and 24 degrees. The cause of the bend is uncertain; however, we favor a scenario where the jet is bent after it departs the galaxy, either by ram pressure due to winds in the intracluster medium or simply by the density gradient in the transition to the intergalactic medium.

We are investigating the kinematics of jets in active galactic nuclei on parsec scales by studying a representative population of sources. This study is being carried out using the Very Long Baseline Array at 15 GHz, with more than 800 images taken since 1994. In this contribution we present an overview of the diversity of kinematics for a complete sample of sources.

We present analysis of the flux and polarization variability of parsec scale radio jets from a dual-frequency, six-epoch, VLBA polarization experiment monitoring 12 blazars. The observations were made at 15 and 22 GHz at bimonthly intervals over 1996. Here we analyze the flux, fractional polarization, and polarization position angle behavior of core regions and jet features, considering both the linear trends of these quantities with time and more rapid fluctuations about the linear trends. The dual frequency nature of the observations allows us to examine spectral evolution, to separate Faraday effects from changes in magnetic field order, and also to deduce empirical estimates for the uncertainties in measuring properties of VLBI jet features (abridged).

We present a detailed recipe for making movies from multi-epoch radio observations of astronomical sources. Images are interpolated linearly in time to create a smooth succession of frames so that a continuous movie can be compiled. Here, we outline the procedure, and draw attention to specific details necessary for making a successful movie. In particular, we discuss the issues pertaining specifically to making polarization movies. The procedure described here has been implemented into scripts in NRAO's AIPS package (Brandeis AIPS Movie Maker -- BAMM) that are available for public use (http://www.astro.brandeis.edu).

We present circular polarization results from a 5 GHz survey of the parsec-scale polarization properties of 40 AGN made with the VLBA. We find 11 circular polarization detections at the 3 sigma level or higher. This nearly quadruples the number of sources detected in circular polarization at VLBI resolution. We find no correlation between fractional linear and circular polarization across our sample. A likely explanation is external Faraday depolarization in the cores of AGN which reduces linear polarization but leaves circular polarization unchanged. In comparing ours and other recent results to observations made 20 years ago, we find that, in five of six cases, sources have the same sign of circular polarization today as they did 20 years ago. This suggests the presence of a long term property of the jets, such as the polarity of a net magnetic flux, which is stable on time-scales much longer than those of individual outbursts.

We review recent observations of circularly polarized radiation from AGN made with the VLBA and with the ATCA. We also discuss briefly the detections of the Galactic sources Sag A* and SS433. The origin of the circular polarization is still an open question in most cases, and we discuss four possible mechanisms. Detectable circular polarization is a common property of quasars, but not of radio galaxies, and is always associated with the compact core. There is growing evidence that the sign of the circular polarization stays the same over at least two or three decades in time, suggesting it is a fundamental property of the jet.

We present proper motions obtained from a dual frequency, six-epoch, VLBA polarization experiment monitoring a sample of 12 blazars. The observations were made at 15 GHz and 22 GHz at bi-monthly intervals over 1996. Ten of the eleven sources for which proper motion could be reliably determined are superluminal. Only J2005+77 has no superluminal components. Three sources (OJ287, J1224+21, and J1512-09) show motion faster than 10h^{-1}c, requiring $\gamma_{pattern}$ of at least 10h^{-1}c (H_0 = 100h km/s/Mpc). We compare our results to those in the literature and find motions outside the previously observed range for four sources. While some jet components exhibit significant non-radial motion, most motion is radial. In at least two sources there are components moving radially at significantly different structural position angles. In five of six sources (3C120, J1224+21, 3C273, 3C279, J1512-09, and J1927+73) that have multiple components with measurable proper motion, the innermost component is significantly slower than the others, suggesting that acceleration occurs in the jet. In the motions of individual components we observe at least one decelerating motion and two ``bending'' accelerations which tend to align their motions with larger scale structure. We also discuss in detail our procedures for obtaining robust kinematical results from multi-frequency VLBI data spanning several epochs.

We present five epochs of 15 GHz VLBA observations of 13 AGN. These observations were specially calibrated to detect parsec scale circular polarization and our calibration techniques are discussed and analyzed in detail. We obtained reliable detections of parsec scale circular polarization in the radio jets of 4 AGN: 3C84, PKS 0528+134, 3C273, and 3C279. For each of these objects our detections are at the level of ~ 0.3-1% local fractional circular polarization. Our observations are consistent across multiple epochs (and different calibration techniques) in the sign and magnitude of the circular polarization observed. 3C273 and 3C279 both undergo core outbursts during our observations and changes in the circular polarization of both sources are correlated with these outbursts. In general, we observe the circular polarization to be nearly coincident with the strong VLBI cores of these objects; however, in 3C84 the circular polarization is located a full milli-arcsecond south of the source peak, and in the 1996.73 epoch of 3C273 the circular polarization is predominately associated with the newly emerging jet component. Circular polarization may be produced as an intrinsic component of synchrotron radiation or by the Faraday conversion of linear to circular polarization. Our single frequency observations do not easily distinguish between these possible mechanisms, but independent of mechanism, the remarkable consistency across epoch of the sign of the observed circular polarization suggests the existence of a long term, stable, uni-directional magnetic field. Single dish observations of 3C273 and 3C279 at 8 GHz by Hodge and Aller suggest that this stability may persist for decades in our frame of observation.

We present a new technique for directly measuring the distances to superluminal radio sources. By comparing the observed proper motions of components in a parsec scale radio jet to their measured Doppler factors, we can deduce the distance to the radio source independent of the standard rungs in the cosmological distance ladder. This technique requires that the jet angle to the line of sight and the ratio of pattern to flow velocities are sufficiently constrained. We evaluate a number of possibilities for constraining these parameters and demonstrate the technique on a well defined component in the parsec scale jet of the quasar 3C279 (z = 0.536). We find an angular size distance to 3C279 of greater than 1.8 (+0.5,-0.3) n^{1/8} Gpc, where n is the ratio of the energy density in the magnetic field to the energy density in the radiating particles in that jet component. For an Einstein-de Sitter Universe, this measurement would constrain the Hubble constant to be H < 65 n^{-1/8} km/s/Mpc at the two sigma level. Similar measurements on higher redshift sources may help discriminate between cosmological models.