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.

Context. Monitoring of BL Lacertae at 15 GHz with the Very Long Baseline Array (VLBA) has revealed a quasi-stationary radio feature in the innermost part of the jet, at 0.26 mas from the radio core. Stationary features are found in many blazars, but they have rarely been explored in detail. Aims. We aim to study the kinematics, dynamics, and brightness of the quasi-stationary feature of the jet in BL Lacertae based on VLBA monitoring with submilliarcsecond resolution (subparsec-scales) over 17 years. Methods. We analysed position uncertainties and flux leakage effects of the innermost quasi-stationary feature and developed statistical tools to distinguish the motions of the stationary feature and the radio core. We constructed a toy model to simulate the observed emission of the quasi-stationary component. Results. We find that trajectories of the quasi-stationary component are aligned along the jet axis, which can be interpreted as evidence of the displacements of the radio core. The intrinsic motions of the core and quasi-stationary component have a commensurate contribution to the apparent motion of the stationary component. During the jet-stable state, the core shift significantly influences the apparent displacements of the stationary component, which shows orbiting motion with reversals. The quasi-stationary component has low superluminal speeds on time scales of months. On time-scales of few years, the apparent mean speeds are subrelativistic, of about 0.15 the speed of light. We find that the brightness profile of the quasi-stationary component is asymmetric along and transverse to the jet axis, and this effect remains unchanged regardless of epoch. Conclusions. Accurate positional determination, a high cadence of observations, and a proper accounting for the core shift are crucial for the measurement of the trajectories and speeds of the quasi-stationary component. (abbrev.)

Observational studies of collimation in jets in active galactic nuclei (AGN) are a key to understanding their formation and acceleration processes. We have performed an automated search for jet shape transitions in a sample of 367 AGN using VLBA data at 15 GHz and 1.4 GHz. This search has found ten out of 29 nearby jets at redshifts z<0.07 with a transition from a parabolic to conical shape, while the full analyzed sample is dominated by distant AGN with a typical z about 1. The ten AGN are UGC00773, NGC1052, 3C111, 3C120, TXS0815-094, Mrk180, PKS1514+00, NGC6251, 3C371, and BL Lac. We conclude that the geometry transition may be a common effect in AGN jets. It can be observed only when sufficient linear resolution is obtained. Supplementing these results with previously reported shape breaks in the nearby AGN 1H0323+342 and M87, we estimate that the break occurs at 10^5-10^6 gravitational radii from the nucleus. We suggest that the jet shape transition happens when the bulk plasma kinetic energy flux becomes equal to the Poynting energy flux, while the ambient medium pressure is assumed to be governed by Bondi accretion. In general, the break point may not coincide with the Bondi radius. The observational data supports our model predictions on the jet acceleration and properties of the break point.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We present RadioAstron Space VLBI imaging observations of the BL Lac object S5 0716+71 made on January 3-4 2015 at a frequency of 22 GHz (wavelength $\lambda=1.3$ cm). The observations were made in the framework of the AGN Polarization Key Science Program. The source was detected on projected space-ground baselines up to 70 833 km (5.6 Earth diameters) for both, parallel hand and cross-hand interferometric visibilities. We have used these detections to obtain a full-polarimetric image of the blazar at an unprecedented angular resolution of 24 $\mu$as, the highest for this source to date. This enabled us to estimate the size of the radio core to be $<12\times5~\mu$as and to reveal a complex structure and a significant curvature of the blazar jet in the inner 100 $\mu$as, which is an indication that the jet viewing angle lies inside the opening angle of the jet conical outflow. Fairly highly (15%) linearly polarized emission is detected in a jet region of 19 $\mu$as in size, located 58 $\mu$as downstream from the core. The highest brightness temperature in the source frame is estimated to be $>2.2\times10^{13}$ K for the blazar core. This implies that the inverse Compton limit must be violated in the rest frame of the source, even for the largest Doppler factor $\delta\thicksim25$ reported for 0716+714.

Rotations of the electric vector position angle (EVPA) in blazars are often close to an integral multiple of 180$^\circ$. There are multiple examples of this in the literature, and our analysis here, of the optical polarization data from the RoboPol monitoring program, strengthens the evidence by showing that $n\pi$ rotations occur more frequently than expected by chance. We explain this with a model consisting of two polarized emission components: a "jet" that is constant in time, and a "burst" that is variable. The EVPA of the combination is $\rm EVPA_{jet}$ at both the beginning and the end of the burst, so the net rotation across the burst must be $n\pi$. Examples are analyzed on the Stokes plane, where the winding number for the Stokes vector of the combination gives the value of $n$. The main conclusion is that the EVPA rotation can be much larger than the physical rotation of the emission region around the axis of the jet, but this requires the EVPAs of the jet and the burst to be nearly orthogonal. A shock-in-jet calculation by Zhang et al. can provide a physical model for our toy model, and in addition automatically gives the needed orthogonality. The model is illustrated with data on OJ287 published by Myserlis et al., and we suggest that the large rapid EVPA rotation seen there might be a phase effect and not representative of a physical rotation.

Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.

Detecting and modelling the reprocessed hard X-ray emission component in the accretion flow, so-called reflection spectrum is a main tool to estimate black hole spins in a wide range of astrophysical black holes regardless of their mass or distance. In this work, we studied the X-ray spectra of the Seyfert I galaxy III Zw 2 using multi-epoch XMM-Newton, NuSTAR and Suzaku observations. The X-ray spectra exhibit a soft-excess below 1 keV and a prominent excess at the location of the broad Fe K{\alpha} line at 6.4 keV. To account for these spectral features, we have fitted the spectra with multiple models including an ionized partially covering absorber and an accretion disk reflection model. To fully resolve the reflection component, we analyzed jointly the XMM-Newton and NuSTAR observations taken in 2017 and archival XMM-Newton data from 2000. Assuming the reflection scenario, the resulting model fits support a rapidly spinning black hole (a > 0.98) in this radio-intermediate source. The X-ray spectra in 2000 and 2017 are remarkably similar with the only difference in the reflection fraction, possibly due to a change in the geometry of the accretion flow. However, the Suzaku observation is markedly different, and we suggest this could be an effect of a jet contribution in the X-ray band, which is supported by the elevated radio flux during this observation.

There are several methods to calculate the radiative and kinetic power of relativistic jets, but their results can differ by one or two orders of magnitude. Therefore, it is necessary to perform a calibration of the jet power, to understand the reasons for these differences (whether wrong hypotheses or intrinsic source variability), and if it is possible to converge to a reliable measurement of this physical quantity. We present preliminary results of a project aimed at calibrating the power of relativistic jets in active galactic nuclei (AGN) and X-ray binaries (XRB). We started by selecting all the AGN associations with known redshift in the Fourth Fermi LAT Gamma-Ray Catalog (4FGL). We then calculated the radiative and/or kinetic powers from available data or we extracted this information from literature. We compare the values obtained for overlapping samples and highlight early conclusions.

