Free Access
Issue
A&A
Volume 551, March 2013
Article Number A94
Number of page(s) 7
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201220612
Published online 01 March 2013

© ESO, 2013

1. Introduction

The class of very high energy (VHE; E > 100 GeV) γ-ray binaries comprises only a handful of known objects in our Galaxy: LS 5039 (Aharonian et al. 2005a), LS I +61 303 (Albert et al. 2006), PSR B1259-63/LS 2883 (Aharonian et al. 2005b) and HESS J0632+057 (Aharonian et al. 2007), the first binary primarily discovered at TeV energies. This class can be extended by two more candidates: Cygnus X-1 (Albert et al. 2007), a stellar-mass black hole binary detected at VHEs at the 4.1σ significance level, and HESS J1018-589 (HESS Collaboration et al. 2012), which shows a point-like component spatially coincident with the GeV binary 1FGL J1018.6-5856 discovered by the Fermi LAT collaboration (Abdo et al. 2010). Only for the source presented in this paper, PSR B1259-63/LS 2883, the compact companion is clearly identified as a pulsar, making it a unique object for the study of the interaction between pulsar and stellar winds and the emission mechanisms in such systems.

PSR B1259-63/LS 2883 was discovered in a high-frequency radio survey devoted to the detection of young, distant, and short-period pulsars (Johnston et al. 1992a,b). It consists of a rapidly rotating pulsar with a spin period of ≃48 ms and a spin-down luminosity of ≃8 × 1035 erg/s in a highly eccentric (e = 0.87) orbit around a massive Be star. The pulsar moves around the companion with a period Porb = 3.4 years (1237 days).

Latest optical observations with VLT UT2 (Negueruela et al. 2011) significantly improved the previously known parameters of the companion star LS 2883. The luminosity of the star is L = 2.3 × 1038 erg s-1. Because of its fast rotation, the star is oblate with an equatorial radius of Req = 9.7 R and a polar radius of Rpole = 8.1 R. This leads to a strong gradient of the surface temperature from Teq ≈ 27 500 K at the equator to Tpole ≈ 34 000 K at the poles. The mass function of the system suggests a mass of the star M ≈ 30  M and an orbital inclination angle iorb ≈ 25° for the smallest neutron star mass of 1.4  M. The optical observations also suggest that the system is located at the same distance as the star association Cen OB1 at d = (2.3 ± 0.4) kpc (Negueruela et al. 2011). The companion Be star features an equatorial disk that is believed to be inclined with respect to the pulsar’s orbital plane (Johnston et al. 1992a; Melatos et al. 1995; Negueruela et al. 2011) in a way that the pulsar crosses the disk twice in each orbit just before (~20 days) and just after (~20 days) the periastron.

Since its discovery in 1992, PSR B1259-63/LS 2883 is constantly monitored by various instruments at all energy bands. The source shows broadband emission and is visible from radio wavelengths up to the VHE regime. The properties of the radio emission differ depending on the distance between the pulsar and the star. Radio observations (Johnston et al. 1999, 2005; Connors et al. 2002) show that when the pulsar is far from the periastron, the observed radio emission consists only of the pulsed component, whose intensity is almost independent of the orbital position. But closer to the periastron, starting at about tp − 100 d, where tp is the time of periastron, the intensity starts to decrease up to the complete disappearance approximately at tp − 20 d. This is followed by an eclipse of the pulsed emission for about 35−40 days when the pulsar is behind the disk. In contrast, a transient unpulsed component appears and sharply rises to a level more than ten times higher than the flux density of the pulsed emission far from the periastron. The unpulsed component is believed to come from synchrotron radiation generated in the shocked wind zone between the relativistic pulsar wind and the stellar disk outflow. After the disk crossing the unpulsed emission shows a slight decrease with another increase around tp + 20 d at the second crossing of the disk. Radio observations around the 2007 periastron passage showed extended unpulsed emission with a total projected extent of ~120 AU and the peak of the emission clearly displaced from the binary system orbit (Moldón et al. 2011). This indicates that a flow of synchrotron-emitting particles, which can travel far away from the system, can be produced in PSR B1259-63/LS 2883. The source was also monitored around the 2010 periastron passage. The pulsed radio emission was monitored with the Parkes telescope, revealing an eclipse of the pulsed signal that lasted from tp − 16 d to tp + 15 d. Radio emission from PSR B1259-63 at frequencies between 1.1 and 10 GHz was observed using the ATCA array in the period from tp − 31 d to tp + 55 d. The detected unpulsed emission around the periastron passage showed a behaviour similar to the one observed during previous observations (Abdo et al. 2011).

