A number of works reported on the existence of a large scale alignment of the polarization plane of extragalactic sources as well as the alignment of radio-sources structural axes. However, both claims and their interpretation remain controversial. For the first time we explore the parsec-scale jets alignments. Additionally, we use archival polarimetric data at different wavelengths in order to compare relative orientations of the jets and the polarization planes of their emission. Using the flux density distribution in very long baseline interferometry (VLBI) radio maps from the Astrogeo database, we determine the parsec-scale jet orientation for the largest sample of active galactic nuclei (AGN) to date. Employing the method of parallel transport and a sample statistics characterizing the jet orientation dispersion among neighbors, we test whether the identified jets are significantly aligned. We show that the parsec-scale jets in our sample do not demonstrate any significant global alignments. Moreover, the jet direction is found to be weakly correlated with the polarization plane direction at different frequencies.

There are several methods to calculate the radiative and kinetic power of relativistic jets, but their results can differ by one or two orders of magnitude. Therefore, it is necessary to perform a calibration of the jet power, to understand the reasons for these differences (whether wrong hypotheses or intrinsic source variability), and if it is possible to converge to a reliable measurement of this physical quantity. We present preliminary results of a project aimed at calibrating the power of relativistic jets in active galactic nuclei (AGN) and X-ray binaries (XRB). We started by selecting all the AGN associations with known redshift in the Fourth Fermi LAT Gamma-Ray Catalog (4FGL). We then calculated the radiative and/or kinetic powers from available data or we extracted this information from literature. We compare the values obtained for overlapping samples and highlight early conclusions.

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.

We present a multi-frequency study of PKS J1222$+$0413 (4C$+$04.42), currently the highest redshift $\gamma$-ray emitting narrow-line Seyfert 1 ($\gamma$-NLS1). We assemble a broad spectral energy distribution (SED) including previously unpublished datasets: X-ray data obtained with the NuSTAR and Neil Gehrels Swift observatories; near-infrared, optical and UV spectroscopy obtained with VLT X-shooter; and multiband radio data from the Effelsberg telescope. These new observations are supplemented by archival data from the literature. We apply physical models to the broadband SED, parameterising the accretion flow and jet emission to investigate the disc-jet connection. PKS J1222$+$0413 has a much greater black hole mass than most other NLS1s, $M_\mathrm{BH}\approx2\times10^{8}$ M$_\odot$, similar to those found in flat spectrum radio quasars (FSRQs). Therefore this source provides insight into how the jets of $\gamma$-NLS1s relate to those of FSRQs.

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

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 third Fermi-LAT catalog (3FGL) presented the data of the first four years of observations from the Fermi Gamma-ray Space Telescope mission. There are 3034 sources, 1010 of which still remain unidentified. Identifying and classifying gamma-ray emitters is of high significance with regard to studying high-energy astrophysics. We demonstrate that optical polarimetry can be an advantageous and practical tool in the hunt for counterparts of the unidentified gamma-ray sources (UGSs). Using data from the RoboPol project, we validated that a significant fraction of active galactic nuclei (AGN) associated with 3FGL sources can be identified due to their high optical polarization exceeding that of the field stars. We performed an optical polarimetric survey within $3\sigma$ uncertainties of four unidentified 3FGL sources. We discovered a previously unknown extragalactic object within the positional uncertainty of 3FGL J0221.2+2518. We obtained its spectrum and measured a redshift of $z=0.0609\pm0.0004$. Using these measurements and archival data we demonstrate that this source is a candidate counterpart for 3FGL J0221.2+2518 and most probably is a composite object: a star-forming galaxy accompanied by AGN. We conclude that polarimetry can be a powerful asset in the search for AGN candidate counterparts for unidentified Fermi sources. Future extensive polarimetric surveys at high galactic latitudes (e.g., PASIPHAE) will allow the association of a significant fraction of currently unidentified gamma-ray sources.

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.

