ABSTRACT We present coordinated multiwavelength observations of the bright, nearby BL Lacertae object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, ...VERITAS, and other collaborations and instruments, providing data from radio to very high energy (VHE) γ-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law, with no evidence for an exponential cutoff or additional hard components up to ∼80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure that relates to the two bumps of the broadband SED. In each bump, the variability increases with energy, which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in situ electron acceleration and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.
Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components that are spatially ...indistinguishable by gamma-ray instruments, but separated by a 10–12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma-ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data obtained by Fermi-LAT and MAGIC are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray source detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broadband emission can be modeled in the framework of a two-zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region.
Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy ...distribution (SED). Aims. The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. Methods. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E> 100 GeV) γ-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. Results. The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2–10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 104<γ< 6 × 105, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. Conclusions. Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population.
ABSTRACT The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 using the MAGIC telescopes. Together with the ...gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability timescale is estimated to be 6.4 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad-line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy.
Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse γ-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus ...cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E ≳ 100 GeV) for a total of 253 h from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6° from the centre, have been detected as point-like VHE γ-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is described well by a power law with a photon index Γ = 3.6 ± 0.2stat ± 0.2syst between 90 GeV and 1200 GeV. We do not detect any diffuse γ-ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes α. For α = 2.2, the CR-to-thermal pressure within the cluster virial radius is constrained to be ≲ 1−2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to ≲ 20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, B0, that depend on the CR distribution. For α = 2.2, B0 ≳ 5−8 μG, which is below the ~25 μG inferred from Faraday rotation measurements, whereas for α ≲ 2.1, the hadronic interpretation of the diffuse radio emission contrasts with our γ-ray flux upper limits independently of the magnetic field strength.
Aims. We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year timescales, which can bring additional insight on the processes responsible for its broadband emission. ...Methods. We observed Mrk 421 in the very high energy (VHE) γ-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 h of effective time after quality cuts. The VHE flux variability is quantified using several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes here for the first time. The 2.3 yr long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Metsähovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year timescales. Results. Mrk 421 was found in different γ-ray emission states during the 2.3 yr long observation period: The flux above 400 GeV spans from the minimum nightly value of (1.3 ± 0.4)×10-11 cm-2 s-1 to the maximum flux, that is about 24 times higher, at (3.1 ± 0.1)×10-10 cm-2 s-1. Flares and different levels of variability in the γ-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Metsähovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. Overall, the fractional variability increases with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and γ-ray photons are either the same, as expected in the synchrotron-self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.
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.
Context. The gamma-ray binary LS I +61°303 is a well-established source from centimeter radio up to very high energy (VHE; E> 100 GeV). The broadband emission shows a periodicity of ~26.5 days, ...coincident with the orbital period. A longer (super-orbital) period of 1667 ± 8 days was proposed from radio variability and confirmed using optical and high-energy (HE; E> 100 MeV) gamma-ray observations. In this paper, we report on a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long-timescale signature in the VHE emission from LS I +61°303. Aims. We focus on the search for super-orbital modulation of the VHE emission, similar to that observed at other energies, and on the search for correlations between TeV emission and an optical determination of the extension of the circumstellar disk. Methods. A four-year campaign has been carried out using the MAGIC telescopes. The source was observed during the orbital phases when the periodic VHE outbursts have occurred (φ = 0.55–0.75, one orbit = 26.496 days). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61°303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past (φ = 0.75–1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. Results. The TeV flux of the periodical outburst in orbital phases φ = 0.5–0.75 was found to show yearly variability consistent with the long-term modulation of ~4.5 years found in the radio band. This modulation of the TeV flux can be well described by a sine function with a best-fit period of 1610 ± 58 days. The complete data, including archival observations, span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star, but this may be affected by the strong, short-timescale (as short as intra-day) variation displayed by the Hα fluxes.
The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula has been ...claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 h of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 × 10-13 TeV cm-2 s-1 for the Geminga pulsar and 3.5 × 10-12 TeV cm-2 s-1 for the surrounding nebula at 50 GeV are the mostconstraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 yr of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.