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

A detailed understanding of how the activity of a galactic nucleus regulates the growth of its host is still missing. To understand the activity and the types of accretion of supermassive black holes in different hosts, it is essential to study radio-optical properties of a large sample of extragalactic sources. In particular, we aim at studying the radio spectral index trends across the optical emission line diagnostic diagrams to search for potential (anti)correlations. To this goal, we combine flux densities from the radio FIRST survey at $1.4\,{\rm GHz}$ (with the flux density range $10\,{\rm mJy} \leq F_{1.4} \leq 1000\,{\rm mJy}$) for 209 SDSS sources at intermediate redshift $(0.04\leq z \leq 0.4)$ with the Effelsberg radiotelescope measurements at $4.85\,{\rm GHz}$ and $10.45\,{\rm GHz}$. The information about the optical emission-line ratios is obtained from the SDSS-DR7 catalogue. Using the Effelsberg data, we were able to infer the two-point radio spectral index distributions for star-forming galaxies, composite galaxies (with a combined contribution to the line emission from the star-formation and AGN activity), Seyferts, and low ionization narrow emission region (LINER) galaxies. While studying the distribution of steep, flat, and inverted sources across optical diagnostic diagrams, we found three distinct classes of radio emitters for our sample: (i) sources with steep radio index, high ionization ratio and high radio loudness, (ii) sources with flat radio index, lower ionization ratio and intermediate radio loudness, (iii) sources with inverted radio index, low ionization ratio and low radio loudness. The classes (i), (ii), (iii) cluster mainly along the transition from Seyfert to LINER sources in the BPT diagram. We interpret these groups as a result of the recurrent nuclear-jet activity.

We present the design and performance of RoboPol, a four-channel optical polarimeter operating at the Skinakas Observatory in Crete, Greece. RoboPol is capable of measuring both relative linear Stokes parameters $q$ and $u$ (and the total intensity $I$) in one sky exposure. Though primarily used to measure the polarization of point sources in the R-band, the instrument features additional filters (B, V and I), enabling multi-wavelength imaging polarimetry over a large field of view (13.6' $\times$ 13.6'). We demonstrate the accuracy and stability of the instrument throughout its five years of operation. Best performance is achieved within the central region of the field of view and in the R band. For such measurements the systematic uncertainty is below 0.1% in fractional linear polarization, $p$ (0.05% maximum likelihood). Throughout all observing seasons the instrumental polarization varies within 0.1% in $p$ and within 1$^\circ$ in polarization angle.

The advent of the Fermi-GST with its unprecedented capability to monitor the entire 4 pi sky within less than 2-3 hours, introduced new standard in time domain gamma-ray astronomy. To explore this new avenue of extragalactic physics the F-GAMMA programme undertook the task of conducting nearly monthly, broadband radio monitoring of selected blazars from January 2007 to January 2015. In this work we release all the light curves at 2.64, 4.85, 8.35, 10.45, 14.6, 23.05, 32, and 43 GHz and present first order derivative data products after all necessary post-measurement corrections and quality checks; that is flux density moments and spectral indices. The release includes 155 sources. The effective cadence after the quality flagging is around one radio SED every 1.3 months. The coherence of each radio SED is around 40 minutes. The released dataset includes more than $4\times10^4$ measurements. The median fractional error at the lowest frequencies (2.64-10.45 GHz) is below 2%. At the highest frequencies (14.6-43 GHz) with limiting factor of the atmospheric conditions, the errors range from 3% to 9%, respectively.

Black holes, anywhere in the stellar-mass to supermassive range, are often associated with relativistic jets. Models suggest that jet production may be a universal process common in all black hole systems regardless of their mass. Although in many cases observations support such hypotheses for microquasars and Seyfert galaxies, little is known on whether boosted blazar jets also comply with such universal scaling laws. We use uniquely rich multiwavelength radio light curves from the F-GAMMA program and the most accurate Doppler factors available to date to probe blazar jets in their emission rest frame with unprecedented accuracy. We identify for the first time a strong correlation between the blazar intrinsic broad-band radio luminosity and black hole mass, which extends over $\sim$ 9 orders of magnitude down to microquasars scales. Our results reveal the presence of a universal scaling law that bridges the observing and emission rest frames in beamed sources and allows us to effectively constrain jet models. They consequently provide an independent method for estimating the Doppler factor, and for predicting expected radio luminosities of boosted jets operating in systems of intermediate or tens-of-solar mass black holes, immediately applicable to cases as those recently observed by LIGO.

LS I +61${^\circ}$303 is a $\gamma$-ray emitting X-ray binary with periodic radio outbursts with time scales of one month. Previous observations have revealed microflares superimposed on these large outbursts with periods ranging from a few minutes to hours. This makes LS I +61${^\circ}$303, along with Cyg X-1, the only TeV emitting X-ray binary exhibiting radio microflares. To further investigate these microflaring activity in LS I +61${^\circ}$303 we observed the source with the 100-m Effelsberg radio telescope at 4.85, 8.35, and 10.45 GHz and performed timing analysis on the obtained data. Radio oscillations of 15 hours time scales are detected at all three frequencies. We also compare the spectral index evolution of radio data to that of the photon index of GeV data observed by Fermi-LAT. We conclude that the observed QPO could result from multiple shocks in a jet.

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

Recent population studies have shown that the variability Doppler factors can adequately describe blazars as a population. We use the flux density variations found within the extensive radio multi-wavelength datasets of the F-GAMMA program, a total of 10 frequencies from 2.64 up to 142.33 GHz, in order to estimate the variability Doppler factors for 58 $\gamma$-ray bright sources, for 20 of which no variability Doppler factor has been estimated before. We employ specifically designed algorithms in order to obtain a model for each flare at each frequency. We then identify each event and track its evolution through all the available frequencies for each source. This approach allows us to distinguish significant events producing flares from stochastic variability in blazar jets. It also allows us to effectively constrain the variability brightness temperature and hence the variability Doppler factor as well as provide error estimates. Our method can produce the most accurate (16\% error on average) estimates in the literature to date.

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.

Radio polarimetry is an invaluable tool to investigate the physical conditions and variability processes in active galactic nuclei (AGN) jets. However, detecting their linear and circular polarization properties is a challenging endeavor due to their low levels and possible depolarization effects. We have developed an end-to-end data analysis methodology to recover the polarization properties of unresolved sources with high accuracy. It has been applied to recover the linear and circular polarization of 87 AGNs measured by the F-GAMMA program from July 2010 to January 2015 with a mean cadence of 1.3 months. Their linear polarization was recovered at four frequencies between 2.64 and 10.45 GHz and the circular polarization at 4.85 and 8.35 GHz. The physical conditions required to reproduce the observed polarization properties and the processes which induce their variability were investigated with a full-Stokes radiative transfer code which emulates the synchrotron emission of modeled jets. The model was used to investigate the conditions needed to reproduce the observed polarization behavior for the blazar 3C 454.3, assuming that the observed variability is attributed to evolving internal shocks propagating downstream.

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

Near-Eddington accretion rates onto low-mass black holes are thought to be a prime driver of the multi-wavelength properties of the narrow-line Seyfert 1 (NLS1) population of active galactic nuclei (AGN). Orientation effects have repeatedly been considered as another important factor involved, but detailed studies have been hampered by the lack of measured viewing angles towards this type of AGN. Here we present multi-epoch, 15 GHz VLBA images (MOJAVE program) of the radio-loud and Fermi/LAT-detected NLS1 galaxy 1H 323+342. These are combined with single-dish, multi-frequency radio monitoring of the source's variability, obtained with the Effelsberg 100-m and IRAM 30-m telescopes, in the course of the F-GAMMA program. The VLBA images reveal 6 components with apparent peeds of ~1 to ~7 c, and one quasi-stationary feature. Combining the obtained apparent jet speed ($\beta_{app}$) and variability Doppler factor ($D_{var}$) estimates together with other methods, we constrain the viewing angle towards 1H 0323+342 to $\theta \leq 4 - 13$ deg. Using literature values of $\beta_{app}$ and $D_{var}$, we also deduce a viewing angle of $\leq$ 8-9 deg towards another radio- and {\gamma}-ray-loud NLS1, namely SBS 0846+513.

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

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

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

