Open Access
Volume 639, July 2020
Article Number A42
Number of page(s) 13
Section Extragalactic astronomy
Published online 07 July 2020

© H. Abdalla et al. 2020

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Blazars are active galactic nuclei (AGNs) with an ultra-relativistic jet pointing towards the Earth. The spectral energy distribution (SED) of blazars exhibits two distinct bumps. The low-energy part (from radio to X-ray) is attributed to synchrotron emission while there is still debate on the emission process responsible for the high-energy bump (from X-ray up to TeV). Synchrotron self-Compton (SSC) models reproduce such emission invoking only leptons. The photons are then produced via synchrotron emission and inverse-Compton scattering. Hadronic blazar models, in which the high-energy component of the blazar SED is ascribed to emission by protons in the jet, or by secondary leptons produced in p − γ interactions, have been widely studied (see e.g. Mannheim 1993; Aharonian 2000; Mücke & Protheroe 2001) as an alternative to leptonic models. These latter convey a certain advantage in that they provide a link between photon, cosmic-ray, and neutrino emission from AGNs, and thus open the multi-messenger path to study AGN jets as cosmic-ray accelerators. Interest in hadronic blazar models has recently increased with the first hint (at 3σ level) of an association of an IceCube high-energy neutrino with the flaring γ-ray blazar TXS 0506+056 (IceCube Collaboration 2018).

To distinguish between the different models, accurate and contemporaneous observations over a wide energy range are of utmost importance. These are possible with the Nuclear Spectroscopic Telescope Array (NuSTAR) launched in 2012, which permits more sensitive studies above 10 keV than previous X-ray missions. Its sensitivity in hard X-rays up to 79 keV enables an examination of the high-energy end of the synchrotron emission even in high-frequency peaked BL Lac (HBL) objects. Such emission is produced by electrons with the highest Lorentz factors, which could be responsible for the γ-ray emission above tens of GeV that can be detected by ground-based facilities such as the High Energy Stereoscopic System (H.E.S.S.).

One of the best-suited objects for joint observations is PKS 2155−304 (z = 0.116, Falomo et al. 1993), a well-known southern object classified as an HBL with HEAO-1 observations in X-rays (Schwartz et al. 1979). The source is a bright and variable γ-ray emitter. Variability with a timescale of about one month was reported in the GeV energy range by the Fermi-Large Area Telescope (LAT; Acero et al. 2015) as well as variations on timescales of approximately one day and rapid flaring events (Cutini 2014, 2013). First detected at TeV energies by Chadwick et al. (1999) in 1996 with the Durham Mark 6 atmospheric Cerenkov telescope, PKS 2155−304 has been regularly observed by H.E.S.S. since the beginning of H.E.S.S. operations, allowing detailed studies of the source variability (H.E.S.S. Collaboration 2017a; Chevalier et al. 2019). The TeV flux of the object exhibits log-normal flux variability behaviour across the whole energy range (H.E.S.S. Collaboration 2017a; Chevalier et al. 2019) making its flux level and variability unpredictable with possible huge flaring events in TeV (Aharonian et al. 2007).

An interesting aspect of this object is the fact that several authors (Zhang 2008; Foschini et al. 2008; Madejski et al. 2016) reported possible contamination of the hard X-ray spectra by the high-energy component (referred to as the hard tail hereafter), but unfortunately, no very high-energy (VHE, E >  100 GeV) data were taken at that time to further constrain the VHE γ-ray flux. Only one multi-wavelength campaign has been conducted so far, using X-ray instruments, Fermi-LAT, and H.E.S.S. (Aharonian et al. 2009). The gathered data were equally well reproduced by either a leptonic model such as the SSC model (Aharonian et al. 2009) or a lepto-hadronic model (Cerruti et al. 2012).

PKS 2155−304 was subsequently the target of a multi-wavelength campaign from June to October 2013 by NuSTAR, H.E.S.S., the Neil Gehrels Swift Observatory, and Fermi-LAT. These instruments observed PKS 2155−304 to provide contemporaneous data for the first time in a very broad energy range, extending from ultra-violet up to TeV γ-rays and yielding a more complete coverage in the X-ray and γ-ray ranges than the previous campaign held in 2008 (Aharonian et al. 2009).

This paper presents the gathered multi-wavelength data from the 2013 campaign and an analysis of these data in Sect. 2. In Sect. 3, the variability of the source and the X-ray spectra are discussed. Section 4 presents the modeling of the data, and Sect. 5 summarizes the findings of this campaign.

2. Data analysis

PKS 2155−304 is an important calibration source in X-rays and was observed during a cross-calibration campaign with other X-ray instruments early in the NuSTAR mission (Madsen et al. 2017). The multi-wavelength observations of the source in April 2013 including NuSTAR, XMM-Newton, and Fermi-LAT were reported by Madejski et al. (2016), and those are denoted “epoch 0” in this paper. Observations of PKS 2155−304 were made as part of the “Principal Investigator” phase of the NuSTAR mission. The aim was to have those observations take place in exact coincidence with observations by the γ-ray observatory H.E.S.S. Because of diverse constraints (technical problems, bad weather, etc.), H.E.S.S., NuSTAR, and Swift only observed PKS 2155−304 simultaneously during four epochs, where each epoch corresponds to observations conducted on a given night (2013-07-17, 2013-08-03, 2013-08-08, and 2013-09-28), labelled as epochs 1, 2, 3, and 4. Both H.E.S.S. and Swift observed the blazar for two additional epochs (2013-06-05 and 2013-06-19; labelled 5 and 6). Epoch 6 is presented for sake of completeness since the Swift data were found to be unusable (see Sect. 2.4). NuSTAR and Swift also observed PKS 2155−304 during three extra epochs (labelled 7, 8, and 9): those are also reported here for the sake of completeness. For each epoch, Fermi-LAT data were analysed and the results are reported in Sect. 2.2. Figure 1 presents the overall light curve derived from all the epochs.

thumbnail Fig. 1.

Multiwavelength light curve of PKS 2155−304 in (from top to bottom) TeV, GeV, X-ray, and UV. The red lines illustrate the epochs mentioned in the text.

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

