We examine the 2008–2016 γ-ray and optical light curves of three bright BL Lac objects, 0716+714, MRK 421, BL Lac, which exhibit large structured variability. We searched for periodicities by using a ...fully Bayesian approach. For two out of three sources investigated, no significant periodic variability was found. In the case of BL Lac, we detected a periodicity of ~680 days. Although the signal related to this is modest, the coincidence of the periods in both gamma and optical bands is indicative of a physical relevance. Taking into consideration previous literature results, possibly related γ-ray and optical periodicities of about one year time scale are proposed in four bright γ-ray blazars out of the ten examined in detail. Compared with results from periodicity search of optical archives of quasars, the presence of quasi-periodicities in blazars may be more frequent by a large factor. This suggests the intriguing possibility that the basic conditions for their observability are related to the relativistic jet in the observer direction, but the overall picture remains uncertain.
Supermassive black holes with masses of millions to billions of solar masses are commonly found in the centers of galaxies. Astronomers seek to image jet formation using radio interferometry but ...still suffer from insufficient angular resolution. An alternative method to resolve small structures is to measure the time variability of their emission. Here we report on gamma-ray observations of the radio galaxy IC 310 obtained with the MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes, revealing variability with doubling time scales faster than 4.8 min. Causality constrains the size of the emission region to be smaller than 20% of the gravitational radius of its central black hole. We suggest that the emission is associated with pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the radio jet.
We present ∼2000 polarimetric and ∼3000 photometric observations of 15 γ-ray bright blazars over a period of 936 days (2008-10-11 to 2012-10-26) using data from the Tuorla blazar monitoring program ...(KVA DIPOL) and Liverpool Telescope (LT) RINGO2 polarimeters (supplemented with data from SkyCamZ (LT) and Fermi-LAT γ-ray data). In 11 out of 15 sources we identify a total of 19 electric vector position angle (EVPA) rotations and 95 flaring episodes. We group the sources into subclasses based on their broad-band spectral characteristics and compare their observed optical and γ-ray properties. We find that (1) the optical magnitude and γ-ray flux are positively correlated, (2) EVPA rotations can occur in any blazar subclass, four sources show rotations that go in one direction and immediately rotate back, (3) we see no difference in the γ-ray flaring rates in the sample; flares can occur during and outside of rotations with no preference for this behaviour, (4) the average degree of polarization (DoP), optical magnitude and γ-ray flux are lower during an EVPA rotation compared with during non-rotation and the distribution of the DoP during EVPA rotations is not drawn from the same parent sample as the distribution outside rotations, (5) the number of observed flaring events and optical polarization rotations are correlated, however we find no strong evidence for a temporal association between individual flares and rotations and (6) the maximum observed DoP increases from ∼10 per cent to ∼30 per cent to ∼40 per cent for subclasses with synchrotron peaks at high, intermediate and low frequencies, respectively.
ABSTRACT OJ 287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts that are predictable in a binary black hole model. The model predicted a major optical ...outburst in 2015 December. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, . The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2% accuracy level, and it opens up the possibility of testing the black hole no-hair theorem with 10% accuracy during the present decade.
In an attempt to constrain and understand the emission mechanism of γ-rays, we perform a cross-correlation analysis of 15 blazars using light curves in millimetre, optical and γ-rays. We use discrete ...correlation function and consider only correlations significant at the 99 per cent level. A strong correlation was found between 37 and 95 GHz with a near-zero time delay in most of the sources, and ∼1 month or longer in the rest. A similar result was obtained between the optical and γ-ray bands. Of the 15 sources, less than 50 per cent showed a strong correlation between the millimetre and γ-ray or millimetre and optical bands. The primary reason for the lack of statistically significant correlation is the absence of a major outburst in the millimetre bands of most of the sources during the 2.5 yr time period investigated in our study. This may indicate that only the long-term variations or large flares are correlated between these bands. The variability of the sources at every waveband was also inspected using fractional rms variability (F
var ). The F
var displays an increase with frequency reaching its maximum in the γ-rays.
Since the launch of the Fermi satellite, BL Lacertae has been moderately active at γ-rays and optical frequencies until 2011 May, when the source started a series of strong flares. The exceptional ...optical sampling achieved by the GLAST-AGILE Support Program of the Whole Earth Blazar Telescope in collaboration with the Steward Observatory allows us to perform a detailed comparison with the daily γ-ray observations by Fermi. Discrete correlation analysis between the optical and γ-ray emission reveals correlation with a time lag of 0 ± 1 d, which suggests cospatiality of the corresponding jet emitting regions. A better definition of the time lag is hindered by the daily gaps in the sampling of the extremely fast flux variations. In general, optical flares present more structure and develop on longer time-scales than corresponding γ-ray flares. Observations at X-rays and at millimetre wavelengths reveal a common trend, which suggests that the region producing the mm and X-ray radiation is located downstream from the optical and γ-ray-emitting zone in the jet. The mean optical degree of polarization slightly decreases over the considered period and in general it is higher when the flux is lower. The optical electric vector polarization angle (EVPA) shows a preferred orientation of about 15°, nearly aligned with the radio core EVPA and mean jet direction. Oscillations around it increase during the 2011-2012 outburst. We investigate the effects of a geometrical interpretation of the long-term flux variability on the polarization. A helical magnetic field model predicts an evolution of the mean polarization that is in reasonable agreement with the observations. These can be fully explained by introducing slight variations in the compression factor in a transverse shock waves model.
Abstract
The X-ray variability of the BL Lacertae source Mrk 501 was studied during 11.5 yr of monitoring with Swift. Here, we report the results of this study pertaining to the epoch of 2014 ...March–October, when our target showed the most powerful and long-lasting X-ray flaring activity. This epoch was characterized by X-ray flares varying in amplitude by factors of 2–5 on time-scales of a few weeks or shorter. We detected 35 instances of the intraday variability, sometimes occurring within the 1 ks observational runs. The X-ray flux was generally correlated with the TeV flux, while the 0.3–300 GeV and optical-UV fluxes did not show a significant correlation. Some notable incidences of more complicated variability patterns could also be recognized, indicating that the high-energy emission in Mrk 501 arose from an emission region more complex than a single zone. The best fits of the 0.3–10 keV spectra were mainly obtained using the logparabola model. Strong spectral variability was detected, affecting the slope but not the curvature of the spectrum. In strong flares, the spectral index became harder than 1.70. The spectral evolution was characterized by a harder-when-brighter behaviour, shifting the peak of the spectral energy distribution by about 20 keV that happens rarely in blazars.
Aims. We present a study of the very high-energy (VHE; E> 100 GeV) γ-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength ...spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z ≥ 0.6 and makes it a promising candidate to be the most distant VHE source. Methods. The source has been observed with the MAGIC telescopes in VHE γ rays for a total observation time of ~33.6 h from 2009 to 2011. A detailed analysis of its γ-ray spectrum and time evolution has been carried out. Moreover, we have collected and analyzed simultaneous and quasi-simultaneous multiwavelength data. Results. The source was marginally detected in VHE γ rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2σ. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 ± 1.2 to 5.0 ± 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009–2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z ~ 0.6.
We present ten years of R-band monitoring data of 31 northern blazars which were either detected at very high-energy (VHE) gamma rays or listed as potential VHE gamma-ray emitters. The data comprise ...11 820 photometric data points in the R-band obtained in 2002–2012. We analyzed the light curves by determining their power spectral density (PSD) slopes assuming a power-law dependence with a single slope β and a Gaussian probability density function (PDF). We used the multiple fragments variance function (MFVF) combined with a forward-casting approach and likelihood analysis to determine the slopes and perform extensive simulations to estimate the uncertainties of the derived slopes. We also looked for periodic variations via Fourier analysis and quantified the false alarm probability through a large number of simulations. Comparing the obtained PSD slopes to values in the literature, we find the slopes in the radio band to be steeper than those in the optical and gamma rays. Our periodicity search yielded one target, Mrk 421, with a significant (p < 5%) period. Finding one significant period among 31 targets is consistent with the expected false alarm rate, but the period found in Mrk 421 is very strong and deserves further consideration.
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 mag 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.