The Exceptional Flare of Mrk 501 in 2014 Combined Observations with H.E.S.S. and FACT
The BL Lac type object Mrk 501 was observed at very high energies (E 100 GeV) in 2014 with the upgraded H.E.S.S. (High Energy Stereoscopic System) phase 2 array. The data collected with the central 28 m telescope allow for a broader energy range extending to lower energies when compared to the one obtained with the four small telescopes alone. A strong flaring event with a flux level comparable to the 1997 historical maximum has been detected as a consequence of target of opportunity observations triggered by alerts from the FACT collaboration. The First G-APD Cherenkov Telescope (FACT) is continuously monitoring bright blazars at TeV energies providing important pre- and post-flare information. For the first time, the data and lightcurves from H.E.S.S. and FACT are compared. These contemporaneous observations allow for a better characterization of the source emission. In a multiwavelength context, more precise correlation studies between VHE and lower energies are possible thanks to the dense sampling of the FACT observations. The hard intrinsic spectrum detected by H.E.S.S. during the flare allows the derivation of strong constraints on the scale of Lorentz invariance violation via the non-detection of EBL opacity modifications and from time-of-flight studies.
]G. Cologna ]N. Chakraborty aff3]A. Jacholkowska ]M. Lorentz aff1]M. Mohamed aff3,aff5]C. Perennes aff6]C. Romoli aff1]S. J. Wagner ]A. Wierzcholska for the H.E.S.S. Collaboration ]D. Dorner for the FACT Collaboration aff10,aff1]O. Kurtanidze
aff1]Landessternwarte, Universität Heidelberg, Königstuhl 12, D 69117 Heidelberg, Germany aff2]Max-Planck-Institut für Kernphysik, P.O. Box 103980, D 69029 Heidelberg, Germany aff3]LPNHE, Université Pierre et Marie Curie Paris 6, Université Denis Diderot Paris 7, CNRS/IN2P3, 4 Place Jussieu, F-75252, Paris Cedex 5, France aff4]18 DSM/Irfu, CEA Saclay, F-91191 Gif-Sur-Yvette Cedex, France aff5]LUTH, Observatoire de Paris, CNRS, Université Paris Diderot, 5 Place Jules Janssen, 92190 Meudon, France aff6]Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland aff7]Insitute of Nuclear Physics, Polish Academy of Sciences, Radzikowskiego 152, PL-31-342 Krakow, Poland aff8]Universität Würzburg, D-97074 Würzburg, Germany aff9]Universität Erlangen-Nürnberg, Physikalisches Institut, Erwin-Rommel-Str. 1, 91058 Erlangen, Germany aff10]Abastumani Observatory, Mt. Kanobili, 0301 Abastumani, Georgia
[cor1]email@example.com, firstname.lastname@example.org, email@example.com,
Mrk 501 ( 0.034) is a blazar belonging to the subclass of the high frequency peaked BL Lac objects (HBL). Its emission is strongly variable in all energy bands from radio up to very high energies (VHE, E100 GeV), where it was discovered as the second extragalactic source in 1995 . Its synchrotron and inverse Compton (IC) peaks are known to move towards higher energies during high flux states. In particular, the synchrotron peak can shift up to two orders of magnitude . As a consequence, a flux dependent spectral hardening has been observed in both X-rays and in the GeV-TeV band [e.g. 2, 3, 4, 5].
Mrk 501 has been extensively monitored over the years, often in multiwavelength (MWL) campaigns in order to characterize its spectral energy distribution (SED). In 1997, a prolonged high state was detected [e.g. 3, 6, 7, 8, 9, 10] and the historical highest flux was measured. Comparable high fluxes have been observed in short events on a few occasions, as in 2005 , 2012  and 2014 . During flares, the VHE flux variability can be as short as a few minutes [11, 13].
2 Observations and Analysis
VHE In 2014, the continuous monitoring of Mrk 501 with FACT  lead to the detection of several high state events with fluxes above a few Crab units (c.u.). In every one of these cases, flare alerts were issued, two of which triggered Target of Opportunity (ToO) observations with H.E.S.S.. A flare with flux comparable to the 1997 historical maximum was detected by both instruments on the night of June 23-24 2014. Mrk 501 was observed in wobble mode with the full HESS II array on several nights between June 19-25 and on July 29-30, and a total of 21 runs (observations with a typical exposure of 28 minutes) were collected. 18 runs pass the standard quality cuts, for a total livetime of 7.7 h (7.1 h when corrected for acceptance). The mean zenith angle is 63.8 causing the energy threshold of the observations to be very high, above 1 TeV. In , only events detected with the CT1-4 telescopes were extracted from the observations. Here, the information coming from all five telescopes are considered. The Model analysis  was used to reduce the stereo data. Very Loose cuts were used in order to obtain the lowest possible energy threshold. This choice does not raise any issue regarding background contamination since the background level is very low. An excess of almost 1900 photons is detected with high significance. Spectral analyses have been carried out on the total dataset and on the flare (four runs in one night) and low state (all runs in the remaining five nights) separately using the forward-folding technique . All results have been cross-checked using a different software and an independent calibration. The FACT monitoring provides valuable pre- and post-flare information. Thus, FACT observations between May 24 and August 4 are presented in this work alongside the H.E.S.S. observations. A total of 204 h in 63 nights have been collected. For the conversion to c.u., data with zenith angle smaller than 45 and only limited amount of moon light have been selected, leaving 144.5 hours of data in 59 nights. Rejecting bad weather data as described in  and excluding all data suffering from calima (the effect of calima on IACT data is described in ), 104 hours of data from 41 nights remain with an average zenith angle of 21.6. In the night of the flare, 4.1 hours of data were recorded. To allow for fast flare alerts, FACT has a quick-look analysis (described in ) providing results with low latency.
HE The Large Area Telescope (LAT) onboard the Fermi satellite provides the data at high energies (HE, 100 MeV E 100 GeV). The observations collected in a four months period centered on the VHE flare have been analysed with the ScienceTool software package version
v9r33p0. Events belonging to the ’Source’ class within 15 from the position of Mrk 501 were selected, and cuts on the zenith angle (100), rocking angle (52) and distance from the Sun (5) were applied. The instrument response functions
P7REP_SOURCE_V15 and the standard Isotropic and Galactic diffuse emission background models iso_source_v05.txt and
gll_iem_v05 rev1.fit were used for the binned maximum-likelihood spectral analysis. Mrk 501 was not sufficiently bright for sufficiently long time in order to derive short time-scale spectra or lightcurves with time bins shorter than one week.
X-rays X-ray data collected with Swift-XRT in WT mode between May 23 and August 6, 2014 (obsIDs: 00035023026-00035023078) have been analysed for this work using the HEASoft software package v. 6.16 (http://heasarc.gsfc.nasa.gov/docs/software/lheasoft) with CALDB v. 20140120. The events were cleaned and calibrated using
xrtpipeline and the data in the energy range 0.3-10 keV with grades 0-2 were analysed. For the spectral fit, the events were grouped in bins with at least 30 counts using the
grappha tool and were fit using XSPEC v. 12.8.2 with a single power-law model and Galactic hydrogen absorption fixed to cm . The lightcurve flux values (Figure 1) were calculated integrating the spectra of single snapshots between 2 and 10 keV. The relation between fluxes and spectral indices shows a clear harder-when-brighter behaviour, as expected from the literature (e.g. ). Most of the X-ray and FACT observations are simultaneous or at least contemporaneous, while X-ray and H.E.S.S. observations are typically offset by 90 minutes in the nights with coverage of both instruments.
Optical Observations in optical R band have been collected with the 70 cm telescope of the Abastunami Observatory (Georgia), which is equipped with an Apogee 6E camera. The data has been analysed with the Daophot II reduction software using an aperture diameter of 10 . Also in this case, nights of contemporaneous data with H.E.S.S. and FACT exist. Notably, strictly simultaneous observations have been taken during the flare peak on June 23, 2014.
3 Mwl Lightcurves
The regular monitoring of Mrk 501 by FACT ensures dense sampling, which is of primary importance in providing valuable pre- and post-flare information, as in the current case. The high sensitivity of H.E.S.S. on the other hand allows for minute dissection of the flare itself, probing timescales down to a few minutes.
In Figure 1, the multiwavelength lightcurves spanning 80 days around the June 23 flare are shown: VHE data from H.E.S.S. (black, E 2 TeV) and FACT (red, E 750 GeV) in the top panel, HE data from Fermi-LAT (0.1-500 GeV) in the second panel, X-ray data from Swift (2-10 keV) in the third panel and optical R-band data from Abastumani in the bottom panel. As can be seen in the top panel, the H.E.S.S. and FACT observations cover the flux rise from less than 1 c.u. before June 19 to 5 c.u. of the peak emission on the night of June 23, as well as the following sharp decline to sub c.u. levels.
Despite the different energy thresholds and not complete simultaneity, the average nightly fluxes recorded with H.E.S.S. and FACT in the same nights correlate well (Figure 2, left). In this period, there was no significant activity either in the optical or at HE. The comparison of the average nightly optical magnitudes with the fluxes measured by FACT (Figure 2, middle) shows no correlation. In particular, the optical emission is stable on the average value during the strictly simultaneous observations of the peak of the flare. The harder-when-brighter behaviour detected in X-rays can in principle explain these findings, since it requires the variability at lower energy to have a lower amplitude. On the other hand, it is known that the source shows activity also in optical. This suggests that at least two zones or mechanisms are necessary to explain the variability in VHE and in optical.
Figure 3 shows the flaring state on June 23 on minute level time scale. The H.E.S.S. and FACT lightcurves have bins of 4 and 10 minutes, respectively. There is evidence of fast variations on timescales of a few minutes at TeV energies. Doubling times below 10 minutes are obtained from the H.E.S.S. data. This favours a lepton induced emission at VHE. It is interesting to note, however, that the correlation expected between X-ray and TeV fluxes, indicative of synchrotron self-Compton (SSC) emission, is violated by the right-most point in Figure 2, right. Again, this suggests that the source state during the 2014 flare is different compared to the rest of the campaign in terms of the physical parameters of the emitting region or of the particles responsible for the -ray emission.
4 TeV SPECTRA
In  and , flux dependent spectral variability has been reported. In those publications, it was determined mostly at energies below 2 TeV. In the present work, the high energy threshold above 1 TeV and the very high maximal energy reached allow for a flux dependent spectral analysis at energies complementary to those of previous publications. Besides the total 2014 dataset, also the flare and the low state nights have been analysed separately.
Power law (PL) spectral fits provide very poor results, in line with the previous findings by, e.g., [3, 9, 11], and with the consideration that the intrinsic emitted spectrum must have been absorbed through the interaction with the extragalactic background light (EBL). In fact, despite Mrk 501 being a close source, a significant optical depth for EBL absorption 1 is expected at the highest energies observed from the source. Unsurprisingly, the curved power law (CPL) and the exponential cut-off power law (ECPL) spectral shapes fit the data satisfactorily. However, being the optical depth significant, the intrinsic spectra will differ from the observed ones. For this reason, fits including EBL absorption from  have been performed.
The spectra are best described by an EBL-absorbed PL. The photon indices are for the flaring state, and for the low state. The respective normalizations at 3.5 TeV are and 10cmsTeV. There is no evidence of intrinsic curvature nor intrinsic cut-off up to 25 TeV. These results show how a flux dependent spectral variability exists for energies extending significantly beyond 10 TeV, up to 20 TeV on a nightly time interval. The fact that the intrinsic spectrum is best described by a simple hard PL during the 2014 flare also indicates that there is no sign of being in the Klein-Nishina regime even at these high energies.
5 Lorentz Invariance Violation Studies
The 2014 flare data set with an energy coverage extending significantly up to 20 TeV has naturally been used to put constraints on Lorentz invariance violation (LIV).
Deviations from Lorentz invariance could appear at an energy scale beyond our current grasp. A generic approach to introduce such effects consists in adding effective terms in the dispersion relation of particles, i.e. for photons
where is the hypothetical energy scale at which Lorentz symmetry could stop being exact, and is the leading order of the perturbation. In some approaches to quantum gravity, is expected to be of the order of Planck energy eV .
Energetic astrophysical -rays can be used to constrain as Equation 1 would lead to non-negligible observational effects. The most direct approach in searching for LIV is to look for energy-dependent time delays in the arrival time of -rays. Constraints obtained with this approach using the 2014 flare data set are described in . These results are comparable with the best time-of-flight constraints for the linear case () derived from AGN observations , and are currently the best time-of-flight constraints for the quadratic case ().
An attractive alternative possibility to test Lorentz invariance takes advantage of the fact that the presence of LIV would affect the kinematics for the pair production of TeV -rays on the EBL, resulting in a modified opacity. In the subluminal case we consider here (minus sign in Equation 1), the threshold energy for the pair creation process is increased. This would reduce the EBL absorption and result in an excess of transparency for -rays at the highest energies with a specific spectral signature [see e.g. 25].
The spectral characteristics of Mrk 501 during flares and its low redshift allow for the detection of the most energetic extragalactic -rays, making this source ideal to investigate LIV through spectral studies, as has already pointed out by, e.g., [26, 27, 28] using the 1997 outburst, or by  discussing the sensitivity of CTA to LIV.
Optical depths using the EBL model of  are computed considering modifications due to LIV for the subluminal case. A log-likelihood spectral fit is performed on the H.E.S.S. 2014 flare dataset assuming an intrinsic power-law and letting the spectral index and the normalization free to vary. Values of are scanned logarithmically in the range of interest for linear and quadratic scenarios. More details can be found in . As the data show no evidence for a high-energy upturn, LIV-free optical depth values are preferred and log-likelihood profiles reach plateaus corresponding to the standard case. This allows to compute exclusion limits on , as summarized in Table 1.
|n=1||eV ()||eV ()||eV ()|
For the linear case, the 5 lower limit is at eV (). It is the best limit obtained using an AGN and confirms the result obtained with time delays using GRB 090510 : the standard photon dispersion relation holds up to the Planck energy scale in the case of linear perturbations. For the quadratic case, the 5 lower limit is at eV, which increases the current limits obtained with both AGNs and GRBs. These strong constraints naturally come from the exceptional spectrum of the 2014 flare data-set where the power law intrinsic emission extends up to 20 TeV with no sign of deviations from standard EBL absorption.
The same LIV analysis using the EBL model of  leads to very similar exclusion limits, showing that the results are not significantly affected by the choice of a specific EBL model. We have considered LIV affecting only photons [as in 28, 29] and not electrons since the constraints on LIV for electrons are very stringent . Relaxing this condition and considering LIV affecting electrons and photons equally amounts to a rescaling of by .
The H.E.S.S. Collaboration acknowledges: The support of the Namibian authorities and of the University of Namibia in facilitating the construction and operation of H.E.S.S. is gratefully acknowledged, as is the support by the German Ministry for Education and Research (BMBF), the Max Planck Society, the German Research Foundation (DFG), the French Ministry for Research, the CNRS-IN2P3 and the Astroparticle Interdisciplinary Programme of the CNRS, the U.K. Science and Technology Facilities Council (STFC), the IPNP of the Charles University, the Czech Science Foundation, the Polish Ministry of Science and Higher Education, the South African Department of Science and Technology and National Research Foundation, and by the University of Namibia. We appreciate the excellent work of the technical support staff in Berlin, Durham, Hamburg, Heidelberg, Palaiseau, Paris, Saclay, and in Namibia in the construction and operation of the equipment. The FACT Collaboration acknowledges: The important contributions from ETH Zurich grants ETH-10.08-2 and ETH-27.12-1 as well as the funding by the German BMBF (Verbundforschung Astro- und Energy Physics; the Laboratory Astroteilchenphysik) and HAP (Helmoltz Alliance for Astro-particle Physics) are gratefully acknowledged. We are thankful for the very valuable contributions from E. Lorenz, D. Renker and G. Viertel during the early phase of the project. We thank the Instituto de Astrofisica de Canarias allowing us to operate the telescope at the Observatorio del Roque de los Muchachos in La Palma, the Max-Planck-Institut fÃ¼r Physik for providing us with the mount of the ormer HEGRA CT3 telescope, and the MAGIC collaboration for their support. G.C. and M.M. are members of the International Max Planck Research School for Astronomy and Cosmic Physics at the University of Heidelberg (IMPRS-HD) and the Heidelberg Graduate School of Fundamental Physics (HGSFP). N.C. acknowledges support from Alexander von Humboldt foundation. The Abastumani Observatory team acknowledges financial support by the Shota Rustaveli NSF through project FR/577/6-320/13.
- J. Quinn et al., \apjl 456, p. L83 (1996).
- E. Pian et al., \apjl 492, L17–L20 (1998).
- A. Djannati-Atai et al., \aap 350, 17–24 (1999).
- A. A. Abdo et al., \apj 727, p. 129 (2011).
- G. Cologna et al., “Spectral characteristics of Mrk 501 during the 2012 and 2014 flaring states,” in 34th International Cosmic Ray Conference (ICRC2015), PoS(ICRC2015)761 (2015) arXiv:1509.04458 .
- M. Catanese et al., \apjl 487, L143–L146 (1997).
- F. W. Samuelson et al., \apjl 501, L17–L20 (1998).
- F. Aharonian et al., \aap 327, L5–L8 (1997).
- F. A. Aharonian et al., \aap 349, 11–28 (1999a).
- F. A. Aharonian et al., \aap 342, 69–86 (1999b).
- J. Albert et al., \apj 669, 862–883 (2007).
- F. Borracci et al., “Magic results and multi-wavelength observations of a mrk 501 flare in june 2012,” in 33rd International Cosmic Ray Conference (ICRC2013), The Astroparticle Physics Conference (2013).
- N. Chakraborty et al., “Rapid variability at very high energies in Mrk 501,” in 34th International Cosmic Ray Conference (ICRC2015), PoS(ICRC2015)872 (2015) arXiv:1509.04893 .
- H. Anderhub et al., Journal of Instrumentation 8, p. 6008P (2013).
- M. de Naurois and L. Rolland, Astroparticle Physics 32, 231–252 (2009).
- F. Piron et al., \aap 374, 895–906 (2001).
- D. Dorner et. al, “FACT - Long-term Monitoring of Bright TeV-Blazars,” in 33rd International Cosmic Ray Conference (ICRC2013), Braz. J. Phys., Vol. 44 (2014), pp. 415–608, arXiv:1311.0478 .
- D. Dorner, K. Nilsson, and T. Bretz, \aap 493, 721–725 (2009).
- D. Dorner et al., “FACT - Monitoring Blazars at Very High Energies,” in 5th Fermi Symposium Proceedings, eConf C141020.1 (2015), pp. 231–236, arXiv:1502.02582 .
- P. M. W. Kalberla et al., \aap 440, 775–782 (2005).
- O. M. Kurtanidze and M. G. Nikolashvili, “Blazar monitoring at Abastumani,” in Blazar Monitoring towards the Third Millennium, Proceedings of the OJ-94 Annual Meeting 1999, edited by C. M. Raiteri, M. Villata, and L. O. Takalo (1999), pp. 25–28.
- A. Franceschini, G. Rodighiero, and M. Vaccari, \aap 487, 837–852 (2008).
- G. Amelino-Camelia, Living Reviews in Relativity 16 (2013).
- H.E.S.S. Collaboration et al., Astroparticle Physics 34, 738–747 (2011).
- U. Jacob and T. Piran, \prd 78, p. 124010 (2008).
- F. W. Stecker and S. L. Glashow, Astroparticle Physics 16, 97–99 (2001).
- J. Biteau and D. A. Williams, \apj 812, p. 60 (2015).
- F. Tavecchio and G. Bonnoli, \aap 585, p. A25 (2016).
- M. Fairbairn et al., \jcap 6, p. 005 (2014).
- M. Lorentz and P. Brun for the H.E.S.S. Collaboration, “Limits on Lorentz invariance violation at the Planck energy scale from H.E.S.S. spectral analysis of the blazar Mrk 501,” in 6th Roma International Conference on AstroParticle Physics (RICAP16), EPJ Web of Conferences (2016) arXiv:1606.08600 .
- V. Vasileiou et al., \prd 87, p. 122001 (2013).
- A. Domínguez et al., \mnras 410, 2556–2578 (2011).
- T. Jacobson, S. Liberati, and D. Mattingly, \nat 424, 1019–1021 (2003).