Free Access
Issue
A&A
Volume 554, June 2013
Article Number A107
Number of page(s) 7
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201321135
Published online 12 June 2013

© ESO, 2013

1. Introduction

Blazars are a composite class of active galactic nuclei (AGN), consisting of BL Lacertae-type objects (BL Lacs) and flat-spectrum radio quasars (FSRQs). They are differentiated by the presence (FSRQs) or the absence (BL Lacs) of strong emission lines in their spectra. The broadband spectra of blazars are dominated by nonthermal emission, characterized by rapid variability (see e.g., Wagner & Witzel 1995) in all frequency regimes, with high and variable polarization in the radio and optical frequency regimes (Aller et al. 2003; Mead et al. 1990). More than three dozen blazars have been detected in VHE γ-rays, the overwhelming majority of which belong to the BL Lac class.

PKS 1510−089 is a FSRQ at a redshift of z = 0.361 (Burbidge & Kinman 1966), with highly polarized radio and optical emission (Stockman et al. 1984). At the milliarcsecond scale individual components in the radio jet show apparent superluminal motion (Homan et al. 2001) as high as 46c (Jorstad et al. 2005) indicating a small inclination angle to the line of sight and high bulk Lorentz factors. Very long baseline interferometry (VLBI) observations of this highly polarized radio jet shows large misalignment between the milliarcsecond and the arcsecond-scale components (Homan et al. 2002). This can be explained by the high bulk Lorentz factor in the jet, which makes a small jet bending appear much larger in the projected orientation seen by an observer.

