We present the results from the timing and spectral study of Mrk\,421 based mainly on the \emph{Swift} data in the X-ray energy range obtained during the time interval 2015\,December--2018\,April. The most extreme X-ray flaring activity on the long-term, daily and intraday timescales was observed during the 2-month period which started in 2017\,December when the 0.3--10\,keV flux exceeded a level of 5$\times$10$^{-9}$erg\,cm$^{-2}$s$^{-1}$, recorded only twice previously. While the TeV-band and X-ray variabilities mostly were correlated, the source often varied in a complex manner in the MeV--GeV and radio--UV energy ranges, indicating that the multifrequency emission of Mrk\,421 could not be always generated in a single zone. The longer-term flares at X-rays and $\gamma$-rays showed a lognormal character, possibly indicating a variability imprint of the accretion disk onto the jet. A vast majority of the 0.3--10\,keV spectra were consistent with the log-parabolic model, showing relatively low spectral curvature and correlations between the different spectral parameters, predicted in the case of the first and second-order Fermi processes. The position of the synchrotron spectral energy distribution (SED) peak showed an extreme variability on diverse timescales between the energies $E_{\rm p}$$<$0.1\,keV and $E_{\rm p}$$>$15\,keV, with 15\% of the spectra peaking at hard X-rays and was related to the peak height as $S_{\rm p}$$\varpropto$$E^{\alpha}_{\rm p}$ with $\alpha$$\sim$0.6, which is expected for the transition from Kraichnan-type turbulence into the \textquotedblleft hard-sphere\textquotedblright~ one. The 0.3--300\,GeV spectra showed the features of the hadronic contribution, jet-star interaction and upscatter in the Klein-Nishina regime in different time intervals.

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $\gamma$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $\gamma$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $\gamma$-ray (VERITAS, MAGIC), high-energy (HE) $\gamma$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets\"ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($\delta \gtrsim 33$) and the size of the emission region ($ \delta^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

Synchrotron self-absorption in active galactic nuclei (AGN) jets manifests itself as a time delay between flares observed at high and low radio frequencies. It is also responsible for the observing frequency dependent change in size and position of the apparent base of the jet, aka the core shift effect, detected with very long baseline interferometry (VLBI). We measure the time delays and the core shifts in 11 radio-loud AGN to estimate the speed of their jets without relying on multi-epoch VLBI kinematics analysis. The 15$-$8 GHz total flux density time lags are obtained using Gaussian process regression, the core shift values are measured using VLBI observations and adopted from the literature. A strong correlation is found between the apparent core shift and the observed time delay. Our estimate of the jet speed is higher than the apparent speed of the fastest VLBI components by the median coefficient of 1.4. The coefficient ranges for individual sources from 0.5 to 20. We derive Doppler factors, Lorentz factors and viewing angles of the jets, as well as the corresponding de-projected distance from the jet base to the core. The results support evidence for acceleration of the jets with bulk motion Lorentz factor $\Gamma\propto R^{0.52\pm0.03}$ on de-projected scales $R$ of 0.5$-$500 parsecs.

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.

Clues to the physical conditions in radio cores of blazars come from measurements of brightness temperatures as well as effects produced by intrinsic opacity. We study the properties of the ultra compact blazar AO 0235+164 with RadioAstron ground-space radio interferometer, multi-frequency VLBA, EVN and single-dish radio observations. We employ visibility modeling and image stacking for deriving structure and kinematics of the source, and use Gaussian process regression to find the relative multi-band time delays of the flares. The multi-frequency core size and time lags support prevailing synchrotron self absorption. The intrinsic brightness temperature of the core derived from ground-based VLBI is close to the equipartition regime value. In the same time, there is evidence for ultra-compact features of the size of less than 10 $\mu$as in the source, which might be responsible for the extreme apparent brightness temperatures of up to $10^{14}$ K as measured by RadioAstron. In 2007--2016 the VLBI components in the source at 43 GHz are found predominantly in two directions, suggesting a bend of the outflow from southern to northern direction. The apparent opening angle of the jet seen in the stacked image at 43 GHz is two times wider than that at 15 GHz, indicating a collimation of the flow within the central 1.5 mas. We estimate the Lorentz factor $\Gamma = 14$, the Doppler factor $\delta=21$, and the viewing angle $\theta = 1.7^\circ$ of the apparent jet base, derive the gradients of magnetic field strength and electron density in the outflow, and the distance between jet apex and the core at each frequency.

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.

The University of Michigan 26-m paraboloid was dedicated to obtaining linear polarization and total flux density observations of blazars from the mid-1960s until June 2012 providing an unprecedented record tracking centimeter-band variability over decades at 14.5, 8.0, and 4.8 GHz for both targeted objects and members of flux-limited samples. In the mid-1970s through the mid-1980s, and during the last decade of the program, observations were additionally obtained of circular polarization for a small sample of radio-bright (S>5Jy), active sources. Key program results include evidence supporting class-dependent differences in the magnetic field geometry of BL Lac and QSO jets, identification of linear polarization changes temporally associated with flux outbursts supporting a shock-in-jet scenario, and determination of the spectral evolution of the Stokes V amplitude and polarity for testing proposed models. Recent radiative transfer modeling during large flares supports a jet scenario with a kinetically-dominated, relativistic flow at parsec scales with embedded turbulent magnetic fields and dynamically-weak ordered components which may be helical; the circular polarization observations are consistent with linear-to-circular mode conversion within this turbulent jet environment.

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous.

We present a pipeline that allows recovering reliable information for all four Stokes parameters with high accuracy. Its novelty relies on the treatment of the instrumental effects already prior to the computation of the Stokes parameters contrary to conventional methods, such as the M\"uller matrix one. The instrumental linear polarization is corrected across the whole telescope beam and significant Stokes $Q$ and $U$ can be recovered even when the recorded signals are severely corrupted. The accuracy we reach in terms of polarization degree is of the order of 0.1-0.2 %. The polarization angles are determined with an accuracy of almost 1$^{\circ}$. The presented methodology was applied to recover the linear and circular polarization of around 150 Active Galactic Nuclei. The sources were monitored from July 2010 to April 2016 with the Effelsberg 100-m telescope at 4.85 GHz and 8.35 GHz with a cadence of around 1.2 months. The polarized emission of the Moon was used to calibrate the polarization angle. Our analysis showed a small system-induced rotation of about 1$^{\circ}$ at both observing frequencies. Finally, we identify five sources with significant and stable linear polarization; three sources remain constantly linearly unpolarized over the period we examined; a total of 11 sources have stable circular polarization degree $m_\mathrm{c}$ and four of them with non-zero $m_\mathrm{c}$. We also identify eight sources that maintain a stable polarization angle over the examined period. All this is provided to the community for polarization observations reference. We finally show that our analysis method is conceptually different from the traditionally used ones and performs better than the M\"uller matrix method. Although it was developed for a system equipped with circularly polarized feeds it can easily be modified for systems with linearly polarized feeds as well.

We studied the pc-scale core shift effect using radio light curves for three blazars, S5 0716+714, 3C 279 and BL Lacertae, which were monitored at five frequencies ($\nu$) between 4.8 GHz and 36.8 GHz using the University of Michigan Radio Astronomical Observatory (UMRAO), the Crimean Astrophysical Observatory (CrAO), and Metsahovi Radio Observatory for over 40 years. Flares were Gaussian fitted to derive time delays between observed frequencies for each flare ($\Delta t$), peak amplitude ($A$), and their half width. Using $A \propto \nu^{\alpha}$ we infer $\alpha$ in the range $-$16.67 to 2.41 and using $\Delta t \propto \nu^{1/k_r}$, we infer $k_r \sim 1$, employed in the context of equipartition between magnetic and kinetic energy density for parameter estimation. From the estimated core position offset ($\Omega_{r \nu}$) and the core radius ($r_{\rm core}$), we infer that opacity model may not be valid in all cases. The mean magnetic field strength at 1 pc ($B_1$) and at the core ($B_{\rm core}$), are in agreement with previous estimates. We apply the magnetically arrested disk model to estimate black hole spins in the range $0.15-0.9$ for these blazars, indicating that the model is consistent with expected accretion mode in such sources. The power law shaped power spectral density has slopes $-$1.3 to $-$2.3 and is interpreted in terms of multiple shocks or magnetic instabilities.

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.

Since the mid-1980s the shock-in-jet model has been the preferred paradigm to explain radio-band flaring in blazar jets. We describe our radiative transfer model incorporating relativistically-propagating shocks, and illustrate how the 4.8, 8, and 14.5 GHz linear polarization and total flux density data from the University of Michigan monitoring program, in combination with the model, constrain jet flow conditions and shock attributes. Results from strong Fermi-era flares in 4 blazars with widely-ranging properties are presented. Additionally, to investigate jet evolution on decadal time scales we analyze 3 outbursts in OT 081 spanning nearly 3 decades and find intrinsic changes attributable to flow changes at a common spatial location, or, alternatively, to a change in the jet segment viewed. The model's success in reproducing these data supports a scenario in which relativistic shocks compress a plasma with an embedded passive, initially-turbulent magnetic field, with additional ordered magnetic field components, one of which may be helical.

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

We observed the prototype blazar, BL Lacertae, extensively in optical and radio bands during an active phase in the period 2010--2013 when the source showed several prominent outbursts. We searched for possible correlations and time lags between the optical and radio band flux variations using multifrequency data to learn about the mechanisms producing variability. During an active phase of BL Lacertae, we searched for possible correlations and time lags between multifrequency light curves of several optical and radio bands. We tried to estimate any possible variability timescales and inter-band lags in these bands. We performed optical observations in B, V, R and I bands from seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in Ukraine, Finland and USA. Significant cross-correlations between optical and radio bands are found in our observations with a delay of cm-fluxes with respect to optical ones of ~250 days. The optical and radio light curves do not show any significant timescales of variability. BL Lacertae showed many optical 'mini-flares' on short time-scales. Variations on longer term timescales are mildly chromatic with superposition of many strong optical outbursts. In radio bands, the amplitude of variability is frequency dependent. Flux variations at higher radio frequencies lead the lower frequencies by days or weeks. The optical variations are consistent with being dominated by a geometric scenario where a region of emitting plasma moves along a helical path in a relativistic jet. The frequency dependence of the variability amplitude supports an origin of the observed variations intrinsic to the source.