Table 1

Analysis results of the H.E.S.S. data for the full observation period as well as for the pre-flare and flare periods.

PSR B1259-63/LS 2883 is very well covered by X-ray observations carried out with various instruments such as ROSAT (Cominsky et al. 1994), ASCA (Kaspi et al. 1995; Hirayama et al. 1999), INTEGRAL (Shaw et al. 2004), and XMM-Newton (Chernyakova et al. 2006). The periastron passage in 2007 was monitored at the same time by Suzaku, Swift, XMM-Newton, and Chandra (Chernyakova et al. 2009). Observations around the 2010 periastron passage were performed by three instruments: Swift, Suzaku, and XMM-Newton (Abdo et al. 2011). Observations confirmed the 1−10 keV light curve shape obtained in previous periastron observations, showing a rapid X-ray brightening that started at about tp − 25 d with a subsequent decrease closer to periastron and a second increase of the X-ray flux after periastron (Abdo et al. 2011). X-ray observations did not show any X-ray pulsed emission from the pulsar. Unpulsed non-thermal radiation from the source appeared to be variable in flux and spectral index. Similarly to radio measurements, the enhancement of the flux occurs shortly before and shortly after the periastron. Unambiguously, the enhancement of the non-thermal emission results from the interaction of the pulsar wind with the circumstellar disk close to the periastron passage.

PSR B1259-63/LS 2883 was observed by H.E.S.S. around the periastron passages in 2004 (Aharonian et al. 2005b) and 2007 (Aharonian et al. 2009), leading to a firm detection on both occasions. In 2004, PSR B1259-63/LS 2883 was observed mostly after the periastron, in 2007 mostly before it. Therefore, the repetitive behaviour of the source, i.e. the recurrent appearance of the source near periastron at the same orbital phase, with the same flux level and spectral shape of the emission, was not precisely confirmed, since the observations covered different orbital phases. However, the similar dependence of the flux on the separation distance between the pulsar and the star for both periastron passages provides a strong indication for the repetitive behaviour (Kerschhaggl 2011). PSR B1259-63/LS 2883 was not detected in observations performed far from periastron in 2005 and 2006, which comprised 8.9 h and 7.5 h of exposure, respectively.

Observations around the 2004 and 2007 periastron passages showed a variable behaviour of the source flux with time. A combined light curve of those two periastron passages indicates two asymmetrical peaks around periastron with a significant decrease of the flux at the periastron itself. Peaks of the TeV emission roughly coincide with the flux enhancement observed in other wavebands as well as with the eclipse of the pulsed radio emission, which indicates the position of the circumstellar disk. This coincidence suggests that the TeV emission from PSR B1259-63/LS 2883 may be connected to the interaction of the pulsar with the disk.

The paper is organised as follows: in Sect. 2, Fermi-LAT observations at the 2010/2011 periastron passage are reviewed. In Sect. 3, the dataset and analysis techniques are described and analysis results are presented. Results are discussed in Sect. 4 and summarised in Sect. 5.

2. Fermi-LAT detection of a post-periastron HE flare

H.E.S.S. observations around the most recent periastron passage, which took place on 15th of December 2010, were performed as part of an extended multiwavelength (MWL) campaign, which included also radio, optical, X-ray and, for the first time, high-energy (HE; E > 100 MeV) observations. This paper is dedicated to the study of the H.E.S.S. results in the context of the HE observations. The detailed study of the MWL emission from the source will be presented in the joint MWL paper, which is currently in preparation.