The broadband spectrum of this source has a synchrotron component that peaks between millimeter and IR wavelengths. Malkan & Moore (1986) report broad emission lines in the spectrum of PKS 1510−089 (confirmed by Tadhunter et al. 1993), as well as a clear UV excess (the “blue bump") on top of the nonthermal continuum. The blue bump is attributed to thermal emission from the accretion disk. The high energy component in the spectrum extends from soft X-rays to GeV γ-rays. PKS 1510−089 has been extensively monitored in X-rays (e.g., Jorstad et al. 2006) and is known to be a bright γ-ray emitter from the EGRET era (Hartman et al. 1999). Kataoka et al. (2008) have shown that the quasi-simultaneous broadband spectral energy distribution of PKS 1510−089 can be well described by an external Compton model with seed photons from a dusty torus. In the soft X-ray band, Suzaku data suggest a hardening in the spectrum, which Kataoka et al. (2008) propose could be due to either a small contribution from the synchrotron self-Compton component, or from bulk-Compton scattered radiation. Fermi-LAT measurements of the high energy (HE, defined as 100 MeV < E < 100 GeV) spectrum of PKS 1510−089 in different flux states during 2008–2009 were presented in Abdo et al. (2010a). The average HE spectrum (derived from the entire data set presented therein) is well described by a log-parabola model, dN/dE ∝ (E/E0)− α − β  ln(E/E0), with the following best-fit values for the three free parameters – the spectral-shape parameters α = 2.23 ± 0.02 and β = 0.09 ± 0.01, and an integral photon flux above 100 MeV of (1.12 ± 0.03) × 10-6 cm-2 s-1, with the reference energy, E0, fixed at 300 MeV. In the multiwavelength data presented by Abdo et al. (2010a) no correlation between the flux variations in the HE and X-ray bands is found. The authors report a positive correlation between the HE and the optical band. They find evidence for a 13 day lag between the HE and optical R-band light curves (with the HE light curve leading). They argue that this behavior can be used to rule out a change in beaming as the main driver for variability in the source. Abdo et al. (2010a) make an estimate for the mass of the black hole in this source of 5.4 × 108M, using a model for the accretion disk temperature profile and the measured UV flux. They deduce a maximum isotropic γ-ray luminosity of ≃ 2 × 1048 erg s-1. Because both these estimates are nearly an order of magnitude smaller than the values they obtain for more distant FSRQs, such as 3C 454.3 or PKS 1502+106, they argue that PKS 1510−089 could be an atypical FSRQ. At the same time, the high ratio between γ-ray and synchrotron luminosities is a typical feature of an FSRQ.

The luminous optical–UV photon fields (broad line emission and the blue bump) in FSRQs can cause substantial absorption of VHE photons by electron-positron pair production (see, e.g., Donea & Protheroe 2003; Liu & Bai 2006; Poutanen & Stern 2010). If the VHE emitting region were to be embedded within the broad-line region (BLR), immersed in the reprocessed accretion-disk emission, the VHE γ-rays may not escape the system. On the other hand, VHE emission has been reported from two FSRQs, viz. 3C 279 (Albert et al. 2008) and PKS 1222+216 (Aleksić et al. 2011).

PKS 1510−089 is a good candidate for a VHE emitting FSRQ because it is a bright GeV blazar with a jet that shows highly relativistic bulk motion (hence beaming effects should be strong). PKS 1510−089 was observed with High Energy Stereoscopic System (H.E.S.S.) in March–April 2009 to search for VHE emission, when it was reported to be flaring in the HE and optical domains. MWL data from this period have been published elsewhere, e.g., Marscher et al. (2010) and Abdo et al. (2010a). The H.E.S.S. data are presented followed by the relevant optical and Fermi-LAT data, used for triggering H.E.S.S. observations and in the discussion section.

2. Observations and analysis results

2.1. H.E.S.S.

Following reports in March 2009 of flaring activity in PKS 1510−089 in the HE domain (D’Ammando et al. 2009; Pucella et al. 2009; Vercellone et al. 2009) as well as in the optical frequencies recorded with the ATOM telescope (an optical flaring state was also independently reported by Villata et al. 2009 and Larionov et al. 2009), it was observed with the H.E.S.S. array (Hinton 2004; Aharonian et al. 2006b) between MJD 54 910 and MJD 54 923. Subsequent observations were performed between MJD 54 948 and MJD 54 950 following more optical flaring and another HE flare that was reported by Cutini et al. (2009). A total of 15.8 h (corrected for dead time) of data passing quality cuts (Aharonian et al. 2006b) were obtained, with zenith angles between 14° and 41°, from all observations.

The data were analyzed using the Model Analysis (de Naurois & Rolland 2009). The relatively high redshift of the source implies that the VHE flux should be strongly attenuated due to the extragalactic background light (EBL) in the optical to IR regime (Nikishov 1962). The amount of EBL extinction increases with the energy of γ-ray photons, resulting in a steepening of the spectrum measured in the VHE band (see, e.g., Salamon & Stecker 1998). Therefore, loose cuts (de Naurois & Rolland 2009), that result in a lower energy threshold, are used in the analysis. A total of 248 γ-ray candidates were recorded from the source direction (Table 1, first row), which corresponds to a firm detection with a statistical significance of 9.2σ (following the method of Li & Ma 1983). The distribution of the squared angular distance of events around the position of PKS 1510−089, determined using the “reflected background” method (Berge et al. 2007), shows a clear excess in the source region as compared to the background region (Fig. 1). The observed excess is consistent with point-like γ-ray emission from PKS 1510−089. The best-fit position of the VHE γ-ray excess is at RA=15h12m52.s2±1.s8stat±1.s3sys\hbox{$\mathrm{RA}=15^\mathrm{h}12^\mathrm{m}52\fs2\pm1\fs8_{\mathrm{stat}}\pm1\fs3_{\mathrm{sys}}$} (J2000), Dec=9621.6±26.5stat±20.0sys\hbox{$\mathrm{Dec}=-9^\circ6\arcmin21\farcs6\pm26\farcs5_{\mathrm{stat}}\pm20\farcs0_{\mathrm{sys}}$} (J2000). This is compatible with the optical position (Andrei et al. 2009) of PKS 1510−089, at a separation of 33.3±22.3\hbox{$33\farcs3\pm22\farcs3$}, within the statistical and pointing errors of H.E.S.S.

Table 1

H.E.S.S. data and analysis results.

thumbnail Fig. 1

Distribution of the squared angular distance of γ-ray candidate events around the position of PKS 1510−089. The angular distribution (in terms of θ2, the square of the angular distance between the source position and the reconstructed arrival direction) of events around the position of PKS 1510−089 is shown. On-source (signal+background) events are shown as hatched histogram, whereas the off-source (background) events are shown as points with error bars. The on-source region is defined by a \hbox{$\theta^2\leqslant0.0125$} deg2, shown by the vertical dotted line. The dashed curve is the PSF modeled as a Gaussian, and a constant scaled to the average background level.

thumbnail Fig. 2

Multiwavelength light curves of PKS 1510−089 for the period between MJD 54 910 to MJD 54 951 in terms of integral fluxes in the respective bands along with the 1σ error bars. Panel a) shows the VHE light curve (one day bins) from H.E.S.S. The horizontal line is obtained by fitting a constant-source model to the data. Panel b) displays the HE light curve derived from Fermi-LAT data. The black open circles show the integrated fluxes in one day bins, whereas a finer binning of 4 h is shown in magenta. The solid line is the average level, and the dashed line is the threshold level used for deriving a flare spectrum during the high state around MJD 54 917. In panel c) the R-band optical fluxes measured with ATOM are shown.

The VHE light curve (in 1 day bins) is shown in the top panel of Fig. 2. The best-fit integral flux (>0.15 TeV), obtained by fitting a constant-flux model to the light curve, is (8.0 ± 1.6) × 10-12 cm-2 s-1 (χ2 = 20.4 for 11 degrees of freedom). The χ2 test for variability, i.e. for a null hypothesis that the flux is constant, yields a p - value ≈ 0.11 for this value of the χ2 test statistic. This is therefore insufficient evidence to claim statistically significant variability (at a 99% confidence level). Because a few bins in the light curve have low statistics (fewer than 10 counts in the source region and/or background region), the appropriate p - value was derived from simulations of a large number (~ 105) of light curves. An average flux of 8.0 × 10-12 cm-2 s-1 was assumed, and the actual exposure times and instrument effective area in each bin were taken into account. All bins could be thus considered in the test for variability.

