This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

The currents carried by the jets of active galactic nuclei (AGNs) can be probed using maps of the Faraday rotation measure (RM), since a jet current will be accompanied by a toroidal magnetic (B) field, which will give rise to a systematic change in the RM across the jet. The aim of this study is to identify new AGNs displaying statistically significant transverse RM gradients across their parsec-scale jets, and to look for overall patterns in the implied directions for the toroidal B-field components and jet currents. We have carried out new analyses of Faraday RM maps derived from previously published 8.1, 8.4, 12.1 and 15.3 GHz data obtained in 2006 on the NRAO VLBA. In a number of important ways, our procedures were identical to those of the original authors, but several other key aspects of the new imaging and analysis differ from the original methods. Our new analysis has substantially increased the number of AGNs known to display transverse RM gradients. The collected data on parsec and kiloparsec scales indicate that the current typically flows inward along the jet axis and outward in a more extended region surrounding the jet, typical to the current structure of a coaxial cable, accompanied by a self-consistent system of nested helical B fields, whose toroidal components give rise to the observed transverse RM gradients. These new results make it possible for the first time to conclusively demonstrate the existence of a preferred direction for the toroidal B-field components - and therefore of the currents - of AGN jets. Discerning the origin of this current-field system is of cardinal importance for understanding the physical mechanisms leading to the formation of the intrinsic jet B field, which likely plays an important role in the propagation and collimation of the jets; one possibility is the action of a "cosmic battery".

A number of groups have recently been active in searching for gradients in the ob- served Faraday rotation measure (RM) across jets of Active Galactic Nuclei (AGNs) on various scales and estimating their reliability. Such RM structures provide direct evidence for the presence of an azimuthal magnetic field component, which may be associated with a helical jet magnetic field, as is expected based on the results of many theoretical studies. We present new parsec-scale RM maps of 4 AGNs here, and analyze their transverse RM structures together with those for 5 previously published RM maps. All these jets display transverse RM gradients with significances of at least 3{\sigma}. This is part of an ongoing effort to establish how common transverse RM gradients that may be associated with helical or toroidal magnetic fields are in AGNs on parsec scales.

The formation of relativistic jets in active galactic nuclei (AGN) is related to accretion on to their central supermassive black holes, and magnetic fields are believed to play a central role in launching, collimating and accelerating the jet streams from very compact regions out to kiloparsec or megaparsec scales. In the presence of helical or toroidal magnetic fields threading the AGN jets and their immediate vicinity, gradients in the observed Faraday rotation measures are expected due to the systematic change in the line-of-sight component of the magnetic field across the jet. We have analysed total intensity, linear polarization, fractional polarization and Faraday rotation maps based on Very Long Baseline Array data obtained at four wavelengths in the 18-22 cm range for six AGN (OJ 287, 3C 279, PKS 1510-089, 3C 345, BL Lac and 3C 454.3). These observations typically probe projected distances out to tens of parsecs from the observed core, and are well suited for Faraday rotation studies due to the relatively long wavelengths used and the similarity of the structures measured at the different wavelengths. We have identified statistically significant, monotonic, transverse Faraday rotation gradients across the jets of four of these six sources, as well as a tentative transverse Faraday rotation gradient across the jet of OJ 287, providing evidence for the presence of toroidal magnetic fields, which may be one component of helical magnetic fields associated with these AGN jets.

Radio polarisation images of the jets of Active Galactic Nuclei (AGN) can provide a deep insight into the launching and collimation mechanisms of relativistic jets. However, even at VLBI scales, resolution is often a limiting factor in the conclusions that can be drawn from observations. The Maximum Entropy Method (MEM) is a deconvolution algorithm that can outperform the more common CLEAN algorithm in many cases, particularly when investigating structures present on scales comparable to or smaller than the nominal beam size with "super-resolution". A new implementation of the MEM suitable for single- or multiple-wavelength VLBI polarisation observations has been developed and is described here. Monte Carlo simulations comparing the performances of CLEAN and MEM at reconstructing the properties of model images are presented; these demonstrate the enhanced reliability of MEM over CLEAN when images of the fractional polarisation and polarisation angle are constructed using convolving beams that are appreciably smaller than the full CLEAN beam. The results of using this new MEM software to image VLBA observations of the AGN 0716+714 at six different wavelengths are presented, and compared to corresponding maps obtained with CLEAN. MEM and CLEAN maps of Stokes $I$, the polarised flux, the fractional polarisation and the polarisation angle are compared for convolving beams ranging from the full CLEAN beam down to a beam one-third of this size. MEM's ability to provide more trustworthy polarisation imaging than a standard CLEAN-based deconvolution when convolving beams appreciably smaller than the full CLEAN beam are used is discussed.

Helical magnetic fields embedded in the jets of active galactic nuclei (AGNs) are required by the broad range of theoretical models that advocate for electromagnetic launching of the jets. In most models, the direction of the magnetic field is random, but if the axial field is generated by a Cosmic Battery generated by current in the direction of rotation in the accretion disk, there is a correlation between the directions of the spin of the AGN accretion disk and of the axial field, which leads to a specific direction for the axial electric current, azimuthal magnetic field, and the resulting observed transverse Faraday-rotation (FR) gradient across the jet, due to the systematic change in the line-of-sight magnetic field. We consider new observational evidence for the presence of a nested helical magnetic-field structure such as would be brought about by the operation of the Cosmic Battery, and make predictions about the expected behavior of transverse FR gradients observed on decaparsec and kiloparsec scales. The collected results can be understood if the dominant azimuthal field on parsec scales corresponds to an axial electric current flowing inward along the jet, whereas the (weaker) dominant azimuthal field on kiloparsec scales corresponds to a outward-flowing current in the outer sheath of the jet and/or an extended disk wind.

A monotonic, statistically significant gradient in the observed Faraday Rotation Measure (RM) across the jet of an Active Galactic Nucleus (AGN) reflects a corresponding gradient in the electron density and/or line-of-sight magnetic (B) field. Such gradients may indicate the presence of a toroidal B field component, possibly associated with a helical jet B field. Although transverse RM gradients have been reported across a number of parsec-scale AGN jets, the same is not true on kiloparsec scales, suggesting that other (e.g. random) B-field components usually dominate on these larger scales. We have identified an extended, monotonic transverse RM gradient across the Northern lobe of a previously published Very Large Array (kiloparsec-scale) RM image of 5C4.114. We reanalyzed these VLA data in order to determine the significance of this RM gradient. The RM gradient across the Northern kiloparsec-scale lobe structure of 5C4.114 has a statistical significance of about 4sigma. There is also a somewhat less prominent monotonic transverse RM gradient across the Southern jet/lobe (significance ~ 3sigma). Other parts of the RM distribution observed across the source are patchy and show no obvious order.This suggests that we are observing a random RM component associated with the foreground material in the cluster in which the radio source is located and through which it is viewed, superposed on a more ordered RM component that arises in the immediate vicinity of the AGN jets. We interpret the transverse RM gradient as reflecting the systematic variations of the line-of-sight component of a helical or toroidal B field associated with the jets of 5C4.114. These results suggest that the helical field that arises due to the joint action of the rotation of the central black hole and its accretion disc and the jet outflow can survive to distances of thousands of parsec from the central engine.

