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

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

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 on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

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.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We present the results of extensive modeling of the spectral energy distributions (SEDs) of the closest blazar (z=0.031) Markarian 421 (Mrk 421) during a giant outburst in February 2010. The source underwent rapid flux variations in both X-rays and very high energy (VHE) gamma-rays as it evolved from a low-flux state on 2010 February 13-15 to a high-flux state on 2010 February 17. During this period, the source exhibited significant spectral hardening from X-rays to VHE gamma-rays while exhibiting a "harder when brighter" behavior in these energy bands. We reproduce the broadband SED using a time-dependent multi-zone leptonic jet model with radiation feedback. We find that an injection of the leptonic particle population with a single power-law energy distribution at shock fronts followed by energy losses in an inhomogeneous emission region is suitable for explaining the evolution of Mrk 421 from low- to high-flux state in February 2010. The spectral states are successfully reproduced by a combination of a few key physical parameters, such as the maximum $\&$ minimum cutoffs and power-law slope of the electron injection energies, magnetic field strength, and bulk Lorentz factor of the emission region. The simulated light curves and spectral evolution of Mrk 421 during this period imply an almost linear correlation between X-ray flux at 1-10 keV energies and VHE gamma-ray flux above 200 GeV, as has been previously exhibited by this source. Through this study, a general trend that has emerged for the role of physical parameters is that, as the flare evolves from a low- to a high-flux state, higher bulk kinetic energy is injected into the system with a harder particle population and a lower magnetic field strength.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

We present the first LOFAR observations of the radio jet in the quasar 4C+19.44 (a.k.a. PKS 1354+19) obtained with the long baselines. The achieved resolution is very well matched to that of archival Jansky Very Large Array (JVLA) observations at higher radio frequencies as well as the archival X-ray images obtained with {\it Chandra}. We found that, for several knots along the jet, the radio flux densities measured at hundreds of MHz lie well below the values estimated by extrapolating the GHz spectra. This clearly indicates the presence of spectral curvature. Radio spectral curvature has been already observed in different source classes and/or extended radio structures and it has been often interpreted as due to intrinsic processes, as a curved particle energy distribution, rather than absorption mechanisms ({ Razin-Tsytovich} effect, free-free or synchrotron self absorption to name a few). Here we discuss our results according to the scenario where particles undergo stochastic acceleration mechanisms also in quasar jet knots.

The quasar 3C454.3 underwent a uniquely-structured multi-frequency outburst in June 2016. The blazar was observed in the optical $R$ band by several ground-based telescopes in photometric and polarimetric modes, at $\gamma$-ray frequencies by the \emph{Fermi}\ Large Area Telescope, and at 43 GHz with the Very Long Baseline Array. The maximum flux density was observed on 2016 June 24 at both optical and $\gamma$-ray frequencies, reaching $S^\mathrm{max}_\mathrm{opt}=18.91\pm0.08$ mJy and $S_\gamma^\mathrm{max} =22.20\pm0.18\times10^{-6}$ ph cm$^{-2}$ s$^{-1}$, respectively. The June 2016 outburst possessed a precipitous decay at both $\gamma$-ray and optical frequencies, with the source decreasing in flux density by a factor of 4 over a 24-hour period in $R$ band. Intraday variability was observed throughout the outburst, with flux density changes between 1 and 5 mJy over the course of a night. The precipitous decay featured statistically significant quasi-periodic micro-variability oscillations with an amplitude of $\sim 2$-$3\%$ about the mean trend and a characteristic period of 36 minutes. The optical degree of polarization jumped from $\sim3\%$ to nearly 20\% during the outburst, while the position angle varied by $\sim120\degr$. A knot was ejected from the 43 GHz core on 2016 Feb 25, moving at an apparent speed $v_\mathrm{app}=20.3c\pm0.8c$. From the observed minimum timescale of variability $\tau_\mathrm{opt}^\mathrm{min}\approx2$ hr and derived Doppler factor $\delta=22.6$, we find a size of the emission region $r\lesssim2.6\times10^{15}$ cm. If the quasi-periodic micro-variability oscillations are caused by periodic variations of the Doppler factor of emission from a turbulent vortex, we derive a rotational speed of the vortex $\sim0.2c$.

Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to shed light on the crucial physical processes giving rise to the jet formation, as well as to its extraordinary radiation output up to gamma-ray energies. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm-VLBI polarimetric observations, reaching an unprecedented resolution (~50 microarcsec). We also investigate the variability and physical process occurring in the source during the observing period which coincides with a very active state of the source till high-energies. The Faraday rotation analysis between 3 and 7mm shows a gradient in rotation measure with a maximum value of ~6X10^4 rad/m^2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7mm EVPAs orientation is different when the component is exiting the core or crossing a stationary feature at ~0.1 mas. The interaction between the superluminal component and a recollimation shock at ~0.1 mas could have triggered the multi-wavelengths flares. The variability Doppler factor associated with such interaction is large enough to explain the high energy emission, as we infer from the analysis of gamma-ray and X-ray data, and it is in agreement with the Doppler factor obtained to explain the extraordinary optical flare by Raiteri et al.(2017).

The mechanisms producing fast variability of the $\gamma$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $\gamma$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $\gamma$-ray bands. In this paper we present the MAGIC VHE $\gamma$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $\gamma$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $\gamma$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $\gamma$-rays is very similar to two other observed VHE $\gamma$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $\gamma$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $\gamma$-ray flare, but all have some tension with the multiwavelength observations.

We report results from a deep polarization imaging of the nearby radio galaxy 3C$\,$84 (NGC$\,$1275). The source was observed with the Global Millimeter VLBI Array (GMVA) at 86$\,$GHz at an ultra-high angular resolution of $50\mu$as (corresponding to 250$R_{s}$). We also add complementary multi-wavelength data from the Very Long Baseline Array (VLBA; 15 & 43$\,$GHz) and from the Atacama Large Millimeter/submillimeter Array (ALMA; 97.5, 233.0, and 343.5$\,$GHz). At 86$\,$GHz, we measure a fractional linear polarization of $\sim2$% in the VLBI core region. The polarization morphology suggests that the emission is associated with an underlying limb-brightened jet. The fractional linear polarization is lower at 43 and 15$\,$GHz ($\sim0.3-0.7$% and $<0.1$%, respectively). This suggests an increasing linear polarization degree towards shorter wavelengths on VLBI scales. We also obtain a large rotation measure (RM) of $\sim10^{5-6}~{\rm rad/m^{2}}$ in the core at $\gtrsim$43$\,$GHz. Moreover, the VLBA 43$\,$GHz observations show a variable RM in the VLBI core region during a small flare in 2015. Faraday depolarization and Faraday conversion in an inhomogeneous and mildly relativistic plasma could explain the observed linear polarization characteristics and the previously measured frequency dependence of the circular polarization. Our Faraday depolarization modeling suggests that the RM most likely originates from an external screen with a highly uniform RM distribution. To explain the large RM value, the uniform RM distribution, and the RM variability, we suggest that the Faraday rotation is caused by a boundary layer in a transversely stratified jet. Based on the RM and the synchrotron spectrum of the core, we provide an estimate for the magnetic field strength and the electron density of the jet plasma.

We report the results of monitoring of the radio galaxy 3C 120 with the Neil Gehrels Swift Observatory, Very Long Baseline Array, and Mets\"ahovi Radio Observatory. The UV-optical continuum spectrum and R-band polarization can be explained by a superposition of an inverted-spectrum source with a synchrotron component containing a disordered magnetic field. The UV-optical and X-ray light curves include dips and flares, while several superluminal knots appear in the parsec-scale jet. The recovery time of the second dip was longer at UV-optical wavelengths, in conflict with a model in which the inner accretion disk (AD) is disrupted during a dip and then refilled from outer to inner radii. We favor an alternative scenario in which occasional polar alignments of the magnetic field in the disk and corona cause the flux dips and formation of shocks in the jet. Similar to observations of Seyfert galaxies, intra-band time lags of flux variations are longer than predicted by the standard AD model. This suggests that scattering or some other reprocessing occurs. The 37 GHz light curve is well correlated with the optical-UV variations, with a ~20-day delay. A radio flare in the jet occurred in a superluminal knot 0.14 milliarcseconds downstream of the 43 GHz "core," which places the site of the preceding X-ray/UV/optical flare within the core 0.5-1.3 pc from the black hole. The inverted UV-optical flare spectrum can be explained by a nearly mono-energetic electron distribution with energy similar to the minimum energy inferred in the TeV gamma-ray emitting regions of some BL Lacertae objects.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

We analyze the linear polarization of the relativistic jet in BL Lacertae object OJ~287 as revealed by multi-epoch Very Long Baseline Array (VLBA) images at 43 GHz and monitoring observations at optical bands. The electric-vector position angle (EVPA) of the optical polarization matches that at 43 GHz at locations that are often in the compact millimeter-wave "core" or, at other epochs, coincident with a bright, quasi-stationary emission feature $\sim0.2$~milliarcsec ($\sim$0.9~pc projected on the sky) downstream from the core. This implies that electrons with high enough energies to emit optical synchrotron and $\gamma$-ray inverse Compton radiation are accelerated both in the core and at the downstream feature, the latter of which lies $\geq10$~pc from the central engine. The polarization vector in the stationary feature is nearly parallel to the jet axis, as expected for a conical standing shock capable of accelerating electrons to GeV energies.

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

We use a combination of high-resolution very long baseline interferometry (VLBI) radio and multi-wavelength flux density and polarization observations to constrain the physics of the dissipation mechanism powering the broadband flares in 3C 279 during an episode of extreme flaring activity in 2013-2014. Six bright flares superimposed on a long-term outburst are detected at $\gamma$-ray energies. Four of the flares have optical and radio counterparts. The two modes of flaring activity (faster flares sitting on top of a long term outburst) present at radio, optical, $\gamma$-ray frequencies are missing in X-rays. X-ray counterparts are only observed for two flares. The first three flares are accompanied by ejection of a new VLBI component (NC2) suggesting the 43 GHz VLBI core as the site of energy dissipation. Another new component, NC3, is ejected after the last three flares, which suggests that the emission is produced upstream from the core (closer to the black hole). The study therefore indicates multiple sites of energy dissipation in the source. An anti-correlation is detected between the optical percentage polarization (PP) and optical/$\gamma$-ray flux variations, while the PP has a positive correlation with optical/$\gamma$-rays spectral indices. Given that the mean polarization is inversely proportional to the number of cells in the emission region, the PP vs. optical/$\gamma$-ray anti-correlation could be due to more active cells during the outburst than at other times. In addition to the turbulent component, our analysis suggests the presence of a combined turbulent and ordered magnetic field, with the ordered component transverse to the jet axis.