We consider the high radio frequency (15 GHz - 353 GHz) properties and variability of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing us to relate our results to the general population. We find that >6.0% of our parent sample (>15.1% if only cool-core clusters are considered) contain a radio-source at 150 GHz of at least 3mJy (~1x10^23 W/Hz at our median redshift of z~0.13). Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component (Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increasing to >8.5% if only cool-core clusters are considered. We see little evidence for strong variability at 15 GHz on short (week-month) time-scales although we see variations greater than 20% at 150 GHz over 6-month times-frames for 4 of the 23 sources with multi-epoch observations. Much more prevalent is long-term (year-decade time-scale) variability, with average annual amplitude variations greater than 1% at 15 GHz being commonplace. There is a weak trend towards higher variability as the peak of the GPS-like component occurs at higher frequency. We demonstrate the complexity that is seen in the radio spectra of BCGs and discuss the potentially significant implications of these high-peaking components for Sunyaev-Zel'dovich cluster searches.

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

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto \nu^\alpha$, $\alpha$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $\Omega_{r\nu} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $\nu$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $\alpha$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

PKS 0521-36 is an Active Galactic Nucleus (AGN) with uncertain classification. We investigate the properties of this source from radio to gamma rays. The broad emission lines in the optical and UV bands and steep radio spectrum indicate a possible classification as an intermediate object between broad-line radio galaxies (BLRG) and steep spectrum radio quasars (SSRQ). On pc-scales PKS 0521-36 shows a knotty structure similar to misaligned AGN. The core dominance and the gamma-ray properties are similar to those estimated for other SSRQ and BLRG detected in gamma rays, suggesting an intermediate viewing angle with respect to the observer. In this context the flaring activity detected from this source by Fermi-LAT between 2010 June and 2012 February is very intriguing. We discuss the gamma-ray emission of this source in the framework of the structured jet scenario, comparing the spectral energy distribution (SED) of the flaring state in 2010 June with that of a low state. We present three alternative models corresponding to three different choices of the viewing angles theta_v = 6, 15, and 20 deg. We obtain a good fit for the the first two cases, but the SED obtained with theta_v = 15 deg if observed at a small angle does not resemble that of a typical blazar since the synchrotron emission should dominate by a large factor (about 100) the inverse Compton component. This suggests that a viewing angle between 6 and 15 deg is preferred, with the rapid variability observed during gamma-ray flares favouring a smaller angle. However, we cannot rule out that PKS 0521-36 is the misaligned counterpart of a synchrotron-dominated blazar.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

As part of a program to identify the physical conditions in the jets of gamma-ray-flaring blazars detected by Fermi, including the role of shocks in the production of high-energy flaring, we obtained 4 years of 3-frequency, centimeter-band total flux density and linear polarization monitoring observations of the radio-bright blazar S5 0716+714 with the University of Michigan 26-m paraboloid. Light curves constructed from these data exhibit a series of rapid, high-amplitude, centimeter-band total flux density outbursts, and changes in the linear polarization consistent with the passage of shocks during the gamma-ray flaring. The observed spectral evolution of the radio-band flares, in combination with radiative transfer simulations incorporating propagating shocks, was used to constrain the shock and jet flow conditions in the parsec-scale regions of the jet. Eight forward-moving, transverse shocks with unusually-strong shock compression factors, a very fast Lorentz factor of the shocks of 77, a bulk Lorentz factor of the flow of 20, a viewing angle of 12 degrees, and an intrinsic opening angle of the radio jet of 5.2 degrees were identified.

We analyze the shock-in-jet models for the gamma-ray flaring blazars 0420-014, OJ 287, and 1156+295 presented in Aller et al. (2014, Paper I), quantifying how well the modeling constrains internal properties of the flow (low energy spectral cutoff, partition between random and ordered magnetic field), the flow dynamics (quiescent flow speed and orientation), and the number and strength of the shocks responsible for radio-band flaring. We conclude that well-sampled, multifrequency polarized flux light curves are crucial for defining source properties. We argue for few, if any, low energy particles in these flows, suggesting no entrainment and efficient energization of jet material, and for approximate energy equipartition between the random and ordered magnetic field components, suggesting that ordered field is built by non-trivial dynamo action from the random component, or that the latter arises from a jet instability that preserves the larger-scale, ordered flow. We present evidence that the difference between orphan radio-band (no gamma-ray counterpart) and non-orphan flares is due to more complex shock interactions in the latter case.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

To investigate parsec-scale jet flow conditions during GeV gamma-ray flares detected by the Fermi Large Angle Telescope, we obtained centimeter-band total flux density and linear polarization monitoring observations from 2009.5 through 2012.5 with the 26-meter Michigan radio telescope for a sample of core-dominated blazars. We use these data to constrain radiative transfer simulations incorporating propagating shocks oriented at an arbitrary angle to the flow direction in order to set limits on the jet flow and shock parameters during flares temporally associated with gamma-ray flares in 0420-014, OJ 287, and 1156+295; these AGN exhibited the expected signature of shocks in the linear polarization data. Both the number of shocks comprising an individual radio outburst (3-4) and the range of the compression ratios of the individual shocks (0.5-0.8) are similar in all three sources; the shocks are found to be forward-moving with respect to the flow. While simulations incorporating transverse shocks provide good fits for 0420-014 and 1156+295, oblique shocks are required for modeling the OJ~287 outburst, and an unusually-low value of the low energy cutoff of the radiating particles' energy distribution is also identified. Our derived viewing angles and shock speeds are consistent with independent VLBA results. While a random component dominates the jet magnetic field, as evidenced by the low fractional linear polarization, to reproduce the observed spectral character requires that a significant fraction of the magnetic field energy is in an ordered axial component. Both straight and low pitch angle helical field lines are viable scenarios.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

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)

The broad-line radio galaxy 3C 111 exhibited a major flux density outburst in 2007. Here, we present imaging and preliminary kinematic results of the jet, based on three millimetre-VLBI observations (86 GHz) using the Global Millimeter VLBI Array (GMVA) covering one year just after the radio flare. The GMVA data allow us to study this outburst with unprecedented image fidelity at highest (sub-parsec) resolution. On these scales, the outburst is resolved into a complex series of plasma components forming an intriguing bent structure. Within 1 mas from the jet base, ejections vary in position angle and components move with an apparent velocity of ~3.7c, significantly slower than the maximum velocity observed with cm-VLBI on scales beyond 1 mas.

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.

As part of a program to investigate jet flow conditions during GeV gamma-ray flares detected by Fermi, we are using UMRAO multi-frequency, centimeter-band total flux density and linear polarization monitoring observations to constrain radiative transfer models incorporating propagating shocks orientated at an arbitrary angle to the flow direction. We describe the characteristics of the model, illustrate how the data are used to constrain the models, and present results for three program sources with diverse characteristics: PKS 0420-01, OJ 287, and 1156+295. The modeling of the observed spectral behavior yields information on the sense, strength and orientation of the shocks producing the radio-band flaring; on the energy distribution of the radiating particles; and on the observer's viewing angle with respect to the jet independent of VLBI data. We present evidence that, while a random component dominates the jet magnetic field, a distinguishing feature of those radio events with an associated gamma-ray flare is the presence of a weak but non-negligible ordered magnetic field component along the jet axis.

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

We report 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 present a study of the central engine in the broad-line radio galaxy 3C120 using a multi-epoch analysis of a deep XMM-Newton observation and two deep Suzaku pointings (in 2012). In order to place our spectral data into the context of the disk-disruption/jet-ejection cycles displayed by this object, we monitor the source in the UV/X-ray bands, and in the radio band. We find three statistically acceptable spectral models, a disk-reflection model, a jet-model and a jet+disk model. Despite being good descriptions of the data, the disk-reflection model violates the radio constraints on the inclination, and the jet-model has a fine-tuning problem, requiring a jet contribution exceeding that expected. Thus, we argue for a composite jet+disk model. Within the context of this model, we verify the basic predictions of the jet-cycle paradigm, finding a truncated/refilling disk during the Suzaku observations and a complete disk extending down to the innermost stable circular orbit (ISCO) during the XMM-Newton observation. The idea of a refilling disk is further supported by the detection of the ejection of a new jet knot approximately one month after the Suzaku pointings. We also discover a step-like event in one of the Suzaku pointings in which the soft band lags the hard band. We suggest that we are witnessing the propagation of a disturbance from the disk into the jet on a timescale set by the magnetic field.

We carried out intra-day variability (IDV) observations from August 2005 to January 2010 with the Urumqi 25m radio telescope for a dozen IDV sources including the quasar 0917+624. This target exhibited pronounced centimeter-band, intra-day variability during the 1980s--1990s, but its strong IDV phase ceased in 2000. The source showed no IDV in the majority of the Urumqi observing sessions, although weak IDV activity was detected in some. Multifrequency UMRAO data for 0917+624 show that the spectral index is steeper during the weak and non-IDV phases than during the strong IDV phase, supporting the idea that the size of the scintillating component may be enlarged in the weak/non IDV phases.

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

We have used single-dish centimeter-band, multi-frequency linear polarization and total flux density monitoring observations obtained by UMRAO at 14.5, 8, and 4.8 GHz during gamma-ray flaring to constrain propagating shock-in-jet models; this procedure allows us to probe jet conditions at or near the presumed parsec-scale gamma-ray emission site within the relativistic jet. Results are presented for the blazar 0420-014 during the gamma-ray flare detected by the LAT which peaked in late January 2010. Our work identifies the passage of three forward-moving transverse shocks during the radio-band flaring, a shock Lorentz factor of 12, and an observer's viewing angle of 2 degrees and sets limits on the energy distribution of the radiating particles.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

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.

We discuss the parsec-scale structural variability of the extragalactic jet 3C 111 related to a major radio flux density outburst in 2007. The data analyzed were taken within the scope of the MOJAVE, UMRAO, and F-GAMMA programs, which monitor a large sample of the radio brightest compact extragalactic jets with the VLBA, the University of Michigan 26 m, the Effelsberg 100 m, and the IRAM 30 m radio telescopes. The analysis of the VLBA data is performed by fitting Gaussian model components in the visibility domain. We associate the ejection of bright features in the radio jet with a major flux-density outburst in 2007. The evolution of these features suggests the formation of a leading component and multiple trailing components.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