The H.E.S.S. array is located in the Khomas Highland, in Namibia (23°16′18″ S, 16°30′01″ E), at an altitude of 1800 m above sea level. Now in its second phase, H.E.S.S., is an array of five imaging Atmospheric Cherenkov telescopes. Four of the telescopes (CT1-4) have segmented optical reflectors of 12 m in diameter consisting of 382 mirrors (Bernlöhr et al. 2003) and cameras composed of 960 photomultipliers. Together these form the array of the H.E.S.S. phase I. The second phase started in September 2012 with the addition of a 28 m diameter telescope (CT5) with a camera of 2048 photomultipliers in the centre of the array. The system operates either in Stereo mode, requiring the detection of an air shower by at least two telescopes (Funk et al. 2004; Holler et al. 2015), or in Mono mode in which the array triggers on events detected only with CT5.

PKS 2155−304 was observed by the full H.E.S.S. phase II array during the present observational campaign. Table 1 gives the date of each observation and the results of the analysis described in the following sections. To ensure good data quality, each observation of 28 min had to pass standard quality criteria (Aharonian et al. 2006). For two nights (2013-08-03 and 2013-09-28; epochs 2 and 4), these criteria were not met by the four 12 m telescopes. Therefore, only CT5 Mono observations are available for these nights.

Table 1.

H.E.S.S. observations of PKS 2155−304.

Data for each night were analysed independently using the Model analysis (de Naurois & Rolland 2009) adapted for the five-telescope array (“Stereo analysis” hereafter). In this case, Loose cuts (with a threshold of 40 photo-electrons) were used to lower the energy threshold. For the Mono analysis, standard cuts (threshold of 60 photo-electrons) were applied to minimize systematic uncertainties.

The spectra obtained at each epoch were extracted using a forward-folding method described in Piron et al. (2001). For each night, a power-law model was used of the form ϕdec(E/Edec)−Γ, where Edec is the decorrelation energy. Table 1 lists the parameters providing the best fits to the data above an energy threshold Eth. This threshold is defined as the energy where the acceptance is 10% of the maximal acceptance.

For completeness, the spectra averaged over the epochs 1, 3, 5, and 6 (Stereo mode observations) and over epochs 2 and 4 (Mono mode observations) were computed separately. Above 200 GeV, both measurements are compatible with each other, with an integrated flux of (4.86 ± 0.30)×10−6 ph cm−2 s−1 for the Stereo mode observations and (2.59 ± 0.38)×10−6 ph cm−2 s−1 for Mono mode observations. All the H.E.S.S. data were analysed together by combining the Stereo and Mono mode observations (see Holler et al. 2015), allowing us to compute an averaged spectrum (see Table 1). The integrated flux above 200 GeV measured for this combined analysis is (3.12 ± 0.47)×10−12 ph cm−2 s−1 TeV−1. A cross check with a different analysis chain (Parsons & Hinton 2014) was performed and yields similar results.

2.2. Fermi-LAT data analysis and results

The Fermi-LAT is a γ-ray pair conversion telescope (Atwood et al. 2009) that is sensitive to γ-rays above 20 MeV. The bulk of LAT observations are performed in an all-sky survey mode ensuring a coverage of the full sky every 3 h.

Data and software used in this work (Fermitools) are publicly available from the Science Support Center1. Events within 10° around the radio coordinates of PKS 2155−304 (region of interest, ROI) and passing the SOURCE selection (Ackermann et al. 2012) were considered corresponding to event class 128 and event type 3 and a maximum zenith angle of 90°. Further cuts on the energy (100 MeV < E < 500 GeV) were made, which remove the events with poor energy resolution. To ensure a significant detection of PKS 2155−304, time windows of 3 days centred on the campaign nights (Table 1) were considered to extract the spectral parameters. To analyse LAT data, P8R3_SOURCE_V2 instrumental response functions (irfs) were used. In the fitting procedure, FRONT and BACK events (Atwood et al. 2009) were treated separately.

The Galactic and extragalactic background models designed for the PASS 8 irfs denoted gll_iem_v07.fits (Acero et al. 2016) and iso_P8R3_SOURCE_V2_v1.txt were used in the sky model, which also contains all the sources of the fourth general Fermi catalogue (4FGL, The Fermi-LAT Collaboration 2020) within the ROI plus 2° to take into account the large point spread function (PSF) of the instrument especially at low energy.

An unbinned maximum likelihood analysis (Mattox et al. 1996) implemented in the gtlike tool2 was used to find the best-fit spectral parameters of each epoch. Models other than the power-law reported here do not significantly improve the fit quality. Table 2 shows the results of the analysis. We note that for epoch 1 with a test statistic (TS) below 25 (≈5σ), a flux upper limit was derived assuming a spectral index of Γ = 1.753.

Table 2.

Fermi-LAT observations of PKS 2155−304.

All the uncertainties presented in this section are statistical only. The most important source of systematic uncertainties in the LAT results is the uncertainty on the effective area, all other systematic effects are listed on the FSSC website4.

2.3. NuSTAR data analysis and results

The NuSTAR satellite developed in the NASA Explorer program features two multilayer-coated telescopes that focus the reflected X-rays onto pixellated CdZnTe focal plane modules and provide an image of a point source with a half-power diameter of ∼1′ (see Harrison et al. 2013, for more details). The advantage of NuSTAR over other X-ray missions is its broad bandpass, 3–79 keV with a spectral resolution of ∼1 keV.