Ever since the discovery by the Fermi mission that active galactic nuclei (AGN) produce copious amounts of high-energy emission, its origin has remained elusive. Using high-frequency radio interferometry (VLBI) polarization imaging, we could probe the magnetic field topology of the compact high-energy emission regions in blazars. A case study for the blazar 3C 279 reveals the presence of multiple gamma-ray emission regions. Pass 8 Fermi-Large Area Telescope (LAT) data are used to investigate the flux variations in the GeV regime; six gamma-ray flares were observed in the source during November 2013 to August 2014. We use the 43 GHz VLBI data to study the morphological changes in the jet. Ejection of a new component (NC2) during the first three gamma-ray flares suggests the VLBI core as the possible site of the high-energy emission. A delay between the last three flares and the ejection of a new component (NC3) indicates that high-energy emission in this case is located upstream of the 43 GHz core (closer to the black hole).

We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $\gamma$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $\gamma$-ray data, together with Swift-XRT, Swift-UVOT, and ground-based GASP-WEBT optical observations, including polarization information, to study the source variability and the overall spectral energy distribution during the $\gamma$-ray flare. The $\gamma$-ray flaring data, compared with as yet unpublished simultaneous optical data which allow to set constraints on the big blue bump disk luminosity, show very high Compton dominance values of $\sim 100$, with a ratio of $\gamma$-ray to optical emission rising by a factor of three in a few hours. The multi-wavelength behavior of the source during the flare challenges one-zone leptonic theoretical models. The new observations during the June 2015 flare are also compared with already published data and non-simultaneous historical 3C 279 archival data.

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

The very-high-energy (VHE, $\gtrsim 100$ GeV) $\gamma$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $\gamma$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $\gamma$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL).

We present total and linearly polarized 3 mm Global mm-VLBI Array images of a sample of blazars and radio galaxies from the VLBA-BU-BLAZAR 7 mm monitoring program designed to probe the innermost regions of active galactic nuclei (AGN) jets and locate the sites of gamma-ray emission observed by the Fermi-LAT. The lower opacity at 3 mm and improved angular resolution, on the order of 50 microarcseconds, allow us to distinguish features in the jet not visible in the 7 mm VLBA data. We also compare two different methods used for the calibration of instrumental polarisation and we analyze the resulting images for some of the sources in the sample.

Blazars exhibit flares across the entire electromagnetic spectrum. Many $\gamma$-ray flares are highly correlated with flares detected at longer wavelengths; however, a small subset appears to occur in isolation, with little or no correlated variability at longer wavelengths. These "orphan" $\gamma$-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. Macdonald et al. have developed the Ring of Fire model to explain the origin of orphan $\gamma$-ray flares from within blazar jets. In this model, electrons contained within a blob of plasma moving relativistically along the spine of the jet inverse-Compton scatter synchrotron photons emanating off of a ring of shocked sheath plasma that enshrouds the jet spine. As the blob propagates through the ring, the scattering of the ring photons by the blob electrons creates an orphan $\gamma$-ray flare. This model was successfully applied to modeling a prominent orphan $\gamma$-ray flare observed in the blazar PKS 1510$-$089. To further support the plausibility of this model, Macdonald et al. presented a stacked radio map of PKS 1510$-$089 containing the polarimetric signature of a sheath of plasma surrounding the spine of the jet. In this paper, we extend our modeling and stacking techniques to a larger sample of blazars: 3C 273, 4C 71$.$01, 3C 279, 1055$+$018, CTA 102, and 3C 345, the majority of which have exhibited orphan $\gamma$-ray flares. We find that the model can successfully reproduce these flares, while our stacked maps reveal the existence of jet sheaths within these blazars.

We analyze the parsec-scale jet kinematics from 2007 June to 2013 January of a sample of $\gamma$-ray bright blazars monitored roughly monthly with the Very Long Baseline Array at 43~GHz. In a total of 1929 images, we measure apparent speeds of 252 emission knots in 21 quasars, 12 BL~Lacertae objects (BLLacs), and 3 radio galaxies, ranging from 0.02$c$ to 78$c$; 21\% of the knots are quasi-stationary. Approximately 1/3 of the moving knots execute non-ballistic motions, with the quasars exhibiting acceleration along the jet within 5~pc (projected) of the core, and knots in the BLLacs tending to decelerate near the core. Using apparent speeds of components and timescales of variability from their light curves, we derive physical parameters of 120 superluminal knots, including variability Doppler factors, Lorentz factors, and viewing angles. We estimate the half-opening angle of each jet based on the projected opening angle and scatter of intrinsic viewing angles of knots. We determine characteristic values of physical parameters for each jet and AGN class based on the range of values obtained for individual features. We calculate intrinsic brightness temperatures of the cores, $T_{\rm b,int}^{\rm core}$, at all epochs, finding that the radio galaxies usually maintain equipartition conditions in the cores, while $\sim$30\% of $T_{\rm b,int}^{\rm core}$ measurements in the quasars and BLLacs deviate from equipartition values by a factor $>$10. This probably occurs during transient events connected with active states. In the Appendix we briefly describe the behavior of each blazar during the period analyzed.

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

PKS1510-089 is a flat spectrum radio quasar located at a redshift of 0.36. It is one of only a few such sources detected in very-high-energy (VHE, >100 GeV) gamma rays. Though PKS1510-089 is highly variable at GeV energies, until recently no variability has been observed in the VHE band. In 2015 May PKS1510-089 showed a high state in optical and in the GeV range. A VHE gamma-ray flare was detected with MAGIC at that time, showing the first instance of VHE gamma-ray flux variability on the time scale of days in this source. We will present the MAGIC results from this observation, discuss their temporal and spectral properties in the multi-wavelength context and present modelling of such emission in the external Compton scenario.

We present arc-second-resolution data in the radio, IR, optical and X-ray for 4C+19.44 (=PKS 1354+195), the longest and straightest quasar jet with deep X-ray observations. We report results from radio images with half to one arc-second angular resolution at three frequencies, plus HST and Spitzer data. The Chandra data allow us to measure the X-ray spectral index in 10 distinct regions along the 18 arcsec jet and compare with the radio index. The radio and X-ray spectral indices of the jet regions are consistent with a value of $\alpha =0.80$ throughout the jet, to within 2 sigma uncertainties. The X-ray jet structure to the south extends beyond the prominent radio jet and connects to the southern radio lobe, and there is extended X-ray emission in the direction of the unseen counter jet and coincident with the northern radio lobe. This jet is remarkable since its straight appearance over a large distance allows the geometry factors to be taken as fixed along the jet. Using the model of inverse Compton scattering of the cosmic microwave background (iC/CMB) by relativistic electrons, we find that the magnetic field strengths and Doppler factors are relatively constant along the jet. If instead the X-rays are synchrotron emission, they must arise from a population of electrons distinct from the particles producing the radio synchrotron spectrum.

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

We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations.

S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13{\sigma} significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to construct the broad-band spectral energy distribution of S5 0716+714 during its brightest outburst. In this work we will present the preliminary analysis of MAGIC and Fermi-LAT data of the flaring activity in January and February 2015 for the HE (0.1 < HE < 300 GeV) and VHE band, together with radio (Mets\"ahovi, OVRO, VLBA, Effelsberg), sub-millimeter (SMA), optical (Tuorla, Perkins, Steward, AZT-8+ST7, LX-200, Kanata), X-ray and UV (Swift-XRT and UVOT), in the same time-window and discuss the time variability of the multiwavelength light curves during this impressive outburst.

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

The high radio frequency polarization imaging of non-thermal emission from active galactic nuclei (AGN) is a direct way to probe the magnetic field strength and structure in the immediate vicinity of supermassive black holes (SMBHs) and is crucial in testing the jet-launching scenario. To explore the the magnetic field configuration at the base of jets in blazars, we took advantage of the full polarization capabilities of the Global Millimeter VLBI Array (GMVA). With an angular resolution of $\sim$50 micro-arcseconds ($\mu$as) at 86 GHz, one could resolve scales up to $\sim$450 gravitational radii (for a 10$^9$ solar mass black hole at a redshift of 0.1). We present here the preliminary results of our study on the blazar BL~Lac. Our results suggest that on sub-mas scales the core and the central jet of BL Lac are significantly polarized with two distinct regions of polarized intensity. We also noted a great morphological similarity between the 7mm/3mm VLBI images at very similar angular resolution.

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models.

The Gamma-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features including an apparent (~100 deg) position angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps together with cm/mm total intensity and Gamma-ray observations from Fermi/LAT from 2008-2014, we aimed to determine the location of Gamma-ray emission and to explain the inner-mas structural changes. Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and allow us to observe above the synchrotron self-absorption peak frequency. The jet was spectrally decomposed at multiple locations along the jet. From this we derived estimates of the magnetic field. How the field decreases down the jet allowed an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics we attempt to locate the site of Gamma-ray activity, radio flares and spectral changes. Strong Gamma-ray flares appeared to originate from either the "core" region, a downstream stationary feature, or both, with Gamma-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be >1.6 G in the "core" and >0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the "core" as >6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.

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

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.

We present the first polarimetric space VLBI imaging observations at 22 GHz. BL Lacertae was observed in 2013 November 10 with the RadioAstron space VLBI mission, including a ground array of 15 radio telescopes. The instrumental polarization of the space radio telescope is found to be within 9%, demonstrating the polarimetric imaging capabilities of RadioAstron at 22 GHz. Ground-space fringes were obtained up to a projected baseline distance of 7.9 Earth's diameters in length, allowing us to image the jet in BL Lacertae with a maximum angular resolution of 21 $\mu$as, the highest achieved to date. We find evidence for emission upstream of the radio core, which may correspond to a recollimation shock at about 40 $\mu$as from the jet apex, in a pattern that includes other recollimation shocks at approximately 100 $\mu$as and 250 $\mu$as from the jet apex. Polarized emission is detected in two components within the innermost 0.5 mas from the core, as well as in some knots 3 mas downstream. Faraday rotation analysis, obtained from combining RadioAstron 22 GHz and ground-based 15 GHz and 43 GHz images, shows a gradient in rotation measure and Faraday corrected polarization vector as a function of position angle with respect to the core, suggesting that the jet in BL Lacertae is threaded by a helical magnetic field. The intrinsic de-boosted brightness temperature in the unresolved core exceeds $3\!\times\!10^{12}$ K, suggesting at the very least departure from equipartition of energy between the magnetic field and radiating particles.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright $\gamma$-ray outburst detected by the {\it Fermi} Large Area Telescope in September-October 2012 when the source reached a flux of F$_{>100~\mathrm{MeV}} =5.2\pm0.4\times10^{-6}$ photons cm$^{-2}$ s$^{-1}$. At the same time the source displayed an unprecedented optical and NIR outburst. We study the evolution of the parsec scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity Very Long Baseline Array images at 43 GHz, covering the observing period from June 2007 to June 2014. We find that the $\gamma$-ray outburst is coincident with flares at all the other frequencies and is related to the passage of a new superluminal knot through the radio core. The powerful $\gamma$-ray emission is associated with a change in direction of the jet, which became oriented more closely to our line of sight ($\theta\sim$1.2$^{\circ}$) during the ejection of the knot and the $\gamma$-ray outburst. During the flare, the optical polarized emission displays intra-day variability and a clear clockwise rotation of EVPAs, which we associate with the path followed by the knot as it moves along helical magnetic field lines, although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We locate the $\gamma$-ray outburst a short distance downstream of the radio core, parsecs from the black hole. This suggests that synchrotron self-Compton scattering of near-infrared to ultraviolet photons is the probable mechanism for the $\gamma$-ray production.

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings.

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.

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

Blazars exhibit flares across the electromagnetic spectrum. Many $\gamma$-ray flares are highly correlated with flares detected at optical wavelengths; however, a small subset appears to occur in isolation, with little or no variability detected at longer wavelengths. These "orphan" $\gamma$-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. We present numerical calculations of the time-variable emission of a blazar based on a proposal by Marscher et al. (2010) to explain such events. In this model, a plasmoid ("blob") propagates relativistically along the spine of a blazar jet and passes through a synchrotron-emitting ring of electrons representing a shocked portion of the jet sheath. This ring supplies a source of seed photons that are inverse-Compton scattered by the electrons in the moving blob. The model includes the effects of radiative cooling, a spatially varying magnetic field, and acceleration of the blob's bulk velocity. Synthetic light curves produced by our model are compared to the observed light curves from an orphan flare that was coincident with the passage of a superluminal knot through the inner jet of the blazar PKS 1510$-$089. In addition, we present Very Long Baseline Array polarimetric observations that point to the existence of a jet sheath in PKS 1510$-$089, thus providing further observational support for the plausibility of our model. An estimate of the bolometric luminosity of the sheath within PKS 1510$-$089 is made, yielding $L_{\rm sh} \sim 3 \times 10^{45} ~ \rm erg ~ \rm s^{-1}$. This indicates that the sheath within PKS 1510$-$089 is potentially a very important source of seed photons.

The Global mm-VLBI Array (GMVA) is a network of 14 3\,mm and 7\,mm capable telescopes spanning Europe and the United States, with planned extensions to Asia. The array is capable of sensitive maps with angular resolution often exceeding 50\,$\mu$as. Using the GMVA, a large sample of prominent $\gamma$-ray blazars have been observed approximately 6 monthly from later 2008 until now. Combining 3\,mm maps from the GMVA with near-in-time 7\,mm maps from the VLBA-BU-BLAZAR program and 2\,cm maps from the MOJAVE program, we determine the sub-pc morphology and high frequency spectral structure of $\gamma$-ray blazars. The magnetic field strength can be estimated at different locations along the jet under the assumption of equipartition between magnetic field and relativistic particle energies. Making assumptions on the jet magnetic field configuration (e.g. poloidal or toroidal), we can estimate the separation of the mm-wave "core" and the jet base, and estimate the strength of the magnetic field there. The results of this analysis show that on average, the magnetic field strength decreases with a power-law $B \propto r^{-n}$, $n=0.3 \pm 0.2$. This suggests that on average, the mm-wave "core" is $\sim 1-3$\,pc downstream of the de-projected jet apex and that the magnetic field strength is of the order $B_{\rm{apex}} \sim 5-20$\,kG, broadly consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks (MAD)).