The presence of a helical magnetic field threading the jet of an Active Galactic Nucleus (AGN) should give rise to a gradient in the observed Faraday rotation measure (RM) across the jet, due to the associated systematic change in the line-of-sight magnetic field. Reports of observations of transverse RM gradients across AGN jets have appeared in the literature starting from 2002, but concerns were raised about the resolution required for these gradients to be reliable, and there was a lack of a full understanding of the best approach to accurate estimation of the uncertainties of local RM values. These questions have now been resolved by recent Monte Carlo simulations carried out by various groups, enabling both a verification of previously published results and reliable analyses of new data. We consider here RM gradients across the jet structures of 15 AGN, some previously published in the refereed literature but without a correct and complete error analysis, and some published for the first time here, all of which have monotonic transverse RM gradients with significances of at least 3 sigma.

One possible origin for polarization structures across jets of Active Galactic Nuclei (AGNs) with a central "spine" of orthogonal magnetic field and a "sheath" of longitudinal magnetic field along one or both edges of the jet is the presence of a helical jet magnetic field. Simultaneous Very Long Baseline Array (VLBA) polarization observations of AGN displaying partial or full spine--sheath polarization structures were obtained at 4.6, 5.0, 7.9, 8.9, 12.9 and 15.4 GHz, in order to search for additional evidence for helical jet magnetic fields, such as transverse Faraday rotation gradients (due to the systematic change in the line-of-sight magnetic-field component across the jet). Results for eight sources displaying monotonic transverse Faraday rotation gradients with significances $\geq 3\sigma$ are presented here. Reversals in the directions of the transverse RM gradients with distance from the core or with time are detected in three of these AGNs. These can be interpreted as evidence for a nested helical magnetic field structure, with different directions for the azimuthal field component in the inner and outer regions of helical field. The results presented here support the idea that many spine--sheath polarization structures reflect the presence of helical magnetic fields being carried by these jets.

Mapping the relativistic jets emanating from AGN requires the use of a deconvolution algorithm to account for the effects of missing baseline spacings. The CLEAN algorithm is the most commonly used algorithm in VLBI imaging today and is suitable for imaging polarisation data. The Maximum Entropy Method (MEM) is presented as an alternative with some advantages over the CLEAN algorithm, including better spatial resolution and a more rigorous and unbiased approach to deconvolution. We have developed a MEM code suitable for deconvolving VLBI polarisation data. Monte Carlo simulations investigating the performance of CLEAN and the MEM code on a variety of source types are being carried out. Real polarisation (VLBA) data taken at multiple wavelengths have also been deconvolved using MEM, and several of the resulting polarisation and Faraday rotation maps are presented and discussed.

Both the emission properties and evolution of Active Galactic Nuclei (AGN) radio jets are dependent on the magnetic (B) fields that thread them. Faraday Rotation measurements are a very important way of investigating these B fields, and can provide information on the orientation and structure of the B field in the immediate vicinity of the jet; for example, a toroidal or helical B field component should give rise to a systematic gradient in the observed Faraday rotation across the jet, as well as characteristic intensity and polarization profiles. However, real observed radio images have finite resolution, usually expressed via convolution with a Gaussian beam whose size corresponds to the central lobe of the point source response function. This will tend to blur the transverse structure of the jet, raising the question of how well resolved a jet must be in the transverse direction in order to reliably detect transverse structure. We present the results of simulated Faraday rotation images designed to directly investigate the effect of finite resolution on observed transverse Faraday rotation measure structures.

We are investigating a new approach to modelling uncertainties in individual pixels of Very Long Baseline Interferometry (VLBI) intensity (Stokes I, Q, U) images. Comparison of distributions of our calculated uncertainties for model sources with the results of Monte Carlo simulations shows that our method correctly reproduces the overall level and pattern of uncertainties in intensity images for model sources that are not too compact. Refinement of the approach to better reproduce uncertainties in compact sources is being studied. This new approach should ultimately provide a means to correctly estimate pixel-based uncertainties in the Stokes parameters I, Q and U, including the effect of correlations between values in different pixels.

This paper describes a method of fitting total intensity and polarization profiles in VLBI images of astrophysical jets to profiles predicted by a theoretical model. As an example, the method is used to fit profiles of the jet in the Active Galactic Nucleus Mrk501 with profiles predicted by a model in which a cylindrical jet of synchrotron plasma is threaded by a magnetic field with helical and disordered components. This fitting yields model Stokes Q profiles that agree with the observed profiles to within the 1-2 \sigma uncertainties; the I model and observed profiles are overall not in such good agreement, with the model I profiles being generally more symmetrical than the observed profiles. Consistent fitting results are obtained for profiles derived from 6cm VLBI images at two distances from the core, and also for profiles obtained for different wavelengths at a single location in the VLBI jet. The most striking success of the model is its ability to reproduce the spine-sheath polarization structure observed across the jet. Using the derived viewing angle in the jet rest frame, \delta' approximately 83 degrees, together with a superluminal speed reported in the literature, \beta apparent = 3.3, yields a solution for the viewing angle and velocity of the jet in the observer's frame \delta degrees and \beta approximately 0.96. Although these results for Mrk501 must be considered tentative, the combined analysis of polarization profiles and apparent component speeds holds promise as a means of further elucidating the magnetic field structures and other parameters of parsec-scale AGN jets.

The idea that systematic Faraday Rotation gradients across the parsec-scale jets of AGNs can reveal the presence of helical magnetic (B) fields has been around since the early 1990s. These gradients are taken to be due to the systematic variation of the line of sight-B-field across the jet. We present here the parsec-scale Faraday Rotation distributions for the BL Lac objects 0716+714 and 1749+701, based on polarization data obtained with the Very Long Baseline Array (VLBA) at two wavelengths near each of the 2cm, 4cm and 6cm bands (0716+714) and at four wavelengths in the range 18-22 cm (1749+701). The Rotation Measure (RM) maps for both these sources indicate systematic gradients across their jets, as expected if these jets have helical B fields. The significance of these transverse RM gradients is > 3 sigma in all cases. We present the results of Monte Carlo simulations directly demonstrating the possibility of observing such transverse RM gradients even if the intrinsic jet structure is much narrower than the observing beam. We observe an intriguing new feature in these sources, a reversal in the direction of the gradient in the jet as compared to the gradient in the core region. This provides new evidence to support models in which field lines emerging from the central region of the accretion disk and closing in the outer region of the accretion disk are both "wound up" by the differential rotation of the disk. The net observed RM gradient will essentially be the sum effect of two regions of helical field, one nested inside the other. The direction of the net RM gradient will be determined by whether the inner or outer helix dominates the RM integrated through the jet, and RM gradient reversals will be observed if the inner and outer helical fields dominate in different regions of the jet. This potentially provides new insights about the overall configuration of the jet B fields.

Long-lived, stable jets are observed in a wide variety of systems, from protostars, through Galactic compact objects to active galactic nuclei (AGN). Magnetic fields play a central role in launching, accelerating, and collimating the jets through various media. The termination of jets in molecular clouds or the interstellar medium deposits enormous amounts of mechanical energy and momentum, and their interactions with the external medium, as well, in many cases, as the radiation processes by which they are observed, are intimately connected with the magnetic fields they carry. This review focuses on the properties and structures of magnetic fields in long-lived jets, from their launch from rotating magnetized young stars, black holes, and their accretion discs, to termination and beyond. We compare the results of theory, numerical simulations, and observations of these diverse systems and address similarities and differences between relativistic and non-relativistic jets in protostellar versus AGN systems. On the observational side, we focus primarily on jets driven by AGN because of the strong observational constraints on their magnetic field properties, and we discuss the links between the physics of these jets on all scales.

