Modulation of the high energy γ-ray flux from PSR B1259-63/SS2883 due to the orbital variation of the maximum energy of accelerated electrons

The inverse Compton (IC) scattering of ultrarelativistic electrons accelerated at the pulsar wind termination shock is generally believed to be responsible for TeV gamma-ray signal recently reported from the binary system PSR B1259-63/SS2883. In such a system the acceleration takes place in the presence of a dense radiation field provided by a companion Be2-type star. Thus it is natural to expect an orbital phase dependence of the acceleration efficiency in the system. The HESS collaboration reported the tendency of reduction of TeV γ-rays around the periastron. In this paper we study a possible explanation of this effect by the “early” (sub-TeV) cutoffs in the energy spectrum of accelerated electrons due to the enhanced rate of Compton losses close to the periastron.     

The cross-correlation between galaxies and groups: probing the galaxy distribution in and around dark matter haloes

Some massive binaries should contain energetic pulsars which inject relativistic leptons from their inner magnetospheres and/or pulsar wind regions. If the binary system is compact enough, then these leptons can initiate inverse Compton (IC) e^± pair cascades in the anisotropic radiation field of a massive star. γ-rays can be produced in the IC cascade during its development in a pulsar wind region and above a shock in a massive star wind region where the propagation of leptons is determined by the structure of a magnetic field around the massive star. For a binary system with specific parameters, we calculate phase-dependent spectra and fluxes of γ-rays escaping as a function of the inclination angle of the system and for different assumptions on injection conditions of the primary leptons (their initial spectra and location of the shock inside the binary). We conclude that the features of γ-ray emission from such massive binaries containing energetic pulsars should allow us to obtain important information on the acceleration of particles by the pulsars, and on interactions of a compact object with the massive star wind. Predicted γ-ray light curves and spectra in the GeV and TeV energy ranges from such binary systems within our Galaxy and Magellanic Clouds should be observed by future AGILE and GLAST satellites and low-threshold Cherenkov telescopes, such as MAGIC, HESS, VERITAS or CANGAROO III.     

Modelling interaction of relativistic and non-relativistic winds in binary system PSR B1259−63/SS2883 – I. Hydrodynamical limit

In this paper, we present a detailed hydrodynamical study of the properties of the flow produced by the collision of a pulsar wind with the surrounding in a binary system. This work is the first attempt to simulate interaction of the ultrarelativistic flow (pulsar wind) with the non-relativistic stellar wind. Obtained results show that the wind collision could result in the formation of an 'unclosed' (at spatial scales comparable to the binary system size) pulsar wind termination shock even when the stellar wind ram pressure exceeds significantly the pulsar wind kinetic pressure. Moreover, the post-shock flow propagates in a rather narrow region, with very high bulk Lorentz factor (γ∼ 100). This flow acceleration is related to adiabatic losses which are purely hydrodynamical effects. Interestingly, in this particular case, no magnetic field is required for formation of the ultrarelativistic bulk outflow. The obtained results provide a new interpretation for the orbital variability of radio, X-ray and gamma-ray signals detected from binary pulsar system PSR B1259−63/SS2883.     

X-ray observations of PSR B1259−63 near the 2007 periastron passage

PSR B1259−63 is a 48-ms radio pulsar in a highly eccentric 3.4-yr orbit with a Be star SS 2883. Unpulsed γ-ray, X-ray and radio emission components are observed from the binary system. It is likely that the collision of the pulsar wind with the anisotropic wind of the Be star plays a crucial role in the generation of the observed non-thermal emission. The 2007 periastron passage was observed in unprecedented details with Suzaku , Swift , XMM–Newton and Chandra missions. We present here the results of this campaign and compare them with previous observations. With these data we are able, for the first time, to study the details of the spectral evolution of the source over a 2-month period of the passage of the pulsar close to the Be star. New data confirm the pre-periastron spectral hardening, with the photon index reaching a value smaller than 1.5, observed during a local flux minimum. If the observed X-ray emission is due to the inverse Compton (IC) losses of the 10-MeV electrons, then such a hard spectrum can be a result of Coulomb losses, or can be related to the existence of the low-energy cut-off in the electron spectrum. Alternatively, if the X-ray emission is a synchrotron emission of very high-energy electrons, the observed hard spectrum can be explained if the high-energy electrons are cooled by IC emission in Klein–Nishina regime. Unfortunately, the lack of simultaneous data in the TeV energy band prevents us from making a definite conclusion on the nature of the observed spectral hardening and, therefore, on the origin of the X-ray emission.     