Table 3 provides the details of individual NuSTAR pointings: this includes the amount of on-source time (after screening for the South Atlantic Anomaly passages and Earth occultation) and mean net (background-subtracted) count rates. After processing the raw data with the NuSTAR Data Analysis Software (NuSTARDAS) package v1.3.1 (with the script nupipeline), the source data were extracted from a region of 45″ radius centred on the centroid of X-ray emission, while the background was extracted from a 1.5′ radius region roughly 5′ southwest of the source location located on the same chip. The choice of these parameters is dictated by the size of the point-spread function of the mirror. However, the derived spectra depend very weakly on the sizes of the extraction regions. The spectra were subsequently binned to have at least 30 total counts per re-binned channel. Spectral channels corresponding nominally to the 3–60 keV energy range were considered, in which the source was robustly detected. The resulting spectral data were fitted with a power law modified by the Galactic absorption with a column density of 1.7 × 1020 atoms cm−2 (Dickey & Lockman 1990) using XSPEC v12.8.2. The standard instrumental response matrices and effective area were derived using the ftool nuproducts. The alternate NH measurement by Kalberla et al. (2005) of 1.4 × 1020 cm−2 was tested, and the best-fit spectral parameters of the source were entirely consistent with results obtained using Dickey & Lockman (1990) values. Data for both NuSTAR detectors were fitted simultaneously, allowing an offset of the normalisation factor for the focal plane module B (FPMB) with respect to module FPMA. Regardless of the adopted models, the normalisation offset was less than 5%. The resulting fit parameters are given in Table 3. More complex models for fitting to the datasets obtained during joint NuSTAR and Swift-XRT pointings were considered, and those are discussed in Sect. 3.2.

Table 3.

Summary of the NuSTAR observations of PKS 2155−304.

The source exhibited significant variability in one of the pointings on August 26 (epoch 8); the NuSTAR X-ray count rate for the FPMA module dropped by almost a factor of two in 25 ks clock time (Fig. 2). This was observed independently by both NuSTAR modules. The other NuSTAR observations showed only modest variability, with a nominal min-to-max amplitude of less than 20% of the mean count rate. Such variability is not uncommon in HBL-type BL Lac objects and has been seen in previous observations of PKS 2155−304 (see, e.g. Zhang 2008). More recently, rapid X-ray variability was seen in PKS 2155−304 when it was simultaneously observed by many X-ray instruments (Madsen et al. 2017). Other HBL-type blazars exhibit similar variability; recent examples are Mkn 421 (Baloković et al. 2016) and Mkn 501 (Furniss et al. 2015).

thumbnail Fig. 2.

Light curve of PKS 2155−304 as seen by the FPMA module of NuSTAR during the observation 60002022012 (epoch 8). The energy range is 3–60 keV and the plotted data are not background subtracted. However, the background rate is always lower than 0.03 counts per second and the background was steady (within 5%) throughout the observation. Each point corresponds to data taken over roughly one orbit, during the time indicated by the red markers.

2.4. Swift-XRT data analysis and results

The details of the Swift X-ray Telescope (XRT, Burrows et al. 2005) observations used here are listed in Table 4. The observations were taken simultaneously (or as close as possible) with the H.E.S.S. and NuSTAR observations. During this campaign, Swift observed the source nine times, but for one of the pointings (corresponding to epoch 6, archive sequence 00030795110), applying standard data quality cuts resulted in no useful source data (the source was outside of the nominal Window Timing – WT – window). Two Swift-XRT observations (sequences 0080280006 and -08) were close in time and were performed during a single NuSTAR observation. Because these observations have consistent fluxes and spectra, they were added together as Swift-XRT data for epoch 7.

Table 4.

Summary of the Swift-XRT observations of PKS 2155−304.

All Swift-XRT observations were carried out using the WT readout mode. The data were processed with the XRTDAS software package (version 3.4.0) developed at Space Science Data Center (SSDC5) and distributed by HEASARC within the HEASoft package (version 6.22.1). Event files were calibrated and cleaned with standard filtering criteria with the xrtpipeline task using the calibration files available in the Swift CALDB (v. 20171113). The average spectrum was extracted from the summed and cleaned event file. Events for the spectral analysis were selected within a circle of 20 pixels (∼46″) radius, which encloses about 80% of the PSF, centred on the source position. The background was extracted from a nearby circular region of 20 pixels radius. The ancillary response files (ARFs) were generated with the xrtmkarf task applying corrections for PSF losses and CCD defects using the cumulative exposure map. The latest response matrices (version 15) available in the Swift CALDB were used. Before the spectral fitting, the 0.4–10 keV source spectra were binned to ensure a minimum of 30 counts per bin. The data extending to the last bin with 30 counts were used, which is typically ∼5 keV.

The spectrum of each Swift-XRT observation was fitted with a simple power law with a Galactic absorption column of 1.7 × 1020 atoms cm−2 using the XSPEC v12.8.2 package. The resulting mean count rates, power law indices, and corresponding 2–10 keV model fluxes are also included in Table 4. No variability was found in individual observations in this energy range.

2.5. Spectral fitting of X-ray data and the search for the hard X-ray “tail”

The results of the individual spectral fits of the Swift-XRT and NuSTAR data are given in Tables 3 and 4, respectively. However, because PKS 2155−304 exhibited complex X-ray spectral structure measured in the joint XMM-Newton plus NuSTAR observation in April 2013 (Madejski et al. 2016), here a joint fit to the lower-energy Swift-XRT and the higher-energy NuSTAR data was performed to investigate the need for such increasingly complex models. Since the source is highly variable, only the strictly simultaneous Swift-XRT and NuSTAR data sets were paired. To account for possible effects associated with variability or imperfect Swift-XRT-to-NuSTAR cross-calibration, the normalisations of the models for the two detectors were allowed to vary, but the difference was in no case greater than 20%, consistent with the findings of Madsen et al. (2017), with the exception of the August 26 observation (epoch 8) where NuSTAR revealed significant variability (see note in Sect. 2.3).

To explore the spectral complexity similar to that seen in April 2013, the following models were considered6: (1) PL: a simple power-law model; and (2) LP: a log-parabola model. The resulting joint spectral fits are given in Table 5.

Table 5.

Joint NuSTAR and Swift-XRT observations of PKS 2155−304.

In four observations (epochs 1, 3, 4 and 9), the model consisting of a simple PL absorbed by the Galactic column fits the data well: no deviation from a simple power-law model is required. However, for epochs 2, 7, and 8, a significant improvement (Δχ2 >  20 for one extra parameter) of the fit quality is found by adopting the LP model. Thus, at these epochs, the spectrum steepens with energy. In conclusion, there are not only spectral index changes from one observation epoch to another, but there is also a significant change of the spectral curvature from one observation to another. Bhatta et al. (2018), using only NuSTAR data, reported results on the same observations and also found a change in the spectral shape for epoch 8 but not for epochs 2 or 7. These latter authors also reported a hardening for epochs 1, 3, and 4, but one which is not significant when compared to a PL fit.

