Pulse profile stability of the Crab pulsar

We present an X-ray timing analysis of the Crab pulsar, PSR B0531+21, using archival RXTE data. We have investigated the stability of the Crab pulse profile, in soft (2-20 keV) and hard (30-100 keV) X-ray energies, over the last decade of RXTE operation. The analysis includes measurement of the separation between the two pulse peaks and the intensity and widths of the two peaks. We did not find any significant time dependency in the pulse shape. The two peaks have been stable in phase, intensity and width f...     

Orbital evolution measurement of the accreting millisecond X-ray pulsar SAX J1808.4-3658

We present results from a pulse timing analysis of the accretion-powered millisecond X-ray pulsar SAX J1808.4-3658 using X-ray data obtained during four outbursts of this source. Extensive observations were made with the proportional counter array of the Rossi X-ray Timing Explorer (RXTE) during the four outbursts that occurred in 1998, 2000, 2002 and 2005. Instead of measuring the arrival times of individual pulses or the pulse arrival time delay measurement that is commonly used to determine the orbital parameters of binary pulsars, we have determined the orbital ephemeris during each observation by optimizing the pulse detection against a range of trial ephemeris values. The source exhibits a significant pulse shape variability during the outbursts. The technique used by us does not depend on the pulse profile evolution, and is therefore, different from the standard pulse timing analysis. Using 27 measurements of orbital ephemerides during the four outbursts spread over more than 7 years and more than 31,000 binary orbits, we have derived an accurate value of the orbital period of 7249.156862(5) s (MJD = 50915) and detected an orbital period derivative of (3.14 ± 0.21) × 10⁻¹² s s⁻¹. We have included a table of the 27 mid-eclipse time measurements of this source that will be valuable for further studies of the orbital evolution of the source, especially with ASTROSAT. We point out that the measured rate of orbital period evolution is considerably faster than the most commonly discussed mechanisms of orbital period evolution like mass transfer, mass loss from the companion star and gravitational wave radiation. The present time scale of orbital period change, 73 Myr is therefore likely to be a transient high value of period evolution and similar measurements during subsequent outbursts of SAX J1808.4-3658 will help us to resolve this.     

The steady spin-down rate of 4U 1907+09

Using X-ray data from the Rossi X-ray Timing Explorer , we report the pulse timing results of the accretion-powered, high-mass X-ray binary pulsar 4U 1907+09, covering a time-span of almost two years. We measured three new pulse periods in addition to the previously measured four pulse periods. We are able to connect pulse arrival times in phase for more than a year. The source has been spinning down almost at a constant rate, with a spin-down rate of     for more than 15 yr. Residuals of pulse arrival times yield a very low level of random-walk noise, with a strength of ∼     on a time-scale of 383 d, which is 40 times lower than that of the high-mass X-ray binary pulsar Vela X-1. The noise strength is only a factor of 5 greater than that of the low-mass X-ray binary pulsar 4U 1626−67. The low level of the timing noise and the very stable spin-down rate of 4U 1907+09 make this source unique among the high-mass X-ray binary pulsars, providing another example, in addition to 4U 1626−67, of long-term quiet spin down from an accreting source. These examples show that the extended quiet spin-down episodes observed in the anomalous X-ray pulsars 1RXS J170849.0−400910 and 1E 2259+586 do not necessarily imply that these sources are not accreting pulsars.     

Search for orbital motion of the pulsar 4U 1626–67: Candidate for a neutron star with a supernova fall-back accretion disk

We report here results from a new search for orbital motion of the accretion powered X-ray pulsar 4U 1626–67 using two different analysis techniques. X-ray light curve obtained with the Proportional Counter Array of the Rossi X-ray Timing Explorer during a long observation carried out in February 1996, was used in this work. The spin period and the local period derivative were first determined from the broad 2–60 keV energy band light curve and these were used for all subsequent timing analysis. In the first technique, the orbital phase dependent pulse arrival times were determined for different trial orbital periods in the range of 500 to 10,000 s. We have determined a 3σ upper limit of 13 lt-ms on the projected semimajor axis of the orbit of the neutron star for most of the orbital period range, while in some narrow orbital period ranges, covering about 10% of the total orbital period range, it is 20lt-ms. In the second method, we have measured the pulse arrival times at intervals of 100 s over the entire duration of the observation. The pulse arrival time data were used to put an upper limit on any periodic arrival time delay using the Lomb-Scargle periodogram. We have obtained a similar upper limit of 10 lt-ms using the second method over the orbital period range of 500–10,000 s. This puts very stringent upper limits for the mass of the compact object except for the unlikely case of a complete face-on orientation of the binary system with respect to our line-of-sight. In the light of this measurement and the earlier reports, we discuss the possibility of this system being a neutron star with a supernovae fall-back accretion disk.     

