Application of ‘CLEAN’ in the power spectral analysis of non-uniformly sampled pulsar timing data

Spectral analysis of the residual pulsearrival times of pulsars is a useful tool in understanding the nature of the underlying processes that may be responsible for the timing noise observed from pulsars. Power spectra of pulsar timing residuals may be described by one or a combination of powerlaws. As these spectra are expected to be very steep, it is important to ensure a high dynamic range in the estimation of the spectrum. This is difficult in practice since one is, in general, dealing with timing measurements made at unevenly placed epochs. In this paper, we present a technique based on, ‘CLEAN’ to obtain high dynamic range spectra from unevenly sampled data. We compare the performance of this technique with other techniques including some that were used earlier for estimation of power spectra of pulsar timing residuals.     

Rotational irregularities in pulsars — A review

Whilst considerable data have been gathered on the rotational behaviour of pulsars since their discovery, the mechanisms behind the rotational irregularities observed in these objects are still not well understood. This paper describes the observational aspects of these irregularities, namely glitches and timing noise. In addition, it reviews and assesses the theoretical models that have been proposed to explain the irregularities.     

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.     

Methods of searching for planets around pulsars

We analyse different methods of searching for planets around neutron stars by timing observations of pulsars. To this end, we study a few interesting models describing TOA residual variations that are observed, or could be observed, and which can mimic planets. We carry out a detailed theoretical analysis of the behaviour of these methods in the situations mentioned. We show that it is very helpful to look at these phenomena as some kind of quasi-periodic variations of residuals of time of arrival of pulsar pulses. We demonstrate that such a model-independent approach leads to promising conclusions that can be useful when analysing timing observations of pulsars to find planets or to prove that observed phenomena are of planetary origin.     

On measuring the gravitational-wave background using Pulsar Timing Arrays

The long-term precise timing of Galactic millisecond pulsars holds great promise for measuring the long-period (months to years) astrophysical gravitational waves. Several gravitational-wave observational programs, called Pulsar Timing Arrays (PTA), are being pursued around the world.
Here, we develop a Bayesian algorithm for measuring the stochastic gravitational-wave background (GWB) from the PTA data. Our algorithm has several strengths: (i) it analyses the data without any loss of information; (ii) it trivially removes systematic errors of known functional form, including quadratic pulsar spin-down, annual modulations and jumps due to a change of equipment; (iii) it measures simultaneously both the amplitude and the slope of the GWB spectrum and (iv) it can deal with unevenly sampled data and coloured pulsar noise spectra. We sample the likelihood function using Markov Chain Monte Carlo simulations. We extensively test our approach on mock PTA data sets and find that the algorithm has significant benefits over currently proposed counterparts. We show the importance of characterizing all red noise components in pulsar timing noise by demonstrating that the presence of a red component would significantly hinder the detection of the GWB.
Lastly, we explore the dependence of the signal-to-noise ratio on the duration of the experiment, number of monitored pulsars and the magnitude of the pulsar timing noise. These parameter studies will help formulate observing strategies for the PTA experiments.     

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.     

Evidence for a binary companion to the central compact object 1E 1207.4-5209

Unique among neutron stars, 1E 1207.4-5209 is an X-ray pulsar with a spin period of 424 ms that contains at least two strong absorption features in its energy spectrum. This neutron star is positionally coincident with the supernova remnant PKS 1209-51/52 and has been identified as a member of the growing class of radio-quiet compact central objects in supernova remnants. From previous observations with Chandra and XMM-Newton, it has been found that the 1E 1207.4-5209 is not spinning down monotonically as is common for young, isolated pulsars. The spin frequency history requires either strong, frequent glitches, the presence of a fall-back disk, or a binary companion. Here, we report on a sequence of seven XMM-Newton observations of 1E 1207.4-5209 performed during a 40 day window between 2005 June 22 and July 31. Due to unanticipated variance in the phase measurements during the observation period that was beyond the statistical uncertainties, we could not identify a unique phase-coherent timing solution. The three most probable timing solutions give frequency time derivatives of +0.9, ?2.6, and +1.6×10^?12 Hz s^?1 (listed in descending order of significance). We conclude that the local frequency derivative during our XMM-Newton observing campaign differs from the long-term spin-down rate by more than an order of magnitude. This measurement effectively rules out glitch models for 1E 1207.4-5209. If the long-term spin frequency variations are caused by timing noise, the strength of the timing noise in 1E 1207.4-5209 is much stronger than in other pulsars with similar period derivatives. Therefore, it is highly unlikely that the spin variations are caused by the same physical process that causes timing noise in other isolated pulsars. The most plausible scenario for the observed spin irregularities is the presence of a binary companion to 1E 1207.4-5209. We identified a family of orbital solutions that are consistent with our phase-connected timing solution, archival frequency measurements, and constraints on the companions mass imposed by deep IR and optical observations.     

Statistical analyses for NANOGrav 5-year timing residuals

In pulsar timing, timing residuals are the differences between the observed times of arrival and predictions from the timing model. A comprehensive timing model will produce featureless residuals, which are presumably composed of dominating noise and weak physical effects excluded from the timing model(e.g. gravitational waves). In order to apply optimal statistical methods for detecting weak gravitational wave signals, we need to know the statistical properties of noise components in the residuals. In this paper we utilize a variety of non-parametric statistical tests to analyze the whiteness and Gaussianity of the North American Nanohertz Observatory for Gravitational Waves(NANOGrav) 5-year timing data, which are obtained from Arecibo Observatory and Green Bank Telescope from 2005 to2010. We find that most of the data are consistent with white noise; many data deviate from Gaussianity at different levels, nevertheless, removing outliers in some pulsars will mitigate the deviations.     

