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

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.

The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $\gamma$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a statistical significance of 7.8 standard deviations above background. The observed VHE flux is variable on monthly time scales, with an elevated flux seen in 2018 observations. The VHE emission during this elevated state is well-characterized by a power-law spectrum with a photon index $\Gamma = 2.20 \pm 0.27$ and flux F($>315$ GeV) = ($7.6\pm 1.2_{\mathrm stat} \pm 2.3_{\mathrm syst})\times 10^{-13}$ cm$^{-2}$ s$^{-1}$, or approximately 0.7% of the Crab Nebula flux above the same threshold. 3C 264 ($z = 0.0217$) is the most distant radio galaxy detected at VHE, and the elevated state is thought to be similar to that of the famously outbursting jet in M 87. Consequently, extensive contemporaneous multi-wavelength data were acquired in 2018 at the time of the VHE high state. An analysis of these data, including VLBA, VLA, HST, Chandra and Swift observations in addition to the VERITAS data, is presented, along with a discussion of the resulting spectral energy distribution.

We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

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.

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.

We have conducted a multiwavelength survey of 42 radio loud narrow-1ine Seyfert 1 galaxies (RLNLS1s), selected by searching among all the known sources of this type and omitting those with steep radio spectra. We analyse data from radio frequencies to X-rays, and supplement these with information available from online catalogs and the literature in order to cover the full electromagnetic spectrum. This is the largest known multiwavelength survey for this type of source. We detected 90% of the sources in X-rays and found 17% at gamma rays. Extreme variability at high energies was also found, down to timescales as short as hours. In some sources, dramatic spectral and flux changes suggest interplay between a relativistic jet and the accretion disk. The estimated masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than those of blazars, while the accretion luminosities span a range from $\sim 0.01$ to $\sim 0.49$ times the Eddington limit, similar to those of quasars. The distribution of the calculated jet power spans a range from $\sim 10^{42.6}$ to $\sim 10^{45.6}$ erg s$^{-1}$, generally lower than quasars and BL Lac objects, but partially overlapping with the latter. Once normalised by the mass of the central black holes, the jet power of the three types of active galactic nuclei are consistent with each other, indicating the scalability of the jet. Despite the observational differences, the central engine of RLNLS1s is apparently quite similar to that of blazars. The historical difficulties in finding radio-loud narrow-line Seyfert 1 galaxies might be due to their low power and to intermittent jet activity.

Jets are a ubiquitous part of the accretion process, seen in a wide variety of objects ranging from active galaxies (AGN) to X-ray binary stars and even newly formed stars. AGN jets are accelerated by the supermassive black hole of their host galaxy by a coupling between the magnetic field and inflowing material. They are the source for many exciting phenomena and can profoundly influence the larger galaxy and surrounding cluster. This White Paper points out what advances can be achieved in the field by new technologies, concentrating on the zone where jets are accelerated to relativistic speeds and collimated. The ngVLA and new space VLBI missions will give higher angular resolution, sensitivity and fidelity in the radio, penetrating this zone for additional objects and allowing us to resolve fundamental questions over the physics of jet acceleration and collimation. Interferometry in other bands would allow us to probe directly flaring components. We also emphasize the need for polarimetry, which is essential to revealing the role and configuration of magnetic fields.

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.

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.

Jetted plasma outflows from active galactic nuclei (AGN) represent the most energetic phenomena in the known universe, and play a key role in regulating galaxy formation through feedback processes. The JVLA and VLBA have played an indispensable role in understanding the physics of these powerful jets and their environments, via high angular resolution full polarization imaging and astrometric studies. By bridging the current interferometric gap between the JVLA and VLBA with intermediate baselines, the ngVLA offers exciting new opportunities to explore the intermediate regions downstream of the high Lorentz factor pc-scale AGN jets imaged by VLBI, where entrainment, deceleration, collimation and particle acceleration all take place. The ability to image exceedingly faint radio emission on scales of 10s to 100s of milliarcseconds will lead to new breakthroughs in resolving some of the most pressing questions regarding the formation, structure, and evolution of AGN jets. These are in turn crucial for a more complete understanding of the formation of supermassive black holes and galaxies in the early universe and their subsequent evolution.

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.

The last decade has witnessed a steadily increasing number of observational studies concerning the rare class of radio loud narrow-line Seyfert 1 galaxies, of which several hundred are currently known. According to the current AGN paradigm, the low black hole masses and high accretion rates of narrow-line Seyfert 1 galaxies (NLSy1) should make them unlikely to launch jets, and indeed the vast majority ($\sim 90\%$) are very radio weak. The remainder, however, display a wide range of radio power, from $\simeq 10^{21}$- $10^{28}$ W/Hz. In this review I discuss recent radio imaging surveys that suggest there are three main classes of NLSy1, which cannot be easily distinguished by the standard radio-loudness parameter alone: (i) radio-weak NLSy1s without jets, (ii) mildly radio-loud NLSy1s that are a mixture of star-forming and jet-dominant AGN, and (iii) very radio loud NLSy1s with extreme properties similar to powerful jet-dominated blazars. I present updated kinematics information from the MOJAVE survey on six of the latter sources (all detected in gamma-rays by Fermi), indicating high bulk Lorentz factors and small viewing angles in three cases. Studies with the JVLA have shown that the jets of radio loud NLSy1s are likely lower- power versions of classical radio galaxies, with typical lengths of less than 10 kpc, although two very radio-loud NLSy1s have de-projected sizes of several hundred kpc. I discuss the challenges of reconciling the heterogeneous radio properties of NLSy1s with their strict optical line criteria, and near-term prospects for the discovery of larger numbers of radio-loud NLSy1s.

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 report the first results of a survey on 74 narrow-line Seyfert 1 galaxies (NLS1s) carried out in 2015 with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A-configuration. So far, this is the largest survey aimed to image the radio continuum of NLS1s. We produced radio maps in order to compare the general properties of three different samples of objects: radio-quiet NLS1s (RQNLS1s), steep-spectrum radio-loud NLS1s (S-NLS1s), and flat-spectrum radio-loud NLS1s (F-NLS1s). We find that the three classes correspond to different radio morphologies, with F-NLS1s being more compact, and RQNLS1s often showing diffuse emission on kpc scales. We also find that F-NLS1s might be low-luminosity and possibly young blazars, and that S-NLS1s are part of the parent population of F-NLS1s. Dedicated studies to RQNLS1s are needed in order to fully understand their role in the unification pictures.

Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.

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.

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 the first results of a radio survey of 79 narrow-line Seyfert 1 (NLS1) carried out with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A configuration aimed at studying the radio properties of these sources. We report the detection of extended emission in one object: Mrk 783. This is intriguing, since the radio-loudness parameter R of this object is close to the threshold between radio-quiet and radio-loud active galactic nuclei (AGN). The galaxy is one of the few NLS1 showing such an extended emission at z < 0.1. The radio emission is divided in a compact core component and an extended component, observed on both sides of the nucleus and extending from 14 kpc south-east to 12 kpc north-west. There is no sign of a collimated jet, and the shape of the extended component is similar to those of some Seyfert galaxies. The properties of the emission are compatible with a relic produced by the intermittent activity cycle of the AGN.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

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.

Relativistic jets are the most energetic manifestation of the active galactic nucleus (AGN) phenomenon. AGN jets are observed from the radio through gamma-rays and carry copious amounts of matter and energy from the sub-parsec central regions out to the kiloparsec and often megaparsec scale galaxy and cluster environs. While most spatially resolved jets are seen in the radio, an increasing number have been discovered to emit in the optical/near-IR and/or X-ray bands. Here we discuss a spectacular example of this class, the 3C 111 jet, housed in one of the nearest, double-lobed FR II radio galaxies known. We discuss new, deep Chandra and HST observations that reveal both near-IR and X-ray emission from several components of the 3C 111 jet, as well as both the northern and southern hotspots. Important differences are seen between the morphologies in the radio, X-ray and near-IR bands. The long (over 100 kpc on each side), straight nature of this jet makes it an excellent prototype for future, deep observations, as it is one of the longest such features seen in the radio, near-IR/optical and X-ray bands. Several independent lines of evidence, including the X-ray and broadband spectral shape as well as the implied velocity of the approaching hotspot, lead us to strongly disfavor the EC/CMB model and instead favor a two-component synchrotron model to explain the observed X-ray emission for several jet components. Future observations with NuSTAR, HST, and Chandra will allow us to further constrain the emission mechanisms.

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)