Global VLBI (EVN+VLBA) polarization observations at 5 and 8.4 GHz of ten high redshift (z > 3) quasars are presented. The core and jet brightness temperatures are found through modelling the self-calibrated uv-data with Gaussian components, which provide reliable estimates of the flux density and size of individual components. The observed high core brightness temperatures (median $T_{\rm b,\,core}=4\times10^{11}$ K) are consistent with Doppler boosted emission from a relativistic jet orientated close to the line-of-sight. This can also explain the dramatic jet bends observed for some of our sources since small intrinsic bends can be significantly amplified due to projection effects in a highly beamed relativistic jet. We also model-fit the polarized emission and, by taking the minimum angle separation between the model-fitted polarization angles at 5 and 8.4 GHz, we calculate the minimum inferred Faraday rotation measure (RM$_{\rm min}$) for each component. We also calculate the minimum intrinsic RM in the rest frame of the AGN (RM$_{\rm min}^{\rm intr}$ = RM$_{\rm min} (1+z)^2$), first subtracting the integrated (presumed foreground) RM in those cases where we felt we could do this reliably. The resulting mean core $|$RM$_{\rm min}^{\rm intr}|$ is 5580 rad m$^{-2}$, with a standard deviation of 3390 rad m$^{-2}$, for four high-z quasars for which we believe we could reliably remove the foreground RM. We find relatively steep core and jet spectral index values, with a median core spectral index of -0.3 and a median jet spectral index of -1.0. Comparing our results with RM observations of more nearby Active Galactic Nuclei at similar emitted frequencies does not provide any significant evidence for dependence of the quasar nuclear environment with redshift.

Previous research showed most BL Lac objects and some quasars have aligned VLBI-core and optical polarizations, although some of the AGNs also showed no obvious relationship between VLBI-core and optical polarization angles. This may indicate that some AGNs have co-spatial regions of optical and radio emission, while others do not. Another possibility is that some of the VLBI cores had Faraday rotations of several tens of thousand of rad/m^2, which were not properly fit using the three-frequency data due to n\pi ambiguities in the observed polarization angles, leading to incorrect subtraction of the effects of the core Faraday rotation, and so incorrect intrinsic radio polarization angles \chi_0. We obtained additional 12+15+22+24+43GHz plus optical observations for 8 of 40 AGNs previously considered. Our results indicate that, although some VLBI radio cores have comparatively high rotation measures, this alone cannot explain the misalignments found between the radio core and optical VLBI polarization angles, \Delta\chi = |\chi_opt - \chi_0|. Comparison between \Delta\chi and (i) the orientation of \chi_0 relative to the jet direction, (ii) the degree of polarization of the core, (iii) a depolarization factor, (iv) the core rotation measures and (v) the core magnetic fields 1pc from the jet base do not yield evidence for any correlations between these properties. There is some evidence that the maximum observed \Delta\chi tends to decrease as the core-region magnetic field increases, suggesting that misalignments in \Delta\chi could be associated in part with relatively low core magnetic fields. Thus, although the overall distribution of \Delta\chi for all 40 sources in our sample does show a significant peak at \Delta\chi\sim0, it remains unclear what distinguishes these AGN cores from those showing appreciable misalignment between optical and VLBI-core polarization position angles.

We present maps of the sign of the Faraday Rotation measure [sign(RM)] obtained from multi-frequency radio observations on the Very Long Baseline Array (VLBA). Many of the Active Galactic Nuclei (AGN) considered have B-field structures with a central "spine" of B-field orthogonal to the jet and/or a longitudinal B-field near one or both edges of the jet. This structure can plausibly be interpreted as being caused by a helical/toroidal jet magnetic field. Faraday Rotation is a rotation of the plane of polarization that occurs when the polarized radiation passes through a magnetized plasma. The sign of the RM is determined by the direction of the line-of-sight B-Field in the region causing the Faraday Rotation, and an ordered toroidal or helical magnetic field associated with an AGN jet will thus produce a distinctive bilateral distribution of positive and negative RMs across the jet. We present and discuss sign(RM) maps and their possible interpretation regarding the magnetic field geometries for several sources.

Both the emission properties and evolution of Active Galactic Nuclei (AGN) radio jets are dependent on the magnetic fields that thread them. Faraday Rotation gradients are a very important way of investigating these magnetic fields, and can provide information on the orientation and structure of the magnetic field in the immediate vicinity of the jet; for example, a toroidal or helical field component should give rise to a systematic gradient in the observed Faraday rotation across the jet, as well as characteristic intensity and polarization profiles. However, real observed radio images have finite resolution, usually expressed via convolution with a Gaussian beam whose size corresponds to the central lobe of the point source response function. This will tend to blur transverse structure in the jet profile, raising the question of how well resolved a jet must be in the transverse direction in order to reliably detect transverse structure associated with a helical jet magnetic field. We present results of simulated intensity, polarization and Faraday rotation images designed to directly and empirically investigate the effect of finite resolution on observed transverse jet structures.

The Maximum Entropy Method (MEM) for the deconvolution of radio interferometry images is mathematically well based and presents a number of advantages over the usual CLEAN deconvolution, such as appreciably higher resolution. The application of MEM for polarisation imaging remains relatively little studied. CLEAN and MEM intensity and polarisation techniques are discussed in application to recently obtained 18cm VLBA polarisation data for a sample of Active Galactic Nuclei.

We discuss the opacity in the core regions of active galactic nuclei observed with Very Long Baseline Interferometry (VLBI), and describe a new method for deriving the frequency-dependent shifts of the VLBI core from the frequency-dependent time lags of flares observed with single-dish observations. Application of the method to the core shifts of the quasar 3C 345 shows a very good agreement between the core shifts directly measured from VLBI observations and derived from flares in the total flux density using the proposed method. The frequency-dependent time lags of flares can be used to derive physical parameters of the jets, such as distance from the VLBI core to the base of the jet and the magnetic fields in the core region. Our estimates for 3C 345 indicate core magnetic fields ~0.1 G and magnetic field at 1 pc ~0.4 G.

A radio study of a carefully selected sample of 20 Seyfert galaxies that are matched in orientation-independent parameters, which are measures of intrinsic active galactic nuclei (AGN) power and host galaxy properties is presented to test the predictions of the unified scheme hypothesis. Our sample sources have core flux densities greater than 8 mJy at 5 GHz on arcsec-scales due to the feasibility requirements. These simultaneous pc-scale and kpc-scale radio observations reveal (i) that Seyfert 1 and Seyfert 2 galaxies have equal tendency to show compact radio structures on mas-scales, (ii) the distributions of pc-scale and kpc-scale radio luminosities are similar for both Seyfert 1 and Seyfert 2 galaxies, (iii) no evidence for relativistic beaming in Seyfert galaxies, (iv) similar distributions of source spectral indices in spite of the fact that Seyferts show nuclear radio flux density variations, and (v) the distributions of projected linear size for Seyfert 1 and Seyfert 2 galaxies are not significantly different as would be expected in the unified scheme. The latter could be mainly due to a relatively large spread in the intrinsic sizes. We also find that a starburst alone cannot power these radio sources. Finally, an analysis of the kpc-scale radio properties of the CfA Seyfert galaxy sample shows results consistent with the predictions of the unified scheme.