Scintillation Observations of Strong Northern Pulsars

Scintillation of pulsar radio emission provides information about the interstellar medium along the path to the pulsar and the velocities of pulsars. It also affects the precision of pulse timing observations. Using a pulsar timing system developed at the Urumqi Astronomical Observatory25 m telescope, we observed diffractive scintillation dynamic spectra for several strong northern pulsars. This paper introduces the observing system and discusses the observational results. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

The role of FAST in pulsar timing arrays

The Five-hundred-meter Aperture Spherical radio Telescope(FAST) will become one of the world-leading telescopes for pulsar timing array(PTA) research. The primary goals for PTAs are to detect(and subsequently study) ultra-low-frequency gravitational waves, to develop a pulsar-based time standard and to improve solar system planetary ephemerides. FAST will have the sensitivity to observe known pulsars with significantly improved signal-to-noise ratios and will discover a large number of currently unknown pulsars. We describe how FAST will contribute to PTA research and show that jitter-and timing-noise will be the limiting noise processes for FAST data sets. Jitter noise will limit the timing precision achievable over data spans of a few years while timing noise will limit the precision achievable over many years.     

Low-frequency flux-density measurements and spectral features of millisecond pulsars

Radio flux-density measurements for a large sample of millisecond pulsars at a low frequency of 102 MHz are presented. Using higher frequency measurements, we construct their spectra in the frequency range from 102 MHz to 4.8 GHz, the widest one studied to date. The spectra of millisecond and normal pulsars have been found to differ. The spectra of millisecond pulsars have no low-frequency turnover typical of normal pulsars. The absence of a low-frequency turnover in the spectrum suggests that the emitting regions of millisecond and normal pulsars differ in geometry, which we interpret by deviation of the magnetic field from a dipole one or by compactness of the emitting region.     

射电脉冲星观测资料及其启示

在本文中先简要介绍了８０年代以来射电脉冲星方面的观测进展，然后详细介绍了射电脉冲星的第一手观测资料，包括空间分布，传播特性，时间特性，脉冲轮廓与偏振以及频谱特性等，并详细综述观测资料对我们提供的信息及给我们的启示，为理论模型的建立和检验提供了基础。     

Millisecond pulsar velocities

We present improved timing parameters for 13 millisecond pulsars (MSPs), including nine new proper motion measurements. These new proper motions bring to 23 the number of MSPs with measured transverse velocities. In light of these new results we present and compare the kinematic properties of MSPs with those of ordinary pulsars. The mean transverse velocity of MSPs was found to be 85±13 km s⁻¹, a value consistent with most models for the origin and evolution of MSPs and approximately a factor of 4 lower than that of ordinary pulsars. We also find that, in contrast to young ordinary pulsars, the vast majority of which are moving away from the Galactic plane, almost half of the MSPs are moving towards the plane. This near-isotropy would be expected of a population that has reached dynamic equilibrium. Accurate measurements of MSP velocities have allowed us to correct their measured spin-down rates for Doppler acceleration effects, and thereby derive their intrinsic magnetic field strengths and characteristic ages. We find that close to half of our sample of MSPs have a characteristic age comparable to or greater than the age of the Galactic disc.     

Dispersion measure variations and their effect on precision pulsar timing

We present an analysis of the variations seen in the dispersion measures (DMs) of 20-ms pulsars observed as part of the Parkes Pulsar Timing Array project. We carry out a statistically rigorous structure function analysis for each pulsar and show that the variations seen for most pulsars are consistent with those expected for an interstellar medium characterized by a Kolmogorov turbulence spectrum. The structure functions for PSRs J1045−4509 and J1909−3744 provide the first clear evidence for a large inner scale, possibly due to ion–neutral damping. We also show the effect of the solar wind on the DMs and show that the simple models presently implemented into pulsar timing packages cannot reliably correct for this effect. For the first time we clearly show how DM variations affect pulsar timing residuals and how they can be corrected in order to obtain the highest possible timing precision. Even with our presently limited data span, the residuals (and all parameters derived from the timing) for six of our pulsars have been significantly improved by correcting for the DM variations.     

毫秒脉冲星定时研究进展

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

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.     

Fe emission and ionized excess absorption in the luminous quasar 3C 109 with XMM–Newton

The Parkes High-Latitude pulsar survey covers a region of the sky enclosed by Galactic longitudes 220° < l < 260° and Galactic latitudes | b | < 60°. The observations have been performed using the 20-cm multibeam receiver on the Parkes 64-m radio telescope. A total of 6456 pointings of 265 s each have been collected. The system adopted provided a sensitivity limit, for long-period pulsars with 5 per cent duty cycles, of ∼0.5 mJy. Data analysis resulted in the detection of 42 pulsars of which 18 were new discoveries. Four of these belong to the class of the millisecond – or recycled – pulsars; three of these four are in binary systems. The double pulsar system J0737−3039 is among those and has been presented elsewhere. Here, we discuss the other discoveries and provide timing parameters for the objects for which we have a phase-connected solution.     

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.     

Research on Ensemble Pulsar Time Based on Observed Data

The applications of pulsar timing to several important aspects are introduced and a skeletal definition of the pulsar time scale is given. By taking advantage of the observational timing data of millisecond pulsars obtained at the Parkes Astronomical Observatory, Australia, the ensemble pulsar time based on 4 millisecond pulsars is built up, the comparison of its stability σz, with that of the atomic time is made, and finally a few of important factors which affect the ensemble pulsar time, as well as their influences on the pulsar timing accuracy and application are analyzed.