We present centimeter-band total flux density and linear polarization light curves illustrating the signature of shocks during radio band outbursts associated in time with gamma-ray flares detected by the Fermi LAT. The general characteristics of the spectral evolution during these events is well-explained by new radiative transfer simulations incorporating propagating oblique shocks and assuming an initially turbulent magnetic field. This finding supports the idea that oblique shocks in the jet are a viable explanation for activity from the radio to the gamma-ray band in at least some gamma-ray flares.

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.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

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

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

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

The `shock in jet' model for cm-waveband blazar variability is revisited, allowing for arbitrary shock orientation with respect to the jet flow direction, and both random and ordered magnetic field. It is shown that oblique shocks can explain events with swings in polarization position angle much less than the 90 deg. associated with transverse structures, while retaining the general characteristics of outbursts, including spectral behavior and level of peak percentage polarization. Models dominated by a force-free, minimum energy magnetic field configuration (essentially helical) display a shallow rise in percentage polarization and frequency dependent swing in polarization position angle not in agreement with the results of single-dish monitoring observations, implying that the field is predominantly random in the quiescent state. Outbursts well-explained by the `shock in jet' model are present during gamma-ray flaring in several sources, supporting the idea that shock events are responsible for activity from the radio to gamma-ray bands.

We present the results of extensive multi-frequency monitoring of the radio galaxy 3C 111 between 2004 and 2010 at X-ray (2.4--10 keV), optical (R band), and radio (14.5, 37, and 230 GHz) wave bands, as well as multi-epoch imaging with the Very Long Baseline Array (VLBA) at 43 GHz. Over the six years of observation, significant dips in the X-ray light curve are followed by ejections of bright superluminal knots in the VLBA images. This shows a clear connection between the radiative state near the black hole, where the X-rays are produced, and events in the jet. The X-ray continuum flux and Fe line intensity are strongly correlated, with a time lag shorter than 90 days and consistent with zero. This implies that the Fe line is generated within 90 light-days of the source of the X-ray continuum. The power spectral density function of X-ray variations contains a break, with steeper slope at shorter timescales. The break timescale of 13 (+12,-6) days is commensurate with scaling according to the mass of the central black hole based on observations of Seyfert galaxies and black hole X-ray binaries (BHXRBs). The data are consistent with the standard paradigm, in which the X-rays are predominantly produced by inverse Compton scattering of thermal optical/UV seed photons from the accretion disk by a distribution of hot electrons --- the corona --- situated near the disk. Most of the optical emission is generated in the accretion disk due to reprocessing of the X-ray emission. The relationships that we have uncovered between the accretion disk and the jet in 3C 111, as well as in the FR I radio galaxy 3C 120 in a previous paper, support the paradigm that active galactic nuclei and Galactic BHXRBs are fundamentally similar, with characteristic time and size scales proportional to the mass of the central black hole

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

The origin of the high-energy emission of blazars is still a matter of debate. To investigate the emission mechanism of extragalactic outflows and to pin down the location of the emission, we have constructed a broadband spectral energy distribution (SED) database covering from the radio to the gamma-ray band for the complete MOJAVE sample, which consists of 135 relativistically beamed AGN with well-studied parsec-scale jets. Typically, the broadband SEDs of blazars shows a double-humped profile. It is believed that the lower-energy hump is due to synchrotron emission from the radio jet, and the higher-energy hump is generated by i) inverse-Compton upscattered seed photons (leptonic), ii) proton-induced shower (hadronic). Combining the results of high-resolution VLBI observations and the gamma-ray properties of the MOJAVE sources, we attempt to reveal the origin of the high-energy emission in relativistic jets, and search for correlations between VLBI and high-energy properties.

The properties of blazar variability in the radio band are studied using the unique combination of temporal resolution from single dish monitoring and spatial resolution from VLBA imaging; such measurements, now available in all four Stokes parameters, together with theoretical simulations, identify the origin of radio band variability and probe the characteristics of the radio jet where the broadband blazar emission originates. Outbursts in total flux density and linear polarization in the optical-to-radio bands are attributed to shocks propagating within the jet spine, in part based on limited modeling invoking transverse shocks; new radiative transfer simulations allowing for shocks at arbitrary angle to the flow direction confirm this picture by reproducing the observed centimeter-band variations observed more generally, and are of current interest since these shocks may play a role in the gamma-ray flaring detected by Fermi. Recent UMRAO multifrequency Stokes V studies of bright blazars identify the spectral variability properties of circular polarization for the first time and demonstrate that polarity flips are relatively common. All-Stokes data are consistent with the production of circular polarization by linear-to-circular mode conversion in a region that is at least partially self-absorbed. Detailed analysis of single-epoch, multifrequency, all-Stokes VLBA observations of 3C 279 support this physical picture and are best explained by emission from an electron-proton plasma.

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

The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) monitored BL Lacertae in 2008-2009 at radio, near-IR, and optical frequencies. During this period, high-energy observations were performed by XMM-Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to interpret the source broad-band variability. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, as well as moderate flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the line of sight, producing an increase of the Doppler beaming factor. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability.

The quasar B0605-085 (OH 010) shows a hint for probable periodical variability in the radio total flux-density light curves. We study the possible periodicity of B0605-085 in the total flux-density, spectra and opacity changes in order to compare it with jet kinematics on parsec scales. We have analyzed archival total flux-density variability at ten frequencies (408 MHz, 4.8 GHz, 6.7 GHz, 8 GHz, 10.7 GHz, 14.5 GHz, 22 GHz, 37 GHz, 90 GHz, and 230 GHz) together with the archival high-resolution very long baseline interferometry data at 15 GHz from the MOJAVE monitoring campaign. Using the Fourier transform and discrete autocorrelation methods we have searched for periods in the total flux-density light curves. In addition, spectral evolution and changes of the opacity have been analyzed. We found a period in multi-frequency total flux-density light curves of 7.9+-0.5 yrs. Moreover, a quasi-stationary jet component C1 follows a prominent helical path on a similar time scale of 8 years. We have also found that the average instantaneous speeds of the jet components show a clear helical pattern along the jet with a characteristic scale of 3 mas. Taking into account average speeds of jet components, this scale corresponds to a time scale of about 7.7 years. Jet precession can explain the helical path of the quasi-stationary jet component C1 and the periodical modulation of the total flux-density light curves. We have fitted a precession model to the trajectory of the jet component C1, with a viewing angle phi=2.6+-2.2 degrees, aperture angle of the precession cone Omega=23.9+-1.9 degrees and fixed precession period (in the observers frame) P = 7.9 yrs.

In the 1990s a comparison of sparse EGRET measurements with single-dish flux density monitoring from the Metsahovi and UMRAO programs established a temporal connection between the onset of flaring at radio band and the occurrence of gamma-ray activity. Correlations between the emergence of new VLBI components from the core, flares in linearly polarized radio flux, and gamma-ray activity in bright EGRET-detected blazars supported a picture in which the gamma-ray and the radio band emission arises in the same shocked region of the jet, with the high energy emission produced via inverse Compton scattering by the synchrotron-emitting electrons in the jet. Quantitative tests of this scenario, however, were hampered by insufficient temporal sampling of the data and the simple nature of the models adopted. The extensive data from Fermi coupled with the wealth of well-sampled radio band data from old as well as new programs such as the F-GAMMA project now permit statistical studies for large numbers of sources, including weak HBLs, and detailed analyses of individual highly-active class members. I summarize progress in understanding the origin of the gamma-ray emission using these new measurements. I focus on three areas: attempts to isolate the physical site of the high energy emission using time delay information; investigation of the emission process using the characteristics of the variability; and quantitative tests of the shock model picture using high-time-sampled multifrequency linear polarization data, VLBP imaging, and new models of propagating oblique relativistic shocks incorporating detailed radiative transfer calculations.

We are constructing the broadband SED catalog of the MOJAVE sample from the radio to the gamma-ray band using MOJAVE, Swift UVOT/XRT/BAT, and Fermi/LAT data, in order to understand the emission mechanism of extragalactic outflows and to investigate the site of high-energy emission in AGN. Since the launch of Fermi gamma-ray Space Telescope in August 2008, two thirds of the MOJAVE sources have been detected by Fermi/LAT. Combining the results of high-resolution VLBI, X-ray, and gamma-ray observations of the jet-dominated AGN sample, we want to pin down the origin of high-energy emission in relativistic jets. Here we present our overall project and preliminary results for 6 selected sources.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

Context. We present evidence that parsec-scale jets in BL Lac objects may be significantly distinct in kinematics from their counterparts in quasars. We argued this previously for the BL lac sources 1803+784 and 0716+714, report here a similar pattern for another well-known BL Lac object, PKS 0735+178, whose nuclear jet is found to exhibit kinematics atypical of quasars. Aims. A detailed study of the jet components' motion reveals that the standard AGN paradigm of apparent superluminal motion does not always describe the kinematics in BL Lac objects. We study 0735+178 here to augment and improve the understanding of the peculiar motions in the jets of BL Lac objects as a class. Methods. We analyzed 15 GHz VLBA (Very Long Baseline Array) observations (2cm/MOJAVE survey) performed at 23 epochs between 1995.27 and 2008.91. Results. We found a drastic structural mode change in the VLBI jet of 0735+178, between 2000.4 and 2001.8 when its twice sharply bent trajectory turned into a linear shape.We further found that this jet had undergone a similar transition sometime between December 1981 and June 1983. A mode change, occurring in the reverse direction (between mid-1992 and mid-1995) has already been reported in the literature. These structural mode changes are found to be reflected in changed kinematical behavior of the nuclear jet, manifested as an apparent superluminal motion and stationarity of the radio knots. In addition, we found the individual mode changes to correlate in time with the maxima in the optical light curve. The last two transitions occurred before a (modest) radio flare. The behavior of this pc-scale jet appears to favor a scenario involving non-ballistic motions of the radio knots, produced by the precession of a continuous jet within the ambient medium.

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 present results from monitoring the multi-waveband flux, linear polarization, and parsec-scale structure of the quasar PKS 1510-089, concentrating on eight major gamma-ray flares that occurred during the interval 2009.0-2009.5. The gamma-ray peaks were essentially simultaneous with maxima at optical wavelengths, although the flux ratio of the two wavebands varied by an order of magnitude. The optical polarization vector rotated by 720 degrees during a 5-day period encompassing six of these flares. This culminated in a very bright, roughly 1 day, optical and gamma-ray flare as a bright knot of emission passed through the highest-intensity, stationary feature (the "core") seen in 43 GHz Very Long Baseline Array images. The knot continued to propagate down the jet at an apparent speed of 22c and emit strongly at gamma-ray energies as a months-long X-ray/radio outburst intensified. We interpret these events as the result of the knot following a spiral path through a mainly toroidal magnetic field pattern in the acceleration and collimation zone of the jet, after which it passes through a standing shock in the 43 GHz core and then continues downstream. In this picture, the rapid gamma-ray flares result from scattering of infrared seed photons from a relatively slow sheath of the jet as well as from optical synchrotron radiation in the faster spine. The 2006-2009.7 radio and X-ray flux variations are correlated at very high significance; we conclude that the X-rays are mainly from inverse Compton scattering of infrared seed photons by 20-40 MeV electrons.