Observations of the binary system PSR B1259-63/LS 2883 at HEs were performed using the Large Area Telescope (LAT) on board of Fermi. The data taken around the periastron passage were analysed by two independent working groups (Abdo et al. 2011; Tam et al. 2011), yielding similar results for the flaring period (see below), although there are some discrepancies related to the first detection period close to the periastron passage. Those differences do not affect the conclusions drawn in this paper, however. The source was detected close to periastron with a very low photon flux above 100 MeV of about (1−2) × 10-7 cm-2 s-1. After the initial detection, the flux decreased and the source was below the detection threshold until 14th of January, tp + 30 d, when a spectacular flare was detected, which lasted for about seven weeks with an average flux that was about 10 times higher than the flux detected close to periastron (Abdo et al. 2011). The highest day-averaged flux during the flare almost reached the estimate of the spin-down luminosity of the pulsar, which indicates a close to 100% efficiency of the conversion of the pulsar’s rotational energy into γ-rays. The HE emission around the periastron as a function of time significantly differs from the two-peak light curves observed in other wavebands. The flare is not coincident with the post-periastron peak in radio, X-rays, and VHE γ-rays. It is also much brighter than the GeV emission detected close to the periastron passage (Abdo et al. 2011; Tam et al. 2011).

3. H.E.S.S. observations and analysis

3.1. The H.E.S.S. instrument

H.E.S.S. (High Energy Stereoscopic System) is an array of four 13 m diameter imaging atmospheric Cherenkov telescopes located in the Khomas Highland, Namibia, at an altitude of 1800 m above sea level (Hinton 2004). The telescopes are optimised for the detection of VHE γ-rays in the range from 100 GeV to several tens of TeV by imaging Cherenkov light emitted by charged particles in extensive air showers. The total field of view of H.E.S.S. is 5°. The angular resolution of the system is ≲0.1° and the average energy resolution is about 15%. The H.E.S.S. array is capable of detecting point sources with a flux of 1% of the Crab nebula flux at a significance level of 5σ in 25 h when observing at low zenith angles (Aharonian et al. 2006).

3.2. Data set and analysis techniques

PSR B1259-63/LS 2883 observations were scheduled to cover the post-periastron period from January to March 2011. The source was not visible for H.E.S.S. before and at the periastron passage. The observations resulted in a rather small dataset due to unfavourable weather conditions. The collected data correspond to 6 h of livetime after the standard quality selection procedure (Aharonian et al. 2006). These data were taken in five nights, namely January 9/10, 10/11, 13/14, 14/15 and 15/16. The observations were performed at a relatively high average zenith angle of 48° and with a mean offset angle of 0.55° from the test region centred at αJ2000 = 13h02m48s, δJ2000 = −63°50′09″ (Wang et al. 2004).

The data were analysed using the model analysis1 technique with standard cuts (de Naurois & Rolland 2009). The test region was a priori defined as a circle with radius 0.1° (i.e. θ2 < 0.01°2, where θ is defined as the angular distance between the γ-ray event and the nominal target position), which is the standard size for point-like sources. The reflected region background technique was used for the background subtraction. The analysis results were cross-checked with an alternative analysis chain2 using a standard Hillas reconstruction (Aharonian et al. 2006) method for γ/hadron separation and an independent calibration of the raw data. Both analysis chains yielded consistent results.

thumbnail Fig. 1

Overall differential energy spectrum of the VHE γ-ray emission from PSR B1259-63/LS 2883 for the whole observation period from 9th to 16th of January 2011. The solid line denotes the spectral fit with a simple power law. The green band represents the 1σ confidence interval. Points are derived for the minimum significance of 1.5σ per bin. Points’ error bars represent 1σ errors.

3.3. Energy spectrum

The source was detected at an 11.5σ level (Li & Ma 1983) (see Table 1). A spectral analysis of the detected excess events shows that the differential energy spectrum of photons is consistent with a simple power law dN/dE = N0(E/1  TeV)−Γ with a flux normalisation at 1 TeV of N0 = (1.95 ± 0.32stat ± 0.39syst) × 10-12 TeV-1 cm-2 s-1 and a spectral index Γ = 2.92 ± 0.30stat ± 0.20syst (see Fig. 1 and Table 1) with a fit probability of 0.64. The integral flux above 1 TeV averaged over the whole observation period is F(E > 1  TeV) = (1.01 ± 0.18stat ± 0.20sys) × 10-12 cm-2 s-1.

3.4. Light curves

To check for variability of the source a light curve was produced on a night-by-night basis assuming the photon spectral index obtained in the spectral fit (Fig. 2). The spectral index was fixed at the value obtained in the spectral analysis of the total data because of the low statistics for each individual night. The light curve is consistent with a constant resulting in a mean flux of (0.77 ± 0.13) × 10-12 cm-2 s-1 (horizontal line in Fig. 2) with χ2/NDF = 6.35/4 (corresponds to the probability of 0.17; NDF is the number of degrees of freedom), yielding no evidence for variability in the seven-nights observation period. For each individual night the source is detected at a statistical significance level >3σ except for the last point, whose significance is only 1.5σ.