A third model consisting of one log-parabola plus a second hard power law with spectral index ΓHT (LPHT)7 was also tested. The model adds a generally harder high-energy “tail” (HT) to the softer log-parabola component. A notable feature is the absence of such a HT in any of the observations (see Sect. 3.2). Therefore, an upper limit on the 20–40 keV flux was computed assuming ΓHT = 2.

2.6. Swift-UVOT data analysis and results

The Ultraviolet/Optical Telescope (UVOT; Burrows et al. 2005) on board Swift also observed PKS 2155−304 during Swift pointings and measured the UV and optical emission in the bands V (500–600 nm), B (380–500 nm), U (300–400 nm), UVW1 (220–400 nm), UVM2 (200–280 nm) and UVW2 (180–260 nm). The values of Schlafly & Finkbeiner (2011) were used to correct for the Galactic absorption8.

The photon count-to-flux conversion is based on the UVOT calibration (Sect. 11 of Poole et al. 2008). A power-law spectral index ΓUV was derived for each epoch and is reported in Table 6. The results presented in this work do not provide evidence for spectral variability in the UV energy range.

Table 6.

Swift-UVOT observations of PKS 2155−304.

3. Discussion

3.1. Flux state and variability in γ-rays

During the observation campaign, PKS 2155−304 was found in a low flux state, in the H.E.S.S. energy range, ϕ(E >  200GeV) = (11.6 ± 1.3)×10−12 ph cm−2 s−1, a factor of approximately five lower than during the 2008 campaign (ϕ(E >  200GeV) = (57.6 ± 1.8)×10−12 ph cm−2 s−1Aharonian et al. 2009); see Fig. 3. The average flux above 200 GeV measured by H.E.S.S. during 9 years of observations (ϕ(E >  200GeV) = (51.0 ± 4.1)×10−12 ph cm−2 s−1, H.E.S.S. Collaboration 2017a) is also more than four times higher than that reported here (for the entire campaign). We note that even lower flux values have been measured over the last 10 years (see Fig. 1 of H.E.S.S. Collaboration 2017a). The source exhibits a harder spectrum (Γ ≈ 2.8) with respect to the H.E.S.S. phase I measurement (Γ ≈ 3.4, Aharonian et al. 2009; H.E.S.S. Collaboration 2017a). This is consistent with the results of H.E.S.S. Collaboration (2017b) and likely to be due to the lower energy threshold achieved with CT5.

thumbnail Fig. 3.

Spectral energy distribution of PKS 2155−304 for each epoch considered in this work. For epoch 0, the red points are directly extracted from Madejski et al. (2016). In the other plots, the purple points are UVOT data, orange are XRT data, and yellow are the NuSTAR data. In γ-rays, the green points and contours are the Fermi-LAT results and H.E.S.S. results are in blue. The black upper limits refer to the hard-tail component (see text) and are used to constrain the inverse-Compton part of the SSC model (black line). The grey points are the data from the 2008 observation campaign (Aharonian et al. 2009) shown for comparison. Black points are the radio data from Abdo et al. (2010) and Liuzzo et al. (2013). The dashed blue line is the synchrotron emission and the orange line is the IC emission. Both are from the SSC calculation, and the black dashed line is the sum of both.

The Fermi-LAT flux averaged over the nine epochs was lower than the flux measured in the 3FGL, (12.6 ± 0.4)×10−8 ph cm−2 s−1, and lower than in 2008 by a factor of approximately two. Similar results were found by H.E.S.S. Collaboration (2017b) showing that the source was in a low flux state in 2013. With a flux of (8 ± 2)×10−8 ph cm−2 s−1 in the 100 MeV–300 GeV energy range, epoch 0 is not different from the epochs reported here. The 2−10 keV X-ray flux was found to be a factor of between approximately three and four lower than in 2008 (Aharonian et al. 2009); see Fig. 3. Only at two epochs (3 and 4), was the 2−10 keV flux measured by NuSTAR lower than the one measured at epoch 0 (1.1 × 10−11 erg cm−2 s−1) and the fluxes of epochs 1, 2, 7, 8, and 9 were higher. The only noticeable difference is at lower energies with the observed optical flux measured by Swift-UVOT: at epoch 0, the flux was higher than that measured in all the other epochs (see Table 6).

3.2. Broad-band X-ray spectrum

In the energy range from 0.3 to 10 keV, the spectrum is usually assumed to be the high-energy end of the synchrotron emission. Indeed, the measured spectral index of PKS 2155−304 in the X-ray regime is generally in agreement with the value expected for a HBL, for which a power-law spectral index, Γ, is typically steeper than 2 (“soft component” hereafter). Nevertheless a single power law is too simple a representation of the spectrum when measured with sensitive instruments affording a good signal-to-noise ratio. As already pointed out by Perlman et al. (2005), the soft X-ray spectra of HBLs are well represented as gradually steepening functions towards higher energies. In the data presented here, the spectral index measured by Swift-XRT is always harder than the one measured by NuSTAR. A Kolmogorov–Smirnov test was performed on both Swift-XRT and NuSTAR spectral index distributions. This test rejects the hypothesis that they are sampled from the same distribution with a p-value of 3%. This suggests that such steepening takes place for PKS 2155−304.

At the end of the X-ray spectrum (roughly above a few keV), Urry & Mushotzky (1982) observed PKS 2155−304 above an energy of a few keV with the HEAO A1 instrument, and Zhang (2008) reported a hard excess in two XMM-Newton observations (confirmed by Foschini et al. 2008 using the same observations). The XMM-Newton observations fit with a broken power-law showed a spectral hardening of ΔΓ = 0.1−0.3 with a break energy of 3–5 keV. Both works interpreted this as a possible contamination of the synchrotron spectra by inverse-Compton emission.

