Context:Blazars are the rarest and most powerful active galactic nuclei, playing a crucial and growing role in today multi-frequency and multi-messenger astrophysics. Current blazar catalogs, however, are incomplete and particularly depleted at low Galactic latitudes. Aims: We aim at augmenting the current blazar census to build a catalog of blazar candidates with homogeneous sky coverage that can provide candidate counterparts to unassociated gamma-ray sources, sources of high-energy neutrino emission, and ultra-high energy cosmic rays. Methods: Starting from the ALMA Calibrator Catalog we built a catalog of 1580 blazar candidates (ALMA Blazar Candidates, ABC) for which we collect multi-wavelength information. We also compared ABC sources with existing blazar catalogs. Results: The ABC catalogue fills the lack of low Galactic latitude sources in current blazar catalogues. ABC sources are significantly dimmer than known blazars in Gaia g band, and they appear bluer in SDSS and WISE colors. The majority of ABC sources (~ 90%) have optical spectra that classify them as QSO, while the remaining sources resulted galactic objects. ABC sources are similar in X-rays to known blazar, while in gamma-rays they are on average dimmer and softer, indicating a significant contribution of FSRQ sources. Making use of WISE colours, we classified 715 ABC sources as candidate gamma-ray blazar of different classes. Conclusions: We built a new catalogue of 1580 candidate blazars with a rich multi-wavelength data-set, filling the lack of low Galactic latitude sources in current blazar catalogues. This will be particularly important to identify the source population of high energy neutrinos or ultra-high energy cosmic rays. The data collected by the upcoming LSST surveys will provide a key tool to investigate the possible blazar nature of these sources.

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

4C 71.07 is a high-redshift blazar whose optical radiation is dominated by quasar-like nuclear emission. We here present the results of a spectroscopic monitoring of the source to study its unbeamed properties. We obtained 24 optical spectra at the Nordic Optical Telescope (NOT) and William Herschel Telescope (WHT) and 3 near-infrared spectra at the Telescopio Nazionale Galileo (TNG). They show no evidence of narrow emission lines. The estimate of the systemic redshift from the H$\beta$ and H$\alpha$ broad emission lines leads to $z_{\rm sys}=2.2130 \pm 0.0004$. Notwithstanding the nearly face-on orientation of the accretion disc, the high-ionization emission lines present large broadening as well as noticeable blueshifts, which increase with the ionizing energy of the corresponding species. This is a clear indication of strong ionized outflows. Line broadening and blueshift appear correlated. We applied scaling relationships to estimate the mass of the supermassive black hole from the Balmer and C IV lines, taking into account the prescriptions to correct for outflow. They give $M_{\rm BH} \sim 2 \times 10^9 \, M_\odot$. We derived an Eddington luminosity $L_{\rm Edd} \sim 2.5 \times 10^{47} \rm \, erg \, s^{-1}$ $\sim L_{\rm disc}$, and a broad line region luminosity $L_{\rm BLR} \sim 1.5 \times 10^{46} \rm \, erg \, s^{-1}$. The line fluxes do not show significant variability in time. In particular, there is no line reaction to the jet flaring activity detected in 2015 October and November. This implies that the jet gives no contribution to the photoionization of the broad line region in the considered period.

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 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 BL Lac type object 1ES 2344+51.4 (redshift $z=0.044$) was one of the first sources to be included in the extreme high-peaked BL Lac (EHBL) family. EHBLs are characterised by a broadband spectral energy distribution (SED) featuring the synchrotron peak above $\sim 10^{17}$ Hz. From previous studies of 1ES 2344+51.4 in the very-high-energy (VHE, $>$100 GeV) gamma-ray range, its inverse Compton (IC) peak is expected around 200 GeV. 1ES 2344+51.4 was first detected in the VHE range by Whipple in 1995 during a very bright outburst showing around 60 % of the flux of the Crab Nebula above 350 GeV. In 1996, during another flare in the X-ray band, observations with Beppo-SAX revealed a large 0.1-10 keV flux variability on timescales of a few hours and an impressive frequency shift of the synchrotron peak to above $10^{18}$ Hz. Later on, this extreme behaviour of the source motivated several multiwavelength campaigns, during most of which the source appeared to be in a low state and showing no clear signs of "extremeness". In August 2016, FACT detected 1ES 2344+51.4 in a high state and triggered multiwavelength observations. The VHE observations show a flux level similar to the historical maximum of 1995. The combination of MAGIC, FACT, and Fermi-LAT spectra provides an unprecedented characterisation of the IC peak. It is the first time that simultaneous HE and VHE data are presented for this object during a flaring episode. We find an atypically hard spectrum in the VHE $\gamma$-rays as well as a hard X-ray spectrum, revealing a renewed extreme behaviour.

The object 4C 71.07 is a high-redshift blazar whose spectral energy distribution shows a prominent big blue bump and a strong Compton dominance. We present the results of a two-year multiwavelength campaign led by the Whole Earth Blazar Telescope (WEBT) to study both the quasar core and the beamed jet of this source. The WEBT data are complemented by ultraviolet and X-ray data from Swift, and by gamma-ray data by Fermi. The big blue bump is modelled by using optical and near-infrared mean spectra obtained during the campaign, together with optical and ultraviolet quasar templates. We give prescriptions to correct the source photometry in the various bands for the thermal contribution, in order to derive the non-thermal jet flux. The role of the intergalactic medium absorption is analysed in both the ultraviolet and X-ray bands. We provide opacity values to deabsorb ultraviolet data, and derive a best-guess value for the hydrogen column density through the analysis of X-ray spectra. We estimate the disc and jet bolometric luminosities, accretion rate, and black hole mass. Light curves do not show persistent correlations among flux changes at different frequencies. We study the polarimetric behaviour and find no correlation between polarisation degree and flux, even when correcting for the dilution effect of the big blue bump. Similarly, wide rotations of the electric vector polarisation angle do not seem to be connected with the source activity.

We present a unique method that allows the LSST to scan the sky for stellar variability on short timescales. The operational component of the strategy requires LSST to take star trail images. The image processing component uses deep learning to sift for transient events on timescales down to 10 ms. We advocate for enabling this observing mode with LSST, as coupling this capability with the LSST's tremendous 319.5 m$^2$deg$^2$ etendue will produce the first wide area optical survey of the universe on these timescales. We explain how these data will advance both planned lines of investigation and enable new research in the areas of stellar flares, cataclysmic variables, active galactic nuclei, Kuiper Belt objects, gamma-ray bursts, and fast radio bursts.

The aim of this white paper is to discuss the observing strategies for the LSST Wide-Fast-Deep that would improve the study of blazars (emission variability, census, environment) and Fast Radio Bursts (FRBs). For blazars, these include the adoption of: i) a reference filter to allow reconstruction of a well-sampled light curve not affected by colour changes effects; ii) two snapshots/visit with different exposure times to avoid saturation during flaring states; iii) a rolling cadence to get better-sampled light curves at least in some time intervals. We also address the potential importance of Target of Opportunity (ToO) observations of blazar neutrino sources, and the advantages of a Minisurvey with a star trail cadence (see white paper by David Thomas et al.) for both the blazar science and the detection of possible very fast optical counterparts of FRBs.

The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $\gamma$-loud blazar whose optical emission is dominated by the thermal radiation from accretion disc. 4C $+$71.07 has been detected in outburst twice by the AGILE $\gamma$-ray satellite during the period end of October - mid November 2015, when it reached a $\gamma$-ray flux of the order of $F_{\rm E>100\,MeV} = (1.2 \pm 0.3)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ and $F_{\rm E>100\,MeV} = (3.1 \pm 0.6)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, respectively, allowing us to investigate the properties of the jet and of the emission region. We investigated its spectral energy distribution by means of almost simultaneous observations covering the cm, mm, near-infrared, optical, ultra-violet, X-ray and $\gamma$-ray energy bands obtained by the GASP-WEBT Consortium, the Swift and the AGILE and Fermi satellites. The spectral energy distribution of the second $\gamma$-ray flare (the one whose energy coverage is more dense) can be modelled by means of a one-zone leptonic model, yielding a total jet power of about $4\times10^{47}$ erg s$^{-1}$. During the most prominent $\gamma$-ray flaring period our model is consistent with a dissipation region within the broad-line region. Moreover, this class of high-redshift, large-mass black-hole flat-spectrum radio quasars might be good targets for future $\gamma$-ray satellites such as e-ASTROGAM.

PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5{\sigma}) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broad-band emission is compatible with the EC scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR.

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.

The detection of gamma-ray emission from narrow-line Seyfert 1 galaxies (NLSy1) has challenged the idea that large black hole (BH) masses ($\ge$10$^8$ M$_{\odot}$) are needed to launch relativistic jets. We present near-infrared imaging data of the gamma-ray-emitting NLSy1 PKS 1502+036 obtained with the Very Large Telescope. Its surface brightness profile, extending to $\sim$ 20 kpc, is well described by the combination of a nuclear component and a bulge with a Sersic index $n$ = 3.5, which is indicative of an elliptical galaxy. A circumnuclear structure observed near PKS 1502+036 may be the result of galaxy interactions. A BH mass of about $\sim 7 \times 10^{8}$ M$_{\odot}$ has been estimated by the bulge luminosity. The presence of an additional faint disc component cannot be ruled out with the present data, but this would reduce the BH mass estimate by only $\sim$ 30%. These results, together with analogous findings obtained for FBQS J1644+2619, indicate that the relativistic jets in gamma-ray-emitting NLSy1 are likely produced by massive black holes at the center of elliptical galaxies.

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.

Blazar emission is dominated by non-thermal radiation from a relativistic jet pointing toward us, therefore undergoing Doppler beaming. This is responsible for flux enhancement and contraction of the variability time scales, so that most blazars appear as luminous sources characterized by noticeable and fast flux changes at all frequencies. The mechanisms producing their unpredictable variability are debated and include injection, acceleration and cooling of particles, with possible intervention of shock waves or turbulence. Changes in the viewing angle of the emitting knots or jet regions have also been suggested to explain flaring events or specific properties such as intraday variability, quasi-periodicities, or the delay of radio flux variations relative to optical changes. However, such a geometric interpretation has not been universally accepted because alternative explanations based on changes of physical conditions can also work in many cases. Here we report the results of optical-to-radio monitoring of the blazar CTA 102 by the Whole Earth Blazar Telescope Collaboration and show that the observed long-term flux and spectral variability is best explained by an inhomogeneous, curved jet that undergoes orientation changes. We propose that magnetohydrodynamic instabilities or rotation of a twisted jet cause different jet regions to change their orientation and hence their relative Doppler factors. In particular, the recent extreme optical outburst (six magnitudes) occurred when the corresponding jet emitting region acquired a minimum viewing angle.