Narrow-line Seyfert 1 galaxies (NLS1s) are active galactic nuclei (AGN) recently identified as a new class of $\gamma$-ray sources. The high energy emission is explained by the presence of a relativistic jet observed at small angles, just like in the case of blazars. When the latter are observed at larger angles they appear as radio-galaxies, but an analogue parent population for beamed NLS1s has not yet been determined. In this work we analyze this problem by studying the physical properties of three different samples of parent sources candidates: steep-spectrum radio-loud NLS1s, radio-quiet NLS1s, and disk-hosted radio-galaxies, along with compact steep-spectrum sources. In our approach, we first derived black hole mass and Eddington ratio from the optical spectra, then we investigated the interaction between the jet and the narrow-line region from the [O III] $\lambda\lambda$4959,5007 lines. Finally, the radio luminosity function allowed us to compare their jet luminosity and hence determine the relations between the samples.

We report about recent updates of broad-band properties of radio-loud narrow-line Seyfert 1 galaxies.

Using high-resolution radio imaging with VLBI techniques, the TANAMI program has been observing the parsec-scale radio jets of southern (declination south of -30{\deg}) gamma-ray bright AGN simultaneously with Fermi/LAT monitoring of their gamma-ray emission. We present the radio and gamma-ray properties of the TANAMI sources based on one year of contemporaneous TANAMI and Fermi/LAT data. A large fraction (72%) of the TANAMI sample can be associated with bright gamma-ray sources for this time range. Association rates differ for different optical classes with all BL Lacs, 76% of quasars and just 17% of galaxies detected by the LAT. Upper limits were established on the gamma-ray flux from TANAMI sources not detected by LAT. This analysis led to the identification of three new Fermi sources whose detection was later confirmed. The gamma-ray and radio luminosities are related by $L_\gamma \propto L_r^{0.89+-0.04}$. The brightness temperatures of the radio cores increase with the average gamma-ray luminosity, and the presence of brightness temperatures above the inverse Compton limit implies strong Doppler boosting in those sources. The undetected sources have lower gamma/radio luminosity ratios and lower contemporaneous brightness temperatures. Unless the Fermi/LAT-undetected blazars are strongly gamma-ray-fainter than the Fermi/LAT-detected ones, their gamma-ray luminosity should not be significantly lower than the upper limits calculated here.

We report on a multiwavelength survey of a sample of 42 flat-spectrum radio-loud narrow-line Seyfert 1 galaxies (RLNLS1s). This is the largest known sample of this type of active galactic nucleus (AGN) to date. We found that 17% of sources were detected at high-energy gamma rays (E>100 MeV), and 90% at X-rays (0.3-10 keV). The masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than the values of blazars. The disk luminosities are about 1-49% of the Eddington value, with one outlier at 0.3%, comparable with the luminosities observed in flat-spectrum radio quasars (FSRQs). The jet powers are $\sim 10^{42-46}$ erg s$^{-1}$, comparable with BL Lac Objects, yet relatively smaller than FSRQs. However, once renormalized by the mass of the central black hole, the jet powers of RLNLS1s, BL Lacs, and FSRQs are consistent each other, indicating the scalability of the jets. We found episodes of extreme variability at high energies on time scales of hours. In some cases, dramatic spectral and flux changes are interpreted as the interplay between the relativistic jet and the accretion disk. We conclude that, despite the distinct observational properties, the central engines of RLNLS1s are similar to those of blazars.

The aim of this work is to understand the nature of the parent population of beamed narrow-line Seyfert 1 galaxies (NLS1s), by studying the physical properties of three parent candidates samples: steep-spectrum radio-loud NLS1s, radio-quiet NLS1s and disk-hosted radio-galaxies. In particular, we focused on the black hole mass and Eddington ratio distribution and on the interactions between the jet and the narrow-line region.

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.

We use the complete MOJAVE 1.5 Jy sample of active galactic nuclei (AGN) to examine the gamma-ray detection statistics of the brightest radio-loud blazars in the northern sky. We find that 23% of these AGN were not detected above 0.1 GeV by the Fermi LAT during the 4-year 3FGL catalog period partly because of an instrumental selection effect, and partly due to their lower Doppler boosting factors. Blazars with synchrotron peaks in their spectral energy distributions located below $10^{13.4}$ Hz also tend to have high-energy peaks that lie below the 0.1 GeV threshold of the LAT, and are thus less likely to be detected by Fermi. The non-detected AGN in the 1.5 Jy sample also have significantly lower 15 GHz radio modulation indices and apparent jet speeds, indicating that they have lower than average Doppler factors. Since the effective amount of relativistic Doppler boosting is enhanced in gamma-rays (particularly in the case of external inverse-Compton scattering), this makes them less likely to appear in the 3FGL catalog. Based on their observed properties, we have identified several bright radio-selected blazars that are strong candidates for future detection by Fermi.

We present multiwavelength imaging observations of PKS 1045-188, 8C 1849+670, and PKS 2216-038, three radio-loud active galactic nuclei from the MOJAVE-Chandra Sample that straddle the Fanaroff-Riley (FR) boundary between low- and high-power jets. These hybrid sources provide an excellent opportunity to study jet emission mechanisms and the influence of the external environment. We used archival VLA observations, and new Hubble and Chandra observations to identify and study the spectral properties of five knots in PKS 1045-188, two knots in 8C 1849+670, and three knots in PKS 2216-038. For the seven X-ray visible knots, we constructed and fit the broadband spectra using synchrotron and inverse Compton/cosmic microwave background (IC/CMB) emission models. In all cases, we found that the lack of detected optical emission ruled out the X-ray emission from the same electron population that produces radio emission. All three sources have high total extended radio power, similar to that of FR II sources. We find this is in good agreement with previously studied hybrid sources, where high-power hybrid sources emit X-rays via IC/CMB and the low-power hybrid sources emit X-rays via synchrotron emission. This supports the idea that it is total radio power rather than FR morphology that determines the X-ray emission mechanism. We found no significant asymmetries in the diffuse X-ray emission surrounding the host galaxies. Sources PKS 1045-188 and 8C 1849+670 show significant differences in their radio and X-ray termination points, which may result from the deceleration of highly relativistic bulk motion.

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.

Flat-spectrum radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s) are a recently discovered class of $\gamma$-ray emitting Active Galactic Nuclei (AGN), that exhibit some blazar-like properties which are explained with the presence of a relativistic jet viewed at small angles. When blazars are observed at larger angles they appear as radio-galaxies, and we expect to observe an analogue parent population for beamed NLS1s. However, the number of known NLS1s with the jet viewed at large angles is not enough. Therefore, we tried to understand the origin of this deficit. Current hypotheses about the nature of parent sources are steep-spectrum radio-loud NLS1s, radio-quiet NLS1s and disk-hosted radio-galaxies. To test these hypotheses we built three samples of candidate sources plus a control sample, and calculated their black hole mass and Eddington ratio using their optical spectra. We then performed a Kolmogorov-Smirnov statistical test to investigate the compatibility of our different samples with a beamed population. Our results indicate that, when the inclination angle increases, a beamed source appears as a steep-spectrum radio-loud NLS1, or possibly even as a disk-hosted radio-galaxy with low black hole mass and high Eddington ratio. Further investigations, involving larger complete samples and observations at radio frequency, are needed to understand the incidence of disk-hosted radio-galaxies in the parent population, and to assess whether radio-quiet NLS1s can play a role, as well.

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.

We have discovered kiloparsec-scale extended radio emission in three narrow-line Seyfert 1 galaxies (NLS1s) in sub-arcsecond resolution 9 GHz images from the Karl G. Jansky Very Large Array (VLA). We find all sources show two-sided, mildly core-dominated jet structures with diffuse lobes dominated by termination hotspots. These span 20-70 kpc with morphologies reminiscent of FR II radio galaxies, while the extended radio luminosities are intermediate between FR I and FR II sources. In two cases the structure is linear, while a $45^{\circ}$ bend is apparent in the third. Very Long Baseline Array images at 7.6 GHz reveal parsec-scale jet structures, in two cases with extended structure aligned with the inner regions of the kiloparsec-scale jets. Based on this alignment, the ratio of the radio core luminosity to the optical luminosity, the jet/counter-jet intensity and extension length ratios, and moderate core brightness temperatures ($\lesssim10^{10}$ K), we conclude these jets are mildly relativistic ($\beta\lesssim0.3$, $\delta\sim1$-$1.5$) and aligned at moderately small angles to the line of sight (10-15$^{\circ}$). The derived kinematic ages of $\sim10^6$-$10^7$ y are much younger than radio galaxies but comparable to other NLS1s. Our results increase the number of radio-loud NLS1s with known kiloparsec-scale extensions from seven to ten and suggest that such extended emission may be common, at least among the brightest of these sources.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

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.

Results from Chandra-HST-VLA observations of 13 hybrid sources are presented. Data from ten sources in the literature are analysed along with new data from three hybrid blazars belonging to the MOJAVE sample. Studies of such hybrid sources displaying both FRI and FRII jet characteristics could provide the key to resolving the long-standing Fanaroff-Riley dichotomy issue. A majority of the 13 hybrid sources show FRII-like total radio powers, i.e., they are "hybrid" in radio morphology but not in total radio power. VLBI observations of ten of the 13 sources show that the X-ray jet is on the same side as the one-sided VLBI jet. X-rays are therefore emitted from relativistically-boosted approaching jets. This is consistent with the X-ray emission being IC/CMB in origin in the majority of sources. It is not completely clear from our study that asymmetries in the surrounding medium can create hybrid sources. Hybrid radio morphologies could also be indicative of intrinsically asymmetric jets.

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

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

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

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)

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

We present the preliminary results of a survey performed with Swift to observe a sample of radio-loud Narrow-Line Seyfert 1 Galaxies (RLNLS1s). Optical-to-X-ray data from Swift are complemented with gamma-ray observations from Fermi/LAT and radio measurements available in the literature. The comparison with a sample of bright Fermi blazars indicates that RLNLS1s seem to be the low-power tail of the distribution.

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

We report the detection of a new TeV gamma-ray source, VER J0521+211, based on observations made with the VERITAS imaging atmospheric Cherenkov telescope array. These observations were motivated by the discovery of a cluster of >30GeV photons in the first year of Fermi-LAT observations. VER J0521+211 is relatively bright at TeV energies, with a mean photon flux of 1.93 +/- 0.13_stat +/- 0.78_sys 10^-11 cm-2 s-1 above 0.2 TeV during the period of the VERITAS observations. The source is strongly variable on a daily timescale across all wavebands, from optical to TeV, with a peak flux corresponding to ~0.3 times the steady Crab Nebula flux at TeV energies. Follow-up observations in the optical and X-ray bands classify the newly-discovered TeV source as a BL Lac-type blazar with uncertain redshift, although recent measurements suggest z=0.108. VER J0521+211 exhibits all the defining properties of blazars in radio, optical, X-ray, and gamma-ray wavelengths.

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)

[Abridged] Context. To construct and interpret the spectral energy distribution (SED) of BL Lacertae objects, simultaneous broad-band observations are mandatory. Aims. We present the results of a dedicated multi-wavelength study of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 1ES 2344+514 by means of a pre-organised campaign. Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were complemented by Mets\"ahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. Results. 1ES 2344+514 was found at very low flux states in both X-rays and very high energy gamma rays. Variability was detected in the low frequency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was possibly caused by shock acceleration characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over time scales of eleven years, contrary to previous findings. There appears to be no significant correlation between the 15 GHz and R-band monitoring light curves. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE energies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. Consequently it remains unclear whether a so-called quiescent state was found in this campaign.

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.

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 report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

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 present the first collective evidence that Fermi-detected jets of high kinetic power (L_kin) are dominated by inverse Compton emission from upscattered external photons. Using a sample with a broad range in orientation angle, including radio galaxies and blazars, we find that very high power sources (L_kin > 10^45.5 erg s^{-1}) show a significant increase in the ratio of inverse Compton to synchrotron power (Compton dominance) with decreasing orientation angle, as measured by the radio core dominance and confirmed by the distribution of superluminal speeds. This increase is consistent with beaming expectations for external Compton (EC) emission, but not for synchrotron self Compton (SSC) emission. For the lowest power jets (L_kin < 10^43.5 erg s^{-1}), no trend between Compton and radio core dominance is found, consistent with SSC. Importantly, the EC trend is not seen for moderately high power flat spectrum radio quasars with strong external photon fields. Coupled with the evidence that jet power is linked to the jet speed (Kharb et al. 2010), this finding suggests that external photon fields become the dominant source of seed photons in the jet comoving frame only for the faster and therefore more powerful jets.

In recent work, we have identified two sub-populations of radio-loud AGN which appear to be distinguished by jet structure, where low-efficiency accreting systems produce `weak' jets which decelerate more rapidly than the `strong' jets of black holes accreting near the Eddington limit. The two classes are comprised of: (1) The weak jet sources, corresponding to FR I radio galaxies, having a decelerating or spine-sheath jet with velocity gradients, and (2) The strong jet sources, having fast, collimated jets, and typically displaying strong emission lines. The dichotomy in the \nu_peak-L_peak plane can be understood as a `broken power sequence' in which jets exist on one branch or the other based on the particular accretion mode. We suggest that the intrinsic kinetic power (as measured by low-frequency, isotropic radio emission), the orientation, and the accretion rate of the SMBH system are the the fundamental axes needed for unification of radio-loud AGN by studying a well-characterized sample of several hundred Fermi-detected jets. Finally, we present very recent findings that the most powerful strong jets produce gamma-rays by external Compton rather than SSC emission, placing the dissipation region in these strong jets at a radius inside the BLR and/or molecular torus.

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.

In examining a select sample of over 200 blazars of known jet kinetic power (L_kin) and well-characterized SEDs, we found (Meyer et al., 2011) that Intermediate synchrotron-peaking (ISP) blazars may have lower gamma-ray output than high synchrotron-peaking (HSP) blazars of similar L_kin, consistent with our hypothesis that ISP blazars are less-beamed versions of HSP blazars, rather than a distinct population. Further, by using the radio core dominance as a measure of relative beaming, we find that gamma-ray luminosity depends on beaming in a consistent way for blazars ranging over all jet kinetic powers (10^42 - 10^46 ergs/s). We re-examine the gamma-ray properties of this core sample of blazars using the 1-year LAT catalog (Abdo et al., 2010). We find that for weak jets, the ratio of inverse Compton to synchrotron emission remains constant with increased beaming, consistent with an SSC model for the jet emission, while the most powerful jets show a strong increase in Compton dominance with orientation, consistent with an external Compton (EC) emission model.

We present results from deep (70 ks) Chandra ACIS observations and Hubble Space Telescope ACS F475W observations of two highly optically polarized quasars belonging to the MOJAVE blazar sample, viz., PKS B0106+013 and 1641+399 (3C345). These observations reveal X-ray and optical emission from the jets in both sources. X-ray emission is detected from the entire length of the 0106+013 radio jet, which shows clear bends or wiggles - the X-ray emission is brightest at the first prominent kpc jet bend. A picture of a helical kpc jet with the first kpc-scale bend representing a jet segment moving close(r) to our line of sight, and getting Doppler boosted at both radio and X-ray frequencies, is consistent with these observations. The X-ray emission from the jet end however peaks at about 0.4" (~3.4 kpc) upstream of the radio hot spot. Optical emission is detected both at the X-ray jet termination peak and at the radio hot spot. The X-ray jet termination peak is found upstream of the radio hot spot by around 0.2" (~1.3 kpc) in the short projected jet of 3C345. HST optical emission is seen in an arc-like structure coincident with the bright radio hot spot, which we propose is a sharp (apparent) jet bend instead of a terminal point, that crosses our line of sight and consequently has a higher Doppler beaming factor. A weak radio hot spot is indeed observed less than 1" downstream of the bright radio hot spot, but has no optical or X-ray counterpart. By making use of the pc-scale radio and the kpc-scale radio/X-ray data, we derive constraints on the jet Lorentz factors (Gamma_jet) and inclination angles (theta): for a constant jet speed from pc- to kpc-scales, we obtain a Gamma_jet of ~70 for 0106+013, and ~40 for 3C345. On relaxing this assumption, we derive a Gamma_jet of ~2.5 for both the sources. Upper limits on theta of ~13 degrees are obtained for the two quasars. (ABRIDGED)