In this work we present a polarization analysis at radio frequencies of Markarian 421, one of the closest (z=0.03) TeV blazars. The observations were obtained, both in total and in polarized intensity, with the Very Long Baseline Array (VLBA) at 15, 24, and 43 GHz throughout 2011, with one observation per month (for a total of twelve epochs). We investigate the magnetic field topology and the polarization structure on parsec scale and their evolution with time. We detect polarized emission both in the core and in the jet region, and it varies with frequency, location and time. In the core region we measure a mean fractional polarization of about 1-2%, with a peak of about 4% in March at 43 GHz; the polarization angle is almost stable at 43 GHz, but it shows significant variability in the range 114-173 deg at 15 GHz. In the jet region the polarization properties show a more stable behavior; the fractional polarization is about 16% and the polarization angle is nearly perpendicular to the jet axis. The higher EVPA variability observed at 15 GHz is due both to a variable Faraday rotation effect and to opacity. The residual variability observed in the intrinsic polarization angle, together with the low degree of polarization in the core region, could be explained with the presence of a blend of variable cross-polarized subcomponents within the beam.

We present Fermi-LAT and multi-frequency, multi-epoch VLBA data for the TeV blazar Mrk421. We collected the data during a long and intensive multi-frequency campaign in 2011. We study the gamma-ray light curve, the photon index evolution and their connection to the radio data on sub-parsec scales, including total intensity, polarized flux density, polarization angle, spectral index, and rotation measure both in the core and the jet region. The VLBA data were obtained at 15 and 24 GHz for 12 epochs and at 43 GHz for 23 epochs, thus providing the best temporal and spatial coverage in the radio band ever achieved for a TeV blazar. We provide significant constraints on the jet Doppler factor, the presence of proper motion, the magnetic field configuration, and an intriguing connection between variability in the radio data and the gamma-ray light curve: the total intensity and polarized core emission reach a peak simultaneously to the main gamma-ray peak, followed by a rotation of the polarization angle at low frequency. Opacity-related, long wavelength polarization swings are also detected later in the year, possibly related to secondary peaks in the gamma-ray light curve, setting constraints on the physics of the gamma-ray zone.

The analysis of $\gamma$-ray flux variability along with the parsec-scale jet kinematics suggests that the high-energy radiation in the BL Lac object S5 0716+714 has a significant correlation with the mm-VLBI core flux density and with the local orientation of the inner jet flow. For the first time in any blazar, we report a significant correlation between the $\gamma$-ray flux variations and the variations in the local orientation of the jet outflow (position angle). We find that the $\gamma$-ray flux variations lead the 7~mm VLBI core flux variations by 82$\pm$32~days, which suggests that the high-energy emission in S5 0716+714 is coming from a region located 3.8$\pm$1.9~parsecs closer to the central black hole than the "core" seen on the mm-VLBI images. The results imply a strong physical and casual connection between $\gamma$-ray emission and the inner jet morphology in the source.

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.

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

We report the discovery of Faraday rotation toward radio source 3C 84, the active galactic nucleus in NGC1275 at the core of the Perseus Cluster. The rotation measure (RM), determined from polarization observations at wavelengths of 1.3 and 0.9 mm, is (8.7 +/- 2.3) x 10^5 radians/m^2, among the largest ever measured. The RM remained relatively constant over a 2 year period even as the intrinsic polarization position angle wrapped through a span of 300 degrees. The Faraday rotation is likely to originate either in the boundary layer of the radio jet from the nucleus, or in the accretion flow onto the central black hole. The accretion flow probably is disk-like rather than spherical on scales of less than a parsec, otherwise the RM would be even larger.

We analyzed data in polarized intensity obtained with the Very Long Baseline Array (VLBA) at twelve epochs (one observation per month from January to December 2011) at 15, 24, and 43 GHz. For the absolute orientation of the electric vector position angles (EVPA) we used the D-terms method. We also used gamma-ray data from the Fermi Large Area Telescope on weekly time bins throughout 2011. The source shows polarized emission, and its properties vary with time, frequency, and location along the jet. The core mean polarization fraction is generally between 1% and 2%, with a 4% peak at 43 GHz in March; the polarization angle is variable, mainly at 15 GHz, where it changes frequently, and less so at 43 GHz, where it oscillates in the range 114 - 173 deg. The jet polarization properties are more stable, with a fractional polarization of around 16% and a polarization angle nearly perpendicular to the jet axis. The average flux and photon index at gamma-ray energies are (17.7+-0.5)x10^8 ph cm^-2 s^-1 and the photon index is 1.77+-0.02. The gamma-ray light curve shows variability, with a main peak that appears to be associated with the peak in the core polarized emission at 43 GHz, as well as with the total intensity light curve. A discrete correlation function analysis yields a correlation coefficient of 0.54 at zero delay, with a significance level above 99.7%. We accurately determine the polarization properties of Mrk 421, both in the core and in the jet region. The radio and gamma-ray light curves are correlated. The observed EVPA variability at 15 GHz is partly due to opacity and partly to a variable Faraday rotation effect. To explain the residual variability of the intrinsic polarization angle and the low degree of polarization in the core region, we invoke a blend of variable cross-polarized subcomponents with different polarization properties within the beam.

Using millimeter-very long baseline interferometry (VLBI) observations of the BL Lac object S5 0716+714 from August 2008 to September 2013, we investigate variations in the core flux density and orientation of the sub-parsec scale jet i.e. position angle. The gamma-ray data obtained by the Fermi-LAT (Large Area Telescope) are used to investigate the high-energy flux variations over the same time period. For the first time in any blazar, we report a significant correlation between the gamma-ray flux variations and the position angle (PA) variations in the VLBI jet. The cross-correlation analysis also indicates a positive correlation such that the mm-VLBI core flux density variations are delayed with respect to the gamma-ray flux by 82$\pm$32 days. This suggests that the high-energy emission is coming from a region located $\geq$(3.8$\pm$1.9) parsecs upstream of the mm-VLBI core (closer to the central black hole). These results imply that the observed inner jet morphology has a strong connection with the observed gamma-ray flares.

The blazar 1156+295 was active at gamma-ray energies, exhibiting three prominent flares during the year 2010. Here, we present results using the combination of broadband (X-ray through mm single dish) monitoring data and radio band imaging data at 43 GHz on the connection of gamma-ray events to the ejections of superluminal components and other changes in the jet of 1156+295. The kinematics of the jet over the interval 2007.0-2012.5 using 43 GHz Very Long Baseline Array observations, reveal the presence of four moving and one stationary component in the inner region of the blazar jet. The propagation of the third and fourth components in the jet corresponds closely in time to the active phase of the source in gamma rays. We briefly discuss the implications of the structural changes in the jet for the mechanism of gamma-ray production during bright flares. To localise the gamma-ray emission site in the blazar, we performed the correlation analysis between the 43 GHz radio core and the gamma-ray light curve. The time lag obtained from the correlation constrains the gamma-ray emitting region in the parsec-scale jet.