Detailed investigation of broadband flux variability in the blazar 3C 273 allows us to probe the location and size of emission regions and their physical conditions. We report the results on correlation studies of the flaring activity observed between 2008 and 2012. The observed broadband variations were investigated using the structure function and the discrete correlation function, and power spectral density analysis (PSD) methods. The PSD analysis showed that the optical/IR light curve slopes are consistent with the slope of white noise processes, while, the PSD slopes at radio, X-ray and gamma-ray energies are consistent with red-noise processes. The flux variations at gamma-ray and mm-radio bands are found to be significantly correlated. Using the estimated time lag of (110\pm27) days between gamma-ray and radio light curves, we constrained the location of the gamma-ray emission region at a de-projected distance of 1.2\pm0.9 pc from the jet apex. Flux variations at X-ray bands were found to have a significant correlation with variations at both radio and \gamma-rays energies. The correlation between X-rays and gamma-rays light curves suggests presence of two components responsible for the X-ray emission. A negative time lag of -(50\pm20) days, where the X-rays are leading the emission, suggests X-rays are emitted closer to the jet apex from a compact region at a distance of ~(0.5\pm0.4) pc from the jet apex. A positive time lag of (110\pm20) days suggests jet-base origin of the other X-ray component at ~(4--5)~pc from the jet apex. The flux variations at radio frequencies were found to be well correlated with each other such that the variations at higher frequencies are leading the lower frequencies, which could be expected in the standard shock-in-jet model.

Context. The origin of blazar variability, seen from radio up to gamma rays, is still a heavily debated matter and broadband flares offer a unique testbed towards a better understanding of these extreme objects. Such an energetic outburst was detected by Fermi/LAT in 2008 from the blazar PKS 1502+106. The outburst was observed from gamma rays down to radio frequencies. Aims. Through the delay between flare maxima at different radio frequencies, we study the frequency-dependent position of the unit-opacity surface and infer its absolute position with respect to the jet base. This nuclear opacity profile enables the magnetic field tomography of the jet. We also localize the gamma-ray emission region and explore the mechanism producing the flare. Methods. The radio flare of PKS 1502+106 is studied through single-dish flux density measurements at 12 frequencies in the range 2.64 to 226.5 GHz. To quantify it, we employ both a Gaussian process regression and a discrete cross-correlation function analysis. Results. We find that the light curve parameters (flare amplitude and cross-band delays) show a power-law dependence on frequency. Delays decrease with frequency, and the flare amplitudes increase up to about 43 GHz and then decay. This behavior is consistent with a shock propagating downstream the jet. The self-absorbed radio cores are located between about 10 and 4 pc from the jet base and their magnetic field strengths range between 14 and 176 mG, at the frequencies 2.64 to 86.24 GHz. Finally, the gamma-ray active region is located at (1.9 +/- 1.1) pc away from the jet base.

We present the outline and first results of a project using the synergies of the long term blazar radiomillimetre monitoring program F-GAMMA, the continued scanning of the millimetre-submillimetre sky by the Planck satellite, together with several dedicated observing programs at the Effelsberg 100m telescope, to obtain a data sample unprecedented in both time resolution and frequency span.

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

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

Context. Blazars are among the most energetic objects in the Universe. In 2008 August, Fermi/LAT detected the blazar PKS 1502+106 showing a rapid and strong gamma-ray outburst followed by high and variable flux over the next months. This activity at high energies triggered an intensive multi-wavelength campaign covering also the radio, optical, UV, and X-ray bands indicating that the flare was accompanied by a simultaneous outburst at optical/UV/X-rays and a delayed outburst at radio bands. Aims: In the current work we explore the phenomenology and physical conditions within the ultra-relativistic jet of the gamma-ray blazar PKS 1502+106. Additionally, we address the question of the spatial localization of the MeV/GeV-emitting region of the source. Methods: We utilize ultra-high angular resolution mm-VLBI observations at 43 and 86 GHz complemented by VLBI observations at 15 GHz. We also employ single-dish radio data from the F-GAMMA program at frequencies matching the VLBI monitoring. Results: PKS 1502+106 shows a compact core-jet morphology and fast superluminal motion with apparent speeds in the range 5--22 c. Estimation of Doppler factors along the jet yield values between ~7 up to ~50. This Doppler factor gradient implies an accelerating jet. The viewing angle towards the source differs between the inner and outer jet, with the former at ~3 degrees and the latter at ~1 degree, after the jet bends towards the observer beyond 1 mas. The de-projected opening angle of the ultra-fast, magnetically-dominated jet is found to be (3.8 +/- 0.5) degrees. A single jet component can be associated with the pronounced flare both at high-energies and in radio bands. Finally, the gamma-ray emission region is localized at less than 5.9 pc away from the jet base.

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

Aims: Our aim is to explore the broad-band radio continuum spectrum of LSI+61303 during its outbursts by employing the available set of secondary focus receivers of the Effelsberg 100 m telescope. Methods: The clear periodicity of the system LSI+61303 allowed observations to be scheduled covering the large radio outburst in March-April 2012. We observed LSI+61303 on 14 consecutive days at 2.6, 4.85, 8.35, 10.45, 14.3, 23, and 32 GHz with a cadence of about 12 hours followed by two additional observations several days later. Based on these observations we obtained a total of 24 quasi-simultaneous broad-band radio spectra. Results: During onset, the main flare shows an almost flat broad-band spectrum, most prominently seen on March 27, 2012, where - for the first time - a flat spectrum (alpha=0.00+/-0.07, S nu^alpha) is observed up to 32 GHz (9 mm wavelength). The flare decay phase shows superimposed 'sub-flares' with the spectral index oscillating between -0.4 and -0.1 in a quasi-regular fashion. Finally, the spectral index steepens during the decay phase, showing optically thin emission with values alpha $\sim$ -0.5 to -0.7. Conclusions: The radio characteristics of LSI+61303 compare well with those of the microquasars XTE J1752-223 and Cygnus X-3. In these systems the flaring phase is actually also composed of a sequence of outbursts with clearly different spectral characteristics: a first outburst with a flat/inverted spectrum followed by a bursting phase of optically thin emission.

As a radio-loud narrow-line Seyfert 1 galaxy (NLS1) detected by Fermi/LAT in GeV $\gamma$-rays, 1H 0323+342 is a remarkable Active Galactic Nucleus (AGN) showing properties characteristic of both NLS1s and blazars. Here we present results of simultaneous X-ray and UV/optical monitoring observations on 1H 0323+342 taken with the UV/Optical Telescope (UVOT) and X-ray Telescope (XRT) onboard the Swift satellite over six years from 2006. Overall, the object showed statistically correlated variations in both the UV and X-ray bands on timescales of years as well as on timescales of days. A deep Suzaku observation reveals X-ray variability on timescales as short as a few tens of thousand seconds, and an X-ray spectrum typical of Seyfert galaxies. The broad-band spectral energy distribution, for which the data of UV and X-ray observations taken on 2009 July 26-27 were used, can be well modeled with a simple one-zone leptonic jet model plus accretion disk/corona emission. The latter is predominantly responsible for the UV/optical and X-ray (0.3-10 keV) emission and their observed variations. The correlated UV-X-ray variability on the timescale of days is consistent with reprocessing of the X-ray radiation by the accretion disk. The shortest timescale and large normalized excess variance of the X-ray variability detected with Suzaku suggest a relatively small black hole mass of the order of $10^7M_{\odot}$, consistent with the estimation based on the broad H$\beta$ line in the optical.

The IDV observations of S4 0917+624 were carried out monthly, from August 2005 to January 2010, with the Urumqi 25m radio telescope at 4.8 GHz. The quasar S4 0917+624 exhibits only very weak or no IDV during our 4.5 year observing interval. Prior to the year 2000, the source S4 0917+624 was one of the most prominent IDV sources. Our new data indicate that the previous strong IDV has ceased. We analyzed the long-term VLBI structural variability using Gaussian model-fitting. From this we obtained the flux densities and the deconvolved sizes of core and inner-jet components of the source. We studied the properties such as core fraction, angular size, spectral index, and brightness temperature of VLBI core for S4 0917+624, as well as the time delay between 5 and 15 GHz variations, and compared them with the IDV properties of S4 0917+624. The source shows ejection of several jet components that are suspected to have partially reduced the IDV amplitude of S4 0917+624. However, during 2005-2006, the VLBI core size was comparable to the size before the year 2000, but no strong IDV was detected in the period, suggesting that the quenching effect due to source size changes may not be responsible for the lack of strong IDV after the year 2000. The refractive scattering properties for the strong IDV phase of S4 0917+624 before the year 2000 are discussed. The disappearance of strong IDV in S4 0917+624 after the year 2000 is a mystery and cannot be explained via the quenching effect by changes in the observable VLBI structure. However, it may be caused by changes in the interstellar medium, i.e. by interstellar weather, which induces changes in the scintillation pattern on timescales of several years. Further coordinated multi-frequency observations will be required to distinguish between the effect of source-intrinsic variability and changing properties of the interstellar medium.

