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.     

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.     

The 1997 periastron passage of the binary pulsar PSR B1259−63

We report here on multifrequency radio observations of the pulsed emission from PSR B1259−63 around the time of the closest approach (periastron) to its B2e companion star. There was a general increase in the dispersion measure (DM) and scatter-broadening of the pulsar, and a decrease in the flux density towards periastron although fluctuation in these parameters were seen on time-scales as short as minutes. The pulsed emission disappeared 16 d prior to periastron and remained undetectable until 16 d after periastron.
The observations are used to determine the parameters of the wind from the Be star. We show that a simple model, in which the wind density varies with radius as r ⁻², provides a good fit to the data. The wind is highly turbulent with an outer scale of ≤10¹⁰ cm and an inner scale perhaps as small as 10⁴ cm, a mean density of ∼10⁶ cm⁻³ and a velocity of ∼2000 km s⁻¹ at a distance of ∼50 stellar radii. We find a correlation between DM variations and the pulse scattering times, suggesting that the same electrons are responsible for both effects.     

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.      

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.     

γ-rays from binary system with energetic pulsar and Be star with aspherical wind: PSR B1259−63/SS2883

At least one massive binary system containing an energetic pulsar, PSR B1259−63/SS2883, has been recently detected in the TeV γ-rays by the HESS telescopes. These γ-rays are likely produced by particles accelerated in the vicinity of the pulsar and/or at the pulsar wind shock, in comptonization of soft radiation from the massive star. However, the process of γ-ray production in such systems can be quite complicated due to the anisotropy of the radiation field, complex structure of the pulsar wind termination shock and possible absorption of produced γ-rays which might initiate leptonic cascades. In this paper, we consider in detail all these effects. We calculate the γ-ray light curves and spectra for different geometries of the binary system PSR B1259−63/SS2883 and compare them with the TeV γ-ray observations. We conclude that the leptonic inverse-Compton model, which takes into account the complex structure of the pulsar wind shock due to the aspherical wind of the massive star, can explain the details of the observed γ-ray light curve.     

One large glitch in PSR B1737-30 detected with the TMRT

One large glitch was detected in PSR B1737-30 using data spanning from MJD 57999 to 58406 obtained with the newly built Shanghai Tian Ma Radio Telescope(TMRT). The glitch took place at the time around MJD 58232.4 when the pulsar underwent an increase in the rotation frequency of △v about 1.38 × 10^-6 Hz, corresponding to a fractional step change of △v/v～ 8.39 × 10^-7. Post-glitch v gradually decreased to the pre-glitch value. The frequency derivative was observed to undergo a step change of about-9×10^-16 s^-2. Since July 1987, there have been 36 glitches already reported in PSR B1737-30 including this one. According to our analysis, the glitch size distribution is well described by a power law with an index of 1.13. The distribution of the interval between two adjacent glitches(waiting time △T) follows a Poisson probability density function. For PSR B1737-30, the interval is prone to be long after a large glitch. However, no correlation is found between glitch size and the interval since the previous glitch.     

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.     

Giant radio pulses from the millisecond pulsar PSR B1937+21 at 327 MHz

Seven giant radio pulses were recorded from the millisecond pulsar PSR B1937+21 during ≈8.1 min observation by the Ooty Radio Telescope (ORT) at 326.5 MHz. Although sparse, these observations support most of the giant pulse behaviour reported at higher radio frequencies (430 to 2380 MHz). Within the main component of the integrated profile, they are emitted only in a narrow (≲47 μs) window of pulse phase, close to its peak. This has important implications for doing super-high precision timing of PSR B1937+21 at low radio frequencies.     

A study of subpulse phase correlation across nulls in PSR B0031−07

The correlation of subpulse phases across nulls is investigated in the radio pulsar PSR B0031−07, using 29 849 periods of high-quality data obtained with the Ooty Radio Telescope (ORT) which operates at 327 MHz. Assuming that the turn-off and turn-on subpulse phases (the phase of the subpulse in the last period before the null and that in the first period after the null, respectively) are independent random variables, the expected distribution of their difference (i.e. the total drift) is inconsistent with the observed distribution for null transitions within the same drift mode; this implies a correlation of subpulse phase across nulls. However, this correlation decreases with null duration for both the dominant drift modes. Substantial drifting occurs during short nulls (one to four periods); the drift rate during the short nulls appears to be constant for a class A transition, whereas it decreases with null duration for class B transitions. These results, together with the reported behaviour of PSR B1944+17 and PSR B0809+74, seem to imply different time-scales for phase correlation in different pulsars.     

Observations of six glitches in PSR B1737−30

Six glitches have been recently observed in the rotational frequency of the young pulsar PSR B1737−30 (J1740−3015) using the 25-m Nanshan telescope of Urumqi Observatory. With a total of 20 glitches in 20 yr, it is one of the most frequently glitching pulsars of the ∼1750 known pulsars. Glitch amplitudes are very variable with fractional increases in rotation rate ranging from 10⁻⁹ to 10⁻⁶. Interglitch intervals are also very variable, but no relationship is observed between interval and the size of the preceding glitch. There is a persistent increase in     , opposite in sign to that expected from slowdown with a positive braking index, which may result from changes in the effective magnetic dipole moment of the star during the glitch.     

The light curve of the companion to PSR B1957+20

We present a new analysis of the light curve for the secondary star in the eclipsing binary millisecond pulsar system PSR B1957+20. Combining previous data and new data points at minimum from the Hubble Space Telescope , we have 100 per cent coverage in the R -band. We also have a number of new K _s-band data points, which we use to constrain the infrared magnitude of the system. We model this with the eclipsing light-curve (ELC) code. From the modelling with the ELC code we obtain colour information about the secondary at minimum light in BVRI and K . For our best-fitting model we are able to constrain the system inclination to 65°± 2° for pulsar masses ranging from 1.3 to  1.9  M_⊙  . The pulsar mass is unconstrained. We also find that the secondary star is not filling its Roche lobe. The temperature of the unirradiated side of the companion is in agreement with previous estimates and we find that the observed temperature gradient across the secondary star is physically sustainable.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.