thumbnail Fig. 2

Integrated photon flux above 1 TeV for individual observation nights. The solid horizontal line indicates the fit of a constant to the distribution. The flare start date is indicated by the dashed vertical line.

For a comparison with the GeV flare (see Sect. 4.2) the whole dataset was divided into two datasets: “pre-flare” (tp + 26 d to tp + 29 d) and “flare” (tp + 30 d to tp + 32 d). These two datasets were analysed independently and revealed similar fluxes and significance levels (see Table 1). A spectral analysis of the two datasets shows that both spectra are consistent with a simple power law, yielding similar values of the spectral index (see Table 1). The two spectral indices are consistent with the one obtained for the total dataset. These results are discussed in Sect. 4.2.

thumbnail Fig. 3

Integrated photon flux above 1 TeV as a function of the time with respect to the periastron passage indicated with the vertical dashed line. The corresponding orbital phases (mean anomaly) are shown on the upper horizontal axis. The data from the 2004 (blue squares) (Aharonian et al. 2005b), 2007 (red triangles) (Aharonian et al. 2009), and 2011 (green circles) observation campaigns are shown. For the 2004 and 2011 data the flux is shown in daily bins while for the 2007 data the flux is shown in monthly bins for clarity.

3.5. Re-analysis of the 2004 data

For the data taken around the 2004 periastron passage the energy spectrum had been measured only up to ~3 TeV (Aharonian et al. 2005b) while for the much smaller dataset of 2011 observations the spectrum was measured up to >4 TeV and for the comparable dataset of 2007 observations up to >10 TeV (Aharonian et al. 2009) using more advanced analysis techniques with a better understanding of weak fluxes. To fully compare observation results around different periastron passages (see below), the 2004 data were re-analysed with the current analysis techniques, the same as used for the analysis of the 2011 data described above. The re-analysis results are compatible with the published ones. The differential energy spectrum measured up to ~10 TeV is consistent with a power law with a spectral index Γ = 2.8 ± 0.1stat ± 0.2syst and a flux normalisation at 1 TeV of N0 = (1.29 ± 0.08stat ± 0.26syst) × 10-12 TeV-1 cm-2 s-1. The new analysis of the 2004 data is therefore compatible with the published results when extrapolated above 3 TeV.

4. Discussion

4.1. Comparison with previous H.E.S.S. observations

In Fig. 3, the integrated photon flux above 1 TeV as a function of time with respect to periastron (indicated by the dashed vertical line) is shown. The light curve compiles the data from all three periastron observation campaigns spanning from 100 days before to 100 days after the periastron. The observed flux from the 2010/2011 observation campaign is compatible with the flux detected in 2004 at the similar orbital phases. Observation periods from 2004 and 2007 were separated in time with respect to the periastron position, i.e., observations in 2004 were performed mainly after and in 2007 mainly before the periastron. Therefore, it was impossible to directly confirm the repetitive behaviour of the source by comparing observations of PSR B1259-63/LS 2883 at the same orbital phases. In this perspective, although the 2011 observations do not exactly overlap with the orbital phases of previous studies, they cover the gap in the 2004 data post-periastron light curve and the integrated flux follows the shape of the light curve, yielding a stronger evidence for the repetitive behaviour of the source.

The spectral shape of the VHE γ-ray emission from PSR B1259-63/LS 2883 around the 2010/2011 periastron passage is similar to what was observed during previous periastron passages (Fig. 4). The photon index of 2.92 ± 0.25stat ± 0.2syst inferred from the 2011 data is well compatible with previous results. The spectrum measured for the 2011 data can be resolved only up to ~4 TeV, which is explained by a very low statistics at higher energies due to a short exposure of the source.

thumbnail Fig. 4

Differential energy spectra of the VHE γ-ray emission from PSR B1259-63/LS 2883 for the data collected around the 2004 (blue squares), 2007 (red triangles), and 2010/2011 (green circles) periastron passages. For the 2004 data the spectrum presented in this paper is shown. The 2007 spectrum is extracted from Aharonian et al. (2009).

4.2. Search for the equivalent “GeV Flare” in the H.E.S.S. data

