A Comparison of the Four Algorithms of Ensemble Pulsar Time

The pulsar time defined by a single pulsar is affected by many kinds of noise sources. Its short-term and long-term degrees of stability are both not good enough. In order to weaken the influence of these noise sources on the single pulsar time, an appropriate algorithm can be adopted to make a synthesis of many single pulsar times, then the ensemble pulsar time is obtained, thereby increasing the long-term degree of stability of the ensemble pulsar time. In this article four kinds of algorithms of the ensemble pulsar time are introduced, i.e., the classical weighting algorithm, wavelet analysis algorithm, Wiener filtering algorithm and Wiener filtering algorithm in wavelet domain. These four algorithms are respectively applied to the timing residuals obtained from the observation of two millisecond pulsars, PSR B1855+09 and PSR B1937+21 made at the Arecibo Astronomical Observatory, and comparisons are carried out.     

四种综合脉冲星时算法比较

由单颗脉冲星定义的脉冲星时受多种噪声源的影响,其短期和长期稳定度都不够好.为了削弱这些噪声源对单脉冲星时的影响,可以采取合适的算法对多个单脉冲星时进行综合得到综合脉冲星时,从而提高综合脉冲星时的长期稳定度.文中介绍4种综合脉冲星时算法:经典加权算法、小波分析算法、维纳滤波算法和小波域中的维纳滤波算法,将这4种算法分别应用于Arecibo天文台对两颗毫秒脉冲星PSR B1855+09和PSRB1937+21观测得到的计时残差并作出比较.     

综合脉冲星时的小波分析算法

由单脉冲星定义的脉冲星时受几种噪声源的影响,为削弱这些影响得到一种更稳定的时间尺度,采取将多颗脉冲星定义的单脉冲星时进行综合的分析方法得到综合脉冲星时．用两种方法:经典加权算法和小波分解算法,对PSRB1855 09和PSR B1937 21两颗脉冲星进行综合计算并做出比较,经典加权算法无法兼顾脉冲星不同频率上的稳定度,而用小波分析的方法对PSR B1855 09和PSR B1937 21两颗脉冲星进行综合,能够有效抑制不同频率上噪声的影响,达到更好的效果．     

基于双谱滤波的综合脉冲星时算法研究

毫秒脉冲星的自转频率非常稳定,提供了一种独立的基于遥远自然天体并能持续数百万乃至数十亿年的时间基准,具有稳定性强、运行时间长、服务范围广等特点.为了减弱毫秒脉冲星计时观测中各种高斯噪声对脉冲星时的影响,研究了一种基于双谱滤波的综合脉冲星时构建算法,处理分析了国际脉冲星计时阵(International Pul-sar Timing Array,IPTA)最新发布的4颗毫秒脉冲星(PSR J0437-4715、J0613-0200、J1713+0747和J1909-3744)的观测数据,分析了不同时间尺度综合脉冲星时的稳定性,并与构成国际原子时(International Atomic Time,TAI)的4家授时单位原子钟稳定性进行了比较.结果表明:双谱滤波算法能够较好地抑制观测噪声,提高综合脉冲星时的稳定性.相比于经典加权算法,综合脉冲星时1 yr、10 yr稳定度从7.77×10^-14、8.56×10^-16分别提高到1.50×10^-14、3.50×10^-16,单脉冲星时稳定性的提升也类似.同时发现,综合脉冲星时稳定性在5 yr及以上时间尺度上优于原子钟稳定性,可用于改善当前原子时的长期稳定性.     

Detection of giant radio pulses from the pulsar PSR B0656+14

Giant pulses (GPs) have been detected from the pulsar PSR B0656 + 14. A pulse that is more intense than the average pulse by a factor of 120 is encountered approximately once in 3000 observed periods of the pulsar. The peak flux density of the strongest pulse, 120 Jy, is a factor of 630 higher than that of the average pulse. The GP energy exceeds the energy of the average pulse by up to a factor of 110, which is comparable to that for other known pulsars with GPs, including the Crab pulsar and the millisecond pulsar PSR B1937+21. The giant pulses are a factor of ～6 narrower than the average pulse and are clustered at the head of the average pulse. PSR B0656+14 along with PSR B0031-07, PSR B1112+50, and PSR J1752+2359 belong to a distinctive group of pulsars in which GPs have been detected without any extremely strong magnetic field on the light cylinder.     

A possible association of a long-period pulsar with a supernova remnant?

In this paper we argue that the PSR B 1849+00 and SNR G33.6+0.1 are possibly associated, because (1) there is only a small probability of chance superposition in this sky region; (2) although the estimates of the distance and the age of the two objects contain a high degree of uncertainty, they are not inconsistent; (3) the implied transverse velocity of the pulsar is also quite reasonable according to the measured distribution of pulsar velocities. The association could be verified by measuring the proper motion of PSR B 1849+00, which is estimated to be about 10 mas·yr^–1.     

Searches for Planets Around Neutron Stars