More recently, and with the increased energy range provided by NuSTAR, Madejski et al. (2016) also measured a hard tail in the X-ray spectrum of PKS 2155−304 (April 2013 observations, epoch 0). Using a broken power-law model, they found a flattening spectrum with a spectral break of ΔΓ >  1 around 10 keV. During that observation, the source was found in a very low flux state (with the 2–10 keV flux of 1.1 × 10−11 erg cm−2 s−1), even lower than the flux reported by Zhang (2008) and Foschini et al. (2008). Jointly fitting the strictly simultaneous XMM-Newton data with the NuSTAR data, a more complete picture emerged, with a log-parabola describing the soft (E <  5 keV) spectrum, and a hard tail which can be described as an additional power law.

Regarding the observations presented in this work, adding an extra hard tail (LPHT model) does not significantly improve the χ2. However, it is important to note that the flux of the object during the April 2013 pointing was relatively low, and the observations were fairly long (about four times longer than any single pointing during the campaign reported here). As noted by Madejski et al. (2016), the hard tail becomes more easily detectable only during low-flux states of the softer, low-energy spectral component.

To detect a possible hard tail in the data set of the present campaign, a simultaneous spectral fit of all data sets was performed. Due to the spectral variability of the soft, low-energy component (Table 3), stacking (or just summing) all spectra simultaneously is inappropriate. Instead, a simultaneous fit of seven individual datasets from Epochs 1, 2, 3, 4, 7, 8, and 9 was considered, allowing the spectral parameters of the soft component (described as a log-parabola) to vary independently. Each epoch was described by a LPHT model (see Sect. 2.5), and with common normalisation of the hard tail for all data sets9. Formally, the fit returns zero flux for the hard-tail component. The 99% confidence upper limit of 1.8 × 10−4 ph keV−1 cm−2 on the normalisation of this component (at χ2 + 2.7) corresponds to a 20–40 keV flux limit of 2.5 × 10−13 erg cm−2 s−1. The normalisation of this hard tail in the data from epoch 0 is 8 × 10−4 ph keV−1 cm−2 (corresponding to a 20–40 keV flux of 12.0 × 10−13 erg cm−2 s−1), or more than four times higher than the upper limit measured during the other epochs. In conclusion, the hard tail is also variable on a timescale of months, but no conclusions on the shorter timescales from the presented NuSTAR data can be drawn.

We note that the source does exhibit a similar flux level in X-rays with respect to the April 2013 data set while in optical, the flux is significantly lower. In an SSC framework, this photon field might be scattered by low-energy electrons to produce hard X-ray photons accounting for the hard tail visible in epoch 0. Nevertheless, when the Fermi measurement is extrapolated down towards the NuSTAR energy range, it always overshoots the X-ray measurement. This could be due to a lack of statistics in the LAT range preventing the detection of spectral curvature such as that reported in the 3FGL catalogue, since only 3 days of data were used in each epoch. The extrapolation of the 3FGL spectrum of PKS 2155−304 does not violate the upper limits derived here on the hard-tail component but cannot reproduce epoch 0.

4. SED modelling

4.1. Leptonic modelling: one zone synchrotron self-Compton

Modelling of blazar SEDs was performed with a one-zone SSC model by Band & Grindlay (1985). The emission zone is considered to be a sphere of radius R filled with a magnetic field B and moving at relativistic speed with a Lorentz factor Γ. In this zone, the emitting particle distribution follows a broken power law:


where N is density of electrons at γ = 1, p1 and p2 are the indices of the electron distribution, and γb is the break energy.

The modelling was performed on the epochs presented in this work (1–5) with UV, X-ray, GeV, and TeV data. Radio data from Abdo et al. (2010) and Liuzzo et al. (2013) were taken from the NED10. The radio emission could originate from another location in the jet, or from the emission zone, and is therefore considered as an upper limit in the model. Historical data taken between 10−2 eV and 1 eV (infra-red range) are found to be quite stable in time with variation of less than a factor of two. Such data were collected using Vizier11 and shown in the SEDs.

For each epoch, a mathematical minimisation (Nelder & Mead 1965) was performed to find the model parameters R, B, N, log(γmin), log(γb), and log(γmax) that best fit the data. The values of p1 and p2 were constrained by the UV and X-ray data, respectively, and were not allowed to vary freely in the fitting procedure. Given the little spectral variability found in UV and GeV, p1 was set to 2.5 and p2 = 2 ⋅ ΓX-ray − 1 (Rybicki & Lightman 1986). The minimisation was performed using a Markov chain Monte Carlo (MCMC) implemented in the emcee python package (Foreman-Mackey et al. 2013). For epochs 1–4, the upper limit on the hard-tail flux (Table 5) is taken into account by forcing the inverse-Compton (IC) component of the model to be below this limit. The resulting parameters are given in Table 7 with their corresponding realisations in Fig. 3.

Table 7.

Model parameters for each epoch.

The model parameters are consistent with previous studies by Kataoka et al. (2000), Foschini et al. (2007), Katarzynski et al. (2008), and Aharonian et al. (2009). As in these previous studies, as well as for other BL Lac objects (e.g. Mrk 421 (Abdo et al. 2011a), Mrk 501 (Abdo et al. 2011b), SHBL J001355.9–185406 (H.E.S.S. Collaboration 2013), etc.), the obtained model is far from equipartition. Even with a very low flux state in the present modelling, particles carry at least ten times more energy density than the magnetic field.

The data from epochs 1–5 are well reproduced by the simple SSC calculation presented here. In contrast to Gaur et al. (2017) for this object or Chen (2017) for Mrk 421, there is no need to invoke a second component to reproduce the SED without over-predicting the radio flux. The main difference is that the hard tail above ≈10 keV seen in the previous observations is not observed in the present data set.

The SSC model was applied to the data of epoch 0 and results are also presented in Table 7. The contemporaneous data are well reproduced. The main difference in the modelling parameters between epoch 0 and the campaign presented in this work lies in the values of γmin. For epoch 0, having log(γmin) = 0 allows a greater inverse-Compton contribution in the X-ray band, making the X-ray tail detectable by NuSTAR. This is also in agreement with the observed decrease in the optical flux in epochs 1–5. Indeed a higher value of γmin decreases the number of electrons emitting in this energy range. We also note that the archival radio data are in disagreement with the modelling of epoch 0, which predict an overly high flux in that energy range. The values obtained for different parameters are not equally well constrained. The shape of the electron distribution (γmin, γbreak and γmax) is quite robust with small errors. Other parameters like the B-field or the size of the emitting region remain poorly known and are indeed different from the model presented in Madejski et al. (2016).

4.2. Emergence of a hadronic component in hard X-rays?

Following the detection of a γ-ray flare from TXS 0506+056 coinciding with a high-energy neutrino (IceCube Collaboration 2018), several authors have independently shown that, while pure hadronic models cannot reproduce the multi-messenger dataset, a scenario in which the photon emission is dominated by an SSC component with a subdominant hadronic component is viable (see, e.g. Ansoldi et al. 2018; Cerruti et al. 2018; Gao et al. 2019; Keivani et al. 2018). The hadronic component emerges in the hard-X-rays as synchrotron radiation by secondary leptons produced via the Bethe–Heitler pair-production channel in this scenario. With this result in mind, it was investigated whether the hardening seen in the NuSTAR data of PKS 2155−304 could be due to subdominant hadronic emission. Starting from the simple SSC model for epoch 0 (see Table 7), a population of relativistic protons was added. It was assumed that pp = pe, 1 (i.e. protons and electrons share the same acceleration mechanism, resulting in the same injection spectral index) and that the maximum proton Lorentz factor γp, max is determined by equating acceleration and cooling timescales. The proton distribution was normalised such that the hadronic component emerges in hard X-rays. For additional details on the hadronic code used see Cerruti et al. (2015). Another change in the SSC part of the model was the increase of the value of log(γmin) to 3.3 in order to avoid overshooting the radio emission.

The key parameter is the power in protons Lp required to provide the observed photon flux, because a very well-known drawback of hadronic blazar models is that they often require proton powers well above the Eddington luminosity LEdd of the super-massive black hole which powers the AGN. For the case of PKS 2155−304, if pp = 2.5, γp, min = 1 and log γp, max = 8.0, then Lp = 5.6 × 1050 erg s−1 is needed, which is around 1000LEdd for a black hole mass of , making this scenario unrealistic. This result is very sensitive to the exact shape of the proton distribution, especially at low Lorentz factors (which cannot be constrained by the data). Moreover, Lp is lower if the proton distribution is harder, or if γp, min >  1. As an example, if pp = 2.0 and γp, min = 1000, then Lp = 6.6 × 1047 erg s−1, which is of the same order of magnitude as LEdd. For this scenario, the hadronic photon emission is shown in Fig. 4, and emerges in X-rays as the emission by Bethe–Heitler pairs, and at VHE as a photo-meson cascade. The model predicts an expected neutrino rate in IceCube of νrate = 0.01 yr−1, which is compatible with the non-detection of PKS 2155−304 by IceCube (computed using the IC effective area12 for a declination of −30°).

thumbnail Fig. 4.

Same as Fig. 3 but for epoch 0 only. The blue and orange dashed lines are the synchrotron and inverse-Compton emission as in Fig. 3. The green line is the emission from Bethe–Heitler pair-production and the red line is that from the photo-meson cascade. The sum of all these components is given by the black dashed line. The black continuous line is the predicted neutrino spectrum.

5. Conclusions

PKS 2155−304 was observed contemporaneously for the first time by Swift, NuSTAR, Fermi-LAT, and H.E.S.S. The source was found in a low flux state in all wavelengths during epochs 1–9. The source flux is lower than during the campaign carried out in 2008.

For each epoch, no hard tail was detected in the X-ray spectra, contrary to what was seen at epoch 0. The computation of an upper limit on the 20–40 keV flux of such a hard tail for each observation and for the full data set shows that this component is variable on a timescale of a few months. For epochs 1–5, the SED is well reproduced by a one-zone SSC model. Such a model fails to reproduce the epoch 0 data due to the required value of the γmin parameter. A low value of γmin is mandatory to reproduce the hardening in X-rays but in return produces an overly high flux in the radio band with respect to the archival measurements.

The emergence of the variable X-ray hard tail cannot be explained by a one-zone SSC model. Several authors have proposed a multi-zone model to tackle this issue, and Gaur et al. (2017) in particular used a spine or layer jet structure. In such a structured jet, synchrotron photons of the slow layer are Comptonised by the electrons of the fast spine to produce the hard X-ray tail. The results presented here would imply that the layer producing the hard tail is variable over a timescale of months. Such a result is in agreement with the model parameters of Gaur et al. (2017). Nevertheless, the variability timescale derived from the model parameters of these latter authors cannot reproduce variability of the source on a timescale of days, as the model was not designed to reproduce such variability.

Here, the possibility of having a lepto-hadronic radiation component was explored. The same parameters as for the SSC model but with log(γmin) = 3.3 were used to reproduce a large part of the SED. The hard tail was successfully reproduced by the hadronic emission. Nevertheless for such a model to be in agreement with the Eddington luminosity of the super-massive black hole, the proton distribution has to be harder (pe = 2.0) than the electron distribution (pe = 2.3) together with γp, min >  1000 and/or have a low-energy cut-off γp, min >  1. In the framework of the lepto-hadronic model, the hard-X-ray emission associated with Bethe–Heitler pair production is independent and is not directly associated with the electron-synchrotron and the SSC components. The detection of the hard-X-ray tail during only one of the NuSTAR observations could therefore be explained by a sudden increase in the hadronic injection. The origin of the hard tail is still uncertain but this feature could help to disentangle different classes of emission models for PKS 2155−304 and blazars in general.


An unbinned analysis is recommended for small time bins


This value has been taken a priori and close to the index found in this work.


Models are corrected for Galactic absorption.


The formula for this LPHT model is ϕ ∝ E−Γ − β ⋅ log(E) + E−ΓHT.


See with a reddening ratio Av/E(B − V) = 3.1 and E(B − V) = 0.022.