Detailed studies of relativistic jets in active galactic nuclei (AGN) require high-fidelity imaging at the highest possible resolution. This can be achieved using very long baseline interferometry (VLBI) at radio frequencies, combining worldwide (global) VLBI arrays of radio telescopes with a space-borne antenna on board a satellite. We present multiwavelength images made of the radio emission in the powerful quasar S5 0836+710, obtained using a global VLBI array and the antenna Spektr-R of the RadioAstron mission of the Russian Space Agency, with the goal of studying the internal structure and physics of the relativistic jet in this object. The RadioAstron observations at wavelengths of 18cm, 6cm, and 1.3cm are part of the Key Science Program for imaging radio emission in strong AGN. The internal structure of the jet is studied by analyzing transverse intensity profiles and modeling the structural patterns developing in the flow. The RadioAstron images reveal a wealth of structural detail in the jet of S5 0836+710 on angular scales ranging from 0.02mas to 200mas. Brightness temperatures in excess of $10^{13}$\,K are measured in the jet, requiring Doppler factors of $\ge 100$ for reconciling them with the inverse Compton limit. Several oscillatory patterns are identified in the ridge line of the jet and can be explained in terms of the Kelvin-Helmholtz (KH) instability. The oscillatory patterns are interpreted as the surface and body wavelengths of the helical mode of the KH instability. The interpretation provides estimates of the jet Mach number and of the ratio of the jet to the ambient density, which are found to be $M_\mathrm{j}\approx 12$ and $\eta\approx 0.33$. The ratio of the jet to the ambient density should be conservatively considered an upper limit because its estimate relies on approximations.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

The largest Key Science Program of the RadioAstron space VLBI mission is a survey of active galactic nuclei (AGN). The main goal of the survey is to measure and study the brightness of AGN cores in order to better understand the physics of their emission while taking interstellar scattering into consideration. In this paper we present detection statistics for observations on ground-space baselines of a complete sample of radio-strong AGN at the wavelengths of 18, 6, and 1.3 cm. Two-thirds of them are indeed detected by RadioAstron and are found to contain extremely compact, tens to hundreds of $\mu$as structures within their cores.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

Imaging relativistic jets in active galactic nuclei (AGN) at angular resolution significantly surpassing that of the ground-based VLBI at centimetre wavelengths is one of the key science objectives of the RadioAstron space-VLBI mission. There are three RadioAstron imaging key science programs that target both nearby radio galaxies and blazars, with one of the programs specifically focusing on polarimetry of the jets. The first images from these programs reach angular resolution of a few tens of microarcseconds and reveal unprecedented details about the jet collimation profile, magnetic field configuration, and Kelvin-Helmholtz instabilities along the flow in some of the most studied AGN (3C 84, BL Lac, 3C 273, S5 0836+710). Here we give an overview of the goals and strategy of these three ongoing programs, highlight their early results, and discuss the challenges of space-VLBI imaging.

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.

We studied the polarization behavior of the quasar 3C 273 over the 1 mm wavelength band at ALMA with a total bandwidth of 7.5GHz across 223 to 243 GHz at 0.8 arcsec resolution, corresponding to 2.1 kpc at the distance of 3C 273. With these observations we were able to probe the optically thin polarized emission close to the jet base, and constrain the magnetic field structure. We computed the Faraday rotation measure using simple linear fitting and Faraday rotation measure synthesis. In addition, we modeled the broadband behavior of the fractional Stokes Q and U parameters (qu-fitting). The systematic uncertainties in the polarization observations at ALMA were assessed through Monte Carlo simulations. We find the unresolved core of 3C 273 to be 1.8% linearly polarized. We detect a very high rotation measure (RM) of ~5.0 x10^5 rad/m^2 over the 1 mm band when assuming a single polarized component and an external RM screen. This results in a rotation of >40 deg of the intrinsic electric vector position angle, which is significantly higher than typically assumed for millimeter wavelengths. The polarization fraction increases as a function of wavelength, which according to our qu-fitting could be due to multiple polarized components of different Faraday depth within our beam, or internal Faraday rotation. With our limited wavelength coverage, we cannot distinguish between the cases, and additional multifrequency and high angular resolution observations are needed to determine the location and structure of the magnetic field of the Faraday active region. Comparing our RM estimate with values obtained at lower frequencies, the RM increases as a function of observing frequency, following a power law with an index of ~2.0 consistent with a sheath surrounding a conically expanding jet. We also detect ~0.2% circular polarization, although further observations are needed to confirm this result.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Ultra-light dark matter may consist of axion-like particles with masses below 10^(-19) eV. Two-photon interactions of these particles affect the polarization of radiation propagating through the dark matter. Coherent oscillations of the Bose condensate of the particles induce periodic changes in the plane of polarisation of emission passing through the condensate. We estimate this effect and analyze MOJAVE VLBA polarization observations of bright downstream features in the parsec-scale jets of active galaxies. Through the non-observation of periodic polarization changes, we are able to constrain the photon coupling of the ultra-light dark-matter axion-like particles at the level of <~ 10^(-12)/GeV for masses between ~5*10^(-23) eV and ~1.2*10^(-21) eV.

We present an analysis of the core linear polarization properties of 387 parsec-scale active galactic nuclei (AGN) jets. Using 15 GHz VLBA data, we revisit the conclusions of the first paper in this series with multi-epoch measurements and more detailed analysis of a larger AGN sample which spans a broader range of synchrotron peak frequencies. Each AGN has been observed for at least five epochs between 1996 and 2017. We find that BL Lac objects have core electric vector position angles (EVPAs) which tend towards alignment with the local jet direction; compared to flat spectrum radio quasars (FSRQs), their EVPAs are also less variable over time. The AGN cores which are most fractionally polarized and least variable in polarization have EVPAs that are closely aligned with the local jet direction; they also have low variability in EVPA. These results support the popular model of a standing transverse shock at the base of the jet which collimates the jet magnetic field perpendicular to the jet direction, increasing the fractional polarization and leading to greater polarization stability over time. High-synchrotron-peaked (HSP) BL Lac objects form a low luminosity, low fractional polarization population. The five narrow-line Seyfert 1 galaxies in our sample have low fractional polarization and large EVPA-jet misalignments. Although AGN detected at gamma-rays are thought to be more Doppler boosted than non-detected AGN, we find no significant differences in fractional polarization based on detection by Fermi-LAT; the gamma-loud AGN are, however, more variable in core EVPAs.