The Global mm-VLBI Array (GMVA) is the highest angular resolution imaging interferometer currently available as a common user facility. It is capable of angular resolutions on the order of 40 microarcseconds. Currently 14 stations in the United States and Europe participate in global 3 mm VLBI observations. The GMVA is used for continuum and spectroscopic imaging, probing the central regions of active galaxies and the origin of jets as these regions are typically not observable at longer wavelengths due to synchrotron self-absorption. In early 2012, fringes were detected to the three stations of the Korean VLBI Network (KVN), opening the possibility of extending the baseline coverage of the VLBI array to the East. In these proceedings, we will present recent images from a monitoring program of gamma-ray blazars using the GMVA, including the sources 3C454.3 and 0235+164, and an update of its current status and abilities.

We study the multifrequency emission and spectral properties of the quasar 3C 279. We observed 3C 279 in very high energy (VHE, E>100GeV) gamma rays, with the MAGIC telescopes during 2011, for the first time in stereoscopic mode. We combine these measurements with observations at other energy bands: in high energy (HE, E>100MeV) gamma rays from Fermi-LAT, in X-rays from RXTE, in the optical from the KVA telescope and in the radio at 43GHz, 37GHz and 15GHz from the VLBA, Mets\"ahovi and OVRO radio telescopes and optical polarisation measurements from the KVA and Liverpool telescopes. During the MAGIC observations (February to April 2011) 3C 279 was in a low state in optical, X-ray and gamma rays. The MAGIC observations did not yield a significant detection. These upper limits are in agreement with the extrapolation of the HE gamma-ray spectrum, corrected for extragalactic background light absorption, from Fermi-LAT. The second part of the MAGIC observations in 2011 was triggered by a high activity state in the optical and gamma-ray bands. During the optical outburst the optical electric vector position angle rotatated of about 180 degrees. There was no simultaneous rotation of the 43GHz radio polarisation angle. No VHE gamma rays were detected by MAGIC, and the derived upper limits suggest the presence of a spectral break or curvature between the Fermi-LAT and MAGIC bands. The combined upper limits are the strongest derived to date for the source at VHE and below the level of the previously detected flux by a factor 2. Radiation models that include synchrotron and inverse Compton emissions match the optical to gamma-ray data, assuming an emission component inside the broad line region (BLR) responsible for the high-energy emission and one outside the BLR and the infrared torus causing optical and low-energy emission. We interpreted the optical polarisation with a bent trajectory model.

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 present $\gamma$-ray, X-ray, ultraviolet, optical, and near-infrared light curves of 33 $\gamma$-ray bright blazars over four years that we have been monitoring since 2008 August with multiple optical, ground-based telescopes and the Swift satellite, and augmented by data from the Fermi Gamma-ray Space Telescope and other publicly available data from Swift. The sample consists of 21 flat-spectrum radio quasars (FSRQs) and 12 BL Lac objects (BL Lacs). We identify quiescent and active states of the sources based on their $\gamma$-ray behavior. We derive $\gamma$-ray, X-ray, and optical spectral indices, $\alpha_\gamma$, $\alpha_X$, and $\alpha_o$, respectively ($F_\nu\propto\nu^\alpha$), and construct spectral energy distributions (SEDs) during quiescent and active states. We analyze the relationships between different spectral indices, blazar classes, and activity states. We find (i) significantly steeper $\gamma$-ray spectra of FSRQs than for BL Lacs during quiescent states, but a flattening of the spectra for FSRQs during active states while the BL Lacs show no significant change; (ii) a small difference of $\alpha_X$ within each class between states, with BL Lac X-ray spectra significantly steeper than in FSRQs; (iii) a highly peaked distribution of X-ray spectral slopes of FSRQs at $\sim-$0.60, but a very broad distribution of $\alpha_X$ of BL Lacs during active states; (iv) flattening of the optical spectra of FSRQs during quiescent states, but no statistically significant change of $\alpha_o$ of BL Lacs between states; and (v) a positive correlation between optical and $\gamma$-ray spectral slopes of BL Lacs, with similar values of the slopes. We discuss the findings with respect to the relative prominence of different components of high-energy and optical emission as the flux state changes.

We present the results of optical (R band) photometric and polarimetric monitoring and Very Long Baseline Array (VLBA) imaging of the blazar S4 0954+658, along with Fermi and gamma;-ray data during a multi-waveband outburst in 2011 March-April. After a faint state with a brightness level R ~17.6 mag registered in the first half of January 2011, the optical brightness of the source started to rise and reached ~14.8 mag during the middle of March, showing flare-like behavior. The most spectacular case of intranight variability was observed during the night of 2011 March 9, when the blazar brightened by ~0.7 mag within ~7 hours. During the rise of the flux the position angle of optical polarization rotated smoothly over more than 300$\deg$. At the same time, within 1$\sigma$ uncertainty a new superluminal knot appeared with an apparent speed of 19.0$\pm$0.3 c. We have very strong evidence for association of this knot with the multi-waveband outburst in 2011 March-April. We also analyze the multi-frequency behavior of S4 0954+658 during a number of minor outbursts from August 2008 to April 2012. We find some evidence of connections between at least two more superluminal ejecta and near-simultaneous optical flares.

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

We report the results obtained for the AGN Markarian 421 by model-fitting the data in the visibility plane, studing the proper motion of jet components, the light curve, and the spectral index of the jet features. We compare the radio data with optical light curves obtained at the Steward Observatory, considering also the optical polarization information. Mrk 421 has a bright nucleus and a one-sided jet extending towards the north-west for a few parsecs. The model-fits show that brightness distribution is well described using 6-7 circular Gaussian components, four of which are reliably identified at all epochs; all components are effectively stationary except for component D, at ~0.4 mas from the core, whose motion is however subluminal. Analysis of the light curve shows two different states, with the source being brighter and more variable in the first half of 2011 than in the second half. The highest flux density is reached in February. A comparison with the optical data reveals an increase of the V magnitude and of the fractional polarization simultaneous with the enhancement of the radio activity.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent \gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each \gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the \gamma-ray outbursts, although the \gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the \gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the \gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that \gamma-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells.

The extragalactic radio source 1803+784 is a BL Lac object that shows rapid variability and strong linear polarization. Very long baseline interferometry observations reveal that the core possesses a distinctive structure in polarized intensity with two maxima on axis and two minima symmetrically placed on either side. The approximately radial pattern of electric field polarization rods is reminiscent of the results obtained earlier by Cawthorne (2006) for conical shocks, but, individually, these do not reproduce the main features of the polarized intensity images. In numerical simulations and experiments, these shocks occur in pairs and help to stabilize jets as they adjust to changes in environment. Here, the polarization resulting from such structures is investigated using an approximate, analytical approach, by making some simple assumptions about the nature of the flow between two such shock waves. For fairly small viewing angles, it is found that a reasonable representation of the core polarization of 1803+784 can be obtained. The similarity between the observed and model polarization supports the view that the core structure in 1803+784 represents a recollimation shock, and that such shock waves may be responsible for the first disturbance and hence brightening of the quiescent flow in astrophysical jets.

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 report results from 5-day VLBI observations of the well-known quasar 3C 279 at 1.3 mm (230 GHz) in 2011. The measured nonzero closure phases on triangles including stations in Arizona, California and Hawaii indicate that the source structure is spatially resolved. We find an unusual inner jet direction at scales of $\sim$1 parsec extending along the northwest-southeast direction (PA = $127^{\circ}\pm3^{\circ}$), as opposed to other (previously) reported measurements on scales of a few parsecs showing inner jet direction extending to the southwest. The 1.3 mm structure corresponds closely with that observed in the central region of quasi-simultaneous super-resolution VLBA images at 7 mm. The closure phase changed significantly on the last day when compared with the rest of observations, indicating that the inner jet structure may be variable on daily timescales. The observed new direction of the inner jet shows inconsistency with the prediction of a class of jet precession models. Our observations indicate a brightness temperature of $\sim 8\times10^{10}$ K in the 1.3 mm core, much lower than that at centimeter wavelengths. Observations with better uv coverage and sensitivity in the coming years will allow the discrimination between different structure models and will provide direct images of the inner regions of the jet with 20--30 $\mu$as (5--7 light months) resolution.

Global VLBI imaging at millimeter and sub-millimeter wavelength overcomes the opacity barrier of synchrotron self-absorption in AGN and opens the direct view into sub-pc scale regions not accessible before. Since AGN variability is more pronounced at short millimeter wavelength, mm-VLBI can reveal structural changes in very early stages after outbursts. When combined with observations at longer wavelength, global 3mm and 1mm VLBI adds very detailed information. This helps to determine fundamental physical properties at the jet base, and in the vicinity of super-massive black holes at the center of AGN. Here we present new results from multi-frequency mm-VLBI imaging of OJ287 during a major outburst. We also report on a successful 1.3mm VLBI experiment with the APEX telescope in Chile. This observation sets a new record in angular resolution. It also opens the path towards future mm-VLBI with ALMA, which aims at the mapping of the black hole event horizon in nearby galaxies, and the study of the roots of jets in AGN.

We present the results of optical (R band) photometric and polarimetric monitoring and Very Long Baseline Array (VLBA) imaging of the blazar S5 0716+714 along with Fermi gamma-ray data during a multi-waveband outburst in 2011 October. We analyze total and polarized intensity images of the blazar obtained with the VLBA at 43 GHz during and after the outburst. Monotonic rotation of the linear polarization vector at a rate of >50 degrees per night coincided with a sharp maximum in gamma-ray and optical flux. At the same time, within the uncertainties, a new superluminal knot appeared, with an apparent speed of ~21c. The general multi-frequency behavior of the outburst can be explained within the framework of a shock wave propagating along a helical path in the blazar's jet.

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 in 2008. 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 7mm 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 at >12 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.

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

Relativistic jets in AGN in general, and in blazars in particular, are the most energetic and among the most powerful astrophysical objects known so far. Their relativistic nature provides them with the ability to emit profusely at all spectral ranges from radio wavelengths to gamma-rays, as well as to vary extremely at time scales from hours to years. Since the birth of gamma-ray astronomy, locating the origin of gamma-ray emission has been a fundamental problem for the knowledge of the emission processes involved. Deep and densely time sampled monitoring programs with the Fermi Gamma-ray Space Telescope and other facilities at most of the available spectral ranges (including millimeter interferometric imaging and polarization measurements wherever possible) are starting to shed light for the case of blazars. After a short review of the status of the problem, we summarize two of our latest results -obtained from the comprehensive monitoring data compiled by the Boston University Blazar monitoring program - that locate the GeV flaring emission of the BL Lac objects AO 0235+164 and OJ287 within the jets of these blazars, at >12 parsecs from the central AGN engine.