Radio-to-TeV γ-ray emission from PSR B1259–63

We discuss the implications of the recent X-ray and TeV γ-ray observations of the PSR B1259–63 system (a young rotation powered pulsar orbiting a Be star) for the theoretical models of interaction of pulsar and stellar winds. We show that previously considered models have problems to account for the observed behaviour of the system. We develop a model in which the broad band emission from the binary system is produced in result of collisions of GeV–TeV energy protons accelerated by the pulsar wind and interacting with the stellar disk. In this model the high energy γ-rays are produced in the decays of secondary neutral pions, while radio and X-ray emission are synchrotron and inverse Compton emission produced by low-energy (≤100 MeV) electrons from the decays of secondary charged π ^± mesons. This model can explain not only the observed energy spectra, but also the correlations between TeV, X-ray and radio emission components.      

Probing pulsar winds using inverse Compton scattering

We investigate the effects of inverse Compton scattering by electrons and positrons in the unshocked winds of rotationally-powered binary pulsars. This process can scatter low energy target photons to produce gamma rays with energies from MeV to TeV. The binary radio pulsars PSR B1259−63 and PSR J0045−73 are both in close eccentric orbits around bright main sequence stars which provide a huge density of low energy target photons. The inverse Compton scattering process transfers momentum from the pulsar wind to the scattered photons, and therefore provides a drag which tends to decelerate the pulsar wind. We present detailed calculations of the dynamics of a pulsar wind which is undergoing inverse Compton scattering, showing that the deceleration of the wind of PSR B1259−63 due to ‘inverse Compton drag' is small, but that this process may confine the wind of PSR J0045−73 before it attains pressure balance with the outflow of its companion star. We calculate the spectra and light curves of the resulting inverse Compton emission from PSR B1259−63 and show that if the size of the pulsar wind nebula is comparable to the binary separation, then the γ-ray emission from the unshocked wind may be detectable by atmospheric Cherenkov detectors or by the new generation of satellite-borne γ-ray detectors such as INTEGRAL and GLAST. This mechanism may therefore provide a direct probe of the freely-expanding regions of pulsar winds, previously thought to be invisible.     

Non-pulsed Emission from the Binary System B1259-63

PSR B1259-63 is the only known binary system with a radio pulsar from which the non-pulsed radio and X-ray emission was detected. The companion star in this system is a Be star SS 2883. A rapidly rotating radio pulsar is expected to produce a wind of relativistic particles. Be stars are known to produce highly asymmetric mass loss. Due to the interaction of the pulsar wind and the Be star wind the system of two shocks between the pulsar and the Be star forms. In this paper we show that the observed non-pulsed radio emission from the system is a result of the synchrotron emission of the relativistic particles in the outflow beyond the shock wave and that the non-pulsed X-ray emission is due to the inverse Compton scattering of the Be star photons on this particles. This revised version was published online in July 2006 with corrections to the Cover Date.     

Timing models for the long orbital period binary pulsar PSR B1259–63

The pulsar PSR B1259–63 is in a highly eccentric 3.4-yr orbit with the Be star SS 2883. Timing observations of this pulsar, made over a 7-yr period using the Parkes 64-m radio telescope, cover two periastron passages, in 1990 August and 1994 January. The timing data cannot be fitted by the normal pulsar and Keplerian binary parameters. A timing solution including a (non-precessing) Keplerian orbit and timing noise (represented as a polynomial of fifth order in time) provides a satisfactory fit to the data. However, because the Be star probably has a significant quadrupole moment, we prefer to interpret the data by a combination of timing noise, dominated by a cubic phase term, and ω^. and x ^. terms. We show that the ω^. and x ^. terms are likely to be a result of a precessing orbit caused by the quadrupole moment of the tilted companion star. We further rule out a number of possible physical effects which could contribute to the timing data of PSR B1259–63 on a measurable level.     