We have obtained a smooth time series for the Electric Vector Position Angle (EVPA) of the blazar OJ 287 at centimeter wavelengths, by making $\pm n\pi$ adjustments to archival values from 1974 to 2016. The data display rotation reversals in which the EVPA rotates counter-clockwise (CCW) for 180 deg and then rotates clockwise (CW) by a similar amount. The time scale of the rotations is a few weeks to a year, and the scale for a double rotation, including the reversal, is one to three years. We have seen four of these events in 40 years. A model consisting of two successive outbursts in polarized flux density, with EVPAs counter-rotating, superposed on a steady polarized jet, can explain many of the details of the observations. Polarization images support this interpretation. The model can also help to explain similar events seen at optical wavelengths. The outbursts needed for the model can be generated by the super-magnetosonic jet model of Nakamura et al. (2010) and Nakamura and Meier (2014), which requires a strong helical magnetic field. This model produces forward and reverse pairs of fast and slow MHD waves, and the plasma inside the two fast/slow pairs rotates around the jet axis, but in opposite directions.

Understanding the launching, acceleration, and collimation of jets powered by active galactic nuclei remains an outstanding problem in relativistic astrophysics. This is partly because observational tests of jet formation models suffer from the limited angular resolution of ground-based very long baseline interferometry that has thus far been able to probe the transverse jet structure in the acceleration and collimation zone of only two sources. Here we report radio interferometric observations of 3C 84 (NGC 1275), the central galaxy of the Perseus cluster, made with an array including the orbiting radio telescope of the RadioAstron mission. The obtained image transversely resolves the edge-brightened jet in 3C 84 only 30 microarcseconds from the core, which is ten times closer to the central engine than what has been possible in previous ground-based observations, and it allows us to measure the jet collimation profile from ~ 100 to ~10000 gravitational radii from the black hole. The previously found, almost cylindrical jet profile on scales larger than a few thousand r_g is now seen to continue at least down to a few hundred r_g from the black hole and we find a broad jet with a transverse radius larger than about 250 r_g at only 350 r_g from the core. If the bright outer jet layer is launched by the black hole ergosphere, it has to rapidly expand laterally on scales smaller than 100 r_g. If this is not the case, then this jet sheath is likely launched from the accretion disk.

Combined with very-long-baseline interferometry measurements, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL~Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL~Lacertae observed by VERITAS, with a rise time of $\sim$2.3~hr and a decay time of $\sim$36~min. The peak flux above 200 GeV is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photon} \;\text{m}^{-2}\; \text{s}^{-1}$ measured with a 4-minute-binned light curve, corresponding to $\sim$180\% of the flux which is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in VLBA observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models which invoke relativistic plasma passing stationary shocks.

We used 15 GHz multi-epoch Very Long Baseline Array (VLBA) polarization sensitive observations of 484 sources within a time interval 1996--2016 from the MOJAVE program, and also from the NRAO data archive. We have analyzed the linear polarization characteristics of the compact core features and regions downstream, and their changes along and across the parsec-scale active galactic nuclei (AGN) jets. We detected a significant increase of fractional polarization with distance from the radio core along the jet as well as towards the jet edges. Compared to quasars, BL Lacs have a higher degree of polarization and exhibit more stable electric vector position angles (EVPAs) in their core features and a better alignment of the EVPAs with the local jet direction. The latter is accompanied by a higher degree of linear polarization, suggesting that compact bright jet features might be strong transverse shocks, which enhance magnetic field regularity by compression.

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.

Observations of fast TeV $\gamma$-ray flares from blazars reveal the extreme compactness of emitting regions in blazar jets. Combined with very-long-baseline radio interferometry measurements, they probe the structure and emission mechanism of the jet. We report on a fast TeV $\gamma$-ray flare from BL Lacertae observed by VERITAS, with a rise time of about 2.3 hours and a decay time of about 36 minutes. The peak flux at $>$200 GeV measured with the 4-minute binned light curve is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photons} \;\text{m}^{-2}\, \text{s}^{-1}$, or $\sim$180% the Crab Nebula flux. Variability in GeV $\gamma$-ray, X-ray, and optical flux, as well as in optical and radio polarization was observed around the time of the TeV $\gamma$-ray flare. A possible superluminal knot was identified in the VLBA observations at 43 GHz. The flare constrains the size of the emitting region, and is consistent with several theoretical models with stationary shocks.

RadioAstron is a 10 m orbiting radio telescope mounted on the Spektr-R satellite, launched in 2011, performing Space Very Long Baseline Interferometry (SVLBI) observations supported by a global ground array of radio telescopes. With an apogee of about 350 000 km, it is offering for the first time the possibility to perform {\mu}as-resolution imaging in the cm-band. We present observations at 22 GHz of 3C 273, performed in 2014, designed to reach a maximum baseline of approximately nine Earth diameters. Reaching an angular resolution of 0.3 mas, we study a particularly low-activity state of the source, and estimate the nuclear region brightness temperature, comparing with the extreme one detected one year before during the RadioAstron early science period. We also make use of the VLBA-BU-BLAZAR survey data, at 43 GHz, to study the kinematics of the jet in a 1.5-year time window. We find that the nuclear brightness temperature is two orders of magnitude lower than the exceptionally high value detected in 2013 with RadioAstron at the same frequency (1.4x10^13 K, source-frame), and even one order of magnitude lower than the equipartition value. The kinematics analysis at 43 GHz shows that a new component was ejected 2 months after the 2013 epoch, visible also in our 22 GHz map presented here. Consequently this was located upstream of the core during the brightness temperature peak. These observations confirm that the previously detected extreme brightness temperature in 3C 273, exceeding the inverse Compton limit, is a short-lived phenomenon caused by a temporary departure from equipartition. Thus, the availability of interferometric baselines capable of providing {\mu}as angular resolution does not systematically imply measured brightness temperatures over the known physical limits for astrophysical sources.

We present 15 GHz stacked VLBA images of 373 jets associated with active galactic nuclei (AGN) having at least five observing epochs within a 20 yr time interval 1994-2015 from the MOJAVE programme and/or its precursor, the 2 cm VLBA Survey. These data are supplemented by 1.4 GHz single-epoch VLBA observations of 135 MOJAVE AGNs to probe larger scale jet structures. The typical jet geometry is found to be close to conical on scales from hundreds to thousands of parsecs, while a number of galaxies show quasi-parabolic streamlines on smaller scales. A true jet geometry in a considerable fraction of AGNs appears only after stacking epochs over several years. The jets with significant radial accelerated motion undergo more active collimation. We have analysed total intensity jet profiles transverse to the local jet ridgeline and derived both apparent and intrinsic opening angles of the flows, with medians of $21.5\deg$ and $1.3\deg$, respectively. The Fermi LAT-detected gamma-ray AGNs in our sample have, on average, wider apparent and narrower intrinsic opening angle, and smaller viewing angle than non LAT-detected AGNs. We have established a highly significant correlation between the apparent opening angle and gamma-ray luminosity, driven by Doppler beaming and projection effects.