We recently argued (Meyer 2011) that the collective properties of radio loud active galactic nuclei point to the existence of two families of sources, one of powerful sources with single velocity jets and one of weaker jets with significant velocity gradients in the radiating plasma. These families also correspond to different accretion modes and therefore different thermal and emission line intrinsic properties: powerful sources have radiatively efficient accretion disks, while in weak sources accretion must be radiatively inefficient. Here, after we briefly review of our recent work, we present the following findings that support our unification scheme: (i) along the broken sequence of aligned objects, the jet kinetic power increases. (ii) in the powerful branch of the sequence of aligned objects the fraction of BLLs decreases with increasing jet power. (iii) for powerful sources, the fraction of BLLs increases for more un-aligned objects, as measured by the core to extended radio emission. Our results are also compatible with the possibility that a given accretion power produces jets of comparable kinetic power.

We report about the multiwavelength campaign on the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by high activity as measured by Fermi/LAT. The peak luminosity in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^48 erg/s, a level comparable to the most powerful blazars. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

Aims. We use a sample of 83 core-dominated active galactic nuclei (AGN) selected from the MOJAVE (Monitoring of Jets in AGN with VLBA Experiments) radio-flux-limited sample and detected with the Fermi Large Area Telescope (LAT) to study the relations between non-simultaneous radio, optical, and gamma-ray measurements. Methods. We perform a multi-band statistical analysis to investigate the relations between the emissions in different bands and reproduce these relations by modeling of the spectral energy distributions of blazars. Results. There is a significant correlation between the gamma-ray luminosity and the optical nuclear and radio (15 GHz) luminosities of blazars. We report a well defined positive correlation between the gamma-ray luminosity and the radio-optical loudness for quasars and BL Lacertae type objects (BL Lacs). A strong positive correlation is found between the radio luminosity and the gamma-ray-optical loudness for quasars, while a negative correlation between the optical luminosity and the gamma-ray-radio loudness is present for BL Lacs. Modeling of these correlations with a simple leptonic jet model for blazars indicates that variations of the accretion disk luminosity (and hence the jet power) is able to reproduce the trends observed in most of the correlations. To reproduce all observed correlations, variations of several parameters, such as the accretion power, jet viewing angle, Lorentz factor, and magnetic field of the jet, are required.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

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

We revisit the concept of a blazar sequence that relates the synchrotron peak frequency ({\nu}peak) in blazars with synchrotron peak luminosity (Lpeak, in {\nu}L{\nu}) using a large sample of radio-loud AGN. We present observational evidence that the blazar sequence is formed from two populations in the synchrotron {\nu}peak - Lpeak plane, each forming an upper edge to an envelope of progressively misaligned blazars, and connecting to an adjacent group of radio galaxies having jets viewed at much larger angles to the line of sight. When binned by jet kinetic power (Lkin; as measured through a scaling relationship with extended radio power), we find that radio core dominance decreases with decreasing synchrotron Lpeak, revealing that sources in the envelope are generally more misaligned. We find population-based evidence of velocity gradients in jets at low kinetic powers (~ 10^42-10^44.5 erg/s), corresponding to FR I radio galaxies and most BL Lacs. These low jet power 'weak jet' sources, thought to exhibit radiatively inefficient accretion, are distinguished from the population of non-decelerating, low synchrotron-peaking (LSP) blazars and FR II radio galaxies ('strong' jets) which are thought to exhibit radiatively efficient accretion. The two-population interpretation explains the apparent contradiction of the existence of highly core-dominated, low-power blazars at both low and high synchrotron peak frequencies, and further implies that most intermediate synchrotron peak (ISP) sources are not intermediate in intrinsic jet power between LSP and high synchrotron-peaking (HSP) sources, but are more misaligned versions of HSP sources with similar jet powers.

The Chandra X-ray observatory has proven to be a vital tool for studying high-energy emission processes in jets associated with Active Galactic Nuclei (AGN).We have compiled a sample of 27 AGN selected from the radio flux-limited MOJAVE (Monitoring of Jets in AGN with VLBA Experiments) sample of highly relativistically beamed jets to look for correlations between X-ray and radio emission on kiloparsec scales. The sample consists of all MOJAVE quasars which have over 100 mJy of extended radio emission at 1.4 GHz and a radio structure of at least 3" in size. Previous Chandra observations have revealed X-ray jets in 11 of 14 members of the sample, and we have carried out new observations of the remaining 13 sources. Of the latter, 10 have Xray jets, bringing the overall detection rate to ~ 78%. Our selection criteria, which is based on highly compact, relativistically beamed jet emission and large extended radio flux, thus provides an effective method of discovering new X-ray jets associated with AGN. The detected X-ray jet morphologies are generally well correlated with the radio emission, except for those displaying sharp bends in the radio band. The X-ray emission mechanism for these powerful FR II (Fanaroff-Riley type II) jets can be interpreted as inverse Compton scattering off of cosmic microwave background (IC/CMB) photons by the electrons in the relativistic jets. We derive viewing angles for the jets, assuming a non-bending, non-decelerating model, by using superluminal parsec scale speeds along with parameters derived from the inverse Compton X-ray model. We use these angles to calculate best fit Doppler and bulk Lorentz factors for the jets, as well as their possible ranges, which leads to extreme values for the bulk Lorentz factor in some cases. When both the non-bending and non-decelerating assumptions are relaxed [abridged]

We report on a multiwavelength campaign on the radio-loud Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.5846) performed in 2010 July-September and triggered by a high-energy gamma-ray outburst observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. The peak flux in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^-6 ph cm^-2 s^-1, corresponding to an observed luminosity of 10^48 erg s^-1. Although the source was too close to the Sun position to organize a densely sampled follow-up, it was possible to gather some multiwavelength data that confirmed the state of high activity across the sampled electromagnetic spectrum. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

We study the correlations between the Very Long Baseline Array radio emission at 15 GHz, extended emission at 151 MHz, and optical nuclear emission at 5100 AA for a complete sample of 135 compact jets. We use the partial Kendall's tau correlation analysis to check the link between radio properties of parsec-scale jets and optical luminosities of host AGN. We find a significant positive correlation for 99 quasars between optical nuclear luminosities and total radio (VLBA) luminosities of unresolved cores at 15 GHz originated at milliarcseconds scales. For 18 BL Lacs, the optical continuum emission correlates with the radio emission of the jet at 15 GHz. We suggest that the radio and optical emission are beamed and originate in the innermost part of the sub--parsec-scale jet in quasars. Analysis of the relation between the apparent speed of the jet and the optical nuclear luminosity at 5100 AA supports the relativistic beaming model for the optical emission generated in the jet, and allows the peak values of the intrinsic optical luminosity of the jet and its Lorentz factor to be estimated for the populations of quasars, BL Lacs, and radio galaxies. The radio-loudness of quasars is found to increase at high redshifts, which can be a result of lower efficiency of the accretion in AGN having higher radio luminosities. A strong positive correlation is found between the intrinsic kinetic power of the jet and the apparent luminosities of the total and the unresolved core emission of the jet at 15 GHz. This correlation is interpreted in terms of intrinsically more luminous parsec-scale jet producing more luminous extended structure which is detectable at low radio frequencies, 151 MHz. A possibility that the low frequency radio emission is relativistically beamed in superluminal AGN and therefore correlates with radio luminosity of the jet at 15 GHz can not be ruled out (abridged).

We report the detection of a non-zero time delay between radio emission measured by the VLBA at 15.4 GHz and gamma-ray radiation (gamma-ray leads radio) registered by the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope for a sample of 183 radio and gamma-ray bright active galactic nuclei (AGNs). For the correlation analysis we used 100 MeV - 100 GeV gamma-ray photon fluxes, taken from monthly binned measurements from the first Fermi LAT catalog, and 15.4 GHz radio flux densities from the MOJAVE VLBA program. The correlation is most pronounced if the core flux density is used, strongly indicating that the gamma-ray emission is generated within the compact region of the 15 GHz VLBA core. Determining the Pearson's r and Kendall's tau correlation coefficients for different time lags, we find that for the majority of sources the radio/gamma-ray delay ranges from 1 to 8 months in the observer's frame and peaks at about 1.2 months in the source's frame. We interpret the primary source of the time delay to be synchrotron opacity in the nuclear region.