We report on long-term multi-wavelength monitoring of blazar Mrk 421 by the GASP-WEBT collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the optical, X-ray, and \gamma-ray fluxes is very variable. The \gamma-ray flux variations show fair correlation with the optical ones starting from 2012. We analyse spectropolarimetric data and find wavelength-dependence of the polarisation degree (P), which is compatible with the presence of the host galaxy, and no wavelength-dependence of the electric vector polarisation angle (EVPA). Optical polarimetry shows a lack of simple correlation between P and flux and wide rotations of the EVPA. We build broad-band spectral energy distributions with simultaneous near-infrared and optical data from the GASP-WEBT and UV and X-ray data from the Swift satellite. They show strong variability in both flux and X-ray spectral shape and suggest a shift of the synchrotron peak up to a factor $\sim 50$ in frequency. The interpretation of the flux and spectral variability is compatible with jet models including at least two emitting regions that can change their orientation with respect to the line of sight.

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.

The discovery of $\gamma$-ray emission from radio-loud narrow-line Seyfert 1 (NLSy1) galaxies has questioned the need for large black hole masses (> 10$^8$ M$_{\odot}$) to launch relativistic jets. We present near-infrared data of the $\gamma$-ray-emitting NLSy1 FBQS J1644+2619 that were collected using the camera CIRCE (Canarias InfraRed Camera Experiment) at the 10.4-m Gran Telescopio Canarias to investigate the structural properties of its host galaxy and to infer the black hole mass. The 2D surface brightness profile is modelled by the combination of a nuclear and a bulge component with a S\'ersic profile with index $n$ = 3.7, indicative of an elliptical galaxy. The structural parameters of the host are consistent with the correlations of effective radius and surface brightness against absolute magnitude measured for elliptical galaxies. From the bulge luminosity, we estimated a black hole mass of (2.1$\pm$0.2) $\times$10$^8$ M$_{\odot}$, consistent with the values characterizing radio-loud active galactic nuclei.

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.

We examine the 2008-2016 $\gamma$-ray and optical light curves of three bright BL Lac objects, 0716+714, MRK 421, BL Lac, which exhibit large structured variability. We searched for periodicities by using a fully Bayesian approach. For two out of three sources investigated no significant periodic variability was found. In the case of BL Lac we detected a periodicity of ~ 680 days. Although the signal related to this is modest, the coincidence of the periods in both gamma and optical bands is indicative of a physical relevance. Considering previous literature results, possibly related $\gamma$-ray and optical periodicities of about one year time scale are proposed in 4 bright $\gamma$-ray blazars out of the 10 examined in detail. Comparing with results from periodicity search of optical archives of quasars, the presence of quasi-periodicities in blazars might be more frequent by a large factor. This suggests the intriguing possibility that the basic conditions for their observability are related to the relativistic jet in the observer direction, but the overall picture remains uncertain.

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.

In 2015 July 29 - September 1 the satellite XMM-Newton pointed at the BL Lac object PG 1553+133 six times, collecting data for 218 hours. During one of these epochs, simultaneous observations by the Swift satellite were requested to compare the results of the X-ray and optical-UV instruments. Optical, near-infrared and radio monitoring was carried out by the Whole Earth Blazar Telescope (WEBT) collaboration for the whole observing season. We here present the results of the analysis of all these data, together with an investigation of the source photometric and polarimetric behaviour over the last three years. The 2015 EPIC spectra show slight curvature and the corresponding light curves display fast X-ray variability with a time scale of the order of 1 hour. In contrast to previous results, during the brightest X-ray states detected in 2015 the simple log-parabolic model that best-fits the XMM-Newton data also reproduces reasonably well the whole synchrotron bump, suggesting a peak in the near-UV band. We found evidence of a wide rotation of the polarization angle in 2014, when the polarization degree was variable, but the flux remained almost constant. This is difficult to interpret with deterministic jet emission models, while it can be easily reproduced by assuming some turbulence of the magnetic field.

Blazars are the most luminous and variable AGNs, and thus excellent probes of accretion and emission processes close to the central engine. We focus on PKS 1510-089 ($z=0.36$), one of the brightest gamma-ray sources in the Fermi LAT catalog, to study its complex multi-wavelength variability. PKS 1510-089 was observed twice in hard X-rays with the IBIS instrument onboard INTEGRAL during the flares of Jan 2009 and Jan 2010, and simultaneously with Swift and NOT, in addition to the constant Fermi monitoring. The optical polarization was measured in several bands on 18 Jan 2010 at the NOT. Using our and archival data we constructed historical light curves at gamma-to-radio wavelengths covering nearly 20 years and applied variability tests. We assembled SEDs in 2009 and 2010 and compared them with those at two previous epochs and with a model based on synchrotron and inverse Compton (IC) radiation. The SED modeling suggests that the physical quantities that undergo the largest variations are the total power injected into the emitting region and the random Lorentz factor of the electron distribution cooling break, that are higher in the higher gamma-ray states. This suggests a correlation of the injected power with enhanced activity of the acceleration mechanism. The cooling likely takes place at a much smaller distance ($\sim$0.03 pc) than the BLR radius. The emission at a few hundred GeV can be reproduced with IC scattering of highly relativistic electrons off FIR photons at $\sim$0.2 pc, presumably in a dusty torus. DCF analysis between the long-term optical and gamma-ray light curves yields a good correlation with no measurable delay. Our time analysis of the RXTE PCA and Fermi LAT light curves reveals no obvious (quasi-)periodicities, up to the maximum time scale (a few years) probed by the light curves, which are severely affected by red noise.

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $\gamma$-ray band with a statistical significance of 5.9 $\sigma$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $\gamma$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

The discovery by the Large Area Telescope on board Fermi of variable gamma-ray emission from radio-loud narrow-line Seyfert 1 (NLSy1) galaxies revealed the presence of a possible third class of Active Galactic Nuclei (AGN) with relativistic jets in addition to blazars and radio galaxies. Considering that NLSy1 are usually hosted in spiral galaxies, this finding poses intriguing questions about the nature of these objects and the formation of relativistic jets. We report on a systematic investigation of the gamma-ray properties of a sample of radio-loud NLSy1, including the detection of new objects, using 7 years of Fermi-LAT data with the new Pass 8 event-level analysis. In addition we discuss the radio-to-very-high-energy properties of the gamma-ray emitting NLSy1, their host galaxy, and black hole mass in the context of the blazar scenario and the unification of relativistic jets at different scales.

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.

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.

In a previous paper, we proposed a new method to select low-power BL Lacs (LPBLs) based on mid-infrared emission and flux contrast through the Ca II spectral break; that study led to the selection of a complete sample formed by 34 LPBLs with 0.05<z<=0.15 and radio luminosities spanning the range log(L_r) = 39.2-41.5 [erg/s]. We now assemble the broadband spectral energy distributions (SEDs) of these sources to investigate their nature and compare them with brighter BL Lacs. We find that the ratios between the X-ray and radio luminosities range from ~20 to ~30000 and that the synchrotron peak frequencies span a wide energy interval, from log(nu_peak)~13.5 to ~20 [Hz]. This indicates a broad variety of SED shapes and a mixture of BL Lac flavors. Indeed, although the majority of our LPBLs are high-energy peaked BL Lacs (HBLs), we find that a quarter of them are low-energy peaked BL Lacs (LBLs), despite the fact that the sample is biased against the selection of LBLs. The analysis of the median LPBL SED confirms disagreement with the blazar sequence at low radio luminosities. Furthermore, if we limit the sample to the LBLs subsample, we find that their median SED shape is essentially indistinguishable from that of the most luminous BL Lacs. We conclude that the observed radio power is not the main driving parameter of the multiwavelength properties of BL Lacs.

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.

A multifrequency campaign on the BL Lac object PG 1553+113 was organized by the Whole Earth Blazar Telescope (WEBT) in 2013 April-August, involving 19 optical, two near-IR, and three radio telescopes. The aim was to study the source behaviour at low energies during and around the high-energy observations by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes in April-July. We also analyse the UV and X-ray data acquired by the Swift and XMM-Newton satellites in the same period. The WEBT and satellite observations allow us to detail the synchrotron emission bump in the source spectral energy distribution (SED). In the optical we found a general bluer-when-brighter trend. The X-ray spectrum remained stable during 2013, but a comparison with previous observations suggests that it becomes harder when the X-ray flux increases. The long XMM-Newton exposure reveals a curved X-ray spectrum. In the SED, the XMM-Newton data show a hard near-UV spectrum, while Swift data display a softer shape that is confirmed by previous HST-COS and IUE observations. Polynomial fits to the optical-X-ray SED show that the synchrotron peak likely lies in the 4-30 eV energy range, with a general shift towards higher frequencies for increasing X-ray brightness. However, the UV and X-ray spectra do not connect smoothly. Possible interpretations include: i) orientation effects, ii) additional absorption, iii) multiple emission components, and iv) a peculiar energy distribution of relativistic electrons. We discuss the first possibility in terms of an inhomogeneous helical jet model.

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.

Among active galactic nuclei, BL Lac objects show extreme properties that have been interpreted as the effect of relativistic beaming on the emission from a plasma jet oriented close to the line of sight. The Doppler amplification of the jet emission makes them ideal targets for studying jet physics. In particular, low-power BL Lacs (LPBL) are very interesting because they probe the jet formation and emission processes at the lowest levels of accretion. However, they are difficult to identify since their emission is swamped by the radiation from the host galaxy in most observing bands. In this paper we propose a new LPBL selection method based on the mid-infrared emission, in addition to the traditional optical indices. We considered the radio-selected sample of Best & Heckman (2012, MNRAS, 421, 1569) and cross-matched it with the WISE all-sky survey. In a new diagnostic plane including the W2-W3 color and the Dn(4000) index, LPBL are located in a region scarcely populated by other sources. By filtering objects with small emission line equivalent width, we isolated 36 LPBL candidates up to redshift 0.15. Their radio luminosity at 1.4 GHz spans the range log L_r = 39.2-41.5 [erg/s]. Considering the completeness of our sample, we analyzed the BL Lac luminosity function (RLF), finding a dramatic paucity of LPBL with respect to the extrapolation of the RLF toward low power. This requires a break in the RLF located at log L_r~40.6 [erg/s]. The consequent peak in the BL Lacs number density is possibly the manifestation of a minimum power required to launch a relativistic jet.

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 an analysis of the multiwavelength behaviour of the blazar OJ 248 at z = 0.939 in the period 2006-2013. We use low-energy data (optical, near-infrared, and radio) obtained by 21 observatories participating in the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from the Swift (optical-UV and X-rays) and Fermi (gamma-rays) satellites, to study flux and spectral variability and correlations among emissions in different bands. We take into account the effect of absorption by the Damped Lyman Alpha intervening system at z = 0.525. Two major outbursts were observed in 2006-2007 and in 2012-2013 at optical and near-IR wavelengths, while in the high-frequency radio light curves prominent radio outbursts are visible peaking at the end of 2010 and beginning of 2013, revealing a complex radio-optical correlation. Cross-correlation analysis suggests a delay of the optical variations after the gamma-ray ones of about a month, which is a peculiar behaviour in blazars. We also analyse optical polarimetric and spectroscopic data. The average polarization percentage P is less than 3 per cent, but it reaches about 19 per cent during the early stage of the 2012-2013 outburst. A vague correlation of P with brightness is observed. There is no preferred electric vector polarisation angle and during the outburst the linear polarization vector shows wide rotations in both directions, suggesting a complex behaviour or structure of the jet and possible turbulence. The analysis of 140 optical spectra acquired at the Steward Observatory reveals a strong Mg II broad emission line with an essentially stable flux of 6.2 e-15 erg cm-2 s-1 and a full width at half-maximum of 2053 km s-1.

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

