Binary System PSR B1259-63 Research Articles

New insight into the origin of the GeV flare in the binary system PSR B1259-63/LS 2883 from the 2017 periastron passage

ABSTRACT PSR B1259-63 is a gamma-ray binary system hosting a radio pulsar orbiting around an O9.5Ve star, LS 2883, with a period of ∼3.4 yr. The interaction of the pulsar wind with the LS 2883 outflow leads to unpulsed broad-band emission in the radio, X-rays, GeV, and TeV domains. While the radio, X-ray, and TeV light curves show rather similar behaviour, the GeV light curve appears very different with a huge outburst about a month after a periastron. The energy release during this outburst seems to significantly exceed the spin-down luminosity of the pulsar and both the GeV light curve and the energy release vary from one orbit to the next. In this paper, we present for the first time the results of optical observations of the system in 2017, and also reanalyse the available X-ray and GeV data. We present a new model in which the GeV data are explained as a combination of the bremsstrahlung and inverse Compton emission from the unshocked and weakly shocked electrons of the pulsar wind. The X-ray and TeV emission is produced by synchrotron and inverse Compton emission of energetic electrons accelerated on a strong shock arising due to stellar/pulsar winds collision. The brightness of the GeV flare is explained in our model as a beaming effect of the energy released in a cone oriented, during the time of the flare, in the direction of the observer.

The geometric distance and binary orbit of PSR B1259–63

The pulsar/massive star binary system PSR B1259-63 / LS 2883 is one of the best-studied gamma-ray binaries, a class of systems whose bright gamma-ray flaring can provide important insights into high-energy physics. Using the Australian Long Baseline Array we have conducted very long baseline interferometric observations of PSR B1259-63 over 4.4 years, fully sampling the 3.4-year orbital period. From our measured parallax of $0.38\pm0.05$ mas we use a Bayesian approach to infer a distance of $2.6^{+0.4}_{-0.3}$ kpc. We find that the binary orbit is viewed at an angle of $154\pm3$ degrees to the line of sight, implying that the pulsar moves clockwise around its orbit as viewed on the sky. Taking our findings together with previous results from pulsar timing observations, all seven orbital elements for the system are now fully determined. We use our measurement of the inclination angle to constrain the mass of the stellar companion to lie in the range 15-31$M_{\odot}$. Our measured distance and proper motion are consistent with the system having originated in the Cen OB1 association and receiving a modest natal kick, causing it to have moved $\sim$8 pc from its birthplace over the past $\sim3\times10^5$ years. The orientation of the orbit on the plane of the sky matches the direction of motion of the X-ray synchrotron-emitting knot observed by the Chandra X-ray Observatory to be moving away from the system.

On the timing behaviour of PSR B1259–63 under the propeller torque from a transient accretion disc

The $\gamma$-ray pulsar binary system PSR B1259-63 flares in GeV after each periastron. The origin of those flares is still under debate. Recently, [2017ApJ...844..114Y] proposed a mechanism that might explain the GeV flares. In that model, a transient accretion disc is expected to be formed from the matter which was gravity-captured by the neutron star from the main sequence companion's circumstellar disc. The transient accretion disc exerts a spin down torque on the neutron star (propeller effect), which might be traceable via pulsar timing observation of PSR B1259-63. In this paper, we phenomenologically consider the propeller effect with a parameter $\chi$, which describes the coupling between the disc matter and the neutron star. Comparing the expected timing residuals against the recent observation in [2014MNRAS.437.3255S], we conclude that the angular momentum transfer is very weak (with the coupling parameter $\chi\le10^{-4}$).

Gamma–Gamma Absorption in the γ-ray Binary System PSR B1259-63/LS 2883