HESS J1616−508: likely to be powered by PSR J1617−5055

HESS J1616−508 is one of the brightest emitters in the TeV sky. Recent observations with the IBIS/ISGRI telescope onboard the INTEGRAL spacecraft have revealed that a young, nearby and energetic pulsar, PSR J1617−5055, is a powerful emitter of soft γ-rays in the 20–100 keV domain. In this paper, we present an analysis of all available data from the INTEGRAL , Swift , BeppoSAX and XMM–Newton telescopes with a view to assessing the most likely counterpart to the High Energy Stereoscopic System (HESS) source. We find that the energy source that fuels the X/γ-ray emissions is derived from the pulsar, both on the basis of the positional morphology, the timing evidence and the energetics of the system. Likewise the 1.2 per cent of the pulsar's spin-down energy loss needed to power the 0.1–10 TeV emission is also fully consistent with other HESS sources known to be associated with pulsars. The relative sizes of the X/γ-ray and very high energy sources are consistent with the expected lifetimes against synchrotron and Compton losses for a single source of parent electrons emitted from the pulsar. We find that no other known object in the vicinity could be reasonably considered as a plausible counterpart to the HESS source. We conclude that there is good evidence to assume that the HESS J1616−508 source is driven by PSR J1617−5055 in which a combination of synchrotron and inverse-Compton processes combine to create the observed morphology of a broad-band emitter from keV to TeV energies.     

TeV gamma-rays from accreting magnetars in massive binaries

This paper focuses on neutron stars (NS) of the magnetar type inside massive binary systems. We determine the conditions under which the matter from the stellar wind can penetrate the inner magnetosphere of the magnetar. At a certain distance from the NS surface, the magnetic pressure can balance the gravitational pressure of the accreting matter, creating a very turbulent, magnetized transition region. It is suggested that this region provides good conditions for the acceleration of electrons to relativistic energies. These electrons lose energy due to the synchrotron process and inverse Compton (IC) scattering of the radiation from the nearby massive stellar companion, producing high-energy radiation from X-rays up to ∼TeV γ-rays. The primary γ-rays can be further absorbed in the stellar radiation field, developing an IC  e^±  pair cascade. We calculate the synchrotron X-ray emission from primary electrons and secondary  e^±  pairs and the IC γ-ray emission from the cascade process. It is shown that quasi-simultaneous observations of the TeV γ-ray binary system LSI +61 303 in the X-ray and TeV γ-ray energy ranges can be explained with such an accreting magnetar model.     

The multi-band non-thermal emission from the supernova remnant RX J1713.7−3946

Non-thermal X-rays and very high energy (VHE) γ-rays have been detected from the supernova remnant (SNR) RX J1713.7−3946, and the recent observations with the Suzaku satellite clearly reveal a spectral cut-off in the X-ray spectrum, which directly relates to the cut-off of the energy spectrum of the parent electrons. However, whether the origin of the VHE γ-rays from the SNR is hadronic or leptonic is still in debate. We studied the multi-band non-thermal emission from RX J1713.7−3946 based on a semi-analytical approach towards the non-linear shock acceleration process by including the contribution of the accelerated electrons to the non-thermal radiation. The results show that the multi-band observations on RX J1713.7−3946 can be well explained in the model with appropriate parameters, and the TeV γ-rays have hadronic origin, i.e. they are produced via proton–proton interactions as the relativistic protons accelerated by the shock collide with the ambient matter.     

Gamma-ray binaries

Recent observations have shown that some compact stellar binaries radiate the highest energy light in the universe. The challenge has been to determine the nature of the compact object and whether the very high energy gamma-rays are ultimately powered by pulsar winds or relativistic jets. Multiwavelength observations have shown that one of the three gamma-ray binaries known so far, PSR B1259−63, is a neutron star binary and that the very energetic gamma-rays from this source and from another gamma-ray binary, LS I +61 303, may be produced by the interaction of pulsar winds with the wind from the companion star. At this time it is an open question whether the third gamma-ray binary, LS 5039, is also powered by a pulsar wind or a microquasar jet, where relativistic particles in collimated jets would boost the energy of the wind from the stellar companion to TeV energies. I.F. Mirabel is on leave from CEA, France.     