High-resolution radio observations are ideal for constraining the value of physical parameters in the inner regions of active-galactic-nucleus jets and complement results on multiwavelength (MWL) observations. This study is part of a wider multifrequency campaign targeting the nearby TeV blazar Markarian 421 (z=0.031), with observations in the sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and gamma rays (Fermi-LAT, MAGIC, VERITAS). We investigate the jet's morphology and any proper motions, and the time evolution of physical parameters such as flux densities and spectral index. The aim of our wider multifrequency campaign is to try to shed light on questions such as the nature of the radiating particles, the connection between the radio and gamma-ray emission, the location of the emitting regions and the origin of the flux variability. We consider data obtained with the Very Long Baseline Array (VLBA) over twelve epochs (one observation per month from January to December 2011) at 15 GHz and 24 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch. The structure of Mrk 421 is dominated by a compact (~0.13 mas) and bright component, with a one-sided jet detected out to ~10 mas. We identify 5-6 components in the jet that are consistent with being stationary during the 12-month period studied here. Measurements of the spectral index agree with those of other works: they are fairly flat in the core region and steepen along the jet length. Significant flux-density variations are detected for the core component. From our results, we draw an overall scenario in which we estimate a viewing angle 2{\deg} < theta < 5{\deg} and a different jet velocity for the radio and the high-energy emission regions, such that the respective Doppler factors are {\delta}r ~3 and {\delta}h.e. ~14.

We present multiwavelength data of the blazar 3C 454.3 obtained during an extremely bright outburst from November 2010 through January 2011. These include flux density measurements with the Herschel Space Observatory at five submillimeter-wave and far-infrared bands, the Fermi Large Area Telescope at gamma-ray energies, Swift at X-ray, ultraviolet (UV), and optical frequencies, and the Submillimeter Array at 1.3 mm. From this dataset, we form a series of 52 spectral energy distributions (SEDs) spanning nearly two months that are unprecedented in time coverage and breadth of frequency. Discrete correlation anlaysis of the millimeter, far-infrared, and gamma-ray light curves show that the variations were essentially simultaneous, indicative of co-spatiality of the emission, at these wavebands. In contrast, differences in short-term fluctuations at various wavelengths imply the presence of inhomegeneities in physical conditions across the source. We locate the site of the outburst in the parsec-scale core, whose flux density as measured on 7 mm Very Long Baseline Array images increased by 70 percent during the first five weeks of the outburst. Based on these considerations and guided by the SEDs, we propose a model in which turbulent plasma crosses a conical standing shock in the parsec-scale region of the jet. Here, the high-energy emission in the model is produced by inverse Compton scattering of seed photons supplied by either nonthermal radiation from a Mach disk, thermal emission from hot dust, or (for X-rays) synchrotron radiation from plasma that crosses the standing shock. For the two dates on which we fitted the model SED to the data, the model corresponds very well to the observations at all bands except at X-ray energies, where the spectrum is flatter than observed.

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.

We have observed 3C~279 in a gamma-ray flaring state in November 2008. We construct quasi-simultaneous spectral energy distributions (SEDs) of the source for the flaring period of 2008 and during a quiescent period in May 2010. Data have been compiled from observations with Fermi, Swift, RXTE, the VLBA, and various ground-based optical and radio telescopes. The objective is to comprehend the correspondence between the flux and polarization variations observed during these two time periods by carrying out a detailed spectral analyses of 3C~279 in the internal shock scenario, and gain insights into the role of intrinsic parameters and interplay of synchrotron and inverse Compton radiation processes responsible for the two states. As a first step, we have used a multi-slice time-dependent leptonic jet model, in the framework of the internal shock scenario, with radiation feedback to simulate the SED of 3C~279 observed in an optical high state in early 2006. We have used physical jet parameters obtained from the VLBA monitoring to guide our modeling efforts. We briefly discuss the effects of intrinsic parameters and various radiation processes in producing the resultant SED.

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 a comparison of the parameters of accretion disc outflows and the jet of the broad-line radio galaxy 3C 111 on sub-pc scales. We make use of published X-ray observations of ultra-fast outflows (UFOs) and new 43GHz VLBA images to track the jet knots ejection. We find that the superluminal jet coexists with the mildly relativistic outflows on sub-pc scales, possibly indicating a transverse stratification of a global flow. The two are roughly in pressure equilibrium, with the UFOs potentially providing additional support for the initial jet collimation. The UFOs are much more massive than the jet, but their kinetic power is probably about an order of magnitude lower, at least for the observations considered here. However, their momentum flux is equivalent and both of them are powerful enough to exert a concurrent feedback impact on the surrounding environment. A link between these components is naturally predicted in the context of MHD models for jet/outflow formation. However, given the high radiation throughput of AGNs, radiation pressure should also be taken into account. From the comparison with the long-term 2-10keV RXTE light curve we find that the UFOs are preferentially detected during periods of increasing flux. We also find the possibility to place the UFOs within the known X-ray dips-jet ejection cycles, which has been shown to be a strong proof of the disc-jet connection, in analogue with stellar-mass black holes. However, given the limited number of observations presently available, these relations are only tentative and additional spectral monitoring is needed to test them conclusively.

We present a preliminary analysis of new high resolution radio observations of the nearby TeV blazar Markarian 421 (z=0.031). This study is part of an ambitious multifrequency campaign, with observations in sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and gamma rays (Fermi-LAT, MAGIC, VERITAS). In this manuscript we consider only data obtained with the Very Long Baseline Array (VLBA) at seven epochs (one observation per month from January to July 2011) at 15 and 23.8 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch. We identified 5-6 components which are consistent with being stationary during the 6-month period reported here. The aim is to try to shed light on questions such as the nature of radiating particles, the connection between radio and gamma-ray emission, the location of emitting regions and the origin of the flux variability.

We analyze total and polarized intensity images of the quasar 3C 454.3 obtained monthly with the VLBA at 43 GHz within the ongoing Boston U. monitoring program of gamma-ray blazars started in June 2007. The data are supplemented by VLBA observations performed during intense campaigns of 2 week duration when the quasar was observed 3 times per campaign. We find a strong increase of activity in the parsec-scale jet of the quasar during high gamma-ray states in December 2009, April 2010, and November 2010. We detect new superluminal knots, K09 and K10, associated with the autumn 2009 and 2010 outbursts, respectively, and compare their kinematic parameters. We analyze optical polarimetric behavior along with polarization parameters of the parsec-scale jet and outline similarities and differences in polarization properties across wavelengths. The results of the analysis support the conclusions that the optical polarized emission is produced in a region located in the vicinity of the mm-wave core of the jet of the quasar, and that the gamma-ray outbursts occur when a superluminal disturbance passes through the core.

Analysis of comprehensive monitoring of 34 gamma-ray bright quasars, BL Lac objects, and radio galaxies reveals a close connection between events in the millimeter-wave emission imaged with the VLBA at 43 GHz and flares at gamma-ray and lower frequencies. Roughly 2/3 of the flares are coincident with the appearance of a new superluminal knot and/or a flare in the millimeter-wave "core'" located parsecs from the central engine. This presents a theoretical challenge to explain how the gamma-ray flux can often be variable on intra-day time-scales. Possible answers to this include very narrow opening angles of the jet, small volume filling factors of the highest energy electrons, chaotic magnetic fields, and turbulent velocity fields relative to the mean jet flow.

We present the results from an ultra-high-resolution 7mm Very Long Baseline Array (VLBA) study of the relativistic jet in the BL Lacertae object OJ287 from 1995 to 2011 containing 136 total intensity images. Analysis of the image sequence reveals a sharp jet-position-angle swing by >100 deg. during [2004,2006], as viewed in the plane of the sky, that we interpret as the crossing of the jet from one side of the line of sight to the other during a softer and longer term swing of the inner jet. Modulating such long term swing, our images also show for the first time a prominent erratic wobbling behavior of the innermost ~0.4mas of the jet with fluctuations in position angle of up to ~40 deg. over time scales ~2yr. This is accompanied by highly superluminal motions along non-radial trajectories, which reflect the remarkable non-ballistic nature of the jet plasma on these scales. The erratic nature and short time scales of the observed behavior rules out scenarios such as binary black hole systems, accretion disk precession, and interaction with the ambient medium as possible origins of the phenomenon on the scales probed by our observations, although such processes may cause longer-term modulation of the jet direction. We propose that variable asymmetric injection of the jet flow; perhaps related to turbulence in the accretion disk; coupled with hydrodynamic instabilities, leads to the non-ballistic dynamics that cause the observed non-periodic changes in the direction of the inner jet.

We locate the gamma-ray and lower frequency emission in flares of the BL Lac object AO 0235+164 at >12pc in the jet of the source from the central engine. We employ time-dependent multi-spectral-range flux and linear polarization monitoring observations, as well as ultra-high resolution (~0.15 milliarcsecond) imaging of the jet structure at lambda=7mm. The time coincidence in the end of 2008 of the propagation of the brightest superluminal feature detected in AO 0235+164 (Qs) with an extreme multi-spectral-range (gamma-ray to radio) outburst, and an extremely high optical and 7mm (for Qs) polarization degree provides strong evidence supporting that all these events are related. This is confirmed at high significance by probability arguments and Monte-Carlo simulations. These simulations show the unambiguous correlation of the gamma-ray flaring state in the end of 2008 with those in the optical, millimeter, and radio regime, as well as the connection of a prominent X-ray flare in October 2008, and of a series of optical linear polarization peaks, with the set of events in the end of 2008. The observations are interpreted as the propagation of an extended moving perturbation through a re-collimation structure at the end of the jet's acceleration and collimation zone.

We introduce MAPCAT, a long-term observing program for "Monitoring of AGN with Polarimetry at the Calar Alto Telescopes". Multi-spectral-range studies are critical to understand some of the most relevant current problems of high energy astrophysics of blazars such as their high energy emission mechanisms and the location of their gamma-ray emission region through event associations across the spectrum. Adding multi-spectral-range polarimetry allows for even more reliable identification of polarized flares across the spectrum in these kind of objects, as well as for more accurate modeling of their magnetic field. As part of a major international effort to study the long term multi-spectral range polarimetric behavior of blazars, MAPCAT uses -since mid 2007- CAFOS on the 2.2m Telescope at the Calar Alto Observatory (Almeria, Spain) to obtain monthly optical (R-band) photo-polarimetric measurements of a sample of 34 of the brightest gamma-ray, optical, and radio-millimeter blazars accessible from the northern hemisphere.

