Binary systems represent an important population of very high energy (VHE) sources. Presently this category comprises four objects [1, 2, 3, 4, 5]. At least one of them, PSR B1259-63/SS2883, is a binary pulsar system consisting of a pulsar in an elliptic orbit around a massive B2e optical star . The multiwavelength observations of the system [10, 8, 11, 7, 13, 9, 12] indicate correlations between different energy bands . One of the interesting features of the high energy emission of PSR B1259-63/SS2883 is the behavior of the gamma-ray light-curve which deviates from the earlier theoretical predictions [15, 16]. The understanding of the nature/origin of the high energy radiation of this object is a quite complex issue. It requires not only proper modeling of the acceleration and radiation process in the context of multiwavelength properties of the source, but also detailed studies of the hydrodynamics of the interaction of the pulsar and stellar winds. Indeed, the particle acceleration regime and formation of energy distribution of relativistic particles after termination of the pulsar wind strongly depend on the dynamics of the flow. In this regard one should note that the recent attempts at explanation of the gamma-ray light-curve of PSR B1259-63/SS2883 [14, 17] contain assumptions which are closely linked to the hydrodynamics of the system. For example, the density changes along current lines determine the adiabatic loss rate in the flow. The flow bulk motion is also responsible for the particle advection (escape) from the system. The nonradiative losses caused by these processes may have a strong impact on the energy spectrum and light-curve of gamma-rays . Therefore, any self-consistent treatment of the high energy radiation of the system requires hydrodynamic calculations of key parameters characterizing the interaction of two winds. This concerns, in particular, the magnetic field. The magnetic field is expected to be nonhomogenous in the pulsar wind nebulae . The calculations of the spatial distribution of the strength of the magnetic field generally requires a proper magnetohydrodynamical (MHD) approach. However, in the case when the the impact of the magnetic field on the dynamics of plasma is not dramatic, it can be treated hydrodynamically within a “froze-in” approach. Below we assume the influence of the magnetic field to be weak and we discuss the impact of hydrodynamics characterizing the collision of pulsar and stellar winds on the high energy radiation of the system. The results are based on the recent hydrodynamic calculations of Bogovalov et al.  conducted through two different approaches using a relativistic code in the “pulsar zone” and a nonrelativistic code in the “optical star” region. It is assumed that initially both the pulsar and stellar winds expand radially.
2 Hydrodynamics of interaction of relativistic and nonrelativistic winds
The discussion below concerns the regions of linear scales comparable to the size of the binary system (because the variable synchrotron and IC radiation components are predominantly generated in this region). The results are based on the hydrodynamical calculations reported by Bogovalov et al. . The properties of the steady-state flow basically depend only on one parameter, namely on the wind ram pressure ratio :
Here is the pulsar spindown luminosity; and are the stellar mass loss rate and the stellar outflow velocity, respectively (both winds are assumed to be isotropic). In the particular case of PSR B1259-63/SS2883 the parameter is expected to be within – . Depending on the value of , the structure of the flow may significantly change. It is demonstrated in Figs. 1 and 2 for and (the color represents the flow bulk Lorentz factor).
Figures 1 and 2 show that at distances comparable to the size of the binary system, the bulk Lorentz factor of the shocked flow may be as high as . Because of Doppler boosting, this should result in a significant anisotropy of nonthermal emission of relativistic particles. Moreover, the large bulk Lorentz factors imply significant adiabatic losses. Indeed, the full energy per particle is conserved along the flow line (here and are the enthalpy and particle densities). Thus the acceleration of the flow is reduced to the transfer of the thermal energy to the bulk motion, i.e. to adiabatic losses. Below we compare the adiabatic loss time to other relevant time-scales, the particle escape time and IC cooling time. Note that all these time-scales depend, unlike the overall structure of the flow, on the separation distance between stars.
As long as the post-shock flow velocity is close to the speed of light , the advection time scale can be estimated as
It is convenient to express all times in the unit of . In Figure 3 we show the IC cooling time for three different separations, which correspond to , and days to the periastron passage. The calculations are performed for the following parameters: temperature of the optical star K; eccentricity of the orbit ; periastron separation cm. As it follows from Fig. 3, the inverse Compton cooling time , i.e. the radiative cooling time is much larger than the particle escape time, especially at large separation distances (obviously, , due to decrease of target photon density).
Although it is not possible to obtain, from first principles, the adiabatic loss rate, the -dependence of the adiabatic cooling time can be defined. Interestingly, the adiabatic loss time has a -dependence similar to . Indeed, since adiabatic losses, which occur between two points on the current line, depend only on the ratio of densities in these points, one obtains
since the traveling time between these two points is proportional to . This allows us to determine the -dependence of the ratio of IC to nonradiative cooling time:
Thus, in a case of collision of two isotropic winds small separation distances are more preferable for VHE gamma-ray production. Although, we note here that it does not necessarily imply higher observable fluxes, due to the anisotropy related to the Dopper boosting of the produced emission.
Generally, the adiabatic loss rate is not homogenous in the flow. In Figure 4 we show typical dependences of adiabatic cooling time, obtained with numerical simulation, along current lines for two cases: and (please note that sharp peaks/deeps in this figure are numerical artifacts). These calculations were performed along current lines close to the contact discontinuity in the “pulsar zone”.
The calculations show that for the flow structure predicted by the hydrodynamical calculations, one should expect very fast (compared to the IC cooling time) adiabatic losses over the entire orbit. This should result in a strong suppression of nonthermal radiation of accelerated particles.
The results presented above demonstrate the importance of hydrodynamics for the calculation of the high energy radiation of particles accelerated in the post-shock flow. Namely, (i) the adiabatic cooling time is significantly shorter than IC cooling time; (ii) the bulk Lorentz factor of the flow is rather large () already at distances, comparable to the size of the binary system. Both these effects have direct implication for interpretation of observations of high energy gamma- and X-ray emission components from the binary pulsar system PSR B1259-63/SS2883. In particular, the derived dependence in Eq. (4) (the ratio of IC cooling time to the nonradiative loss time as a function of the separation distance) implies that the HD simulations do not support (at least in the particular case of two isotropic winds colliding) the assumption of a sharp increase in the role of adiabatic losses towards periastron as hypothesized in ref.  for explanation of the observed deep in the gamma-ray light-curve of PSR B1259-63/SS2883 around periastron. On the other hand, the obtained absolute timescales of the adiabatic losses (compared to the Compton cooling times) show that the VHE signal should be significantly suppressed, especially at large separation distances. The HESS observation of PSR B1259-63/SS2883  supports this conclusion: the signal from the source was detected only during epochs close to the periastron passage.
The effect of the flow bulk acceleration strongly modifies the relationship between the synchrotron X-ray and inverse Compton gamma-ray fluxes produced by the same population of relativistic electrons. Indeed, as long as the magnetic field is frozen in the plasma, its strength is determined by the structure of the flow and may significantly vary across the flow. The plasma bulk motion affects the IC radiation production as well, but in a different way. Namely, in the co-moving frame, where electrons are distributed isotropically, the target radiation field will be Doppler boosted. This implies that hydrodynamical effects can lead to a significant deviation from the standard relations between the X-ray and VHE gamma-ray fluxes (see e.g. [20, 21]).
Moreover, due to the large bulk Lorentz factors (for details see ), the observed fluxes should have strong orbital phase dependence. Indeed, the direction of the post shock flow varies with the motion of the pulsar along the orbit around the star. This implies significant changes of the Doppler factor . Namely, for small viewing angles (when the flow bulk velocity is directed towards the observer); and for large viewing angles. Such a geometry is expected, in particular, at the periastron passage. Correspondingly, this will have a strong impact on the light-curve of nonthermal radiation of electrons, where typically . Interestingly, in such a case one should expect a direct correlation between different energy bands (from radio to VHE), even though the electrons responsible for these energy intervals of electromagnetic radiation may not have the same origin. This conclusion is supported by multiwavelength observation of PSR B1259-63/SS2883. Moreover, this effect could be the reason for the deep in the VHE light-curve close to the periastron as reported by the HESS collaboration .
- F. Aharonian et al. (HESS Collaboration), A&A 442, 1 (2005).
- F. Aharonian et al. (HESS Collaboration), Science 309, 746 (2005).
- F. Aharonian et al. (HESS Collaboration), A&A 460, 743 (2006).
- J. Albert et al. (MAGIC Collaboration), Science 312, 1771A (2006).
- J. Albert et al. (MAGIC Collaboration), ApJ 665, L51 (2007).
- S. Johnston et al., ApJ 387L, 37 (1992).
- L. Cominsky et al., ApJ 427, 978 (1994).
- T.W. Connors et al., MNRAS 336, 1201 (2002).
- M. Hirayama et al., ApJ 521, 718 (1999).
- S. Johnston et al., MNRAS 302, 277 (1999).
- S. Johnston et al., MNRAS 358, 1069 (2005).
- M. Chernyakova et al., MNRAS 367, 1201C (2006).
- V.M. Kaspi et al., ApJ 453, 424 (1995).
- A. Neronov and M. Chernyakova, Ap&SS 309, 253N (2007).
- J. Kirk et al., Astroparticle Physics 10, 31 (1999).
- A. Kawachi et al., ApJ 607, 949 (2004).
- D. Khangulyan et al., MNRAS 380, 320 (2007).
- C.F. Kennel and F.V. Coroniti, ApJ 283, 710 (1984).
- S.V. Bogovalov et al., arXiv:0710.1961v1.
- F. Aharonian et al., MNRAS 291, 162 (1997).
- D. Khangulyan and F. Aharonian, “High Energy Gamma-Ray Astronomy”, AIP Conference Proceedings Vol.745, pp.359-364 (2005).