Abstract The observed TeV light curve from the γ-ray binary PSR B1259-63/LS 2883 shows a decrease in the flux at periastron that has not been fully explained by emission mechanisms alone. This observed decrease can, however, be explained by γγ absorption due to the stellar and disk photons. We calculate the γγ absorption in PSR B1259-63/LS 2883 taking into account photons from both the circumstellar disk and star, assuming that the γ-rays originate at the position of the pulsar. The γγ absorption due to the circumstellar disk photons produces a ≈14% decrease in the flux, and there is a total decrease of ≈52% (>1 TeV) within a few days before periastron, accompanied by a hardening of the γ-ray photon index. While the γγ absorption alone is not sufficient to explain the full complexity of the H.E.S.S. γ-ray light curve, it results in a significant decrease in the predicted flux, which is coincident with the observed decrease. In addition, we have calculated an upper limit on the γγ absorption, assuming that the emission is produced at the apex of the bow shock. Future observations with CTA during the 2021 periastron passage may be able to confine the location of the emission based on the degree of γγ absorption, as well as measure the hardening of the spectrum around periastron.

A SEARCH FOR VERY HIGH ENERGY GAMMA RAYS FROM THE MISSING LINK BINARY PULSAR J1023+0038 WITH VERITAS

ABSTRACT The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259–63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than ∼2 G before the disappearance of the radio pulsar and greater than ∼10 G afterward.

Multiwavelength monitoring and X-ray brightening of Be X-ray binary PSR J2032+4127/MT91 213 on its approach to periastron

The radio and gamma-ray pulsar PSR J2032+4127 was recently found to be in a decades-long orbit with the Be star MT91 213, with the pulsar moving rapidly towards periastron. This binary shares many similar characteristics with the previously unique binary system PSR B1259-63/LS 2883. Here, we describe radio, X-ray, and optical monitoring of PSR J2032+4127/MT91 213. Our extended orbital phase coverage in radio, supplemented with Fermi gamma-ray data, allows us to update and refine the orbital period to 45-50 yr and time of periastron passage to 2017 November. We analyse archival and recent Chandra and Swift observations and show that PSR J2032+4127/MT91 213 is now brighter in X-rays by a factor of ~70 since 2002 and ~20 since 2010. While the pulsar is still far from periastron, this increase in X-rays is possibly due to collisions between pulsar and Be star winds. Optical observations of the Halpha emission line of the Be star suggest that the size of its circumstellar disc may be varying by ~2 over timescales as short as 1-2 months. Multiwavelength monitoring of PSR J2032+4127/MT91 213 will continue through periastron passage, and the system should present an interesting test case and comparison to PSR B1259-63/LS 2883.

Optical spectroscopy of PSR B1259-63/LS 2883 during the 2014 periastron passage with the Southern African Large Telescope

The gamma-ray binary system PSR B1259-63/LS 2883 went through periastron in May 2014. We report on the optical spectroscopic monitoring of the system from 33 d before to 78 d after periastron, undertaken with the Southern African Large Telescope (SALT). The H{\alpha} and He I ({\lambda}6678) lines exhibit an orbital variation around periastron, with the line strengths reaching a maximum ~13 d after periastron. The line strength is weaker than observed around the previous periastron in 2010. There is also a marked change in the line-strength and asymmetry around the first disc crossing. These observations are consistent with the disruption of the circumstellar disc around periastron due to the interaction with the pulsar.

Multi-wavelength observations of the binary system PSR B1259−63/LS 2883 around the 2014 periastron passage

We report on the results of the extensive multi-wavelength campaign from optical to GeV gamma-rays of the 2014 periastron passage of PSR B1259-63, which is a unique high-mass gamma-ray emitting binary system with a young pulsar companion. Observations demonstrate the stable nature of the post-periastron GeV flare and prove the coincidence of the flare with the start of rapid decay of the H$\alpha$ equivalent width, usually interpreted as a disruption of the Be stellar disk. Intensive X-ray observations reveal changes in the X-ray spectral behaviour happening at the moment of the GeV flare. We demonstrate that these changes can be naturally explained as a result of synchrotron cooling of monoenergetic relativistic electrons injected into the system during the GeV flare.

Variability of VHEγ–ray emission from the binary PSR B1259–63/LS 2883

We examine changes of the $\gamma$-ray intensity observed from the direction of the binary system PSR B1259-63/LS 2883 during campaigns around its three periastron passages. A simple and straightforward method is applied to the published data obtained with the Imaging Atmospheric Cherenkov Technique. Regardless of many issues of the detection process, the method works only with numbers of very high energetic photons registered in the specified regions. Within the realm of this scheme, we recognized changes attributable to the variations of the intrinsic source activity at high levels of significance.