We present multi-frequency radio observations from the Very Long Baseline Array (VLBA) of selected AGN that seem to have a B-field structure with a central "spine" of B-field orthogonal to the jet and a longitudinal B-field near one or both edges of the jet. Two explanations for this structure have been discussed in the literature: shocks making the central orthogonal field combined with a jet-medium interaction causing the longitudinal "sheath", or both components produced by a helical jet magnetic field. One way to investigate this structure is to look for gradients in the Faraday Rotation across the jet. We will discuss results for 0333+321, 1150+812 and 2037+511 providing evidence for the latter picture.

We are in the process of obtaining VLBA polarisation data for the 135 MOJAVE-I Active Galactic Nuclei at four wavelengths in the 18-22cm band. These observations will enable studies of the evolution of the intensity and magnetic-field structures of these AGN jets as they propagate from parsec to kiloparsec scales, as well as studies of the thermal plasma present in the vicinity of the jets on these scales, manifest via Faraday rotation. A wide range of other multi-wavelength studies can also be carried out using these data. Preliminary results for selected sources from the first 3 of 9 observing sessions will be presented. We aim to have 18-cm intensity and polarisation images available via the MOJAVE website within 18 months after the last observing session.

We present polarization measurements at 8.4, 22, and 43 GHz made with the VLA of a complete sample of extragalactic sources stronger than 1 Jy in the 5-year WMAP catalogue and with declinations north of -34 degrees. The observations were motivated by the need to know the polarization properties of radio sources at frequencies of tens of GHz in order to subtract polarized foregrounds for future sensitive Cosmic Microwave Background (CMB) experiments. The total intensity and polarization measurements are generally consistent with comparable VLA calibration measurements for less-variable sources, and within a similar range to WMAP fluxes for unresolved sources. A further paper will present correlations between measured parameters and derive implications for CMB measurements.

We use observations of six "blazars" with the Very Long Baseline Array (VLBA), at eight frequencies (4.6, 5.1, 7.9, 8.9, 12.9, 15.4, 22.2, 43.1 GHz), to investigate the frequency-dependent position of their VLBI cores ("core-shift") and their overall jet spectral distribution. By cross-correlating the optically thin jet emission, we are able to accurately align the multi-frequency images of three of the jets (1418+546, 2007+777, 2200+420), whose core-shifts and spectra we find consistent with the equipartition regime of the Blandford & Konigl conical jet model, where the position of the radio core from the base of the jet follows $r_{core}\propto\nu^{-1}$. For the jet of 0954+658, we align the higher frequency images using our lower frequency measurements assuming equipartition in the radio core from 4.6-43 GHz. The jet emission of the other two sources in our sample (1156+295, 1749+096) is too sparse for our alignment technique to work. Using our measured core-shifts, we calculate equipartition magnetic field strengths of the order of 10's to 100's of mG in the radio cores of these four AGN from 4.6-43 GHz. Extrapolating our results back to the accretion disk and black hole jet-launching distances, we find magnetic field strengths consistent with those expected from theoretical models of magnetically powered jets.

The origin of cosmic magnetic (B) fields remains an open question. It is generally believed that very weak primordial B fields are amplified by dynamo processes, but it appears unlikely that the amplification proceeds fast enough to account for the fields presently observed in galaxies and galaxy clusters. In an alternative scenario, cosmic B fields are generated near the inner edges of accretion disks in Active Galactic Nuclei (AGNs) by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with photons. While dynamo processes show no preference for the polarity of the (presumably random) seed field that they amplify, this alternative mechanism uniquely relates the polarity of the poloidal B field to the angular velocity of the accretion disk, resulting in a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry came about by chance being less than 1%. This lends support to the hypothesis that the Universe is seeded by B fields that are generated in AGN via this mechanism and subsequently injected into intergalactic space by the jet outflows.

Several multi-frequency polarization studies have shown the presence of systematic Faraday Rotation gradients across the parsec-scale jets of Active Galactic Nuclei (AGN), taken to be due to the systematic variation of the line-of-sight component of a helical magnetic (B) field across the jet. Other studies have confirmed the presence and sense of these gradients in several sources, thus providing evidence that these gradients persist over time and over large distances from the core. However, we find surprising new evidence for a reversal in the direction of the Faraday Rotation gradient across the jet of B1803+784, for which multi-frequency polarization observations are available at four epochs. At our three epochs and the epoch of Zavala & Taylor (2003), we observe transverse Rotation Measure (RM) gradients across the jet, consistent with the presence of a helical magnetic field wrapped around the jet. However, we also observe a "flip" in the direction of the gradient between June 2000 and August 2002. Although the origins of this phenomena are not entirely clear, possibly explanations include (i) the sense of rotation of the central supermassive black hole and accretion disc has remained the same, but the dominant magnetic pole facing the Earth has changed from North to South; (ii) a change in the direction of the azimuthal B field component as a result of torsional oscillations of the jet; and (iii) a change in the relative contributions to the observed rotation measures of the "inner" and "outer" helical fields in a magnetic-tower model. Although we cannot entirely rule out the possibility that the observed changes in the RM distribution are associated instead with changes in the thermal-electron distribution in the vicinity of the jet, we argue that this explanation is unlikely.

Parsec-scale multi-wavelength VLBA polarization observations can be used to study the magnetic-field structures of Active Galactic Nuclei (AGN) based on Faraday Rotation (FR) gradients. A number of transverse FR gradients have been found, and interpreted as corresponding to helical magnetic fields wrapped around the jets; the gradients reflect the systematic change in the line-of-sight component of a toroidal or helical magnetic field across the jet (e.g Gabuzda, Murray & Cronin 2004). Our observations of a sample of BL Lac objects at six wavelengths near 2, 4 and 6 cm have also revealed a previously undetected phenomena: these transverse gradients sometimes change their direction with distance from the core. We have observed this behaviour in at least five sources, which display gradients in their VLBI core region opposite to those in the jet. We suggest that this may be linked to magnetic tower models. In magnetic tower models, the field lines go outward with the jet and return and close in the accretion disk (or vice versa); differential rotation of the accretion disk winds up the inner and outer field lines into two helices (the inner helix "nested" in the outer helix). The total observed FR gradient is a sum of the effect of these two helical fields. It may be that gradients detected relatively far from the core correspond to the outer helix, while gradients detected in the core region correspond to dominance of the inner helix. This provides tentative evidence for the unification of helical magnetic fields and magnetic tower models, which could provide crucial new information for understanding AGN jets. Further VLBI studies with resolution sufficient to reliably detect these gradients in the cm-wavelength core and inner jet will be important for further investigations of this phenomena.

The parsec-scale Faraday rotation measure (RM) distribution of six "blazars" is investigated using multi-frequency (4.6--43 GHz) polarization observations taken on 2006 July 2 with the VLBA. Analysis of the RM provides the direction of the line-of-sight (LoS) magnetic field component, as well as the intrinsic 2-D polarization distribution on the plane of the sky. Our results show that the magnitude of the core RM increases systematically with frequency, and is well described by a power-law, where |RM_{core}| \propto \nu^a. Our measured values of $a$ vary from 0.9 to 3.8, providing information on the assumed power-law fall-off in the electron density with distance from the central engine for each source. RM gradients were detected across the jets of three sources, supporting the presence of helical magnetic fields in a sheath or boundary layer surrounding their jets. We find a bi-modal distribution of the intrinsic jet polarization orientation; either aligned or orthogonal to the jet direction. A helical magnetic field geometry can neatly explain both the bi-model distribution of the jet polarization orientation and the ordered polarization structure detected on these scales. In half the sources, we find that the core RM changes sign with distance from the central engine. We provide an explanation for this by considering a boundary layer of Faraday rotating material threaded by a helical magnetic field, where bends in the relativistic jet or accelerating/decelerating flows give rise to changes in the dominant LoS components of the magnetic field, which in turn gives rise to different signs of the RM. (abridged)