We present a multi-frequency, dense radio monitoring program of the blazar OJ287 using the 100m Effelsberg radio telescope. We analyze the evolution in total flux density, linear and circular polarization to study the jet structure and its magnetic field geometry. The total flux density is measured at nine bands from 2.64 GHz to 43 GHz, the linear polarization parameters between 2.64 GHz and 10.45 GHz, and the circular polarization at 4.85 GHz and 8.35 GHz. The mean cadence is 10 days. Between MJD 57370 and 57785, OJ287 showed flaring activity and complex linear and circular polarization behavior. The radio EVPA showed a large clockwise (CW) rotation by ~340$^{\circ}$ with a mean rate of -1.04 $^{\circ}$/day. Based on concurrent VLBI data, the rotation seems to originate within the jet core at 43 GHz (projected size $\le$ 0.15 mas or 0.67 pc). Moreover, optical data show a similar monotonic CW EVPA rotation with a rate of about -1.1 $^{\circ}$/day which is superposed with shorter and faster rotations of about 7.8 $^{\circ}$/day. The observed variability is consistent with a polarized emission component propagating on a helical trajectory within a bent jet. We constrained the helix arc length to 0.26 pc and radius to $\le$ 0.04 pc as well as the jet bending arc length projected on the plane of the sky to $\le$ 1.9-7.6 pc. A similar bending is observed in high angular resolution VLBI images at the innermost jet regions. Our results indicate also the presence of a stable polarized emission component with EVPA (-10$^{\circ}$) perpendicular to the large scale jet, suggesting dominance of the poloidal magnetic field component. Finally, the EVPA rotation begins simultaneously with an optical flare and hence the two might be physically connected. That optical flare has been linked to the interaction of a secondary SMBH with the inner accretion disk or originating in the jet of the primary.

We have monitored the R-band optical linear polarisation of ten jetted NLSy1 galaxies with the aim to quantify their variability and search for candidate long rotation of the polarisation plane. In all cases for which adequate datasets are available we observe significant variability of both the polarisation fraction and angle. In the best-sampled cases we identify candidate long rotations of the polarisation plane. We present an approach that assesses the probability that the observed phenomenology is the result of pure noise. We conclude that although this possibility cannot be excluded it is much more likely that the EVPA undergoes an intrinsic evolution. We compute the most probable parameters of the intrinsic event which forecasts events consistent with the observations. In one case we find that the EVPA shows a preferred direction which, however, does not imply any dominance of a toroidal or poloidal component of the magnetic field at those scales.

Narrow line Seyfert 1 galaxies (NLSy1s) constitute the AGN subclass associated with systematically smaller black hole masses. A few radio loud ones have been detected in MeV -- GeV energy bands by Fermi and evidence for the presence of blazar-like jets has been accumulated. In this study we wish to quantify the temporal behaviour of the optical polarisation, fraction and angle, for a selected sample of radio loud NLSy1s. We also search for rotations of the polarisation plane similar to those commonly observed in blazars. We have conducted R-band optical polarisation monitoring of a sample of 10 RL NLSy1s 5 of which have been previously detected by Fermi. The dataset includes observations with the RoboPol, KANATA, Perkins and Steward polarimeters. In the cases where evidences for long rotations of the polarisation plane are found, we carry out numerical simulations to assess the probability that they are caused by intrinsically evolving EVPAs instead of observational noise. Even our moderately sampled sources show indications of variability, both in polarisation fraction and angle. For the four best sampled objects in our sample we find multiple periods of significant polarisation angle variability. In the two best sampled cases, namely J1505+0326 and J0324+3410, we find indications for three long rotations. We show that although noise can induce the observed behaviour, it is much more likely that the apparent rotation is caused by intrinsic evolution of the EVPA. To our knowledge this is the very first detection of such events in this class of sources. In the case of the largest dataset (J0324+3410) we find that the EVPA concentrates around a direction which is at 49.3\degr to the 15-GHz radio jet implying a projected magnetic field at an angle of 40.7\degr to that axis.

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

We present the first results of an ongoing intra-day variability (IDV) flux density monitoring program of 107 blazars, which were selected from a sample of RadioAstron space very long baseline interferometry (VLBI) targets. The~IDV observations were performed with the Effelsberg 100-m radio telescope at 4.8\,GHz, focusing on the statistical properties of IDV in a relatively large sample of compact active galactic nuclei (AGN). We investigated the dependence of rapid ($<$3 day) variability on various source properties through a likelihood approach. We found that the IDV amplitude depends on flux density and that fainter sources vary by about a factor of 3 more than their brighter counterparts. We also found a significant difference in the variability amplitude between inverted- and flat-spectrum radio sources, with the former exhibiting stronger variations. $\gamma$-ray loud sources were found to vary by up to a factor 4 more than $\gamma$-ray quiet ones, with 4$\sigma$ significance. However a galactic latitude dependence was barely observed, which suggests that it is predominantly the intrinsic properties (e.g., angular size, core-dominance) of the blazars that determine how they scintillate, rather than the directional dependence in the interstellar medium (ISM). We showed that the uncertainty in the VLBI brightness temperatures obtained from the space VLBI data of the RadioAstron satellite can be as high as $\sim$70\% due to the presence of the rapid flux density variations. Our statistical results support the view that IDV at centimeter wavelengths is predominantly caused by interstellar scintillation (ISS) of the emission from the most compact, core-dominant region in an AGN.