thumbnail Fig. 3

VHE spectrum of PKS 1510−089 measured with the H.E.S.S. instrument, during March 2009. The solid line is the best-fit power law obtained using the forward-folding method (see text for details). The butterfly is the 68% confidence band, and the points with errorbars (1σ statistical errors) are the energy flux. Arrows denote the 99% C.L. upper limits.

To get the best spectral reconstruction, an additional quality criterion was applied, only those data that were taken with the full array of 4 telescopes were accepted (see Table 1, second row). A total of 6.7 h of good quality 4-telescope data were taken during all observations, yielding 159 γ-ray candidates from the source direction and a statistical significance of 8.2σ. The energy spectrum is derived using a forward-folding technique (Piron et al. 2001). The analysis threshold, Ethr ≈ 0.15 TeV, is given by the energy at which the effective area falls to 10% of its maximum value. For these observations the maximum of the measured differential rate is also at this energy. The likelihood maximization for a power-law hypothesis, dN/dE = N0(E/E0)− Γ, in the energy range 0.15 TeV–1.0 TeV, yields a spectral index of Γ = 5.4  ±  0.7stat ± 0.3sys and a normalization constant of N0 = (1.1  ±  0.2stat ± 0.2sys) × 10-10  cm-2 s-1 TeV-1 at the decorrelation energy, E0 = 0.18  TeV (equivalent to a χ2 of 10.3 with 7 degrees of freedom). The spectral slope is steep compared to other VHE detected blazars, for example, cf. Γ = 4.11 ± 0.68stat for 3C 279 (Albert et al. 2008). The reconstructed H.E.S.S. spectrum is shown in Fig. 3. The integral flux, I(0.15 TeV < E < 1.0 TeV) = (1.0 ± 0.2stat ± 0.2sys) × 10-11  cm-2 s-1 corresponds to ≈3% of the Crab Nebula flux in the same band. Fitting a broken-power-law or a power law with an exponential cutoff does not give a better fit. A spectrum obtained from the entire H.E.S.S. data (first row in Table 1) is compatible with the spectrum derived from the high-quality 4-telescope data, within the statistical errors. The results have been cross-checked using a different analysis method based on multivariate analysis technique (see Ohm et al. 2009, and the references therein), and an independent calibration procedure (Aharonian et al. 2006b).

2.2. Fermi-LAT

It is desirable to compare the VHE spectrum to the Fermi-LAT (Atwood et al. 2009) spectrum during the contemporaneous period. While the Fermi-LAT team has published a spectrum for this source from this period, the integration period used (~1 month) is much larger than the shortest flaring timescale seen by Fermi (Abdo et al. 2010a). The intention here is to obtain a spectrum for the highest state within the Fermi-LAT flares (as defined in Abdo et al. 2010a) that is contemporaneous with HESS observations. Thus the Fermi-LAT data were analyzed using the publicly available Fermi Science Tools1 (v9r23p1-fssc-20111006) and the P7SOURCE_V6 instrument response functions. The light curve over a contemporaneous period as the H.E.S.S. observations is produced by a binned likelihood analysis retaining photons (the source class events) with energies between 200 MeV and 100 GeV from a region of interest of a radius of 10° around the position of PKS 1510−089. All sources from the Fermi-LAT two-years point source catalog (2FGL; Nolan et al. 2012) within an angular distance of 15° of PKS 1510−089 were modeled simultaneously. Pass 7 models2 of the Galactic and extragalactic backgrounds were used. For these two diffuse backgrounds, the normalizations are treated as free parameters in the likelihood analysis. As mentioned before the Fermi-LAT spectrum of this source is best described by a curved log parabola model (Abdo et al. 2010a). Therefore, a log parabola model is used for the source in the likelihood analysis, using the gtlike tool, with the normalization and spectral parameters left free. The light curve derived between MJD 54 909.5 and MJD 54 951.5 is shown in Fig. 2, middle panel. The average (daily binned) integral flux between 200 MeV–100 GeV is (1.26 ± 0.03) × 10-6 cm-2 s-1. Two flares are evident, one centered around MJD 54 916, and the second centered around MJD 54 948. Whereas the H.E.S.S. observations during the first HE flare had good coverage and data quality, the H.E.S.S. observations of the second HE flare commenced only after the HE flare had peaked. A finer binning of 4 h was used around MJD 54 916 to precisely trace the development of the flare.