We present a comprehensive 5-43 GHz VLBA study of the blazar 3C 273 initiated after an onset of a strong $\gamma$-ray flare in this source. We have analyzed the kinematics of new-born components, light curves, and position of the apparent core to pinpoint the location of the $\gamma$-ray emission. Estimated location of the $\gamma$-ray emission zone is close to the jet apex, 2 pc to 7 pc upstream from the observed 7 mm core. This is supported by ejection of a new component. The apparent core position was found to be inversely proportional to frequency. The brightness temperature in the 7 mm core reached values up to at least $10^{13}$ K during the flare. This supports the dominance of particle energy density over that of magnetic field in the 7 mm core. Particle density increased during the radio flare at the apparent jet base, affecting synchrotron opacity. This manifested itself as an apparent core shuttle along the jet during the 7 mm flare. It is also shown that a region where optical depth decreases from $\tau\sim1$ to $\tau<<1$ spans over several parsecs along the jet. The jet bulk flow speed estimated at the level of 12c on the basis of time lags between 7 mm light curves of stationary jet features is 1.5 times higher than that derived from VLBI apparent kinematics analysis.

Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.

1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $\gamma$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Mets\"ahovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission.

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)

High-resolution millimeter wave very-long-baseline interferometry (mm-VLBI) is an ideal tool for probing the structure at the base of extragalactic jets in detail. The TeV blazar Mrk 501 is one of the best targets among BL Lac objects for studying the nature of off-axis jet structures because it shows different jet position angles at different scales. The aim of this study is to investigate the properties of the off-axis jet structure through high-resolution mm-VLBI images at the jet base and physical parameters such as kinematics, flux densities, and spectral indices. We performed Very Long Baseline Array (VLBA) observations over six epochs from 2012 February to 2013 February at 43 GHz. Quasi-simultaneous Global Millimeter VLBI Array (GMVA) observations at 86 GHz were performed in May 2012. We discover a new jet component at the northeast direction from the core in all the images at 43 and 86 GHz. The new component shows the off-axis location from the persistent jet extending to the southeast. The 43 GHz images reveal that the scattering of the positions of the NE component is within ~0.2 mas. The 86 GHz data reveals a jet component located 0.75 mas southeast of the radio core. We also discuss the spectral indices between 43 and 86 GHz, where the northeast component has steeper spectral index and the southeast component has comparable or flatter index than the radio core does.

Following the discovery of high-energy (HE; $E>10\,{\rm MeV}$) and very-high-energy (VHE; $E>100\,{\rm GeV}$) $\gamma$-ray emission from the low-frequency-peaked BL~Lac (LBL) object AP Librae, its electromagnetic spectrum is studied over 60 octaves in energy. Contemporaneous data in radio, optical and UV together with the $\gamma$-ray data are used to construct the most precise spectral energy distribution of this source. The data have been found to be modeled with difficulties with single zone homogeneous leptonic synchrotron self-Compton (SSC) radiative scenarios due to the unprecedented width of the high-energy component when compared to the lower-energy component. The two other LBL objects also detected at VHE appear to have similar modeling difficulties. Nevertheless, VHE $\gamma$ rays produced in the extended jet could account for the VHE flux observed by H.E.S.S.

Opacity-driven shifts of the apparent VLBI core position with frequency (the "core shift" effect) probe physical conditions in the innermost parts of jets in active galactic nuclei. We present the first detailed investigation of this effect in the brightest gamma-ray blazar 3C454.3 using direct measurements from simultaneous 4.6-43 GHz VLBA observations, and a time lag analysis of 4.8-37 GHz lightcurves from the UMRAO, CrAO, and Metsahovi observations in 2007-2009. The results support the standard Konigl model of jet physics in the VLBI core region. The distance of the core from the jet origin r_c(nu), the core size W(nu), and the lightcurve time lag DT(nu) all depend on the observing frequency nu as r_c(nu)~W(nu)~ DT(nu)~nu^-1/k. The obtained range of k=0.6-0.8 is consistent with the synchrotron self-absorption being the dominating opacity mechanism in the jet. The similar frequency dependence of r_c(nu) and W(nu) suggests that the external pressure gradient does not dictate the jet geometry in the cm-band core region. Assuming equipartition, the magnetic field strength scales with distance r as B = 0.4(r/1pc)^-0.8 G. The total kinetic power of electron/positron jet is about 10^44 ergs/s.

We present results from 43 GHz (VLBA, six epochs from 2012.2 to 2013.2) and 86 GHz (GMVA, one epoch in 2012.4) observations toward the basis of the jet in the TeV Blazar Mrk 501. The 43-GHz data analysis reveals a new feature located northeast of the radio core, with a flux density of several tens of mJy, perpendicularly to the jet axis. The 86-GHz image shows the jet feature located 0.75 mas southeast of the radio core, which is consistent with the previous result. The location of Gaussian model for 0.75 mas feature does not coincide with those for the jets in the 43-GHz image, however, a distribution of emission is found. We also discuss the spectral indices of the core, the northeast feature, and the jet feature between 43 GHz and 86 GHz, which show flat-to-steep, steep, and flat-to-invert, respectively.

We present the analysis of the radio jet evolution of the radio galaxy 3C 120 during a period of prolonged gamma-ray activity detected by the Fermi satellite between December 2012 and October 2014. We find a clear connection between the gamma-ray and radio emission, such that every period of gamma-ray activity is accompanied by the flaring of the mm-VLBI core and subsequent ejection of a new superluminal component. However, not all ejections of components are associated with gamma-ray events detectable by Fermi. Clear gamma-ray detections are obtained only when components are moving in a direction closer to our line of sight.This suggests that the observed gamma-ray emission depends not only on the interaction of moving components with the mm-VLBI core, but also on their orientation with respect to the observer. Timing of the gamma-ray detections and ejection of superluminal components locate the gamma-ray production to within almost 0.13 pc from the mm-VLBI core, which was previously estimated to lie about 0.24 pc from the central black hole. This corresponds to about twice the estimated extension of the broad line region, limiting the external photon field and therefore suggesting synchrotron self Compton as the most probable mechanism for the production of the gamma-ray emission. Alternatively, the interaction of components with the jet sheath can provide the necessary photon field to produced the observed gamma-rays by Compton scattering.

The radio galaxy 3C84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGN) and one of the best laboratories to study the radio properties of subparsec scale jets. We discuss here the past and present activity of the nuclear region within the central 1pc and the properties of subparsec-sized components C1, C2 and C3. We compare these results with the high resolution space-VLBI image at 5GHz obtained with the RadioAstron satellite and we shortly discuss the possible correlation of radio emission with the gamma-ray emission.

We used 15 GHz VLBA observations of 366 sources having at least 5 epochs within a time interval 1995-2013 from the MOJAVE program and/or its predecessor, the 2 cm VLBA Survey. For each source we produced a corresponding stacked image averaging all available epochs for a better reconstruction of the cross section of the flow. We have analyzed jet profiles transverse to the local jet ridge line and derived both apparent and intrinsic opening angles of the parsec-scale outflows. The sources detected by the Fermi Large Area Telescope (LAT) during the first 24 months of operation show wider apparent jet opening angle and smaller viewing angles on a very high level of significance supporting our early findings. Analyzing transverse shapes of the outflows we found that most sources have conical jet geometry at parsec scales, though there are also sources that exhibit active jet collimation.

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.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of radio-loud Active Galactic Nuclei. For a summary, we refer to the paper.