We introduce the TANAMI program (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) which is monitoring an initial sample of 43 extragalactic jets located south of -30 degrees declination at 8.4 GHz and 22 GHz since 2007. All aspects of the program are discussed. First epoch results at 8.4 GHz are presented along with physical parameters derived therefrom. We present first epoch images for 43 sources, some observed for the first time at milliarcsecond resolution. Parameters of these images as well as physical parameters derived from them are also presented and discussed. These and subsequent images from the TANAMI survey are available at http://pulsar.sternwarte.uni-erlangen.de/tanami/ We obtain reliable, high dynamic range images of the southern hemisphere AGN. All the quasars and BL Lac objects in the sample have a single-sided radio morphology. Galaxies are either double-sided, single-sided or irregular. About 28% of the TANAMI sample has been detected by LAT during its first three months of operations. Initial analysis suggests that when galaxies are excluded, sources detected by LAT have larger opening angles than those not detected by LAT. Brightness temperatures of LAT detections and non-detections seem to have similar distributions. The redshift distributions of the TANAMI sample and sub-samples are similar to those seen for the bright gamma-ray AGN seen by LAT and EGRET but none of the sources with a redshift above 1.8 have been detected by LAT.

Magnetically-driven non-stationary acceleration of jets in AGNs results in the leading parts of the flow been accelerated to much higher Lorentz factors than in the case of steady state acceleration with the same parameters. The higher Doppler-boosted parts of the flow may dominate the high energy emission of blazar jets. We suggest that highly variable GeV and TeV emission in blazars is produced by the faster moving leading edges of highly magnetized non-stationary ejection blobs, while the radio data trace the slower-moving bulk flow. Model predictions compare favorably with the latest Fermi gamma-ray and MOJAVE radio VLBI results.

A bright feature 80 pc away from the core in the powerful jet of M87 shows highly unusual properties. Earlier radio, optical and X-ray observations have shown that this feature, labeled HST-1, is superluminal, and is possibly connected with the TeV flare detected by HESS in 2005. It has been claimed that this feature might have a blazar nature, due to these properties. To examine the possible blazar-like nature of HST-1, we analyzed lambda 2 cm VLBA archival data from dedicated full-track observations and the 2 cm survey/MOJAVE VLBI monitoring programs obtained between 2000 and 2009. Applying VLBI wide-field imaging techniques, the HST-1 region was imaged at milliarcsecond resolution. Here we present the first 2 cm VLBI detection of this feature in observations from early 2003 to early 2007, and analyze its evolution over this time. Using the detections of HST-1, we find that the projected apparent speed is 0.61 +/- 0.31 c. A comparison of the VLA and VLBA flux densities of this feature indicate that is mostly resolved on molliarcsecond scales. This feature is optically thin between lambda 2 cm and lambda 20 cm. We do not find evidence of a blazar nature for HST-1.

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

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.

We present a summary of the observation strategy of TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry), a monitoring program to study the parsec-scale structure and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Australian Long Baseline Array (LBA) and the trans-oceanic antennas Hartebeesthoek, TIGO, and O'Higgins. TANAMI is focusing on extragalactic sources south of -30 degrees declination with observations at 8.4 GHz and 22 GHz every ~2 months at milliarcsecond resolution. The initial TANAMI sample of 43 sources has been defined before the launch of the Fermi Gamma Ray Space Telescope to include the most promising candidates for bright gamma-ray emission to be detected with its Large Area Telescope (LAT). Since November 2008, we have been adding new sources to the sample, which now includes all known radio- and gamma-ray bright AGN of the Southern Hemisphere. The combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN and the nature of relativistic jets.

We present the results of a study on the 1.4 GHz kpc-scale radio emission in the complete flux density limited MOJAVE sample, comprising 135 radio-loud AGNs. While extended emission is detected in the majority of the sources, about 7% of the sources exhibit only radio core emission. Many BL Lacs exhibit extended radio power and kpc-scale morphology typical of powerful FRII jets, while a substantial number of quasars possess radio powers intermediate between FRIs and FRIIs. This poses challenges to the simple radio-loud unified scheme, which links BL Lacs to FRIs and quasars to FRIIs. We find a significant correlation between extended radio emission and pc-scale jet speeds: the more radio powerful sources possess faster jets. This indicates that the 1.4 GHz (or low frequency) radio emission is indeed related to jet kinetic power. Various properties such as extended radio power and apparent pc-scale jet speeds vary smoothly between different blazar subclasses, suggesting that, at least in terms of radio jet properties, the distinction between quasars and BL Lac objects, at an emission-line equivalent width of 5 Angstrom is essentially an arbitrary one. Based on the assumption that the extended radio luminosity is affected by the kpc-scale environment, we define the ratio of extended radio power to absolute optical magnitude as a proxy for environmental effects. Trends with this parameter suggest that the pc-scale jet speeds and the pc-to-kpc jet misalignments are not affected by the large-scale environment, but are more likely to depend upon factors intrinsic to the AGN, or its local pc-scale environment. We suggest that some of the extremely misaligned MOJAVE blazar jets could be "hybrid" morphology sources, with an FRI jet on one side and an FRII jet on the other. (Abridged)

The TANAMI (Tracking AGN with Austral Milliarcsecond Interferometry) and associated programs provide comprehensive radio monitoring of extragalactic gamma-ray sources south of declination -30 degrees. Joint quasi-simultaneous observations between the Fermi Gamma-ray Space Telescope and ground based observatories allow us to discriminate between competing theoretical blazar emission models. High resolution VLBI observations are the only way to spatially resolve the sub-parsec level emission regions where the high-energy radiation originates. The gap from radio to gamma-ray energies is spanned with near simultaneous data from the Swift satellite and ground based optical observatories. We present early results from the TANAMI program in the context of this panchromatic suite of observations.

TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) is a monitoring program to study the parsec-scale structures and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Long Baseline Array and associated telescopes. Extragalactic jets south of -30 degrees declination are observed at 8.4 GHz and 22 GHz every two months at milliarcsecond resolution. The initial TANAMI sample is a hybrid radio and gamma-ray selected sample since the combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

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.

The radio-loud active galactic nucleus in M 87 hosts a powerful jet fueled by a super-massive black hole in its center. A bright feature 80 pc away from the M 87 core has been reported to show superluminal motions, and possibly to be connected with a TeV flare observed around 2005. To complement these studies and to understand the nature of this feature, we analyzed 2 cm VLBI data from 15 observing runs between 2000 and 2009. This feature is successfully detected at the milli-Jansky level from 2003 to 2007. Our detections show that its milli-arcsecond structure appears to be extended with a steep spectrum, and no compact or rapidly moving features are observed. Our results do not favor a blazar scenario for this feature.

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.

We report the discovery of high-energy (E>100 MeV) gamma-ray emission from NGC 1275, a giant elliptical galaxy lying at the center of the Perseus cluster of galaxies, based on observations made with the Large Area Telescope (LAT) of the Fermi Gamma ray Space Telescope. The positional center of the gamma-ray source is only ~3' away from the NGC 1275 nucleus, well within the 95% LAT error circle of ~5'.The spatial distribution of gamma-ray photons is consistent with a point source. The average flux and power-law photon index measured with the LAT from 2008 August 4 to 2008 December 5 are F_gamma = (2.10+-0.23)x 10^{-7} ph (>100 MeV) cm^{-2} s^{-1} and Gamma = 2.17+-0.05, respectively. The measurements are statistically consistent with constant flux during the four-month LAT observing period.Previous EGRET observations gave an upper limit of F_gamma < 3.72x 10 ^{-8} ph (>100 MeV) cm^{-2} s^{-1} to the gamma-ray flux from NGC 1275. This indicates that the source is variable on timescales of years to decades, and therefore restricts the fraction of emission that can be produced in extended regions of the galaxy cluster. Contemporaneous and historical radio observations are also reported. The broadband spectrum of NGC 1275 is modeled with a simple one-zone synchrotron/synchrotron self-Compton model and a model with a decelerating jet flow.

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

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

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

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