Analysis of the Precision of Pulsar Time Clock Modeltwo

Millisecond pulsars have a very high rotation stability, which can be applied to many research fields, such as the establishment of the pulsar time standard, the detection of gravitational wave, the spacecraft navigation by using X-ray pulsars and so on. In this paper, we employ two millisecond pulsars PSR J0437-4715 and J1713+0743, which are observed by the International Pulsar Timing Array (IPTA), to analyze the precision of pulsar clock parameter and the prediction accuracy of pulse time of arrival (TOA). It is found that the uncertainty of spin frequency is 10^?15 Hz, the uncertainty of the first derivative of spin frequency is 10^?23 s^?2, and the precision of measured rotational parameters increases by one order of magnitude with the accumulated observational data every 4～5 years. In addition, the errors of TOAs within 4.8 yr which are predicted by the clock model established by the 10 yr data of J0437-4715 are less than 1 μs. Therefore, one can use the pulsar time standard to calibrate the atomic clock, and make the atomic time deviate from the TT (Terrestrial Time) less than 1 μs within 4.8 yr.     

A method for judging decay or growth of the magnetic field of pulsar

This paper provides a method for judging growth or decay of the magnetic field of pulsar by using pulse period P, or frequency ν, and its first and second derivatives $ \dot P,\ddot P $ or $ \dot v,\ddot v $ . The author uses this method to judge the growth or decay of the magnetic field of Crab pulsar. The judged result for Crab pulsar is that the magnetic field of Crab pulsar is growing now, but it is not decaying. The result corresponds with the actual case of Crab pulsar.     

Maximum-likelihood TOA Estimation of X-ray Pulsar Signals on the Basis of Poison Model

The cycle-stationary Poison model for the photon arrival of X-ray pulsars and the estimation of the pulse's time of arrival (TOA) are discussed. Based on this model, the maximum likelihood estimation (MLE) of the TOA, as well as the Cramer-Rao boundary (CRB), are presented, especially the approximate formulation of the likelihood function and the Cramer-Rao boundary are derived for the low-R_SN (signal-to-noise ratio) cases. By using the analytical pulse profile, we have made the Monte-Carlo simulation on the TOA estimation of PSR B1821-24, discussed the estimate error for different observation times and signal-to-noise ratios, and presented the corresponding R_SN-thresholds. The results show that this method of data analysis can estimate effectively the timing accuracy of X-ray pulsar pulses and help to evaluate the corresponding performances in other applications.     

The accretion powered spin-up of GRO J1750−27

The timing properties of the 4.45 s pulsar in the Be X-ray binary system GRO J1750−27 are examined using hard X-ray data from INTEGRAL and Swift during a type II outburst observed during 2008. The orbital parameters of the system are measured and agree well with those found during the last known outburst of the system in 1995. Correcting the effects of the Doppler shifting of the period, due to the orbital motion of the pulsar, leads to the detection of an intrinsic spin-up that is well described by a simple model including     and     terms of  −7.5 × 10⁻¹⁰ s s⁻¹  and  1 × 10⁻¹⁶ s s⁻²  , respectively. The model is then used to compare the time-resolved variation of the X-ray flux and intrinsic spin-up against the accretion torque model of Ghosh & Lamb; this finds that GRO J1750−27 is likely located 12–22 kpc distant and that the surface magnetic field of the neutron star is  ∼2 × 10¹²  G. The shape of the pulse and the pulsed fraction shows different behaviour above and below 20 keV, indicating that the observed pulsations are the convolution of many complex components.     

Short time scale pulse stability of the Crab pulsar in the optical band

The fine structure and the variations of the optical pulse shape and phase of the Crab pulsar are studied on various time scales. The observations have been carried out on 4-m William Hershel and 6-m BTA telescopes with APD photon counter, photomultiplier based 4-channel photometer and PSD based panoramic spectrophotopolarimeter with 1 μs time resolution in 1994, 1999, 2003 and 2005–2006 years. The upper limit on the pulsar precession on Dec 2, 1999 is placed in the 10 s–2 hr time range. The evidence of a varying from set to set fine structure of the main pulse is found in the 1999 and 2003 years data. No such fine structure is detected in the integral pulse shape of 1994, 1999 and 2003 years. The drastic change of the pulse shape in the 2005–2006 years set is detected along with the pulse shape variability and quasi-periodic phase shifts. This work has been supported by the Russian Foundation for Basic Research (grant No 04-02-17555), Russian Academy of Sciences (program “Evolution of Stars and Galaxies”), by the Russian Science Support Foundation, and by INTAS (grant No 04-78-7366).     