HIGH-ENERGY OBSERVATIONS OF PSR B1259–63/LS 2883 THROUGH THE 2014 PERIASTRON PASSAGE: CONNECTING X-RAYS TO THE GeV FLARE

The binary system PSR B1259-63/LS 2883 is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital phase-dependent variability in these frequencies. The first detection of GeV gamma-rays from the system was made around the 2010 periastron passage. In this Letter, we present an analysis of X-ray and gamma-ray data obtained by the Swift/XRT, NuSTAR/FPM, and Fermi/LAT, through the recent periastron passage which occurred on 2014 May 4. While PSR B1259-63/LS 2883 was not detected by the LAT before and during this passage, we show that the GeV flares occurred at a similar orbital phase as in early 2011, thus establishing the repetitive nature of the post-periastron GeV flares. Multiple flares each lasting for a few days have been observed and short-term variability is seen as well. We also found X-ray flux variation contemporaneous with the GeV flare for the first time. A strong evidence of the keV-to-GeV connection came from the broadband high-energy spectra, which we interpret as synchrotron radiation from the shocked pulsar wind.

Multifrequency Behaviour of the Gamma-Ray Binary System PSR B1259-63: Modelling the FERMI Flare

This paper presents a brief overview of the multifrequency properties of the gamma-ray binary system PSR B1259-63 from radio to very high energy gamma-rays. A summary is also presented of the various models put forward to explain the Fermi "flare" detected in 2011. Initial results are presented of a new turbulence driven model to explain the GeV observations.

Pair cascades in the disk environment of the binary system PSR B1259-63/LS 2883

PSR B1259-63/LS 2883 is a very high energy (VHE; E>100 GeV) γ-ray emitting binary consisting of a 48 ms pulsar orbiting around a Be star with a period of ∼3.4 years. The Be star features a circumstellar disk which is inclined with respect to the orbit in such a way that the pulsar crosses it twice every orbit. The circumstellar disk provides an additional field of target photons which may contribute to inverse Compton scattering and γγ-absorption, leaving a characteristic imprint in the observed spectrum of the high energy emission. At GeV energies, the source was detected for the first time during the previous periastron passage which took place on December 15, 2010. The Fermi Large Area Telescope (LAT) reported a spectacular and unexpected γ-ray flare occurring around 30 days after periastron and lasting for about 7 weeks. In this paper, we study the signatures of Compton-supported, VHE γ-ray induced pair cascades in the circumstellar disc of the Be star and their possible contribution to the GeV flux. We show that cascade emission generated in the disk cannot be responsible for the GeV flare, but it might explain the GeV emission observed close to periastron. We also show that the γγ-absorption in the disk might explain the observed TeV light curve.

Multiwavelength observations of the binary system PSR B1259−63/LS 2883 around the 2010–2011 periastron passage

We report on broad multi-wavelength observations of the 2010-2011 periastron passage of the gamma-ray loud binary system PSR B1259-63. High resolution interferometric radio observations establish extended radio emission trailing the position of the pulsar. Observations with the Fermi Gamma-ray Space Telescope reveal GeV gamma-ray flaring activity of the system, reaching the spin-down luminosity of the pulsar, around 30 days after periastron. There are no clear signatures of variability at radio, X-ray and TeV energies at the time of the GeV flare. Variability around periastron in the H$\alpha$ emission line, can be interpreted as the gravitational interaction between the pulsar and the circumstellar disk. The equivalent width of the H$\alpha$ grows from a few days before periastron until a few days later, and decreases again between 18 and 46 days after periastron. In near infrared we observe the similar decrease of the equivalent width of Br$\gamma$ line between the 40th and 117th day after the periastron. For the idealized disk, the variability of the H$\alpha$ line represents the variability of the mass and size of the disk. We discuss possible physical relations between the state of the disk and GeV emission under assumption that GeV flare is directly related to the decrease of the disk size.

H.E.S.S. observations of the binary system PSR B1259-63/LS 2883 around the 2010/2011 periastron passage

