The origin of high-energy emission in blazars jets (i.e., leptonic versus hadronic) has been a long-standing matter of debate. Here, we focus on one variant of hadronic models where proton synchrotron radiation accounts for the observed steady $\gamma$-ray blazar emission. Using analytical methods, we derive the minimum jet power ($P_{j,\min}$) for the largest blazar sample analyzed to date (145 sources), taking into account uncertainties of observables and jet's physical parameters. We compare $P_{j,\min}$ against three characteristic energy estimators for accreting systems, i.e., the Eddington luminosity, the accretion disk luminosity, and the power of the Blandford-Znajek process, and find that $P_{j,\min}$ is about 2 orders of magnitude higher than all energetic estimators for the majority of our sample. The derived magnetic field strengths in the emission region require either large amplification of the jet's magnetic field (factor of 30) or place the $\gamma$-ray production site at sub-pc scales. The expected neutrino emission peaks at $\sim 0.1-10$ EeV, with typical peak neutrino fluxes $\sim 10^{-4}$ times lower than the peak $\gamma$-ray fluxes. We conclude that if relativistic hadrons are present in blazar jets, they can only produce a radiatively subdominant component of the overall spectral energy distribution of the blazar's steady emission.

In this paper, we propose a new approach to determining cosmological distances to active galactic nuclei (AGN) via light travel-time arguments, which can be extended from nearby sources to very high redshift sources. The key assumption is that the variability seen in AGN is constrained by the speed of light and therefore provides an estimate of the linear size of an emitting region. This can then be compared with the angular size measured with very long baseline interferometry (VLBI) in order to derive a distance. We demonstrate this approach on a specific well studied low redshift (z = 0.0178) source 3C84 (NGC 1275), which is the bright radio core of the Perseus Cluster. We derive an angular diameter distance including statistical errors of $D_{A} = 72^{+5}_{-6}$ Mpc for this source, which is consistent with other distance measurements at this redshift. Possible sources of systematic errors and ways to correct for them are discussed.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

X-ray polarization should provide new probes of magnetic field geometry and acceleration physics near the base of blazar jets, but near-future missions will have limited sensitivity. We thus use existing lower energy data and X-ray variability measurements in the context of a basic synchro-Compton model to predict the X-ray polarization level and the probability of detection success for individual sources, listing the most attractive candidates for an IXPE campaign. We find that, as expected, several high-peak blazars such as Mrk 421 can be easily measured in 100 ksec exposures. Most low peak sources should only be accessible to triggered campaigns during bright flares. Surprisingly, a few intermediate peak sources can have anomalously high X-ray polarization and thus are attractive targets.

Even with several thousand {\it Fermi}-LAT blazar detections, the $\gamma$-ray emission mechanism is poorly understood. We explore correlated optical/$\gamma$-ray flux variations for 178 {\it Fermi}-LAT blazars regularly monitored by KAIT, SMARTS, and the Steward Observatory. Out of the 178 sources, 121 show a measurable ($>1\sigma$) discrete correlation function peak. Using the derived time-lags and Bayesian block light-curve decompositions, we measure the fraction of common and orphan flares between the two bands. After accounting for sampling and sensitivity limitations we quantify for the first time the true orphan flare rates of optical and $\gamma$-ray flares: 54.5\% of optical and 20\% of $\gamma$-ray flares are orphan events. Both the intraband temporal relation and the small orphan $\gamma$-ray flare fraction point toward leptonic processes as the likely mechanism for the high-energy emission. Motivated to discriminate between synchrotron self-Compton and external-Compton dominance in individual sources, we use the flux-flux variations to determine the slope $m$ of the ${\rm log}\,f_{\rm opt} - {\rm log}\,f_{\gamma}$ dependence. The slope distribution suggests a bimodal population with high and intermediate synchrotron-peak objects showing larger $m$ than low synchrotron-peak objects. We find that $m$ is naturally decreased through pollution from the orphan (typically optical) flares and develop a method to statistically recover, given the sources' measured orphan flare rate, the intrinsic $m$. While source classes show composite behavior, the majority of BL Lac objects favor $m=2$, indicating a synchrotron self-Compton origin for the $\gamma$-rays. No preference for either $m$ is found in flat spectrum radio quasars.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15.

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.

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

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

Relativistic effects dominate the emission of blazar jets complicating our understanding of their intrinsic properties. Although many methods have been proposed to account for them, the variability Doppler factor method has been shown to describe the blazar populations best. We use a Bayesian hierarchical code called {\it Magnetron} to model the light curves of 1029 sources observed by the Owens Valley Radio Observatory's 40-m telescope as a series of flares with an exponential rise and decay, and estimate their variability brightness temperature. Our analysis allows us to place the most stringent constraints on the equipartition brightness temperature i.e., the maximum achieved intrinsic brightness temperature in beamed sources which we found to be $\rm \langle T_{eq}\rangle=2.78\times10^{11}K\pm26\%$. Using our findings we estimated the variability Doppler factor for the largest sample of blazars increasing the number of available estimates in the literature by almost an order of magnitude. Our results clearly show that $\gamma$-ray loud sources have faster and higher amplitude flares than $\gamma$-ray quiet sources. As a consequence they show higher variability brightness temperatures and thus are more relativistically beamed, with all of the above suggesting a strong connection between the radio flaring properties of the jet and $\gamma$-ray emission.

Blazars are known for their energetic multiwavelength flares from radio wavelengths to high-energy $\gamma$-rays. In this work, we study radio, optical, and $\gamma$-ray light curves of 145 bright blazars spanning up to 8~yr, to probe the flaring activity and interband correlations. Of these, 105 show $>1\sigma$ correlations between one or more wavebands, 26 of which have a $>3\sigma$ correlation in at least one wavelength pair, as measured by the discrete correlation function. The most common and strongest correlations are found between the optical and $\gamma$-ray bands, with fluctuations simultaneous within our $\sim 30$~d resolution. The radio response is usually substantially delayed with respect to the other wavelengths with median time lags of $\sim 100$--160~d. A systematic flare identification via Bayesian block analysis provides us with a first uniform sample of flares in the three bands, allowing us to characterise the relative rates of multiband and "orphan" flares. Multiband flares tend to have higher amplitudes than "orphan" flares.

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

Context: Understanding the acceleration mechanism of astrophysical jets has been a cumbersome endeavor from both the theoretical and observational perspective. Although several breakthroughs have been achieved in recent years, on all sides, we are still missing a comprehensive model for the acceleration of astrophysical jets. Aims: In this work we attempt to construct a simple toy model that can account for several observational and theoretical results and allow us to probe different aspects of blazar jets usually inaccessible to observations. Methods: We used the toy model and Lorentz factor estimates from the literature to constrain the black hole spin and external pressure gradient distributions of blazars. Results: Our results show that (1) the model can reproduce the velocity, spin and external pressure gradient of the jet in M87 inferred independently by observations; (2) blazars host highly spinning black holes with 99% of BL Lac objects and 80% of flat spectrum radio quasars having spins a>0.6; (3) the dichotomy between BL Lac objects and Flat Spectrum Radio Quasars could be attributed to their respective accretion rates. Using the results of the proposed model, we estimated the spin and external pressure gradient for 75 blazars.

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.

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.

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.

Blazars are known to show periods of quiescence followed by outbursts visible throughout the electromagnetic spectrum. We present a novel maximum likelihood approach to capture this bimodal behavior by examining blazar radio variability in the flux-density domain. We separate quiescent and flaring components of a source's light curve by modeling its flux-density distribution as a series of "off" and "on" states. Our modeling allows us to extract information regarding the flaring ratio, duty cycle, and the modulation index in the "off"-state, in the "on"-state, as well as throughout the monitoring period of each blazar. We apply our method to a flux-density-limited subsample from the Owens Valley Radio observatory's 15 GHz blazar monitoring program, and explore differences in the variability characteristics between BL Lacs and FSRQs as well as between $\gamma$-ray detected and non-detected sources. We find that: (1) BL Lacs are more variable and have relatively larger outbursts than the FSRQs, (2) unclassified blazar candidates in our sample show similar variability characteristics as the FSRQs, and (3) $\gamma$-ray detected differ from the $\gamma$-ray non-detected sources in all their variability properties, suggesting a link between the production of $\gamma$-rays and the mechanism responsible for the radio variability. Finally, we fit distributions for blazar flaring ratios, duty cycles, and on- and off- modulation indices that can be used in population studies of variability-dependent blazar properties.

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.

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.

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)

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

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.

In a companion paper we presented a statistical model for the blazar population, consisting of distributions for the unbeamed radio luminosity function and the Lorentz factor distribution of each of the BL Lac and Flat Spectrum Radio Quasar (FSRQ) classes. Our model has been optimized so that it reproduces the MOJAVE distributions of apparent speeds and redshifts when the appropriate flux limit is applied and a uniform distribution of jet viewing angles is assumed for the population. Here we use this model to predict the Doppler factor distribution for various flux-limited samples for which Doppler factors have been estimated in a variety of ways (equipartition, variability + equipartition, inverse Compton dominance) on a blazar-by-blazar basis. By comparing the simulated and data-estimated Doppler factor distributions in each case, we evaluate the different methods of estimating blazar Doppler factors. We find that the variability Doppler factors assuming equipartition are the ones in the best agreement with our statistical model, whereas the inverse Compton Doppler factor method is only suitable for FSRQs. In the case of variability Doppler factors, we find that while random errors are relatively low ($\sim 30\%$), uncertainties are dominated by systematic effects. In the case of inverse Compton Doppler factors, random errors appear to dominate, but are significantly larger ($\sim 60\%$).

Observations of blazar jets are shrouded in relativistic effects, thus hindering our understanding of their intrinsic properties and dominant physical processes responsible for their generation, evolution,radiation and particle emission. In this work we focus on extracting information about timescales in the jet rest frame using a population-modeling approach. We employ Monte Carlo simulations to derive a simple population model for the intrinsic unbeamed luminosities and the Lorentz $\Gamma$ of blazar jets that adequately describe the observed redshift and apparent superluminal speed distributions for flux-limited blazar samples. We derive separate models for BL Lacs and Flat Spectrum Radio Quasars. We then use these models to compute the predicted distribution of Doppler factors in each blazar class, and address the following questions: (a) What is the relativistically induced spread in observed timescales (e.g., event duration, time lags) in a flux-limited sample of relativistic jets? (b) Could differences between BL Lacs and FSRQs observed in the time domain be attributed to differences in beaming between the two populations? (c) Is there a statistically preferred amount of beaming in a flux-limited sample? How large are statistical deviations from that preferred value? We use our findings to propose promising approaches in phenomenological studies of timescales in blazar jets.

In this study, we consider the problem of geometry of curved jets. We work with a simplified two-dimensional kinematical model assuming that the jet velocity is varying. We prescribe the velocity variation and study two cases concerning parsec and kilo-parsec scale jets by applying some simplified numerical simulations. For parsec scale jets, we use some general expressions and derive the shape of jets for certain numerical cases. For kilo-parsec scale jets, we study the cases of 4C53.37 and 1159+583 of Abell 2220 and Abell 1446 galaxy clusters, respectively, in an effort to test our model and compare our numerical results with observational data.