Relativistic jets in active galactic nuclei represent one of the most powerful phenomena in the Universe. They form in the surroundings of the supermassive black holes as a by-product of accretion onto the central black hole in active galaxies. The flow in the jets propagates at velocities close to the speed of light. The distance between the first part of the jet that is visible in radio images (core) and the black hole is still a matter of debate. Only very-long-baseline interferometry observations resolve the innermost compact regions of the radio jet. Those can access the jet base, and combining data at different wavelenghts, address the physical parameters of the outflow from its emission. We have performed an accurate analysis of the frequency-dependent shift of the VLBI core location for a multi-wavelength set of images of the blazar CTA 102 including data from 6 cm down to 3 mm. The measure of the position of the central black hole, with mass $\sim 10^{8.93}\,M_\odot$, in the blazar CTA 102 reveals a distance of $\sim 8\times10^4$ gravitational radii to the 86 GHz core, in agreement with similar measures obtained for other blazars and distant radio galaxies, and in contrast with recent results for the case of nearby radio galaxies, which show distances between the black hole and the radio core that can be two orders of magnitude smaller.

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.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

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 report that active galactic nucleus (AGN) jets are causally connected on parsec scales, based on 15 GHz Very Long Baseline Array (VLBA) data from a sample of 133 AGN jets. This result is achieved through a new method for measuring the product of the jet Lorentz factor and the intrinsic opening angle Gamma*theta_j from measured apparent opening angles in flux density limited samples of AGN jets. The Gamma*theta_j parameter is important for jet physics because it is related to the jet-frame sidewise expansion speed and causal connection between the jet edges and its symmetry axis. Most importantly, the standard model of jet production requires that the jet be causally connected with its symmetry axis, implying that Gamma*theta_j < 1. When we apply our method to the MOJAVE flux density limited sample of radio loud objects, we find Gamma*theta_j = 0.2, implying that AGN jets are causally connected. We also find evidence that AGN jets viewed very close to the line of sight effectively have smaller intrinsic opening angles compared with jets viewed more off-axis, which is consistent with Doppler beaming and a fast inner spine/slow outer sheath velocity field. Notably, gamma-ray burst (GRB) jets have a typical Gamma*theta_j that is two orders of magnitude higher, suggesting that different physical mechanisms are at work in GRB jets compared to AGN jets. A useful application of our result is that a jet's beaming parameters can be derived. Assuming Gamma*theta_j is approximately constant in the AGN jet population, an individual jet's Doppler factor and Lorentz factor (and therefore also its viewing angle) can be determined using two observable quantities: apparent jet opening angle and the apparent speed of jet components.

Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of both random and deterministic cases. Our analysis indicates that the EVPA variation is likely dominated by a random process in the low brightness state of the jet and by a deterministic process in the flaring state.

In September 2012, the high-synchrotron-peaked (HSP) blazar Markarian 421 underwent a rapid wideband radio flare, reaching nearly twice the brightest level observed in the centimeter band in over three decades of monitoring. In response to this event we carried out a five epoch centimeter- to millimeter-band multifrequency Very Long Baseline Array (VLBA) campaign to investigate the aftermath of this emission event. Rapid radio variations are unprecedented in this object and are surprising in an HSP BL Lac object. In this flare, the 15 GHz flux density increased with an exponential doubling time of about 9 days, then faded to its prior level at a similar rate. This is comparable with the fastest large-amplitude centimeter-band radio variability observed in any blazar. Similar flux density increases were detected up to millimeter bands. This radio flare followed about two months after a similarly unprecedented GeV gamma-ray flare (reaching a daily E>100 MeV flux of (1.2 +/- 0.7)x10^(-6) ph cm^(-2) s^(-1)) reported by the Fermi Large Area Telescope (LAT) collaboration, with a simultaneous tentative TeV detection by ARGO-YBJ. A cross-correlation analysis of long-term 15 GHz and LAT gamma-ray light curves finds a statistically significant correlation with the radio lagging ~40 days behind, suggesting that the gamma-ray emission originates upstream of the radio emission. Preliminary results from our VLBA observations show brightening in the unresolved core region and no evidence for apparent superluminal motions or substantial flux variations downstream.

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)

The nearby active galaxy IC 310 (z=0.019), located in the Perseus cluster of galaxies is a bright and variable multi-wavelength emitter from the radio regime up to very high gamma-ray energies above 100 GeV. Very recently, a blazar-like compact radio jet has been found by parsec-scale VLBI imaging. Along with the unusually flat gamma-ray spectrum and variable high-energy emission, this suggests that IC 310 is the closest known blazar and therefore a key object for AGN research. As part of an intense observing program at TeV energies with the MAGIC telescopes, an exceptionally bright flare of IC 310 was detected in November 2012 reaching a flux level of up to >0.5 Crab units above 300 GeV. We have organized a multi-wavelength follow-up program, including the VLBA, Effelsberg 100 m, KVA, Swift, INTEGRAL, Fermi/LAT, and the MAGIC telescopes. We present preliminary results from the multi-wavelength follow-up program with the focus on the response of the jet to this exceptional gamma-ray flare.

The temporal and spatial spectral evolution of the jets of AGN can be studied with multi-frequency, multi-epoch VLBI observations. The combination of both, morphological and spectral parameters can be used to derive source intrinsic physical properties such as the magnetic field and the non-thermal particle density. In the first two papers of this series, we analyzed the single-dish light curves and the VLBI kinematics of the blazar CTA 102 and suggested a shock-shock interaction between a traveling and a standing shock wave as a possible scenario to explain the observed evolution of the component associated to the 2006 flare. In this paper we investigate the core-shift and spectral evolution to test our hypothesis of a shock-shock interaction. We used 8 multi-frequency VLBA observations to analyze the temporal and spatial evolution of the spectral parameters during the flare. We observed CTA 102 between May 2005 and April 2007 using the VLBA at six different frequencies spanning from 2 - 86 GHz. After the calibrated VLBA images were corrected for opacity, we performed a detailed spectral analysis. From the derived values we estimated the magnetic field and the density of the relativistic particles. The detailed analysis of the opacity shift reveals that the position of the jet core is proportional to nu^-1 with some temporal variations. The value suggests possible equipartition between magnetic field energy and particle kinetic energy densities at the most compact regions. From the variation of the physical parameters we deduced that the 2006 flare in CTA 102 is connected to the ejection of a new traveling feature (t=2005.9) and the interaction between this shock wave and a stationary structure around 0.1 mas from the core. The source kinematics together with the spectral and structural variations can be described by helical motions in an over-pressured jet.