For the BL Lac object S5 0716+714 conflicting apparent velocities of jet component motion have been reported in the literature. This object is an intra-day variable source and suited to investigate a possible correlation between kinematic properties and flux-density variability on different timescales. We (re-)analyze 50 VLBI observations obtained with the VLBA at 5 different frequencies (5 to 43 GHz) between 1992.73 and 2006.32. We analyze the jet component motion in detail taking care not only to account for motion in the radial but also in the orthogonal direction. We study the evolution of the jet ridge line and search for correlations between radio band light curves and the kinematic properties of the jet components. We present an alternative kinematic scenario for jet component motion in S5 0716+714 with apparent stationarity of jet components (with regard to their core separation) with time. Jet components, however, do seem to move significantly non-radially with regard to their position angle and in a direction perpendicular to the major axis of the jet. We discuss a possible correlation between the long-term radio flux-density variability and the kinematics. We conclude that a geometric contribution to the origin of the long-term variability might not be negligible. Subluminal motion has been reported for most of the TeV blazars. Our analysis also confirms this for S5 0716+714. This result increases the number of TeV blazars showing apparent subluminal motion to 7.

We present a kinematic analysis of jet component motion in the VLBI jet of the BL Lac object S5 1803+784, which does not reveal long-term outward motion for most of the components. Understanding the complex kinematic phenomena can possibly provide insights into the differences between quasars and BL Lac objects. The blazar S5 1803+784 has been studied with VLBI at $\nu$ =1.6, 2.3, 5, 8.4, and 15 GHz between 1993.88 and 2005.68 in 26 observing runs. We (re)analyzed the data and present Gaussian model-fits. We collected the already published kinematic information for this source from the literature and re-identified the components according to the new scenario presented in this paper. Altogether, 94 epochs of observations have been investigated. A careful study of the long-term kinematics reveals a new picture for component motion in S5 1803+784. In contrast to previously discussed motion scenarios, we find that the jet structure within 12 mas of the core can most easily be described by the coexistence of several bright jet features that remain on the long-term at roughly constant core separations (in addition to the already known {\it stationary} jet component $\sim$ 1.4 mas) and one faint component moving with an apparent superluminal speed ($\sim$ 19c, based on 3 epochs). While most of the components maintain long-term roughly constant distances from the core, we observe significant, smooth changes in their position angles. We report on an evolution of the whole jet ridge line with time over the almost 12 years of observations. The width of the jet changes periodically with a period of $\sim$ 8 to 9 years. We find a correlation between changes in the position angle and maxima in the total flux-density. We present evidence for a geometric origin of the phenomena and discuss possible models.

We describe and present initial results from a Fermi cycle 2 program designed to monitor the behavior of the centimeter-band linear polarization and total flux density emitted by gamma-ray-bright blazars during flaring. The goal of the program is to identify changes in the magnetic field structure in the radio jet associated with gamma-ray flaring and ultimately to test whether gamma-ray flaring is associated with the onset of shocks in the radio jet. Light curves illustrating radio band variability patterns are shown for sample program sources.

We present analysis of both the short-term optical and long-term multiwavelength variability of CTA 102. In 2004, this object was observed in an intense optical flaring state. Extensive R-band microvariability observations were carried out during this high state. In 2005, we obtained several weeks of contemporaneous radio, optical, and X-ray observations of CTA 102. These observations recorded distinct flaring activity in all three wavebands. Subsequent analysis revealed that this object may appear redder when in a brighter optical state, and that the X-ray, optical, and radio activity do not appear to be correlated. The shape of the observed spectral energy distributions suggests that both synchrotron-related and external inverse Compton processes may contribute to the X-ray emission. Our results are also compared with other results on this object and archival microvariability observations. It appears that more rapid, dramatic microvariability events occur when CTA 102 is in an elevated optical flux state.

We report results from an intensive multiwavelength campaign on the intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a strong outburst of very high energy gamma-ray emission in June 2008. The very high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with a flux F(>200 GeV) = (5.7+-0.6)x10^-11 cm-2s-1, about three times brighter than during the discovery of gamma-ray emission from W Com by VERITAS in 2008 March. The initial detection of this flare by VERITAS at energies above 200 GeV was followed by observations in high energy gamma-rays (AGILE, E>100 MeV), and X-rays (Swift and XMM-Newton), and at UV, and ground-based optical and radio monitoring through the GASP-WEBT consortium and other observatories. Here we describe the multiwavelength data and derive the spectral energy distribution (SED) of the source from contemporaneous data taken throughout the flare.

Radio-bright BL Lacertae objects (BLOs) are typically variable and exhibit prominent flaring. We use a sample of 24 BLOs to get a clear idea of their flaring behavior and to find possible commonalities in their variability patterns. Our goal was to compare the results given by computational time scales and the observed variability parameters determined directly from the flux curves. Also, we wanted to find out if the BLO flares adhere to the generalized shock model. We use long-term monitoring data from 4.8, 8, 14.5, 22, 37, 90 and 230 GHz. The structure function, discrete correlation function and Lomb-Scargle periodogram time scales, calculated in a previous study, are analyzed in more detail. We determine flare durations, rise and decay times, absolute and relative peak fluxes from the monitoring data. We find that BLOs demonstrate a wide range of variability behavior. BLOs include sources with fast and strong variability, such as OJ 287, PKS 1749+096 and BL Lac, but also sources with more rolling fluctuations like PKS 0735+178. The most extreme flares can last for up to 13 years or have peak fluxes of approximately 12 Jy in the observer's frame. When the Doppler boosting effect is taken into account, the peak flux of a flare does not depend on the duration of the flare. A rough analysis of the time lags and peak flux evolution indicates that BLO flares in the mm - cm wavelengths are high-peaking, i.e., are in the adiabatic stage. Thus, the results concur with the generalized shock model.

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

We present the results of extensive multi-frequency monitoring of the radio galaxy 3C 120 between 2002 and 2007 at X-ray, optical, and radio wave bands, as well as imaging with the Very Long Baseline Array (VLBA). Over the 5 yr of observation, significant dips in the X-ray light curve are followed by ejections of bright superluminal knots in the VLBA images. Consistent with this, the X-ray flux and 37 GHz flux are anti-correlated with X-ray leading the radio variations. This implies that, in this radio galaxy, the radiative state of accretion disk plus corona system, where the X-rays are produced, has a direct effect on the events in the jet, where the radio emission originates. The X-ray power spectral density of 3C 120 shows a break, with steeper slope at shorter timescale and the break timescale is commensurate with the mass of the central black hole based on observations of Seyfert galaxies and black hole X-ray binaries. These findings provide support for the paradigm that black hole X-ray binaries and active galactic nuclei are fundamentally similar systems, with characteristic time and size scales linearly proportional to the mass of the central black hole. The X-ray and optical variations are strongly correlated in 3C 120, which implies that the optical emission in this object arises from the same general region as the X-rays, i.e., in the accretion disk-corona system. We numerically model multi-wavelength light curves of 3C 120 from such a system with the optical-UV emission produced in the disk and the X-rays generated by scattering of thermal photons by hot electrons in the corona. From the comparison of the temporal properties of the model light curves to that of the observed variability, we constrain the physical size of the corona and the distances of the emitting regions from the central BH.

We present the result of our extensive intranight optical monitoring of the well known low-energy peaked BL Lac (LBL) object PKS 0735+178. This long-term follow-up consists of $R$-band monitoring for a minimum duration of $\sim 4$ hours, on 17 nights spanning 11 years (1998-2008). Using the CCD as an N-star photometer, a detection limit of around 1% was attained for the intra-night optical variability (INOV). Remarkably, an INOV amplitude of $\geq 3%$ on hour-like time scale was not observed on any of the 17 nights, even though the likelihood of a typical LBL showing such INOV levels in a single session of $\ga 4$ hours duration is known to be high ($\sim50%$). Our observations have thus established a peculiar long-term INOV quiescence of this radio-selected BL Lac object. Moreover, the access to unpublished optical monitoring data of similarly high sensitivity, acquired in another programme, has allowed us to confirm the same anomalous INOV quiescence of this LBL all the way back to 1989, the epoch of its historically largest radio outburst. Here, we present observational evidence revealing the very unusual INOV behaviour of this classical BL Lac object and discuss this briefly in the context of its other known exceptional properties.

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.

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 present an analysis of the 101 ks, 2007 Suzaku spectrum of the LINER galaxy NGC 1052. The 0.5-10 keV continuum is well-modeled by a power-law modified by Galactic and intrinsic absorption, and it exhibits a soft, thermal emission component below 1 keV. Both a narrow core and a broader component of Fe K emission centered at 6.4 keV are robustly detected. While the narrow line is consistent with an origin in material distant from the black hole, the broad line is best fit empirically by a model that describes fluorescent emission from the inner accretion disk around a rapidly rotating black hole. We find no evidence in this observation for Comptonized reflection of the hard X-ray source by the disk above 10 keV, however, which casts doubt on the hypothesis that the broad iron line originates in the inner regions of a standard accretion disk. We explore other possible scenarios for producing this spectral feature and conclude that the high equivalent width (EW ~ 185 keV) and full-width-half-maximum velocity of the broad iron line (v ~ 0.37c) necessitate an origin within d ~ 8 gravitational radii of the hard X-ray source. Based on the confirmed presence of a strong radio jet in this galaxy nucleus, the broad iron line may be produced in dense plasma near the base of the jet, implying that emission mechanisms in the centralmost portions of active galactic nuclei are more complex than previously thought.

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)

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

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