The RoboPol instrument and the relevant program was developed in order to conduct a systematic study of the optical polarisation variability of blazars. Driven by the discovery that long smooth rotations of the optical polarisation plane can be associated with the activity in other bands and especially in gamma rays, the program was meant to investigate the physical mechanisms causing them and quantify the optical polarisation behaviour in blazars. Over the first three nominal observing seasons (2013, 2014 and 2015) RoboPol detected 40 rotations in 24 blazars by observing a gamma-ray-loud and gamma-ray-quite unbiassed sample of blazars, providing a reliable set of events for exploring the phenomenon. The obtain datasets provided the ground for a systematic quantification of the variability of the optical polarisation in such systems. In the following after a brief review of the discoveries that relate to the gamma-ray loudness of the sources we move on to discuss a simple jet model that explains the observed dichotomy in terms of polarisation between gamma-ray-loud and quite sources and the dependence of polarisation and the stability of the polarisation angle on the synchrotron peak frequency.

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.

We use results of our 3 year polarimetric monitoring program to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by {\em RoboPol} is analysed together with the gamma-ray data provided by {\em Fermi}-LAT. We confirm that polarization plane rotations are indeed related to the closest gamma-ray flares in blazars and the time lags between these events are consistent with zero. Amplitudes of the rotations are anticorrelated with amplitudes of the gamma-ray flares. This is presumably caused by higher relativistic boosting (higher Doppler factors) in blazars that exhibit smaller amplitude polarization plane rotations. Moreover, the time scales of rotations and flares are marginally correlated.

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.

Blazar population models have shown that the inverse-Compton and variability Doppler factor estimates yield consistent results at the population level for flat spectrum radio quasars (FSRQs). The two methods, however, are inconsistent when compared on a source-by-source basis. In this work, we attempt to understand the source of the discrepancy by tracing the potential sources of systematic and statistical error for the inverse-Compton Doppler factors. By eliminating these sources of error, we provide stronger constrains on the value of the Doppler factor in blazar jets. We re-estimate the inverse-Compton Doppler factor for 11 sources that meet certain criteria for their synchrotron peak frequency and the availability of Doppler factor estimates in the literature. We compare these estimates with the average of two different estimates of the variability Doppler factor obtained using various datasets and methodologies to identify any discrepancies and, in each case, trace their sources in the methodology or assumptions adopted. We identify three significant sources of error for the inverse-Compton Doppler factors: a) contamination of the X-ray flux by non-synchrotron self-Compton emission; b) radio observations at frequencies other than the synchrotron turnover frequency; c) non-simultaneity between radio and X-ray observations. We discuss key aspects in the correct application of the inverse-Compton method in light of these potential errors. We are able to constrain the Doppler factor of 3C273, 3C345, 3C454.3, PKS1510-089, and PKS1633+382 effectively, since all available estimates from both methods converge to the same values for these five sources.

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.

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

Narrow-line Seyfert 1 (NLS1) galaxies provide us with unique insights into the drivers of AGN activity under extreme conditions. Given their low black hole (BH) masses and near-Eddington accretion rates, they represent a class of galaxies with rapidly growing supermassive BHs in the local universe. Here, we present the results from our multi-frequency radio monitoring of a sample of {\gamma}-ray loud NLS1 galaxies ({\gamma}NLS1s), including systems discovered only recently, and featuring both the nearest and the most distant {\gamma}NLS1s known to date. We also present high-resolution radio imaging of 1H 0323+342, which is remarkable for its spiral or ring-like host. Finally, we present new radio data of the candidate {\gamma}-emitting NLS1 galaxy RX J2314.9+2243, characterized by a very steep radio spectrum, unlike other {\gamma}NLS1s.

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 RoboPol program has been monitoring the $R$-band linear polarisation parameters of an unbiased sample of 60 gamma-ray-loud blazars and a "control" sample of 15 gamma-ray-quite ones. The prime drive for the program has been the systematic study of the temporal behaviour of the optical polarisation and particularly the potential association of smooth and long rotations of the polarisation angle with flaring activity at high energies. Here we present the program and discuss a list of selected topics from our studies of the first three observing seasons (2013--2015) both in the angle and in the amplitude domain.