Blazars are astrophysical sources whose emission is dominated by non-thermal processes, typically interpreted as synchrotron and inverse Compton emission. Although the general picture is rather robust and consistent with observations, many aspects are still unexplored. Polarimetric monitoring can offer a wealth of information about the physical processes in blazars. Models with largely different physical ingredients can often provide almost indistinguishable predictions for the total flux, but usually are characterized by markedly different polarization properties. We explore, with a pilot study, the possibility to derive structural information about the emitting regions of blazars by means of a joint analysis of rapid variability of the total and polarized flux at optical wavelengths. Short timescale (from tens of seconds to a couple of minutes) optical linear polarimetry and photometry for two blazars, BL Lacertae and PKS 1424+240, was carried out with the PAOLO polarimeter at the 3.6m Telescopio Nazionale Galileo. Several hours of almost continuous observations were obtained for both sources. Our intense monitoring allowed us to draw strongly different scenarios for BL Lacertae and PKS 1424+240, with the former characterized by intense variability on time-scales from hours to a few minutes and the latter practically constant in total flux. Essentially the same behavior is observed for the polarized flux and the position angle. The variability time-scales turned out to be as short as a few minutes, although involving only a few percent variation of the flux. The polarization variability time-scale is generally consistent with the total flux variability. Total and polarized flux appear to be essentially uncorrelated. However, even during our relatively short monitoring, different regimes can be singled out. (abridged)

We constrain the behavior of the radio luminosity function (RLF) of two classes of active galactic nuclei (AGN) namely AGN of low radio power (LRP) and BL Lac objects. The extrapolation of the observed steep RLFs to low power predicts a space density of such objects that exceeds that of the sources that can harbor them and this requires a break to a shallower slope. For LRP AGN we obtain P_br,LRP > 10^20.5 W/Hz at 1.4 GHz to limit their density to be smaller than that of elliptical galaxies with black hole masses M_BH > 10^7.5 solar masses. By combining this value with the limit derived by the observations the break must occur at P_br,LRP~10^20.5-10^21.5 W/Hz. For BL Lacs we find P_br,BLLAC > 10^23.3 W/Hz otherwise they would outnumber the density of weak-lined and compact radio sources, while the observations indicate P_br,BLLAC < 10^24.5 W/Hz. In the framework of the AGN unified model a low luminosity break in the RLF of LRP AGN must correspond to a break in the RLF of BL Lacs. The ratio between P_br,LRP and P_br,BLLAC is ~10^3, as expected for a jet Doppler factor of ~10.

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

We report on multifrequency observations performed during 2012 December-2013 August of the first narrow-line Seyfert 1 galaxy detected in gamma rays, PMN J0948+0022 ($z$ = 0.5846). A gamma-ray flare was observed by the Large Area Telescope on board Fermi during 2012 December-2013 January, reaching a daily peak flux in the 0.1-100 GeV energy range of (155 $\pm$ 31) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring period triggered Swift and VERITAS observations in addition to radio and optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the gamma-ray flare, reaching a record flux for this source from optical to gamma rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F (> 0.2 TeV) < 4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$. We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.

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

Multiwavelength variability of blazars offers indirect insight into their powerful engines and on the mechanisms through which energy is propagated from the centre down the jet. The BL Lac object Mkn 421 is a TeV emitter, a bright blazar at all wavelengths, and therefore an excellent target for variability studies. Mkn 421 was observed by INTEGRAL and Fermi-LAT in an active state on 16-21 April 2013. Well sampled optical, soft, and hard X-ray light curves show the presence of two flares. The average flux in the 20-100 keV range is 9.1e-11 erg/s/cm2 (~4.5 mCrab) and the nuclear average apparent magnitude, corrected for Galactic extinction, is V ~12.2. In the time-resolved X-ray spectra (3.5-60 keV), which are described by broken power laws and, marginally better, by log-parabolic laws, we see a hardening that correlates with flux increase, as expected in refreshed energy injections in a population of electrons that later cool via synchrotron radiation. The hardness ratios between the JEM-X fluxes in two different bands and between the JEM-X and IBIS/ISGRI fluxes confirm this trend. During the observation, the variability level increases monotonically from the optical to the hard X-rays, while the large LAT errors do not allow a significant assessment of the MeV-GeV variability. The cross-correlation analysis during the onset of the most prominent flare suggests a monotonically increasing delay of the lower frequency emission with respect to that at higher frequency, with a maximum time-lag of about 70 minutes, that is however not well constrained. The spectral energy distributions from the optical to the TeV domain are satisfactorily described by homogeneous models of blazar emission based on synchrotron radiation and synchrotron self-Compton scattering, except in the state corresponding to the LAT softest spectrum and highest flux.

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.

The infrared properties of blazars can be studied from the statistical point of view with the help of sky surveys, like that provided by the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All Sky Survey (2MASS). However, these sources are known for their strong and unpredictable variability, which can be monitored for a handful of objects only. In this paper we consider the 28 blazars (14 BL Lac objects and 14 flat-spectrum radio quasars, FSRQs) that are regularly monitored by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) since 2007. They show a variety of infrared colours, redshifts, and infrared-optical spectral energy distributions (SEDs), and thus represent an interesting mini-sample of bright blazars that can be investigated in more detail. We present near-IR light curves and colours obtained by the GASP from 2007 to 2013, and discuss the infrared-optical SEDs. These are analysed with the aim of understanding the interplay among different emission components. BL Lac SEDs are accounted for by synchrotron emission plus an important contribution from the host galaxy in the closest objects, and dust signatures in 3C 66A and Mkn 421. FSRQ SEDs require synchrotron emission with the addition of a quasar-like contribution, which includes radiation from a generally bright accretion disc, broad line region, and a relatively weak dust torus.

We report on radio-to-gamma-ray observations during 2011 May-September of PMN J0948+0022, the first narrow-line Seyfert 1 (NLSy1) galaxy detected in gamma-rays by Fermi-LAT. Strong variability was observed in gamma-rays, with two flaring periods peaking on 2011 June 20 and July 28. The variability observed in optical and near-infrared seems to have no counterpart in gamma-rays. This different behaviour could be related to a bending and inhomogeneous jet or a turbulent extreme multi-cell scenario. The radio spectra showed a variability pattern typical of relativistic jets. The XMM spectrum shows that the emission from the jet dominates above 2 keV, while a soft X-ray excess is evident in the low-energy part of the X-ray spectrum. Models where the soft emission is partly produced by blurred reflection or Comptonisation of the thermal disc emission provide good fits to the data. The X-ray spectral slope is similar to that found in radio-quiet NLSy1, suggesting that a standard accretion disc is present, as expected from the high accretion rate. Except for the soft X-ray excess, unusual in jet-dominated AGNs, PMN J0948+0022 shows all characteristics of the blazar class.

Since the launch of the Fermi satellite, BL Lacertae has been moderately active at gamma-rays and optical frequencies until May 2011, when the source started a series of strong flares. The exceptional optical sampling achieved by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) in collaboration with the Steward Observatory allows us to perform a detailed comparison with the daily gamma-ray observations by Fermi. Discrete correlation analysis between the optical and gamma-ray emission reveals correlation with a time lag of 0 +- 1 d, which suggests cospatiality of the corresponding jet emitting regions. A better definition of the time lag is hindered by the daily gaps in the sampling of the extremely fast flux variations. In general, optical flares present more structure and develop on longer time scales than corresponding gamma-ray flares. Observations at X-rays and at millimetre wavelengths reveal a common trend, which suggests that the region producing the mm and X-ray radiation is located downstream from the optical and gamma-ray-emitting zone in the jet. The mean optical degree of polarisation slightly decreases over the considered period and in general it is higher when the flux is lower. The optical electric vector polarisation angle (EVPA) shows a preferred orientation of about 15 deg, nearly aligned with the radio core EVPA and mean jet direction. Oscillations around it increase during the 2011-2012 outburst. We investigate the effects of a geometrical interpretation of the long-term flux variability on the polarisation. A helical magnetic field model predicts an evolution of the mean polarisation that is in reasonable agreement with the observations. These can be fully explained by introducing slight variations in the compression factor in a transverse shock waves model.