The radio loud galaxy NGC 1052 is being studied in an intensive multi-band campaign including X-ray brigthness monitoring and spectroscopic observations, single-dish radio brightness monitoring at centimetre wavelengths, and a high-frequency very-long-baseline interferometry monitoring program. Here we present a progress report on our studies from this program. The final goal of our observations is to relate the findings from the high-resolution radio images with the observed variations in the X-ray regime, to address the accretion processes and their relationship with the radio jet activity.

We report the detection by the AGILE (Astro-rivelatore Gamma a Immagini LEggero) satellite of an intense gamma-ray flare from the source AGL J1511-0909, associated with the powerful quasar PKS 1510-089, during ten days of observations from 23 August to 1 September 2007. During the observation period, the source was in optical decrease following a flaring event monitored by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). The simultaneous gamma-ray, optical, and radio coverage allows us to study the spectral energy distribution and the theoretical models based on the synchrotron and inverse Compton (IC) emission mechanisms. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector and the hard X-ray imager Super-AGILE sensitive in the 30 MeV - 50 GeV and 18 - 60 keV bands, respectively. Between 23 and 27 August 2007, AGILE detected gamma-ray emission from PKS 1510-089 when this source was located about 50 degrees off-axis, with an average flux of (270 +/- 65) x 10^{-8} photons cm^{-2} s^{-1} for photon energy above 100 MeV. In the following period, 28 August - 1 September, after a satellite re-pointing, AGILE detected the source at about 35 degrees off-axis, with an average flux (E > 100 MeV) of (195 +/- 30) x 10^{-8} photons cm^{-2} s^{-1}. No emission was detected by Super-AGILE, with a 3-sigma upper limit of 45 mCrab in 200 ksec. The spectral energy distribution is modelled with a homogeneous one-zone synchrotron self Compton (SSC) emission plus contributions by external photons: the SSC emission contributes primarily to the X-ray band, whereas the contribution of the IC from the external disc and the broad line region match the hard gamma-ray spectrum observed.

A large number of AGN have been monitored for nearly 30 years at 22, 37 and 87 GHz in Mets\"ahovi Radio Observatory. These data were combined with lower frequency 4.8, 8.0 and 14.5 GHz data from the University of Michigan Radio Astronomy Observatory, higher frequency data at 90 and 230 GHz from SEST, and supplementary higher frequency data from the literature to study the long-term variability of a large sample of AGN. Both the characteristics of individual flares from visual inspection and statistically-determined variability timescales as a function of frequency and optical class type were determined. Based on past behaviour, predictions of sources expected to exhibit large flares in 2008--2009 appropriate for study by GLAST and other instruments are made. The need for long-term data for properly understanding source behaviour is emphasised.

We present the results of extensive multi-waveband monitoring of the blazar 3C~279 between 1996 and 2007 at X-ray energies (2-10 keV), optical R band, and 14.5 GHz, as well as imaging with the Very Long Baseline Array (VLBA) at 43 GHz. In all bands the power spectral density corresponds to "red noise" that can be fit by a single power law over the sampled time scales. Variations in flux at all three wavebands are significantly correlated. The time delay between high and low frequency bands changes substantially on time scales of years. A major multi-frequency flare in 2001 coincided with a swing of the jet toward a more southerly direction, and in general the X-ray flux is modulated by changes in the position angle of the jet near the core. The flux density in the core at 43 GHz--increases in which indicate the appearance of new superluminal knots--is significantly correlated with the X-ray flux. We decompose the X-ray and optical light curves into individual flares, finding that X-ray leads optical variations (XO) in 6 flares, the reverse occurs in 3 flares (OX), and there is essentially zero lag in 4 flares. Upon comparing theoretical expectations with the data, we conclude that (1) XO flares can be explained by gradual acceleration of radiating electrons to the highest energies; (2) OX flares can result from either light-travel delays of the seed photons (synchrotron self-Compton scattering) or gradients in maximum electron energy behind shock fronts; and (3) events with similar X-ray and optical radiative energy output originate well upstream of the 43 GHz core, while those in which the optical radiative output dominates occur at or downstream of the core.

We present an update of 3C 273's database hosted by the ISDC, completed with data from radio to gamma-ray observations over the last 10 years. We use this large data set to study the multiwavelength properties of this quasar,especially focussing on its variability behaviour. We study the amplitude of the variations and the maximum variability time scales across the broad-band spectrum and correlate the light curves in different bands, specifically with the X-rays, to search for possible connections between the emission at different energies. 3C 273 shows variability at all frequencies, with amplitudes and time scales strongly depending on the energy and being the signatures of the different emission mechanisms. The variability properties of the X-ray band imply the presence of either two separate components (possibly a Seyfert-like and a blazar-like) or at least two parameters with distinct timing properties to account for the X-ray emission below and above ~20 keV. The dominant hard X-ray emission is most probably not due to electrons accelerated by the shock waves in the jet as their variability does not correlate with the flaring millimeter emission, but seems to be associated to long-timescale variations in the optical. This optical component is consistent with being optically thin synchrotron radiation from the base of the jet and the hard X-rays would be produced through inverse Compton processes (SSC and/or EC) by the same electron population. We show evidence that this synchrotron component extends from the optical to the near-infrared domain, where it is blended by emission of heated dust that we find to be located within about 1 light-year from the ultraviolet source.

We have studied the flare characteristics of 55 AGN at 8 different frequency bands between 4.8 and 230 GHz. Our extensive database enables us to study the various observational properties of flares in these sources and compare our results with theoretical models. We visually extracted 159 individual flares from the flux density curves and calculated different parameters, such as the peak flux density and duration, in all the frequency bands. The selection of flares is based on the 22 and 37 GHz data from Mets\"ahovi Radio Observatory and 90 and 230 GHz data from the SEST telescope. Additional lower frequency 4.8, 8, and 14.5 GHz data are from the University of Michigan Radio Observatory. We also calculated variability indices and compared them with earlier studies. The observations seem to adhere well to the shock model, but there is still large scatter in the data. Especially the time delays between different frequency bands are difficult to study due to the incomplete sampling of the higher frequencies. The average duration of the flares is 2.5 years at 22 and 37 GHz, which shows that long-term monitoring is essential for understanding the typical behaviour in these sources. It also seems that the energy release in a flare is independent of the duration of the flare.

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.

The TeV blazar PKS 2155--304 was the subject of an intensive 2 week optical and near-infrared observing campaign in 2004 August with the CTIO 0.9m telescope. During this time, simultaneous X-ray data from RXTE were also obtained. We compare the results of our observations to the results from two previous simultaneous multiwavelength campaigns on PKS 2155-304. We conclude that the correlation between the X-ray and UV/optical variability is strongest and the time lag is shortest (only a few hours) when the object is brightest. As the object becomes fainter, the correlations are weaker and the lags longer, increasing to a few days. Based on the results of four campaigns, we find evidence for a linear relationship between the mean optical brightness and lag time of X-ray and UV/optical events. Furthermore, we assert that this behavior, along with the different multiwavelength flare lag times across different flux states is consistent with a highly relativistic shock propagating down the jet producing the flares observed during a high state. In a quiescent state, the variability is likely to be due to a number of factors including both the jet and contributions outside of the jet, such as the accretion disk.

During the period 1966.5 - 2006.2 the 15GHz and 8GHz lightcurves of 3C454.3 (z=0.859) show a qsasi-periodicity of ~12.8 yr (~6.9 yr in the rest frame of the source) with a double-bump structure. This periodic behaviour is interpreted in terms of a rotating double-jet model in which the two jets are created from the black holes in a binary system and rotate with the period of the orbital motion. The periodic variations in the radio fluxes of 3C454.3 are suggested to be mainly due to the lighthouse effects (or the variation in Doppler boosting) of the precessing jets which are caused by the orbital motion. In addition, variations in the mass-flow rates accreting onto the black holes may be also involved.

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 present a study of variability time scales in a large sample of Active Galactic Nuclei at several frequencies between 4.8 and 230 GHz. We investigate the differences of various AGN types and frequencies and correlate the measured time scales with physical parameters such as the luminosity and the Lorentz factor. Our sample consists of both high and low polarization quasars, BL Lacertae objects and radio galaxies. The basis of this work is the 22 GHz, 37 GHz and 87 GHz monitoring data from the Metsahovi Radio Observatory spanning over 25 years. In addition,we used higher 90 GHz and 230 GHz frequency data obtained with the SEST-telescope between 1987 and 2003. Further lower frequency data at 4.8 GHz, 8 GHz and 14.5 GHz from the University of Michigan monitoring programme have been used. We have applied three different statistical methods to study the time scales: The structure function, the discrete correlation function and the Lomb-Scargle periodogram. We discuss also the differences and relative merits of these three methods. Our study reveals that smaller flux density variations occur in these sources on short time scales of 1-2 years, but larger outbursts happen quite rarely, on the average only once in every 6 years. We do not find any significant differences in the time scales between the source classes. The time scales are also only weakly related to the luminosity suggesting that the shock formation is caused by jet instabilities rather than the central black hole.

We present the results of a multi-frequency analysis of the structural variability in the parsec-scale jet of the blazar S5 1803+784. More than 90 epochs of observations at 6 frequencies from 1.6 GHz up to 22 GHz have been combined and analyzed. We discuss an alternative jet model for the source. In contrast to previously discussed motion scenarios for S5 1803+784, we find that the jet structure within 12 mas of the core can most easily be described by seven "oscillating" jet features. We find that the parameters of jet features, such as core separation, position angle and flux density, change in a periodic way with a timescale of about 4 years. We also find evidence for a correlation between these parameters and the total flux density variability. We suggest a scenario incorporating a periodic form of motion (e.g. rotation, precession), with a non-negligible geometrical contribution to explain the observational results.

15 VLBA observations of the BL Lac object 0716+714, performed between 1994.67 and 2006.40 at 15 GHz, have been analyzed. Part of the data (12 epochs) were obtained within the Mojave-survey project and reanalyzed by us. We present a new motion scenario for jet model component motion in this source which suggests that no longterm-outward motion but rather an oscillation of components around an average core separation is taking place. Although no significant outward motion in the core separation can be found, motion with regard to the position angle is observed. We give lower limits for the derived apparent motions and Doppler factors and compare them with values published in the literature. We find a relation between the total flux-density evolution at 14.5 GHz (UMRAO data) and the position angle changes observed for the jet component closest to the core. We suggest a significant geometric contribution to the longterm flux-density variability in 0716+714.

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

