Long-term lightcurves from combined unified very high energy $γ$-ray data
Key Words.:Gamma-rays: observations – Galaxies: active – Galaxies: individual: Mrk 421 – Galaxies: individual: Mrk 501
Context: Very high-energy (VHE, E 100 GeV) -ray data are a valuable input for multi-wavelength and multi-messenger (e.g. combination with neutrino data) studies.
Aims: We aim at the conservation and homogenization of historical, current, and future VHE -ray-data on active galactic nuclei (AGN).
Methods: We have collected lightcurve data taken by major VHE experiments since 1991 and combined them into long-term lightcurves for several AGN, and now provide our collected datasets for further use. Due to the lack of common data formats in VHE -ray astronomy, we have defined relevant datafields to be stored in standard data formats. The time variability of the combined VHE lightcurve data was investigated, and correlation with archival X-ray data collected by RXTE/ASM tested.
Results: The combination of data on the prominent blazar Mrk 421 from different experiments yields a lightcurve spanning more than a decade. From this combined dataset we derive an integral baseline flux from Mrk 421 that must be lower than 33 % of the Crab Nebula flux above 1 TeV. The analysis of the time variability yields log-normal flux variations in the VHE-data on Mrk 421.
Conclusions: Existing VHE data contain valuable information concerning the variability of AGN and can be an important ingredient for multi-wavelength or multi-messenger studies. In the future, upcoming and planned experiments will provide more data from many transient objects, and the interaction of VHE astronomy with classical astronomy will intensify. In this context a unified and exchangeable data format will become increasingly important. Our data collection is available at the url: http://nuastro-zeuthen.desy.de/magic_experiment/projects/light_curve_archive/index_eng.html.
The broad-band emission observed from active galactic nuclei (AGN) spans the complete electromagnetic spectrum from radio to VHE (very high-energy, E 100 GeV) -rays. Since the discovery of the blazar Mrk 421 in the VHE regime (Punch et al. 1992), many new detections of AGNs by different experiments have been reported (see, e.g. Horns 2008; Hinton & Hofmann 2009, for reviews). Strong flux variations on different timescales were observed from many AGNs. Variability in the VHE -ray regime was measured down to the minute timescale (MAGIC Collaboration et al. 2008; Aharonian et al. 2007). Many AGN lightcurves were produced by different -ray experiments from observations of flaring states, dedicated monitoring programs, or from joint campaigns with other experiments. Long-term -ray lightcurves are valuable to address open questions concerning AGNs, such as possible periodicities (see Thieler et al. 2010, an analysis using the data collection of the present publication), the log-normality of their flux state distribution, the nature of the radiation mechanism of AGNs, or the estimation of their dutycycle.
A log-normal behavior can be indicative of a multiplicative process of the underlying mechanism governing the variability of the object. In the case of AGNs, this could be evidence for a connection of the observed emission to accretion disk activity (see Giebels & Degrange 2009, and references therein). Previously, log-normal flux variations were reported in the X-ray band from two objects: the narrow-line Seyfert 1 galaxy IRAS 13244-3809 (Gaskell 2004) and BL Lac (Giebels & Degrange 2009). In VHE data, log-normal variability was observed from a high flux state of the BL Lac object PKS 2155-304 (Degrange et al. 2008).
The nature of the AGN radiation mechanism in the VHE regime, i.e. whether the observed radiation has a leptonic or hadronic origin, still remains ambiguous. The detection of neutrinos from these objects would prove the existence of a hadronic component. In the framework of multi-messenger strategies (e.g., combining electromagnetic with neutrino data), the phenomenology of lightcurves in the electromagnetic wavelength band can give valuable input. The Neutrino triggered Target of Opportunity (NToO) program (Bernardini 2005; Ackermann et al. 2008) is based on the idea of neutrino events from a variable object (single events or multiplets) that are used to trigger VHE monitoring of the same object. Coincidences between neutrino triggers and -ray flux high states, which occur more often than expected from random coincidence with atmospheric neutrinos, would be evidence of a hadronic component of the VHE -ray signal and a cosmic origin of the neutrinos. Long-term -ray lightcurves can be used to estimate the AGN dutycycle, which is an important input parameter for such analyses.
2 Lightcurve data
Considering the heterogeneous nature of historical and present VHE data comprising different file and content formats from different experiments, a common data format is desirable. Moreover, as described above, VHE astronomy has already started to interface strongly with different fields of classical astronomy. In view of an effective exchange and diffusion of VHE data within the astronomy and astroparticle community, our storage strategy is to use an ascii file format (SLF) and the widely used FITS and VOTable file formats (see following section).
The use of a standard data format will also become increasingly important in the framework of the next-generation Cherenkov telescope systems, such as the upcoming CTA (Cherenkov Telescope Array) experiment (see e.g. Martinez 2008). In the CTA era, VHE astronomy will intensify multiwavelength (and multimessenger) interactions with other fields. Furthermore, a standard data format is essential for running an experiment like an observatory, and make it open to the whole astronomy community.
Publicly available lightcurve data from 1992 until today were collected from the Whipple (Kerrick et al. 1995; Schubnell et al. 1996; Buckley et al. 1996; Maraschi et al. 1999a), HEGRA (Aharonian et al. 1999b, a; Krawczynski et al. 2001b; Aharonian et al. 2001, 2002, 2003, 2004; Kestel 2002), CAT (Piron 2000; Piron et al. 2001), H.E.S.S. (Aharonian et al. 2005, 2006), MAGIC (Albert et al. 2008; Donnarumma et al. 2009), and VERITAS (Rebillot et al. 2006; Donnarumma et al. 2009) experiments. We are also working on collecting more data from these and other experiments such as the Crimean GT-48 Cherenkov telescope (e.g. Neshpor et al. 2007), the Tibet air-shower array (e.g. Amenomori et al. 2003), or the Patchmarhi Cherenkov Telescope Array (e.g. Gupta et al. 2008).
2.1 Data collection and formatting
Different formats and standards are used by different experiments. We use simple directly usable ascii tables and the standard astronomy file formats FITS and VOTable to store the collected, combined data. In the present work, the types and units of the FITS/VOTable datafields (i.e. columns) follow a simple lightcurve format (hereafter SLF) given in Table 1. The data are also provided in ascii file form following the SLF definition. These files are referred to as slf-files. The tables contain rows of (typically nightwise) integral flux-measurements and further information necessary for subsequent analyses. The units of the datafields are defined within the FITS/VOTable files and follow the conventions defined in Table 1 within the slf-files (ascii).
|1||mjd_mid_exp||MJD||center of observation exposure|
|2||mjd_start||MJD||start of observation|
|3||mjd_end||MJD||end of observation|
|4||int_flux||Crab units||Observed integral flux|
|5||sigma_int_flux_stat||Crab units||Statistical error on|
|6||sigma_int_flux_sys||Crab units||Systematical error on|
|7||alpha||none||Differential spectral index|
|8||sigma_alpha_stat||none||Statistical error on|
|9||sigma_alpha_sys||none||Systematical error on|
|10||e_thr||TeV||Energy threshold in TeV|
|11||e_cut||TeV||Cutoff energy in TeV|
|13||Duration||duration||hours||Duration of observation|
|14||Reference||reference||none||Reference string (pref. ADS format)|
|15||Flux-flag||fflag||none||flux-measurement: ’’, upper limit: ’’|
|16||additional entry 1||add_entry1||e.g., special parametrizations|
|17||additional entry 2||add_entry2|
The SLF datafield format definition provides the means to effectively combine heterogeneous datasets. In case further information than included in columns 1 to 12 are necessary, additional entries are foreseen. For example, to define a spectral shape different from the standard pure powerlaw, the spectral parameterization together with the additional parameters must be defined as comments within the data files. In case of lightcurve data only containing information on a few datafields, the flexibility of the binary FITS and VOTable file format allows minimizing the data volume. For example, if only the observation date and the flux value are known for the total lightcurve data of an object, the binary file can only contain and , along with the definition of units used. While this example represents the simplest case of a homogeneous dataset, in practice, the heterogeneous lightcurve data combinations from different experiments can include information on a specific datafield for one subset of the combined data (i.e. for one experiment), and not for another subset. In this (most common) case, when information is not available for datafields of some of the measurements, the following conventions for default datafield values apply.
If no start and end times are known, then is set to the same value as , i.e. the middle of the exposure time, and to -1. In this case the duration might still be given.
In some cases only the MJD of the observations day is known. Then, is the MJD of the observation day.
Most energy spectra can be described by a pure powerlaw (differential spectral index ) or a powerlaw with an exponential cutoff (cutoff energy ). Additional fields are foreseen for cases with different parameterizations of the energy spectrum. If such a different parameterization is given, a defining formula containing the relevant fields has to be written in the comment section of the lightcurve file.
The flux-flag (fflag) is a one-character flag indicating whether the integral flux is an upper limit (fflag = ’’) or a flux measurement (fflag = ’=’).
In the present work, we chose Crab Nebula flux units for the observed integral flux above the energy threshold .
For any other entry, the value -1 means that no information is available.
2.2 Lightcurve combination
For a combination of lightcurves from different experiments, the measured integral flux values must be converted to a common energy threshold. For the results in this paper, we chose an energy threshold of 1 TeV. The conversion of integral fluxes in Crab units to the same energy threshold requires knowledge of the energy spectrum of the considered object and of the Crab Nebula as measured by the same experiment. The differential energy spectrum of the Crab Nebula in the energy range covered by most VHE -ray experiments is described well by a pure powerlaw of the form in the energy range from 100 GeV to few tens of TeV. Deviations from the pure powerlaw form at both ends of this energy range (as expected in the framework of the inverse Compton scenario) were not taken into account here. Extrapolation to energies below 100 GeV requires a full parameterization of the spectral energy distribution of the Crab Nebula (e.g. Aharonian et al. 2004). An integral flux in Crab units above the energy threshold is given by
where and are the differential energy spectra of the source and of the Crab Nebula (in units of photons TeVcms). The corresponding integral fluxes above a given energy threshold are denoted with . In the case of a pure powerlaw differential energy spectrum , the conversion of an observed integral flux to a given threshold energy is given by
where is the energy threshold of the observation.
The combined day-wise, integral flux lightcurve of the BL Lac object Mrk 421 above 1 TeV is shown in Figure 1. This lightcurve includes all data we have collected so far, covering an uprecedented 17-year time-span from 1992 to 2008.
In Figure 2, the day-wise integral flux lightcurve of the BL Lac object Mrk 501 is shown.
2.3 Systematic errors
The inhomogeneity of data combined from different experiments and the inclusion of data from the pioneering time of VHE -ray astronomy induces systematic errors that are not easily evaluated. Combination of data from many different experiments, different time periods, partly inaccurate flux normalizations (i.e. -ray rate measurements) and energy thresholds can lead to systematic errors in the relative flux normalizations. When combining data, the flux normalization depends on assumptions on the spectral shape and the energy threshold. Here, limiting factors can be experimental uncertainties (especially in old data) induced by poor knowledge of atmospheric conditions, aging detector effects, and intrinsic spectral variability of the considered objects. We convert flux and rate measurements to Crab Nebula flux units using Crab Nebula data taken by the same experiment as close in time as possible, therewith reducing systematic uncertainties induced by seasonal or instrumental variations.
The error on the reconstructed spectral shape results in an error on the flux state when converting to different thresholds. For example, when converting integral fluxes from 1 TeV to 100 GeV, typical systematic errors on the spectral index of 0.1 (pure powerlaw case) lead to a relative systematic error on the integral flux of 30 % due to this extrapolation. The datafields contained in our FITS files provide all information needed to estimate these systematic errors for conversion to different energy thresholds. A conversion over a narrower energy range than used in this example limits this systematic error to less than 30 %.
The flux units used in some publications are given in uncalibrated units (e.g. counts per minute), i.e. as an instrument-specific -ray count-rate. The energy threshold is also unknown in some cases. For flux calibration of lightcurves from the Whipple experiment before 1998, we used Crab Nebula data taken by the same experiment as close in time as possible (see Hillas et al. 1998, and references therein). However, due to dependencies on zenith angle and weather conditions, normalization of rate measurements is very uncertain. We estimate that the systematic error induced when including rate measurements is below 40 %. This form of relative systematic error can be avoided when restricting analysis to flux measurements, by restricting to data from a single experiment, or by intercalibrating the different measurements when observations from different experiments overlap.
Significant evidence of a correlated gamma- and X-ray emission of blazars was presented in earlier studies (see e.g. Takahashi et al. 1996; Buckley et al. 1996; Maraschi et al. 1999b; Krawczynski et al. 2001a; Neshpor et al. 2007). Usually the reported correlation is linear, but in a few cases a quadratic relation between the fluxes in both bands was found (see e.g. Krawczynski et al. 2000). Such a correlation provides us with essential information on the underlying acceleration and emission processes and is especially valuable for variability studies. Very often the gamma- X-ray correlation is interpreted as a strong argument in favor of the so-called synchrotron-Compton jet emission models in which the same population of ultrarelativistic electrons is responsible for production of both X-rays via synchrotron radiation and TeV -rays via inverse Compton scattering (Katarzyński et al. 2005). However, it can also be accommodated in the hadronic framework in particular in models which assume that the observed -ray emission is a result of interactions of accelerated protons and ambient gas or low-frequency radiation (Aharonian 2000; Mücke et al. 2003).
Using contemporanous X-ray data extracted from the RXTE/ASM database web interface at
In Figure 3, the distribution of the VHE flux states of Mrk 421
(also see Figure 1) is shown. Each entry of this histogram represents a snapshot flux measurement of Mrk 421 as given by the combined dataset. All individual observed integral flux values were converted to flux values in units of the Crab Nebula flux and normalized to a common energy threshold of 1 TeV. The observed energy spectra of the Crab Nebula and Mrk 421 as observed by the individual experiments were taken into account. The overall distribution can be described by an exponential law (left panel) above a flux of a few tenths of Crab. An exponential shape might indicate that the measured flux states mainly reflect a stochastic outburst state of the object. As shown in the right panel of Figure 3, a better fit to the data is obtained when using the sum of a Gaussian and a log-normal distribution (Aitchison & Brown 1957; Limpert et al. 2001), as given by
The data were further divided into time intervals of equal length, each comprising 40 flux observations with statistical errors . The excess variance (Vaughan et al. 2003) was calculated for each time interval. The excess variance is a measure of the Poisson noise corrected rms value of the fluxes within the corresponding time interval. Figure 4 shows as a function of the average flux within the interval. A clear proportionality is seen with . A log-normally shaped distribution with a proportionality of the excess variance with the flux are evidence of log-normal flux variations (Aitchison & Brown 1957).
Below a flux level of few tenths of Crab, one might expect a low baseline flux level that can be described by a Gaussian. At these low fluxes, however, detector sensitivity threshold effects become increasingly important, i.e. lower fluxes cannot be detected significantly within the short observation windows. Therefore, the mean of the Gaussian fit (0.33 Crab flux units) must be treated as an absolute upper limit on the integral baseline flux above 1 TeV from Mrk 421.
More information such as the variability type (e.g. red noise / white noise / blue noise) or the duty-cycle of the object can be extracted from the flux state distributions. Previously, first steps in this direction were taken by Tluczykont et al. (2007); however, further investigations are necessary using simulations and will be the topic of a subsequent publication.
With increasing exposure achieved by the major current generation experiments more data will soon be available. As discussed previously, one important aspect for long-term variability studies will be to carry out unbiased (random) observations to avoid systematic selection effects.
4 Summary & outlook
VHE -ray flux measurements of AGN have been collected from observations of
different experiments since 1992.
For the first time these data were combined into single long-term lightcurves,
and are provided for further analysis in the standard FITS file format.
The BL Lac object Mrk 421 yields the most extended dataset with a combined
lightcurve spanning more than a decade.
The collected data is publicly available for download
The observed flux states averaged over 12 h of Mrk 421 and Mrk 501 above 1 TeV are consistent within statistical and systematic errors with a linear correlation with daily averaged RXTE/ASM X-ray count rates, with correlation coefficients of 0.65 and 0.68 respectively.
The long-term flux state distribution of Mrk 421 can be well described by the sum of a Gaussian and a log-normal distribution. This behavior can be interpreted by the observation of a baseline flux and a stochastic flux state population governed by an underlying multiplicative process. In the framework of this interpretation, the combined dataset on Mrk 421 allowed setting an absolute upper limit on the baseline flux of Mrk 421 of 0.33 Crab flux units above 1 TeV.
Unbiased data from monitoring campaigns by ongoing experiments (for a review see Weekes 2008) and data from future experiments such as CTA will further extend the long-term coverage of variable objects. A unified standard data format will become increasingly important in the future, when more high-sensitivity VHE data becomes available, and when interfacing with other fields of classical astronomy and astroparticle physics (multi-wavelength and multi-messenger) intensifies. We intend to store these future data as described in this work in standard data formats using the introduced SLF field structure definition.
Acknowledgements.We thank all colleagues who provided data in electronic form to us. MT thanks Dr. B. Giebels and Prof. Dr. D. Horns for valuable discussions. We thank Dr. B. Pence for helpful comments on an ealier draft. This research has made use of NASA’s Astrophysics Data System Bibliographic Services.
- Ackermann, M., Bernardini, E., Galante, N., & et al. 2008, in International Cosmic Ray Conference, Vol. 3, International Cosmic Ray Conference, 1257–1260
- Aharonian, F., Akhperjanian, A., Barrio, J., et al. 2001, ApJ, 546, 898
- Aharonian, F., Akhperjanian, A., Beilicke, M., et al. 2002, A&A, 393, 89
- Aharonian, F., Akhperjanian, A., Beilicke, M., et al. 2004, ApJ, 614, 897
- Aharonian, F., Akhperjanian, A., Beilicke, M., et al. 2003, A&A, 410, 813
- Aharonian, F., Akhperjanian, A. G., Aye, K., et al. 2005, A&A, 437, 95
- Aharonian, F., Akhperjanian, A. G., Barrio, J. A., et al. 1999a, A&A, 349, 29
- Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2007, ApJ, 664, L71
- Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 457, 899
- Aharonian, F. A. 2000, New A, 5, 377
- Aharonian, F. A., Akhperjanian, A. G., Barrio, J. A., et al. 1999b, A&A, 342, 69
- Aitchison, J. & Brown, J. A. C. 1957, The Lognormal Distribution, cambridge University Press: Cambridge
- Albert, J., Aliu, E., Anderhub, H., et al. 2008, ApJ, 674, 1037
- Amenomori, M., Ayabe, S., Cui, S. W., et al. 2003, ApJ, 598, 242
- Bernardini, E. 2005, in Proceedings Towards A Network Of Cherenkov Detectors VII, ed. B. Degrange, Palaiseau, 307–315
- Buckley, J. H., Akerlof, C. W., Biller, S., et al. 1996, ApJ, 472, L9+
- Degrange, B., Superina, G., Giebels, B., & Volpe, F. 2008, in Blazar Variability across the Electromagnetic Spectrum
- Donnarumma, I., Vittorini, V., Vercellone, S., et al. 2009, ApJ, 691, L13
- Gaskell, C. M. 2004, ApJ, 612, L21
- Giebels, B. & Degrange, B. 2009, A&A, 503, 797
- Gupta, A. C., Acharya, B. S., Bose, D., Chitnis, V. R., & Fan, J. 2008, Chinese J. Astron. Astrophys., 8, 395
- Hillas, A. M., Akerlof, C. W., Biller, S. D., et al. 1998, ApJ, 503, 744
- Hinton, J. A. & Hofmann, W. 2009, ARA&A, 47, 523
- Horns, D. 2008, in Reviews in Modern Astronomy, ed. S. Röser, Vol. 20, 167–+, arXiv:astro-ph/0808.3744
- Katarzyński, K., Ghisellini, G., Tavecchio, F., et al. 2005, A&A, 433, 479
- Kerrick, A. D., Akerlof, C. W., Biller, S. D., et al. 1995, ApJ, 438, L59
- Kestel, M. 2002, phD thesis, Technische Universität München, Germany
- Krawczynski, H., Coppi, P. S., Maccarone, T., & Aharonian, F. A. 2000, A&A, 353, 97
- Krawczynski, H., Coppi, P. S., Maccarone, T., & Aharonian, F. A. 2001a, X-ray Astronomy: Stellar Endpoints, AGN, and the Diffuse X-ray Background, 599, 694
- Krawczynski, H., Sambruna, R., Kohnle, A., et al. 2001b, ApJ, 559, 187
- Limpert, E., Stahel, W., & Abbt, M. 2001, BioScience, 51 (5), 341â352
- MAGIC Collaboration, Albert, J., Aliu, E., et al. 2008, Physics Letters B, 668, 253
- Maraschi, L., Fossati, G., Tavecchio, F., et al. 1999a, ApJ, 526, L81
- Maraschi, L., Fossati, G., Tavecchio, F., et al. 1999b, Astroparticle Physics, 11, 189
- Martinez, M. 2008, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 824–829
- Mücke, A., Protheroe, R. J., Engel, R., Rachen, J. P., & Stanev, T. 2003, Astroparticle Physics, 18, 593
- Neshpor, Y. I., Eliseev, V. S., Jogolev, N. A., et al. 2007, Bulletin Crimean Astrophysical Observatory, 103, 16
- Piron, F. 2000, PhD thesis, AA(Université de Paris-Sud (Paris XI))
- Piron, F., Djannati-Atai, A., Punch, M., et al. 2001, A&A, 374, 895
- Punch, M., Akerlof, C. W., Cawley, M. F., et al. 1992, Nature, 358, 477
- Rebillot, P. F., Badran, H. M., Blaylock, G., et al. 2006, ApJ, 641, 740
- Schubnell, M. S., Akerlof, C. W., Biller, S., et al. 1996, ApJ, 460, 644
- Takahashi, T., Tashiro, M., Madejski, G., et al. 1996, ApJ, 470, L89+
- Thieler, A. et al. 2010, private communication, DPG conference, Bonn
- Tluczykont, M., Shayduk, M., Kalekin, O., & Bernardini, E. 2007, in Journal of Physics: Conference Series, ed. T. M. F. Halzen, Vol. 60, 318–320, TeV Particle Astrophysics II Workshop
- Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271
- Weekes, T. C. 2008, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 3–17