In an SSC or lepto-hadronic scenario, one would expect the hard X-ray tail to be the low-energy counterpart of the Fermi spectra. The approach made here with the assumption of a constant normalisation for the tail is more conservative than using the γ-ray spectral results.


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 Helmholtz Association, the Alexander von Humboldt Foundation, the French Ministry of Higher Education, Research and Innovation, the Centre National de la Recherche Scientifique (CNRS/IN2P3 and CNRS/INSU), the Commissariat à l’énergie atomique et aux énergies alternatives (CEA), the U.K. Science and Technology Facilities Council (STFC), the Knut and Alice Wallenberg Foundation, the National Science Centre, Poland grant no. 2016/22/M/ST9/00382, the South African Department of Science and Technology and National Research Foundation, the University of Namibia, the National Commission on Research, Science & Technology of Namibia (NCRST), the Austrian Federal Ministry of Education, Science and Research and the Austrian Science Fund (FWF), the Australian Research Council (ARC), the Japan Society for the Promotion of Science and by the University of Amsterdam. We appreciate the excellent work of the technical support staff in Berlin, Zeuthen, Heidelberg, Palaiseau, Paris, Saclay, Tübingen and in Namibia in the construction and operation of the equipment. This work benefited from services provided by the H.E.S.S. Virtual Organisation, supported by the national resource providers of the EGI Federation. The Fermi-LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase from the following agencies is also gratefully acknowledged: the Istituto Nazionale di Astrofisica in Italy and and the Centre National d’Etudes Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. This work was supported under NASA Contract No. NNG08FD60C and made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. We thank the NuSTAR Operations, Software, and Calibration teams for support with the execution and analysis of these observations. This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center (ASDC, Italy) and the California Institute of Technology (USA). This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research made use of Enrico, a community-developed Python package to simplify Fermi-LAT analysis (Sanchez & Deil 2013). This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France. The original description of the VizieR service was published in A&AS 143, 23. This work has been done thanks to the facilities offered by the Université Savoie Mont Blanc MUST computing center. M. Cerruti has received financial support through the Postdoctoral Junior Leader Fellowship Programme from la Caixa Banking Foundation, grant n. LCF/BQ/LI18/11630012. M. B. gratefully acknowledges financial support from NASA Headquarters under the NASA Earth and Space Science Fellowship Program (grant NNX14AQ07H), and from the Black Hole Initiative at Harvard University, which is funded through a grant from the John Templeton Foundation.


  1. Abdo, A. A., Ackermann, M., Agudo, I., et al. 2010, ApJ, 716, 30 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011a, ApJ, 736, 131 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011b, ApJ, 727, 129 [NASA ADS] [CrossRef] [Google Scholar]
  4. Acero, F., Ackermann, M., Ajello, M., et al. 2015, ApJS, 218, 23 [Google Scholar]
  5. Acero, F., Ackermann, M., Ajello, M., et al. 2016, ApJS, 223, 26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Ackermann, M., Ajello, M., Albert, A., et al. 2012, ApJS, 203, 4 [NASA ADS] [CrossRef] [Google Scholar]
  7. Aharonian, F. A. 2000, New Astron., 5, 377 [NASA ADS] [CrossRef] [Google Scholar]
  8. Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 457, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2007, ApJ, 664, L71 [Google Scholar]
  10. Aharonian, F., Akhperjanian, A. G., Anton, G., et al. 2009, ApJ, 696, L150 [NASA ADS] [CrossRef] [Google Scholar]
  11. Ansoldi, S., Antonelli, L. A., Arcaro, C., et al. 2018, ApJ, 863, L10 [NASA ADS] [CrossRef] [Google Scholar]
  12. Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 [NASA ADS] [CrossRef] [Google Scholar]
  13. Baloković, M., Paneque, D., Madejski, G., et al. 2016, ApJ, 819, 156 [NASA ADS] [CrossRef] [Google Scholar]
  14. Band, D. L., & Grindlay, J. E. 1985, ApJ, 298, 128 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bernlöhr, K., Carrol, O., Cornils, R., et al. 2003, Astropart. Phys., 20, 111 [NASA ADS] [CrossRef] [Google Scholar]
  16. Bhatta, G., Mohorian, M., & Bilinsky, I. 2018, A&A, 619, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165 [NASA ADS] [CrossRef] [Google Scholar]
  18. Cerruti, M., Zech, A., Boisson, C., & Inoue, S. 2015, MNRAS, 448, 910 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cerruti, M., Zech, A., Boisson, C., & Inoue, S. 2012, in 5th International Meeting on High Energy Gamma-Ray Astronomy, AIP Conf. Ser., 1505, 635 [Google Scholar]
  20. Cerruti, M., Zech, A., Boisson, C., et al. 2018, MNRAS, 483, L12 [NASA ADS] [CrossRef] [Google Scholar]
  21. Chadwick, P. M., Lyons, K., McComb, T. J. L., et al. 1999, ApJ, 513, 161 [NASA ADS] [CrossRef] [Google Scholar]
  22. Chen, L. 2017, ApJ, 842, 129 [CrossRef] [Google Scholar]
  23. Chevalier, J., Sanchez, D. A., Serpico, P. D., Lenain, J.-P., & Maurin, G. 2019, MNRAS, 484, 749 [CrossRef] [Google Scholar]
  24. Cutini, S. 2013, ATel, 4755 [Google Scholar]
  25. Cutini, S. 2014, ATel, 6148 [Google Scholar]
  26. de Naurois, M., & Rolland, L. 2009, Astropart. Phys., 32, 231 [NASA ADS] [CrossRef] [Google Scholar]
  27. Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  28. Falomo, R., Pesce, J. E., & Treves, A. 1993, ApJ, 411, L63 [NASA ADS] [CrossRef] [Google Scholar]
  29. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
  30. Foschini, L., Ghisellini, G., Tavecchio, F., et al. 2007, ApJ, 657, L81 [NASA ADS] [CrossRef] [Google Scholar]
  31. Foschini, L., Treves, A., Tavecchio, F., et al. 2008, A&A, 484, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Funk, S., Hermann, G., Hinton, J., et al. 2004, Astropart. Phys., 22, 285 [NASA ADS] [CrossRef] [Google Scholar]
  33. Furniss, A., Noda, K., Boggs, S., et al. 2015, ApJ, 812, 65 [NASA ADS] [CrossRef] [Google Scholar]
  34. Gao, S., Fedynitch, A., Winter, W., & Pohl, M. 2019, Nat. Astron., 3, 88 [NASA ADS] [CrossRef] [Google Scholar]
  35. Gaur, H., Chen, L., Misra, R., et al. 2017, ApJ, 850, 209 [CrossRef] [Google Scholar]
  36. Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJ, 770, 103 [NASA ADS] [CrossRef] [Google Scholar]
  37. H.E.S.S. Collaboration (Abramowski, A., et al.) 2013, A&A, 554, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. H.E.S.S. Collaboration (Abdalla, H., et al.) 2017a, A&A, 598, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. H.E.S.S. Collaboration (Abdalla, H., et al.) 2017b, A&A, 600, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Holler, M., Berge, D., van Eldik, C., et al. 2015, ArXiv e-prints [arXiv:1509.02902] [Google Scholar]
  41. IceCube Collaboration (Aartsen, M. G., et al.) 2018, Science, 361, eaat1378 [NASA ADS] [Google Scholar]
  42. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Kataoka, J., Takahashi, T., Makino, F., et al. 2000, ApJ, 528, 243 [NASA ADS] [CrossRef] [Google Scholar]
  44. Katarzynski, K., Lenain, J. P., Zech, A., Boisson, C., & Sol, H. 2008, MNRAS, 390, 371 [NASA ADS] [CrossRef] [Google Scholar]
  45. Keivani, A., Murase, K., Petropoulou, M., et al. 2018, ApJ, 864, 84 [NASA ADS] [CrossRef] [Google Scholar]
  46. Liuzzo, E., Falomo, R., Treves, A., et al. 2013, AJ, 145, 73 [CrossRef] [Google Scholar]
  47. Madejski, G. M., Nalewajko, K., Madsen, K. K., et al. 2016, ApJ, 831, 142 [NASA ADS] [CrossRef] [Google Scholar]
  48. Madsen, K. K., Beardmore, A. P., Forster, K., et al. 2017, AJ, 153, 2 [NASA ADS] [CrossRef] [Google Scholar]
  49. Mannheim, K. 1993, A&A, 269, 67 [NASA ADS] [Google Scholar]
  50. Mattox, J. R., Bertsch, D. L., Chiang, J., et al. 1996, ApJ, 461, 396 [NASA ADS] [CrossRef] [Google Scholar]
  51. Mücke, A., & Protheroe, R. J. 2001, Astropart. Phys., 15, 121 [NASA ADS] [CrossRef] [Google Scholar]
  52. Nelder, J., & Mead, R. 1965, Comput. J., 7, 308 [Google Scholar]
  53. Parsons, R. D., & Hinton, J. A. 2014, Astropart. Phys., 56, 26 [NASA ADS] [CrossRef] [Google Scholar]
  54. Perlman, E. S., Madejski, G., Georganopoulos, M., et al. 2005, ApJ, 625, 727 [NASA ADS] [CrossRef] [Google Scholar]
  55. Piron, F., Djannati-Atai, A., Punch, M., et al. 2001, A&A, 374, 895 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Poole, T. S., Breeveld, A. A., Page, M. J., et al. 2008, MNRAS, 383, 627 [Google Scholar]
  57. Rybicki, G. B., & Lightman, A. P. 1986, Radiative Processes in Astrophysics (Germany: Wiley-VCH), 400 [Google Scholar]
  58. Sanchez, D. A., & Deil, C. 2013, in Proceedings of the 33rd International Cosmic Ray Conference (ICRC 2013) [Google Scholar]
  59. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  60. Schwartz, D. A., Griffiths, R. E., Schwarz, J., Doxsey, R. E., & Johnston, M. D. 1979, ApJ, 229, L53 [NASA ADS] [CrossRef] [Google Scholar]
  61. The Fermi-LAT Collaboration 2020, ApJS, 247, 33 [NASA ADS] [CrossRef] [Google Scholar]
  62. Urry, C. M., & Mushotzky, R. F. 1982, ApJ, 253, 38 [NASA ADS] [CrossRef] [Google Scholar]
  63. Zhang, Y. H. 2008, ApJ, 682, 789 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