We study the synchrotron flaring behaviour of the blazar 3C279 based on an extensive dataset covering 10 years of monitoring at 19 different frequencies in the radio-to-optical range. The properties of a typical outburst are derived from the observations by decomposing the 19 lightcurves into a series of self-similar events. This analysis is achieved by fitting all data simultaneously to a succession of outbursts defined according to the shock-in-jet model of Marscher & Gear (1985). We compare the derived properties of the synchrotron outbursts in 3C279 to those obtained with a similar method for the quasar 3C273. It is argued that differences in the flaring behaviour of these two sources are intrinsic to the sources themselves rather than being due to orientation effects. We also compare the start times and flux densities of our modelled outbursts with those measured from radio components identified in Very Long Baseline Interferometry (VLBI) images. We find VLBI counterparts for most of our model outbursts, although some high-frequency peaking events are not seen in the radio maps. Finally, we study the link between the appearance of a new synchrotron component and the EGRET gamma-ray state of the source at 10 different epochs. We find that an early-stage shock component is always present during high gamma-ray states, while in low gamma-ray states the time since the onset of the last synchrotron outburst is significantly longer. This statistically significant correlation supports the idea that gamma-ray flares are associated with the early stages of shock components propagating in the jet. We note, however, that the shock wave is already beyond the broad line region during the gamma-ray flaring.

In this letter we present the results of a wavelet analysis of the radio light curve of the BL Lac Object PKS B0048-097 from the University of Michigan Radio Astronomy Observatory monitoring program at 8GHz during twenty-five years, from 1979 to 2004. The results show a remarkable periodicity of 450-470 days in the early 1980s that changed to a ~585 day periodicity in the late 1980s to early 1990s. A less pronounced ~400 day periodicity is found after \~1995. Very-long-baseline interferometry imaging at 15GHz shows dramatic structural changes in the usually unresolved source between two epochs, 1995.57 and 2002.38. The pronounced northward directed jet seen in the 2002 image differs by more than 90 degrees in direction from the source structure found in earlier epochs. These findings make PKS B0048-097 a primary target for multi-wavelength observations and intensive radio monitoring to decipher the blazar-variability phenomenon.

BL Lacertae has been the target of several observing campaigns by the Whole Earth Blazar Telescope (WEBT) collaboration and is one of the best studied blazars at all accessible wavelengths. A recent analysis of the optical and radio variability indicates that part of the radio variability is correlated with the optical light curve. Here we present an analysis of a huge VLBI data set including 108 images at 15, 22, and 43 GHz obtained between 1995 and 2004. The aim of this study is to identify the different components contributing to the single-dish radio light curves. We obtain separate radio light curves for the VLBI core and jet and show that the radio spectral index of single-dish observations can be used to trace the core variability. Cross-correlation of the radio spectral index with the optical light curve indicates that the optical variations lead the radio by about 100 days at 15 GHz. By fitting the radio time lags vs. frequency, we find that the power law is steeper than expected for a freely expanding conical jet in equipartition with energy density decreasing as the square of the distance down the jet as in the K\"onigl model. There is a section of the compact radio jet where the emission is weak such that flares propagating down the jet are bright first in the core region with a secondary increase in flux about 1.0 mas from the core. This illustrates the importance of direct imaging to the interpretation of multi-wavelength light curves that can be affected by several distinct components at any given time. We discuss how the complex behaviour of the light curves and correlations can be understood within the framework of a precessing helical jet model.

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

PG 1553+11 was the target of a coordinated three week multiwavelength campaign during 2003 April and May. A significant X-ray flare was observed during the second half of this campaign. Although no optical flare was recorded during the X-ray campaign, optical observations obtained immediately prior to the campaign displayed a higher flux than that recorded during the campaign. An optical flare was observed a few days after the end of the X-ray campaign and may be related to the X-ray flare. Radio observations were made at three frequencies, with no significant changes in flux detected near the times of the optical and X-ray flares. The spectral energy distributions and flux ratios in different wavebands observed for this object are compared to other X-ray selected blazars to demonstrate how PG 1553+11 is an extreme member of this group.

Aims. The aim of this work is to identify and study spectral features in the quasar 3C 273 usually blended by its strong jet emission. Method. A historic minimum in the sub-millimetre emission of 3C 273 triggered coordinated multi-wavelength observations in June 2004. X-ray observations from the INTEGRAL, XMM-Newton and RXTE satellites are complemented by ground-based optical, infrared, millimetre and radio observations. The overall spectrum is used to model the infrared and X-ray spectral components. Results. Three thermal dust emission components are identified in the infrared. The dust emission on scales from 1 pc to several kpc is comparable to that of other quasars, as expected by AGN unification schemes. The observed weakness of the X-ray emission supports the hypothesis of a synchrotron self-Compton origin for the jet component. There is a clear soft-excess and we find evidence for a very broad iron line which could be emitted in a disk around a Kerr black hole. Other signatures of a Seyfert-like X-ray component are not detected.

We report on a multiwavelength campaign on the TeV gamma-ray blazar Markarian (Mrk) 421 performed during December 2002 and January 2003. These target of opportunity observations were initiated by the detection of X-ray and TeV gamma-ray flares with the All Sky Monitor (ASM) on board the Rossi X-ray Timing Explorer (RXTE) and the 10 m Whipple gamma-ray telescope.The campaign included observational coverage in the radio (University of Michigan Radio Astronomy Observatory), optical (Boltwood, La Palma KVA 0.6m, WIYN 0.9m), X-ray (RXTE pointed telescopes), and TeV gamma-ray (Whipple and HEGRA) bands. At TeV energies, the observations revealed several flares at intermediate flux levels, peaking between 1 and 1.5 times the flux from the Crab Nebula. While the time averaged spectrum can be fitted with a single power law of photon index Gamma =2.8, we find some evidence for spectral variability. Confirming earlier results, the campaign reveals a rather loose correlation between the X-ray and TeV gamma-ray fluxes. In one case, a very strong X-ray flare is not accompanied by a comparable TeV gamma-ray flare. Although the source flux was variable in the optical and radio bands, the sparse sampling of the optical and radio light curves does not allow us to study the correlation properties in detail. We present a simple analysis of the data with a synchrotron-self Compton model, emphasizing that models with very high Doppler factors and low magnetic fields can describe the data.

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

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 report results from {\it XMM-Newton} observations of thirteen X-ray bright BL Lacertae objects, selected from the {\it Einstein} Slew Survey sample. The survey was designed to look for evidence of departures of the X-ray spectra from a simple power law shape (i.e., curvature and/or line features), and to find objects worthy of deeper study. Our data are generally well fit by power-law models, with three cases having hard ($\Gamma<2; dN/dE \propto E^{-\Gamma}$) spectra that indicate synchrotron peaks at $E \gsim 5$ keV. Previous data had suggested a presence of absorption features in the X-ray spectra of some BL Lacs. In contrast, none of these spectra show convincing examples of line features, either in absorption or emission, suggesting that such features are rare amongst BL Lacs, or, more likely, artifacts caused by instrumental effects. We find significant evidence for intrinsic curvature (steepening by $d\Gamma / d({\rm log} E) = 0.4 \pm 0.15$) in fourteen of the seventeen X-ray spectra. This cannot be explained satisfactorily via excess absorption, since the curvature is essentially constant from $0.5-6$ keV, an observation which is inconsistent with the modest amounts of absorption that would be required. We use the {\it XMM-Newton} Optical Monitor data with concurrent radio monitoring to derive broadband spectral energy distributions and peak frequency estimates. From these we examine models of synchrotron emission and model the spectral curvature we see as the result of episodic particle acceleration.

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

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

A very detailed monitoring of a radio flare in the Seyfert I galaxy III Zw 2 with the VLA and the VLBA is presented. The relative astrometry in the VLBA observations was precise on a level of a few microarcseconds. Spectral and spatial evolution of the source are closely linked and these observations allowed us to study in great detail a textbook example of a synchrotron self-absorbed jet. We observe a phase where the jet gets frustrated, without expansion and no spectral evolution. Then the jet breaks free and starts to expand with apparent superluminal motion. This expansion is accompanied by a strong spectral evolution. The results are a good confirmation of synchrotron theory and equipartition for jets.

We report results from {\it XMM-Newton} observations of thirteen X-ray bright BL Lacertae objects, selected from the {\it Einstein} Slew Survey sample. The spectra are generally well fit by power-law models, with four objects having hard ($\alpha<1; F_\nu \propto \nu^{-\alpha}$) spectra that indicates synchrotron peaks at $>5$ keV. None of our spectra show line features, indicating that soft X-ray absorption ``notches'' must be rare amongst BL Lacs, rather than common or ubiquitous as had previously been asserted. We find significant curvature in most of the spectra. This curvature is almost certainly intrinsic, as it appears nearly constant from 0.5 to 6 keV, an observation which is inconsistent with the small columns seen in these sources.

We have carried out a multifrequency analysis of the radio variability of blazars, exploiting the data obtained during the extensive monitoring programs carried out at the University of Michigan Radio Astronomy Observatory (UMRAO, at 4.8, 8, and 14.5 GHz) and at the Metsahovi Radio Observatory (22 and 37 GHz). Two different techniques detect, in the Metsahovi light curves, evidences of periodicity at both frequencies for 5 sources (0224+671, 0945+408, 1226+023, 2200+420, and 2251+158). For the last three sources consistent periods are found also at the three UMRAO frequencies and the Scargle (1982) method yields an extremely low false-alarm probability. On the other hand, the 22 and 37 GHz periodicities of 0224+671 and 0945+408 (which were less extensively monitored at Metsahovi and for which we get a significant false-alarm probability) are not confirmed by the UMRAO database, where some indications of ill-defined periods about a factor of two longer are retrieved. We have also investigated the variability index, the structure function, and the distribution of intensity variations of the most extensively monitored sources. We find a statistically significant difference in the distribution of the variability index for BL Lac objects compared to flat-spectrum radio quasars (FSRQs), in the sense that the former objects are more variable. For both populations the variability index steadily increases with increasing frequency. The distribution of intensity variations also broadens with increasing frequency, and approaches a log-normal shape at the highest frequencies. We find that variability enhances by 20-30% the high frequency counts of extragalactic radio-sources at bright flux densities, such as those of the WMAP and Planck surveys.