Aims: 3C345 was recently identified as a gamma-ray emitter, based on the first 20 months of Fermi-LAT data and optical monitoring. In this paper, a connection between the gamma-ray and optical variability of 3C345 and properties of its parsec-scale radio emission is investigated. Methods: The Fermi-LAT data of 3C345, covering an energy range of 0.1-300 GeV, were combined with 32 Very Long Baseline Array observations of the object made at 43.2 GHz in the period of January 2008 - March 2010. Results: The VLBA data reveal morphology and kinematics of the flow on scales of up to ~5 milliarcseconds (mas; deprojected linear distances of 380 parsecs). The brightness temperature, T_b(r), measured along the jet first decreases with distance proportional to r^-(0.95 +/-0.69) and later exhibits a break at ~0.3 mas, with T_b(r) proportional to r^-(4.11 +/-0.85) at larger separations. Variations of the gamma-ray, optical and parsec-scale radio emission show a similar long-term trend persistent during the entire VLBA monitoring period. The gamma-ray and optical variations on shorter time scales are related to structural changes in the jet on scales of ~0.3 mas (~23 parsecs, deprojected), with the gamma-ray and optical flares possibly related to the evolution of four distinct superluminal components identified in the flow. Conclusions: The observations indicate that both the quiescent and flaring components of the gamma-ray emission are produced in a region of the jet of ~23 pc in extent. This region may mark the Compton-loss dominated zone of the flow and its large extent may favor the synchrotron self-Compton mechanism for gamma-ray production in the relativistic jet of the quasar 3C345.

We present multiwavelength spectral analyses of two Fermi-LAT blazars, OJ 287 and 3C 279, that are part of the Boston University multiwaveband polarization program. The data have been compiled from observations with Fermi, RXTE, the VLBA, and various ground-based optical and radio telescopes. We simulate the dynamic spectral energy distributions (SEDs) within the framework of a multi-slice, time-dependent leptonic jet model for blazars, with radiation feedback, in the internal shock scenario. We use the physical jet parameters obtained from the VLBA monitoring to guide our modeling efforts. We discuss the role of intrinsic parameters and the interplay between synchrotron and inverse Compton radiation processes responsible for producing the resultant SEDs.

We compare the \gamma-ray light curves of the blazars, constructed with data provided by the Fermi Large Area Telescope, with flux and polarization variations in the VLBI core and bright superluminal knots obtained via monthly monitoring with the Very Long Baseline Array at 43 GHz. For all blazars in the sample that exhibit a high \gamma-ray state on time scales from several weeks to several months, an increase of the total flux in the mm-wave core is contemporaneous with the \gamma-ray activity (more than a third of the sample). Here we present the results for quasars with the most extreme \gamma-ray behavior (3C 454.3, 3C 273, 3C 279, 1222+216, and 1633+382). The sources show that in addition to the total flux intensity behavior, a maximum in the degree of polarization in the core or bright superluminal knot nearest to the core coincides with the time of a \gamma-ray peak to within the accuracy of the sampling of the radio data. These argue in favor of location of many of \gamma-ray outbursts in blazars outside of the broad line region, either in the vicinity or downstream of the mm-wave VLBI core.

We combine the Fermi-LAT light curve of the BL Lacertae type blazar OJ287 with time-dependent multi-waveband flux and linear polarization observations and submilliarcsecond-scale polarimetric images at lambda=7mm to locate the gamma-ray emission in prominent flares in the jet of the source >14pc from the central engine. We demonstrate a highly significant correlation between the strongest gamma-ray and millimeter-wave flares through Monte Carlo simulations. The two reported gamma-ray peaks occurred near the beginning of two major millimeter-wave outbursts, each of which is associated with a linear polarization maximum at millimeter wavelengths. Our very long baseline array observations indicate that the two millimeter-wave flares originated in the second of two features in the jet that are separated by >14pc. The simultaneity of the peak of the higher-amplitude gamma-ray flare and the maximum in polarization of the second jet feature implies that the gamma-ray and millimeter-wave flares are cospatial and occur >14pc from the central engine. We also associate two optical flares, accompanied by sharp polarization peaks, with the two gamma-ray events. The multi-waveband behavior is most easily explained if the gamma-rays arise from synchrotron self Compton scattering of optical photons from the flares. We propose that flares are triggered by interaction of moving plasma blobs with a standing shock.

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 present optical photopolarimetric observations of the BL Lac object S4 0954+658 obtained with the 70-cm telescope in Crimea, 40-cm telescope in St.Petersburg, and 1.8-m Perkins telescope at Lowell Observatory (Flagstaff, Az). After a faint state with a brightness level R ~17.6 mag registered in the first half of January 2011, the optical brightness of the source started to rise and reached ~14.8 mag during the middle of March, showing flare-like behavior. The most spectacular case of intranight variability was observed during the night of 2011 March 9, when the blazar brightened by ~0.7 mag within ~7 hours. During the rise of the flux the position angle of optical polarization rotated smoothly over more than 200 degrees. S4 0954+658 is a gamma-ray blazar with gamma-ray flux of (5{\pm}3)x10^{-10} phot/cm^2/s according to the Fermi 11-month Catalog Extragalactic Sources. Our analysis of contemporaneous Fermi LAT data does not show any sign of increased gamma-ray activity above the detection threshold except for an elevated flux on 2011 March 5, JD2455626, coincident with the local optical maximum.

We report on the location of the gamma-ray emission region in flares of the BL Lacertae object OJ287 at >14pc from the central supermassive black hole. We employ data from multi-spectral range (total flux and linear polarization) monitoring programs combined with sequences of ultra-high-resolution 7mm VLBA images. The correlation between the brightest gamma-ray and mm flares is found to be statistically significant. The two gamma-ray peaks, detected by Fermi-LAT, that we report here happened at the rising phase of two exceptionally bright mm flares accompanied by sharp linear polarization peaks. The VLBA images show that these mm flares in total flux and polarization degree occurred in a jet region at >14pc from the innermost jet region. The time coincidence of the brighter gamma-ray flare and its corresponding mm linear polarization peak evidences that both the gamma-ray and mm outbursts occur >14pc from the central black hole. We find two sharp optical flares occurring at the peak times of the two reported gamma-ray flares. This is interpreted as the gamma-ray flares being produced by synchrotron self-Compton scattering of optical photons from the flares triggered by the interaction of moving knots with a stationary conical shock in the jet.

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at gamma-ray energies <=20 MeV. Recently a bright gamma-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640+396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 - April 2010) to associate and identify the dominating gamma-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this gamma-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing gamma-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No gamma-ray emission from B3 1640+396 was detected.

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.

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 source of Faraday rotation in the jet of the radio galaxy 3C120 is analyzed through Very Long Baseline Array observations carried out between 1999 and 2007 at 86, 43, 22, 15, 12, 8, 5, 2, and 1.7 GHz. Comparison of observations from 1999 to 2001 reveals uncorrelated changes in the linear polarization of the underlying jet emission and the Faraday rotation screen: while the rotation measure (RM) remains constant between approximately 2 and 5 mas from the core, the RM-corrected electric vector position angles (EVPAs) of two superluminal components are rotated by almost 90 degrees when compared to other components moving through similar jet locations. On the other hand, the innermost 2 mas experiences a significant change in RM -- including a sign reversal -- but without variations in the RM-corrected EVPAs. Similarly, observations in 2007 reveal a double sign reversal in RM along the jet, while the RM-corrected EVPAs remain perpendicular to the jet axis. Although the observed coherent structure and gradient of the RM along the jet supports the idea that the Faraday rotation is produced by a sheath of thermal electrons that surrounds the emitting jet, the uncorrelated changes in the RM and RM-corrected EVPAs indicate that the emitting jet and the source of Faraday rotation are not closely connected physically and have different configurations for the magnetic field and/or kinematical properties. Furthermore, the existence of a region of enhanced RM whose properties remain constant over three years requires a localized source of Faraday rotation, favoring a model in which a significant fraction of the RM originates in foreground clouds.

A possible source of $\gamma$-ray photons observed from the jets of blazars is inverse Compton scattering by relativistic electrons of infrared seed photons from a hot, dusty torus in the nucleus. We use observations from the Spitzer Space Telescope to search for signatures of such dust in the infrared spectra of four $\gamma$-ray bright blazars, the quasars 4C 21.35, CTA102, and PKS 1510$-$089, and the BL Lacertae object ON231. The spectral energy distribution (SED) of 4C 21.35 contains a prominent infrared excess indicative of dust emission. After subtracting a non-thermal component with a power-law spectrum, we fit a dust model to the residual SED. The model consists of a blackbody with temperature $\sim1200$ K, plus a much weaker optically thin component at $\sim660$ K. The total luminosity of the thermal dust emission is $7.9\pm0.2 \times 10^{45}$ erg s$^{-1}$. If the dust lies in an equatorial torus, the density of IR photons from the torus is sufficient to explain the $\gamma$-ray flux from 4C 21.35 as long as the scattering occurs within a few parsecs of the central engine. We also report a tentative detection of dust in the quasar CTA102, in which the luminosity of the infrared excess is $7 \pm 2 \times 10^{45}$ erg s$^{-1}$. However, in CTA102 the far-IR spectra are too noisy to detect the $10 \mu$m silicate feature. Upper limits to the luminosity from thermal emission from dust in PKS 1510-089, and ON231, are, $2.3\times10^{45}$, and $6.6\times10^{43}$ erg s$^{-1}$, respectively. These upper limits do not rule out the possibility of inverse Compton up-scattering of IR photons to $\gamma$-ray energies in these two sources. The estimated covering factor of the hot dust in 4C 21.35, 22%, is similar to that of non-blazar quasars; however, 4C 21.35 is deficient in cooler dust.

We are leading a comprehensive multi-waveband monitoring program of 34 gamma-ray bright blazars designed to locate the emission regions of blazars from radio to gamma-ray frequencies. The "maps" are anchored by sequences of images in both total and polarized intensity obtained with the VLBA at an angular resolution of ~ 0.1 milliarcseconds. The time-variable linear polarization at radio to optical wavelengths and radio to gamma-ray light curves allow us to specify the locations of flares relative to bright stationary features seen in the images and to infer the geometry of the magnetic field in different regions of the jet. Our data reveal that some flares occur simultaneously at different wavebands and others are only seen at some of the frequencies. The flares are often triggered by a superluminal knot passing through the stationary "core" on the VLBA images. Other flares occur upstream or even parsecs downstream of the core.