The absence of the flux enhancement during the GeV flare at radio and X-ray wavebands indicates that the GeV flare may be created by physical processes different from the those responisble for the emission at other wavelengths. The VHE post-periastron data obtained with H.E.S.S. around the 2004 periastron passage do not show any evidence of a flux outburst at orbital phases at which the GeV flare is observed. However, the H.E.S.S. observations around the 2004 periastron passage do not comprise the orbital phase when the GeV flare starts. Moreover, to compare H.E.S.S. 2004 data with the GeV flare observed after the 2010 periastron passage, one has to assume that the GeV flare is a periodic phenomenon, which may not be the case. The H.E.S.S. data taken between 9th and 16th of January in 2011 provide a three-day overlap in time with the GeV flare. Therefore, it is possible to directly study any flux enhancement in the VHE band on the time scale of the HE flare. To improve the sensitivity of the variability search the whole period of the H.E.S.S observations was divided into two time periods of almost equal length: before (“pre-flare”) and during (“flare”) the HE flare (see Sect. 3.4). The pre-flare and flare dataset analysis results are presented in Table 1.

To search for variability, the flux as a function of time was fitted with a constant, which resulted in a mean flux of (0.91 ± 0.18)  ×  10-12 cm-2 s-1 (black horizontal dashed line in Fig. 5 top). The fit has a χ2-to-NDF ratio of 0.73/1, which corresponds to a χ2 probability of 0.39, showing no indication for a flux enhancement. Note that the spectral parameters obtained by an independent fit of each of the two periods have been used here.

thumbnail Fig. 5

(Top) Integrated photon fluxes above 1 TeV for the pre-flare and flare periods (see text) are shown as black filled boxes. The dashed horizontal line shows the best fit with a constant. The HE data points above 0.1 GeV as reported by Abdo et al. (2011) are shown as red filled circles. The flare start date is indicated by the dashed vertical line. The left axis indicates the units for the VHE flux and the right (red) axis denotes the units for the HE flux. (Bottom) The spectral energy distribution of the HE-VHE emission. For the HE emission the overall flare spectrum is shown as reported by Abdo et al. (2011). Marking of the data points is the same as in the top panel. Solid lines denote the fit of the Fermi data only with the power law with exponential cut-off (red) and the fit of the H.E.S.S. data only with the power law (black). The dashed black line denotes the fit of the Fermi (excluding upper limits) and H.E.S.S. data together with the power law.

If one assumes that HE and VHE emission are created according to the same scenario, i.e. the same acceleration and radiation processes and sites, then a flux enhancement of the same magnitude as observed at HEs should be also seen at VHEs. To investigate this hypothesis, the flare coefficient κ is introduced as the ratio of the fluxes during the flare period and the pre-flare period. The ratio of the HE (E > 0.1 GeV) flux averaged over the three-day interval between (tp + 30 d) and (tp + 32 d) to the upper limit on the HE pre-flare emission (see Fig. 5 top) yields a lower limit on the HE emission flare coefficient κHE ≥ 9.2. An upper limit on the VHE flare coefficient can be estimated using the profile likelihood method. The likelihood function is defined as a product of two Gaussian distributions of the pre-flare and flare flux measurements φ1 and φ2 correspondingly, stating that the flare measurement φ2 varies around ˜κ˜φ1\hbox{$\tilde{\kappa} \tilde{\phi_1}$}, where the tilde denotes the true value for a parameter. The profile likelihood λ is then built as a function of κ: λ(κ)=L(φ̂1|φ1,φ2)L(φ̂1,κ̂|φ1,φ2)·\begin{equation} \lambda(\kappa) = \frac{L(\hat{\phi}_{1}, \kappa | \phi_{1}, \phi_{2})}{L(\hat{\phi}_{1}, \hat{\kappa}| \phi_{1}, \phi_{2})}\cdot \end{equation}(1)Where the hat denotes the maximum likelihood estimate for a parameter. The variable −2log λ follows a χ2 distribution with one degree of freedom, which allows one to calculate the 99.7 % confidence level (equivalent of 3σ) upper limit of κ99.7% < 3.5. The obtained upper limit is lower than the observed lower limit on κHE.

The statistical tests presented above give two main results:

  • A flare of similar magnitude as observed in the HE band can befirmly rejected in the VHE band at the same orbital phase.

  • There is no significant difference between the pre-flare and the flare flux in the VHE band.

