A fast, very-high-energy γ -ray flare from BL Lacertae during a period of multi-wavelength activity in June 2015

The mechanisms producing fast variability of the γ-ray emission in active galactic nuclei (AGNs) are under debate. The MAGIC telescopes detected a fast, very-high-energy (VHE, E > 100 GeV) γ-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of (1.5 ± 0.3) × 10-10 photons cm-2 s-1 and halving time of 26 ± 8 min. The MAGIC observations were triggered by a high state in the optical and high-energy (HE, E > 100 MeV) γ-ray bands. In this paper we present the MAGIC VHE γ-ray data together with multi-wavelength data from radio, optical, X-rays, and HE γ rays from 2015 May 1 to July 31. Well-sampled multi-wavelength data allow us to study the variability in detail and compare it to the other epochs when fast, VHE γ-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays, and HE γ-rays is very similar to two other observed VHE γ-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43 GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the γ-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE γ-ray flare, but all have some tension with the multi-wavelength observations

The broad-band properties of the intermediate synchrotron peaked BL Lac S2 0109+22 from radio to VHE gamma-rays

The Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes observed S2 0109+22 in 2015 July during its flaring activity in high-energy gamma-rays observed by Fermi-Large Area Telescope. We analyse the MAGIC data to characterize the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lacertae (BL Lac) objects. We study the multifrequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multifrequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of 5.3 sigma. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as z = 0.36 +/- 0.07. The SSC parameters describing the SED are rather typical for blazars

Detection of persistent VHE gamma-ray emission from PKS 1510-089 by the MAGIC telescopes during low states between 2012 and 2017

Context. PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. To date, very high-energy (VHE, >100 GeV) emission has been observed from this source either during long high states of optical and GeV activity or during short flares.Aims. We search for low-state VHE gamma-ray emission from PKS 1510-089. We characterize and model the source in a broadband context, which would provide a baseline over which high states and flares could be better understood.Methods. PKS 1510-089 has been monitored by the MAGIC telescopes since 2012. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray, and gammaray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modeled within an external Compton scenario with a stationary emission region through which plasma and magnetic fields are flowing. We also perform the emissio-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models.Results. The MAGIC telescopes collected 75 hr of data during times when the Fermi-LAT flux measured above 1 GeV was below 3 x10(-8)cm(-2) s(-1) , which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.50-sigma ) VHE gamma-ray emission at the level of (4.27 +/- 0.61(stat)) x 10(-12) cm(-2) s(-1) above 150 GeV, a factor of 80 lower than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broadband emission is compatible with the external Compton scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma ray state of a FSRQ. For the used model of the BLR, the 95% confidence level on the location of the emission region allows us to place it at a distance >74% of the outer radius of the BLR

Measurement of the extragalactic background light using MAGIC and Fermi-LAT gamma-ray observations of blazars up to z=1

We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03-0.944, obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible, and infrared bands. Major contributors to the EBL are the light emitted by stars through the history of the Universe, and the fraction of it that was absorbed by dust in galaxies and re-emitted at longer wavelengths. The EBL can be studied indirectly through its effect on very high energy photons that are emitted by cosmic sources and absorbed via gamma gamma interactions during their propagation across cosmological distances. We obtain estimates of the EBL density in good agreement with state-of-the-art models of the EBL production and evolution. The 1 sigma upper bounds, including systematic uncertainties, are between 13 per cent and 23 per cent above the nominal EBL density in the models. No anomaly in the expected transparency of the Universe to gamma-rays is observed in any range of optical depth. We also perform a wavelength-resolved EBL determination, which results in a hint of an excess of EBL in the 0.18-0.62 mu m range relative to the studied models, yet compatible with them within systematics

Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and γ-ray activity in 2015-2016

We report a characterization of the multiband flux variability and correlations of the nearby (z = 0.031) blazar Markarian 421 (Mrk 421) using data from Metsahovi, Swift, Fermi-LAT, MAGIC, FACT, and other collaborations and instruments from 2014 November till 2016 June. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F->1 TeV 0.1 TeV) gamma-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HRkeV and HRTeV show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 d centred at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 d centred at a time lag of 43(-6)(+9) d. This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a lognormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process

Unraveling the Complex Behavior of Mrk 421 with Simultaneous X-Ray and VHE Observations during an Extreme Flaring Activity in 2013 April*

We report on a multiband variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 19. The study uses, among others, data from GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), Swift, Nuclear Spectroscopic Telescope Array (NuSTAR), Fermi Large Area Telescope, Very Energetic Radiation Imaging Telescope Array System (VERITAS), and Major Atmospheric Gamma Imaging Cherenkov (MAGIC). The large blazar activity and the 43 hr of simultaneous NuSTAR and MAGIC/VERITAS observations permitted variability studies on 15 minute time bins over three X-ray bands (3-7 keV, 7-30 keV, and 30-80 keV) and three very-high-energy (VHE; >0.1 TeV) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV, and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all of the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-rays and the VHE gamma-rays. The three VHE bands and the three X-ray bands are positively correlated with no time lag, but the strength and characteristics of the correlation change substantially over time and across energy bands. Our findings favor multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer