Abstract
We use results of our 3 yr polarimetric monitoring programme 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 RoboPol is analysed together with the gamma-ray data provided by 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 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 3s level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band luminosity 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 randomness 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 preferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.
Context. 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. Aims. We wish to ...understand if there are differences between the BL Lac objects that have been detected with the current-generation TeV instruments and those objects that have not yet been detected. Methods. We used a maximum-likelihood method to investigate the optical polarization fraction and its variability in these sources. In order to study the polarization position angle variability, we calculated the time derivative of the electric vector position angle (EVPA) change. We also studied the spread in the Stokes Q/I−U/I plane and rotations in the polarization plane. Results. 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 seem 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. Conclusions. 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.
ABSTRACT
We report results from our ongoing project MOMO (Multiwavelength Observations and Modelling of OJ 287). In this latest publication of a sequence, we combine our Swift UVOT–XRT and Effelsberg ...radio data (2.6–44 GHz) between 2019 and 2022.04 with public SMA data and gamma-ray data from the Fermi satellite. The observational epoch covers OJ 287 in a high state of activity from radio to X-rays. The epoch also covers two major events predicted by the binary supermassive black hole (SMBH) model of OJ 287. Spectral and timing analyses clearly establish: a new UV–optical minimum state in 2021 December at an epoch where the secondary SMBH is predicted to cross the disc surrounding the primary SMBH; an overall low level of gamma-ray activity in comparison to pre-2017 epochs; the presence of a remarkable, long-lasting UV–optical flare event of intermediate amplitude in 2020–2021; a high level of activity in the radio band with multiple flares; and particularly a bright, ongoing radio flare peaking in 2021 November that may be associated with a gamma-ray flare, the strongest in 6 yr. Several explanations for the UV–optical minimum state are explored, including the possibility that a secondary SMBH launches a temporary jet, but the observations are best explained by variability associated with the main jet.
Abstract
Based on our dedicated Swift monitoring program, MOMO, OJ 287 is one of the best-monitored blazars in the X-ray–UV–optical regime. Here, we report results from our accompanying, dense, ...multifrequency (1.4–44 GHz) radio monitoring of OJ 287 between 2015 and 2022 covering a broad range of activity states. Fermi
γ
-ray observations were added. We characterize the radio flux and spectral variability in detail, including discrete correlation function and other variability analyses, and discuss its connection with the multiwavelength emission. Deep fades of the radio and optical–UV fluxes are found to occur every 1–2 yr. Further, it is shown that a precursor flare of thermal bremsstrahlung predicted by one of the binary supermassive black hole (SMBH) models of OJ 287 was absent. We then focus on the nature of the extraordinary, nonthermal, 2016/2017 outburst that we initially discovered with Swift. We interpret it as the latest of the famous optical double-peaked outbursts of OJ 287, favoring binary scenarios that do not require a highly precessing secondary SMBH.
We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the 2013 May–June 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-mag of 17 and as low as 0.035 down to 18 mag. 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 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.
Context. With the current study we aim at understanding the properties of radio emission and the assumed jet from four radio-loud and γ-ray-loud narrow-line Seyfert 1 galaxies that have been detected ...by Fermi. These are Seyfert 1 galaxies with emission lines at the low end of the FWHM distribution. Aims. The ultimate goal is twofold: first we investigate whether a relativistic jet is operating at the source producing the radio output, and second, we quantify the jet characteristics to understand possible similarities with and differences from the jets found in typical blazars. Methods. We relied on the most systematic monitoring of radio-loud and γ-ray-detected narrow-line Seyfert 1 galaxies in the cm and mm radio bands conducted with the Effelsberg 100 m and IRAM 30 m telescopes. It covers the longest time-baselines and the most radio frequencies to date. This dataset of multi-wavelength, long-term radio light-curves was analysed from several perspectives. We developed a novel algorithm to extract sensible variability parameters (mainly amplitudes and time scales) that were then used to compute variability brightness temperatures and the corresponding Doppler factors. The jet powers were computed from the light curves to estimate the energy output and compare it with that of typical blazars. The dynamics of radio spectral energy distributions were examined to understand the mechanism causing the variability. Results. The length of the available light curves for three of the four sources in the sample allowed a firm understanding of the general behaviour of the sources. They all display intensive variability that appears to be occurring at a pace rather faster than what is commonly seen in blazars. The flaring events become progressively more prominent as the frequency increases and show intensive spectral evolution that is indicative of shock evolution. The variability brightness temperatures and the associated 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 recursive flaring activity. Finally, in one case, the sudden disappearance of a γ-ray flare below some critical frequency in our band needs a more detailed investigation of the possible mechanism causing the evolution of broadband events. Conclusions. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are powered by relativistic jets, such as intensive variability, spectral evolution across the different bands following evolutionary paths explained by travelling shocks, and Doppler factors indicating mildly relativistic speeds.
ABSTRACT
At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over ...a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-interval between successive observations – i.e. the cadence – is short enough to ensure that the correct modulo 180° value is selected. This optimal cadence depends on the maximum intrinsic EVPA rotation speed in blazars, which is currently not known. In this paper, we address the following questions for the RoboPol sample: What range of rotation speeds for rotations greater than 90° can we expect? What observation cadence is required to detect such rotations? Have rapid rotations been missed in EVPA rotation studies thus far? What fraction of data is affected by the ambiguity? And how likely are detected rotations affected by the ambiguity? We answer these questions with three seasons of optical polarimetric observations of a statistical sample of blazars sampled weekly with the RoboPol instrument and an additional season with daily observations. We model the distribution of EVPA changes on time-scales from 1–30 d and estimate the fraction of changes exceeding 90°. We show that at least daily observations are necessary to measure ${\gt}96{{\ \rm per\ cent}}$ of optical EVPA variability in the RoboPol sample of blazars correctly and that intraday observations are needed to measure the fastest rotations that have been seen thus far.
Context.
In February 2017 the blazar OJ 287, one of the best super-massive binary-black-hole-system candidates, was detected for the first time at very high energies (VHEs;
E
> 100 GeV) with the ...ground-based
γ
-ray observatory VERITAS.
Aims.
Very high energy
γ
rays are thought to be produced in the near vicinity of the central engine in active galactic nuclei. For this reason, and with the main goal of providing useful information for the characterization of the physical mechanisms connected with the observed teraelectronvolt flaring event, we investigate the parsec-scale source properties by means of high-resolution very long baseline interferometry observations.
Methods.
We use 86 GHz Global Millimeter-VLBI Array (GMVA) observations from 2015 to 2017 and combine them with additional multiwavelength radio observations at different frequencies from other monitoring programs. We investigate the source structure by modeling the brightness distribution with two-dimensional Gaussian components in the visibility plane.
Results.
In the GMVA epoch following the source VHE activity, we find a new jet feature (labeled K) at ∼0.2 mas from the core region and located in between two quasi-stationary components (labeled S1 and S2). Multiple periods of enhanced activity are detected at different radio frequencies before and during the VHE flaring state.
Conclusions.
Based on the findings of this work, we identify as a possible trigger for the VHE flaring emission during the early months of 2017 the passage of a new jet feature through a recollimation shock (represented by the model-fit component S1) in a region of the jet located at a de-projected distance of ∼10 pc from the radio core.