Although the continua of radio-loud Active Galactic Nuclei (AGN) are typically dominated by synchrotron radiation over virtually the entire spectrum, it is not clear whether the radio and higher-frequency emission originate in the same or different parts of the jet. Several different radio--optical correlations based on polarization data have been found recently, suggesting that the optical and radio polarization may be closely related, and that the corresponding emission regions may be cospatial (Gabuzda et. al2006, Jorstad et al. 2007, D'Arcangelo et al. 2007) Our joint analysis of optical and VLBA polarization data for a sample of about 40 AGNs shows that, after correction for the inferred VLBA core Faraday rotations, most BL Lac objects and some quasars have aligned VLBA-core and optical polarizations, although many quasars also show no obvious relationship between their VLBA-core and optical polarization angles. This may indicate that not all AGNs have cospatial regions of optical and radio emission in their jets. However, another possibility is that some of the 7mm-2cm VLBA cores have Faraday rotations of the order of several tens of thousand of rad/m^2, which were not properly fit by our three-frequency data due to n*pi ambiguities in the observed polarization angles, leading to incorrect subtraction of the effects of the core Faraday rotation, and so incorrect "zero-wavelength" radio polarization angles. The possibility of such high core Faraday rotations is supported by the results of the parsec-scale Faraday-rotation studies of Zavala & Taylor (2004) and Jorstad et al. (2007).

We present the results of a parsec-scale polarization study of three FRI radio galaxies - 3C66B, 3C78 and 3C264 - obtained with the Very Long Baseline Array at 5, 8 and 15 GHz. Parsec-scale polarization has been detected in a large number of beamed radio-loud active galactic nuclei, but in only a handful of the relatively unbeamed radio galaxies. We report here the detection of parsec-scale polarization at one or more frequencies in all three FRI galaxies studied. We detect Faraday rotation measures of the order of a few hundred rad/m^2 in the nuclear jet regions of 3C78 and 3C264. In 3C66B polarization was detected at 8 GHz only. A transverse rotation measure gradient is observed across the jet of 3C78. The inner-jet magnetic field, corrected for Faraday rotation, is found to be aligned along the jet in both 3C78 and 3C264, although the field becomes orthogonal further from the core in 3C78. The RM values in 3C78 and 3C264 are similar to those previously observed in nearby radio galaxies. The transverse RM gradient in 3C78, the increase in the degree of polarization at the jet edge, the large rotation in the polarization angles due to Faraday rotation and the low depolarization between frequencies, suggests that a layer surrounding the jet with a sufficient number of thermal electrons and threaded by a toroidal or helical magnetic field is a good candidate for the Faraday rotating medium. This suggestion is tentatively supported by Hubble Space Telescope optical polarimetry but needs to be examined in a greater number of sources.

Accurate measurement of the frequency-dependent shift of the self-absorbed radio core is required for multi-frequency analysis of VLBI data since absolute positional information is lost as a result of phase self-calibration. We use the cross-correlation technique of Croke & Gabuzda (2008) on the optically thin jet emission to align our VLBA images. Our results are consistent with those obtained from the phase-referencing method, as well as alignment by model-fitted optically thin jet components. Physical parameters of the compact jet regions, such as the magnetic field strength (B_core) and the distance of the radio core to the jet origin (r_core), can be calculated from these measurements. For the source Mrk 501, we find a magnetic field strength of 0.15\pm0.04 G in the 8.4-GHz core at a distance of 0.8\pm0.2 pc from the base of the jet. By extrapolating our 4.6 to 15.4 GHz results for BL Lac (2200+420), we estimate magnetic field strengths of the order of 1 G in the millimetre VLBI core. Using our core-shift measurement between 1.6 and 4.8 GHz for 1803+784, we find B_core(4.8 GHz) = 0.11\pm0.02 G and r_core(4.8 GHz) = 20\pm5 pc. The phase-referencing observations of this source at 8.4 and 43 GHz by Jim\'enez-Monferrer et al. (2008) imply B_core(43 GHz) = 1.0\pm0.4 G and r_core(43 GHz) = 2.0\pm0.9 pc.

Many important techniques for investigating the properties of extragalactic radio sources, such as spectral-index and rotation-measure mapping, involve the comparison of images at two or more frequencies. In the case of radio interferometric data, this can be done by comparing the CLEAN maps obtained at the different frequencies. However, intrinsic differences in images due to the frequency dependence of the radio emission can be distorted by additional differences that arise due to source variability (if the data to be compared is obtained at different times), image misalignment, and the frequency dependence of the sensitivity to weak emission and the angular resolution provided by the observations (the resolution of an interferometer depends on the lengths of its baselines in units of the observing wavelength). These effects must be corrected for as best as possible before multi-frequency data comparison techniques can be applied. We consider the origins for the afore-mentioned factors, outline the standard techniques used to overcome these difficulties, and describe in detail a technique developed by us, based on the cross-correlation technique widely used in other fields, to correct for misalignments between maps at different frequencies.

We present the results of parsec-scale circular polarization (CP) measurements based on VLBA data for a number of radio-bright, core-dominated AGN obtained simultaneously at 15, 22 and 43 GHz. The degrees of CP $m_c$ for the VLBI core region at 15 GHz are similar to values reported earlier at this wavelength, with typical values of a few tenths of a percent.We find that $m_c$ as often rises as falls with increasing frequency between 15 and 22 GHz, while $m_c$ at 43 GHz is in all cases higher than at 22 and 15 GHz. This behaviour seems contrary to expectations, since the degree of CP from both synchrotron radiation and Faraday conversion of linear to circular polarization should {\em decrease} towards higher frequencies if the source is homogeneous. The increase in $m_c$ at 43 GHz may be due to the presence of regions of both positive and negative CP with different frequency dependences on small scales within the core region; alternatively, it may be associated with the intrinsic inhomogeneity of a Blandford-K\"onigl-like jet. In several objects, the detected CO appears to be near, but not coincident with the core, although further observations are needed to confirm this. We find several cases of changes in sign with frequency, most often between 22 and 43 GHz. We find tentative evidence for transverse structure in the CP of 1055+018 and 1334$-$127 that is consistent with their being generated by either the synchrotron mechanism or Faraday conversion in a helical magnetic field. Our results confirm the earlier finding that the sign of the CP at a given observing frequency is generally consistent across epochs separated by several years or more, suggesting stability of the magnetic field orientation in the innermost jets.

We are using "broadband" (4.6 to 43 GHz) multi-frequency VLBA polarization observations of compact AGN to investigate the 3-D structure of their jet magnetic (B) fields. Observing at several frequencies, separated by short and long intervals, enables reliable determination of the distribution of Faraday Rotation, and thereby the intrinsic B field structure. Transverse Rotation Measure (RM) gradients were detected in the jets of 0954+658 and 1418+546, providing evidence for the presence of a helical B field surrounding the jet. The RM in the core regions of 2200+420 (BL Lac), 0954+658 and 1418+546 display different signs in different frequency-intervals (on different spatial scales); we suggest an explanation for this in terms of modest bends in a helical B field surrounding their jets. In future, polarization observations with a combination of VSOP-2 at 8, 22 and 43 GHz and ground arrays at frequencies with corresponding resolution will help map out the distributions of Faraday rotation, spectral index and the 3-D B field structure both across the jet and closer to the central engine, providing strong constraints for any jet B field models.

Multi-frequency (4.6, 5, 5.5, 8, 8.8, 13, 15, 22 & 43 GHz) polarization observations of 6 "blazars" were obtained on the American Very Long Baseline Array (VLBA) over a 24-hr period on 2 July 2006. Observing at several frequencies, separated by short and long intervals, enabled reliable determination of the distribution of Faraday Rotation on a range of scales. In all cases the magnitude of the RM increases in the higher frequency observations, implying that the electron density and/or magnetic field strength is increasing as we get closer to the central engine. After correcting for Faraday rotation, the polarization orientation in the jet is either parallel or perpendicular to the jet direction. A transverse Rotation Measure (RM) gradient was detected in the jet of 0954+658, providing evidence for the presence of a helical magnetic field surrounding the jet. For three of the sources (0954+658, 1418+546, 2200+420), the sign of the RM in the core region changes in different frequency-intervals, indicating that the line-of-sight component of the magnetic field is changing with distance from the base of the jet. We suggest an explanation for this in terms of bends in a relativistic jet surrounded by a helical magnetic field; where there is no clear evidence for pc-scale bends, the same effect can be explained by an accelerating/decelerating jet.

Homan & Lister (2006) have recently published circular-polarization (CP) detections for 34 objects in the MOJAVE sample - a set of bright, compact AGN being monitored by the Very Long Baseline Array at 15 GHz. We report the detection of 15-GHz parsec-scale CP in two more AGN (3C345 and 2231+114), and confirm the MOJAVE detection of CP in 1633+382. It is generally believed that the most likely mechanism for the generation of this CP is Faraday conversion of linear polarization to CP. A helical jet magnetic-field (B-field) geometry can facilitate this process - linearly polarized emission from the far side of the jet is converted to CP as it passes through the magnetised plasma at the front side of the jet on its way toward the observer. In this case, the sign of the generated CP is essentially determined by the pitch angle and helicity of the helical B field. We have determined the pitch-angle regimes and helicities of the helical jet B fields in 8 AGN for which parsec-scale CP has been detected, and used them to predict the expected CP signs for these AGN if the CP is generated via conversion in these helical fields. We have obtained the intriguing result that our predictions agree with the observed signs in all eight cases, provided that the longitudinal B-field components in the jets correspond to South magnetic poles. This clearly non-random pattern demonstrates that the observed CP in AGN is directly associated with the presence of helical jet B fields. These results suggest that helical B fields are ubiquitous in AGN jets.

We present very long baseline interferometry polarization images of an X-ray selected sample of BL Lacertae objects belonging to the first High Energy Astronomy Observatory (HEAO-1) and the ROSAT-Green Bank (RGB) surveys. These are primarily high-energy-peaked BL Lacs (HBLs) and exhibit core-jet radio morphologies on pc-scales. They show moderately polarized jet components, similar to those of low-energy-peaked BL Lacs (LBLs). The fractional polarization in the unresolved cores of the HBLs is, on average, lower than in the LBLs, while the fractional polarizations in the pc-scale jets of HBLs and LBLs are comparable. However a difference is observed in the orientation of the inferred jet magnetic fields -- while LBL jets are well-known to preferentially exhibit transverse magnetic fields, the HBL jets tend to display longitudinal magnetic fields. Although a `spine-sheath' jet velocity structure, along with larger viewing angles for HBLs could produce the observed magnetic field configuration, differences in other properties of LBLs and HBLs, such as their total radio power, cannot be fully reconciled with the different-angle scenario alone. Instead it appears that LBLs and HBLs differ intrinsically, perhaps in the spin rates of their central black holes.

Preliminary VLBA polarisation results on 6 ``blazars'' from 6.5-cm to 7-mm are presented here. Observing at several different wavelengths, separated by short and long intervals, enabled reliable information about the magnetic (B) field structure to be obtained and for the effect of Faraday Rotation to be determined and corrected. For all sources the magnitude of the core Rotation Measure (RM) derived from the shorter wavelength data was greater than that derived from the longer wavelength data, consistent with a higher electron density and/or B-field strength closer to the central engine. A transverse RM gradient was detected in the jet of 0954+658, providing evidence for the presence of a helical B-field surrounding the jet. The RM in the core region of 2200+420 (BL Lac) displays sign changes in different wavelength intervals (on different spatial scales); we suggest an explanation for this in terms of modest bends in a helical B-field surrounding the jet.

Several multi-frequency polarization studies have shown the presence of systematic Faraday Rotation gradients across the parsec-scale jets of Active Galactic Nuclei (AGN), taken to be due to the systematic variation of the line-of-sight component of a helical magnetic field across the jet. Other studies have confirmed the presence and sense of these gradients in several sources, thus providing evidence that these gradients persist over time and over large distances from the core. However, we find surprising new evidence for a reversal in the direction of the Faraday Rotation gradient across the jet of 1803+784, for which multi-wavelength polarization observations are available at four epochs. At all four epochs, we observe transverse Rotation Measure (RM) gradients across the jet, consistent with the presence of a helical magnetic field wrapped around the jet. However, we also observe a ''flip'' in the direction of the gradient between June 2000 and August 2002. Although the origins of this phenomena are not understood, one way to interpret this change is if the sense of rotation of the central supermassive black hole and accretion disc has remained the same, but the dominant magnetic pole facing the Earth has changed from North to South.

Previously, a multi-wavelength (2cm, 4cm and 6cm) polarization study by Gabuzda, Murray and Cronin (2004) showed systematic Faraday Rotation gradients across the parsec-scale jets of several BL Lac objects, interpreted as evidence for helical magnetic (B) fields - the gradients were taken to be due to the systematic variation of the line-of-sight B field across the jet. We present here new results for the parsec-scale Faraday Rotation distributions for eight additional BL Lac objects, based on polarization data obtained with the Very Long Baseline Array (VLBA) at two wavelengths near each of the 2cm, 4cm and 6cm bands. The Rotation Measure (RM) maps for all these sources indicate gradients across their jets, as expected if these jets have helical B fields. Such gradients are also detected in the cores of several sources. This provides evidence that these gradients are present in appreciable regions of the jets and are not isolated phenomena. We also observe reversals in the RM gradient in the core region as compared to the gradient in the jet in at least three sources.

Relatively few Very Long Baseline Interferometry (VLBI) polarization observations have been carried out at 18 cm. The importance of such observations lies in their ability to reveal information about the jet magnetic ({\bf B}) field structure and the environment of the jet on scales intermediate between those probed by higher-frequency VLBI and connected-element interferometers such as the Very Large Array. We have obtained polarization observations of 34 BL Lac objects with the Very Long Baseline Array (VLBA), at 4 separate wavelengths in the 18-20 cm band. The 18-cm jets typically extend to tens of parsecs. In some cases, the decaparsec jet is a continuation of the jet on smaller scales, while in others, we see appreciable bending. We have constructed Faraday rotation-measure maps and used them to study the jet {\bf B} field structures and distribution of thermal plasma around the jets. The Faraday rotation is typically large at these wavelengths, and knowledge of the rotation-measure distribution is essential to derive the {\bf B} field structures of the jets. The high sensitivity of these observations to Faraday rotation makes them an effective tool for studies of possible interactions between the jets and the media through which they propagate.

The combined data of the University of Michigan Radio Astronomy Observatory and Metsahovi Radio Observatory provide us with radio light curves for Active Galactic Nuclei monitored by both observatories from 4.8 to 37 GHz covering time intervals up to ~25 years. We consider here such composite light curves for four gamma-ray blazars that have been nearly continuously monitored at both observatories: 0458-020, 0528+134, 1730-130 and 2230+114. We have decomposed the most prominent outbursts in the light curves of these four blazars into individual components using Gaussian model fitting, and estimated the epochs, amplitudes, and half-widths of these components as functions of frequency. We attempt to distinguish "core outbursts", which show frequency-dependent time delays and are associated with brightening of the core, from "jet outbursts", which appear nearly synchronously at all frequencies and are accompanied by the emergence of new jet components and their subsequent evolution. Available 43 GHz VLBA images allow us to identify only one pure core outburst (in 2230+114) and one pure jet outburst (0458-020). Most of the outbursts analyzed are mixed, in the sense that they display frequency-dependent time delays (i.e., they are optically thick) and are associated with the eventual emergence of new jet components. The maxima of the jet and mixed outbursts probably correspond to epochs when newly ejected components become fully optically thin. These epochs are also marked by a significant increase in the angular velocities of the ejected components. There is evidence that the outbursts in 2230+114 repeat every 8.0+-0.3 years, with the positions of individual sub-outbursts being preserved from one quasi-periodic eight-year cycle to another, even though their amplitudes vary by more than a factor of two.

The TeV-emitting BL Lac object Mkn 421 was observed with very long baseline interferometry (VLBI) at three closely-spaced epochs one-month apart in March-April 1998. The source was also monitored at very-high gamma-ray energies (TeV measurements) during the same period in an attempt to search for correlations between TeV variability and the evolution of the radio morphology on parsec scales. While the VLBI maps show no temporal changes in the Mkn 421 VLBI jet, there is strong evidence of complex variability in both the total and polarized fluxes of the VLBI core of Mkn 421 and in its spectrum over the two-month span of our data. The high-energy measurements indicate that the overall TeV activity of the source was rising during this period, with a gamma-ray flare detected just three days prior to our second VLBI observing run. Although no firm correlation can be established, our data suggest that the two phenomena (TeV activity and VLBI core variability) are connected, with the VLBI core at 22 GHz being the self-absorbed radio counterpart of synchrotron self-Compton (SSC) emission at high energies. Based on the size of the VLBI core, we could derive an upper limit of 0.1 pc (3 x 10**17 cm) for the projected size of the SSC zone. This determination is the first model-free estimate of the size of the gamma-ray emitting region in a blazar.

We present polarimetric 5 GHz to 43 GHz VLBI observations of the BL Lacertae object PKS 0735+178, spanning March 1996 to May 2000. Comparison with previous and later observations suggests that the overall kinematic and structural properties of the jet are greatly influenced by its activity. Time intervals of enhanced activity, as reported before 1993 and after 2000 by other studies, are followed by highly superluminal motion along a rectilinear jet. In contrast the less active state in which we performed our observations, shows subluminal or slow superluminal jet features propagating through a twisted jet with two sharp bends of about 90 deg. within the innermost three-milliarcsecond jet structure. Proper motion estimates from the data presented here allow us to constrain the jet viewing angle to values < 9 deg., and the bulk Lorentz factor to be between 2 and 4.

We present images showing the first detections of polarization on parsec scales in the nuclei of four Fanaroff-Riley type I (low-luminosity) radio galaxies. Observations with Very Long Baseline Interferometry at \lambda 3.6cm reveal the presence of ordered magnetic fields within ~1650 Schwarzchild radii of the putative central supermassive black hole. The relatively high fractional polarization in the pc-scale jets of these galaxies is consistent with the standard scheme unifying low-luminosity radio galaxies with BL Lac objects. This result also suggests that these radio galaxies lack the obscuring tori that apparently depolarize the nuclear emission in the more powerful FRII radio galaxies, and that their supermassive blackholes are poorly fed and/or inefficient radiators.

In 1999, the highly compact and variable BL Lac object AO 0235+164 was identified as the highest brightness temperature active galactic nucleus observed with the VLBI Space Observatory Programme (VSOP), with T_B > 5.8 x 10^{13} K. The sub-milliarcsecond radio structure of this source has been studied with dual-frequency (1.6 and 5 GHz), polarization-sensitive VSOP observations during 2001 and 2002. Here we present the results of this monitoring campaign. At the time of these observations, the source was weakly polarized and characterized by a radio core that is clearly resolved on space-ground baselines.

Frequency-dependent time lags for strong outbursts in four gamma-blazars are determined. The time lags for two adjacent outbursts in 2230+114 are correlated with the outburst amplitudes. There is evidence that bright outbursts in 2230+114 appear with a quasi-period of (8.0 +- 0.3) yr.

Results of multi-frequency VLBA polarization imaging for sources showing evidence of quasi-periodic variations in their total flux density are discussed. Possible origins of the observed behavior are proposed.

The parsec-scale total intensity structures of BL Lac objects share many characteristics with quasars: e.g. well-collimated jets and superluminal motion. However, the jets of BL Lac objects appear to fade much more quickly than those in quasars and on VLA scales many BL Lacs have structures comparable to those of low luminosity (FR I) radio galaxies, in which kpc-scale flow speeds are non-relativistic. The region between parsec and kiloparsec scales is therefore one of transition in BL Lac objects. We have carried out VLBI observations of a small sample of BL Lac objects at frequencies between 2.3 GHz and 327 MHz in order to investigate this transition region. Preliminary results for two of these sources are presented here. Coherent structures are detected in these sources out to a distance of several tens of parsecs. We find evidence for significant changes in jet structure on scales of a few tens of milliarcseconds, both in terms of the orientation and collimation of the jets.

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

We have obtained mas-scale resolution very long baseline interferometry (VLBI) images of a sample of Seyfert 1 and Seyfert 2 galaxies at 5 GHz (wavelength, $\lambda$ = 6 cm). The Seyferts of the two types were chosen to be matched in several orientation-independent properties, primarily in order to rigorously test predictions of the unified scheme. We detected all the 15 objects that we observed. In this paper we describe the observations and data reduction procedures, and present the VLBI radio images as well as simultaneous Very Large Array images that we obtained for these 15 Seyferts.

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

Evidence has been mounting that many of the transverse jet B fields observed in BL Lac objects on parsec scales represent the dominant vvroidal compoent of the intrinsic jet B fields. Such fields could come about, for example, as a result of the "winding up" of an initial "seed" field with a significant longitudinal component by the rotation of the central accreting object. If this is the case, this should give rise to gradients in the rotation measure (RM) across the jets, due to the systematic change in the line-of-sight component of the jet B field. We present evidence for transverse RM gradients in four BL Lac objects, strengthening arguments that the jets of these objects do indeed have toroidal or helical B fields. This underlines the view of the jets as fundamentally electromagnetic structures, and suggests that they may well carry non-zero currents. It also provides a natural means to collimate the jets.

If the transverse jet B fields observed in BL Lac objects on parsec scales represent the toroidal component of the jet B fields, this should give rise to rotation-measure gradients across the jets, due to the systematic change in the line-of-sight component of the jet B field. We have found evidence for such rotation-measure gradients in one-third to one-half of the 1 Jy BL Lac objects we have analyzed. We present and discuss these new results, together with some of their implications for our understanding of the parsec-scale jets of AGN.

Evidence is mounting that many of the transverse jet B fields observed in BL Lac objects on parsec scales represent the dominant toroidal component of the intrinsic jet B fields. If this is the case, this may give rise to rotation-measure (RM) gradients across the jets, due to the systematic change in the line-of-sight component of the jet B field. We have found evidence for such RM gradients in several BL Lac objects. We discuss these new results, together with some of their implications for our understanding of the pc-scale jets of AGN.

We have observed the five radio galaxies NGC1052, NGC4261, Hydra A, Centaurus A and Cygnus A with polarimetric VLBI at 15 GHz to look for possible Faraday rotation in foreground absorbers. Only Cygnus A was found to be polarized at the level of 0.4 %, the other four objects were found to be unpolarized (<0.3 %). We discuss implications for magnetic field structures and physical conditions in the absorbers imposed by the absence of polarized emission at this high frequency.

The results of VLBI total intensity (I) and linear polarization (P) observations at 6cm are presented for ten radio bright BL Lacertae objects. These images complete first-epoch polarization observations for the 1-Jy sample of northern BL Lacertae objects defined by Kuhr and Schmidt. Estimates of superluminal speeds are presented for several sources, bringing the total number of sources in the sample for which such estimates are available to 16. Second epoch observations currently being reduced should yield speed estimates for VLBI features in essentially all the sources in the sample. The jet magnetic fields of these BL Lacertae objects are usually transverse to the local jet direction, but a sizeable minority (about 30%) have VLBI jet components with longitudinal magnetic fields. This may suggest that the conditions in the VLBI jets of BL Lacertae objects are favorable for the formation of relativistic shocks; alternatively, it may be that the toroidal component of the intrinsic jet magnetic field is characteristically dominant in these sources.

We present VLA and first-epoch VLBA observations that are part of a program to study the parsec-scale radio structure of a sample of fifteen high-energy-peaked BL Lacs (HBLs). The sample was chosen to span the range of logarithmic X-ray to radio flux ratios observed in HBLs. As this is only the first epoch of observations, proper motions of jet components are not yet available; thus we consider only the structure and alignment of the parsec- and kiloparsec-scale jets. Like most low-energy-peaked BL Lacs (LBLs), our HBL sample shows parsec-scale, core-jet morphologies and compact, complex kiloparsec-scale morphologies. Some objects also show evidence for bending of the jet 10-20pc from the core, suggesting interaction of the jet with the surrounding medium. Whereas LBLs show a wide distribution of parsec- to kpc-scale jet misalignment angles, there is weak evidence that the jets in HBLs are more well-aligned, suggesting that HBL jets are either intrinsically straighter or are seen further off-axis than LBL jets.

Very Long Baseline Polarimetric observations of the Broad Line Radio galaxy 3C111 performed in July and September of 1996 at 8 and 43 GHz reveal rapidly evolving parsec-scale radio structure after a large millimetre outburst. The B-field geometry is not simple. We present a first analysis of possible Faraday and optical depth effects based on a comparison of the polarization images for the two frequencies.

We compare the systematics of the magnetic field geometry in the "red" low-energy peaked BL Lacs (LBLs) and "blue" high-energy peaked BL Lacs (HBLs) using VLBI polarimetric images. The LBLs are primarily "radio--selected" BL Lacs and the HBLs are primarily "X-ray selected". In contrast to the LBLs, which show predominantly transverse jet magnetic fields, the HBLs show predominantly longitudinal fields. Thus, while the SED peaks of core-dominated quasars, LBLs and HBLs form a sequence of increasing frequency, the magnetic field geometry does not follow an analogous sequence. We briefly investigate possible connections between the observed parsec-scale magnetic field structures and circular polarization measurements in the literature on various spatial scales.

We study two aspects of the differences between Active Galactic Nuclei (AGN) of the two Fanaroff-Riley types in order to investigate the causes of the F-R divide. We (a) contrast the properties of the optical cores in beamed and unbeamed AGN of the two types, incorporating Hubble Space Telescope measurements of the unbeamed objects, and (b) contrast the nuclear magnetic field geometry of the beamed AGN of the two types using Very Long Baseline Polarimetry.

We present multi-epoch 8.4 and 43 GHz Very Long Baseline Array images of the BL Lac object 0735+178. The images confirm the presence of a twisted jet with two sharp apparent bends of 90$^{\circ}$ within two milliarcseconds of the core, resembling a helix in projection. The observed twisted geometry could be the result of precession of the jet inlet, but is more likely produced by pressure gradients in the external medium through which the jet propagates. Quasi-stationary components are observed at the locations of the 90$^{\circ}$ bends, possibly produced by differential Doppler boosting. Identification of components across epochs, since the earliest VLBI observations of this source in 1979.2, proves difficult due to the sometimes large time gaps between observations. One possible identification suggests the existence of superluminal components following non--ballistic trajectories with velocities up to $11.6\pm 0.6 h_{65}^{-1} c$. However, in images obtained after mid-1995, components show a remarkable tendency to cluster near several jet positions, suggesting a different scenario in which components have remained nearly stationary in time at least since mid-1995. Comparison with the earlier published data, covering more than 19 years of observations, suggests a striking qualitative change in the jet trajectory sometime between mid-1992 and mid-1995, with the twisted jet structure with stationary components becoming apparent only at the later epochs. This would require a re-evaluation of the physical parameters estimated for 0735+178, such as the observing viewing angle, the plasma bulk Lorentz factor, and those deduced from these.

We present the results of simultaneous VLBA polarisation observations of the BL Lacertae object 0820+225 at 5, 8, and 15 GHz, together with earlier images at 5 GHz. This source has an unusually rich total intensity and polarisation structure compared to other objects with comparable redshifts. The magnetic field in the inner part of the complex and highly twisted VLBI jet is transverse, as is typical of BL Lacertae objects, but becomes roughly longitudinal further from the core, possibly due to shear. Although the integrated rotation measure of 0820+225 is modest, the rotation-measure distribution on parsec scales is non-uniform, and clearly shows regions where the rotation measure is substantially higher than the integrated value.

We present a new centimeter polarimetric VLBI image of the BL Lac object 0735+178. This source exhibits one of the most pronounced curvatures observed in jets of AGNs, with two sharp apparent bends of 90 degrees within the inner 2 milliarcseconds from the core. Through the analysis of the data gathered over the past decades we study whether this curvature is produced by precession of the jet with ejection of ballistic components, or a precession in such a way that components' velocity vector are always parallel to the jet axis. These possibilities are also studied by comparison with 3D hydrodynamic relativistic simulations of precessing jets.

We present results from VLBA+Effelsberg 1.7 GHz, VLBA 8.4 GHz and VLBA 15 GHz observations of the JVAS gravitational lens system B1030+074. The VLBA+Effelsberg 1.7 GHz data with 3 mas resolution reveal a detailed 60 mas jet structure in the strong A component of the lens system and maybe a first hint of the corresponding jet structure in the faint component B. As we go to the higher resolution 8.4 GHz and 15 GHz data we can see more of the inner structure of the A core, which is resolved, while the B core still remains unresolved.

We present two epoch polarimetric images of the BL Lac object 0735+178 obtained with the Very Long Baseline Array at 22 and 43 GHz. These images provide the highest resolution observations of this source to date, and show a twisted jet with two sharp apparent bends of 90 deg. within two milliarcseconds of the core. The magnetic field appears to smoothly follow one of the bends in the jet, suggesting that this structure may be the result of a precessing nozzle in the jet of 0735+178. Proper motion information is not conclusive as to whether the twisted geometry is the result of knots moving ballistically or following the direction of the jet axis.