We report on the first detection of the theoretically-predicted rare phenomenon of multiple parsec-scale imaging of an active galactic nucleus induced by refractive effects due to localized foreground electron density enhancements, e.g., in an AU-scale plasma lens(es) in the ionized component of the Galactic interstellar medium. We detected multiple imaging in the low galactic latitude (b=-2 deg) quasar 2023+335 from the 15.4 GHz MOJAVE observations when the source was undergoing an ESE. While the parsec-scale jet of the source normally extends along PA -20 deg, in the 28 May 2009 and 23 July 2009 images a highly significant multi-component pattern of secondary images is stretched out nearly along the constant galactic latitude line with a local PA 40 deg, indicating that the direction of relative motion of the plasma lens is close to orbital. Weaker but still detectable imaging patterns at similar position angles are sporadically manifest at several other epochs. Modeling the ESE that occurred in early 2009 and lasted ~0.14 yr, we determined that the foreground screen has a double-lens structure, with proper motion (~6.8 mas/yr), and angular size (~0.27 mas). We also found that the angular separation between the two brightest sub-images roughly follows a wavelength-squared dependence expected from plasma scattering. Furthermore, by analyzing archival non-simultaneous VLBA observations covering a wide frequency range from 1.4 to 86 GHz, we found that the scattered angular size of the VLBI core follows a $\nu^{-1.89}$ dependence, implying the presence of a turbulent, refractive dominated scattering screen that has a confined structure or is truncated transverse to the line of sight toward 2023+335.

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.

The 2006 radio flare in the blazar CTA102 is the largest ever reported for this source. The analysis of the single--dish light curves revealed a possible shock-shock interaction as a scenario for the 2006 outburst. In order to confirm this hypothesis we analyzed eight multi-frequency (2 GHz - 86 GHz) VLBI observations covering the 2006 flare. We detected a radio component ejected around 2005.9 and a stationary feature around 0.1 mas away from the core. We also performed a core-shift and spectral analysis on the fully calibrated VLBI maps. The frequency dependent position of the core slightly deviates from the typical $\nu^{-1}$-dependence during the flaring state. This behaviour may reflect a mismatch between magnetic energy density and energy density of the relativistic particles. From the core-shift corrected maps we obtained the spectral parameters and computed the temporal and spatial evolution of the jet intrinsic parameters such as the magnetic field and the particle density. The overall picture of the source could be best interpreted as an over-pressured jet perturbed by traveling shock waves.

We present a new method to measure Gamma*theta_j in flux-limited samples of active galactic nuclei (AGN) jets, where Gamma is the bulk Lorentz factor and theta_j is the jet's half-opening angle. The Gamma*theta_j parameter is physically important for models of jet launching, and also determines the effectiveness of jet instabilities and magnetic reconnection. We measure Gamma*theta_j by analyzing the observed distribution of apparent opening angles in very long baseline interferometry (VLBI) flux-limited samples of jets, given some prior knowledge of the active galactic nuclei (AGN) radio luminosity function. We then apply this method to the MOJAVE flux-limited sample of radio loud objects and find Gamma*theta_j = 0.1 +- 0.03, which implies that AGN jets are subject to a variety of physical processes that require causal connection.

We performed multifrequency multiepoch Very Long Baseline Array (VLBA) observations of the blazar CTA 102 during its 2006 radio flare, the strongest ever reported for this source. These observations provide an excellent opportunity to investigate the evolution of the physical properties of blazars, especially during these flaring events. We want to study the kinematic changes in the source during the strong radio outburst in April 2006 and test the assumption of a shock-shock interaction. This assumption is based on the analysis and modeling of the single-dish observations of CTA\,102 (Paper I). In this paper we study the kinematics of CTA 102 at several frequencies using VLBI observations. From the modeled jet features we derived estimates for the evolution of the physical parameters, such as the particle density and the magnetic field. Furthermore, we combined our observations during the 2006 flare with long-term VLBA monitoring of the source at 15 GHz and 43 GHz. We cross-identified seven features throughout our entire multifrequency observations and find evidence of two possible recollimation shocks around 0.1 mas (deprojected 18 pc at a viewing angle of 2.6 degrees) and 6.0 mas (deprojected 1 kpc) from the core. The 43 GHz observations reveal a feature ejected at epoch $t_\mathrm{ej}=2005.9\pm0.2$, which could be connected to the 2006 April radio flare. Furthermore, this feature might be associated with the traveling component involved in the possible shock-shock interaction, which gives rise to the observed double peak structure in the single-dish light curves reported in Paper I.

We have investigated a frequency-dependent shift in the absolute position of the optically thick apparent origin of parsec-scale jets ("core shift" effect) to probe physical conditions in ultra-compact relativistic outflows in AGN. We used multi-frequency Very Long Baseline Array (VLBA) observations of 191 sources carried out in 12 epochs in 2006 within the MOJAVE program. The observations were performed at 8.1, 8.4, 12.1, and 15.4 GHz. We implemented a method of determining the core shift vector based on (i) image registration by two-dimensional normalized cross-correlation and (ii) model-fitting the source brightness distribution to take into account a non-zero core component offset from the phase center. The 15.4-8.1, 15.4-8.4, and 15.4-12.1 GHz core shift vectors are derived for 163 sources, and have median values of 0.128, 0.125, and 0.088 mas, respectively, compared to the typical measured errors of 0.050, 0.051, 0.035 mas. The effect occurs predominantly along the jet direction, with departures smaller than 45 deg from the median jet position angle in over 80% of the cases. Despite the moderate ratio of the observed frequencies (<2), core shifts significantly different from zero (>2sigma) are detected for about 55% of the sources. These shifts are even better aligned with the jet direction, deviating from the latter by less than 30 deg in over 90% of the cases. There is an indication that the core shift decreases with increasing redshift. Magnetic fields in the jet at a distance of 1 parsec from the central black hole, calculated from the obtained core shifts, are found to be systematically stronger in quasars (median B1~0.9 G) than those in BL Lacs (median B1~0.4 G). We also constrained the absolute distance of the core from the apex of the jet at 15 GHz as well as the magnetic field strength in the 15 GHz core region.

We have investigated the differences in apparent opening angles between the parsec-scale jets of the active galactic nuclei (AGN) detected by the Fermi Large Area Telescope (LAT) during its first 24 months of operations and those of non-LAT-detected AGN. We used 15.4 GHz VLBA observations of 215 sources from the 2 cm VLBA MOJAVE program. The apparent opening angles were determined by analyzing transverse jet profiles from the data in the image plane by using stacked images constructed from all available MOJAVE epochs for a given source. We confirm our earlier result based on the first three months of scientific operations of the LAT. The apparent opening angles of \gamma-ray bright AGN are preferentially larger than those of \gamma-ray weak sources, suggesting smaller viewing angles to the \gamma-ray bright AGN. Intrinsic opening angles for BL Lacs are wider than those in quasars.

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.

