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.

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

Significant positional offsets of the value from 1 mas to more than 10 mas were found previously between radio (VLBI) and optical (Gaia) positions of active galactic nuclei (AGN). They happen preferentially parallel to the parsec-scale jet direction. AGN with VLBI-to-Gaia offsets pointed downstream the jet are found to have favourably higher optical polarization, as expected if extended optical jets dominate in the emission and shift the Gaia centroid away from the physical nucleus of the source. Upstream offsets with the suggested domination of accretion disks manifest themselves through the observed low optical polarization. Direction of linear optical polarization is confirmed to preferentially align with parsec-scale jets in AGN with dominant jets consistent with a toroidal magnetic field structure. Our findings support the disk-jet interpretation of the observed positional offsets. These results call on an intensification of AGN optical polarization monitoring programs in order to collect precious observational data. Taken together with the continued VLBI and Gaia observations, they will allow researchers to reconstruct detailed models of the disk-jet system in AGN on parsec scales.

Polarization measurements done using Imaging Polarimeters such as the Robotic Polarimeter are very sensitive to the presence of artefacts in images. Artefacts can range from internal reflections in a telescope to satellite trails that could contaminate an area of interest in the image. With the advent of wide-field polarimetry surveys, it is imperative to develop methods that automatically flag artefacts in images. In this paper, we implement a Convolutional Neural Network to identify the most dominant artefacts in the images. We find that our model can successfully classify sources with 98\% true positive and 97\% true negative rates. Such models, combined with transfer learning, will give us a running start in artefact elimination for near-future surveys like WALOP.

After three years of polarimetric monitoring of blazars, the RoboPol project has uncovered several key characteristics of polarimetric rotations in the optical for these most variable sources. The most important of these is that polarization properties of the synchrotron emission in the optical appear to be directly linked with gamma-ray activity. In this paper, we discuss the evidence for this connection, as well as the broader features of polarimetric behavior in blazars that are key in making progress with theoretical modeling of blazar emission.

The Fermi gamma-ray space telescope has revolutionized our view of the gamma-ray sky and the high energy processes in the Universe. While the number of known gamma-ray emitters has increased by orders of magnitude since the launch of Fermi, there is an ever increasing number of, now more than a thousand, detected point sources whose low-energy counterpart is to this day unknown. To address this problem, we combined optical polarization measurements from the RoboPol survey as well as other discriminants of blazars from publicly available all-sky surveys in machine learning (random forest and logistic regression) frameworks that could be used to identify blazars in the Fermi unidentified fields with an accuracy of >95%. Out of the potential observational biases considered, blazar variability seems to have the most significant effect reducing the predictive power of the frameworks to ~80-85%. We apply our machine learning framework to six unidentified Fermi fields observed using the RoboPol polarimeter. We identified the same candidate source proposed by Mandarakas et al. for 3FGL J0221.2+2518.

Galactic dust emission is polarized at unexpectedly high levels, as revealed by Planck. The origin of the observed $\simeq 20\%$ polarization fractions can be identified by characterizing the properties of optical starlight polarization in a region with maximally polarized dust emission. We measure the R-band linear polarization of 22 stars in a region with a submillimeter polarization fraction of $\simeq 20%$. A subset of 6 stars is also measured in the B, V and I bands to investigate the wavelength dependence of polarization. We find that starlight is polarized at correspondingly high levels. Through multiband polarimetry we find that the high polarization fractions are unlikely to arise from unusual dust properties, such as enhanced grain alignment. Instead, a favorable magnetic field geometry is the most likely explanation, and is supported by observational probes of the magnetic field morphology. The observed starlight polarization exceeds the classical upper limit of $\left[p_V/E\left(B-V\right)\right]_{\rm max} = 9$%mag$^{-1}$ and is at least as high as 13%mag$^{-1}$ that was inferred from a joint analysis of Planck data, starlight polarization and reddening measurements. Thus, we confirm that the intrinsic polarizing ability of dust grains at optical wavelengths has long been underestimated.

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.

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

The 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 availability of large datasets with stellar distance and polarization information will enable a tomographic reconstruction of the (plane-of-the-sky-projected) interstellar magnetic field in the near future. We demonstrate the feasibility of such a decomposition within a small region of the diffuse ISM. We combine measurements of starlight (R-band) linear polarization obtained using the RoboPol polarimeter with stellar distances from the second Gaia data release. The stellar sample is brighter than 17 mag in the R band and reaches out to several kpc from the Sun. HI emission spectra reveal the existence of two distinct clouds along the line of sight. We decompose the line-of-sight-integrated stellar polarizations to obtain the mean polarization properties of the two clouds. The two clouds exhibit significant differences in terms of column density and polarization properties. Their mean plane-of-the-sky magnetic field orientation differs by 60 degrees. We show how our tomographic decomposition can be used to constrain our estimates of the polarizing efficiency of the clouds as well as the frequency dependence of the polarization angle of polarized dust emission. We also demonstrate a new method to constrain cloud distances based on this decomposition. Our results represent a preview of the wealth of information that can be obtained from a tomographic map of the ISM magnetic field.

PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for inflationary gravitational waves, providing unique information regarding line-of-sight integration effects. Optical polarization observations of a large number of stars at known distances, tracing the same dust that emits polarized microwaves, can map the magnetic field between them. The Gaia mission is measuring distances to a billion stars, providing an opportunity to produce a tomographic map of Galactic magnetic field directions, using optical polarization of starlight. Such a map will not only boost CMB polarization foreground removal, but it will also have a profound impact in a wide range of astrophysical research, including interstellar medium physics, high-energy astrophysics, and evolution of the Galaxy. Taking advantage of the novel technology implemented in our high-accuracy Wide-Area Linear Optical Polarimeters (WALOPs) currently under construction at IUCAA, India, we will engage in a large-scale optopolarimetric program that can meet this challenge: a survey of both northern and southern Galactic polar regions targeted by CMB experiments, covering over 10,000 square degrees, which will measure linear optical polarization of over 360 stars per square degree (over 3.5 million stars, a 1000-fold increase over the state of the art). The survey will be conducted concurrently from the South African Astronomical Observatory in Sutherland, South Africa in the southern hemisphere, and the Skinakas Observatory in Crete, Greece, in the north.

Current models that explain giant (type II) X-ray outbursts in Be/X-ray binaries (BeXB), are based on the idea of highly distorted disks. They are believed to occur when a misaligned and warped disk becomes eccentric, allowing the neutron star to capture a large amount of material. The BeXB 4U 0115+63 underwent two major outbursts in 2015 and 2017. Our aim is to investigate whether the structural changes in the disk expected during type II outbursts can be detected through optical polarimetry. We present the first optical polarimetric observations and new optical spectra of the BeXB 4U 0115+63 covering the period 2013-2017. We study in detail the shape of the H$\alpha$ line profile and the polarization parameters before, during, and after the occurrence of a type II X-ray outburst. We find significant changes in polarization degree and polarization angle and highly distorted line profiles during the 2017 X-ray outburst. The degree of polarization decreased by $\sim$ 1%, while the polarization angle, which is supposed to be related with the disk orientation, first increased by $\sim 10^{\circ}$ in about two months and then decreased by a similar amount and on a similar timescale once the X-ray activity ceased.We interpret the polarimetric and spectroscopic variability as evidence for the presence of a warped disk.

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

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

We conduct a small-scale pathfinding survey designed to identify the average polarization properties of the diffuse ISM locally at the lowest dust content regions. We perform deep optopolarimetric surveys within three $\sim 15' \times 15'$ regions located at $b > 48^\circ$, using the RoboPol instrument. The observed samples of stars are photometrically complete to $\sim$16 mag in the R-band. The selected regions exhibit low dust emission at 353 GHz and low total reddening compared to the majority of high-latitude sightlines. We measure the level of systematic uncertainty for all observing epochs and find it to be 0.1\% in fractional linear polarization, $p$. The majority of individual stellar measurements are non-detections. However, our survey strategy enables us to locate the mean fractional linear polarization $p_{mean}$ in each of the three regions. The region with lowest dust content yields $p_{mean}=(0.054 \pm 0.038) \%$, not significantly different from zero. We find significant detections for the remaining two regions of: $p_{mean}=(0.113 \pm 0.036) \%$ and $p_{mean}=(0.208 \pm 0.044) \%$. Using a Bayesian approach we provide upper limits on the intrinsic spread of the small-scale distributions of $q$ and $u$. At the detected $p_{mean}$ levels, the determination of the systematic uncertainty is critical for the reliability of the measurements. We verify the significance of our detections with statistical tests, accounting for all sources of uncertainty. Using publicly available HI emission data, we identify the velocity components that most likely account for the observed $p_{mean}$ and find their morphologies to be misaligned with the orientation of the mean plane-of-sky magnetic field at a spatial resolution of 10$\arcmin$. We find indications that the standard upper envelope of $p$ with reddening underestimates the maximum $p$ at very low E(B-V) ($\leq 0.01$ mag).

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

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.

We identify rotations of the polarization angle in a sample of blazars observed for three seasons with the RoboPol instrument. A simplistic stochastic variability model is tested against this sample of rotation events. The model is capable of producing samples of rotations with parameters similar to the observed ones, but fails to reproduce the polarization fraction at the same time. Even though we can neither accept nor conclusively reject the model, we point out various aspects of the observations that are fully consistent with a random walk process.

We report on the variation in the optical polarization of the blazar PKS 1749+096 observed in 2008--2015. The degree of polarization (PD) tends to increase in short flares having a time-scale of a few days. The object favors a polarization angle (PA) of $40^\circ$--$50^\circ$ at the flare maxima, which is close to the position angle of the jet ($20^\circ$--$40^\circ$). Three clear polarization rotations were detected in the negative PA direction associated with flares. In addition, a rapid and large decrease in the PA was observed in the other two flares, while another two flares showed no large PA variation. The light curve maxima of the flares possibly tend to lag behind the PD maxima and color-index minima. The PA became $-50^\circ$ to $-20^\circ$ in the decay phase of active states, which is almost perpendicular to the jet position angle. We propose a scenario to explain these observational features, where transverse shocks propagate along curved trajectories. The favored PA at the flare maxima suggests that the observed variations were governed by the variations in the Doppler factor, $\delta$. Based on this scenario, the minimum viewing angle of the source, $\theta_\mathrm{min}=4.8^\circ$--$6.6^\circ$, and the location of the source, $\Delta r\gtrsim 0.1$pc, from the central black hole were estimated. In addition, the acceleration of electrons by the shock and synchrotron cooling would have a time-scale similar to that of the change in $\delta$. The combined effect of the variation in $\delta$ and acceleration/cooling of electrons is probably responsible for the observed diversity of the polarization variations in the flares.

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.

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

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

In any flux-density limited sample of blazars, the distribution of the timescale modulation factor $\Delta t'/\Delta t$, which quantifies the change in observed timescales compared to the rest-frame ones due to redshift and relativistic compression follows an exponential distribution with a mean depending on the flux limit of the sample. In this work we produce the mathematical formalism that allows us to use this information in order to uncover the underlining rest-frame probability density function of measurable timescales of blazar jets. We extensively test our proposed methodology using a simulated FSRQ population with a 1.5 Jy flux-density limit in the simple case (where all blazars share the same intrinsic timescale), in order to identify limits of applicability and potential biases due to observational systematics and sample selection. We find that for monitoring with time intervals between observations longer than $\sim$30\% of the intrinsic timescale under investigation the method loses its ability to produce robust results. For time intervals of $\sim$3\% of the intrinsic timescale the error of the method is as low as 1\% in recovering the intrinsic rest-frame timescale. We applied our method to rotations of the optical polarization angle of blazars observed by RoboPol. We found that the intrinsic timescales of the longest-duration rotation event in each blazar follows a narrow distribution, well-described by a normal distribution with mean 87 days and standard deviation 5 days. We discuss possible interpretations of this result.

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.

The main goal of this work is to perform a detailed study of the optical variability of the Be/X-ray binary SAX J2239.3+6116. We obtained multi-colour BVRI photometry and polarimetry and 4000-7000 A spectroscopy. The optical counterpart to SAX J2239.3+6116 is a V=14.8 B0Ve star located at a distance of ~4.9 kpc. The interstellar reddening in the direction of the source is E(B-V)=1.70 mag. The monitoring of the Halpha line reveals a slow long-term decline of its equivalent width since 2001. The line profile is characterized by a stable double-peak profile with no indication of large-scale distortions. Although somewhat higher than predicted by the models, the optical polarization is consistent with electron scattering in the circumstellar disk. We attribute the long-term decrease in the intensity of the Halpha line to the dissipation of the circumstellar disk of the Be star. The longer variability timescales observed in SAX J2239.3+6116 compared to other Be/X-ray binaries may be explained by the wide orbit of the system.

We present the first linear polarimetric survey of white dwarfs (WDs). Our sample consists of WDs of DA and DC spectral types in the SDSS r magnitude range from 13 to 17. We performed polarimetric observations with the RoboPol polarimeter attached to the 1.3-m telescope at the Skinakas Observatory. We have 74 WDs in our sample, of which almost all are low polarized WDs with polarization degree (PD) smaller than 1%, while only 2 have PD higher than 1%. There is an evidence that on average the isolated WDs of DC type have higher PD (with median PD of 0.78%) than the isolated DA type WDs (with median PD of 0.36%). On the other hand, the median PD of isolated DA type WDs is almost the same, i.e. 0.36% as the median PD of DA type white dwarfs in binary systems with red dwarfs (dM type), i.e. 0.33%. This shows, as expected, that there is no contribution to the PD from the companion if the WD companion is the red dwarf, which is the most common situation for WDs binary systems. We do not find differences in the polarization degree between magnetic and non-magnetic WDs. Because 97% of WDs in our sample have PD lower than 1%, they can be used as faint zero--polarized standard star in the magnitude range from 13 up to 17 of SDSS r filter. They cover the Northern sky between 13 hour to 23 hour in right ascension and from -11 degrees to 78 degrees in declination. Additionally, we found that for low extinction values (< 0.04) the best model that describes the dependence of PD on E(B-V) is given by the equation: PD_{max, ISM}[%] = 0.65 x E(B-V)^{0.12}.

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)

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.

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

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

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

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

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

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.

Polarization in classical Be stars results from Thomson scattering of the unpolarized light from the Be star in the circumstellar disc. Theory and observations agree that the maximum degree of polarization from isolated Be stars is < 4%. We report on the first optical polarimetric observations of the Be/X-ray binary EXO\,2030+375. We find that the optical (R band) light is strongly linearly polarized with a degee of polarization of 19%, the highest ever measured either in a classical or Be/X-ray binary. We argue that the interstellar medium cannot account for this high polarization degree and that a substantial amount must be intrinsic to the source. We propose that it may result from the alignment of non-spherical ferromagnetic grains in the Be star disc due to the strong neutron star magnetic field.

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

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.

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

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

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.

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.

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

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

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

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

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

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

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

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

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

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

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

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

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

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

We present results from monitoring the multi-waveband flux, linear polarization, and parsec-scale structure of the quasar PKS 1510-089, concentrating on eight major gamma-ray flares that occurred during the interval 2009.0-2009.5. The gamma-ray peaks were essentially simultaneous with maxima at optical wavelengths, although the flux ratio of the two wavebands varied by an order of magnitude. The optical polarization vector rotated by 720 degrees during a 5-day period encompassing six of these flares. This culminated in a very bright, roughly 1 day, optical and gamma-ray flare as a bright knot of emission passed through the highest-intensity, stationary feature (the "core") seen in 43 GHz Very Long Baseline Array images. The knot continued to propagate down the jet at an apparent speed of 22c and emit strongly at gamma-ray energies as a months-long X-ray/radio outburst intensified. We interpret these events as the result of the knot following a spiral path through a mainly toroidal magnetic field pattern in the acceleration and collimation zone of the jet, after which it passes through a standing shock in the 43 GHz core and then continues downstream. In this picture, the rapid gamma-ray flares result from scattering of infrared seed photons from a relatively slow sheath of the jet as well as from optical synchrotron radiation in the faster spine. The 2006-2009.7 radio and X-ray flux variations are correlated at very high significance; we conclude that the X-rays are mainly from inverse Compton scattering of infrared seed photons by 20-40 MeV electrons.

We use optical photometric and polarimetric data of BL Lacertae that cover a period of 22 years to study the variability of the source. The long-term observations are employed for establishing parameters of a stochastic model consisting of the radiation from a steady polarized source and a number of variable components with different polarization parameters, proposed by Hagen-Thorn et al. earlier. We infer parameters of the model from the observations using numerical simulations based on a Monte Carlo method, with values of each model parameter selected from a Gaussian distribution. We determine the best set of model parameters by comparing model distributions to the observational ones using the chi-square criterion. We show that the observed photometric and polarimetric variability can be explained within a model with a steady source of high polarization, ~40%, and with direction of polarization parallel to the parsec scale jet, along with 10+-5 sources of variable polarization.

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