Aim. In this paper we present very high energy (VHE; E>100 GeV) data from the \gamma-ray binary system PSR B1259-63/LS 2883 taken around its periastron passage (15th of December 2010) with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes. We aim to search for a possible TeV counterpart of the GeV flare detected by the Fermi LAT. In addition, we aim to study the current periastron passage in the context of previous observations taken at similar orbital phases, testing the repetitive behavior of the source. Methods. Observations at VHE were conducted with H.E.S.S. from 9th to 16th of January 2011. The total dataset amounts to around 6 h of observing time. Results. The source is detected in the 2011 data at a significance level of 11.5\sigma\ revealing an averaged integral flux above 1 TeV of (1.01 \pm 0.18_{stat} \pm 0.20_{sys}) \times 10^{-12} cm^{-2}s^{-1}. The differential energy spectrum follows a power-law shape with a spectral index \Gamma = 2.92 \pm 0.30_{stat} \pm 0.20_{sys} and a flux normalisation at 1 TeV of N_{0} = 1.95 \pm 0.32_{stat} \pm 0.39_{sys}) \times 10^{-12} TeV^{-1} cm^{-2} s^{-1}. The measured lightcurve does not show any evidence for variability of the source on the daily scale. Conclusions. The measured integral flux and the spectral shape of the 2011 data are compatible with the results obtained around previous periastron passages. The absence of variability in the H.E.S.S. data indicates that the GeV flare observed by Fermi LAT in the time period covered also by H.E.S.S. observations originates in a different physical scenario than the TeV emission. Additionaly, new results compared to those obtained in the observations which were performed in 2004 at a similar orbital phase, further support the hypothesis of the repetitive behavior of the source.

Anisotropic inverse Compton scattering of photons from the circumstellar disc in PSR B1259−63

The gamma-ray binary system PSR B1259-63 consists of a 48 ms pulsar orbiting a Be star. The system is particularly interesting because it is the only gamma-ray binary system where the nature of the compact object is known. The non-thermal radiation from the system is powered by the spin-down luminosity of the pulsar and the unpulsed radiation originates from the stand-off shock front which forms between the pulsar and stellar wind. The Be star/optical companion in the system produces an excess infrared flux from the associated circumstellar disc. This infrared excess provides an additional photon source for inverse Compton scattering. We discuss the effects of the IR excess near periastron, for anisotropic inverse Compton scattering and associated gamma-ray production. We determine the infrared excess from the circumstellar disc using a modified version of a curve of growth method, which takes into account the changing optical depth through the circumstellar disc during the orbit. The model is constrained using archive data and additional mid-IR observations obtained with the VLT during January 2011. The inverse Compton scattering rate was calculated for three orientations of the circumstellar disc. The predicted gamma-ray light curves show that the disc contribution is a maximum around periastron and not around the disc crossing epoch. This is a result of the disc being brightest near the stellar surface. Additional spectroscopic and near-infrared observations were obtained of the system and these are discussed in relation to the possibility of shock heating during the disc crossing epoch.

ANISOTROPIC SCATTERING FROM THE CIRCUMSTELLAR DISC IN PSR B1259-63

The gamma-ray binary system PSR B1259-63 has recently passed through periastron and has been of particular interest as it was observed by Fermi near the December 2010 periastron passage. The system has been detected at very high energies with H.E.S.S. The most probable production mechanism is inverse Compton scattering between target photons from the optical companion and disc, and relativistic electrons in the pulsar wind. We present results of a full anisotropic inverse Compton scattering model of the system, taking into account the IR excess from the extended circumstellar disc around the optical companion.

Modelling the interaction between relativistic and non-relativistic winds in the binary system PSR B1259−63/SS2883★- II. Impact of the magnetization and anisotropy of the pulsar wind

In this paper, we present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with its environment in binary system. We compare the impact of both the magnetic field and the wind anisotropy to the benchmark case of a purely hydrodynamical (HD) interaction of isotropic winds, which has been studied in detail by Bogovalov et al. (2008). We consider the interaction in axisymmetric approximation, i.e. the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star and the effects related to the pulsar orbiting are neglected. The impact of the magnetic field is studied for the case of weak magnetization (with magnetization parameter $\sigma<0.1$), which is consistent with conventional models of pulsar winds. The effects related to anisotropy in pulsar winds are modeled assuming that the kinetic energy flux in a non-magnetized pulsar wind is strongly anisotropic, with the minimum at the pulsar rotation axis and the maximum in the perpendicular direction. We show that, although both considered effects change the shape of the region occupied by the terminated pulsar wind, their impact appears to be small. In particular, for the magnetization of the pulsar wind below 0.1, the magnetic field pressure remains well below the plasma pressure in the post-shock region. Thus, in the case of interaction of a pulsar with the stellar wind environment (opposite to the case of plerions, i.e. the pulsar interaction with the interstellar medium, when the magnetic field becomes dynamically important independently on the wind magnetization) the HD approach represents a feasible approximation for numerical modelling.