A bright feature 100 pc away from the core in the powerful jet of M 87 shows mysterious properties. Earlier radio, optical and X-ray observations have shown that this feature, labelled HST-1, is superluminal, and is possibly connected with the TeV flare detected by HESS in 2005. To examine the possible blazar-like nature of HST-1, we analyzed 2 cm VLBA data from dedicated full-track observations and the 2 cm survey/MOJAVE VLBI monitoring programs observed from 2000 to 2008. Applying wide-field imaging techniques, the HST-1 region was imaged at milliarcsecond resolutions. Here we present the first 15 GHz VLBI detection of this feature and discuss the connection between our radio findings and the TeV detection.

We show that using either the method of Page & Carrera or the well-known $1/V_a$ method to construct the binned luminosity function (LF) of a flux limited sample of Active Galactic Nuclei (AGN) can produce an artificial flattening (or steepening in the case of negative evolution) of the binned LF for bins intersected by the flux cutoff of the sample. This effect is more pronounced for samples with steep and strongly evolving parent LFs but is still present even for non-evolving LFs. As a result of this distortion of the true LF, fitting a model LF to binned data may lead to errors in the estimation of the parameters and may even prompt the erroneous use of broken power law functions. We compute the expected positions of apparent breaks in the binned LF. We show that these spurious breaks in the binned LFs can be avoided if the binning is done in the flux--redshift plane instead of the typically used luminosity--redshift plane. Binning in the flux--redshift plane can be used in conjunction with the binning in the luminosity--redshift plane to test for real breaks in the binned LFs and to identify the features that are the result of binning biases. We illustrate this effect for most typical forms of luminosity dependence and redshift evolution and show how the proposed method helps address this problem. We also apply this method to the MOJAVE AGN sample and show that it eliminates an apparent break in the binned LF.

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

We re-examine the differences observed in the pc-scale magnetic field geometry of high and low optical polarization Quasars (HPQs, LPRQs) using the MOJAVE sample. We find that, as previously reported, HPQ jets exhibit predominantly transverse B fields while LPRQ jets tend to display longitudinal B fields. We attempt to understand these results along with the different B field geometry observed in the low and high energy peaked BL Lacs (LBLs, HBLs) using a simple picture wherein the spinning central black holes in these AGNs influence the speed and strength of the jet components (spine, sheath). Higher spin rates in HPQs compared to LPRQs and LBLs compared to HBLs could explain the different total radio powers, VLBI jet speeds, and the observed B field geometry in these AGN classes.

The HEGRA gamma-ray source TeV J2032+4130 is considered the prototypical 'dark accelerator', since it was the first TeV source detected with no firm counterparts at lower frequencies. The Whipple collaboration observed this source in 2003-5 and the emission hotspot appears displaced about 9 arcminutes to the northeast of the HEGRA position, though given the large positional uncertainties the HEGRA and Whipple positions are consistent. Here we report on Westerbork Synthesis Radio Telescope (WSRT), Very Large Array (VLA), Chandra and INTEGRAL data covering the locations of the Whipple and HEGRA hotspots. We confirm a dual-lobed radio source (also see Marti et al., 2007) coincident with the Whipple hotspot, as well as a weak, partially non-thermal shell-like object, with a location and morphology very similar to the HEGRA source, in our WSRT and mosaicked VLA datasets, respectively. Due to its extended nature, it is likely that the latter structure is a more plausible counterpart of the reported very high energy (VHE) gamma-ray emissions in this region. If so, TeV J2032+4130 may not be a 'dark accelerator' after all. Further observations with the new generation of imaging Cherenkov telescopes are needed to pin down the precise location and morphology of the TeV emission region and thus clear up the confusion over its possible lower frequency counterparts.

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

(Abridged) We use a complete sample of active galactic nuclei (AGN) selected on the basis of relativistically beamed 15 GHz radio flux density to derive the parent radio luminosity function (RLF) of bright radio-selected blazar cores. We use a maximum likelihood method to fit a beamed RLF to the observed data and thereby recover the parameters of the intrinsic (unbeamed) RLF. We analyze two subsamples of the MOJAVE sample: the first contains only objects of known FR II class, with a total of 103 sources, and the second subsample adds 24 objects of uncertain FR class for a total of 127 sources. Both subsamples exclude four known FR I radio galaxies and two gigahertz-peaked spectrum sources. We obtain good fits to both subsamples using a single power law intrinsic RLF with pure density evolution function. We find that a previously reported break in the observed MOJAVE RLF actually arises from using incomplete bins (because of the luminosity cutoff) across a steep and strongly evolving RLF, and does not reflect a break in the intrinsic RLF. The derived space density of the parent population of the FR II sources from the MOJAVE sample (with L>1.3e25 W/Hz) is approximately 1600/Gpc^3.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. Approximately twenty-five percent of the mission time was dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz. This paper, the fifth in the series, presents images and models for the remaining 140 sources not included in Paper III, which contained 102 sources. For most sources, the plots of the uv-coverage, the visibility amplitude versus uv-distance, and the high resolution image are presented. Model fit parameters to the major radio components are determined, and the brightness temperature of the core component for each source is calculated. The brightness temperature distributions for all of the sources in the VSOP AGN survey are discussed.

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

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

Observations made with the VLBA have led to fundamental advances in our understanding of how radio jets in AGN evolve from parsec-scales out to distances exceeding several hundred kiloparsecs. In this review I discuss current models of young radio source evolution, as well as the observational evidence for a rapid change in jet properties on scales of ~1 kpc. A central topic of current debate is the relative importance of intermittent jet fueling versus jet-environment interactions in causing a drop-off in powerful radio sources at this critical evolutionary stage. Recent 3-D hydrodynamical jet simulations suggest that dense environments and cloud collisions can temporarily stifle, but not completely halt powerful relativistic jets. Several VLBA studies of jet-ISM interactions in both blazars and weak Seyfert jets have indicated that collimated outflows are indeed possible in dense environments. At present, the bulk of the evidence favors intermittent AGN accretion as the dominant factor in determining the evolutionary path of large numbers of AGN jets.

MOJAVE is a long term VLBA program to investigate the kinematics and polarization evolution of a complete sample of 133 active galactic nuclei selected on the basis of compact, relativistically beamed jet emission at 15 GHz. By fitting to the apparent distributions of superluminal speed and jet luminosity, we can constrain the Lorentz factor distribution and intrinsic luminosity function of the radio-selected blazar parent population. These low-energy peaked blazars formed a significant fraction of all EGRET detections, and should figure prominently in the GLAST source catalog. Using simple models, we investigate the predicted distribution of GLAST blazars in the gamma-ray/radio flux density plane, and describe an extension of the MOJAVE survey that will provide extensive parsec-scale jet information in complete regions of this plane. We find that if a population of intrinsically radio bright yet gamma-ray weak blazars exists, its signal will be largely wiped out by the large gamma-ray flux scatter associated with Doppler beaming.

We report on multi-frequency linear polarization monitoring of 15 active galactic nuclei containing highly relativistic jets with apparent speeds from $\sim$4 $c$ to $>40c$. The measurements were obtained at optical, 1 mm, and 3 mm wavelengths, and at 7 mm with the Very Long Baseline Array. The data show a wide range in degree of linear polarization among the sources, from $<$1% to $>$30%, and interday polarization variability in individual sources. The polarization properties suggest separation of the sample into three groups with low, intermediate, and high variability of polarization in the core at 7 mm : LVP, IVP, and HVP, respectively. The groups are partially associated with the common classification of active galactic nuclei as radio galaxies and quasars with low optical polarization (LVP), BL Lacertae objects (IVP), and highly optically polarized quasars (HVP). Our study investigates correlations between total flux, fractional polarization, and polarization position angle at the different wavelengths. We interpret the polarization properties of the sources in the sample through models in which weak shocks compress turbulent plasma in the jet. The differences in the orientation of sources with respect to the observer, jet kinematics, and abundance of thermal matter external to the jet near the core can account for the diversity in the polarization properties. The results provide strong evidence that the optical polarized emission originates in shocks, most likely situated between the 3 mm and 7 mm VLBI cores. They also support the idea that the 1 mm core lies at the edge of the transition zone between electromagnetically dominated and turbulent hydrodynamical sections of the jet.