The analysis of the radio light curves of the blazar CTA102 during its 2006 flare revealed a possible interaction between a standing shock wave and a traveling one. In order to better understand this highly non-linear process, we used a relativistic hydrodynamic code to simulate the high energy interaction and its related emission. The calculated synchrotron emission from these simulations showed an increase in turnover flux density, $S_{m}$, and turnover frequency, $\nu_{m}$, during the interaction and decrease to its initial values after the passage of the traveling shock wave.

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.

Context: The blazar CTA 102 (z=1.037) underwent a historical radio outburst in April 2006. This event offered a unique chance to study the physical properties of the jet. Aims: We used multifrequency radio and mm observations to analyze the evolution of the spectral parameters during the flare as a test of the shock-in-jet model under these extreme conditions. Methods: For the analysis of the flare we took into account that the flaring spectrum is superimposed on a quiescent spectrum. We reconstructed the latter from archival data and fitted a synchrotron self-absorbed distribution of emission. The uncertainties of the derived spectral parameters were calculated using Monte Carlo simulations. The spectral evolution is modeled by the shock-in-jet model, and the derived results are discussed in the context of a geometrical model (varying viewing angle) and shock-shock interaction. Results: The evolution of the flare in the turnover frequency-turnover flux density plane shows a double peak structure. The nature of this evolution is dicussed in the frame of shock-in-jet models. We discard the generation of the double peak structure in the turnover frequency-turnover flux density plane purely based on geometrical changes (variation of the Doppler factor). The detailed modeling of the spectral evolution favors a shock-shock interaction as a possible physical mechanism behind the deviations from the standard shock-in-jet model.

The radio light curve and spectral evolution of the blazar CTA 102 during its 2006 outburst can be rather well explained by the standard shock-in-jet model. The results of a pixel-to-pixel spectral analysis of multi-frequency VLBI images, together with kinematics derived from the MOJAVE survey lead to the picture of an over-pressured jet with respect to the ambient medium. The interaction of a traveling shock wave with a standing one is a possible scenario which could explain the observed spectral behaviour

We have found evidence for interaction between a standing and a traveling shock in the jet of the blazar CTA 102. Our result is based in the study of the spectral evolution of the turnover frequency-turnover flux density plane. The radio/mm light curves were taken during a major radio outburst in April 2006.

The blazar CTA 102 underwent a major radio flare in April 2006. We used several 15 GHz VLBI observations from the MOJAVE program to investigate the influence of this extreme event on jet kinematics. The result of modeling and analysis lead to the suggestion of an interaction between traveling and standing shocks 0.2 mas away from the VLBI core.

Single-zone synchrotron self-Compton and external Compton models are widely used to explain broad-band Spectral Energy Distributions (SEDs) of blazars from infrared to gamma-rays. These models bear obvious similarities to the homogeneous synchrotron cloud model which is often applied to explain radio emission from individual components of parsec-scale radio jets. The parsec-scale core, typically the brightest and most compact feature of blazar radio jet, could be the source of high-energy emission. We report on ongoing work to test this hypothesis by deriving the physical properties of parsec-scale radio emitting regions of twenty bright Fermi blazars using dedicated 5-43 GHz VLBA observations and comparing these parameters to results of SED modeling.

We have compared the parsec-scale jet linear polarization properties of the Fermi LAT-detected and non-detected sources in the complete flux-density-limited (MOJAVE-1) sample of highly beamed AGN. Of the 123 MOJAVE sources, 30 were detected by the LAT during its first three months of operation. We find that during the era since the launch of Fermi, the unresolved core components of the LAT-detected jets have significantly higher median fractional polarization at 15 GHz. This complements our previous findings that these LAT sources have higher apparent jet speeds, brightness temperatures and Doppler factors, and are preferentially found in higher activity states.

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.

The spatially resolved broad-band spectroscopy with Very Long Baseline Interferometry (VLBI) is one of the few methods that can probe the physical conditions inside blazar jets. We report on measurements of the magnetic field strength in parsec-scale radio structures of selected bright Fermi blazars, based on fitting the synchrotron spectrum to VLBA images made at seven frequencies in a 4.6 -- 43.2 GHz range. Upper limits of B <= 10^-2 -- 10^2 G (observer's frame) could be placed on the magnetic field strength in 13 sources. Hard radio spectra (-0.5 <= a <= +0.1, S_nu ~ nu^a) observed above the synchrotron peak may either be an indication of a hard energy spectrum of the emitting electron population or result from a significant inhomogeneity of the emitting region.

We have investigated the differences in apparent opening angles between the parsec-scale jets of the active galactic nuclei (AGN) detected by the Fermi Large Area Telescope (LAT) during its first three months of operations and those of non-LAT-detected AGN. We used 15.4 GHz VLBA observations of sources from the 2 cm VLBA MOJAVE program, a subset of which comprise the statistically complete flux density limited MOJAVE sample. We determined the apparent opening angles by analyzing transverse jet profiles from the data in the image plane and by applying a model fitting technique to the data in the (u,v) plane. Both methods provided comparable opening angle estimates. The apparent opening angles of gamma-ray bright blazars are preferentially larger than those of gamma-ray weak sources. At the same time, we have found the two groups to have similar intrinsic opening angle distributions, based on a smaller subset of sources. This suggests that the jets in gamma-ray bright AGN are oriented at preferentially smaller angles to the line of sight resulting in a stronger relativistic beaming. The intrinsic jet opening angle and bulk flow Lorentz factor are found to be inversely proportional, as predicted by standard models of compact relativistic jets. If a gas dynamical jet acceleration model is assumed, the ratio of the initial pressure of the plasma in the core region P_0 to the external pressure P_ext lies within the range 1.1 to 34.6, with a best fit estimate of P_0/P_ext=2.

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)

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.

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 observed the blazar \object{0716+714} with the VLBA during its active state in 2003-2004. In this paper we discuss multi-frequency analysis of the inner jet (first 1 mas) kinematics. The unprecedentedly dense time sampling allows us to trace jet components without misidentification and to calculate the component speeds with good accuracy. In the smooth superluminal jet we were able to identify and track three components over time moving outwards with relatively high apparent superluminal speeds (8.5-19.4 $c$), which contradicts the hypothesis of a stationary oscillating jet in this source. Component ejections occur at a relatively high rate (once in two months), and they are accompanied by mm-continuum outbursts. Superluminal jet components move along wiggling trajectories, which is an indication of actual helical motion. Fast proper motion and rapid decay of the components suggest that this source should be observed with the VLBI at a rate of at least once in one or two months in order to trace superluminal jet components without confusion.

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.