Detection of blazar pc scale counter-jets is difficult, but it can provide invaluable insight into the relativistic effects, radiative processes and the complex mechanisms of jet production, collimation and accelation in blazars. We build on recent populations models (optimized using the MOJAVE apparent velocity and redshift distributions) in order to derive the distribution of jet-to-counter-jet ratios and the flux densities of the counter-jet at different frequencies, in an effort to set minimum sensitivity limits required for existing and future telescope arrays in order to detect these elusive counter-jets. We find that: for the BL Lacs $5\%$ of their counter-jets have a flux-density higher than 100mJy, $15\%$ are higher than 10 mJy, and $32\%$ have higher flux-density than 1 mJy, whereas for the FSRQs $8\%$ have a flux-density higher than 10mJy, $17\%$ are higher than 1 mJy, and $32\%$ are higher than 0.1 mJy (at 15 GHz). Future telescopes like the SKA and newly operating like e-MERLIN and JVLA may detect up to $99\%$ of the BL Lac and $77\%$ of the FSRQ counter-jets. Sources with both low apparent velocity and a low Doppler factor make prime candidates for counter-jet detection. Combining our findings with literature values we have identified five such counter-jet detection candidates. Finally, we discuss possible effects beyond relativistic deboosting that may complicate the detection of counter-jets and that need to be accounted for in the interpretation of detections.

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3{\sigma} level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

We investigate the optical polarization properties of high-energy BL Lac objects using data from the RoboPol blazar monitoring program and the Nordic Optical Telescope. We wish to understand if there are differences in the BL Lac objects that are detected with the current-generation TeV instruments compared to those that have not yet been detected. The mean polarization fraction of the TeV-detected BL Lacs is 5% while the non-TeV sources show a higher mean polarization fraction of 7%. This difference in polarization fraction disappears when the dilution by the unpolarized light of the host galaxy is accounted for. The TeV sources show somewhat lower fractional polarization variability amplitudes than the non-TeV sources. Also the fraction of sources with a smaller spread in the Q/I - U/I -plane and a clumped distribution of points away from the origin, possibly indicating a preferred polarization angle, is larger in the TeV than in the non-TeV sources. These differences between TeV and non-TeV samples seems to arise from differences between intermediate and high spectral peaking sources instead of the TeV detection. When the EVPA variations are studied, the rate of EVPA change is similar in both samples. We detect significant EVPA rotations in both TeV and non-TeV sources, showing that rotations can occur in high spectral peaking BL Lac objects when the monitoring cadence is dense enough. Our simulations show that we cannot exclude a random walk origin for these rotations. These results indicate that there are no intrinsic differences in the polarization properties of the TeV-detected and non-TeV-detected high-energy BL Lac objects. This suggests that the polarization properties are not directly related to the TeV-detection, but instead the TeV loudness is connected to the general flaring activity, redshift, and the synchrotron peak location. (Abridged)

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.

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

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.

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $\mu m$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

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 measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

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

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.

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.

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.

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.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

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.

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

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.

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.

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.

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.

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.

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.

The high-frequency peaked blazars and known TeV emitters 1ES 1011+496, Mrk 180 and 1ES 2344+514 have been observed in the course of multi-wavelength campaigns in 2008, covering the frequency bands from radio up to TeV energies. For all three sources, these coordinated observations represent the first of their kind. We will present and discuss the campaigns, resulting light curves and spectral energy distributions.

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.

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 describe an upgrade to the Cosmic Background Imager instrument to increase its surface brightness sensitivity at small angular scales. The upgrade consisted of replacing the thirteen 0.9-m antennas with 1.4-m antennas incorporating a novel combination of design features, which provided excellent sidelobe and spillover performance for low manufacturing cost. Off-the-shelf spun primaries were used, and the secondary mirrors were oversized and shaped relative to a standard Cassegrain in order to provide an optimum compromise between aperture efficiency and low spillover lobes. Low-order distortions in the primary mirrors were compensated for by custom machining of the secondary mirrors. The secondaries were supported on a transparent dielectric foam cone to minimize scattering. The antennas were tested in the complete instrument, and the beam shape and spillover noise contributions were as expected. We demonstrate the performance of the telescope and the inter-calibration with the previous system using observations of the Sunyaev-Zel'dovich effect in the cluster Abell 1689. The enhanced instrument has been used to study the cosmic microwave background, the Sunyaev-Zel'dovich effect and diffuse Galactic emission.

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.

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

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

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

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 discovery by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope of high-energy gamma-ray emission from the peculiar quasar PMN J0948+0022 (z=0.585+/-0.001). Contrary to the expectations, the optical spectrum of this quasar shows only narrow lines [FWHM(Hbeta) ~ 1500 km/s] and the typical characteristics of narrow-line Seyfert 1 type galaxies. However, the strong radio emission and the flat spectrum suggest the presence of a relativistic jet, which can now be confirmed by the detection of MeV-GeV photons. PMN J0948+0022 is therefore the first radio-loud narrow-line Seyfert 1 quasar to be detected at gamma-rays and the third type of gamma-ray emitting AGN, after blazars and radiogalaxies.