We present an X-ray image of the BL Lacertae object OJ287 revealing a long jet, curved by 55 degrees and extending 20", or 90 kpc from the nucleus. This de-projects to >1 Mpc based on the viewing angle on parsec scales. Radio emission follows the general X-ray morphology but extends even farther from the nucleus. The upper limit to the isotropic radio luminosity, ~2E24 W/Hz, places the source in the Fanaroff-Riley 1 (FR 1) class, as expected for BL Lac objects. The spectral energy distribution indicates that the extended X-ray emission is from inverse Compton scattering of cosmic microwave background photons. In this case, the derived magnetic field is B ~ 5 microGauss, the minimum electron energy is 7-40 m_e c^2, and the Doppler factor is delta ~ 8 in a knot 8" from the nucleus. The minimum total kinetic power of the jet is 1-2E45 erg/s. Upstream of the bend, the width of the X-ray emission in the jet is about half the projected distance from the nucleus. This implies that the highly relativistic bulk motion is not limited to an extremely thin spine, as has been proposed previously for FR 1 sources. The bending of the jet, the deceleration of the flow from parsec to kiloparsec scales, and the knotty structure can all be caused by standing shocks inclined by ~7 degrees to the jet axis. Moving shocks resulting from major changes in the flow properties can also reproduce the knotty structure, but such a model does not explain as many of the observational details.

3C 345 is one of the archetypical active galactic nuclei, showing structural and flux variability on parsec scales near a compact unresolved radio core. During the last 2 years, the source has been undergoing a period of high activity visible in the broad spectral range, from radio through high-energy bands. We have been monitoring parsec-scale radio emission in 3C 345 during this period at monthly intervals, using the VLBA at 15, 24, and 43 GHz. Our radio observations are compared with gamma-ray emission detected by Fermi-LAT in the region including 3C 345 (1FGL J1642.5+3947). Three distinct gamma-ray events observed in this region are associated with the propagation of relativistic plasma condensations inside the radio jet of 3C 345. We report on evidence for the gamma-rays to be produced in a region of the jet of up to 40 pc (de-projected) in extent. This suggests the synchrotron self-Compton process as the most likely mechanism for production of gamma-rays in the source.

We combine time-dependent multi-waveband flux and linear polarization observations with sub-milliarcsecond-scale polarimetric images at lambda=7mm of the BL Lacertae-type blazar OJ287 to locate the gamma-ray emission in prominent flares in the jet of the source >14pc from the central engine. We demonstrate a highly significant correlation between the strongest gamma-ray and millimeter-wave flares through Monte-Carlo simulations. The two reported gamma-ray peaks occurred near the beginning of two major mm-wave outbursts, each of which is associated with a linear polarization maximum at millimeter wavelengths. Our Very Long Baseline Array observations indicate that the two mm-wave flares originated in the second of two features in the jet that are separated by >14 pc. The simultaneity of the peak of the higher-amplitude gamma-ray flare and the maximum in polarization of the second jet feature implies that the gamma-ray and mm-wave flares are co-spatial and occur >14 pc from the central engine. We also associate two optical flares, accompanied by sharp polarization peaks, with the two gamma-ray events. The multi-waveband behavior is most easily explained if the gamma-rays arise from synchrotron self-Compton scattering of optical photons from the flares. We propose that flares are triggered by interaction of moving plasma blobs with a standing shock. The gamma-ray and optical emission is quenched by inverse Compton losses as synchrotron photons from the newly shocked plasma cross the emission region. The mm-wave polarization is high at the onset of a flare, but decreases as the electrons emitting at these wavelengths penetrate less polarized regions.

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.

Comprehensive VLBI and multi-waveband monitoring indicate that a single superluminal knot can cause a number of gamma-ray flares at different locations. However, the often very rapid variability timescale is a challenge to theoretical models when a given flare (perhaps the majority of those observed) is inferred from observations to lie near the 43 GHz core, parsecs from the central engine. We present some relevant observational results, using the BL Lac object AO 0235+164 as an example. We propose a turbulent cell model leading to a frequency-dependent filling factor of the emission region. This feature of the model can provide a solution to the timescale dilemma and other characteristics of blazar emission.

We analyze the behavior of the parsec-scale jet of the quasar 3C~454.3 during pronounced flaring activity in 2005-2008. Three major disturbances propagated down the jet along different trajectories with Lorentz factors $\Gamma>$10. The disturbances show a clear connection with millimeter-wave outbursts, in 2005 May/June, 2007 July, and 2007 December. High-amplitude optical events in the $R$-band light curve precede peaks of the millimeter-wave outbursts by 15-50 days. Each optical outburst is accompanied by an increase in X-ray activity. We associate the optical outbursts with propagation of the superluminal knots and derive the location of sites of energy dissipation in the form of radiation. The most prominent and long-lasting of these, in 2005 May, occurred closer to the black hole, while the outbursts with a shorter duration in 2005 Autumn and in 2007 might be connected with the passage of a disturbance through the millimeter-wave core of the jet. The optical outbursts, which coincide with the passage of superluminal radio knots through the core, are accompanied by systematic rotation of the position angle of optical linear polarization. Such rotation appears to be a common feature during the early stages of flares in blazars. We find correlations between optical variations and those at X-ray and $\gamma$-ray energies. We conclude that the emergence of a superluminal knot from the core yields a series of optical and high-energy outbursts, and that the mm-wave core lies at the end of the jet's acceleration and collimation zone.

We present 1.7, 5, 15, 22 and 43 GHz polarimetric multi-epoch VLBA observations of the radio galaxy 3C 120. The higher frequency observations reveal a new component, not visible before April 2007, located 80 mas from the core (which corresponds to a deprojected distance of 140 pc), with a brightness temperature about 600 times higher than expected at such distances. This component (hereafter C80) is observed to remain stationary and to undergo small changes in its brightness temperature during more than two years of observations. A helical shocked jet model - and perhaps some flow acceleration - may explain the unusually high Tb of C80, but it seems unlikely that this corresponds to the usual shock that emerges from the core and travels downstream to the location of C80. It appears that some other intrinsic process in the jet, capable of providing a local burst in particle and/or magnetic field energy, may be responsible for the enhanced brightness temperature observed in C80, its sudden appearance in April 2007, and apparent stationarity.

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.

As part of our comprehensive long-term multi-waveband monitoring of 34 blazars, we followed the activity in the jet of the blazar PKS 1510-089 during major outbursts during the first half of 2009. The most revealing event was a two-month long outburst that featured a number of gamma-ray flares. During the outburst, the position angle of optical linear polarization rotated by about 720 degrees, which implies that a single emission feature was responsible for all of the flares during the outburst. At the end of the rotation, a new superluminal knot (~ 22c) passed through the "core" seen on 43 GHz VLBA images at essentially the same time as an extremely sharp, high-amplitude gamma-ray and optical flare occurred. We associate the entire multi-flare outburst with this knot. The ratio of gamma-ray to synchrotron integrated flux indicates that some of the gamma-ray flares resulted from inverse Compton scattering of seed photons outside the ultra-fast spine of the jet. Because many of the flares occurred over time scales of days or even hours, there must be a number of sources of IR-optical-UV seed photons -- probably synchrotron emission -- surrounding the spine, perhaps in a slower sheath of the jet.

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.

We perform monthly total and polarized intensity imaging of a sample of $\gamma$-ray blazars (33 sources) with the Very Long Baseline Array (VLBA) at 43 GHz with the high resolution of 0.1 milliarcseconds. From Summer 2008 to October 2009 several of these blazars triggered Astronomical Telegrams due to a high $\gamma$-ray state detected by the Fermi Large Area Telescope (LAT): AO 0235+164, 3C 273, 3C 279, PKS 1510-089, and 3C 454.3. We have found that 1) $\gamma$-ray flares in these blazars occur during an increase of the flux in the 43 GHz VLBI core; 2) strong $\gamma$-ray activity, consisting of several flares of various amplitudes and durations (weeks to months), is simultaneous with the propagation of a superluminal knot in the inner jet, as found previously for BL Lac (Marscher et al. 2008); 3) coincidence of a superluminal knot with the 43 GHz core precedes the most intense $\gamma$-ray flare by 36$\pm$24 days. Our results strongly support the idea that the most dramatic $\gamma$-ray outbursts of blazars originate in the vicinity of the mm-wave core of the relativistic jet. These results are preliminary and should be tested by future monitoring with the VLBA and Fermi.

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 explore the variability and cross-frequency correlation of the flux density and polarization of the blazar OJ287, using imaging at 43 GHz with the Very Long Baseline Array, as well as optical and near-infrared polarimetry. The polarization and flux density in both the optical waveband and the 43 GHz compact core increased by a small amount in late 2005, and increased significantly along with the near-IR polarization and flux density over the course of 10 days in early 2006. Furthermore, the values of the electric vector position angle (EVPA) at the three wavebands are similar. At 43 GHz, the EVPA of the blazar core is perpendicular to the flow of the jet, while the EVPAs of emerging superluminal knots are aligned parallel to the jet axis. The core polarization is that expected if shear aligns the magnetic field at the boundary between flows of disparate velocities within the jet. Using variations in flux density, percentage polarization, and EVPA, we model the inner jet as a spine-sheath system. The model jet contains a turbulent spine of half-width 1.2 degrees and maximum Lorentz factor of 16.5, a turbulent sheath with Lorentz factor of 5, and a boundary region of sheared field between the spine and sheath. Transverse shocks propagating along the fast, turbulent spine can explain the superluminal knots. The observed flux density and polarization variations are then compatible with changes in the direction of the inner jet caused by a temporary change in the position of the core if the spine contains wiggles owing to an instability. In addition, we can explain a stable offset of optical and near-IR percentage polarization by a steepening of spectral index with frequency, as supported by the data.

We present some results of the deepest extragalactic survey performed by the INTEGRAL satellite. The fraction of very absorbed AGN is quite large. The sharp decrease in the absorption fraction with X-ray luminosity observed at lower-energy X-rays is not observed. The current lack of truly Compton-thick objects, with an upper limit of 14% to the size of this population, is just compatible with recent modeling of the cosmic X-ray background. We also study the prospects for a future hard X-ray serendipitous survey with Simbol-X. We show that Simbol-X will easily detect a large number of serendipitous AGN, allowing us to study the evolution of AGN up to redshifts about 2, opening the door to the cosmological study of hard X-ray selected AGN, which is barely possible with existing satellites like Swift and INTEGRAL.

We present results from a multi-year monitoring campaign of the broad-line radio galaxy 3C 120, using the Rossi X-ray Timing Explorer (RXTE) for nearly five years of observations. Additionally, we present coincident optical monitoring using data from several ground-based observatories. Both the X-ray and optical emission are highly variable and appear to be strongly correlated, with the X-ray emission leading the optical by 28 days. The X-ray power density spectrum is best fit by a broken power law, with a low-frequency slope of -1.2, breaking to a high-frequency slope of -2.1, and a break frequency of log nu_b=-5.75 Hz, or 6.5 days. This value agrees well with the value expected based on 3C 120's mass and accretion rate. We find no evidence for a second break in the power spectrum. Combined with a moderately soft X-ray spectrum (Gamma=1.8) and a moderately high accretion rate (mdot / mdot_Edd ~ 0.3), this indicates that 3C 120 fits in wellwith the high/soft variability state found in most other AGNs. Previous studies have shown that the spectrum has a strong Fe K-alpha line, which may be relativistically broadened. The presence of this line, combined with a power spectrum similar to that seen in Seyfert galaxies, suggests that the majority of the X-ray emission in this object arises in or near the disk, and not in the jet.