We review the methodology of searches for planet-mass bodies around neutron stars observable as radio pulsars and discuss the current status of these searches. PSR B1257 + 12, a 6.2-millisecond pulsar, remains the only neutron star accompanied by confirmed planets. It is possible that there is a fourth distant planet in the 1257+12 system. The best of the other candidates for pulsar planets under consideration is a distant, possibly Jovian-mass companion to PSR B1620-26, a 11-millisecond pulsar in the globular cluster M4. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Stokes phase portraits of descattered pulse profiles of a few pulsars

The observing signals from pulsar are always influenced by the interstellar medium (ISM) scattering. In the lower frequency observation, the intensity profiles are broadened and the plane of polarization angle (PPA) curves are flattened by the scattering effect of the ISM. So before we analyze the scattered signal, we should take a proper approach to clear scattering effect from it. Observing data and simulation have shown that the Stokes phase portraits I–U, I–Q and Q–U are also distorted by the ISM scattering. In this paper, a simulation is held to demonstrate a scattering and a descattering of the Stokes phase portraits of a single pulse profile of a pulsar. As a realization of the simulation method, this paper has studied the descattering of Stokes phase portraits of lower frequency observation of PSR B1356?60, PSR B1831?03, PSR B1859+03, PSR B1946+35.     

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.     

The effect of pulse profile evolution on pulsar dispersion measure

In an earlier paper, based on simultaneous multifrequency observations with the Giant Metrewave Radio Telescope (GMRT), we reported the variation of pulsar dispersion measures (DMs) with frequency. A few different explanations are possible for such frequency dependence, and a possible candidate is the effect of pulse shape evolution on the DM estimation technique. In this paper we describe extensive simulations we have done to investigate the effect of pulse profile evolution on pulsar DM estimates. We find that it is only for asymmetric pulse shapes that the DM estimate is significantly affected due to profile evolution with frequency. Using multifrequency data sets from our earlier observations, we have carried out systematic analyses of PSR B0329+54 and PSR B1642−03. Both these pulsars have central core-dominated emission which does not show significant asymmetric profile evolution with frequency. Even so, we find that the estimated DM shows significant variation with frequency for these pulsars. We also report results from new, simultaneous multifrequency observations of PSR B1133+16 carried out using the GMRT in phased array mode. This pulsar has an asymmetric pulse profile with significant evolution with frequency. We show that in such a case, amplitude of the observed DM variations can be attributed to profile evolution with frequency. We suggest that genuine DM variations with frequency could arise due to propagation effects through the interstellar medium and/or the pulsar magnetosphere.     

Optimal filters for the construction of the ensemble pulsar time

An algorithm of the ensemble pulsar time based on the optimal Wiener filtration method has been constructed. This algorithm allows the separation of the contributions to the post-fit pulsar timing residuals of the atomic clock and the pulsar itself. Filters were designed using the cross- and auto-covariance functions of the timing residuals. The method has been applied to the timing data of millisecond pulsars PSR B1855+09 and B1937+21 and allowed the filtering out of the atomic-scale component from the pulsar data. Direct comparison of the terrestrial time TT(BIPM06) and the ensemble pulsar time PT_ens revealed that the fractional instability of TT(BIPM06)−PT_ens is equal to  σ _z = (0.8 ± 1.9) × 10⁻¹⁵  . Based on the  σ _z   statistics of TT(BIPM06)−PT_ens, a new limit of the energy density of the gravitational wave background was calculated to be equal to  Ω_g h ²∼ 3 × 10⁻⁹  .     

Hyperfast pulsars as the remnants of massive stars ejected from young star clusters

Recent proper motion and parallax measurements for the pulsar PSR B1508+55 indicate a transverse velocity of  ∼1100 km s⁻¹  , which exceeds earlier measurements for any neutron star. The spin-down characteristics of PSR B1508+55 are typical for a non-recycled pulsar, which implies that the velocity of the pulsar cannot have originated from the second supernova disruption of a massive binary system. The high velocity of PSR B1508+55 can be accounted for by assuming that it received a kick at birth or that the neutron star was accelerated after its formation in the supernova explosion. We propose an explanation for the origin of hyperfast neutron stars based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star which attained its peculiar velocity (similar to that of the pulsar) in the course of a strong dynamical three- or four-body encounter in the core of dense young star cluster. To check this hypothesis, we investigated three dynamical processes involving close encounters between: (i) two hard massive binaries, (ii) a hard binary and an intermediate-mass black hole (IMBH) and (iii) a single stars and a hard binary IMBH. We find that main-sequence O-type stars cannot be ejected from young massive star clusters with peculiar velocities high enough to explain the origin of hyperfast neutron stars, but lower mass main-sequence stars or the stripped helium cores of massive stars could be accelerated to hypervelocities. Our explanation for the origin of hyperfast pulsars requires a very dense stellar environment of the order of  10⁶– 10⁷ stars pc⁻³  . Although such high densities may exist during the core collapse of young massive star clusters, we caution that they have never been observed.     

The possible companions of young radio pulsars