Blazars are among the most variable objects in the universe. They feature energetic jets of plasma that launch from the cores of these active galactic nuclei (AGN), triggering activity from radio up to gamma-ray energies. Spatial localization of the region of their MeV/GeV emission is a key question in understanding the blazar phenomenon. The flat spectrum radio quasar (FSRQ) PKS 1502+106 has exhibited extreme and correlated, radio and high-energy activity that triggered intense monitoring by the Fermi-GST AGN Multi-frequency Monitoring Alliance (F-GAMMA) program and the Global Millimeter VLBI Array (GMVA) down to $\lambda$3 mm (or 86 GHz), enabling the sharpest view to date towards this extreme object. Here, we report on preliminary results of our study of the gamma-ray loud blazar PKS 1502+106, combining VLBI and single dish data. We deduce the critical aspect angle towards the source to be $\theta_{\rm c} = 2.6^{\circ}$, calculate the apparent and intrinsic opening angles and constrain the distance of the 86 GHz core from the base of the conical jet, directly from mm-VLBI but also through a single dish relative timing analysis. Finally, we conclude that gamma rays from PKS 1502+106 originate from a region between ~1-16 pc away from the base of the hypothesized conical jet, well beyond the bulk of broad-line region (BLR) material of the source.

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

We present a high-frequency very long baseline interferometry (VLBI) kinematical study of the BL Lac object S5 0716+714 over the time period of September 2008 to October 2010. The aim of the study is to investigate the relation of the jet kinematics to the observed broadband flux variability. We find significant non-radial motions in the jet outflow of the source. In the radial direction, the highest measured apparent speed is \sim37 c, which is exceptionally high, especially for a BL Lac object. Patterns in the jet flow reveal a roughly stationary feature \sim0.15 mas downstream of the core. The long-term fits to the component trajectories reveal acceleration in the sub-mas region of the jet. The measured brightness temperature, T_{B}, follows a continuous trend of decline with distance, T_B \propto r_{jet}^{-(2.36\pm0.41)}, which suggests a gradient in Doppler factor along the jet axis. Our analysis suggest that a moving disturbance (or a shock wave) from the base of the jet produces the high-energy (optical to \gamma-ray) variations upstream of the 7 mm core, and then later causes an outburst in the core. Repetitive optical/\gamma-ray flares and the curved trajectories of the associated components suggest that the shock front propagates along a bent trajectory or helical path. Sharper \gamma-ray flares could be related to the passage of moving disturbances through the stationary feature. Our analysis suggests that the \gamma-ray and radio emission regions have different Doppler factors.

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

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 studied the radio emission from four radio-loud and gamma-ray-loud narrow-line Seyfert 1 galaxies. The goal was to investigate whether a relativistic jet is operating at the source, and quantify its characteristics. We relied on the most systematic monitoring of such system in the cm and mm radio bands which is conducted with the Effelsberg 100 m and IRAM 30 m telescopes and covers the longest time-baselines and the most radio frequencies to date. We extract variability parameters and compute variability brightness temperatures and Doppler factors. The jet powers were computed from the light curves to estimate the energy output. The dynamics of radio spectral energy distributions were examined to understand the mechanism causing the variability. All the sources display intensive variability that occurs at a pace faster than what is commonly seen in blazars. The flaring events show intensive spectral evolution indicative of shock evolution. The brightness temperatures and Doppler factors are moderate, implying a mildly relativistic jet. The computed jet powers show very energetic flows. The radio polarisation in one case clearly implies a quiescent jet underlying the recursive flaring activity. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are powered by relativistic jets.

Blazars are among the most powerful extragalactic objects, as a sub-class of active galactic nuclei. They launch relativistic jets and their emitted radiation shows strong variability across the entire electro-magnetic spectrum. The mechanisms producing the variability are still controversial and different models have been proposed to explain the observed variations in multi-frequency blazar light curves.We investigate the capabilities of the classical shock-in-jet model to explain and reconstruct the observed evolution of flares in the turnover frequency turnover flux density plane and their frequency-dependent light curve parameters. With a detailed parameter space study we provide the framework for future, detailed comparisons of observed flare signatures with the shock-in-jet scenario. Based on the shock model we compute synthetic single-dish light curves at different radio frequencies (2.6 to 345 GHz) and for different physical conditions in a conical jet (e.g. magnetic field geometry and Doppler factor). From those we extract the slopes of the different energy loss stages within the $\nu_\mathrm{m}$-$S_\mathrm{m}$ plane and deduce the frequency-dependence of different light curve parameters such as flare amplitude, time scale and cross-band delays. The evolution of the Doppler factor along the jet has the largest influence on the evolution of the flare and on the frequency-dependent light curve parameters. The synchrotron stage can be hidden in the Compton or in the adiabatic stage, depending mainly on the evolution of the Doppler factor, which makes it difficult to detect its signature in observations. In addition, we show that the time lags between different frequencies can be used as an efficient tool to better constrain the physical properties of these objects.

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.

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

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