H.E.S.S. observations of PKS 2155−304.

Table 2.

Fermi-LAT observations of PKS 2155−304.

Table 3.

Summary of the NuSTAR observations of PKS 2155−304.

Table 4.

Summary of the Swift-XRT observations of PKS 2155−304.

Table 5.

Joint NuSTAR and Swift-XRT observations of PKS 2155−304.

Table 6.

Swift-UVOT observations of PKS 2155−304.

Table 7.

Model parameters for each epoch.

All Figures

thumbnail Fig. 1.

Multiwavelength light curve of PKS 2155−304 in (from top to bottom) TeV, GeV, X-ray, and UV. The red lines illustrate the epochs mentioned in the text.

In the text
thumbnail Fig. 2.

Light curve of PKS 2155−304 as seen by the FPMA module of NuSTAR during the observation 60002022012 (epoch 8). The energy range is 3–60 keV and the plotted data are not background subtracted. However, the background rate is always lower than 0.03 counts per second and the background was steady (within 5%) throughout the observation. Each point corresponds to data taken over roughly one orbit, during the time indicated by the red markers.

In the text
thumbnail Fig. 3.

Spectral energy distribution of PKS 2155−304 for each epoch considered in this work. For epoch 0, the red points are directly extracted from Madejski et al. (2016). In the other plots, the purple points are UVOT data, orange are XRT data, and yellow are the NuSTAR data. In γ-rays, the green points and contours are the Fermi-LAT results and H.E.S.S. results are in blue. The black upper limits refer to the hard-tail component (see text) and are used to constrain the inverse-Compton part of the SSC model (black line). The grey points are the data from the 2008 observation campaign (Aharonian et al. 2009) shown for comparison. Black points are the radio data from Abdo et al. (2010) and Liuzzo et al. (2013). The dashed blue line is the synchrotron emission and the orange line is the IC emission. Both are from the SSC calculation, and the black dashed line is the sum of both.

In the text
thumbnail Fig. 4.

Same as Fig. 3 but for epoch 0 only. The blue and orange dashed lines are the synchrotron and inverse-Compton emission as in Fig. 3. The green line is the emission from Bethe–Heitler pair-production and the red line is that from the photo-meson cascade. The sum of all these components is given by the black dashed line. The black continuous line is the predicted neutrino spectrum.

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.