We report the discovery by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope of high-energy gamma-ray emission from the peculiar quasar PMN J0948+0022 (z=0.5846). The optical spectrum of this object exhibits rather narrow Hbeta (FWHM(Hbeta) ~ 1500 km s^-1), weak forbidden lines and is therefore classified as a narrow-line type I quasar. This class of objects is thought to have relatively small black hole mass and to accrete at high Eddington ratio. The radio loudness and variability of the compact radio core indicates the presence of a relativistic jet. Quasi simultaneous radio-optical-X-ray and gamma-ray observations are presented. Both radio and gamma-ray emission (observed over 5-months) are strongly variable. The simultaneous optical and X-ray data from Swift show a blue continuum attributed to the accretion disk and a hard X-ray spectrum attributed to the jet. The resulting broad band spectral energy distribution (SED) and, in particular, the gamma-ray spectrum measured by Fermi are similar to those of more powerful FSRQ. A comparison of the radio and gamma-ray characteristics of PMN J0948+0022 with the other blazars detected by LAT shows that this source has a relatively low radio and gamma-ray power, with respect to other FSRQ. The physical parameters obtained from modelling the SED also fall at the low power end of the FSRQ parameter region discussed in Celotti & Ghisellini (2008). We suggest that the similarity of the SED of PMN J0948+0022 to that of more massive and more powerful quasars can be understood in a scenario in which the SED properties depend on the Eddington ratio rather than on the absolute power.

We present the results of the flux density measurements at 4.85 GHz and 10.45 GHz of a sample of 5998 NVSS radio sources with the Effelsberg 100 m telescope. The initial motivation was the need to identify the NVSS radio sources that could potentially contribute significant contaminating flux in the frequency range at which the Cosmic Background Imager experiment operated. An efficient way to achieve this challenging goal has been to compute the high frequency flux density of those sources by extrapolating their radio spectrum. This is determined by the three-point spectral index measured on the basis of the NVSS entry at 1.4 GHz and the measurements at 4.85 GHz and 10.45 GHz carried out with the 100 m Effelsberg telescope. These measurements are important since the targeted sample probes the weak part of the flux density distribution, hence the decision to make the data available. We present the table with flux density measurements of 3434 sources that showed no confusion allowing reliable measurements, their detection rates, their spectral index distribution and an interpretation which explains satisfactorily the observed uncertainties.

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.

Within the context of investigating possible differences between the mechanisms at play in Radio Loud AGN and those in Radio Quiet ones, we study the spectral characteristics of a selected sample of Intermediate-Luminosity Broad-Line Radio Galaxies in X-rays, optical, IR and radio. Here, we present the radio spectra acquired with the 100-m radio telescope in Effelsberg between 2.6 and 32 GHz. These measurements reveal a large variety of spectral shapes urging for radio imaging that would disclose the source morphology. Such studies could potentially discriminate between different mechanisms.

Gigahertz-Peaked Spectrum (GPS) sources are probably the precursors of local radio galaxies.Existing GPS source samples are small (<200). It is necessary to extend the availabe sample of the Gigahertz-Peaked Spectrum (GPS) and High Frequency Peaker (HFP) sources in order to study their nature with greater details and higher statistical significance. A sample of 214 radio sources, which were extracted from the SPECFIND catalog and show an inverted radio spectrum, were observed quasi-simultaneously at 4.85, 10.45, and 32GHz with the 100-m Effelsberg radio telescope. Using the VLBA calibrator survey (VCS) we have investigated the parsec-scale morphology of the sources. About 45% of the sources in our sample are classified as GPS or HFP candidates. We add 65 new GPS/HFP candidates to existing samples. We confirm the expected tendency that HFP are more compact on milliarcsecond scale than the 'classical' GPS sources, which peak at lower frequencies. The data mining of the SPECFIND database represents a promising tool for the discovery of new GPS/HFP sources.

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.

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