Modelling the thermal absorption and radio spectra of the binary pulsar B1259–63

We present the results of modelling of the radio spectrum evolution and dispersion measure variations of PSR B1259–63, a pulsar in a binary system with Be star LS 2883.We base our model on a hypothesis that the observed variations of the spectrum are caused by thermal free-free absorption occurring in the pulsar surroundings. We reproduce the observed pulsar spectral shapes in order to examine the influence of the stellar wind of LS 2883 and the equatorial disc on the pulsar’s radiation.The simulations of the pulsar’s radio emission and its consequent free-free absorption give us an insight into the impact of stellar wind and equatorial disc of LS 2883 has on the shapes of PSR B1259–63 radio spectra, providing an evidence for the connection between gigahertz-peaked spectra phenomenon and the close environment of the pulsar. Additionally, we supplement our model with an external absorbing medium, which results in a good agreement between simulated and observational data.     

Hard X-ray observations of PSR J1833−1034 and its associated pulsar wind nebula

PSR J1833−1034 and its associated pulsar wind nebula (PWN) have been investigated in depth through X-ray observations ranging from 0.1 to 200 keV. The low-energy X-ray data from Chandra reveal a complex morphology that is characterized by a bright central plerion, no thermal shell and an extended diffuse halo. The spectral emission from the central plerion softens with radial distance from the pulsar, with the spectral index ranging from  Γ= 1.61  in the central region to  Γ= 2.36  at the edge of the PWN. At higher energy, INTEGRAL detected the source in the 17–200 keV range. The data analysis clearly shows that the main contribution to the spectral emission in the hard X-ray energy range is originated from the PWN, while the pulsar is dominant above 200 keV. Recent High Energy Stereoscopic System (HESS) observations in the high-energy gamma-ray domain show that PSR J1833−1034 is a bright TeV emitter, with a flux corresponding to ∼2 per cent of the Crab in 1–10 TeV range. In addition, the spectral shape in the TeV energy region matches well with that in the hard X-rays observed by INTEGRAL . Based on these findings, we conclude that the emission from the pulsar and its associated PWN can be described in a scenario where hard X-rays are produced through synchrotron light of electrons with Lorentz factor  γ∼ 10⁹  in a magnetic field of ∼10 μG. In this hypothesis, the TeV emission is due to inverse-Compton interaction of the cooled electrons off the cosmic microwave background photons. Search for PSR J1833−1034 X-ray pulsed emission, via RXTE and Swift X-ray observations, resulted in an upper limit that is about 50 per cent.     

TeV γ-rays from old supernova remnants

We study the emission from an old supernova remnant (SNR) with an age of around 10⁵ yr and that from a giant molecular cloud (GMC) encountered by the SNR. When the SNR age is around 10⁵ yr, proton acceleration is efficient enough to emit TeV γ-rays both at the shock of the SNR and that in the GMC. The maximum energy of primarily accelerated electrons is so small that TeV γ-rays and X-rays are dominated by hadronic processes,  π⁰  -decay and synchrotron radiation from secondary electrons, respectively. However, if the SNR is older than several 10⁵ yr, there are few high-energy particles emitting TeV γ-rays because of the energy-loss effect and/or the wave-damping effect occurring at low-velocity isothermal shocks. For old SNRs or SNR–GMC interacting systems capable of generating TeV γ-ray emitting particles, we calculated the ratio of TeV γ-ray (1–10 TeV) to X-ray (2–10 keV) energy flux and found that it can be more than  ∼10²  . Such a source showing large flux ratio may be a possible origin of recently discovered unidentified TeV sources.     

Gamma-rays and cosmic rays from a pulsar in Cygnus OB2

We argue that γ-ray sources observed in the direction of the Cygnus OB2 association in the GeV and TeV energy range are due to a pulsar that was created by a supernova a few tens of thousands of years ago. The GeV emission is produced by a middle-aged pulsar, a factor of 2 older than the Vela pulsar. The TeV emission is produced by high-energy hadrons and/or leptons accelerated in pulsar wind nebulae. We suggest, moreover, that the excess of cosmic rays at ∼10¹⁸ eV observed from the direction of the Cygnus region can also be related to the appearance of this very energetic pulsar in the Cyg OB2 association. Some of the relativistic hadrons, captured in strong magnetic fields of a high-density region of Cyg OB2, produce neutrons and γ-rays in collisions with matter. These neutrons can arrive from Cyg OB2, creating an excess of cosmic rays.     

High-energy γ-rays from globular clusters

It is expected that specific globular clusters (GCs) can contain up to a hundred of millisecond pulsars. These pulsars can accelerate leptons at the shock waves originated in collisions of the pulsar winds and/or inside the pulsar magnetospheres. Energetic leptons diffuse gradually through the GC Comptonizing stellar and microwave background radiation. We calculate the GeV–TeV γ-ray spectra for different models of injection of leptons and parameters of the GCs assuming reasonable, of the order of 1 per cent, efficiency of energy conversion from the pulsar winds into the relativistic leptons. It is concluded that leptons accelerated in the GC cores should produce well localized γ-ray sources which are concentric with these GCs. The results are shown for four specific GCs (47 Tuc, Ter 5, M13 and M15), in which significant population of millisecond pulsars have been already discovered. We argue that the best candidates, which might be potentially detected by the present Cherenkov telescopes and the planned satellite telescopes (AGILE, GLAST), are 47 Tuc on the Southern hemisphere, and M13 on the Northern hemisphere. We conclude that detection (or non-detection) of GeV–TeV γ-ray emission from GCs by these instruments put important constraints on the models of acceleration of leptons by millisecond pulsars.     

PSR J1740−3052: a pulsar with a massive companion

We report on the discovery of a binary pulsar, PSR J1740−3052, during the Parkes multibeam survey. Timing observations of the 570-ms pulsar at Jodrell Bank and Parkes show that it is young, with a characteristic age of 350 kyr, and is in a 231-d, highly eccentric orbit with a companion whose mass exceeds 11 M_⊙. An accurate position for the pulsar was obtained using the Australia Telescope Compact Array. Near-infrared 2.2-μm observations made with the telescopes at the Siding Spring observatory reveal a late-type star coincident with the pulsar position. However, we do not believe that this star is the companion of the pulsar, because a typical star of this spectral type and required mass would extend beyond the orbit of the pulsar. Furthermore, the measured advance of periastron of the pulsar suggests a more compact companion, for example, a main-sequence star with radius only a few times that of the Sun. Such a companion is also more consistent with the small dispersion measure variations seen near periastron. Although we cannot conclusively rule out a black hole companion, we believe that the companion is probably an early B star, making the system similar to the binary PSR J0045−7319.     

Very-high-energy gamma radiation associated with the unshocked wind of the Crab pulsar

We show that the relativistic wind of the Crab pulsar, which is commonly thought to be invisible in the region upstream of the termination shock at r r _S∼0.1 pc, in fact could be directly observed through its inverse Compton (IC) γ -ray emission. This radiation is caused by illumination of the wind by low-frequency photons emitted by the pulsar, and consists of two, pulsed and unpulsed , components associated with the non-thermal (pulsed) and thermal (unpulsed) low-energy radiation of the pulsar, respectively. These two components of γ -radiation have distinct spectral characteristics, which depend essentially on the site of formation of the kinetic-energy-dominated wind, as well as on the Lorentz factor and the geometry of propagation of the wind. Thus, the search for such specific radiation components in the spectrum of the Crab Nebula can provide unique information about the unshocked pulsar wind that is not accessible at other wavelengths. In particular, we show that the comparison of the calculated flux of the unpulsed IC emission with the measured γ -ray flux of the Crab Nebula excludes the possibility of formation of a kinetic-energy-dominated wind within 5 light-cylinder radii of the pulsar, R _w5 R _L. The analysis of the pulsed IC emission, calculated under reasonable assumptions concerning the production site and angular distribution of the optical pulsed radiation, yields even tighter restrictions, namely R _w30 R _L.