These two results suggest that the HE flare emission has a different nature than the VHE emission.

This conclusion is also supported by the inconsistency of HE and VHE emission spectra (see Fig. 5 bottom). The joint fit of the Fermi and H.E.S.S. data points with the simple power law (the dashed line on Fig. 5 bottom) results in a fit probability of 0.004 and, hence, fails to explain the combined HE/VHE emission even ignoring the Fermi upper limits, which cannot be taken into account in the fit procedure. Moreover, the Fermi upper limits at 1−100 GeV violate any reasonable model that would be able to explain the HE and VHE emission together. The Fermi spectrum alone is consistent with the the power law with exponential cut-off E2dN/dE = N0(E/0.1GeV)peE/Ecutoff with the index p = 0.16 ± 0.32, the cut-off energy of Ecutoff = 0.5 ± 0.2 GeV and the normalisation N0 = (4 ± 0.4) × 10-10 erg cm-2 s-1. The fit probability is 0.27.

Several models have been proposed to explain the VHE emission from the source. In a hadronic scenario, the VHE γ-ray emission could be produced by the interaction of the ultrarelativistic pulsar wind particles with the dense equatorial disk outflow with subsequent production of π0 pions and hence VHE γ-rays (Kawachi et al. 2004; Neronov & Chernyakova 2007). However, the detection of the source before the expected disk passage in 2007 casts doubts on the hadronic scenario, suggesting that the VHE emission should be created at least partly by leptonic processes (Aharonian et al. 2009). Within the leptonic scenario, VHE emission from PSR B1259-63/LS 2883 is explained by the inverse Compton (IC) scattering of shock-accelerated electrons on stellar photons (Tavani & Arons 1997; Kirk et al. 1999; Dubus 2006; Khangulyan et al. 2007).

Few possible explanations for the nature of the HE flare are discussed in the literature. One of them, suggested by Khangulyan et al. (2012), is based on IC scattering of the unshocked pulsar wind on the stellar and circumstellar disk photons. While the pulsar is inside the disk, the IC scattering of the unshocked pulsar wind is supressed due to the high ram pressure. But immediately after the pulsar escapes the disk in the post-periastron phase, the unshocked pulsar wind zone towards the observer increases significantly, while the density of the circumstellar disk photons is still high enough for efficient IC scattering. Therefore, the enhancement of the HE flux is observed. This is not expected in the pre-periastron phase because the termination shock should expand towards the direction opposite to the observer. This model also predicts the difference between the HE and VHE emission, the latter expected to result from the upscattering of the stellar photons by the electrons accelerated at the termination shock between pulsar and stellar winds. Another explanation of the HE flare can be the Doppler boosting of the radiation created by the shocked pulsar wind (Bogovalov et al. 2008; Dubus et al. 2010; Kong et al. 2012). It is unclear though why the flare is not detected at other wavebands, since the Doppler boosting should also enhance X-ray and VHE γ-ray fluxes. This question, however, can be resolved assuming a specific anisotropy of the pulsar wind and the difference of the emission behaviour in different regions of the termination shock, the isotropic emission in the apex and the beamed emission in the tail of the shock (Kong et al. 2012). In this particular case the HE flare is explained by the Doppler boosting of the synchrotron emission. Abdo et al. (2011) suggested that the flare can also be explained by an anisotropy of the pulsar wind and/or stellar material. The anisotropy of electrons with the highest energies would cause an anisotropy of the synchrotron radiation at high energies. In this interpretation, the HE emission is produced by the synchrotron mechanism. The local increase of the stellar wind density would increase the Bremsstrahlung component, which may also cause the HE flare. Regardless of, which mechanism is responsible for the HE flare, the fact that it is observed only after the periastron indicates either a strong dependency of the HE emission on the geometry of the system, i.e., its configuration with respect to the direction to the observer, or some local perturbation of the stellar material. A detailed interpretation of the HE-VHE emission is beyond the scope of this paper and will be discussed in the joint MWL paper.

5. Summary

The binary system PSR B1259-63/LS 2883 was monitored by H.E.S.S. around the periastron passage on 15th of December 2010. The observed flux and spectral shape agree well with what was measured during previous periastron passages. The observations were performed at similar orbital phases as around the 2004 periastron passage, strengthening the evidence for the repetitive behaviour of the source at VHEs.