The MOJAVE blazar sample consists of the 133 brightest, most compact radio-loud AGN in the northern sky, and is selected on the basis of VLBA 2 cm correlated flux density exceeding 1.5 Jy (2 Jy for declinations south of 0 degrees) at any epoch between 1994 and 2003. Since 1994 we have been gathering VLBA data on the sample to measure superluminal jet speeds and to better understand the parsec-scale kinematics of AGN jets. We have obtained 1.4 GHz VLA-A configuration data on 57 of these sources to investigate whether the extended luminosity of blazars is correlated with parsec-scale jet speed, and also to determine what other parsec-scale properties are related to extended morphology, such as optical emission line strength and gamma$ray emission. We present images and measurements of the kilo-parsec scale emission from the VLA data, which will be used in subsequent statistical studies of the MOJAVE sample.

We carried out a 5-pointing mosaic observation of TeV J2032+4130 at 1.4 and 4.8 GHz with the VLA in April of 2003. The analysis of the 4.8GHz data indicate weak wispy shell-like radio structure(s) which are at least partially non-thermal. The radio data is compatible with one or more young supernova remnants or perhaps the signature of large scale cluster shocks in this region induced by the violent action of the many massive stars in Cyg OB2.

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

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

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

We present total and polarized intensity images of 15 active galactic nuclei obtained with the Very Long Baseline Array at 7 mm at 17 epochs from 1998 March to 2001 April. At some epochs the images are accompanied by nearly simultaneous polarization measurements at 3 mm, 1.35/0.85 mm, and optical wavelengths. Here we analyze the 7 mm images to define the properties of the jets of two radio galaxies, five BL Lac objects, and eight quasars on angular scales $\gtrsim 0.1$ milliarcseconds. We determine the apparent velocities of 106 features in the jets; for many of the features we derive Doppler factors using a new method based on comparison of timescale of decline in flux density with the light-travel time across the emitting region. This allows us to estimate the Lorentz factors, intrinsic brightness temperatures, and viewing angles of 73 superluminal knots, as well as the opening angle of the jet for each source. We analyze the derived physical parameters of the jets. In nine sources we detect statistically meaningful deviations from ballistic motion, with the majority of components accelerating with distance from the core. In six sources we identify jet features with characteristics of trailing shocks that form behind the primary strong perturbations in jet simulations. The apparent speeds of these components increase with distance from the core suggesting an acceleration of the underlying jet.

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

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

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsec resolution using an orbiting 8 m telescope, HALCA, along with global arrays of Earth-based telescopes. Approximately 25% of the observing time is devoted to a survey of compact AGN which are stronger than 1 Jy at 5 GHz-the VSOP AGN Survey. This paper, the third in the series, presents the results from the analysis of the first 102 Survey sources. We present high resolution images and plots of visibility amplitude versus projected baseline length. In addition, model-fit parameters to the primary radio components are listed, and from these the angular size and brightness temperature for the radio cores are calculated. For those sources for which we were able to determine the source frame core brightness temperature, a significant fraction (53 out of 98) have a source frame core brightness temperature in excess of 10^12 K. The maximum source frame core brightness temperature we observed was 1.2 X 10^13 K. Explaining a brightness temperature this high requires an extreme amount of relativistic Doppler beaming. Since the maximum brightness temperature one is able to determine using only ground-based arrays is of the order of 10^12 K, our results confirm the necessity of using space VLBI to explore the extremely high brightness temperature regime.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsecond angular resolution using an orbiting 8-m telescope, HALCA and global arrays of Earth-based telescopes. Approximately 25% of the observing time has been devoted to a survey of compact AGN at 5 GHz which are stronger than 1 Jy -- the VSOP AGN Survey. This paper, the second in a series, describes the data calibration, source detection, self-calibration, imaging and modeling, and gives examples illustrating the problems specific to space VLBI. The VSOP Survey web-site which contains all results and calibrated data is described.

A sample of radio-loud active galactic nuclei (AGN) at 2cm is studied to test the isotropic distribution of radio sources in the sky and their uniform distribution in space. The sample is complete flux-density limits of 1.5Jy for positive declinations and 2Jy for declinations between 0 and -20 degrees. The active galactic nuclei sample comprises of 133 members. Application of the two-dimensional Kolmogorov-Smirnov test shows that there is no significant deviation from the isotropic distribution in the sky, while the generalised V/Vm test shows that the space distribution of AGN is not uniform at high confidence level (99.9%). This is indicative of a strong positive evolution of AGN with cosmic epoch implying that AGN (or jet activity phenomena) were more populous at high redshifts. It is shown that the evolution depends strongly on luminosity: low-luminosity QSOs show a strong positive evolution, while high-luminosity counterparts do not seem to show any evolution with cosmic epoch.

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

The idea that the radio jets in AGN contain material in relativistic motion is supported by many lines of observational evidence, including morphology, brightness temperature estimated with interferometers and with intrinsic variations, interstellar scintillations, X-rays, and superluminal motion. These are largely independent, and taken together make an irrefutable case for relativistic motion.

A revised sample of the 2 cm Very Long Baseline Array (VLBA) survey is studied to test the isotropic distribution of radio sources on the sky and their uniform distribution in space. The revised sample is complete to flux-density limits of 1.5 Jy for positive declinations and 2 Jy for declinations between 0 and -20 degrees. At present the active galactic nuclei sample comprises 122 members. Application of the two-dimensional Kolmogorov-Smirnov (K-S) test shows that there is no significant deviation from the homogenous distribution in the sky, while the V/Vmax test shows that the space distribution of active nuclei is not uniform at high confidence level (99.9%). This is indicative of a strong luminosity and/or density evolution implying that active nuclei (or jet activity phenomena) were more populous at high redshifts, z~2.

We discuss a new VLBA program to investigate structural and polarization changes in the brightest AGN jets in the northern sky. Our study represents a significant improvement over previous surveys in terms of image fidelity, size, and completeness, and will serve to characterize the kinematics of AGN jets and determine how these are related to other source properties. We discuss several preliminary results from our program, including the detection of parsec-scale circularly polarized jet emission, enhanced magnetic field ordering at the sites of apparent bends, and intrinsic differences in the jets of strong- and weak-lined AGN.

The intrinsic luminosity functions of extremely fast jets found in many active galaxies and gamma-ray bursts are difficult to measure since their apparent luminosities are strongly affected by relativistic beaming. Past studies have only provided analytical predictions for the beamed characteristics of populations in which all jets have the same Lorentz factor. However, the jets found in active galaxies are known to span a large range of speeds. Here we derive analytical expressions for the expected Doppler factor distributions and apparent (beamed) luminosity functions of randomly-oriented, two-sided jet populations which have bulk Lorentz factors distributed according to a simple power-law over a range [G1, G2]. We find that if a jet population has a uniform or inverted Lorentz factor distribution, its Doppler factor distribution and beamed luminosity function will be very similar to that of a jet population with Lorentz factors all equal to G2. In particular, the beamed and intrinsic luminosity functions will have the same slope at high luminosities. In the case of a steeply declining Lorentz factor distribution, the slopes at high luminosities will also be identical, but only up to an apparent luminosity equivalent to the upper cutoff of the intrinsic luminosity function. Also, at very high apparent luminosities, the slope will be proportional to the power law slope of the Lorentz factor distribution. At very low luminosities, the form of the apparent luminosity function is very sensitive to the lower cutoff and steepness of the Lorentz factor distribution. We discuss how it should be possible to recover useful information about the intrinsic luminosity functions of relativistic jets by studying flux-limited samples selected on the basis of beamed jet emission.

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

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

The study of compact active galactic nuclei (AGN) that possess convex radio spectra (the gigahertz-peaked spectrum radio sources) offers a unique opportunity to probe both the early evolutionary stages of relativistic AGN jets and their immediate nuclear environments. In this article I trace Ken Kellermann's early investigations of these sources, which played a major role in justifying the development of modern-day VLBI techniques. I describe how our understanding of these AGN has progressed since Kellermann's early discoveries, and discuss several ways in which the current classification scheme can be simplified to reflect intrinsic source characteristics, rather than observer-biased quantities. Finally, I discuss recent results from the VLBA 2 cm survey concerning the relativistic jet kinematics of the two-sided peaked-spectrum sources 4C +12.50 (PKS 1345+125) and OQ 208 (1404+286).

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

