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 quasi-periodic variability found in two blazars included in the Steward Observatory Blazar Monitoring data sample: the BL Lac object 3C 66A and the Flat Spectrum Radio Quasar B2 1633+38. We collect optical photometric and polarimetric data in V and R bands of these sources from different observatories: St. Petersburg University, Crimean Astrophysical Observatory, WEBT-GASP, Catalina Real-Time Transient Survey, Steward Observatory, STELLA Robotic Observatory and Katzman Automatic Imaging Telescope. In addition, an analysis of the $\gamma$-ray light curves from $\textit{Fermi}$-LAT is included. Three methods are used to search for any periodic behaviour in the data: the Z-transform Discrete Correlation Function, the Lomb-Scargle periodogram and the Weighted Wavelet Z-transform. We find evidences of possible quasi-periodic variability in the optical photometric data of both sources with periods of $\sim$3 years for 3C 66A and $\sim$1.9 years for B2 1633+38, with significances between 3$\sigma$ and 5$\sigma$. Only B2 1633+38 shows evidence of this behaviour in the optical polarized data set at a confidence level of 2$\sigma$-4$\sigma$. This is the first reported evidence of quasi-periodic behaviour in the optical light curve of B2 1633+38. Also a hint of quasi-periodic behaviour is found in the $\gamma$-ray light curve of B2 1633+38 with a confidence level $\geqslant$2$\sigma$, while no periodicity is observed for 3C 66A in this energy range. We propose different jet emission models that could explain the quasi-periodic variability and the differences found between these two sources.

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 the results of flux density, spectral index, and polarization intra-night monitoring studies of a sample of eight optically bright blazars, carried out by employing several small to moderate aperture (0.4\,m to 1.5\,m diameter) telescopes fitted with CCDs and polarimeters located in Europe, India, and Japan. The duty cycle of flux variability for the targets is found to be $\sim 45$ percent, similar to that reported in earlier studies. The computed two-point spectral indices are found to be between 0.65 to 1.87 for our sample, comprised of low- and intermediate frequency peaked blazars, with one exception; they are also found to be statistically variable for about half the instances where `confirmed' variability is detected in flux density. In the analysis of the spectral evolution of the targets on hourly timescale, a counter-clockwise loop (soft-lagging) is noted in the flux-spectral index plane on two occasions, and in one case a clear spectral flattening with the decreasing flux is observed. In our data set, we also observe a variety of flux-polarization degree variability patterns, including instances with a relatively straightforward anti-correlation, correlation, or counter-clockwise looping. These changes are typically reflected in the flux-polarization angle plane: the anti-correlation between the flux and polarization degree is accompanied by an anti-correlation between the polarization angle and flux, while the counter-clockwise flux-PD looping behaviour is accompanied by a clockwise looping in the flux-polarization angle representation. We discuss our findings in the framework of the internal shock scenario for blazar sources.

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 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$.

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.

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.

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 UBVRI photometry of the supernova 2017eaw in NGC 6946, obtained in the period from May 14 until December 7, 2017. We derive dates and magnitudes of maximum light in the UBVRI bands and the parameters of the light curves. We discuss colour evolution, extinction and maximum luminosity of SN 2017eaw. Preliminary modeling is carried out, and the results are in satisfactory agreement with the light curves in the UBVRI bands.

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.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

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.

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

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.

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.

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.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

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 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.

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.

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 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.

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.

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

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

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

We present 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.

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 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.

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.

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.

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 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.