Multi-frequency VLBI observations allow studies of the continuum spectrum in the different parts of the parsec scale jets of AGN, providing information on the physical properties of the plasma and magnetic fields in them. Since VLBI networks cannot be scaled, the range of spatial frequencies observed differs significantly between the different observing frequencies, which makes it difficult to obtain a broadband spectrum of the individual emission features in the jet. In this paper we discuss a model-fitting based spectral extraction method, which can significantly relieve this problem. The method uses a priori knowledge of the source structure, measured at high frequencies, to allow at lower frequencies the derivation of the sizes and flux densities of even those emission features that have mutual separations significantly less than the Rayleigh limit at the given frequency. We have successfully used this method in the analysis of 5-86 GHz VLBA data of 3C273. The spectra and sizes of several individual jet features were measured, thus allowing derivation of the magnetic flux density and the energy density of the relativistic electrons in the different parts of the jet. We discuss the results, which include e.g. a detection of a strong gradient in the magnetic field across the jet of 3C273.

We present here the discovery of rapid, large amplitude intraday variability in the compact flat-spectrum radio quasar 1156+295. The detection of 40% flux density variations at 15 GHz on a timescale of only 2.7 hours was serendipitously made when the source was observed with the Very Long Baseline Array as a part of the MOJAVE survey programme on February 5, 2007. Intraday variability on timescales of a few hours or less is rare, and there exist very few sources that show large-amplitude variations on a timescale as short as what is now observed for 1156+295. The shape of the visibility function of the source changes very little during the observation, although the correlated flux density changes by 40%. This suggests that the variability occurs in a single dominant compact component. The observed variability characteristics are consistent with interstellar scintillation in nearby, highly turbulent medium. The rms amplitude of modulation at 15 GHz is unusually large and it implies a rather high scattering measure along the line-of-sight towards 1156+295.

We report a serendipitous detection of rapid, large amplitude flux density variations in the highly core-dominated, flat-spectrum radio quasar 1156+295 during an observing session at the Very Long Baseline Array (VLBA). The source was observed as a part of the MOJAVE survey programme with the VLBA at 15 GHz on February 5, 2007. Large amplitude variability in the correlated flux density, unexplainable in terms of the source structure, was first discovered while processing the data, and later confirmed by calibrating the antenna gains using 24 other sources observed in the experiment. The source shows variations in the correlated flux density as high as 40% on a timescale of only 2.7 hours. This places 1156+295 between the classical IDV sources and the so-called intra-hour variables. The observed variability timescale and the modulation index of 13% are consistent with interstellar scintillation by a nearby, highly turbulent scattering screen. The large modulation index at 15 GHz implies a scattering measure that is atypically high for a high galactic latitude source such as 1156+295.

We present first results from a multifrequency VLBA observations of 3C273 in 2003. The source was observed simultaneously at 5.0, 8.4, 15.3, 22.2, 43.2 and 86.2 GHz, and from this multifrequency data set, spectra of 16 emission features in the parsec scale jet were carefully constructed by using a new model-fitting based method. The measured spectra and sizes of the emission features were used to calculate the magnetic field density and the energy density of the relativistic electrons in the different parts of the parsec scale jet, independent of any equipartition assumption. We measure magnetic field density of an order of 1 Gauss in the core. The magnetic energy density in the core dominates over that of the relativistic electrons, while in the downstream region our data are roughly consistent with an equipartition. A strong gradient in the magnetic field density across the jet width, coincident with a transverse velocity structure at about 1.5 mas from the core, was found: the slower superluminal component B2 on the northern side of the jet has a magnetic field density two orders of magnitude lower than the faster southern components B3 and B4.

We report the discovery of an extremely curved jet in the radio-loud quasar PKS2136+141. Multi-frequency Very Long Baseline Array (VLBA) images show a bending jet making a turn-around of 210 degrees in the plane of the sky, which is, to our knowledge, the largest ever observed change in the position angle of an astrophysical jet. Images taken at six different frequencies, from 2.3 to 43 GHz, reveal a spiral-like trajectory, which is likely a sign of an intrinsic helical geometry. A space-VLBI image, taken with the HALCA satellite at 5 GHz and having comparable resolution to our ground-based 15 GHz data, confirms that the bend is a frequency-independent structure. VLBA monitoring data at 15 GHz, covering eight years of observations, show knots in the jet clearly deviating from ballistic motion, which suggests that the bending may be caused by a growing helical Kelvin-Helmholtz normal mode. The jet appearance suggests a helical wave at a frequency well below the "resonant" frequency of the jet, which indicates that the wave is driven by a periodic perturbation at the base of the jet. We fit the observed structure in the source with a helical twist, and we find that a simple isothermal model with a constant wave speed and wavelength gives a good fit. The measured apparent velocities indicate some degree of acceleration along the jet, which together with an observed change in the apparent half-opening angle of the jet allow us to estimate the changes in the angle between the local jet direction and our line of sight. We suggest that the jet in PKS2136+141 is distorted by a helical Kelvin-Helmholtz normal mode externally driven into the jet (e.g. by precession), and that our line of sight falls within the opening angle of the helix cone.

In this first of a series of papers describing polarimetric multifrequency Very Long Baseline Array (VLBA) monitoring of 3C 273 during a simultaneous campaign with the INTEGRAL gamma-ray satellite in 2003, we present 5 Stokes I images and source models at 7 mm. We show that a part of the inner jet (1-2 milliarcseconds from the core) is resolved in a direction transverse to the flow, and we analyse the kinematics of the jet within the first 10 mas. Based on the VLBA data and simultaneous single-dish flux density monitoring, we determine an accurate value for the Doppler factor of the parsec scale jet, and using this value with observed proper motions, we calculate the Lorentz factors and the viewing angles for the emission components in the jet. Our data indicates a significant velocity gradient across the jet with the components travelling near the southern edge being faster than the components with more northern path. We discuss our observations in the light of jet precession model and growing plasma instabilities.

We report the discovery of an extremely curved jet in the high frequency peaking GPS quasar PKS 2136+141. Our multi-frequency VLBA images show a bending jet making a turn-around of 210 degrees on the plane of the sky, which is, to our knowledge, the largest ever observed change in a structural position angle of an extragalactic jet. Images taken at six different frequencies, from 2 to 43 GHz, beautifully reveal a spiral-like trajectory. We discuss possibilities to constrain the 3-D geometry of the source and suggest that it could be used as a testbed for models describing the bending of the relativistic jets.

Light curves of 3C 279 are presented in optical (R-band), X-rays (RXTE/PCA), and gamma rays (CGRO/EGRET) for 1999 Jan-Feb and 2000 Jan-Mar. During both of those epochs the gamma-ray levels were high, and all three observed bands demonstrated substantial variation, on time scales as short as one day. Correlation analyses provided no consistent pattern, although a rather significant optical/gamma-ray correlation was seen in 1999, with a gamma-ray lag of ~2.5 days, and there are other suggestions of correlations in the light curves. For comparison, correlation analysis is also presented for the gamma-ray and X-ray light curves during the large gamma ray flare in 1996 Feb and the two gamma-bright weeks leading up to it; the correlation at that time was strong, with a gamma-ray/X-ray offset of no more than 1 day.