The narrow-line Seyfert 1 galaxy SBS 0846+513 was first detected by the Large Area Telescope (LAT) on-board Fermi in 2011 June-July when it underwent a period of flaring activity. Since then, as Fermi continues to accumulate data on this source, its flux has been monitored on a daily basis. Two further gamma-ray flaring episodes from SBS 0846+513 were observed in 2012 May and August, reaching a daily peak flux integrated above 100 MeV of (50+/-12)x10^-8 ph/cm^2/s, and (73+/-14)x10^-8 ph/cm^2/s on May 24 and August 7, respectively. Three outbursts were detected at 15 GHz by the Owens Valley Radio Observatory 40-m telescope in 2012 May, 2012 October, and 2013 January, suggesting a complex connection with the gamma-ray activity. The most likely scenario suggests that the 2012 May gamma-ray flare may not be directly related to the radio activity observed over the same period, while the two gamma-ray flaring episodes may be related to the radio activity observed at 15 GHz in 2012 October and 2013 January. The gamma-ray flare in 2012 May triggered Swift observations that confirmed that SBS 0846+513 was also exhibiting high activity in the optical, UV and X-ray bands, thus providing a firm identification between the gamma-ray source and the lower-energy counterpart. We compared the spectral energy distribution (SED) of the flaring state in 2012 May with that of a quiescent state. The two SEDs, modelled as an external Compton component of seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. No significant evidence of thermal emission from the accretion disc has been observed. Interestingly, in the 5 GHz radio luminosity vs. synchrotron peak frequency plot SBS 0846+513 seems to lie in the flat spectrum radio quasar part of the so-called `blazar sequence'.

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 report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival radio-to-X-ray data. The source was not active at gamma-ray energies during the first two years of Fermi operation. A significant increase in activity was observed during 2010 October-2011 August. In particular a strong gamma-ray flare was observed in 2011 June reaching an isotropic gamma-ray luminosity (0.1-300 GeV) of 1.0x10^48 erg/s, comparable to that of the brightest flat spectrum radio quasars, and showing spectral evolution in gamma rays. An apparent superluminal velocity of (8.2+/-1.5)c in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indications of the presence of a relativistic jet as powerful as those of blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

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 investigate the optical spectral properties of the blazar PKS 1222+216 during a period of 3 years. While the continuum is highly variable the broad line emission is practically constant. This supports a scenario in which the broad line region is not affected by jet continuum variations. We thus infer the thermal component of the continuum from the line luminosity and we show that it is comparable with the continuum level observed during the phases of minimum optical activity. The mass of the black hole is estimated through the virial method from the FWHM of MgII, Hbeta, and Halpha broad lines and from the thermal continuum luminosity. This yields a consistent black hole mass value of 6x10^8 solar masses.

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 multiwavelength observations of 3C 454.3 from April 2008 to March 2010. The radio to optical data are mostly from the GASP-WEBT, UV and X-ray data from Swift, and gamma-ray data from the AGILE and Fermi satellites. We improved the calibration of optical-UV data from the UVOT and OM instruments and estimated the Lyalpha flux to disentangle the contributions from different components in this spectral region. The observations reveal prominent variability above 8 GHz. In the optical-UV band, the variability amplitude decreases with increasing frequency due to a steadier radiation from both a broad line region and an accretion disc. The optical flux reaches nearly the same levels in the 2008-2009 and 2009-2010 observing seasons; the mm one shows similar behaviour, whereas the gamma and X-ray flux levels rise in the second period. Two prominent gamma-ray flares in mid 2008 and late 2009 show a double-peaked structure, with a variable gamma/optical flux ratio. The X-ray flux variations seem to follow the gamma-ray and optical ones by about 0.5 and 1 d, respectively. We interpret the multifrequency behaviour in terms of an inhomogeneous curved jet, where synchrotron radiation of increasing wavelength is produced in progressively outer and wider jet regions, which can change their orientation in time. In particular, we assume that the long-term variability is due to this geometrical effect. By combining the optical and mm light curves to fit the gamma and X-ray ones, we find that the gamma (X-ray) emission may be explained by inverse-Comptonisation of synchrotron optical (IR) photons by their parent relativistic electrons (SSC process). A slight, variable misalignment between the synchrotron and Comptonisation zones would explain the increased gamma and X-ray flux levels in 2009-2010, as well as the change in the gamma/optical flux ratio during the outbursts peaks.

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

We report on 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.

During the first 3 years of operation the Gamma-Ray Imaging Detector onboard the AGILE satellite detected several blazars in a high gamma-ray activity: 3C 279, 3C 454.3, PKS 1510-089, S5 0716+714, 3C 273, W Comae, Mrk 421, PKS 0537-441 and 4C +21.35. Thanks to the rapid dissemination of our alerts, we were able to obtain multiwavelength data from other observatories such as Spitzer, Swift, RXTE, Suzaku, INTEGRAL, MAGIC, VERITAS, and ARGO as well as radio-to-optical coverage by means of the GASP Project of the WEBT and the REM Telescope. This large multifrequency coverage gave us the opportunity to study the variability correlations between the emission at different frequencies and to obtain simultaneous spectral energy distributions of these sources from radio to gamma-ray energy bands, investigating the different mechanisms responsible for their emission and uncovering in some cases a more complex behaviour with respect to the standard models. We present a review of the most interesting AGILE results on these gamma-ray blazars and their multifrequency data.

Extragalactic nuclear activity is best explored with observations at high energies, where the most extreme flux and spectral variations are expected to occur, witnessing changes in the accretion flow or in the kinematics of the plasma. In active galactic nuclei of blazar type, these variations are the most dramatic. By following blazar outbursts from their onset and by correlating the observed variations at many different wavelengths we can reconstruct the behavior of the plasma and map out the development of the flare within the jet. The advent of the Fermi satellite has allowed the start of a systematic and intensive monitoring program of blazars. Blazar outbursts are very effectively detected by the LAT instrument in the MeV-GeV domain, and these can be promptly followed up with other facilities. Based on a Fermi LAT detection of a high MeV-GeV state, we have observed the blazar PKS 1502+106 with the INTEGRAL satellite between 9 and 11 August 2008. Simultaneous Swift observations have been also accomplished, as well as optical follow-up at the Nordic Optical Telescope. The IBIS instrument onboard INTEGRAL detected a source at a position inconsistent with the optical coordinates of PKS 1502+106, but consistent with those of the Seyfert 1 galaxy Mkn 841, located at 6.8 arcmin south-west from the blazar, which is therefore responsible for all the hard X-ray flux detected by IBIS. At the location of the blazar, IBIS sets an upper limit of ~10^{-11} erg/s/cm2 on the 15-60 keV flux, that turns out to be consistent with a model of inverse Compton scattering accounting for the soft X-ray and gamma-ray spectra measured by Swift XRT and Fermi LAT, respectively. The gamma-ray spectrum during the outburst indicates substantial variability of the characteristic energy of the inverse Compton component in this blazar.(abridged).

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

The blazar PG 1553+113 is a well known TeV gamma-ray emitter. In this paper, we determine its spectral energy distribution using simultaneous multi-frequency data in order to study its emission processes. An extensive campaign was carried out between March and April 2008, where optical, X-ray, high-energy (HE) gamma-ray, and very-high-energy (VHE) gamma-ray data were obtained with the KVA, Abastumani, REM, RossiXTE/ASM, AGILE and MAGIC telescopes, respectively. This is the first simultaneous broad-band (i.e., HE+VHE) gamma-ray observation, though AGILE did not detect the source. We combine data to derive source's spectral energy distribution and interpret its double peaked shape within the framework of a synchrotron self compton model

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

We have analysed optical spectra of BL Lacertae, the prototype of its class, to verify the presence and possible flux variations of its broad Ha line. We used the spectroscopic information also to investigate the question of its parent population. Four spectra were acquired at the TNG in 2007-2008, when the source was in a relatively faint state. In three cases we were able to measure the broad Ha and several narrow emission lines. A comparison with previous results suggests that the broad Ha has increased by ~50% in ten years, a change not unusual for Broad Lined AGN. We estimated a black hole mass for BL Lac of 4-6 10^8 solar masses from its relation with the bulge luminosity. The virial mass estimated from the spectroscopic data is instead 20-30 times smaller. We suggest that this discrepancy is due to the BL Lac BLR being underluminous. Finally, we addressed the problem of the BL Lac parent population, comparing its isotropic quantities with those of other AGN classes. From the point of view of the narrow emission line spectrum, the source locates close to low-excitation radio galaxies. When also considering its diffuse radio power, an association with FRI radio galaxies is severely questioned due to the lower radio luminosity (at give line luminosity) of BL Lac. The narrow line and radio luminosities of BL Lac instead match those of a sample of miniature radio galaxies, which however do not show a BLR. Yet, if existing, "misaligned BL Lac" objects should have entered that sample. We also rule out the possibility that they have been excluded because of a QSO optical appearance. The observational constraints suggest that BL Lac is caught in a short term transient stage, which does not leave a detectable evolutionary "trace" in the AGN population. We speculate on a scenario that can account for the observed properties. [ABRIDGED]

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.

Since 1997, BL Lacertae has undergone a phase of high optical activity, with the occurrence of several prominent outbursts. Starting from 1999, the Whole Earth Blazar Telescope (WEBT) consortium has organized various multifrequency campaigns on this blazar, collecting tens of thousands of data points. One of the main issues in the analysis of this huge dataset has been the study of colour variability. The massive amount of optical and near-infrared data collected during the campaigns enables us to perform a deep analysis of multiband data, with the aim of understanding the flux variability mechanisms. We use a new approach for the analysis of these data, focusing on the source spectral evolution. We show that the overall behaviour of the BL Lacertae light and colour curves can be explained in terms of changing viewing angle of a moving, discrete emitting region, which causes variable Doppler boosting of the corresponding radiation. A fractal helical structure is suggested to be at the origin of the different time scales of variability.

We report on the second AGILE multiwavelength campaign of the blazar 3C 454.3 during the first half of December 2007. This campaign involved AGILE, Spitzer, Swift,Suzaku,the WEBT consortium,the REM and MITSuME telescopes,offering a broad band coverage that allowed for a simultaneous sampling of the synchrotron and inverse Compton (IC) emissions.The 2-week AGILE monitoring was accompanied by radio to optical monitoring by WEBT and REM and by sparse observations in mid-Infrared and soft/hard X-ray energy bands performed by means of Target of Opportunity observations by Spitzer, Swift and Suzaku, respectively.The source was detected with an average flux of~250x10^{-8}ph cm^-2s^-1 above 100 MeV,typical of its flaring states.The simultaneous optical and gamma-ray monitoring allowed us to study the time-lag associated with the variability in the two energy bands, resulting in a possible ~1-day delay of the gamma-ray emission with respect to the optical one. From the simultaneous optical and gamma-ray fast flare detected on December 12, we can constrain the delay between the gamma-ray and optical emissions within 12 hours. Moreover, we obtain three Spectral Energy Distributions (SEDs) with simultaneous data for 2007 December 5, 13, 15, characterized by the widest multifrequency coverage. We found that a model with an external Compton on seed photons by a standard disk and reprocessed by the Broad Line Regions does not describe in a satisfactory way the SEDs of 2007 December 5, 13 and 15. An additional contribution, possibly from the hot corona with T=10^6 K surrounding the jet, is required to account simultaneously for the softness of the synchrotron and the hardness of the inverse Compton emissions during those epochs.

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

PKS 1510-089 is a powerful Flat Spectrum Radio Quasar at z=0.361 with radiative output dominated by the gamma-ray emission. In the last two years PKS 1510-089 showed high variability over all the electromagnetic spectrum and in particular very high gamma-ray activity was detected by the Gamma-Ray Imaging Detector on board the AGILE satellite with flaring episodes in August 2007 and March 2008. An extraordinary activity was detected in March 2009 with several flaring episodes and a flux that reached 500 x 10^{-8} photons cm^{-2} s^{-1}. Multiwavelength observations of PKS 1510-089 seem to indicate the presence of Seyfert-like features such as little blue bump and big blue bump. Moreover, X-ray observations suggested the presence of a soft X-ray excess that could be a feature of the bulk Comptonization mechanism. We present the results of the analysis of the multiwavelength data collected by GASP-WEBT, REM, Swift and AGILE during these gamma-ray flares and the theoretical implications for the emission mechanisms.

We performed an optical spectroscopic monitoring of the blazar 3C 454.3 from September 2003 to July 2008. Sixteen optical spectra were obtained during different runs, which constitute the first spectroscopic monitoring done in the rest-frame UV region (z=0.859). An overall flux variation of the MgII (2800 A) by a factor ~ 3 was observed, while the corresponding UV continuum (F_cont at 3000 A) changed by a factor ~ 14. The MgII emission lines respond proportionally to the continuum variations when the source is in a low-activity state. In contrast, near the optical outbursts detected in 2005 and 2007, the MgII emission lines showed little response to the continuum flux variations. During the monitored period the UV FeII flux changed by a factor ~ 6 and correlated with F_cont (r = 0.92). A negative correlation between EW(Mg II) and F_cont was found, i.e. the so-called "Intrinsic Baldwin Effect".

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.

In 2007-2008 we carried out a new multiwavelength campaign of the Whole Earth Blazar Telescope (WEBT) on BL Lacertae, involving three pointings by the XMM-Newton satellite, to study its emission properties. The source was monitored in the optical-to-radio bands by 37 telescopes. The brightness level was relatively low. Some episodes of very fast variability were detected in the optical bands. The X-ray spectra are well fitted by a power law with photon index of about 2 and photoelectric absorption exceeding the Galactic value. However, when taking into account the presence of a molecular cloud on the line of sight, the data are best fitted by a double power law, implying a concave X-ray spectrum. The spectral energy distributions (SEDs) built with simultaneous radio-to-X-ray data at the epochs of the XMM-Newton observations suggest that the peak of the synchrotron emission lies in the near-IR band, and show a prominent UV excess, besides a slight soft-X-ray excess. A comparison with the SEDs corresponding to previous observations with X-ray satellites shows that the X-ray spectrum is extremely variable. We ascribe the UV excess to thermal emission from the accretion disc, and the other broad-band spectral features to the presence of two synchrotron components, with their related SSC emission. We fit the thermal emission with a black body law and the non-thermal components by means of a helical jet model. The fit indicates a disc temperature greater than 20000 K and a luminosity greater than 6 x 10^44 erg/s.

Since 2001, the radio quasar 3C 454.3 has undergone a period of high optical activity, culminating in the brightest optical state ever observed, during the 2004-2005 outburst. The Whole Earth Blazar Telescope (WEBT) consortium has carried out several multifrequency campaigns to follow the source behaviour. The GLAST-AGILE Support Program (GASP) was born from the WEBT to provide long-term continuous optical-to-radio monitoring of a sample of gamma-loud blazars, during the operation of the AGILE and GLAST (now known as Fermi GST) gamma-ray satellites. The main aim is to shed light on the mechanisms producing the high-energy radiation, through correlation analysis with the low-energy emission. Thus, since 2008 the monitoring task on 3C 454.3 passed from the WEBT to the GASP, while both AGILE and Fermi detected strong gamma-ray emission from the source. We present the main results obtained by the GASP at optical, mm, and radio frequencies in the 2008-2009 season, and compare them with the WEBT results from previous years. An optical outburst was observed to peak in mid July 2008, when Fermi detected the brightest gamma-ray levels. A contemporaneous mm outburst maintained its brightness for a longer time, until the cm emission also reached the maximum levels. The behaviour compared in the three bands suggests that the variable relative brightness of the different-frequency outbursts may be due to the changing orientation of a curved inhomogeneous jet. The optical light curve is very well sampled during the entire season, which is also well covered by the various AGILE and Fermi observing periods. The relevant cross-correlation studies will be very important in constraining high-energy emission models.

We present optical to very-high energy (VHE) gamma-ray observations of Mrk 421 between 2008 May 24 and June 23. A high-energy (HE) gamma-ray signal was detected by AGILE-GRID during June 9-15, brighter than the average flux observed by EGRET in Mrk 421 by a factor of approx. 1.5. In 20-60 keV X-rays, a large-amplitude 5-day flare (June 9-15) was resolved with a maximum flux of approx. 55 mCrab. SuperAGILE, RXTE/ASM and Swift/BAT data show a clearly correlated flaring structure between soft and hard X-rays, with a high flux/amplitude variability in hard X-rays. Hints of the same flaring behavior is also detected in the simultaneously recorded GASP-WEBT optical data. A target of opportunity observation by Swift near the flare maximum on June 12-13 revealed the highest 2-10 keV flux ever observed (>100 mCrab) and a peak synchrotron energy of approx. 3 keV, a large shift from typical values of 0.5-1 keV. Observations at VHE (E>200 GeV) gamma-rays during June 6-8 show the source flux peaking in a bright state, well correlated with the simultaneous peak in the X-rays. The gamma-ray flare can be interpreted within the framework of the Synchrotron Self Compton model in terms of a rapid acceleration of leptons in the jet.

Since 1997, BL Lacertae has undergone a phase of high optical activity, with the occurrence of several prominent outbursts. Starting from 1999, the Whole Earth Blazar Telescope (WEBT) consortium has organized various multifrequency campaigns on this blazar, collecting tens of thousands of data points. One of the main issues in the study of this huge dataset has been the search for correlations between the optical and radio flux variations, and for possible periodicities in the light curves. The analysis of the data assembled during the first four campaigns (comprising also archival data to cover the period 1968-2003) revealed a fair optical-radio correlation in 1994-2003, with a delay of the hard radio events of ~100 days. Moreover, various statistical methods suggested the existence of a radio periodicity of ~8 years. In 2004 the WEBT started a new campaign to extend the dataset to the most recent observing seasons, in order to possibly confirm and better understand the previous results. In this campaign we have collected and assembled about 11000 new optical observations from twenty telescopes, plus near-IR and radio data at various frequencies. Here, we perform a correlation analysis on the long-term R-band and radio light curves. In general, we confirm the ~100-day delay of the hard radio events with respect to the optical ones, even if longer (~200-300 days) time lags are also found in particular periods. The radio quasi-periodicity is confirmed too, but the "period" seems to progressively lengthen from 7.4 to 9.3 years in the last three cycles. The optical and radio behaviour in the last forty years suggests a scenario where geometric effects play a major role. In particular, the alternation of enhanced and suppressed optical activity (accompanied by hard and soft radio events, respectively) can

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 and analyse the WEBT multifrequency observations of 3C 454.3 in the 2007-2008 observing season, including XMM-Newton observations and near-IR spectroscopic monitoring, and compare the recent emission behaviour with the past one. In the optical band we observed a multi-peak outburst in July-August 2007, and other faster events in November 2007 - February 2008. During these outburst phases, several episodes of intranight variability were detected. A mm outburst was observed starting from mid 2007, whose rising phase was contemporaneous to the optical brightening. A slower flux increase also affected the higher radio frequencies, the flux enhancement disappearing below 8 GHz. The analysis of the optical-radio correlation and time delays, as well as the behaviour of the mm light curve, confirm our previous predictions, suggesting that changes in the jet orientation likely occurred in the last few years. The historical multiwavelength behaviour indicates that a significant variation in the viewing angle may have happened around year 2000. Colour analysis reveals a complex spectral behaviour, which is due to the interplay of different emission components. All the near-IR spectra show a prominent Halpha emission line, whose flux appears nearly constant. The analysis of the XMM-Newton data indicates a correlation between the UV excess and the soft-X-ray excess, which may represent the head and the tail of the big blue bump, respectively. The X-ray flux correlates with the optical flux, suggesting that in the inverse-Compton process either the seed photons are synchrotron photons at IR-optical frequencies or the relativistic electrons are those that produce the optical synchrotron emission. The X-ray radiation would thus be produced in the jet region from where the IR-optical emission comes.

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

[Abridged] We report on a multiwavelength observation of the blazar 3C 454.3 (which we dubbed "crazy diamond") carried out on November 2007 by means of the astrophysical satellites AGILE, INTEGRAL, Swift, the WEBT Consortium, and the optical-NIR telescope REM. 3C 454.3 is detected at a $\sim 19-\sigma$ level during the 3-week observing period, with an average flux above 100 MeV of $F_{\rm E>100MeV} = (170 \pm 13) \times 10^{-8}$ \phcmsec. The gamma-ray spectrum can be fit with a single power-law with photon index $\Gamma_{\rm GRID} = 1.73 \pm 0.16$ between 100 MeV and 1 GeV. We detect significant day-by-day variability of the gamma-ray emission during our observations, and we can exclude that the fluxes are constant at the 99.6% ($\sim 2.9 \sigma$) level. The source was detected typically around 40 degrees off-axis, and it was substantially off--axis in the field of view of the AGILE hard X-ray imager. However, a 5-day long ToO observation by INTEGRAL detected 3C 454.3 at an average flux of about $F_{\rm 20-200 keV} = 1.49 \times 10^{-3}$ \phcmsec with an average photon index of $\Gamma_{\rm IBIS} = 1.75 \pm 0.24$ between 20--200 keV. Swift also detected 3C 454.3 with a flux in the 0.3--10 keV energy band in the range $(1.23-1.40) \times 10^{-2}$ \phcmsec{} and a photon index in the range $\Gamma_{\rm XRT} = 1.56-1.73$. In the optical band, both WEBT and REM show an extremely variable behavior in the $R$ band. A correlation analysis based on the entire data set is consistent with no time-lags between the gamma-ray and the optical flux variations. Our simultaneous multifrequency observations strongly indicate that the dominant emission mechanism between 30 MeV and 30 GeV is dominated by inverse Compton scattering of relativistic electrons in the jet on the external photons from the broad line region.

We report the gamma-ray activity from the intermediate BL Lac S5 0716+714 during observations acquired by the AGILE satellite in September and October 2007. These detections of activity were contemporaneous with a period of intense optical activity, which was monitored by GASP-WEBT. This simultaneous optical and gamma-ray coverage allows us to study in detail the light curves, time lags, gamma-ray photon spectrum, and Spectral Energy Distributions (SEDs) during different states of activity. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE), which are sensitive to the 30 MeV-50 GeV and 18-60 keV energy ranges, respectively. Observations were completed in two different periods, the first between 2007 September 4-23, and the second between 2007 October 24-November 1. Over the period 2007 September 7-12, AGILE detected gamma-ray emission from the source at a significance level of 9.6-sigma with an average flux (E>100 MeV) of (97 +/- 15) x 10^{-8} photons cm^{-2} s^{-1}, which increased by a factor of at least four within three days. No emission was detected by Super-AGILE for the energy range 18-60 keV to a 3-sigma upper limit of 10 mCrab in 335 ksec. In October 2007, AGILE repointed toward S5 0716+714 following an intense optical flare, measuring an average flux of (47 +/- 11) x 10^{-8} photons cm^{-2} s^{-1} at a significance level of 6.0-sigma. The gamma-ray flux of S5 0716+714 detected by AGILE is the highest ever detected for this blazar and one of the most intense gamma-ray fluxes detected from a BL Lac object. The SED of mid-September appears to be consistent with the synchrotron self-Compton (SSC) emission model, but only by including two SSC components of different variabilities.

We report the gamma-ray activity from the Intermediate BL Lac S5 0716+714 during 2007 September-October observations by the AGILE satellite, coincident with a period of intense optical activity of the source monitored by GASP-WEBT. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE) sensitive in the energy range 30 MeV-50 GeV and 18-60 keV respectively, in two different periods: the first between 4 and 23 September 2007, the second between 24 October and 1 November 2007. Over the period 7-12 September, AGILE detected gamma-ray emission from the source at a significance level of 9.6-sigma with an average flux (E>100 MeV) of (97+/-15) x 10^{-8} photons cm^{-2} s^{-1}, increasing by a factor of at least four within three days. No emission was detected by Super-AGILE in the energy range 18-60 keV, with a 3-sigma upper limit of 10 mCrab in 335 ksec. The gamma-ray flux of S5 0716+714 detected by AGILE is the highest ever detected for this blazar and one of the most intense gamma-ray fluxes detected from a BL Lac object. The Spectral Energy Distribution (SED) of mid-September seems to be consistent with the synchrotron self-Compton (SSC) emission model, but only by including two SSC components with different variability. In October 2007 AGILE repointed toward S5 0716+714 following an intense optical flare, measuring an average flux of (47+/-11) x 10^{-8} photons cm^{-2} s^{-1} at a significance level of 6.0-sigma. The gamma-ray flux during both AGILE pointings appears to be highly variable on timescales of 1 day.

The quasar-type blazar 3C 454.3 underwent a phase of high activity in summer and autumn 2007, which was intensively monitored in the radio-to-optical bands by the Whole Earth Blazar Telescope (WEBT). The gamma-ray satellite AGILE detected this source first in late July, and then in November-December 2007. In this letter we present the multifrequency data collected by the WEBT and collaborators during the second AGILE observing period, complemented by a few contemporaneous data from UVOT onboard the Swift satellite. The aim is to trace in detail the behaviour of the synchrotron emission from the blazar jet, and to investigate the contribution from the thermal emission component. Optical data from about twenty telescopes have been homogeneously calibrated and carefully assembled to construct an R-band light curve containing about 1340 data points in 42 days. This extremely well-sampled optical light curve allows us to follow the dramatic flux variability of the source in detail. In addition, we show radio-to-UV spectral energy distributions (SEDs) at different epochs, which represent different brightness levels. In the considered period, the source varied by 2.6 mag in a couple of weeks in the R band. Many episodes of fast (i.e. intranight) variability were observed, most notably on December 12, when a flux increase of about 1.1 mag in 1.5 hours was detected, followed by a steep decrease of about 1.2 mag in 1 hour. The contribution by the thermal component is difficult to assess, due to the uncertainties in the Galactic, and possibly also intrinsic, extinction in the UV band. However, polynomial fitting of radio-to-UV SEDs reveals an increasing spectral bending going towards fainter states, suggesting a UV excess likely due to the thermal emission from the accretion disc.

AIMS: Our goal is to understand the nature of blazars and the mechanisms for the generation of high-energy gamma-rays, through the investigation of the prototypical blazar PKS 2155-304, which shows complex behaviour. METHODS: We analyze simultaneous infrared-to-X-ray observations obtained with XMM-Newton and REM on November 7, 2006, when the source was in a low X-ray state. We perform a comparative analysis of these results with those obtained from previous observations in different brightness states. RESULTS: We found that the peak of the synchrotron emission moved from ultraviolet to optical wavelengths and the X-ray spectrum is best fit with a broken power law model with Gamma_2 ~ 2.4 harder than Gamma_1 ~ 2.6 and a break at about 3.5 keV. This suggests that the soft X-rays (E < 3.5 keV) are related to the high-energy tail of the synchrotron emission, while the hard X-rays (E > 3.5 keV) are from the energy region between the synchrotron and inverse-Compton humps. The different variability at energies below and above the break strengthens this hypothesis. Our results also stress the importance of monitoring this source at both low and high energies to better characterize its variability behaviour.

Since the CGRO operation in 1991-2000, one of the primary unresolved questions about the blazar gamma-ray emission has been its possible correlation with the low-energy (in particular optical) emission. To help answer this problem, the Whole Earth Blazar Telescope (WEBT) consortium has organized the GLAST-AGILE Support Program (GASP) to provide the optical-to-radio monitoring data to be compared with the gamma-ray detections by the AGILE and GLAST satellites. This new WEBT project started in early September 2007, just before a strong gamma-ray detection of 0716+714 by AGILE. We present the GASP-WEBT optical and radio light curves of this blazar obtained in July-November 2007, about various AGILE pointings at the source. We construct NIR-to-UV spectral energy distributions (SEDs), by assembling GASP-WEBT data together with UV data from the Swift ToO observations of late October. We observe a contemporaneous optical-radio outburst, which is a rare and interesting phenomenon in blazars. The shape of the SEDs during the outburst appears peculiarly wavy because of an optical excess and a UV drop-and-rise. The optical light curve is well sampled during the AGILE pointings, showing prominent and sharp flares. A future cross-correlation analysis of the optical and AGILE data will shed light on the expected relationship between these flares and the gamma-ray events.

The blazar AO 0235+164 was claimed to show a quasi-periodic behaviour in the radio and optical bands. Moreover, an extra emission component contributing to the UV and soft X-ray flux was detected, whose nature is not yet clear. A predicted optical outburst was observed in late 2006/early 2007. We here present the radio-to-optical WEBT light curves during the outburst, together with UV data acquired by Swift in the same period. We found the optical outburst to be as strong as the big outbursts of the past: starting from late September 2006, a brightness increase of 5 mag led to the outburst peak in February 19-21, 2007. We also observed an outburst at mm and then at cm wavelengths, with an increasing time delay going toward lower frequencies during the rising phase. Cross-correlation analysis indicates that the 1 mm and 37 GHz flux variations lagged behind the R-band ones by about 3 weeks and 2 months, respectively. These short time delays suggest that the corresponding jet emitting regions are only slightly separated and/or misaligned. In contrast, during the outburst decreasing phase the flux faded contemporaneously at all cm wavelengths. This abrupt change in the emission behaviour may suggest the presence of some shutdown mechanism of intrinsic or geometric nature. The behaviour of the UV flux closely follows the optical and near-IR one. By separating the synchrotron and extra component contributions to the UV flux, we found that they correlate, which suggests that the two emissions have a common origin.

Multicolor (UBVRIJHK) observations of the blazar AO 0235+16 are analyzed. The light curves were compiled at the Turin Observatory from literature data and the results of observations obtained in the framework of the WEBT program (http://www.to.astro/blazars/webt/). The color variability of the blazar was studied in eight time intervals with a sufficient number of multicolor optical observations; JHK data are available for only one of these. The spectral energy distribution (SED) of the variable component remained constant within each interval, but varied strongly from one interval to another. After correction for dust absorption, the SED can be represented by a power law in all cases, providing evidence for a synchrotron nature of the variable component. We show that the variability at both optical and IR wavelengths is associated with the same variable source.

The blazar 3C 454.3 underwent an unprecedented optical outburst in spring 2005. This was first followed by a mm and then by a cm radio outburst, which peaked in February 2006. We report on follow-up observations by the WEBT to study the multiwavelength emission in the post-outburst phase. XMM-Newton observations on July and December 2006 added information on the X-ray and UV fluxes. The source was in a faint state. The radio flux at the higher frequencies showed a fast decreasing trend, which represents the tail of the big radio outburst. It was followed by a quiescent state, common at all radio frequencies. In contrast, moderate activity characterized the NIR and optical light curves, with a progressive increase of the variability amplitude with increasing wavelength. We ascribe this redder-when-brighter behaviour to the presence of a "little blue bump" due to line emission from the broad line region, which is clearly visible in the source SED during faint states. Moreover, the data from the XMM-Newton OM reveal a rise of the SED in the UV, suggesting the existence of a "big blue bump" due to thermal emission from the accretion disc. The X-ray spectra are well fitted with a power-law model with photoelectric absorption, possibly larger than the Galactic one. However, the comparison with previous X-ray observations would imply that the amount of absorbing matter is variable. Alternatively, the intrinsic X-ray spectrum presents a curvature, which may depend on the X-ray brightness. In this case, two scenarios are possible.

We present the results from a study of the long-term optical spectral variations of BL Lacertae, using the long and well-sampled B and R-band light curves of the Whole Earth Blazar Telescope (WEBT) collaboration, binned on time intervals of 1 day. The relation between spectral slope and flux (the spectrum gets bluer as the source flux increases) is well described by a power-law model, although there is significant scatter around the best-fitting model line. To some extent, this is due to the spectral evolution of the source (along well-defined loop-like structures) during low-amplitude events, which are superimposed on the major optical flares, and evolve on time scales of a few days. The "bluer-when-brighter" mild chromatism of the long-term variations of the source can be explained if the flux increases/decreases faster in the B than in the R band. The B and R-band variations are well correlated, with no significant, measurable delays larger than a few days. On the other hand, we find that the spectral variations lead those in the flux light curves by ~ 4 days. Our results can be explained in terms of Doppler factor variations due to changes in the viewing angle of a curved and inhomogeneous emitting jet.

At the end of July 2006, the blazar PKS 2155-304 (z=0.116) underwent a strong outburst observed at TeV energies by HESS (up to 17 Crab flux level at E > 200 GeV). The Swift satellite followed the evolution of the source for about one month. The data analysis - reported in another paper - has shown that, despite the violent activity at TeV energies, the synchrotron energy distribution increased in normalization, but only with a small shift in frequency. In the present work, we fit the broad-band spectrum with a log-parabolic model, to search for indications of intrinsic curvature, which in turn is usually interpreted as a signature of energy-dependent acceleration mechanisms of electrons.

We present all the publicly available data, from optical/UV wavelengths (UVOT) to X-rays (XRT, BAT), obtained from Swift observations of the blazar PKS 2155-304, performed in response to the rapid alert sent out after the strong TeV activity (up to 17 Crab flux level at E > 200 GeV) at the end of July 2006. The X-ray flux increased by a factor of 5 in the 0.3-10 keV energy band and by a factor of 1.5 at optical/UV wavelengths, with roughly one day of delay. The comparison of the spectral energy distribution built with data quasi-simultaneous to the TeV detections shows an increase of the overall normalization with respect to archival data, but only a small shift of the frequency of the synchrotron peak that remains consistent with the values reported in past observations when the TeV activity was much weaker.

New data and results on the optical behaviour of the blazar PKS 0735+178 are presented. In addition the whole historical light curve, and a new photometric calibration of comparison stars are reported. Optical spectral indexes are calculated and studied on years scales, while several methods for time-series analysis are applied to the whole historical series and to each observing season of our data set. This allowed to search and identify optical variability modes, characteristic timescales and the signal power spectrum over 3 decades in time. In the last 10 years the optical flux of PKS 0735+178 exhibited a rather achromatic long-term behaviour and a variability mode resembling the shot-noise. The brightness level was in an intermediate/low state with a mild flaring activity and a superimposition/succession of rapid and slower flares, with no extraordinary/isolated outbursts but, at any rate, characterized by 1 major active phase in 2001. Several mid-term scales (days, weeks) were found, the more common falling into values of about 27-28 days, 50-56 days and 76-79 days. The rapid variability in the historical curve appear to be modulated by a general, slower and rather oscillating trend, where typical timescales of about 4.5, 8.5 and 11-13 years can be identified. This spectral and temporal analysis, accompanying our data publication, suggests the occurrence of distinctive variability signatures at days/weeks scales, that can likely be of transitory nature. On the other hand the possible pseudo-cyclical or multi-component modulations at long times could be more stable, recurrent and correlated to the bimodal radio flux behaviour and the twisted radio structure observed by many years in this blazar.

We performed optical spectroscopic monitoring of the BL Lac object AO 0235+164 (z=0.94) with the VLT and TNG telescopes from Aug. 2003 to Dec. 2004, during an extended WEBT campaign. The flux of this source is both contaminated and absorbed by a foreground galactic system at z=0.524, the stars of which can act as gravitational micro-lenses. In this period the object was in an optically faint, though variable state, and a broad Mg II emission line was visible at all epochs. The spectroscopic analysis reveals an overall variation in the Mg II line flux of a factor 1.9, while the corresponding continuum flux density changed by a factor 4.3. Most likely, the photoionising radiation can be identified with the emission component that was earlier recognised to be present as a UV-soft-X-ray bump in the source SED and that is visible in the optical domain only in very faint optical states. We estimate an upper limit to the BLR size of a few light months from the historical minimum brightness level; from this we infer the maximum amplification of the Mg II line predicted by the microlensing scenario. Unless we have strongly overestimated the size of the BLR, only very massive stars could significantly magnify the broad Mg II emission line, but the time scale of variations due to these (rare) events would be of several years. In contrast, the continuum flux, coming from much smaller emission regions in the jet, could be affected by microlensing from the more plausible MACHO deflectors, with variability time scales of some months.

In spring 2005 the blazar 3C 454.3 was observed in an unprecedented bright state from the near-IR to the hard X-ray frequencies. A mm outburst peaked in June-July 2005, and it was followed by a flux increase at high radio frequencies. In this paper we report on multifrequency monitoring by the WEBT aimed at following the further evolution of the outburst in detail. In particular, we investigate the expected correlation and time delays between the optical and radio emissions in order to derive information on the variability mechanisms and jet structure. A comparison among the light curves at different frequencies is performed by means of visual inspection and discrete correlation function, and the results are interpreted with a simple model taking into account Doppler factor variations of geometric origin. The high-frequency radio light curves show a huge outburst starting during the dimming phase of the optical one and lasting more than 1 year. The first phase is characterized by a slow flux increase, while in early 2006 a major flare is observed. The lower-frequency radio light curves show a progressively delayed and fainter event, which disappears below 8 GHz. We suggest that the radio major peak is not physically connected with the spring 2005 optical one, but it is actually correlated with a minor optical flare observed in October-November 2005. This interpretation involves both an intrinsic and a geometric mechanism. The former is represented by disturbances travelling down the emitting jet, the latter being due to the curved-jet motion, with the consequent differential changes of viewing angles of the different emitting regions.

Context: Being dominated by non-thermal emission from aligned relativistic jets, blazars allow us to elucidate the physics of extragalactic jets, and, ltimately, how the energy is extracted from the central black hole in radio-loud active galactic nuclei. Aims: Crucial information is provided by broad-band spectral energy distributions (SEDs), their trends with luminosity and correlated multi-frequency variability. With this study we plan to obtain a database of contemporaneous radio-to-optical spectra of a sample of blazars, which are and will be observed by current and future high-energy satellites. Methods: Since December 2004 we are performing a monthly multi-frequency radio monitoring of a sample of 35 blazars at the antennas in Medicina and Noto. Contemporaneous near-IR and optical observations for all our observing epochs are organised. Results: Until June 2006 about 4000 radio measurements and 5500 near-IR and optical measurements were obtained. Most of the sources show significant variability in all observing bands. Here we present the multi-frequency data acquired during the first eighteen months of the project, and construct the SEDs for the best-sampled sources.

A huge multiwavelength campaign targeting the blazar AO 0235+164 was organized by the Whole Earth Blazar Telescope (WEBT) in 2003-2005 to study the variability properties of the source. Monitoring observations were carried out at cm and mm wavelengths, and in the near-IR and optical bands, while three pointings by the XMM-Newton satellite provided information on the X-ray and UV emission. We present the data acquired during the second observing season, 2004-2005, by 27 radio-to-optical telescopes. They reveal an increased near-IR and optical activity with respect to the previous season. Increased variability is also found at the higher radio frequencies, down to 15 GHz, but not at the lower ones. The radio (and optical) outburst predicted to peak around February-March 2004 on the basis of the previously observed 5-6 yr quasi-periodicity did not occur. The analysis of the optical light curves reveals now a longer characteristic time scale of 8 yr, which is also present in the radio data. The spectral energy distributions corresponding to the XMM-Newton observations performed during the WEBT campaign are compared with those pertaining to previous pointings of X-ray satellites. Bright, soft X-ray spectra can be described in terms of an extra component, which appears also when the source is faint through a hard UV spectrum and a curvature of the X-ray spectrum. Finally, there might be a correlation between the X-ray and optical bright states with a long time delay of about 5 yr, which would require a geometrical interpretation.

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

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.

(abridged) We performed a homogeneous and systematic analysis of simultaneous X-ray and optical/UV properties of a group of 15 gamma-ray loud AGN, using observations performed with XMM-Newton. The sample is composed of 13 blazars (6 BL Lac and 7 Flat-Spectrum Radio Quasar) and 2 radio galaxies, that are associated with detections at energies >100 MeV. The data for 7 of them are analyzed here for the first time, including the first X-ray observation of PKS 1406-706. Then, the spectral characteristics of the sources in the present sample have been compared with those in previous catalogs of blazars and other AGN, in order to search for any difference or long term change. All the selected sources appear to follow the classic "blazar sequence" and the spectral energy distributions (SED) built with the present X-ray and optical/UV data and completed with historical data, confirm the findings of previous studies on this type of sources. Some sources display interesting features worth noting: four of them, namely AO 0235+164, PKS 1127-145, S5 0836+710 and PKS 1830-211 show the presence of an intervening absorption system along the line of sight, but only the latter is known to be surely gravitationally lensed. AO 0235+164 was detected during an outburst and its SED shows a clear shift of the synchrotron peak. 3C 273 shows a change in the state with respect to the previous BeppoSAX observations, that can be interpreted as an increase of the Seyfert-like component and a corresponding decline of the jet emission. This is consistent with the monitoring at radio wavelengths performed during the same years. PKS 1406-706 is detected with a flux higher than in the past, but with a corresponding low optical flux. Moreover, although it is classified as FSRQ, the SED can be modelled with a simple synchrotron self-Compton model.

We report about a short flare from the blazar NRAO 530 occurred on 17 February 2004 and detected serendipitously by the IBIS/ISGRI detector on board INTEGRAL. In the 20-40 keV energy range, the source, that is otherwise below the detection limit, is detected at a level of ~2 x 10^-10 erg cm^-2 s^-1 during a time interval of less than 2000 s, which is about a factor 2 above the detection threshold. At other wavelengths, only nearly-simultaneous radio data are available (1 observation at 2 cm on 11 February 2004), indicating a moderate increase of the polarization. This appears to be the shortest time variability episode ever detected in a high luminosity blazar at hard X-rays, unless the blazar is contaminated by the presence of an unknown unresolved rapidly varying source.

The radio quasar 3C 454.3 underwent an exceptional optical outburst lasting more than 1 year and culminating in spring 2005. The maximum brightness detected was R = 12.0, which represents the most luminous quasar state thus far observed (M_B ~ -31.4). In order to follow the emission behaviour of the source in detail, a large multiwavelength campaign was organized by the Whole Earth Blazar Telescope (WEBT). Continuous optical, near-IR and radio monitoring was performed in several bands. ToO pointings by the Chandra and INTEGRAL satellites provided additional information at high energies in May 2005. The historical radio and optical light curves show different behaviours. Until about 2001.0 only moderate variability was present in the optical regime, while prominent and long-lasting radio outbursts were visible at the various radio frequencies, with higher-frequency variations preceding the lower-frequency ones. After that date, the optical activity increased and the radio flux is less variable. This suggests that the optical and radio emissions come from two separate and misaligned jet regions, with the inner optical one acquiring a smaller viewing angle during the 2004-2005 outburst. Moreover, the colour-index behaviour (generally redder-when-brighter) during the outburst suggests the presence of a luminous accretion disc. A huge mm outburst followed the optical one, peaking in June-July 2005. The high-frequency (37-43 GHz) radio flux started to increase in early 2005 and reached a maximum at the end of our observing period (end of September 2005). VLBA observations at 43 GHz during the summer confirm the

In this paper we analyse five observations of the BL Lac object AO 0235+16 performed with the Chandra and XMM-Newton satellites during the years 2000-2005. In the February 2002 observation the source is found in a bright state and presents a steep X-ray spectrum, while in all the other epochs it is faint and the spectrum is hard. The soft X-ray spectrum appears to be strongly absorbed, likely by the intervening system at z=0.524, which also absorbs the optical-UV radiation. We find that models that consider spectral curvature are superior to single power law ones in fitting the X-ray spectrum. In particular, we favour a double power law model, which agrees with the assumption of a superposition of two different components in the X-ray domain. Both in the Chandra and in one of the XMM-Newton observations, a tentative detection of the redshifted Fe Kalpha emission line may suggest its origin from the inner part of an accretion disc. Thermal emission from this accretion disc might explain the UV-soft-X-ray bump that appears in the spectral energy distributions, when the X-ray spectra are complemented with the optical-UV data from the Optical Monitor onboard XMM-Newton. More likely, the bump can be interpreted in terms of an additional synchrotron component emitted from an inner region of the jet with respect to that where the lower-energy emission comes from. An inspection of the X-ray light curves reveals that intraday variability occurs only when the source is in a bright state.

Notwithstanding the big efforts devoted to the investigation of the mechanisms responsible for the high-energy (E>100 MeV) gamma-ray emission in active galactic nuclei (AGN), the definite answer is still missing. The X-ray energy band (0.4-10 keV) is crucial for this type of study, since both synchrotron and inverse Compton emission can contribute to the formation of the continuum. Within an ongoing project aimed at the investigation of the gamma-ray emission mechanism acting in the AGN detected by the EGRET telescope onboard CGRO, we firstly focused on the sources for which X-ray and optical/UV data are available in the XMM-Newton public archive. The preliminary results are outlined here.

A multiwavelength campaign on the BL Lac object AO 0235+16 has been organized by the WEBT collaboration during the observing seasons 2003-2004 and 2004-2005. Here we report on the results of the first season, which saw the participation of 24 optical and near-IR telescopes and 4 radio telescopes, as well as the first XMM-Newton pointing, which occurred on January 18-19, 2004. Unpublished data from previous epochs were also collected. In the optical and near-IR the source has been very active in the last 3 years, though rather faint most of the time, with noticeable variations of more than a mag in a few days. In contrast, in the radio bands it appears "quiescent" since early 2000. The X-ray spectra obtained by the three EPIC detectors are well fitted by a power law with extra-absorption at z=0.524. No significant variation of the X-ray flux occurred during the pointing. In contrast, contemporaneous dense radio monitoring with the 100 m telescope at Effelsberg shows a ~2-3% flux decrease in 6-7 hours, which, if intrinsic, would imply a brightness temperature well above the Compton limit. We construct the broad-band spectral energy distribution of January 18-19, 2004 with simultaneous radio data from Effelsberg, optical data from the Nordic Optical Telescope, optical-UV data from the Optical Monitor onboard XMM-Newton, and X-ray data by the EPIC instruments. Particular care was taken to correct data for extinction due to both the Milky Way and the z=0.524 absorber. The resulting SED suggests the existence of a bump in the UV spectral region.

BL Lacertae has been the target of four observing campaigns by the Whole Earth Blazar Telescope (WEBT) collaboration. In this paper we present $UBVRI$ light curves obtained by the WEBT from 1994 to 2002, including the last, extended BL Lac 2001 campaign. A total of about 7500 optical observations performed by 31 telescopes from Japan to Mexico have been collected, to be added to the $\sim 15600$ observations of the BL Lac Campaign 2000. All these data allow one to follow the source optical emission behaviour with unprecedented detail. The analysis of the colour indices reveals that the flux variability can be interpreted in terms of two components: longer-term variations occurring on a few-day time scale appear as mildly-chromatic events, while a strong bluer-when-brighter chromatism characterizes very fast (intraday) flares. By decoupling the two components, we quantify the degree of chromatism inferring that longer-term flux changes imply moving along a $\sim 0.1$ bluer-when-brighter slope in the $B-R$ versus $R$ plane; a steeper slope of $\sim 0.4$ would distinguish the shorter-term variations. This means that, when considering the long-term trend, the $B$-band flux level is related to the $R$-band one according to a power law of index $\sim 1.1$. Doppler factor variations on a "convex" spectrum could be the mechanism accounting for both the long-term variations and their slight chromatism.

The long-term variability of the multiwavelength blazar emission can be interpreted in terms of orientation variations of a helical, inhomogeneous, non-thermally emitting jet, possibly caused by the orbital motion of the parent black hole in a binary system (Villata & Raiteri 1999). The helical-jet model is here applied to explain the quasi-periodic radio-optical light curves and the broad-band spectral energy distributions (SEDs) of the BL Lac object AO 0235+16. Through a suitable choice of the model parameters, the helix rotation can well account for the periodicity of the main radio and optical outbursts and for the corresponding SED variability, while the interspersed minor radio events could be interpreted as due either to some local distortions of the helical structure or to other phenomena contributing to the source emission. In particular, the probable existence of flow instabilities provides a viable interpretation for the non-periodic features.

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

We present the most continuous data base of optical $BVR_{c}I_{c}$ observations ever published on the BL Lacertae object GC 0109+224, collected mainly by the robotic telescope of the Perugia University Observatory in the period November 1994-February 2002. These observations have been complemented by data from the Torino Observatory, collected in the period July 1995-January 1999, and Mt. Maidanak Observatory (December 2000). GC 0109+224 showed rapid optical variations and six major outbursts were observed at the beginning and end of 1996, in fall 1998, at the beginning and at the end of 2000, and at the beginning of 2002. Fast and large-amplitude drops characterized its flux behaviour. The $R_c$ magnitude ranged from 13.3 (16.16 mJy) to 16.46 (0.8 mJy), with a mean value of 14.9 (3.38 mJy). In the periods where we collected multi-filter observations, we analyzed colour and spectral indexes, and the variability patterns during some flares. The long-term behaviour seems approximatively achromatic, but during some isolated outbursts we found evidence of the typical loop-like hysteresis behaviour, suggesting that rapid optical variability is dominated by non-thermal cooling of a single emitting particle population. We performed also a statistical analysis of the data, through the discrete correlation function (DCF), the structure function (SF), and the Lomb-Scargle periodogram, to identify characteristic times scales, from days to months, in the light curves, and to quantify the mode of variability. We also include the reconstruction of the historical light curve and a photometric calibration of comparison stars, to favour further extensive optical monitoring of this interesting blazar.

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

We present UBVRI light curves of BL Lacertae from May 2000 to January 2001, obtained by 24 telescopes in 11 countries. More than 15000 observations were performed in that period, which was the extension of the Whole Earth Blazar Telescope (WEBT) campaign originally planned for July-August 2000. Rapid flux oscillations are present all the time, involving variations up to a few tenths of mag on hour time scales, and witnessing an intense intraday activity of this source. Colour indexes have been derived by coupling the highest precision B and R data taken by the same instrument within 20 min and after subtracting the host galaxy contribution from the fluxes. The 620 indexes obtained show that the optical spectrum is weakly sensitive to the long-term trend, while it strictly follows the short-term flux behaviour, becoming bluer when the brightness increases. Thus, spectral changes are not related to the host galaxy contribution, but they are an intrinsic feature of fast flares. We suggest that the achromatic mechanism causing the long-term flux base-level modulation can be envisaged in a variation of the relativistic Doppler beaming factor, and that this variation is likely due to a change of the viewing angle. Discrete correlation function (DCF) analysis reveals the existence of a characteristic time scale of variability of about 7 h in the light curve of the core WEBT campaign, while no measurable time delay between variations in the B and R bands is found.

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

We report on two BeppoSAX observations of BL Lac (2200+420) performed respectively in June and December 1999, as part of a ToO program to monitor blazars in high states of activity. During both runs the source has been detected up to 100 keV, but it showed quite different spectra: in June it was concave with a very hard component above 5-6 keV (alpha_1 ~ 1.6; alpha_2 ~ 0.15); in December it was well fitted by a single power law (alpha ~ 0.6). During the first BeppoSAX observation BL Lac showed an astonishing variability episode: the 0.3 - 2 keV flux doubled in ~ 20 minutes, while the flux above 4 keV was almost constant. This frequency-dependent event is one of the shortest ever recordered for BL Lac objects and places lower limits on the dimension and the magnetic field of the emitting region and on the energy of the synchrotron radiating electrons. A similar but less extreme behaviour is detected also in optical light curves, that display non-simultaneous, smaller fluctuations of ~ 20% in 20 min. We fit the spectral energy distributions with a homogeneous, one-zone model to constrain the emission region in a very simple but effective SSC + external Compton scenario, highlightning the importance of the location of the emitting region with respect to the Broad Line Region and the relative spectral shape dependence. We compare our data with historical radio to gamma-ray Spectral Energy Distributions.

We present the BeppoSAX observations of 3C 279 performed in 1997 January 13-23 simultaneously with gamma-ray pointings with EGRET on board CGRO and optical observations at the Torino Observatory (1997 January 11-18). ISO data close to this epoch are also available. We compare the derived Spectral Energy Distribution with those obtained at all other epochs with adequate multiwavelength coverage. Simple spectral models fitted to the multiepoch SEDs suggest that changes in intensity can be ascribed mainly to the variation of the bulk Lorentz factor of the plasma in the jet.

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

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

We present the preliminary results of a long BeppoSAX observation of the BL Lac object Mkn501 carried out in June 1999. The source was fainter than found during the BeppoSAX pointings of 1997 and 1998, but is still detected with a good signal-to-noise ratio up to ~40 keV. The X-ray spectrum in the energy range 0.1-40 keV, produced through synchrotron radiation, is steeper than in the previous years, it is clearly curved, and peaks (in nu*F_nu) at ~0.5 keV. This energy is much lower than those at which the synchrotron component was found to peak in 1997 and 1998. Some intraday variability suggests that activity of the source on small time scales accompanies the large long time scale changes of brightness and spectrum.

We present the results of two ToO BeppoSAX observations of BL Lac. During the first observation we detected both the synchrotron and the Compton components. Fast time variability was present, but only for the synchrotron component. During the second observation the spectrum was flatter and only the Compton component was present. Four different SED are presented and can be described by a blob moving along a jet and responsible for the SSC emission. The seed photons for the Compton scattering are the synchrotron photons themself plus possibly external photons coming from the broad line region.

We present BeppoSAX observations of three gamma-ray emitting quasars, namely 0836+710, 1510-089 and 2230+114. The three objects have been detected up to 100 keV showing extremely flat power-law spectra above 2 keV (energy index alpha_{2-10}=0.3-0.5). The soft X-ray spectrum of 0836+710 implies either an absorption column density higher than the galactic one or an intrinsically very hard slope (alpha_{0.1-1}=-0.2) below 1 keV. 1510-089 shows a soft excess, with the low energy spectrum steeper (alpha_{0.1-1}=1.6) than the high energy power-law. The results are discussed in the framework of current Inverse Compton models for the high energy emission of Flat Spectrum Radio Quasars and are used to estimate the physical quantities in the jet emitting region and to shed light on the energy transport mechanisms in jets. Finally we discuss the estimates of the jet luminosity in the context of the Blandford & Znajek mechanism for jet production.

We report first observational results of multifrequency campaigns on the prominent Virgo blazars 3C 273 and 3C 279 which were carried out in January and February 1999. Both blazars are detected from radio to gamma-ray energies. We present the measured X- to gamma-ray spectra of both sources, and for 3C 279 we compare the 1999 broad-band (radio to gamma-ray) spectrum to measured previous ones.

We present $Beppo$SAX observations of the $\gamma $-ray emitting quasars 0836+710, 1510-089 and 2230+114. All the objects have been detected in the PDS up to 100 keV and have extremely flat power-law spectra above 2 keV ($\alpha _x$=0.3--0.5). 0836+710 shows absorption higher than the galactic value and marginal evidence for the presence of the redshifted 6.4 keV Iron line. 1510-089 shows a spectral break around 1 keV, with the low energy spectrum steeper ($\alpha_l$=1.6) than the high energy power-law ($\alpha_h$=0.3). The data are discussed in the light of current Inverse Compton models for the high energy emission.