We present radio-to-optical data taken by the WEBT, supplemented by VLBA and RXTE observations, of 3C 279. Our goal is to use this extensive database to draw inferences regarding the physics of the relativistic jet. We assemble multifrequency light curves with data from 30 ground-based observatories and the space-based instruments, along with linear polarization vs. time in the optical R band. In addition, we present a sequence of 22 images (with polarization vectors) at 43 GHz at resolution 0.15 milliarcsec, obtained with the VLBA. We analyse the light curves and polarization, as well as the spectral energy distributions at different epochs, corresponding to different brightness states. The IR-optical-UV continuum spectrum of the variable component corresponds to a power law with a constant slope of -1.6, while in the 2.4-10 keV X-ray band it varies in slope from -1.1 to -1.6. The steepest X-ray spectrum occurs at a flux minimum. During a decline in flux from maximum in late 2006, the optical and 43 GHz core polarization vectors rotate by ~300 degrees. The continuum spectrum agrees with steady injection of relativistic electrons with a power-law energy distribution of slope -3.2 that is steepened to -4.2 at high energies by radiative losses. The X-ray emission at flux minimum comes most likely from a new component that starts in an upstream section of the jet where inverse Compton scattering of seed photons from outside the jet is important. The rotation of the polarization vector implies that the jet contains a helical magnetic field that extends ~20 pc past the 43 GHz core.

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.

Very long baseline interferometric observations of the radio galaxy 3C 120 show a systematic presence of gradients in Faraday rotation and degree of polarization across and along the jet. These are revealed by the passage of multiple superluminal components throughout the jet as they move out from the core in a sequence of 12 monthly polarimetric observations taken with the VLBA at 15, 22, and 43 GHz. The degree of polarization has an asymmetric profile in which the northern side of the jet is more highly polarized. The Faraday rotation measure is also stratified across the jet width, with larger values for the southern side. Superposed on this structure we find a localized region of high Faraday rotation measure (about 6000 rad/m^2) between approximately 3 and 4 mas from the core. This region of enhanced Faraday rotation may result from the interaction of the jet with the ambient medium, which may also explain the stratification in degree of polarization. The data are also consistent with a helical magnetic field in a two-fluid jet model, consisting of an inner emitting jet and a sheath of nonrelativistic electrons.

We present a sequence of 12 monthly polarimetric 15, 22, and 43 GHz VLBA observations of the radio galaxy 3C 120 revealing a systematic presence of gradients in Faraday rotation and degree of polarization across and along the jet. The degree of polarization increases with distance from the core and toward the jet edges, and has an asymmetric profile in which the northern side of the jet is more highly polarized. The Faraday rotation measure is also stratified across the jet width, with larger values for the southern side. We find a localized region of high Faraday rotation measure superposed on this structure between approximately 3 and 4 mas from the core, with a peak of about 6000 rad/m^2. Interaction of the jet with the external medium or a cloud would explain the confined region of enhanced Faraday rotation, as well as the stratification in degree of polarization and the flaring of superluminal knots when crossing this region. The data are also consistent with a helical field in a two-fluid jet model, consisting of an inner, emitting jet and a sheath containing nonrelativistic electrons. However, this helical magnetic field model cannot by itself explain the localized region of enhanced Faraday rotation. The polarization electric vectors, predominantly perpendicular to the jet axis once corrected for Faraday rotation, require a dominant component parallel to the jet axis (in the frame of the emitting plasma) for the magnetic field in the emitting region.

We present the results of polarimetric ($R$ band) and multicolor photometric ($BVRIJHK$) observations of the blazar AO 0235+16 during an outburst in 2006 December. The data reveal a short timescale of variability (several hours), which increases from optical to near-IR wavelengths; even shorter variations are detected in polarization. The flux density correlates with the degree of polarization, and at maximum degree of polarization the electric vector tends to align with the parsec-scale jet direction. We find that a variable component with a steady power-law spectral energy distribution and very high optical polarization (30-50%) is responsible for the variability. We interpret these properties of the blazar withina model of a transverse shock propagating down the jet. In this case a small change in the viewing angle of the jet, by $\lesssim 1^o$, and a decrease in the shocked plasma compression by a factor of $\sim$1.5 are sufficient to account for the variability.

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

An 11-day monitoring campaign in late 2005 reveals clear correlation in polarization between the optical emission and the region of the intensity peak (the "pseudocore") at the upstream end of the jet in 43 GHz Very Long Baseline Array images in the highly variable quasar PKS 0420-014. The electric-vector position angle (EVPA) of the pseudocore rotated by about 80 degrees in four VLBA observations over a period of nine days, matching the trend of the optical EVPA. In addition, the 43 GHz EVPAs agree well with the optical values when we correct the former for Faraday rotation. Fluctuations in the polarization at both wavebands are consistent with the variable emission arising from a standing conical shock wave that compresses magnetically turbulent plasma in the ambient jet. The volume of the variable component is the same at both wavebands, although only ~20% of the total 43 GHz emission arises from this site. The remainder of the 43 GHz flux density must originate in a separate region with very low polarization. If 0420-014 is a typical case, the nonthermal optical emission from blazars originates primarily in and near the pseudocore rather than closer to the central engine where the flow collimates and accelerates.

We investigate the variations of the magnetic field, Doppler factor, and relativistic particle density along the jet of a quasar at z=0.72. We chose 4C 19.44 for this study because of its length and straight morphology. The 18 arcsec length of the jet provides many independent resolution elements in the Chandra X-ray image. The straightness suggests that geometry factors, although uncertain, are almost constant along the jet. We assume the X-ray emission is from inverse Compton scattering of the cosmic microwave background. With the aid of assumptions about jet alignment, equipartition between magnetic-field and relativistic-particle energy, and filling factors, we find that the jet is in bulk relativistic motion with a Doppler factor about 6 at an angle no more than 10 degrees to the line of sight over deprojected distances about 150--600 kpc from the quasar, and with a magnetic field approximately 10 micro Gauss.

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.

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

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

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

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

Intensive VLBI monitoring programs of jets in AGN are showing the existence of intricate emission patterns, such as upstream motions or slow moving and quasi-stationary componentes trailing superluminal features. Relativistic hydrodynamic and emission simulations of jets are in very good agreement with these observations, proving as a powerful tool for the understanding of the physical processes taking place in the jets of AGN, microquasars and GRBs. These simulations show that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or ambient medium. Both observations and simulations suggest that shock-in-jet models may be an overly simplistic idealization when interpreting the emission structure observed in actual jets.

We present Chandra X-ray (0.2-8 keV) and Very Large Array radio (15 and 5 GHz) images of the $\gamma$-ray bright, superluminal quasar 0827+243. The X-ray jet bends sharply--by ~90 deg., presumably amplified by projection effects--5'' from the core. Only extremely weak radio emission is detected between the nuclear region and the bend. The X-ray continuum spectrum of the combined emission of the knots is rather flat, with a slope of $-0.4\pm 0.2$, while the 5-15 GHz spectra are steeper for knots detected in the radio. These characteristics, as well as non-detection of the jet in the optical band by the Hubble Space Telescope, pose challenges to models for the spectral energy distributions (SEDs) of the jet features. The SEDs could arise from pure synchrotron emission from either a single or dual population of relativistic electrons only if the minimum electron energy per unit mass $\gamma_{\min} \gtrsim 1000$. Alternatively, the X-ray emission could result from inverse Compton scattering of the Cosmic Microwave Background photons by electrons with Lorentz factors as low as $\gamma \sim 15$. In all models, the bulk Lorentz factor of the jet flow $\Gamma\gtrsim 20$ found on parsec scales must continue without substantial deceleration out to 800 kpc (deprojected) from the nucleus. Deceleration does appear to occur at and beyond the sharp bend, such that the flow could be only mildly relativistic at the end of the jet.

We present a 16-month sequence of monthly polarimetric 43 GHz VLBA images of the radio galaxy 3C 120. The images probe the inner regions of the radio jet of this relatively nearby superluminal radio galaxy at a linear resolution of 0.07 $h_{65}^{-1}$ pc ($H_o= 65 h_{65}$ km s$^{-1}$ Mpc$^{-1}$). We follow the motion of a number of features with apparent velocities between 4.01$\pm$0.08 and $5.82\pm 0.13 h_{65}^{-1} c$. A new superluminal knot, moving at $4.29\pm 0.16 h_{65}^{-1} c$, is observed to be ejected from the core at a time coincident with the largest flare ever observed for this source at millimeter wavelengths. Changes in the position angle of this component, as well as a progressive rotation of its magnetic polarization vector, suggest the presence of a twisted (resembling a helix in projection) configuration of the underlying jet magnetic field and jet geometry. We identify several knots that appear in the wake of the new superluminal component, moving at proper motions $\sim 4$ times slower than any of the other moving knots observed in 3C 120. These features have properties similar to those of the ``trailing'' shocks seen in relativistic, time-dependent, hydrodynamical and emission simulations of compact jets. Such trailing compressions are triggered by pinch-mode jet-body instabilities caused by the propagation of a strong perturbation, which we associate with the new strong superluminal component.

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 the results of a program to monitor the structure of the radio emission in 42 $\gamma$-ray bright blazars (31 quasars and 11 BL Lac objects) with the VLBA at 43, 22, and occasionally 15 and 8.4 GHz, over the period from November 1993 to July 1997. We determine proper motions in 33 sources and find that the apparent superluminal motions in $\gamma$-ray sources are much faster than for the general population of bright compact radio sources. This follows the strong dependence of the $\gamma$-ray flux on the level of relativistic beaming for both external-radiation Compton and synchrotron self-Compton emission. There is a positive correlation (correlation coefficient $r$=0.45) between the flux density of the VLBI core and the $\gamma$-ray flux and a moderate correlation (partial correlation coefficient $r$=0.31) between $\gamma$-ray apparent luminosity and superluminal velocities of jet components, as expected if the $\gamma$-ray emission originates in a very compact region of the relativistic jet and is highly beamed. In 43% of the sources the jet bends by more than 20$^\circ$ on parsec scales, which is consistent with amplification by projection effects of modest actual changes in position angle. In 27 of the sources in the sample there is at least one non-core component that appears to be stationary during our observations. Different characteristics of stationary features close to and farther from the core lead us to suggest two different classes of stationary components: those within about 2 milliarcseconds (mas) of the core, probably associated with standing hydrodynamical compressions, and those farther down the jet, which tend to be associated with bends in the jet.