关于脉冲星脉冲到达时间转换方程

较详细地介绍了脉冲星脉冲到达时间（TOA）转换方程,讨论了脉冲星TOA转换方程在航天器导航算法中的具体应用问题。同时,对导航用的脉冲星脉冲TOA转换方程与地面射电计时观测用的脉冲星脉冲TOA转换方程进行了比较研究。     

Mechanism of generation of the emission bands in the dynamic spectrum of the Crab pulsar

We show that the proportionately spaced emission bands in the dynamic spectrum of the Crab pulsar fit the oscillations of the square of a Bessel function whose argument exceeds its order. This function has already been encountered in the analysis of the emission from a polarization current with a superluminal distribution pattern: a current whose distribution pattern rotates (with an angular frequency ω) and oscillates (with a frequency  Ω > ω  differing from an integral multiple of ω) at the same time. Using the results of our earlier analysis, we find that the dependence on frequency of the spacing and width of the observed emission bands can be quantitatively accounted for by an appropriate choice of the value of the single free parameter  Ω/ω  . In addition, the value of this parameter, thus implied by Hankins & Eilek's data, places the last peak in the amplitude of the oscillating Bessel function in question at a frequency  (∼Ω³/ω²)  that agrees with the position of the observed ultraviolet peak in the spectrum of the Crab pulsar. We also show how the suppression of the emission bands by the interference of the contributions from differing polarizations can account for the differences in the time and frequency signatures of the interpulse and the main pulse in the Crab pulsar. Finally, we put the emission bands in the context of the observed continuum spectrum of the Crab pulsar by fitting this broad-band spectrum (over 16 orders of magnitude of frequency) with that generated by an electric current with a superluminally rotating distribution pattern.     

High time-resolution spectroscopic imaging using intensified CCD detectors

The micro‐channel plate intensified CCD photon counting detector developed at University College London has been upgraded to allow time-resolved spectroscopic optical data to be acquired on periodical sources such as pulsars. First observing trials have been carried out, acquiring spectroscopic data on the Crab pulsar. The detector was phase locked to the pulsar period and a temporal resolution of 41.4 μs employed. The phase locking allowed the coaddition of time slices over a large number of pulsar periods to build up quantifiable spectroscopic data when observing in a flux-limited regime.     

A search for soft X-ray radiation from pulsars with the Astronomical Netherlands Satellite

We present the results of a search for X-ray emission in the energy range 0.2–0.28 keV and 1–7 keV from a number of radio pulsars, including Crab, Vela and the binary pulsar PSR 1913+16, using the soft X-ray experiment aboard ANS. Except for the Crab no pulsed flux has been found. From the Vela pulsar we have detected continuous flux in agreement with earlier observations. Upper limits are given.Paper presented at the COSPAR Symposium on Fast Transients in X-and Gamma-Rays, held at Varna, Bulgaria, 29–31 May, 1975.     

Two-frequency timing of the pulsar B1937+21 in Kalyazin and Kashima in 1997–2002

We present the results from our timing of the millisecond pulsar B1937+21, performed jointly since 1997 on two radio telescopes: the RT-64 in Kalyazin (Russia) at a frequency of 0.6GHz and RT-34 in Kashima (Japan) at a frequency of 2.15 GHz. The rms value of the pulse time of arrival (TOA) residuals for the pulsar at the barycenter of the Solar system is 1.8 μs (the relative variation is ≈10^?14 over the observing period). The TOA residuals are shown to be dominated by white phase noise, which allows this pulsar to be used as an independent time scale keeper. The upper limit for the gravitational background energy density Ω_gh² at frequencies ≈6.5 × 10^?9 Hz is estimated to be no higher than 10^?6. Based on the long-term timing of the pulsar, we have improved its parameters and accurately determined the dispersion measure and its time variation over the period 1984–2002, which was, on average, ?0.00114(3) pc cm^?3 yr^?1.     

Search for TeV gamma-rays from Geminga pulsar

We present the results of observation of the Geminga pulsar carried out in the TeV energy band during the 6 year period spanning 2000–2006 using the Pachmarhi Array of Cherenkov Telescopes (PACT). A long stretch of data, new computer codes and the “Tempo” package have been used in the present analysis. We have searched for evidence of pulsed emission of γ-rays from the Geminga pulsar using the post-glitch pulsar elements obtained by Jackson and Halpern from X-ray/γ-ray satellite data. We do not see any significant evidence for pulsed emission from the Geminga pulsar at a threshold energy of 825 GeV. In this paper we present our results on the light curve in the TeV energy band, set an upper limit on the time averaged flux of γ-rays, and compare our results with other ground based observations.     