We present the results of a multi-epoch VLBA study of the powerful radio source 4C +12.50 (PKS 1345+125) at a wavelength of 2 cm. This compact radio source is associated with a hidden quasar whose host galaxy shows signs of a recent merger. It has been classified as a compact symmetric source (CSO) due to its small overall size (~220 pc) and twin-jet morphology, although it also has faint extended emission that may be a relic of previous activity. We report the detection of exceedingly high linear fractional polarization in isolated features of the southern jet (up to 60%), which is highly unusual for a CSO. Given the large amount of gas present in the host galaxy, we would expect significant Faraday depolarization across the whole source, unless the depolarizing gas is fairly clumpy. The southern jet also contains two features that are moving outward from the core at apparent speeds of v/c = 1.0 +/- 0.3 and 1.2 +/- 0.3. These represent the first positive detections of superluminal motion in a CSO, and taken together with the jet to counter-jet flux ratio, suggest an intrinsic flow speed of v/c = 0.84 +/- 0.12. The apparent ridgeline of the jet and counter-jet are consistent with a conical helix of wavelength 280 pc that is the result of Kelvin-Helmholtz instabilities driven by a slow precession of the jet nozzle. A fit to the data implies that the nozzle is precessing around a cone with half-angle 23 degrees, whose axis lies at an angle of 82 degrees to the line of sight. We suggest that the "S"-shaped jet morphologies commonly seen in recent AGN outflows such as 4C +12.50 may simply reflect the fact that their black hole spin axes are still precessing, and have not had sufficient time to align with their accretion disks.

We discuss the morphology and kinematics of five gigahertz-peaked spectrum (GPS) sources that have been observed with the VLBA. We find a wide range of observed properties including core-jet structure, superluminal motion, variability, extended structure, and polarization, all of which appear to deviate from commonly-accepted GPS paradigms (e.g., O'Dea 1998). We suggest that the observed low frequency cutoff in GPS sources may be primarily due to free-free absorption rather than synchrotron self-absorption.

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

We present results from the highest resolution polarization imaging survey of a complete sample of extragalactic radio sources carried out to date. Our sample consists of a statistically complete flat-spectrum subset of 32 AGNs from the Pearson-Readhead survey, which has been studied at a variety of wavelengths and resolutions, including space-VLBI. Our observations were made with the VLBA at 43 GHz, where the relatively higher resolution and weaker opacity effects have allowed us to probe magnetic field structures in the jets of these AGNs much closer to the central engine than in previous studies. At 43 GHz, the bulk of the total intensity and polarized emission in most flat-spectrum AGNs originates from an unresolved core component located at the extreme end of a faint jet. The luminosity of the core is positively correlated with the total source luminosity in soft x-rays, in the optical, and at 5 GHz. The most strongly polarized cores display electric vectors that are roughly aligned with the jet axis, which is consistent with a strong transverse shock that enhances the perpendicular component of the jet magnetic field. Sources with highly polarized cores also tend to have high optical polarizations and flatter overall radio spectra. Approximately half of the jets in our sample display apparently bent morphologies that are suggestive of streaming motions along a helical path. The straightest jets in the sample tend to display slower superluminal speeds than those that are significantly bent. We also confirm that intrinsic differences in the magnetic field properties of BL Lacertae objects and quasars previously seen on scales of tens of milliarcseconds are also present in regions much closer to the base of the jet.

We discuss the impact of special relativistic effects on the observed light curves and variability duty cycles of AGNs. We model the properties of AGN light curves at radio wavelengths using a simulated shot noise process in which the occurrence of major flaring events in a relativistic jet is governed by Poisson statistics. We show that flaring sources whose radiation is highly beamed toward us are able to reach very high flux levels, but will in fact spend most of their time in relatively low flaring states due to relativistic contraction of flare time scales in the observer frame. The fact that highly beamed AGNs do not return to a steady-state quiescent level between flares implies that their weakly beamed counterparts should have highly stable flux densities that result from a superposition of many long-term, low-amplitude flares. The ``apparent'' quiescent flux levels of these weakly beamed AGNs (identified in many unified models as radio galaxies) will be significantly higher than their ''true'' quiescent (i.e., non-flaring) levels. We use Monte Carlo simulations to investigate flux variability bias in the selection statistics of flat-spectrum AGN samples. In the case of the Caltech-Jodrell Flat-spectrum survey, the predicted orientation bias towards jets seen end-on is weakened if the parent population is variable, since the highly beamed sources have a stronger tendency to be found in low flaring states. This effect is small, however, since highly beamed sources are relatively rare, and their fluxes tend to be boosted sufficiently above the survey limit such that they are selected regardless of their flaring level. We find that for larger flat-spectrum AGN surveys with fainter flux cutoffs, variability should not be an appreciable source of selection bias.

We present images from a space-VLBI survey using the facilities of the VLBI Space Observatory Programme (VSOP), drawing our sample from the well-studied Pearson-Readhead survey of extragalactic radio sources. Our survey has taken advantage of long space-VLBI baselines and large arrays of ground antennas, such as the Very Long Baseline Array and European VLBI Network, to obtain high resolution images of 27 active galactic nuclei, and to measure the core brightness temperatures of these sources more accurately than is possible from the ground. A detailed analysis of the source properties is given in accompanying papers. We have also performed an extensive series of simulations to investigate the errors in VSOP images caused by the relatively large holes in the (u,v) plane when sources are observed near the orbit normal direction. We find that while the nominal dynamic range (defined as the ratio of map peak to off-source error) often exceeds 1000:1, the true dynamic range (map peak to on-source error) is only about 30:1 for relatively complex core-jet sources. For sources dominated by a strong point source, this value rises to approximately 100:1. We find the true dynamic range to be a relatively weak function of the difference in position angle (PA) between the jet PA and (u,v) coverage major axis PA. For low signal-to-noise regions typically located down the jet away from the core, large errors can occur, causing spurious features in VSOP images that should be interpreted with caution.

We have performed a multi-dimensional correlation analysis on the observed properties of a statistically complete core-selected sample of compact radio-loud active galactic nuclei, based on data from the VLBI Space Observing Programme (Paper I) and previously published studies. Our sample is drawn from the well-studied Pearson-Readhead (PR) survey, and is ideally suited for investigating the general effects of relativistic beaming in compact radio sources. In addition to confirming many previously known correlations, we have discovered several new trends that lend additional support to the beaming model. These trends suggest that the most highly beamed sources in core-selected samples tend to have a) high optical polarizations; b) large pc/kpc-scale jet misalignments; c) prominent VLBI core components; d) one-sided, core, or halo radio morphology on kiloparsec scales; e) narrow emission line equivalent widths; and f) a strong tendency for intraday variability at radio wavelengths. We have used higher resolution space and ground-based VLBI maps to confirm the bi-modality of the jet misalignment distribution for the PR survey, and find that the sources with aligned parsec- and kiloparsec-scale jets generally have arcsecond-scale radio emission on both sides of the core. The aligned sources also have broader emission line widths. We find evidence that the BL Lacertae objects in the PR survey are all highly beamed, and have very similar properties to the high-optically polarized quasars, with the exception of smaller redshifts. A cluster analysis on our data shows that after partialing out the effects of redshift, the luminosities of our sample objects in various wave bands are generally well-correlated with each other, but not with other source properties.

A significant fraction of compact radio-loud quasars display most of the characteristics of relativistically beamed, high-optical polarization blazars, yet are weakly polarized in the optical regime. We have used the VLBA at 22 and 43 GHz to look for differences in the parsec-scale magnetic field structures of 18 high- and low-optically polarized, compact radio-loud quasars (HPQs and LPRQs, respectively). We find a strong correlation between the polarization level of the unresolved parsec-scale radio core at 43 GHz and overall optical polarization of the source, which suggests a common (possibly co-spatial) origin for the emission at these two wavelengths. The magnetic fields of the polarized 43 GHz radio cores are aligned roughly transverse to the jet axis. Similar orientations are seen in the optical, suggesting that the polarized flux at both wavelengths is due to one or more strong transverse shocks located very close to the base of the jet. In LPRQs, these shocks appear to be weak near the core, and gradually increase in strength down the jet. The LPRQs in our sample tend to have less luminous radio cores than the HPQs, and jet components with magnetic fields predominantly parallel to the flow. The components in HPQ jets, on the other hand, tend to have transverse alignments. These differences cannot be accounted for by a simple model in which HPQs and LPRQs are the same type of object, seen at different angles to the line of sight. A more likely scenario is that LPRQs represent a quiescent phase of blazar activity, in which the inner jet flow does not undergo strong shocks.