To characterize the range of variations in the HE band, spectra were extracted over two different epochs from Fermi-LAT data. A spectrum for the long-term average state was derived from the first two years of Fermi-LAT data and a second spectrum was extracted for the high state centered around MJD 54 916. The region of interest, the sources modeled together in the binned likelihood analysis, and the parametrization of the source and the background were identical to that used for the light curve generation. The long-term average spectrum was derived from the data taken between MJD 54 682–MJD 55 412 (i.e., between 04.08.2008–04.08.2010) for energies above 200 MeV using a binned likelihood analysis. This results in a total test statistic (TS, e.g., see Mattox et al. 1996) of 29 440.8. The likelihood analysis was applied in an iterative way, such that the parameters of the 2FGL sources within 15° of PKS 1510−089, and the spectral-shape parameters of PKS 1510−089 are fixed in successive steps, leaving only the normalizations of the source and the two diffuse background components free to vary in the last step. The best-fit values for the log-parabola model (dN/dE = n0   [E/E0− α − β  ln [E/E0]) are – a normalization of n0 = (1.41 ± 0.02) × 10-9  cm-2 s-1 MeV-1, and slope parameters of α = 2.21 ± 0.03 and β = 0.083 ± 0.011. The choice of the reference energy, E0, does not affect the spectral shape and was thus fixed at 260  MeV. These parameters are consistent with the average spectrum presented in Abdo et al. (2010a). A spectrum in logarithm bins in energy (using a log-parabola source model) is also derived for the same period. The average and the energy-binned spectrum are shown in Fig. 4, left panel. A Fermi-LAT spectrum that is strictly simultaneous to the H.E.S.S. spectrum (that is derived from a total exposure of 6.7 h, with gaps in the observations, as well as varying live times for individual exposures) cannot be constructed, because an integration period of at least a few days of Fermi observations is required to derive a meaningful spectrum. Moreover, the H.E.S.S. data used for spectral analysis were spread over a period of 8 days, MJD 54 910–MJD 54 918, during which large flux variations were seen at the Fermi-LAT energies. Because nearly half of the H.E.S.S. exposure was during the brightest phase of the GeV-flare centered around MJD 54 916, the Fermi-LAT data during this period were used to derive a spectrum representing the brightest HE flux state, contemporaneous to the H.E.S.S. observations. Data taken between MJD 54 914.8 to MJD 54 917.5 were used, where the flux points in the Fermi-LAT light curve (in 4-h bins) were above a threshold flux of 2.0 × 10-6 cm-2 s-1 (see Fig. 2). This threshold flux was chosen to ensure adequate statistics to construct spectral bins reaching at least 10 GeV. The analysis of this data-set yielded a total TS = 2354.0, with the best-fit values for a log-parabola model given by a normalization of n0 = (10.5 ± 0.7) × 10-9  cm-2 s-1 MeV-1, and slope parameters of α = 1.81 ± 0.13 and β = 0.161 ± 0.048 (reference energy fixed at 260  MeV). The corresponding energy-binned spectrum is shown in Fig. 4 (left panel). This is roughly an order of magnitude brighter than the average flux, though it should be noted that the spectrum is biased towards the low-energy bins because of the lack of statistics at higher energies.

thumbnail Fig. 4

HE and VHE spectra of PKS 1510−089 measured with Fermi-LAT and H.E.S.S., and upper limits on the EBL induced opacity. Left: the filled butterfly shows the H.E.S.S. spectrum as a 1σ confidence band; the blue and red butterflies shows the spectrum corrected for EBL absorption, using the low model for the EBL from Gilmore et al. (2009) and the upper-limit EBL model from Aharonian et al. (2006a) respectively. The triangles show the long-term average Fermi-LAT spectrum, with error bars denoting 1σ statistical errors; the solid curve is the log-parabola model obtained from a binned likelihood analysis of these data (see text for details). The open squares show the Fermi-LAT spectrum derived from the integration period between MJD 54 914.8 and MJD 54 917.5 (error bars illustrate 1σ statistical errors). Arrows denote upper limits on the flux for those bins where insufficient statistics prohibits a flux measurement. The dashed curve is obtained by re-normalizing the long-term average Fermi-LAT spectrum to the level of the high-flux state. Right: upper limits on the EBL opacity derived from spectral measurements of PKS 1510−089 compared to published models. The solid line with arrows at both ends shows the maximum optical depth, corresponding to upper limits on the EBL. For comparison the opacity corresponding to the previously published H.E.S.S. upper-limit EBL model Aharonian et al. (2006a), the Franceschini et al. (2008) model and the low model from Gilmore et al. (2009) are shown.

2.3. ATOM

The Automatic Telescope for Optical Monitoring (ATOM), Hauser et al. (2004), is a 1  m class optical telescope, located at the H.E.S.S. site. PKS 1510−089 is regularly observed with ATOM. From the onset of the H.E.S.S. observation campaigns, the frequency of optical observations was increased to four observations in the R-band (~640 nm) per night. The R-band light curve, with 500 s exposure per point, not corrected for galactic extinction, is shown in Fig. 2. The optical light curve shows clear and pronounced variability (a reduced χ2 = 55.3, with 117 degrees of freedom from fitting a constant-source model, and the highest flux deviating more than 20σ away from the period average). A prominent flare centered around MJD 54 917 was observed, with another flare occurring after MJD 54 925 that could not be followed up with H.E.S.S. observations because it happened around a full-Moon phase when H.E.S.S. does not operate. A third brightening was seen around MJD 54 946.

These observations show PKS 1510−089 in a relatively high state compared to its average optical brightness, measured with ATOM, between May 2007 and August 2009. For example, in the R-band, compared to the long-term average flux of 1.862 ± 0.004 mJy, the average during the H.E.S.S. observations was 2.597 ± 0.009 mJy and a maximum measured flux of 5.724 ± 0.150 mJy.

3. Discussion

The VHE γ-ray flux of distant (z ≳ 0.1) blazars suffer significant extinction due to pair-production on the diffuse UV-IR photon field in the intergalactic medium (Nikishov 1962; Gould & Schréder 1966). This phenomenon can be used to derive upper limits on the photon density of the UV-IR part of the EBL (e.g., Aharonian et al. 2006a; Mazin & Raue 2007).

Here upper limits are derived by comparing the VHE spectra measured with H.E.S.S. to an extrapolation of the contemporaneous HE spectral measurements obtained with Fermi-LAT, following the procedures in, e.g., Georganopoulos et al. (2010), Aleksić et al. (2011) and Meyer et al. (2012). Modifying Eq. (2) from Aleksić et al. (2011), to include the effect of systematic errors, the 95% confidence upper limit on the EBL, expressed in terms of the optical depth, τmax, is given by τmax(E)=ln[Fint(E)(1ε)×(Fobs(E)1.64×ΔFobs(E))],\begin{equation} \label{Eq1} \tau_{\mathrm{max}}(E) = \ln \left[\frac{F_{\mathrm{int}}(E)}{(1-\varepsilon)\times(F_{\mathrm{obs}}(E) - 1.64 \times \Delta F_{\mathrm{obs}}(E))} \right], \end{equation}(1)where Fobs(E) is the measured flux, ΔFobs(E) is the corresponding 1σ statistical error, Fint(E) is the assumed unattenuated intrinsic flux and ε is the systematic error expressed as a fraction of the measured flux. For a steep spectrum this systematic error (i.e. the uncertainty in the absolute normalization of the flux) is ~25% (i.e. ε = 0.25). The systematic error is conservatively considered as a net bias in the flux measurements. It is the factor by which Fobs(E) could overestimate the true flux. The estimation of τmax in Eq. (1) therefore involves a correction factor of (1 − ε). The extrapolation of the HE spectrum measured with Fermi-LAT into the VHE regime is considered as an estimate for the unattenuated spectrum. The large variations in the HE light curve during the H.E.S.S. observing period require the need of contemporaneous spectral measurements. As already mentioned, a strictly simultaneous Fermi-LAT spectrum corresponding only to the short exposure of H.E.S.S. (Row 2 of Table 1) cannot be obtained because of low statistics in the Fermi-LAT data. Even the spectrum obtained during the high flux state (integrated over 2.7d, centered around MJD 54 916) does not provide sufficient statistics at the high-energy end. Therefore, a conservative estimate resulting in the highest possible upper limit on the opacity is derived by re-normalizing the long-term average Fermi-LAT spectrum to the level of the high-flux state (dashed curve in Fig.4, left panel). This is justified because the long-term spectrum gives the best available description of the HE spectral shape up to the highest GeV energies, and changes in spectral shape during flares are small3. This scaled-up spectrum is extrapolated to the VHE regime to estimate the unattenuated flux, Fint(E), used in Eq. (1). It is assumed that the spectrum smoothly extends from the HE to the VHE band, without spectral breaks or cutoffs. Any intrinsic spectral breaks resulting in lower fluxes at higher energies would yield a lower opacity from EBL attenuation. These assumptions along with the above mentioned considerations (the inclusion of systematic errors on the H.E.S.S. flux and the selection of the intrinsic spectral shape consistent with the brightest HE flux state) give firm upper limits to the EBL opacity. Other sources of opacity, e.g., absorption due to internal photon fields from the BLR or the dusty torus (Wagner et al. 1995), would decrease the contribution from EBL extinction to the overall opacity and thus the previous statement still holds. The EBL upper limits derived using Eq. (1) are shown in Fig. 4, right panel. The EBL limits from this distant FSRQ and the limits in Aharonian et al. (2006a), that are derived from BL Lac type objects that are relatively nearby, are at a comparable level.

Owing to their higher bolometric luminosities FSRQs are detected up to larger distances than BL Lacs. Furthermore, due to the bright emission lines in the spectrum of these objects, which the BL Lacs lack by definition, the redshifts of FSRQs can be accurately measured. Thus the potential of putting strong constraints on the EBL extinction is very promising with FSRQs. However, in the case of PKS 1510−089, other factors currently offset this potential. The H.E.S.S. spectrum turns out to be steep as expected for a high redshift source. This results in a slightly higher systematic error of 25% on the VHE flux measurements (cf. ~20% for a Crab-nebula-like spectrum, index ~2.4, Aharonian et al. 2006b). Furthermore, this sharply falling spectrum makes it difficult to collect good statistics at energies >0.5 TeV and hence the statistical errors in this energy range are also high. The combined effect of the large statistical and systematic errors in the VHE spectrum (denoted by ΔFobs and ε respectively in Eq. (1)), because of the steep VHE spectral slope of this source, limits the effectiveness of the method used above. Thus only weak EBL limits could be derived. Another aspect that limits the ability to put strong EBL constraints is the uncertainty in the intrinsic spectral shape. For the observations presented here there is insufficient evidence for claiming variability in the VHE band. However, because the unabsorbed HE band is highly variable (time scale of hours) it is difficult to obtain an accurate description of the intrinsic spectral shape (Fint in Eq. (1)) due to the typically multi-day integration time required, hence making it difficult to put stricter EBL constraints. For example, if the estimated Fint at the lowest energies were estimated to be 20% less than the value used, which is perfectly plausible when comparing to the factor of ~4 variations in the Fermi fluxes seen during this period, the EBL UL from this source calculated as above would have been more constraining than the Aharonian et al. (2006a) limits. With more multiwavelength monitoring of bright FSRQs, such as PKS 1510−089 and other FSRQs that have less steep VHE spectrum, it could be possible to obtain a HE and VHE spectrum when the fluxes in neither band varies. This of course is better done with more sensitive Cherenkov telescopes, such as the H.E.S.S. II telescope array, which should provide richer statistics due to their higher sensitivities and larger energy coverage. This should allow putting stronger constraints on the EBL extinction by comparing their quiescent-state HE and VHE spectrum.

The luminous broad line emission in the optical-UV band as well as the thermal UV excess in PKS 1510−089 indicate an intense internal photon field. This can in principle cause substantial absorption of VHE γ-rays within the BLR radius, or due to the reprocessed disk emission from the dusty torus at larger distances from the accretion disk. The H.E.S.S. detection of this source implies a low optical depth in the VHE emitting region. This can be tentatively explained by hypothesizing that the VHE emitting region is in an optically thin part of the jet, presumably far outside the BLR, see e.g., Tanaka et al. (2011) and Tavecchio et al. (2011) for discussion on the FSRQ PKS 1222+216. It should be noted that in the study made in Poutanen & Stern (2010), where a search is made for signature features of BLR absorption in the HE spectrum of several of the brightest FSRQs detected by Fermi, the HE spectrum of PKS 1510−089 is consistent with a negligible amount of absorption due to line emission in BLR clouds. This is consistent with the hypothesis that the γ-ray emitting zone is in an optically thin region. A detailed discussion of the internal opacity is beyond the scope of this work.

4. Summary

VHE emission was detected from the flat-spectrum radio quasar PKS 1510−089 with a statistical significance of 9.2 standard deviations in 15.8 h of H.E.S.S. data taken during March and April 2009. An integral flux, in the energy regime between 0.15–1.0   TeV, of (1.0 ± 0.2stat ± 0.2sys) × 10-11  cm-2 s-1 is measured, which is ≈3% of the Crab Nebula flux. The spectrum is extremely steep with a photon index of Γ = 5.4 ± 0.7stat ± 0.3sys for a power-law hypothesis.

There is insufficient evidence to claim significant variability in the H.E.S.S. data. However, the multifrequency data on PKS 1510−089 during the H.E.S.S. observation period shows clear and pronounced variability in the HE and optical bands. Using both the Fermi-LAT and H.E.S.S. γ-ray spectrum, from 200  MeV to 1  TeV, upper limits on the optical depth due to the EBL are derived. The EBL-limits in this work are comparable with the limits in Aharonian et al. (2006a), that were derived from BL Lac objects that are relatively nearby. FSRQs, such as PKS 1510−089, due to their higher luminosities compared to BL Lacs, can in principle, allow us to probe the EBL density to relatively higher redshifts. With the upcoming H.E.S.S. II telescope array more precise measurements of the spectra of distant blazars are expected, which will allow us to put stronger constraints on the EBL.


3

Abdo et al. (2010a) find a harder when brighter trend in the Fermi-LAT spectrum when the flux above 0.2  GeV is ≳ 2.4 × 10-7 cm-2 s-1. However, the change in the photon index is small, ≲0.2.

Acknowledgments

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 UK 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. We acknowledge the information provided by the Fermi-LAT collaboration on the Fermi-LAT measured fluxes of PKS 1510−089 during the 2009 flaring episodes, which helped scheduling of the H.E.S.S. observations presented here. This work has been partially supported by the International Max Planck Research School (IMPRS) for Astronomy & Cosmic Physics at the University of Heidelberg. We acknowledge the anonymous referee for providing useful suggestions that improved the clarity.

References

  1. Abdo, A. A., Ackermann, M., Ajello, M., et al. (Fermi Collaboration) 2009, ApJ, 700, 597 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abdo, A. A., Ackermann, M., Agudo, I., et al. (Fermi Collaboration) 2010a, ApJ, 721, 1425 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abdo, A. A., Ackermann, M., Ajello, M., et al. (Fermi Collaboration) 2010b, ApJS, 188, 405 [NASA ADS] [CrossRef] [Google Scholar]
  4. Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. (H.E.S.S. Collaboration) 2006a, Nature, 440, 1018 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  5. Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. (H.E.S.S. Collaboration) 2006b, A&A, 457, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Albert, J., Aliu, E., Anderhub, H., et al. (MAGIC Collaboration) 2008, Science, 320, 1752 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  7. Aleksić, J., Antonelli, L. A., Antoranz, P., et al. (MAGIC Collaboration) 2011, ApJ, 730, L8 [NASA ADS] [CrossRef] [Google Scholar]
  8. Aller, M. F., Aller, H. D., & Hughes, P. A. 2003, ApJ, 586, 33 [NASA ADS] [CrossRef] [Google Scholar]
  9. Andrei, A. H., Souchay, J., Zacharias, N., et al. 2009, A&A, 505, 385 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  10. Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 [NASA ADS] [CrossRef] [Google Scholar]
  11. Berge, D., Funk, S., Hinton, J., et al. 2007, A&A, 466, 1219 [Google Scholar]
  12. Burbidge, E. M., & Kinman, T. D. 1966, ApJ, 145, 654 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cutini, S., Hays, E., & the Fermi Large Area Telescope Collaboration 2009, ATel, 2033 [Google Scholar]
  14. D’Ammando, F., Vercellone, S., Tavani, M., et al. 2009, ATel, 1957 [Google Scholar]
  15. de Naurois, M., & Rolland, L. 2009, Astropart. Phys., 32, 231 [NASA ADS] [CrossRef] [Google Scholar]
  16. Donea, A.-C., & Protheroe, R. J. 2003, Astropart. Phys., 18, 377 [Google Scholar]
  17. Franceschini, A., Rodighiero, G., & Vaccari, M. 2008, A&A, 487, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Georganopoulos, M., Finke, J. D., & Reyes, L. C. 2010, ApJ, 714, L157 [NASA ADS] [CrossRef] [Google Scholar]
  19. Gilmore, R., Madau, P., Primack, J. R., et al. 2009, MNRAS, 399, 1694 [Google Scholar]
  20. Gould, R. J., & Schréder, G. P. E. 1966, Phys. Rev. Lett., 16, 252 [Google Scholar]
  21. Hartman, R. C., Bertsch, D. L., Bloom, S. D., et al. 1999, ApJS, 123, 79 [NASA ADS] [CrossRef] [Google Scholar]
  22. Hauser, M., Möllenhoff, C., Pühlhofer, G., et al. 2004, Astron. Nachr., 325, 659 [Google Scholar]
  23. Hinton, J. 2004, New. Astron. Rev., 48, 331 [Google Scholar]
  24. Homan, D. C., Ojha, R., Wardle, J. F. C., et al. 2001, ApJ, 549, 840 [NASA ADS] [CrossRef] [Google Scholar]
  25. Homan, D. C., Wardle, J. F. C., Cheung, C. C., et al. 2002, ApJ, 580, 742 [NASA ADS] [CrossRef] [Google Scholar]
  26. Jorstad, S. G., Marscher, A. P., Lister, M. L., et al. 2005, AJ, 130, 1418 [NASA ADS] [CrossRef] [Google Scholar]
  27. Jorstad, S. G., Marscher, A. P., Aller, M. F., & Balonek, T. J. 2006, AGN Variability from X-Rays to Radio Waves, ASP Conf. Ser., 360, 169 [NASA ADS] [Google Scholar]
  28. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Kataoka, J., Madejski, G., Sikora, M., et al. 2008, ApJ, 672, 787 [NASA ADS] [CrossRef] [Google Scholar]
  30. Larionov, V. M., Villata, M., Raiteri, C. M., et al. 2009, ATel, 1990 [Google Scholar]
  31. Li, T.-P., & Ma, Y.-Q. 1983, ApJ, 272, 317 [NASA ADS] [CrossRef] [Google Scholar]
  32. Liu, H. T., & Bai, J. M. 2006, ApJ, 653, 1089 [NASA ADS] [CrossRef] [Google Scholar]
  33. Malkan, M. A., & Moore, R. L. 1986, ApJ, 300, 216 [NASA ADS] [CrossRef] [Google Scholar]
  34. Mattox, J. R., Bertsh, D. L., Chiang, J., et al. 1996, ApJ, 461, 396 [NASA ADS] [CrossRef] [Google Scholar]
  35. Marscher, A. P., Jorstad, S. G., Larinov, V. M., et al. 2010, ApJ, 710, L126 [NASA ADS] [CrossRef] [Google Scholar]
  36. Mazin, D., & Raue, M. 2007, A&A, 471, 439 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Mead, A. R. G., Ballard, K. R., Brand, P. W. J. L., et al. 1990, A&AS, 83, 183 [NASA ADS] [Google Scholar]
  38. Meyer, M., Raue, M., Mazin, D., et al. 2012, A&A, 542, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Nikishov, A. I. 1962, Sov. Phys. J. Exp. Theor. Phys., 14, 393 [Google Scholar]
  40. Nolan, P. A., Abdo, A. A., Ackermann, M., et al. (Fermi Collaboration) 2012, ApJS, 199, 31 [NASA ADS] [CrossRef] [Google Scholar]
  41. Ohm, S., van Eldik, C., & Egberts, K. 2009, Astropart. Phys., 31, 383 [NASA ADS] [CrossRef] [Google Scholar]
  42. Piron, F., Djannati-Ataï, A., Punch, M., et al. 2001, A&A, 374, 895 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Poutanen, J., & Stern, B. 2010, ApJ, 717, L118 [Google Scholar]
  44. Pucella, G., D’Ammando, F., Tavani, M., et al. 2009, ATel, 1968 [Google Scholar]
  45. Salamon, M. H., & Stecker, F. W. 1998, ApJ, 493, 547 [NASA ADS] [CrossRef] [Google Scholar]
  46. Stockman, H. S., Moore, R. L., & Angel, J. R. P. 1984, ApJ, 279, 485 [NASA ADS] [CrossRef] [Google Scholar]
  47. Tadhunter, C. N., Morganti, R., Di Serego-Alighieri, S., et al. 1993, MNRAS, 263, 999 [NASA ADS] [CrossRef] [Google Scholar]
  48. Tavecchio, F., Becerra-Gonzalez, J., Ghisellini, G., et al. 2011, A&A, 534, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Tanaka, Y. T., Stawarz, Ł., Thompson, D. J., et al. 2011, ApJ, 733, 19 [NASA ADS] [CrossRef] [Google Scholar]
  50. Vercellone, S., D’Ammando, F., Pucella, G., et al. 2009, ATel, 1976 [Google Scholar]
  51. Villata, M., Raiteri, C. M., Larionov, V. M., et al. 2009, ATel, 1988 [Google Scholar]
  52. Wagner, S. J., & Witzel, A. 1995, ARA&A, 33, 163 [NASA ADS] [CrossRef] [Google Scholar]
  53. Wagner, S. J., Camenzind, M., Dreissigacker, O., et al. 1995, A&A, 298, 688 [NASA ADS] [Google Scholar]

All Tables

Table 1

H.E.S.S. data and analysis results.

All Figures

thumbnail Fig. 1

Distribution of the squared angular distance of γ-ray candidate events around the position of PKS 1510−089. The angular distribution (in terms of θ2, the square of the angular distance between the source position and the reconstructed arrival direction) of events around the position of PKS 1510−089 is shown. On-source (signal+background) events are shown as hatched histogram, whereas the off-source (background) events are shown as points with error bars. The on-source region is defined by a \hbox{$\theta^2\leqslant0.0125$} deg2, shown by the vertical dotted line. The dashed curve is the PSF modeled as a Gaussian, and a constant scaled to the average background level.

In the text
thumbnail Fig. 2

Multiwavelength light curves of PKS 1510−089 for the period between MJD 54 910 to MJD 54 951 in terms of integral fluxes in the respective bands along with the 1σ error bars. Panel a) shows the VHE light curve (one day bins) from H.E.S.S. The horizontal line is obtained by fitting a constant-source model to the data. Panel b) displays the HE light curve derived from Fermi-LAT data. The black open circles show the integrated fluxes in one day bins, whereas a finer binning of 4 h is shown in magenta. The solid line is the average level, and the dashed line is the threshold level used for deriving a flare spectrum during the high state around MJD 54 917. In panel c) the R-band optical fluxes measured with ATOM are shown.