H.E.S.S. observations were part of a joint MWL campaign that also included radio, optical, X-ray, and HE observations. A spectacular flare observed at HEs with Fermi LAT overlapped in time with the H.E.S.S. observations. A careful statistical study showed that the HE flare does not have a counterpart at VHEs, indicating that the HE and VHE emissions are produced in different physical scenarios.


1

Paris Analysis software version 0-8-18.

2

H.E.S.S. Analysis Package (HAP) version 11-02-pl04.

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 French Ministry for Research, the CNRS-IN2P3 and the Astroparticle Interdisciplinary Programme of the CNRS, the UK Particle Physics and Astronomy Research Council (PPARC), the IPNP of the Charles University, 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.

References

  1. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJS, 188, 405 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011, ApJ, 736, L11 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abramowski, A., et al. (H.E.S.S. Collaboration) 2012, A&A, 541, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Aharonian, F., Akhperjanian, A. G., Aye, K.-M., et al. 2005a, Science, 309, 746 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  5. Aharonian, F., Akhperjanian, A. G., Aye, K.-M., et al. 2005b, A&A, 442, 1 [Google Scholar]
  6. Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 457, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Aharonian, F. A., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2007, A&A, 469, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Aharonian, F., Akhperjanian, A. G., Anton, G., et al. 2009, A&A, 507, 389 [Google Scholar]
  9. Albert, J., Aliu, E., Anderhub, H., et al. 2006, Science, 312, 1771 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  10. Albert, J., Aliu, E., Anderhub, H., et al. 2007, ApJ, 665, L51 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bogovalov, S. V., Khangulyan, D. V., Koldoba, A. V., Ustyugova, G. V., & Aharonian, F. A. 2008, MNRAS, 387, 63 [NASA ADS] [CrossRef] [Google Scholar]
  12. Chernyakova, M., Neronov, A., Lutovinov, A., Rodriguez, J., & Johnston, S. 2006, MNRAS, 367, 1201 [NASA ADS] [CrossRef] [Google Scholar]
  13. Chernyakova, M., Neronov, A., Aharonian, F., Uchiyama, Y., & Takahashi, T. 2009, MNRAS, 397, 2123 [NASA ADS] [CrossRef] [Google Scholar]
  14. Cominsky, L., Roberts, M., & Johnston, S. 1994, ApJ, 427, 978 [NASA ADS] [CrossRef] [Google Scholar]
  15. Connors, T. W., Johnston, S., Manchester, R. N., & McConnell, D. 2002, MNRAS, 336, 1201 [NASA ADS] [CrossRef] [Google Scholar]
  16. de Naurois, M., & Rolland, L. 2009, Astropart. Phys., 32, 231 [NASA ADS] [CrossRef] [Google Scholar]
  17. Dubus, G. 2006, A&A, 456, 801 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Dubus, G., Cerutti, B., & Henri, G. 2010, A&A, 516, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Hinton, J. A. 2004, New A Rev., 48, 331 [Google Scholar]
  20. Hirayama, M., Cominsky, L. R., Kaspi, V. M., et al. 1999, ApJ, 521, 718 [NASA ADS] [CrossRef] [Google Scholar]
  21. Johnston, S., Lyne, A. G., Manchester, R. N., et al. 1992a, MNRAS, 255, 401 [NASA ADS] [CrossRef] [Google Scholar]
  22. Johnston, S., Manchester, R. N., Lyne, A. G., et al. 1992b, ApJ, 387, L37 [Google Scholar]
  23. Johnston, S., Manchester, R. N., McConnell, D., & Campbell-Wilson, D. 1999, MNRAS, 302, 277 [NASA ADS] [CrossRef] [Google Scholar]
  24. Johnston, S., Ball, L., Wang, N., & Manchester, R. N. 2005, MNRAS, 358, 1069 [NASA ADS] [CrossRef] [Google Scholar]
  25. Kaspi, V. M., Tavani, M., Nagase, F., et al. 1995, ApJ, 453, 424 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kawachi, A., Naito, T., Patterson, J. R., et al. 2004, ApJ, 607, 949 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kerschhaggl, M. 2011, A&A, 525, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Khangulyan, D., Hnatic, S., Aharonian, F., & Bogovalov, S. 2007, MNRAS, 380, 320 [NASA ADS] [CrossRef] [Google Scholar]
  29. Khangulyan, D., Aharonian, F. A., Bogovalov, S. V., & Ribó, M. 2012, ApJ, 752, L17 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kirk, J. G., Ball, L., & Skjaeraasen, O. 1999, Astropart. Phys., 10, 31 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kong, S. W., Cheng, K. S., & Huang, Y. F. 2012, ApJ, 753, 127 [NASA ADS] [CrossRef] [Google Scholar]
  32. Li, T.-P., & Ma, Y.-Q. 1983, ApJ, 272, 317 [NASA ADS] [CrossRef] [Google Scholar]
  33. Melatos, A., Johnston, S., & Melrose, D. B. 1995, MNRAS, 275, 381 [NASA ADS] [Google Scholar]
  34. Moldón, J., Johnston, S., Ribó, M., Paredes, J. M., & Deller, A. T. 2011, ApJ, 732, L10 [NASA ADS] [CrossRef] [Google Scholar]
  35. Negueruela, I., Ribó, M., Herrero, A., et al. 2011, ApJ, 732, L11 [NASA ADS] [CrossRef] [Google Scholar]
  36. Neronov, A., & Chernyakova, M. 2007, Ap&SS, 309, 253 [NASA ADS] [CrossRef] [Google Scholar]
  37. Shaw, S. E., Chernyakova, M., Rodriguez, J., et al. 2004, A&A, 426, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Tam, P. H. T., Huang, R. H. H., Takata, J., et al. 2011, ApJ, 736, L10 [NASA ADS] [CrossRef] [Google Scholar]
  39. Tavani, M., & Arons, J. 1997, ApJ, 477, 439 [NASA ADS] [CrossRef] [Google Scholar]
  40. Wang, N., Johnston, S., & Manchester, R. N. 2004, MNRAS, 351, 599 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Analysis results of the H.E.S.S. data for the full observation period as well as for the pre-flare and flare periods.