The discovery of the M-sigma relation, the local galaxy bimodality, and the link between black-hole and host-galaxy properties, have raised the question whether AGN play a role in galaxy evolution. Several theoretical models implement AGN feedback to explain the observed galaxy luminosity function, and possibly the color and morphological transformation of spiral galaxies into passive ellipticals. To understand the importance of AGN feedback, a study of the AGN populations in the radio-optical domain is crucial. A mass sequence linking star-forming galaxies and AGN has been already noted in previous works, and it is now investigated as possible evolutionary sequence. We observed a sample of 119 intermediate-redshift (0.04<z<0.4) SDSS-FIRST radio emitters with the Effelsberg 100-m telescope at 4.85 and 10.45 GHz and obtained spectral indices. We find indications of spectral index flattening in high-metallicity star-forming galaxies, composite galaxies, and Seyferts. This "flattening sequence" along the [NII]-based emission-line diagnostic diagram is consistent with the hardening of galaxy ionizing field, due to nuclear activity. After combining our data with FIRST measurements at 1.4 GHz, we find that the three-point radio spectra of Seyferts and LINERs show substantial differences, attributable to small radio core components and larger (arcsecond sized) jet/lobe components, respectively. A visual inspection of FIRST images seems to confirm this hypothesis. Galaxies along this sequence are hypothesized to be transitioning from the active star-forming galaxies (blue cloud) to the passive elliptical galaxies (red sequence). This supports the suggestion that AGN play a role in shutting down star-formation, and allow the transition from one galaxy class to the other.

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light-curves for the brightest detected blazars from the first 3 years of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 years of 15GHz observations from the OVRO 40-m monitoring program. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light-curves. Modeling the light-curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high quality data in both bands shows correlations with significance larger than 3-sigma (AO 0235+164), with only two more larger than even 2.25-sigma (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in July-September 2012. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light-curves even when many years of data are used. This should be a caution. In all four sources the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.

Context. Blazars are a subset of active galactic nuclei (AGN) with jets that are oriented along our line of sight. Variability and spectral energy distribution (SED) studies are crucial tools for understanding the physical processes responsible for observed AGN emission. Aims. We report peculiar behaviour in the bright gamma-ray blazar PKS 1424-418 and use its strong variability to reveal information about the particle acceleration and interactions in the jet. Methods. Correlation analysis of the extensive optical coverage by the ATOM telescope and nearly continuous gamma-ray coverage by the Fermi Large Area Telescope is combined with broadband, time-dependent modeling of the SED incorporating supplemental information from radio and X-ray observations of this blazar. Results. We analyse in detail four bright phases at optical-GeV energies. These flares of PKS 1424-418 show high correlation between these energy ranges, with the exception of one large optical flare that coincides with relatively low gamma-ray activity. Although the optical/gamma-ray behaviour of PKS 1424-418 shows variety, the multiwavelength modeling indicates that these differences can largely be explained by changes in the flux and energy spectrum of the electrons in the jet that are radiating. We find that for all flares the SED is adequately represented by a leptonic model that includes inverse Compton emission from external radiation fields with similar parameters. Conclusions. Detailed studies of individual blazars like PKS 1424-418 during periods of enhanced activity in different wavebands are helping us identify underlying patterns in the physical parameters in this class of AGN.

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.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

Radio emission in blazars -- the aligned subset of Active Galactic Nuclei (AGN) -- is produced by synchrotron electrons moving relativistically in their jet's magnetic field. Under the assumption of some degree of uniformity of the field, the emission can be highly polarized -- linearly and circularly. In the radio regime, the observed variability is in most of the cases attributed to flaring events undergoing opacity evolution, i.e. transitions from optically thick to thin emission (or vice versa). These transistions have a specific signature in the polarization parameter space (angle and magnitude) which can be traced with high cadence polarization monitoring and provide us with a unique probe of the microphysics of the emitting region. Here we present the full Stokes analysis of radio emission from blazars observed in the framework of the F-GAMMA program and discuss the case study of PKS\,1510$-$089 which has shown a prominent polarization event around MJD 55900.

The exact location of the gamma-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and gamma-ray (0.1-300 GeV) ~ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and gamma-ray correlations (radio lagging gamma rays) when averaging over the whole sample, (ii) average time delays (source frame: 76+/-23 to 7+/-9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence tau_r,gamma (nu) ~ nu^-1, in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk gamma-ray production region typically located within/upstream of the 3 mm core region (tau_3mm,gamma=12+/-8 days), (iv) mean distances between the region of gamma-ray peak emission and the radio "tau=1 photosphere" decreasing from 9.8+/-3.0 pc (11 cm) to 0.9+/-1.1 pc (2 mm) and 1.4+/-0.8 pc (0.8 mm), (v) 3 mm/gamma-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 times 10^-6), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk gamma-ray production region from the supermassive black hole (SMBH) of ~ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm tau=1 surface at ~ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

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

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

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.

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

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

In September 2012, the high-synchrotron-peaked (HSP) blazar Markarian 421 underwent a rapid wideband radio flare, reaching nearly twice the brightest level observed in the centimeter band in over three decades of monitoring. In response to this event we carried out a five epoch centimeter- to millimeter-band multifrequency Very Long Baseline Array (VLBA) campaign to investigate the aftermath of this emission event. Rapid radio variations are unprecedented in this object and are surprising in an HSP BL Lac object. In this flare, the 15 GHz flux density increased with an exponential doubling time of about 9 days, then faded to its prior level at a similar rate. This is comparable with the fastest large-amplitude centimeter-band radio variability observed in any blazar. Similar flux density increases were detected up to millimeter bands. This radio flare followed about two months after a similarly unprecedented GeV gamma-ray flare (reaching a daily E>100 MeV flux of (1.2 +/- 0.7)x10^(-6) ph cm^(-2) s^(-1)) reported by the Fermi Large Area Telescope (LAT) collaboration, with a simultaneous tentative TeV detection by ARGO-YBJ. A cross-correlation analysis of long-term 15 GHz and LAT gamma-ray light curves finds a statistically significant correlation with the radio lagging ~40 days behind, suggesting that the gamma-ray emission originates upstream of the radio emission. Preliminary results from our VLBA observations show brightening in the unresolved core region and no evidence for apparent superluminal motions or substantial flux variations downstream.

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 1156+295 (4C +29.45) is one of the targets in the Urumqi monitoring program which aimed to search for evidence of annual modulation in the timescales of intra-day variable (IDV) sources. The IDV observations of 1156+295 were carried out nearly monthly from October 2007 to October 2009, with the Urumqi 25m radio telescope at 4.8 GHz. The source has shown prominent IDV of total flux density in most observing sessions with variability timescales of <= 1 day at 4.8 GHz. The estimated IDV timescales seem to follow an annual cycle that can be fitted with an anisotropic interstellar scintillation (ISS) model suggesting that a significant part of the flux density variations is due to ISS. The source underwent a dramatic flare in 2008. We studied the possible consequences of the flare on the IDV of 1156+295 by comparing the changes in its IDV characteristics with the evolution of the 43 GHz Very Long Baseline Array (VLBA) core size of the source. Conclusions: The quasar 1156+295 shows evidence for an annual modulation of its IDV timescales at 4.8 GHz, the ISS-induced IDV timescales and variability strength might be affected by the overall activity state of the source core. More frequent IDV and VLBI measurements are required to confirm the relation between the IDV appearance and the core-size evolution of the source.

We present a \gamma-ray photon flux and spectral variability study of the flat-spectrum radio quasar 3C 273 over a rapid flaring activity period between September 2009 to April 2010. Five major flares are observed in the source during this period. The most rapid flare observed in the source has a flux doubling time of 1.1 hr. The rapid \gamma-ray flares allow us to constrain the location and size of the \gamma-ray emission region in the source. The \gamma \gamma-opacity constrains the Doppler factor, $\delta_{\gamma} \geq$ 10 for the highest energy (15 GeV) photon observed by the {\it Fermi}-Large Area Telescope (LAT). Causality arguments constrain the size of the emission region to 1.6$\times 10^{15}$ cm. The \gamma-ray spectra measured over this period show clear deviations from a simple power law with a break in 1-2 GeV energy range. We discuss possible explanations for the origin of the \gamma-ray spectral breaks. Our study suggests that the \gamma-ray emission region in 3C 273 is located within the broad line region ($<$1.6 pc). The spectral behavior and temporal characteristics of the individual flares indicate the presence of multiple shock scenarios at the base of the jet.

Aims. The observation of gamma -ray flares from blazar 0836+710 in 2011, following a period of quiescence, offered an opportunity to study correlated activity at different wavelengths for a high-redshift (z=2.218) active galactic nucleus. Methods. Optical and radio monitoring, plus Fermi-LAT gamma-ray monitoring provided 2008-2012 coverage, while Swift offered auxiliary optical, ultraviolet, and X-ray information. Other contemporaneous observations were used to construct a broad-band spectral energy distribution. Results. There is evidence of correlation but not a measurable lag between the optical and gamma-ray flaring emission. On the contrary, there is no clear correlation between radio and gamma-ray activity, indicating radio emission regions that are unrelated to the parts of the jet that produce the gamma-rays. The gamma-ray energy spectrum is unusual in showing a change of shape from a power law to a curved spectrum when going from the quiescent state to the active state.

Narrow-Line Seyfert 1 galaxies have been established as a new class of gamma-ray emitting AGN with relatively low black hole masses, but near-Eddington accretion rates. Other extragalactic gamma-ray sources observed so far such as Flat Spectrum Radio Quasars, Radio Galaxies, and BL Lacertae Objects generally exhibit much higher black hole masses and, in the case of BL Lac objects and FRI Radio Galaxies, much lower accretions rates. The multifrequency campaign of 2013 centered on the bright source 1H 0323+342 will provide further insights into the nature of the jets and their gamma ray production mechanisms in a largely unexplored corner of AGN parameter space. Here, we show preliminary results of this campaign and discuss them.

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

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

The detection of gamma rays from a small number of Narrow Line Seyfert 1 galaxies by the LAT instrument onboard Fermi seriously challenged our understanding of AGN physics. Among the most important findings associated with their discovery has been the realisation that smaller-mass black holes seem to be hosted by these systems. Immediately after their discovery a radio multi- frequency monitoring campaign was initiated to understand their jet radio emission. Here the first results of the campaign are presented. The light curves and some first variability analyses are discussed, showing that the brightness temperatures and Doppler factors are moderate. The phenomenologies are typically blazar-like. The frequency domain on the other hand indicates intense spectral evolution and the variability patterns indicate mechanisms similar to those acting in the jets of BL Lacs and FSRQs. Finally, the linear polarisation also reveals the presence of a quiescent, optically thin jet in certain cases.

Densely time sampled multi-frequency flux measurements of the extreme BL Lac object S5 0716+714 over the past three years allow us to study its broad-band variability, and the detailed underlying physics, with emphasis on the location and size of the emitting regions and the evolution with time. We study the characteristics of some prominent mm-/gamma-ray flares in the context of the shock-in-jet model and investigate the location of the high energy emission region. The rapid rise and decay of the radio flares is in agreement with the formation of a shock and its evolution, if a geometrical variation is included in addition to intrinsic variations of the source. We find evidence for a correlation between flux variations at gamma-ray and radio frequencies. A two month time-delay between gamma-ray and radio flares indicates a non-cospatial origin of gamma-rays and radio flux variations in S5 0716+714.

The GeV observations by Fermi-LAT give us the opportunity to characterize the high-energy emission (100 MeV - 300 GeV) variability properties of the BL Lac object S5 0716+714. In this study, we performed flux and spectral analysis of more than 3 year long (August 2008 to April 2012) Fermi-LAT data of the source. During this period, the source exhibits two different modes of flux variability with characteristic timescales of ~75 and ~140 days, respectively. We also notice that the flux variations are characterized by a weak spectral hardening. The GeV spectrum of the source shows a clear deviation from a simple power law, and is better explained by a broken power law. Similar to other bright Fermi blazars, the break energy does not vary with the source flux during the different activity states. We discuss several possible scenarios to explain the observed spectral break.

Combining measurements taken using the Wilkinson Microwave Anisotropy Probe (WMAP) from 2001 to 2008 with measurements taken using Planck from 2009 to 2010, we investigate the long-term flux density variability of extragalactic radio sources selected from the Planck Early Release Compact Source Catalogue. The single-year, single-frequency WMAP maps are used to estimate yearly-averaged flux densities of the sources in the four WMAP bands: Ka (33 GHz), Q (41 GHz), V (61 GHz), and W (94 GHz). We identify 82, 67, 32, and 15 sources respectively as variable at greater than 99% confidence level in these four bands. The amplitudes of variation are comparable between bands, and are not correlated with either the flux densities or the spectral indices of the sources. The number counts of WMAP Ka-band sources are stable from year to year despite the fluctuation caused by individual source variability. Most of our sources show strong correlation in variability between bands. Almost all the sources that show variability are blazars. We have attempted to fit two simple, four-parameter models to the time-series of 32 sources showing correlated variability at multiple frequencies - a long-term flaring model and a rotating-jet model. We find that 19 sources (60%) can be fit with the simple rotating-jet model, and ten of these also fit the simple long-term flaring model. The remaining 13 sources (40%) show more complex variability behaviour that is not consistent with either model. Extended radio galaxies in our sample show no sign of variability, as expected, with the exception of Pictor A for which we report evidence for a millimetre flare lasting between 2002 and 2010.

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

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

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

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

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high.

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.

The BLLac object S4 0954+65 is one of the main targets of the Urumqi monitoring program targeting IntraDay Variable (IDV) sources. Between August 2005 and December 2009, the source was included in 41 observing sessions, carried out at a frequency of 4.8 GHz. The time analysis of the collected light curves, performed by applying both a structure function analysis and a specifically developed wavelet-based algorithm, discovered an annual cycle in the variability timescales, suggesting that there is a fundamental contribution by interstellar scintillation to the IDV pattern of the source. The combined use of the two analysis methods also revealed that there was a dramatic change in the variability characteristics of the source between February and March 2008, at the starting time of a strong outburst phase. The analysis' results suggest that the flaring state of the source coincides with the appearance of multiple timescales in its light curves, indicating that changes in the structure of the relativistically moving emitting region may strongly influence the variability observed on IDV timescales.

So far, no systematic long-term blazar monitoring programs and detailed variability studies exist at sub-mm wavelengths. Here, we present a new sub-mm blazar monitoring program using the APEX 12-m telescope. A sample of about 40 gamma-ray blazars has been monitored since 2007/2008 with the LABOCA bolometer camera at 345 GHz. First light curves, preliminary variability results and a first comparison with the longer cm/mm bands (F-GAMMA program) are presented, demonstrating the extreme variability characteristics of blazars at such short wavelengths.

We aim to search for evidence of annual modulation in the time scales of the BL Lac object S5 0716+714. The intra-day variability (IDV) observations were carried out monthly from 2005 to 2009, with the Urumqi 25m radio telescope at 4.8 GHz. The source has shown prominent IDV as well as long-term flux variations. The IDV time scale does show evidence in favor of an annual modulation, suggesting that the IDV of 0716+714 is dominated by interstellar scintillation. The source underwent a strong outburst phase between mid-2008 and mid-2009; a second intense flare was observed in late 2009, but no correlation between the total flux density and the IDV time scale is found, implying that the flaring state of the source does not have serious implications for the general characteristics of its intra-day variability. However, we find that the inner-jet position angle is changing throughout the years, which could result in an annual modulation noise in the anisotropic ISS model fit. There is also an indication that the lowest IDV amplitudes (rms flux density) correspond to the slowest time scales of IDV, which would be consistent with an ISS origin of the IDV of 0716+714.

The F-GAMMA program is among the most comprehensive programs that aim at understanding the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are presented. (a) We discuss that the variability patterns traced can be classified into two classes: (1) to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic behaviour. We show that a simple two-component model can very well reproduce the observed phenomenologies. (b) We present the cm-to-mm behaviour of three gamma-ray bright Narrow Line Seyfert 1 galaxies over time spans varying between ~1.5 and 3 years and compare their variability characteristics with typical blazars.

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

After a long period of quiescence in \gamma-rays, blazar 0836+710 (4C +71.07) flared in the Spring of 2011. We found only limited multiwavelength coverage of the source. An indication of correlated optical / \gamma-ray variability is not surprising for a Flat Spectrum Radio Quasar (FSRQ) like this one. Radio observations at high frequencies, however, had seen a flare in 2010, well offset from possible \gamma-ray activity. The 2011 \gamma-ray activity comes during a period of rising radio emission, a pattern that has been seen since the EGRET era.

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

Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

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

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

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

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of {\em Fermi} blazars. The current study is concerned with the broad-band radio spectra composed of measurement at ten frequencies between 2.64 and 142 GHz. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. The first four types are dominated by spectral evolution and can be reproduced by a simple two-component system made of the quiescent spectrum of a large scale jet populated with a flaring event evolving according to Marscher & Gear (1985). The last type is characterized by an achromatic change of the broad-band spectrum which must be attributed to a completely different mechanism. Here are presented, the classification, the assumed physical system and the results of simulations that have been conducted.

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

Using the Urumqi 25m radio telescope, sources from the first three-month Fermi-LAT detected AGN catalog with declination >0 were observed in 2009 at 4.8 GHz. The radio flux density appears to correlate with the gamma-ray intensity. Intra-day variability (IDV) observations were performed in March, April and May in 2009 for selected 42 gamma-ray bright blazars, and 60% of them show evident flux variability at 4.8 GHz during the IDV observations, the IDV detection rate is higher than that in previous flat-spectrum AGN samples. The IDV appears more often in the VLBI-core dominant blazars, and the non-IDV blazars show relatively `steeper' spectral indices than the IDV blazars. Pronounced inter-month variability has been found in two BL Lac objects: J0112+2244 and J0238+1636.

Simultaneous multi-wavelength (MW) campaigns are the most promising approaches to investigate the still unrevealed nature of blazars, active galactic nuclei which are variable on all time scales from radio to TeV energies. In 2008, two MW campaigns on the high-frequency peaked blazars Mrk 180 and 1ES 2344+514 have been organised by the MAGIC collaboration. From radio to TeV gamma rays, RATAN-600, Mets\"ahovi, Effelsberg, VLBA (only 1ES 2344+514), IRAM, KVA, Swift, AGILE, Fermi-LAT and MAGIC-I were taking part in these campaigns. Mrk 180 had just been discovered at TeV energies by MAGIC in 2006, whereas 1ES 2344+514 is a known TeV emitter since many years. Due to their rather faint emission particularly at TeV energies, the campaigns represented quite challenging observations. In fact, Mrk 180 has not been investigated until now in MW campaigns, and for 1ES 2344+514 only one campaign including TeV measurements has been reported in literature up to now. In this contribution, we will present detailed MW light curves for both sources and describe the composite wide range spectral energy distributions by theoretical models.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

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

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

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

The Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope provides an unprecedented opportunity to study gamma-ray blazars. To capitalize on this opportunity, beginning in late 2007, about a year before the start of LAT science operations, we began a large-scale, fast-cadence 15 GHz radio monitoring program with the 40-m telescope at the Owens Valley Radio Observatory (OVRO). This program began with the 1158 northern (declination>-20 deg) sources from the Candidate Gamma-ray Blazar Survey (CGRaBS) and now encompasses over 1500 sources, each observed twice per week with a ~4 mJy (minimum) and 3% (typical) uncertainty. Here, we describe this monitoring program and our methods, and present radio light curves from the first two years (2008 and 2009). As a first application, we combine these data with a novel measure of light curve variability amplitude, the intrinsic modulation index, through a likelihood analysis to examine the variability properties of subpopulations of our sample. We demonstrate that, with high significance (7-sigma), gamma-ray-loud blazars detected by the LAT during its first 11 months of operation vary with about a factor of two greater amplitude than do the gamma-ray-quiet blazars in our sample. We also find a significant (3-sigma) difference between variability amplitude in BL Lacertae objects and flat-spectrum radio quasars (FSRQs), with the former exhibiting larger variability amplitudes. Finally, low-redshift (z<1) FSRQs are found to vary more strongly than high-redshift FSRQs, with 3-sigma significance. These findings represent an important step toward understanding why some blazars emit gamma-rays while others, with apparently similar properties, remain silent.

Narrow-Line Seyfert 1 (NLS1) class of active galactic nuclei (AGNs) is generally radio-quiet, but a small percent of them are radio-loud. The recent discovery by Fermi/LAT of high-energy gamma-ray emission from 4 NLS1s proved the existence of relativistic jets in these systems. It is therefore important to study this new class of gamma-ray emitting AGNs. Here we report preliminary results about the observations of the July 2010 gamma-ray outburst of PMN J0948+0022, when the source flux exceeded for the first time 10^-6 ph cm^-2 s^-1 (E > 100 MeV).

Since 2008 the Fermi/LAT instrument has delivered highly time-resolved gamma-ray spectra and detailed variability curves for a steadily increasing number of AGN. For detailed AGN studies the Fermi/LAT data have to be combined with, and accompanied by, dedicated ground- and space-based multi-frequency observations. In this framework, the Fermi AGN team has realized a detailed plan for multi-wavelength campaigns including a large suite of cm/mm/sub-mm band instruments. Many of those campaigns have been triggered, often for sources detected in flaring states. We review here a few interesting results recently obtained during three such campaigns, namely for the flat-spectrum radio quasar 3C 279, the Narrow Line Seyfert 1 PMN J0948+0022 and quasar 3C 454.3.

The F -GAMMA program is a coordinated effort of several observing facilities to understand the AGN/blazar phenomenon via a multi-frequency monitoring approach, especially in the era of Fermi-GST. Some 60 prominent sources are monitored monthly with the Effelsberg 100-m telescope, the IRAM 30-m telescope and more frequently but in a less uniform fashion at the APEX 12-m telescope, covering from 2.64 to 345 GHz. The program has been running since January 2007 and here some of its findings are summarized. (a) There are two major variability patterns that the spectra of sources follow, one spectral-evolution-dominated and one achromatic. (b) The FSRQs show higher brightness temperatures indicative of larger Doppler factors at play and (c) a statistically significant radio-gamma-ray correlation has been found with a method recently suggested by Pavlidou et al. (in prep.).

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

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

The \emph{Fermi} Large Area Telescope (LAT) discovered a new gamma-ray source near the Galactic plane, \object{Fermi J0109+6134}, when it flared brightly in 2010 February. The low Galactic latitude (b =-1.2\degr) indicated that the source could be located within the Galaxy, which motivated rapid multi-wavelength follow-up including radio, optical, and X-ray observations. We report the results of analyzing all 19 months of LAT data for the source, and of X-ray observations with both \emph{Swift} and the \emph{Chandra X-ray Observatory}. We determined the source redshift, z =0.783, using a Keck LRIS observation. Finally, we compiled a broadband spectral energy distribution (SED) from both historical and new observations contemporaneous with the 2010 February flare. The redshift, SED, optical line width, X-ray absorption, and multi-band variability indicate that this new GeV source is a blazar seen through the Galactic plane. Because several of the optical emission lines have equivalent width >5\AA, this blazar belongs in the flat-spectrum radio quasar category.

Radio and gamma-ray emission from active galactic nuclei (AGN) are thought to share a common origin, related to the ejection phenomena in the vicinity of supermassive black holes. Thanks to its sensitivity, surveying capability, and broad energy range, the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope has permitted us to discover and characterize a huge number of extragalactic gamma-ray sources. Similarly to what was found by EGRET, these sources are typically associated with blazars, characterized by significant radio emission and flat spectrum. The radio luminosity distribution is extended over 7 orders of magnitudes, with flat spectrum radio quasars clustered at higher powers and BL Lacs more scattered; the average spectral index is consistent with alpha=0, although a few remarkable sources have alpha>0.5. A comparison of the radio flux density and the gamma-ray photon flux is presented, although claims on its significance require a detailed discussion and Monte Carlo simulations which will be presented in a future paper.

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

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

We present results from the first two years of our fast-cadence 15 GHz gamma-ray blazar monitoring program, part of the F-GAMMA radio monitoring project. Our sample includes the 1158 blazars north of -20 degrees declination from the Candidate Gamma-Ray Blazar Survey (CGRaBS), which encompasses a significant fraction of the extragalactic sources detected by the Fermi Gamma-ray Space Telescope. We introduce a novel likelihood analysis for computing a time series variability amplitude statistic that separates intrinsic variability from measurement noise and produces a quantitative error estimate. We use this method to characterize our radio light curves. We also present results indicating a statistically significant correlation between simultaneous average 15 GHz radio flux density and gamma-ray photon flux.

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

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

Short time-scale flux density variations of flat spectrum radio sources are often explained by the scattering of radio waves in the turbulent, ionized Interstellar Matter of the Milky Way. One of the most convincing observational arguments in favor of this is the annual modulation of the variability time-scale caused by the Earth orbital motion around the Sun. J1128+592 is an IDV source with a possible annual modulation in its variability time-scale. We observed the source in 6 epochs with the VLBA at 5, 8 and 15 GHz in total intensity and polarization. The VLBA observations revealed an east-west oriented core-jet structure. Its position angle agrees with the angle of anisotropy derived from the annual modulation model. No significant long-term structural changes were observed with VLBI on mas-scales, however, compared to archival data, the VLBI core size is expanded. This expansion offers a possible explanation to the observed decrease of the strength of IDV. VLBI polarimetry revealed significant changes in the electric vector position angle and Rotation Measure of the core and jet. Part of the observed RM variability could be attributed to a scattering screen (37 pc distance), which covers the source (core and jet) and which may be responsible for the IDV. Superposition of polarized sub-components below the angular resolution limit may affect the observed RM as well.

The F-GAMMA-project is the coordinated effort of several observatories to understand the AGN phenomenon and specifically blazars via multi-frequency monitoring in collaboration with the {\sl Fermi}-GST satellite since January 2007. The core observatories are: the Effelsberg 100-m, the IRAM 30-m and the OVRO 40-m telescope covering the range between 2.6 and 270 GHz. Effelsberg and IRAM stations do a monthly monitoring of the cm to sub-mm radio spectra of 60 selected blazars whereas the OVRO telescope is observing roughly 1200 objects at 15 GHz with a dense sampling of 2 points per week. The calibration uncertainty even at high frequencies, is of a few percent. 47% of the Effelsberg/Pico Veleta sample is included in the LBAS list. An update of the monitored sample is currently underway.

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

We describe a 4.85 GHz survey of bright, flat-spectrum radio sources conducted with the Effelsberg 100 m telescope in an attempt to improve the completeness of existing surveys, such as CRATES. We report the results of these observations and of follow-up 8.4 GHz observations with the VLA of a subset of the sample. We comment on the connection to the WMAP point source catalog and on the survey's effectiveness at supplementing the CRATES sky coverage.

The analysis of the spectral energy distribution variability at frequencies from radio to TeV is a powerful tool in the investigation of the dynamics, the physics and the structure evolution occurring in the most exotic flavour of active galaxies, the blazars. In particular, the presence of {\sl Fermi-GST} is providing a unique opportunity for such studies delivering $\gamma$-ray data of unprecedented quality. Here we introduce a monitoring program that has been running at the Effelsberg 100m telescope since January 2007, underpinning a broad multi-frequency collaboration of facilities that cover the band from radio to infrared. Sixty one selected blazars are observed monthly between 2.64 GHz and 43 GHz. The calibration accuracy is better than a few percent as it is demonstrated with some preliminary examples.

The BL Lac object S5 0716+71 was observed in a global multi-frequency campaign to search for rapid and correlated flux density variability and signatures of an inverse-Compton (IC) catastrophe during the states of extreme apparent brightness temperatures. The observing campaign involved simultaneous monitoring at radio to IR/optical wavelengths centered around a 500-ks INTEGRAL pointing (November 10-17, 2003). We present the combined analysis and results of the cm- to sub-mm observations including a detailed study of the inter- to intra-day variability and spectral characteristics of 0716+714. We further constrain the variability brightness temperatures (T_B) and Doppler factors (delta) comparing the radio-bands with the hard X-ray emission (3-200 keV). 0716+714 was in an exceptionally high state (outburst) and different (slower) phase of short-term variability. The flux density variability in the cm- to mm-bands is dominated by a correlated, ~4 day time scale amplitude increase of up to ~35% systematically more pronounced towards shorter wavelengths. This contradicts expectations from standard interstellar scintillation (ISS) and suggests a source-intrinsic origin of the variability. The derived lower limits to T_B exceed the 10^12 K IC-limit by up to 3-4 orders of magnitude. Assuming relativistic boosting, we obtain robust and self-consistent lower limits of delta >= 5-33, in good agreement with delta_VLBI obtained from VLBI studies and the IC-Doppler factors delta_IC > 14-16 obtained from the INTEGRAL data. Since a strong contribution from ISS can be excluded and a simultaneous IC catastrophe was not observed, we conclude that relativistic Doppler boosting naturally explains the apparent violation of the theoretical limits within standard synchrotron-self-Compton (SSC) jet models of AGN.

A sample of classical IntraDay Variable (IDV) and IDV candidate sources has been monitored with the Urumqi 25m telescope and the Effelsberg 100m telescope. Aim of the project is to investigate the origin of IntraDay Variability, a phenomenon which has been observed in about 30% of flat spectrum radio quasars. Simultaneous Effelsberg-Urumqi observations demonstrated that the Urumqi antenna, although relatively small in diameter, is well suitable for IDV experiments. A few Urumqi datasets, however, turned out to be affected by a spurious $\sim 24$ hours modulation, an effect which has been removed by means of a new procedure for data reduction. In about 14 months, 12 epochs of observation have been collected, for a total observing time of more than 45 days. The epochs are regularly distributed over the whole year, in order to check for the presence of systematic annual changes in the variability time scales - a crucial test for verifying the consistency of source-extrinsic models of the variability. Preliminary time-analysis of the monitored sources revealed some hint for a slowing down of the characteristic time scales of S5~0716+714, a result that, if confirmed, would be compatible with a source-extrinsic origin of the variability, in contrast to previous IDV studies. No significant modulation of the time scales has been detected for S4~0954+658.

To locate and image the compact emission regions in quasars, which are closely connected to the phenomenon of IntraDay Variability (IDV), space VLBI observations are of prime importance. Here we report on VSOP observations of two prominent IDV sources, the BL Lac objects S5 0716+714. To monitor their short term variability, these sources were observed with VSOP at 5 GHz in several polarisation sensitive experiments, separated in time by one day to six days, in autumn 2000. Contemporaneous flux density measurements with the Effelsberg 100m radio telescope were used to directly compare the single dish IDV with changes of the VLBI images. A clear IDV behaviour in total intensity and linear polarization was observed in 0716+714. Analysis of the VLBI data shows that the variations are located inside the VLBI core component of 0716+714. In good agreement with the single-dish measurements, the VLBI ground array images and the VSOP images, both show a decrease in the total flux density of ~20 mJy and a drop of ~5 mJy in the linear polarization of the VLBI core. No variability was found in the jet. From the variability timescales we estimate a source size of a few micro-arcseconds and brightness temperatures exceeding 10^15 K. Independent of whether the interpretation of the IDV seen in the VLBI core is source intrinsic or extrinsic a lower limit of T_B > 2x10^12 K is obtained by model fitting of the VLBI-core. Our results show that future VSOP2 observations should be accompanied by a single dish monitoring not only to discriminate between source-extrinsic and source-intrinsic effects but to allow also a proper calibration and interpretation of ultra-high resolution VSOP2 images.

For more than a decade from the late 1980s to the late 1990s, the intraday variable (IDV) quasar 0917+624 was known to be strongly variable with amplitude variations of 10-20% within 0.8 to 1.6 days. During that period it showed faster variability in the polarized flux and the polarization angle than in total intensity. However, since 2000 the IDV of this source has almost ceased, attaining a variability index of only less than a few percent. If interstellar scintillation was responsible for the IDV in the past, one possible explanation for the sudden change of the variability mode could be an intrinsic change in the morphology of the source, such as an increase in the size of the structural components that are responsible for the scintillation. New jet components appearing at the jet base could lead to a temporarily larger size of the scintillating component when compared to the scattering size of the interstellar medium, causing quenched scintillation. In order to test this scenario, VLBI observations of 0917+624 from 1999 to 2005 have been analysed. Here, we summarize first results and discuss them in the framework of the quenched scintillation model.

Short time-scale radio variations of compact extragalactic radio quasars and blazars known as IntraDay Variability (IDV) can be explained in at least some sources as a propagation effect; the variations are interpreted as scintillation of radio waves in the turbulent interstellar medium of the Milky Way. One of the most convincing observational arguments in favor of a propagation-induced variability scenario is the observed annual modulation in the characteristic time scale of the variation due to the Earth's orbital motion. So far there are only two sources known with a well-constrained seasonal cycle. Annual modulation has been proposed for a few other less well-documented objects. However, for some other IDV sources source-intrinsic structural variations which cause drastic changes in the variability time scale were also suggested. J1128+592 is a recently discovered, highly variable IDV source. Previous, densely time-sampled flux-density measurements with the Effelsberg 100-m radio telescope (Germany) and the Urumqi 25-m radio telescope (China), strongly indicate an annual modulation of the time scale. The most recent 4 observations in 2006/7, however, do not fit well to the annual modulation model proposed before. In this paper, we investigate a possible explanation of this discrepancy.

The Large Area Telescope (LAT) instrument onboard GLAST offers a tremendous opportunity for future blazar studies. In order to fully benefit from its capabilities and to maximize the scientific return from the LAT, it is of great importance to conduct dedicated multi-frequency monitoring campaigns that will result comprehensive observations. Consequently, we initiated an effort to conduct a GLAST-dedicated, quasi-simultaneous, broad-band flux-density (and polarization) monitoring of potential GLAST blazars with the Effelsberg and OVRO radio telescopes (11cm to 7mm wavelength). Here, we present a short overview of these activities which will complement the multi-wavelengths activities of the GLAST/LAT collaboration towards the 'low-energy' radio bands. Further we will give a brief outlook including the extension of this coordinated campaign towards higher frequencies and future scientific aims.

Short time-scale radio variations of compact extragalactic radio sources, known as IntraDay Variability, can be explained in at least some sources by a source-extrinsic effect, in which the variations are interpreted as scintillation of radio waves caused by the turbulent ISM of the Milky Way. One of the most convincing observational arguments in favour of propagation-induced variability is the so called annual modulation of the characteristic variability time-scale, which is due to the orbital motion of the Earth. Data for the recently discovered and highly variable IDV source J1128+5925 are presented. We study the frequency and time dependence of the IDV in this compact quasar. We measure the characteristic variability time-scale of the IDV throughout the year, and analyze whether the observed changes in the variability time-scale are consistent with annual modulation. We monitored the flux density variability of J1128+5925 with dense time sampling between 2.7 and 10.45GHz with the 100m Effelsberg radio telescope of the MPIfR and with the 25m Urumqi radio telescope. From ten observing sessions, we determine the variability characteristics and time-scales. The observed pronounced changes of the variability time-scale of J1128+5925 are modelled with an anisotropic annual modulation model. The observed frequency dependence of the variation is in good agreement with the prediction from interstellar scintillation. Adopting a simple model for the annual modulation model and using also the frequency dependence of the IDV, we derive a lower limit to the distance of the scattering screen and an upper limit to the scintillating source size. The latter is found to be consistent with the measured core size from VLBI.

In radio bands, the study of compact radio sources can be affected by propagation effects through the Interstellar Medium. These are usually attributed to the presence of turbulent intervening plasma along the line of sight. Here, two of such effects are presented. The line of sight of B 2005+403 passes through the heavily scattered region of Cygnus, which causes substantial angular broadening of the source images obtained at frequencies between 0.6 GHz and 8 GHz. At higher frequencies, however the intrinsic source structure shines through. Therefore, multi-frequency VLBI observations allow to study the characteristics of the intervening material, the source morphology and the "coupling" of them in forming the observed image.

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.

The BL Lac object S5 0716+71 was target of a coordinated and global multi-frequency campaign to search for rapid and correlated variability and signatures of the inverse-Compton catastrophe. Here we present first results obtained from a combined analysis of the cm- to sub-mm observations over a period of seven days aiming at a detailed study of the intra- to inter-day variability characteristics and to obtain constraints on the variability brightness temperatures and Doppler factors comparing the radio data with the high energy emission recorded by INTEGRAL. A more detailed description of the whole cm- to sub-mm observations and our analysis/results will be presented in a forthcoming paper. Our analysis reveals the source to be in a particular short-term variability phase when compared to the past with a correlated >~4 day time scale amplitude increase of up to 35%, which is systematically more pronounced towards higher frequencies. The obtained frequency dependent variability amplitudes and time lags contradict expectations from interstellar scintillation and strongly suggest a source intrinsic origin of this inter-day variability. A 7-day spectral evolution study indicate time-variable synchrotron self-absorption and expansion of the emission region, consistent with standard models. Assuming relativistic boosting, our different estimates of the Doppler factor yield robust lower limits of D_var,IC > 5-22 using the inverse-Compton limit and D_var,eq > 8-33 using the equipartition argument. Although high, these values are in good agreement with Doppler factors obtained from recent VLBI studies and from the inverse-Compton Doppler factors D_IC > 14-16 derived with the X-ray emission seen by INTEGRAL at 3-200 keV.

We present the results of a series of Swift and quasi simultaneous ground-based infra-red observations of the blazar 3C 454.3 carried out in April-May 2005 when the source was 10 to 30 times brighter than previously observed. We found 3C 454.3 to be very bright and variable at all frequencies covered by our instrumentation. The broad-band Spectral Energy Distribution (SED) shows the usual two-bump shape (in nu-nu f(nu) space) with the Infra-red, optical and UV data sampling the declining part of the synchrotron emission that, even during this extremely large outburst, had its maximum in the far-infrared. The X-ray spectral data from the XRT and BAT instruments are flat and due to inverse Compton emission. The remarkable SED observed implies that at the time of the Swift pointings 3C 454.3 was one of the brightest objects in the extragalactic sky with a gamma-ray emission similar or brighter than that of 3C 279 when observed in a high state by EGRET. Time variability in the optical-UV flux is very different from that in the X-ray data: while the first component varied by about a factor two within a single exposure, but remained approximately constant between different observations, the inverse Compton component did not vary on short time-scales but changed by more than a factor of 3 between observations separated by a few days. This different dynamical behaviour illustrates the need to collect simultaneous multi-frequency data over a wide range of time-scales to fully constrain physical parameters in blazars.

We report on a densely time sampled polarimetric flux density monitoring of the BL Lac object S5 0716+71 at 86 GHz and 229 GHz with the IRAM 30 m telescope within a coordinated broad spectral band campaign, centred around an INTEGRAL observation during November 10 to 16, 2003. Our aim was to search for signatures of inverse-Compton "catastrophes". At 86 GHz, making use of a new calibration strategy, we reach a relative rms accuracy of the flux density measurements of 1.2%. At this frequency, S5 0716+71 showed no intra-day variability, but showed remarkable inter-day variability with a flux density increase of 34% during the first four observing days, which can not be explained by source extrinsic causes. The 86 GHz linear polarization fraction of S5 0716+71 was unusually large 15.0+-1.8%. Inter-day variability in linear polarization at 86 GHz, with significance level >~95%; sigma_P/<P>=15% and sigma_chi=6 deg., was also observed. From the emission variations at the synchrotron turnover frequency (~86 GHz) we compute an apparent brightness temperature T_B,app>1.4x10^14K at a redshift of 0.3, which exceeds by two orders of magnitude the inverse-Compton limit. A relativistic correction for T_B,app with a Doppler factor delta > 7.8 brings the observed brightness temperature down to the inverse Compton limit. A more accurate lower limit of delta > 14.0, is obtained from the comparison of the 86 GHz synchrotron flux density and the upper limits for the synchrotron self-Compton flux density obtained from the INTEGRAL observations. The relativistic beaming of the emission by this high Doppler factor explains the non-detection of "catastrophic" inverse-Compton avalanches by INTEGRAL.

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.

The flat-spectrum radio quasar 3C 454.3 is well known to be a highly active and variable source with outbursts occurring across the whole electromagnetic spectrum over the last decades. In spring 2005, 3C 454.3 has been reported to exhibit a strong optical outburst which subsequently triggered multi-frequency observations of the source covering the radio up to gamma-ray bands. Here, we present first results of our near-IR/optical (V, R, I, H band) photometry performed between May 11 and August 5, 2005 with the Rapid Eye Mount (REM) at La Silla in Chile and the Automatic Imaging Telescope (AIT) of the Perugia University Observatory. 3C 454.3 was observed during an exceptional and historical high state with a subsequent decrease in brightness over our 86 days observing period. The continuum spectral behaviour during the flaring and declining phase suggests a synchrotron peak below the near-IR band as well as a geometrical origin of the variations e.g. due to changes in the direction of forward beaming.

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.

The evidence for refractive interstellar scintillation (RISS) being the main cause for rapid intraday variations (Intraday Variability, IDV) in Quasars and BL Lacs has recently become stronger. If IDV is still a complex composition of extrinsic and source intrinsic effects, the intrinsic part of the IDV pattern should show up in the millimeter and sub-millimeter regime due to the frequency dependence of RISS. Hence, observations at higher frequencies are essential in order to exclude RISS as the sole cause of IDV. Here we report on our new attempt to search for rapid variations at much higher frequencies. In addition, the possibility of a direct detection of the postulated scattering screen in front of IDV sources will be discussed. Our recent line observations towards a few IDV sources lead to the first detection of a high latitude molecular cloud in front of an intraday variable radio core.

Intraday Variability of compact extragalactic radio sources can be interpreted as quenched scintillation due to turbulent density fluctuations of the nearby ionized interstellar medium. We demonstrate that the statistical analysis of IDV time series contains both information about sub-structure of the source on the scale of several 10 micro-arcsec and about the turbulent state of the ISM. The source structure and ISM properties cannot be disentangled using IDV observations alone. A comparison with the known morphology of the `local bubble' and the turbulent ISM from pulsar observations constrains possible source models. We further argue that earth orbit synthesis fails for non-stationary relativistic sources and no reliable 2D-Fourier reconstruction is possible.

Intraday variability (IDV) of compact extragalactic radio sources is a complex phenomenon and shows a wavelength dependent mixture of refractive interstellar scintillation (RISS) (dominant at long cm-wavelengths) and source-intrinsic effects (dominant at shorter wavelengths). Detailed investigations of individual sources and new high frequency observations suggest a source-intrinsic contribution to the IDV pattern at least in some sources. However, the sizes of intraday variable sources at cm-wavelength are typically smaller than the scattering size set by the ISM in our galaxy and scintillation must be present, too. We present new IDV observations in different regimes (from cm to sub-mm wavelengths) and show how such multi-frequency study can be used as powerful instrument to describe different aspects of IDV.

Radio variability on timescales from a few hours to several days in extragalactic flat-spectrum radio sources is generally classified as intra-day variability (IDV). The origin of this short term variability is still controversial and both extrinsic and intrinsic mechanisms must be considered and may both contribute to the observed variations. The measured linear and circular polarization of IDV sources constrains the low energy end of the electron population. Any population of cold electrons within sources at or above the equipartition temperature of 10^11 K depolarizes the emission and can be ruled out. Intrinsic shock models are shown to either violate the large fraction of sources displaying IDV or they do not relax the light travel time argument for intrinsic variations. From structure function analysis, we further conclude that interstellar scintillation also leads to tight size estimates unless a very local cloud in the ISM is responsible for IDV.

There is new evidence which identifies seasonal changes of the variability time scale in intraday variable sources with refractive interstellar scintillation effects. Such a RISS-model takes the annual change of the Earth velocity relative to the scattering medium into account. In September 1998 we found a remarkable prolongation of the variability time scale in the IDV-source 0917+624 with only small variations in flux density during a period of 5 days. This was explained as a seasonal effect, in which the velocity vector of the Earth and the interstellar medium nearly canceled. In order to further investigate the applicability of the model for 0917+624, we performed an Effelsberg 6cm-flux monitoring program over the course of one year. Since September 2000, the source appears to be remarkably inactive and yet (May 2001), no return to its normal, faster and stronger variability pattern is observed. Here, our observational results and a possible explanation for the current quiescence are presented.

We summarize results from flux density monitoring campaigns performed with the 100 meter radio-telescope at Effelsberg and the VLA during the past 15 yrs. We briefly discuss some of the statistical properties of the rapid variability from now more than 40 high declination sources, which show Intraday Variability (IDV). In general, IDV is more pronounced for sources with flat radio spectra and compact VLBI structures. The frequency dependence of the variability pattern varies with source and observing time. For 0917+62, we present new VLBI images, which suggest that the variability pattern is modified by the occurrence of new jet components. For 0716+71, we show the first detection of IDV at millimeter wavelengths (32 GHz). For the physical interpretation of the IDV phenomenon, a complex source and frequency dependent superposition of interstellar scintillation and source intrinsic variability should be considered.

Extreme Scattering Events (ESEs) are dramatic variations of the flux density at Gigahertz frequencies caused by ray path distortions within an isolated inhomogeneity ("plasma lens") in the interstellar medium. These events are characterized by a deep flux density minimum in the light curve with, in some cases, surrounding maxima. The variability time scales range from weeks to months. These phenomena show a strong frequency dependence, in which the variability amplitudes increase with wavelength. During an Intraday Variability (IDV) monitoring project (March 2000), a feature resembling an ESE-like event appeared in the variable light curve of 0954+658, however with a time scale of less than two days. We will discuss this effect and its implications for a better description of the interstellar medium.