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.
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.
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.
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.
MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the ...exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ∼ 50GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220GeV is (0.66 ± 0.03)% of Crab Nebula flux in 50h of observations. The angular resolution, defined as the σ of a 2-dimensional Gaussian distribution, at those energies is ≲ 0.07°, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11%–18% in flux normalization and ± 0.15 for the energy spectrum power-law slope.
The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. The telescopes are designed to measure Cherenkov light from air showers ...initiated by gamma rays in the energy regime from around 50GeV to more than 50TeV. The two telescopes were built in 2004 and 2009, respectively, with different cameras, triggers and readout systems. In the years 2011–2012 the MAGIC collaboration undertook a major upgrade to make the stereoscopic system uniform, improving its overall performance and easing its maintenance. In particular, the camera, the receivers and the trigger of the first telescope were replaced and the readout of the two telescopes was upgraded. This paper (Part I) describes the details of the upgrade as well as the basic performance parameters of MAGIC such as raw data treatment, linearity in the electronic chain and sources of noise. In Part II, we describe the physics performance of the upgraded system.
Context. The radio galaxy IC 310 has recently been identified as a γ-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E > 100 GeV) with the MAGIC ...telescopes. Originally classified as a head-tail radio galaxy, the nature of this object is subject of controversy since its nucleus shows blazar-like behavior. Aims. To understand the nature of IC 310 and the origin of the VHE emission, we studied the spectral and flux variability of IC 310 from the X-ray band to the VHE γ-ray regime. Methods. The light curve of IC 310 above 300 GeV has been measured with the MAGIC telescopes from 2009 October to 2010 February. Contemporaneous Fermi-LAT data (2008−2011) in the 10−500 GeV energy range were also analyzed. In the X-ray regime, archival observations from 2003 to 2007 with XMM-Newton, Chandra, and Swift-XRT in the 0.5−10 keV band were studied. Results. The VHE light curve reveals several high-amplitude and short-duration flares. Day-to-day flux variability is clearly present (>5σ). The photon index between 120 GeV and 8 TeV remains at the value Γ ~ 2.0 during both low and high flux states. The VHE spectral shape does not show significant variability, whereas the flux at 1 TeV changes by a factor of ~7. Fermi-LAT detected only eight γ-ray events in the energy range 10 GeV–500 GeV in three years of observation. The measured photon index of Γ = 1.3 ± 0.5 in the Fermi-LAT range is very hard. The X-ray measurements show strong variability in both flux and photon index. The latter varied from 1.76 ± 0.07 to 2.55 ± 0.07. Conclusions. The rapid variability measured in γ-rays and X-rays confirms the blazar-like behavior of IC 310. The multi-TeV γ-ray emission seems to originate from scales of less than 80 Schwarzschild radii (for a black hole mass of 2 × 108 M⊙) within the compact core of its FR I radio jet with orientation angle 10°−38°. The spectral energy distribution resembles that of an extreme blazar, albeit the luminosity is more than two orders of magnitude lower.
The MAGIC stereoscopic system collected 69 hours of Crab Nebula data between October 2009 and April 2011. Analysis of this data sample using the latest improvements in the MAGIC stereoscopic software ...provided an unprecedented precision of spectral and night-by-night light curve determination at gamma rays. We derived a differential spectrum with a single instrument from 50 GeV up to almost 30 TeV with 5 bins per energy decade. At low energies, MAGIC results, combined with Fermi-LAT data, show a flat and broad Inverse Compton peak. The overall fit to the data between 1 GeV and 30 TeV is not well described by a log-parabola function. We find that a modified log-parabola function with an exponent of 2.5 instead of 2 provides a good description of the data (χred2=35/26). Using systematic uncertainties of the MAGIC and Fermi-LAT measurements we determine the position of the Inverse Compton peak to be at (53±3stat+31syst−13syst) GeV, which is the most precise estimation up to date and is dominated by the systematic effects. There is no hint of the integral flux variability on daily scales at energies above 300 GeV when systematic uncertainties are included in the flux measurement. We consider three state-of-the-art theoretical models to describe the overall spectral energy distribution of the Crab Nebula. The constant B-field model cannot satisfactorily reproduce the VHE spectral measurements presented in this work, having particular difficulty reproducing the broadness of the observed IC peak. Most probably this implies that the assumption of the homogeneity of the magnetic field inside the nebula is incorrect. On the other hand, the time-dependent 1D spectral model provides a good fit of the new VHE results when considering a 80 μG magnetic field. However, it fails to match the data when including the morphology of the nebula at lower wavelengths.
We examine the 2008-2016 gamma-ray and optical light curves of three bright BLLac objects, 0716+714, MRK421, BLLac, 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 BLLac, 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 gamma-ray and optical periodicities of about one year time scale are proposed in four bright gamma-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.
Aims. We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it ...was identified as a potential very high energy (VHE; E> 100 GeV) gamma-ray emitter. We investigate whether this object is aVHE emitter, characterize its gamma-ray spectrum, and study the broadband emission within the one-zone synchrotron self-Compton (SSC) scenario, which is commonly used to describe the emission in blazars. Moreover, we also intend to determine the redshift of this object, which is a crucial parameter for its scientific interpretation. Methods. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 h of good quality stereoscopic data. Besides MAGIC, the campaign involved, observations with Fermi-LAT, Swift-XRT/UVOT, the optical telescopes KVA, Goddard Robotic Telescope, Galaxy View observatory, Crimean Astrophysical observatory, St. Petersburg observatory, and the Owens Valley Radio Observatory. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. We characterize the radio to VHE spectral energy distribution and quantify the multiband variability and correlations over short (few days) and long (many months) timescales. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. Results. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3σ (E > 70 GeV) during a 1.3 h long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. Besides the variability on few-day timescales, the long-term monitoring of MAGIC J2001+439 shows that, the gamma-ray, optical, and radio emissions gradually decreased on few-month timescales from 2010 through 2011, indicating that at least some of the radio, optical and gamma-ray emission is produced in a single region by the same population of particles. We also determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 ± 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 ± 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.