GAMMA-RAY SIGNAL FROM THE PULSAR WIND IN THE BINARY PULSAR SYSTEM PSR B1259–63/LS 2883

Binary pulsar systems emit potentially detectable components of gamma ray emission due to Comptonization of the optical radiation of the companion star by relativistic electrons of the pulsar wind, both before and after termination of the wind. The recent optical observations of binary pulsar system PSR B1259-63/LS 2883 revealed radiation properties of the companion star which differ significantly from previous measurements. In this paper we study the implications of these observations for the interaction rate of the unshocked pulsar wind with the stellar photons and the related consequences for fluxes of high energy (HE) and very high energy (VHE) gamma rays. We show that the signal should be strong enough to be detected with Fermi close to the periastron passage, unless the pulsar wind is strongly anisotropic or the Lorentz factor of the wind is smaller than $10^3$ or larger that $10^5$. The higher luminosity of the optical star also has two important implications: (i) attenuation of gamma rays due to photon-photon pair production, and (ii) Compton drag of the unshocked wind. While the first effect has an impact on the lightcurve of VHE gamma rays, the second effect may significantly decrease the energy available for particle acceleration after termination of the wind.

DISCOVERY OF HIGH-ENERGY GAMMA-RAY EMISSION FROM THE BINARY SYSTEM PSR B1259–63/LS 2883 AROUND PERIASTRON WITH FERMI

We report on the discovery of \geq 100 MeV {\gamma} rays from the binary system PSR B1259-63/LS 2883 using the Large Area Telescope (LAT) on board Fermi. The system comprises a radio pulsar in orbit around a Be star. We report on LAT observations from near apastron to ~ 60 days after the time of periastron, tp, on 2010 December 15. No {\gamma}-ray emission was detected from this source when it was far from periastron. Faint {\gamma}-ray emission appeared as the pulsar approached periastron. At ~ tp + 30d, the \geq 100 MeV {\gamma}-ray flux increased over a period of a few days to a peak flux 20-30 times that seen during the pre-periastron period, but with a softer spectrum. For the following month, it was seen to be variable on daily time scales, but remained at ~ 1 - 4 \times 10^-6 cm^-2 s^-1 before starting to fade at ~ tp + 57d. The total {\gamma}-ray luminosity observed during this period is comparable to the spin-down power of the pulsar. Simultaneous radio and X-ray observations of the source showed no corresponding dramatic changes in radio and X-ray flux between the pre-periastron and post-periastron flares. We discuss possible explanations for the observed {\gamma}-ray-only flaring of the source.

DISCOVERY OF GeV γ-RAY EMISSION FROM PSR B1259–63/LS 2883

The binary system PSR B1259-63/LS 2883 consists of a 47.8 ms radio pulsar that orbits the companion Be star with a period of 3.4 years in a highly eccentric orbit. The system is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital-phase-dependent variability in all observed frequencies. Here we report on the discovery of >100 MeV gamma-rays from PSR B1259-63/LS 2883 through the 2010 pariastron passage. Using data collected with the Large Area Telescope aboard Fermi from 33 days before pariastron to 75 days after pariastron, PSR B1259-63/LS 2883 is detected at a significance of 13.6 standard deviations. The gamma-ray light curve is highly variable over the above period, with changing photon index that correlates with gamma-ray flux. In particular, two major flares that occur after the pariastron passage were observed. The onset of gamma-ray emission occurs close to, but not at the same orbital phases as, the two disk passages that occur ~1 month before and ~1 month after the pariastron passage. The fact that the GeV orbital light curve is different from that of the X-ray and TeV light curves strongly suggests that GeV gamma-ray emission originates from a different component. We speculate that the observed GeV flares may be resulting from Doppler boosting effects.