About twenty years of radio observations in five bands (from 4.8 to 37 GHz) of the BL Lac object GC 0109+224 (S2 0109+22, RGB J0112+227), are presented and analysed together with the optical data. Over the past ten years this blazar has exhibited enhanced activity. There is only weak correlation between radio and optical flares delays, usually protracted on longer timescales in the radio with respect to the optical. In some cases no radio flare counterpart was observed for the optical outbursts. The radio variability, characterised by peaks superposition, shows hints of some characteristic timescales (around the 3-4 years), and a fluctuation mode between the flickering and the shot noise. The reconstructed spectral energy distribution, poorly monitored at high energies, is preliminarily parameterised with a synchrotron-self Compton description. The smooth synchrotron continuum, peaked in the near-IR-optical bands, strengthens the hypothesis that this source could be an intermediate blazar. Moreover the intense flux in millimetre bands, and the optical and X-ray brightness, might suggest a possible detectable gamma-ray emission.

Following the detection of strong TeV gamma-ray flares from the BL Lac object 1ES 1959+650 with the Whipple 10 m Cherenkov telescope on May 16 and 17, 2002, we performed intensive Target of Opportunity (ToO) radio, optical, X-ray and TeV gamma-ray observations from May 18, 2002 to August 14, 2002. Observations with the X-ray telescope RXTE and the Whipple and HEGRA gamma-ray telescopes revealed several strong flares, enabling us to sensitively test the X-ray/gamma-ray flux correlation properties. Although the X-ray and gamma-ray fluxes seemed to be correlated in general, we found an ``orphan'' gamma-ray flare that was not accompanied by an X-ray flare. After describing in detail the radio (UMRAO, VLA), optical (Boltwood, Abastumani), X-ray (RXTE) and gamma-ray (Whipple, HEGRA) light curves and Spectral Energy Distributions (SEDs) we present initial modeling of the SED with a simple Synchrotron Self-Compton (SSC) model. With the addition of another TeV blazar with good broadband data, we consider the set of all TeV blazars to begin to look for a connection of the jet properties to the properties of the central accreting black hole thought to drive the jet. Remarkably, the temporal and spectral X-ray and gamma-ray emission characteristics of TeV blazars are very similar, even though the masses estimates of their central black holes differ by up to one order of magnitude.

We introduce an algorithm for applying a cross-wavelet transform to analysis of quasiperiodic variations in a time-series, and introduce significance tests for the technique. We apply a continuous wavelet transform and the cross-wavelet algorithm to the Pearson-Readhead VLBI survey sources using data obtained from the University of Michigan 26-m parabloid at observing frequencies of 14.5, 8.0, and 4.8 GHz. Thirty of the sixty-two sources were chosen to have sufficient data for analysis, having at least 100 data points for a given time-series. Of these thirty sources, a little more than half exhibited evidence for quasiperiodic behavior in at least one observing frequency, with a mean characteristic period of 2.4 yr and standard deviation of 1.3 yr. We find that out of the thirty sources, there were about four time scales for every ten time series, and about half of those sources showing quasiperiodic behavior repeated the behavior in at least one other observing frequency.

Eight optical and four radio observatories have been intensively monitoring the BL Lac object 0716+714 in the last years: 4854 data points have been collected in the UBVRI bands since 1994, while radio light curves extend back to 1978. Many of these data are presented here for the first time. The long-term trend shown by the optical light curves seems to vary with a characteristic time scale of about 3.3 years, while a longer period of 5.5-6 years seems to characterize the radio long-term variations. In general, optical colour indices are only weakly correlated with brightness. The radio flux behaviour at different frequencies is similar, but the flux variation amplitude decreases with increasing wavelength. The radio spectral index varies with brightness (harder when brighter), but the radio fluxes seem to be the sum of two different-spectrum contributions: a steady base level and a harder-spectrum variable component. Once the base level is removed, the radio variations appear as essentially achromatic, similarly to the optical behaviour. Flux variations at the higher radio frequencies lead the lower-frequency ones with week-month time scales. The behaviour of the optical and radio light curves is quite different, the broad radio outbursts not corresponding in time to the faster optical ones and the cross-correlation analysis indicating only weak correlation with long time lags. However, minor radio flux enhancements simultaneous with the major optical flares can be recognized, which may imply that the mechanism producing the strong flux increases in the optical band also marginally affects the radio one.

We present the results of the BeppoSAX observations of two BL Lacs, OQ 530 and S5 0716+714, as part of a ToO program for the simultaneous observation at radio, optical, X-ray and TeV energies. Both sources are detected in the LECS and MECS, with S5 0716+714 visible also in the PDS band, up to about 60 keV. The X-ray spectra of both sources are better fitted by a double power-law model, with a steep soft X-ray component flattening at harder energies, with breaks at 0.3 and 1.5 keV, respectively. The concave shape of the spectra in both objects is consistent with soft X-rays being produced by the synchrotron and harder X-rays by the inverse Compton processes. Also the X-ray variability properties confirm this scenario, in particular for S5 0716+714 our observation shows variations by about a factor 3 over one hour below 3 keV and no variability above. Their simultaneous broad band energy spectral distributions can be successfully interpreted within the frame of a homogeneous synchrotron and inverse Compton model, including a possible contribution from an external source of seed photons with the different spectral states of S5 0716+714 being reproduced by changing the injected power. The resulting parameters are fully consistent with the two sources being intermediate objects within the "sequence" scenario proposed for blazars.

Using UMRAO centimeter-band total flux density and linear polarization monitoring observations of the complete Pearson-Readhead extragalactic source sample obtained between August 1984 and March 2001, we identify the range of variability in extragalactic objects as functions of optical and radio morphological classification, and relate total flux density variations to structural changes in published coeval VLBI maps in selected objects. We have identified preferred orientations of the electric vector of the polarized emission at 14.5 and 4.8 GHz in roughly half of the objects, and compared these with orientations of the flow direction indicated by VLBI morphology. When comparing the distributions of the orientation offsets for the BL Lacs and for the QSOs, we find differences in both range and mean value, in support of intrinsic class differences. In the shock-in-jet scenario, we attribute this to the allowed range of obliquities of shocks developing in the flow relative to the flow direction. The fact that we find longterm stability in EVPA over many events implies that a dominant magnetic field orientation persists; in the core-dominated objects, with small contribution from the underlying quiescent jet, this plausibly suggests that the magnetic field has a longterm memory, with subsequent shock events exhibiting similar EVPA orientation, or, alternatively, the presence of a standing shock in the core. We have looked for systematic, monotonic changes in EVPA which might be expected in the emission from a precessing jet, a model currently invoked for some AGNs; none were identified. Further, we carried out a Scargle periodogram analysis of the total flux density observations, but found no strong evidence for periodicity in any of the sample sources.

We interpret the linear polarization structures observed in extragalactic radio sources, even those oriented at oblique angles to the jet flows, to be due to oblique, relativistic shock fronts in the emitting regions. Many sources exhibit indications of such oblique structures, and the goal of this investigation is to test this hypothesis quantitatively. A selected group of ten highly-variable extragalactic sources were observed with the VLBA at 15 and 43 GHz, on nine epochs spanning a 30-month period; five of these objects were also observed at 8.0 and 22 GHz. The integrated total flux densities and linear polarizations of the selected objects were also observed several times a month at 4.8, 8.0 and 14.5 GHz with the University of Michigan 26-meter telescope. All objects exhibited variability with several exhibiting more than one independent outburst during the period. We show the evolution of the polarized components with time and discuss the relativistic shock parameters required to match the observed polarization structures. Even cases where the magnetic field is apparently oriented along the jet flow can be fit by oblique shock models when relativistic aberration effects are included.

Monitoring results are presented for the Stangellini 1 Jy GHz-peaked-spectrum source sample, illustrating that several members exhibit variability in total flux and/or linear polarization over timescales of order a decade. The variability occurs while the spectrum, based on the integrated fluxes, remains steep and characteristic of a transparent source. Total flux variability is unexpected in view of recent VLBI observations indicating no or hidden cores in several members. However, both the variability, and the detection of circular polarization in one class member, argue for the presence of opacity in some portion of the radio jet.

We present two BeppoSAX observations of BL Lacertae as part of a multiwavelength radio-to-TeV campaign. During the first observation we observe a faint Compton spectrum, while during the second, we detect a synchrotron spectrum with the highest [2-10] keV flux ever measured; above 10 keV an inverse Compton component begin to dominate. The synchrotron flux is very variable with time scales of 1 hr. We describe four different SED shifting the synchrotron peak both in frequency and flux intensity and we sketch a scenario in which a blob moves along a jet and can be located in or outside the BLR. This implies different radiative mechanism at work, SSC or external Compton, producing different high energy spectra.

We have organized an extensive multiwavelength campaign on BL Lacertae in the second half of 2000. Simultaneous or quasi-simultaneous observations were taken at radio frequencies, in the optical - carried out by the WEBT collaboration -, in X-rays - using BeppoSAX and RXTE -, and at VHE gamma-rays with the CAT and HEGRA Cherenkov telescope facilities. In this paper, we are presenting first results from this campaign. The WEBT optical campaign achieved an unprecedented time coverage, virtually continuous over several 10 - 20 hour segments, and revealed intraday variability on time scales of ~ 1.5 hours. The X-ray observations of Nov. 1 - 2, 2000, revealed significant variability on similar time scales, and provided evidence for the synchrotron spectrum extending out to ~ 10 keV during that time. At higher energies, the onset of the hard power-law commonly observed in radio-loud quasars and radio-selected BL Lac objects is seen. From the energy-resolved X-ray variability patterns, we find evidence that electron cooling might not be dominated by one single radiation mechanism; rather, a combination of synchrotron and Compton cooling seems to be required.

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

New optical and radio data on the BL Lacertae object AO 0235+16 have been collected in the last four years by a wide international collaboration, which confirm the intense activity of this source. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) first-light campaign organized in November 1997. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection, Discrete Correlation Function analysis, and Discrete Fourier Transform technique. The major radio outbursts are found to repeat quasi-regularly with a periodicity of about 5.7 years; this period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them.

