Studying millisecond pulsars in X-rays

Millisecond pulsars represent an evolutionarily distinct group among rotation-powered pulsars. Outside the radio band, the soft X-ray range (～0.1–10 keV) is most suitable for studying radiative mechanisms operating in these fascinating objects. X-ray observations revealed diverse properties of emission from millisecond pulsars. For the most of them, the bulk of radiation is of a thermal origin, emitted from small spots (polar caps) on the neutron star surface heated by relativistic particles produced in pulsar acceleration zones. On the other hand, a few other very fast rotating pulsars exhibit almost pure nonthermal emission generated, most probably, in pulsar magnetospheres. There are also examples of nonthermal emission detected from X-ray nebulae powered by millisecond pulsars, as well as from pulsar winds shocked in binary systems with millisecond pulsars as companions. These and other most important results obtained from X-ray observations of millisecond pulsars are reviewed in this paper, as well as results from the search for millisecond pulsations in X-ray flux of the radio-quite neutron star RX J1856.5-3754.     

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

毫秒脉冲星的自转频率非常稳定,提供了一种独立的基于遥远自然天体并能持续数百万乃至数十亿年的时间基准,具有稳定性强、运行时间长、服务范围广等特点.为了减弱毫秒脉冲星计时观测中各种高斯噪声对脉冲星时的影响,研究了一种基于双谱滤波的综合脉冲星时构建算法,处理分析了国际脉冲星计时阵(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及以上时间尺度上优于原子钟稳定性,可用于改善当前原子时的长期稳定性.     

PSR J1453+1902 and the radio luminosities of solitary versus binary millisecond pulsars

We present 3 yr of timing observations for PSR J1453+1902, a 5.79-ms pulsar discovered during a 430-MHz drift-scan survey with the Arecibo telescope. Our observations show that PSR J1453+1902 is solitary and has a proper motion of  8 ±  2  mas yr⁻¹. At the nominal distance of 1.2 kpc estimated from the pulsar's dispersion measure, this corresponds to a transverse speed of  46 ± 11   km s⁻¹  , typical of the millisecond pulsar population. We analyse the current sample of 55 millisecond pulsars in the Galactic disc and revisit the question of whether the luminosities of isolated millisecond pulsars are different from their binary counterparts. We demonstrate that the apparent differences in the luminosity distributions seen in samples selected from 430-MHz surveys can be explained by small-number statistics and observational selection biases. An examination of the sample from 1400-MHz surveys shows no differences in the distributions. The simplest conclusion from the current data is that the spin, kinematic, spatial and luminosity distributions of isolated and binary millisecond pulsars are consistent with a single homogeneous population.     

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 new look at pulsar statistics — Birthrate and evidence for injection

We make a statistical analysis of the periodsP and period-derivativesP of pulsars using a model independent theory of pulsar flow in theP-P diagram. Using the available sample ofP andP values, we estimate the current of pulsars flowing unidirectionally along theP-axis, which is related to the pulsar birthrate. Because of radio luminosity selection effects, the observed pulsar sample is biased towards lowP and highP. We allow for this by weighting each pulsar by a suitable scale factor. We obtain the number of pulsars in our galaxy to be 6.05_−2.80 ^+3.32 × 10⁵ and the birthrate to be 0.048_−0.011 ^+0.014 pulsars yr⁻¹ galaxy⁻¹. The quoted errors refer to 95 per cent confidence limits corresponding to fluctuations arising from sampling, but make no allowance for other systematic and random errors which could be substantial. The birthrate estimated here is consistent with the supernova rate. We further conclude that a large majority of pulsars make their first appearance at periods greater than 0.5 s. This ‘injection’, which runs counter to present thinking, is probably connected with the physics of pulsar radio emission. Using a variant of our theory, where we compute the current as a function of pulsar ‘age’ (1/2P/P), we find support for the dipole braking model of pulsar evolution upto 6 × 10⁶ yr of age. We estimate the mean pulsar braking index to be 3.7_−0.8 ^+0.8.     

Scintillation parameters for 49 pulsars

We report on multi-epoch, multifrequency observations of 64 pulsars with high spectral and time resolution. Scintillation parameters were obtained for 49 pulsars, including 13 millisecond pulsars. Scintillation speeds were derived for all 49, which doubles the number of pulsars with speeds measured in this way. There is excellent agreement between the scintillation speed and proper motion for the millisecond pulsars in our sample using the simple assumption of a mid-placed scattering screen. This indicates that the scaleheight of scattering electrons is similar to that of the dispersing electrons. In addition, we present observations of the Vela pulsar at 14 and 23 GHz, and show that the scintillation bandwidth scales as ν^3.93 over a factor of 100 in observing frequency. We show that for PSR J0742−2822, and perhaps PSR J0837−4135, the Gum nebula is responsible for the high level of turbulence along their lines of sight, contrary to previous indications. There is a significant correlation between the scintillation speeds and the product of the pulsar's period and period derivative for the 'normal' pulsars. However, we believe this to be caused by selection effects both in pulsar detection experiments and in the choice of pulsars used in scintillation studies.     

The gamma-ray conversion efficiency of rotation-powered pulsars

We use a new self-consistent model to derive the conversion efficiency from rotation power to γ-ray power for pulsars (ε_th). Our result indicates that ε_th∝τ^6/7 P ², where τ and P are the characteristic age and period of the pulsar, which shows that although the efficiency increases with the characteristic age of the pulsar, it also depends on the pulsar period. We test our model results with the survey of high-energy γ-rays from pulsars by EGRET. Our model not only successfully explains the efficiency of the confirmed γ-ray pulsars but also explains why the γ-ray efficiency of millisecond pulsars is so low.     

Radio pulsars in Terzan 5

We report on searches of the globular cluster Terzan 5 for low-luminosity and accelerated radio pulsars using the 64-m Parkes radio telescope. One new millisecond pulsar, designated PSR J1748−2446C, was discovered, having a period of 8.44 ms. Timing measurements using the 76-m Lovell radio telescope at Jodrell Bank show that it is a solitary pulsar and lies close to the core of the cluster. We also present the results of timing measurements which show that the longer period pulsar PSR J1748−2444 (formerly known as PSR B1744−24B) lies 10 arcmin from the core of the cluster and is unlikely to be associated with the cluster. We conclude that there are further pulsars to be detected in the cluster.     

Space distribution and birth-rate of pulsars

Using the data of the 321 pulsars so far known the Galactic distribution and the luminosity function of pulsars have been investigated. The total number of pulsars in our Galaxy is found to be 9 × 10⁴. If the mean age of pulsars is 1.8 × 10⁶ years, the birth-rate of pulsars in the Galaxy will be one every 20 years. This rate is not in contradiction with the birth-rate of supernovae.     

Millisecond pulsar radiation properties

Two investigations of millisecond pulsar radiation are discussed: average total intensity pulse morphology and individual pulse to pulse fluctuations. The average emission profiles of millisecond pulsars are compared with those of slower pulsars in the context of polar cap models. In general the full widths of pulsar emission regions continue to widen inversely with periodP as P^-(0.30-0.5) as expected for dipole polar cap models. Many pulse components are very narrow. The period scaling of pulsar profiles -separations and widths -can tell us about the angular distribution of radiating currents. An investigation of individual pulses from two millisecond pulsars at 430 MHz shows erratic pulse to pulse variations similar to that seen in slow pulsars. PSR B1937+21 displays occasional strong pulses that are located in the trailing edge of the average profile with relative flux densities in the range of 100 to 400. These are similar to the giant pulses seen in the Crab pulsar.     

Can the known millisecond pulsars help in the detection of intermediate-mass black holes at the centers of globular clusters?

We consider the possibility of detecting intermediate-mass (10³–10⁴ M _⊙) black holes, whose existence at the centers of globular clusters is expected from optical and infrared observations, using precise pulse arrival timing for the millisecond pulsars in globular clusters known to date. For some of these pulsars closest to the cluster centers, we have calculated the expected delay times of pulses as they pass in the gravitational field of the central black hole. The detection of such a time delay by currently available instruments for the known pulsars is shown to be impossible at a black hole mass of 10³ M _⊙ and very problematic at a black hole mass of 10⁴ M _⊙. In addition, the signal delay will have a negligible effect on the pulsar periods and their first derivatives compared to the current accuracy of their measurements.     

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.     

Arrival-time analysis for a millisecond pulsar

Arrival times from a fast, quiet pulsar can be used to obtain accurate determinations of pulsar parameters. In the case of the millisecond pulsar, PSR 1937 + 214, the remarkably small rms residual to the timing fit indicates that precise measurements of position, proper motion and perhaps even trigonometric parallax will be possible (Backer 1984). The variances in these parameters, however, will depend strongly on the nature of the underlying noise spectrum. We demonstrate that for very red spectrai.e. those dominated by low-frequency noise, the uncertainties can be larger than the present estimates (based on a white-noise model) and can even grow with the observation period. The possibility of improved parameter estimation through prewhitening’ the data and the application of these results to other pulsar observations are briefly discussed. The post-fit rms residual of PSR 1937 + 214 may be used to limit the energy density of a gravitational radiation background at periods of a few months to years. However, fitting the pulsar position and pulse-emission times filters out significant amounts of residual power, especially for observation periods of less than three years. Consequently the present upper bound on the energy density of gravitational waves Ω_g <3 × 10^-4 R _Μs ² , though already more stringent than any other available, is not as restrictive as had been previously estimated. The present limit is insufficient to exclude scenarios which use primordial cosmic strings for galaxy formation, but should improve rapidly with time. On leave from Raman Research Institute, Bangalore 560080, India.     

Precision timing measurements of PSR J1012+5307

We present results and applications of high-precision timing measurements of the binary millisecond pulsar J1012+5307. Combining our radio timing measurements with results based on optical observations, we derive complete 3D velocity information for this system. Correcting for Doppler effects, we derive the intrinsic spin parameters of this pulsar and a characteristic age of 8.6±1.9 Gyr . Our upper limit for the orbital eccentricity of only 8×10⁻⁷ (68 per cent confidence level) is the smallest ever measured for a binary system. We demonstrate that this makes the pulsar an ideal laboratory in which to test certain aspects of alternative theories of gravitation. Our precision measurements suggest deviations from a simple pulsar spin-down timing model, which are consistent with timing noise and the extrapolation of the known behaviour of slowly rotating pulsars.     

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.     

Morphology and characteristics of radio pulsars

This review describes the observational properties of radio pulsars, fast rotating neutron stars, emitting radio waves. After the introduction we give a list of milestones in pulsar research. The following chapters concentrate on pulsar morphology: the characteristic pulsar parameters such as pulse shape, pulsar spectrum, polarization and time dependence. We give information on the evolution of pulsars with frequency since this has a direct connection with the emission heights, as postulated in the radius to frequency mapping (RFM) concept. We deal successively with the properties of normal (slow) pulsars and of millisecond (fast-recycled) pulsars. The final chapters give the distribution characteristics of the presently catalogued 1300 objects.Received: 5 December 2003, Published online: 15 April 2004 Correspondence to: Richard Wielebinski     

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.     

毫秒脉冲星定时研究进展

毫秒脉冲星守时的理论和方法研究己取得重要进展,利用现有2颗毫秒脉冲星约10yr的计时观测资料分析得到的 TAI-PT,其长期稳定度为 2×10－14。采用合适的长期稳定度算法,由多颗毫秒脉冲星计时观测可以建立综合脉冲星时间尺度。它可以成为与原子时系统比较的重要手段,并对原子时长期稳定度的改进做出贡献。介绍了该领域研究的基本状况,重点对毫秒脉冲星守时的理论方法,综合脉冲星时间及与原子时的关系等进行了讨论和评述。对由双星系统内毫秒脉冲星的轨道运动定义的双星脉冲星时也做了介绍。     

The Parkes Southern Pulsar Survey — III. Timing of long-period pulsars

The Parkes survey of the entire southern sky for millisecond and other low-luminosity pulsars has now been completed. The survey system and initial results were described by Manchester et al. and the final results and population studies are described by Lyne et al. A total of 298 pulsars were detected, including 101 new discoveries of which 17 were millisecond pulsars. Here we report on timing observations at Parkes and Jodrell Bank of the 84 long-period (non-millisecond) pulsars discovered in the survey, including timing solutions for 78 of them. Pulse profiles and arrival times were obtained at several different frequencies over intervals of more than a year, yielding a position, period, period derivative and dispersion measure for each pulsar. Pulse profiles at frequencies near 400 MHz and 600 or 1400 MHz are presented for most of the observed pulsars. Significant timing noise was detected for five pulsars and a glitch was observed in the period of PSR J1123–6259.     

Populating the Galaxy with pulsars – I. Stellar and binary evolution

The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad hoc manner. Thus, no complete and direct comparison with observations of the pulsar population within the Galactic disc has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed typical millisecond pulsar formation process. Results also focus on reproducing the observed     diagram for Galactic pulsars and how this precludes short time-scales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass.