We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the datas
et presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was used to process and calibrate the entire dataset, and we make the data publicly available for use by the astronomical community. Average quantities summarising optopolarimetric behaviour (average degree of polarization, polarization variability index) are also provided for each source we have observed and for the time interval we have followed it.
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 ca
ses 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.
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-d
ensity domain. We separate quiescent and flaring components of a sources 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 observatorys 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.
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 w
hether 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.
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 ev
ents, 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.
For a general dark-energy equation of state, we estimate the maximum possible radius of massive structures that are not destabilized by the acceleration of the cosmological expansion. A comparison with known stable structures constrains the equation
of state. The robustness of the constraint can be enhanced through the accumulation of additional astrophysical data and a better understanding of the dynamics of bound cosmic structures.
We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected fro
m the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.
Blazars are powerful, variable emitters from radio to gamma-ray wavelengths. Even though the general picture of synchrotron emission at low energies and inverse Compton at the high energy end is well established, many important aspects of these remar
kable objects are still not well understood. For example, even the location of the gamma-ray emission region is still not clearly established, with some theories locating it close to the black hole/accretion disk while others place it at parsec scales in the radio jet. Since mid-2007 we have carried out a large scale monitoring program at 15 GHz using the OVRO 40 m telescope. We are currently observing about 1700 blazars twice per week. The sample includes all the Fermi-LAT detected blazars north of declination -20 degrees. Here, we study the existence of correlated variability between these two bands for 86 sources bright enough to be detected weekly by LAT. The existence of correlated variability can be interpreted as an indication of a related spatial locations for the radio and gamma-ray emission, making the evaluation of its statistical significance a key goal of our program. A study of the statistical significance of these cross-correlations is presented along with a discussion of the Monte Carlo simulations used to evaluate them. More information about the conditions on the radio emission zone can be obtained through polarization monitoring which tells us about the configuration of the magnetic fields in this region. To study radio polarization variability we are building KuPol, a radio polarization receiver for the 12 to 18 GHz band that will replace the current total power receiver at the OVRO 40 meter telescope.
Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance
any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.