In the text
thumbnail Fig. 3

VHE spectrum of PKS 1510−089 measured with the H.E.S.S. instrument, during March 2009. The solid line is the best-fit power law obtained using the forward-folding method (see text for details). The butterfly is the 68% confidence band, and the points with errorbars (1σ statistical errors) are the energy flux. Arrows denote the 99% C.L. upper limits.

In the text
thumbnail Fig. 4

HE and VHE spectra of PKS 1510−089 measured with Fermi-LAT and H.E.S.S., and upper limits on the EBL induced opacity. Left: the filled butterfly shows the H.E.S.S. spectrum as a 1σ confidence band; the blue and red butterflies shows the spectrum corrected for EBL absorption, using the low model for the EBL from Gilmore et al. (2009) and the upper-limit EBL model from Aharonian et al. (2006a) respectively. The triangles show the long-term average Fermi-LAT spectrum, with error bars denoting 1σ statistical errors; the solid curve is the log-parabola model obtained from a binned likelihood analysis of these data (see text for details). The open squares show the Fermi-LAT spectrum derived from the integration period between MJD 54 914.8 and MJD 54 917.5 (error bars illustrate 1σ statistical errors). Arrows denote upper limits on the flux for those bins where insufficient statistics prohibits a flux measurement. The dashed curve is obtained by re-normalizing the long-term average Fermi-LAT spectrum to the level of the high-flux state. Right: upper limits on the EBL opacity derived from spectral measurements of PKS 1510−089 compared to published models. The solid line with arrows at both ends shows the maximum optical depth, corresponding to upper limits on the EBL. For comparison the opacity corresponding to the previously published H.E.S.S. upper-limit EBL model Aharonian et al. (2006a), the Franceschini et al. (2008) model and the low model from Gilmore et al. (2009) are shown.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.