国家授时中心40米射电望远镜Crab脉冲星周期跃变监测

脉冲星周期跃变是一种罕见的现象,是研究其内部结构的探针。针对2019年2月~12月国家授时中心昊平观测站40 m射电望远镜在脉冲星计时观测中监测Crab脉冲星的数据,采用脉冲星计时方法,用TEMPO2拟合程序进行分析。结果表明,Crab脉冲星在2019年7月23日(MJD 58687)附近发生了一次周期跃变现象,该跃变自转增量为Δv_g=5.33(4)×10^-7Hz,自转变化量为Δv_g/v=17.9(1)×10^-9,并伴随着恢复系数Q~0.88的指数恢复过程。此次Crab脉冲星周期跃变的监测及处理,证实了40 m射电望远镜对脉冲星的监测性能,同时为研究周期跃变的产生机理积累了样本。     

Detecting the errors in solar system ephemeris by pulsar timing

Pulsar timing uses planetary ephemerides to convert the measured pulse arrival time at an observatory to the arrival time at the Solar System barycenter(SSB). Since these planetary ephemerides cannot be perfect, a method of detecting the associated errors based on a pulsar timing array is developed. By using observations made by an array of 18 millisecond pulsars from the Parkes Pulsar Timing Array, we estimated the vector uncertainty from the Earth to the SSB of JPL DE421, which reflects the offset of the ephemeris origin with respect to the ideal SSB, in different piecewise intervals of pulsar timing data, and found consistent results. To investigate the stability and reliability of our method, we divided all the pulsars into two groups. Both groups yield largely consistent results, and the uncertainty of the Earth-SSB vector is several hundred meters, which is consistent with the accuracy of JPL DE421. As an improvement in the observational accuracy, pulsar timing will be helpful to improve the solar system ephemeris in the future.     

X-ray observations of PSR B0355+54 and its pulsar wind nebula

We present X-ray data of the middle-aged radio pulsar PSR B0355+54. The XMM-Newton and Chandra observations show not only emission from the pulsar itself, but also compact diffuse emission extending ∼50″ in the opposite direction to the pulsar’s proper motion. Our analysis also indicates the presence of fainter diffuse emission extending ∼5′ from the point source. The morphology of the diffuse component is similar to the ram-pressure confined pulsar wind nebulae detected for other sources. We find that the compact diffuse component is well-fitted with a power-law, with an index that is consistent with the values found for other pulsar wind nebulae. The core emission from the pulsar can be characterized with a thermal plus power-law fit, with the thermal emission most likely originating in a hot polar cap.     

Broadband observations of the transient X-ray pulsar SAX J2103.5+4545

We investigated the optical, X-ray, and gamma-ray variability of the pulsar SAX J2103.5+4545. Our timing and spectral analyses of the X-ray and gamma-ray emissions from the source using RXTE and INTEGRAL data show that the shape of its spectrum in the energy range 3–100 keV is virtually independent of its intensity and the orbital phase. Based on XMM-Newton data, we accurately (5″) localized the object and determined the optical counterpart in the binary. We placed upper limits on the variability of the latter in the Hα line on time scales of the orbital and pulse periods, respectively.     

Formation of the radio profile components of the Crab pulsar

The induced Compton scattering of radio emission off the particles of the ultrarelativistic electron–positron plasma in the open field line tube of a pulsar is considered. We examine the scattering of a bright narrow radio beam into the background over a wide solid angle and specifically study the scattering in the transverse regime, which holds in a moderately strong magnetic field and gives rise to the scattered component nearly antiparallel to the streaming velocity of the scattering particles. Making use of the angular distribution of the scattered intensity and taking into account the effect of rotational aberration in the scattering region, we simulate the profiles of the backscattered components as applied to the Crab pulsar. It is suggested that the interpulse (IP), the high-frequency interpulse (IP') and the pair of so-called high-frequency components (HFC1 and HFC2) result from the backward scattering of the main pulse (MP), precursor (PR) and low-frequency component (LFC), respectively. The components of the high-frequency profiles, the IP' and HFCs, are interpreted for the first time. The HFC1 and HFC2 are argued to be a single component split by the rotational aberration close to the light cylinder. It is demonstrated that the observed spectral and polarization properties of the profile components of the Crab pulsar as well as the giant pulse phenomenon outside the MP can be explained in terms of our model.