All Figures

thumbnail Fig. 1

Overall differential energy spectrum of the VHE γ-ray emission from PSR B1259-63/LS 2883 for the whole observation period from 9th to 16th of January 2011. The solid line denotes the spectral fit with a simple power law. The green band represents the 1σ confidence interval. Points are derived for the minimum significance of 1.5σ per bin. Points’ error bars represent 1σ errors.

In the text
thumbnail Fig. 2

Integrated photon flux above 1 TeV for individual observation nights. The solid horizontal line indicates the fit of a constant to the distribution. The flare start date is indicated by the dashed vertical line.

In the text
thumbnail Fig. 3

Integrated photon flux above 1 TeV as a function of the time with respect to the periastron passage indicated with the vertical dashed line. The corresponding orbital phases (mean anomaly) are shown on the upper horizontal axis. The data from the 2004 (blue squares) (Aharonian et al. 2005b), 2007 (red triangles) (Aharonian et al. 2009), and 2011 (green circles) observation campaigns are shown. For the 2004 and 2011 data the flux is shown in daily bins while for the 2007 data the flux is shown in monthly bins for clarity.

In the text
thumbnail Fig. 4

Differential energy spectra of the VHE γ-ray emission from PSR B1259-63/LS 2883 for the data collected around the 2004 (blue squares), 2007 (red triangles), and 2010/2011 (green circles) periastron passages. For the 2004 data the spectrum presented in this paper is shown. The 2007 spectrum is extracted from Aharonian et al. (2009).

In the text
thumbnail Fig. 5

(Top) Integrated photon fluxes above 1 TeV for the pre-flare and flare periods (see text) are shown as black filled boxes. The dashed horizontal line shows the best fit with a constant. The HE data points above 0.1 GeV as reported by Abdo et al. (2011) are shown as red filled circles. The flare start date is indicated by the dashed vertical line. The left axis indicates the units for the VHE flux and the right (red) axis denotes the units for the HE flux. (Bottom) The spectral energy distribution of the HE-VHE emission. For the HE emission the overall flare spectrum is shown as reported by Abdo et al. (2011). Marking of the data points is the same as in the top panel. Solid lines denote the fit of the Fermi data only with the power law with exponential cut-off (red) and the fit of the H.E.S.S. data only with the power law (black). The dashed black line denotes the fit of the Fermi (excluding upper limits) and H.E.S.S. data together with the power law.

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.