We examine the coincidence of times of high $\gamma$-ray flux and ejections of superluminal components from the core in EGRET blazars based on a VLBA monitoring program at 22 and 43 GHz from November 1993 to July 1997. In 23 cases of $\gamma$-ray flares for which sufficient VLBA data exist, 10 of the flares (in 8 objects) fall within 1$\sigma$ uncertainties of the extrapolated epoch of zero separation from the core of a superluminal radio component. In each of two sources (0528+134 and 1730-130) two successive $\gamma$-ray flares were followed by the appearance of new superluminal components. We carried out statistical simulations which show that if the number of coincidences $\ge$ 10 the radio and $\gamma$-ray events are associated with each other at >99.999% confidence. Our analysis of the observed behavior, including variability of the polarized radio flux, of the sources before, during, and after the $\gamma$-ray flares suggests that the $\gamma$-ray events occur in the superluminal radio knots. This implies that the $\gamma$-ray flares are caused by inverse Compton scattering by relativistic electrons in the parsec-scale regions of the jet rather than closer to the central engine.

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

So far, all relativistically boosted jets with superluminal motion have only been detected in typical radio galaxies with early type host galaxies. We have now discovered superluminal motion in the Seyfert I galaxy III Zw 2, classified as a spiral. The lower limit for the apparent expansion speed is 1.25 c. Spectral and spatial evolution are closely linked. Before and after this rapid expansion we have seen a period of virtually no expansion with an expansion speed less than 0.04 c. However, at 15 GHz the picture is completely different. III Zw 2 shows slow expansion ($\sim 0.6 c$) during the time of no expansion at 43 GHz and no expansion during the rapid expansion at 43 GHz. The difference between the two frequencies is qualitatively explained by optical-depth effects in an 'inflating- balloon model', describing the evolution of radio lobes on an ultra-compact scale. The stop-and-go behavior could be explained by a jet interacting with a molecular cloud or the molecular torus. Since III Zw 2 is also part of a sample of so- called radio-intermediate quasars (RIQ), it confirms earlier predictions of superluminal motion for this source, based on the argument that RIQs could be relativistically boosted jets in radio-weak quasars and Seyfert galaxies.

AO 0235+164 is a very compact, flat spectrum radio source identified as a BL Lac object at a redshift of z=0.94. It is one of the most violently variable extragalactic objects at both optical and radio wavelengths. The radio structure of the source revealed by various ground-based VLBI observations is dominated by a nearly unresolved compact component at almost all available frequencies. Dual-frequency space VLBI observations of AO 0235+164 were made with the VSOP mission in January-February 1999. The array of the Japanese HALCA satellite and co-observing ground radio telescopes in Australia, Japan, China and South Africa allowed us to study AO 0235+164 with an unprecedented angular resolution at frequencies of 1.6 and 5 GHz. We report on the sub-milliarcsecond structural properties of the source. The 5-GHz observations led to an estimate of T_B > 5.8 x 10^{13} K for the rest-frame brightness temperature of the core, which is the highest value measured with VSOP to date.

So far all relativistically boosted jets with superluminal motion have only been detected in typical radio galaxies with early type host galaxies. We have now discovered superluminal motion in the Seyfert I galaxy III Zw 2, classified as a spiral. Superluminal motion was first inferred from the spectral evolution of the source and then confirmed by VLBI observations. The lower limit for the apparent expansion speed is 1.25 +/- 0.09 c. The fact that the spectral and spatial evolution are closely linked demonstrates that we are dealing with real physical expansion. Prior to this rapid expansion we have seen a period of virtually no expansion with an expansion speed less than 0.04 c. Since III Zw 2 is also part of a sample of so called radio-intermediate quasars (RIQ), it confirms earlier predictions of superluminal motion for this source, based on the argument that RIQs could be relativistically boosted jets in radio-weak quasars and Seyfert galaxies.

We present radio observations at three frequencies and contemporaneous optical monitoring of the peculiar BL Lac object AO 0235+164. During a three-week campaign with the VLA we observed intraday variability in this source and found a distinct peak which can be identified throughout the radio frequencies and tentatively connected to the R-band variations. This event is characterized by unusual properties: its strength increases, and its duration decreases with wavelength, and it peaks earlier at 20cm than at 3.6 and 6cm. We discuss several generic models (a `standard' shock-in-jet model, a precessing beam, free-free-absorption in a foreground screen, interstellar scattering, and gravitational microlensing), and explore whether they can account for our observations. Most attempts at explaining the data on 0235+164 require an extremely small source size, which can be reconciled with the $10^{12}$ K inverse Compton limit only when the Doppler factor of the bulk flow is of order 100. However, none of the models is completely satisfactory, and we suggest that the observed variability is due to a superposition of intrinsic and propagation effects.

III Zw 2 is a spiral galaxy with an optical spectrum and faint extended radio structure typical of a Seyfert galaxy, but also with an extremely variable, blazar-like radio core. We have now discovered a new radio flare where the source has brightened more than twenty-fold within less than two years. A broad-band radio spectrum between 1.4 and 666 GHz shows a textbook-like synchrotron spectrum peaking at 43 GHz, with a self-absorbed synchrotron spectral index +2.5 at frequencies below 43 GHz and an optically thin spectral index -0.75 at frequencies above 43 GHz. The outburst spectrum can be well fitted by two homogenous, spherical components with equipartition sizes of 0.1 and 0.2 pc at 43 and 15 GHz, and with magnetic fields of 0.4 and 1 Gauss. VLBA observations at 43 GHz confirm this double structure and these sizes. Time scale arguments suggest that the emitting regions are shocks which are continuously accelerating particles. This could be explained by a frustrated jet scenario with very compact hotspots. Similar millimeter-peaked spectrum (MPS) sources could have escaped our attention because of their low flux density at typical survey frequencies and their strong variability.

We present multiwaveband observations of a well selected sample of 28 quasars and two radio galaxies with flat radio spectra and strong millimeter wave emission (referred to here as FSRQ's). The observations include multifrequency VLBI measurements, X-ray observations with ROSAT and submillimeter observations with the JCMT. Particularly interesting among many findings is a correlation between the X-ray to millimeter spectral index and fraction of flux density contained in the VLBI core. This tendency toward higher X-ray fluxes from sources with stronger jet emission implies that the knots in the jet are the prominent source of X-rays.

Results from a long-term program to quantify the range of behavior of the cm-wavelength total flux and linear polarization variability properties of a sample of 41 radio-loud BL Lac objects using weekly to tri-monthly observations with the University of Michigan 26-m telescope operating at 14.5, 8.0, and 4.8 GHz are presented; these observations are used to identify class-dependent differences between these BL Lacs and QSOs in the Pearson-Readhead sample. The BL Lacs are found to be more highly variable in total flux density than the QSOs, exhibiting changes that are often nearly-simultaneous and of comparable amplitude at 14.5 and 4.8 GHz in contrast to the behavior in the QSOs and supporting the existence of class-dependent differences in opacity within the parsec-scale jet flows. Structure function analyses of the flux observations quantify that a characteristic timescale is identifiable in only 1/3 of the BL Lacs. The time-averaged fractional linear polarizations are only on the order of a few percent and are consistent with the presence of tangled magnetic fields within the emitting regions. In many sources a preferred long-term orientation of the EVPA is present; when compared with the VLBI structural axis, no preferred position angle difference is identified. The polarized flux typically exhibits variability with timescales of months to a few years and shows the signature of a propagating shock during several resolved outbursts. The observations indicate that the source emission is predominately due to evolving source components and support the occurrence of more frequent shock formation in BL Lac parsec-scale flows than in QSO jets. The differences in variability behavior and polarization between BL Lacs and QSOs can be explained by differences in jet stability.

We present a wide multi-wavelength database of most observations of the quasar 3C 273 obtained during the last 30 years. This database is the most complete set of observations available for an active galactic nucleus (AGN). It contains nearly 20'000 observations grouped together into 70 light curves covering 16 orders of magnitude in frequency from the radio to the gamma-ray domain. The database is constituted of many previously unpublished observations and of most publicly available data gathered in the literature and on the World Wide Web (WWW). It is complete to the best of our knowledge, except in the optical (UBV) domain where we chose not to add all observations from the literature. In addition to the photometric data, we present the spectra of 3C 273 obtained by the International Ultraviolet Explorer (IUE) satellite. In the X-ray domain, we used the spectral fit parameters from the literature to construct the light curves. Apart from describing the data, we show the most representative light curves and the average spectrum of 3C 273. The database is available on the WWW in a homogeneous and clear form and we wish to update it regularly by adding new observations.

We use a continuous wavelet transform to analyze more than two decades of data for the BL~Lac object OJ~287 acquired as part of the UMRAO variability program. We find clear evidence for a persistent modulation of the total flux and polarization with period $\sim 1.66$ years, and for another signal that dominates activity in the 1980s with period $\sim 1.12$ years. The relationship between these two variations can be understood in terms of a `shock-in-jet' model, in which the longer time scale periodicity is associated with an otherwise quiescent jet, and the shorter time scale activity is associated with the passage of a shock; the different periodicities of these two components may reflect different internal conditions of the two flow domains, leading to different wave speeds, or different contractions of a single underlying periodicity, due to the different Doppler factors of the two flow components. We suggest that the modulation arises from a wave driven by some asymmetric disturbance close to the central engine. The periodic behavior in polarization exhibits excursions in $U$ which correspond to a direction $\sim 45^{\circ}$ from the VLBI jet axis. This behavior is not explained by the random walk in the $Q$-$U$ plane which is expected from models in which a pattern of randomly aligned magnetic field elements propagate across the visible portion of the flow, and suggests a small amplitude, cyclic variation in the flow direction in that part of the flow that dominates cm-wavelength emission.

Optical, near-infrared, and radio observations of the BL Lac object PKS2155-304 were obtained simultaneously with a continuous UV/EUV/X-ray monitoring campaign in 1994 May. Further optical observations were gathered throughout most of 1994. The radio, millimeter, and near-infrared data show no strong correlations with the higher energies. The optical light curves exhibit flickering of 0.2-0.3 mag on timescales of 1-2 days, superimposed on longer timescale variations. Rapid variations of ~0.01 mag/min, which, if real, are the fastest seen to date for any BL Lac object. Small (0.2-0.3 mag) increases in the V and R bands occur simultaneously with a flare seen at higher energies. All optical wavebands (UBVRI) track each other well over the period of observation with no detectable delay. For most of the period the average colors remain relatively constant, although there is a tendency for the colors (in particular B-V) to vary more when the source fades. In polarized light, PKS 2155-304 showed strong color dependence and the highest optical polarization (U = 14.3%) ever observed for this source. The polarization variations trace the flares seen in the ultraviolet flux.