We discuss the formation of pulsars with massive companions in eccentric orbits. We demonstrate that the probability for a non-recycled radio pulsar to have a white dwarf as a companion is comparable to that of having an old neutron star as a companion. Special emphasis is given to PSR B1820−11 and PSR B2303+46. Based on population synthesis calculations we argue that PSR B1820−11 and PSR B2303+46 could very well be accompanied by white dwarfs with mass ≳1.1 M_⊙. For PSR B1820−11, however, we cannot exclude the possibility that its companion is a main-sequence star with a mass between ∼0.7 M_⊙ and ∼5 M_⊙.     

The micro-glitch in PSR B1821−24: a case for a strange pulsar?

The single glitch observed in PSR B1821−24, a millisecond pulsar in M28, is unusual on two counts. First, the magnitude of this glitch is at least an order of magnitude smaller  (Δν/ν∼ 10⁻¹¹)  than the smallest glitch observed to date. Secondly, all other glitching pulsars have strong magnetic fields with   B ≳ 10¹¹ G  and are young, whereas PSR B1821−24 is an old recycled pulsar with a field strength of  2.25 × 10⁹ G  . We have earlier suggested that some of the recycled pulsars could actually be strange quark stars. In this work, we argue that the crustal properties of such a strange pulsar are just right to give rise to a glitch of this magnitude, explaining the scarcity of larger glitches in millisecond pulsars.     

Constraint on Parameters of Inverse Compton Scattering Model for PSR B2319+60

Using the multifrequency radio profiles of pulsar PSR B2319+60, two parameters of inverse Compton scattering model, the initial Lorentz factor and the factor of energy loss of relativistic particles are constrained.     

The Radio and X-ray Mode-Switching Pulsar PSR B0943+10

Observations obtained in the last years challenged the widespread notion that rotation-powered neutron stars are steady X-ray emitters. Besides a few allegedly rotation-powered neutron stars that showed ‘magnetar-like’ variability, a particularly interesting case is that of PSR B0943+10. Recent observations have shown that this pulsar, well studied in the radio band where it alternates between a bright and a quiescent mode, displays significant X-ray variations, anticorrelated in flux with the radio emission. The study of such synchronous radio/X-ray mode switching opens a new window to investigate the processes responsible for the pulsar radio and high-energy emission. Here we review the main X-ray properties of PSR B0943+10 derived from recent coordinated X-ray and radio observations.     

环加速间隙下的脉冲星星风制动模型

为了解释间歇脉冲星PSR B1931+24在射电噪比射电宁静状态下更大的自转减慢率和模拟蟹状星云脉冲星的自转演化,建立同时考虑了具有不同加速电势的核区和环区的环加速间隙下的星风制动模型.其中对于PSR B1931+24通过计算得到它的磁场强度和磁倾角,并且预言了其理论制动指数.对于蟹状星云脉冲星,通过计算得到它的磁场强度和磁倾角,还计算得到其制动指数随周期的演化和它在周期-周期导数图上的自转演化.相比于真空加速间隙、外加速间隙等,环加速间隙也同样能够适用于星风制动模型.     

A search for giant pulses in Vela-like pulsars

We have carried out a survey for 'giant pulses' in six young, Vela-like pulsars. In no cases did we find single pulses with flux densities more than 10 times the mean flux density. However, in PSR  B1706–44  we have detected giant micro-pulses very similar to those seen in the Vela pulsar. In PSR  B1706–44  these giant micro-pulses appear on the trailing edge of the profile and have an intrinsic width of ∼1 ms. The cumulative probability distribution of their intensities is best described by a power law. If the power law continues to higher intensities, then  3.7×10⁶  rotations are required to obtain a pulse with 20× the mean pulse flux. This number is similar to the giant pulse rate in PSR B1937+21 and PSR  B1821–24  but significantly higher than that for the Crab.     

Long-term scintillation observations of five pulsars at 1540 MHz

From 2001 January to 2002 June, we monitored PSRs B0329+54, B0823+26, B1929+10, B2020+28 and B2021+51 using the Nanshan 25-m radio telescope of the Urumqi Observatory to study their diffractive interstellar scintillation (DISS). The average interval between observations was about 9 d and the observation duration ranged between 2 and 6 h depending on the pulsar. Wide variations in the DISS parameters were observed over the 18-month data span. Despite this, the average scintillation velocities are in excellent agreement with the proper motion velocities. The average two-dimensional autocorrelation function for PSR B0329+54 is well described by a thin-screen Kolmogorov model, at least along the time and frequency axes. Observed modulation indices for the DISS time-scale and bandwidth and the pulsar flux density are greater than values predicted for a Kolmogorov spectrum of electron density fluctuations. Correlated variations over times that are long compared to the nominal refractive scintillation time are observed, suggesting that larger scale density fluctuations are important. For these pulsars, the scintillation bandwidth as a function of frequency has a power-law index  (∼3.6)  much less than that expected for Kolmogorov turbulence (∼4.4). Sloping fringes are commonly observed in the dynamic spectra, especially for PSR B0329+54. The detected range of fringe slopes are limited by our observing resolution. Our observations are sensitive to larger-scale fringes and hence smaller refractive angles, corresponding to the central part of the scattering disc.     

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.