A Statistical Analysis of Radio Pulsar Timing Noise

We present an analysis of the timing observations on 27 radio pulsars, collected at Hartebeesthoek Radio Astronomy Observatory (HartRAO), with time spans ranging between ～ 9 and 14yr. Our results show that the measured pulsar frequency second derivatives are non-stationary. Both the magnitude and the sign of the ■ values depend upon the choice of epoch and data span. A simple statistical analysis of the observed second time derivative of the pulse frequency (■obs) of a large sample of 391 (25 HartRAO and 366 Jodrell Bank Observatory) pulsars reveals that ■ is only marginally correlated with both the pulsar spin-d own rate ( ■) and the characteristic age (τc). We find correlation coefficients of ～ 0.20a nd-0.30 between the measured braking indices and, respectively, ■ and τc. This result reaffirms earlier conclusions that the braking indices of most radio pulsars, obtained through the standard timing technique, are strongly dominated by sustained random fluctuations in the observed pulse phase.     

脉冲星射电轮廓中高斯成分在主脉冲和中间脉冲的不同分布

研究了高斯辐射成分在可视点所画出轨迹上的分布,这个可视点因脉冲星的转动而作非匀速度运动.通过假设辐射区域围绕磁轴均匀分布,一个高斯辐射成分便对应于可视轨迹划过的一个辐射区域.因为演示辐射区域在可视轨迹上是不均匀的分布,因此高斯成分沿轨迹也是不均匀的,而高斯成分的密度在磁轴与视线距离最近时为最大.高斯成分的分布取决于脉冲星的两个角度:旋转轴和视线之间的夹角,以及磁轴和旋转轴之间的倾角.基于此模型,一个脉冲星平均轮廓中观察到的多个高斯成分便对应于可视轨迹在特定的转动相位范围内的辐射区域.演示了脉冲星旋转的近侧和远侧的相位,分别对应的主脉冲和中间脉冲,两者高斯成分的数量和分布是不同的.而且还发现,沿可视轨迹上的辐射区域总数与围绕磁轴的辐射区域的总数是不同,并且预测的辐射区域数目会因忽略可见点的运动而明显不同.拟合表明脉冲星轮廓的高斯成分的形状和数量可能与实际构成轮廓的成分的形状和数量不同.以PSR B0826–34的辐射为例,并假设辐射来自单一磁极.     

Pulsar braking indices revisited

Using the standard equation for the slowdown of a neutron star, we derive a formula for the braking index via integration rather than the conventional differentiation. The new formula negates the need to measure the second time derivative of the rotation frequency, ν¨ . We show that the method gives similar braking indices for PSR B1509−58 and the Crab pulsar to those already in the literature. We point out that our method is useful for obtaining the braking indices of moderate-aged pulsars without the need for long, phase-connected timing solutions. We applied the method to 20 pulsars and discuss the implications of the results. We find that virtually all the derived braking indices are dominated by the effects of (unseen) glitches, the recovery from which corrupts the value of ν˙ . However, any real, large, positive braking index has implications for magnetic field decay and offers support to recent models of pulsar evolution.     

The Effect of Lunisolar Precession on Pulsar Observations

In the measurements of pulsars’ periods and their rates of change there exists the effect of lunisolar precession. Via the analysis of this effect, we find that it can affect the accurate measurements of the periods of certain pulsars. An investigation of the periods of 1771 pulsars reveals that the measured periods of 81 pulsars suffer the influence of lunisolar precession in the range of the accuracy of measurements. With the improvement of the precision of measurements of pulsars’ periods and their change rates, this effect should get the corresponding corrections.     

Pulsar radio luminosity laws revisited

A luminosity law for radio pulsars was proposed by Stollman, based on features of the RudermanSutherland model. We analyse that proposition in the light of the most recent data (the catalogue of Taylor, Manchester & Lyne), and perform a set of statistical tests to verify the reliability of the results. We find little evidence favouring Stollman's law, and point out a possible incompleteness of the data sets, without conclusive evidence for or against the RudermanSutherland model.     

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.     

Crab脉冲星巨脉冲辐射等待时间分布的非稳态泊松模拟

利用中国科学院新疆天文台南山观测站26m射电望远镜, 在中心频率1556MHz, 对Crab脉冲星(PSR B0531+21)进行了长达12.6h的连续观测, 观测带宽为512MHz, 时间分辨率为32μs, 研究了巨脉冲辐射的等待时间分布特征. 观测共探测到2097个信噪比大于10的巨脉冲, 对应的流量密度大于100Jy. 巨脉冲的爆发率表现为高度的间歇性, 在较短的时间内具有较高的爆发率, 在相对长的宁静期内巨脉冲的爆发率较低, 尤其是中间脉冲相位内的巨脉冲爆发. 相邻两个巨脉冲的等待时间分布表现为幂律分布特征, 可以用一个非稳态的泊松过程进行模拟, 这表明巨脉冲的爆发是一种独立的随机事件. 此外, 主脉冲和中间脉冲相位上的巨脉冲具有不同的等待时间分布特征, 这意味着脉冲星不同磁极的巨脉冲辐射机制可能是不同的. 这些观测结果对于理解脉冲星的射电辐射机制具有重要意义.     

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.     

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

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

A Large Glitch in the Crab Pulsar

1 INTRODUCTIONA pulsar glitch is a phenomenon in which the pulse frequency has a sudden increase, typicallywith a fractional ~litude an/u in the range of 10--8 -- 10--6. Glitches are Unpredictable butoccur at intervals of a few y6ars in many young pulsars. Coincident with the glitch, there is oftenan increajse in the magnitude of the frequency derivative followed by thlee types of recovery:almost no recoveryt partial recovery and total recovery on timescales from days to years. It isbelie…     

Circular Polarization in Pulsar Integrated Profiles： Updates

We update the systematic studies of circular polarization in integrated pulse profiles by Han et al. Data of circular polarization profiles are compiled. Sense reversals can occur in core or cone components, or near the intersection between components. The correlation between the sense of circular polarization and the sense of position angle variation for conal-double pulsars is confirmed with a much large database. Circular polarization of some pulsars has clear changes with frequency. Circular polarization of millisecond pulsars is marginally different from that of normal pulsars.     

脉冲星数据比对分析和可视化系统设计与实现

脉冲星数据比对分析和可视化系统(PSRDB,URL:http://www.psrdb.net/),由FAST(Five-hundred-meter Aperture Spherical Radio Telescope)早期科学数据中心团队为快速开展脉冲星候选体比对分析和数据管理研发.通过前端数据提交页面,接收和维护来自FAST及其他研究机构的候选体数据.目前,PSRDB已收录自1967年人类发现第1颗脉冲星以来所有公开文献发表的2811颗脉冲星样本,并采集了当前主要巡天项目尚未正式发表的源和候选体,如FAST多科学目标同时扫描巡天(CRAFTS)候选体数据.基于入库基础数据,利用位置、周期、色散等参数进行比对分析,辅助科研工作者在线检索匹配已知星表数据,最后将检索匹配、比对分析结果生成图表供进一步分析.目前,PSRDB已被应用于FAST脉冲星搜寻和候选体数据管理.未来,PSRDB可在新源认证、后随观测、观测计划制定和原始数据处理流程设计等方面提供数据和工具支撑.     

Pulsar Radio and Gamma-Ray Emission

A total of eight gamma-ray sources are identified with pulsars and these include some of the strongest gamma-ray sources in the sky. About 20 of the unidentified gamma-ray sources are very likely to be associated with currently known pulsars and there is little doubt that many of the others, at least those at low Galactic latitudes, will ultimately be identified with pulsars. How many of these and future gamma-ray detections will be detectable at radio wavelengths depends on the details of the radio and gamma-ray beaming. There is good evidence that the radio beams in young and millisecond pulsars are very wide, implying that most gamma-ray pulsars will be detectable in the radio band.     

The Parkes Multibeam Pulsar Survey – VI. Discovery and timing of 142 pulsars and a Galactic population analysis

We present the discovery and follow-up observations of 142 pulsars found in the Parkes 20-cm multibeam pulsar survey of the Galactic plane. These new discoveries bring the total number of pulsars found by the survey to 742. In addition to tabulating spin and astrometric parameters, along with pulse width and flux density information, we present orbital characteristics for 13 binary pulsars which form part of the new sample. Combining these results from another recent Parkes multibeam survey at high Galactic latitudes, we have a sample of 1008 normal pulsars which we use to carry out a determination of their Galactic distribution and birth rate. We infer a total Galactic population of  30 000 ± 1100  potentially detectable pulsars (i.e. those beaming towards us) having 1.4-GHz luminosities above 0.1 mJy kpc². Adopting the Tauris & Manchester beaming model, this translates to a total of  155 000 ± 6000  active radio pulsars in the Galaxy above this luminosity limit. Using a pulsar current analysis, we derive the birth rate of this population to be  1.4 ± 0.2  pulsars per century. An important conclusion from our work is that the inferred radial density function of pulsars depends strongly on the assumed distribution of free electrons in the Galaxy. As a result, any analyses using the most recent electron model of Cordes & Lazio predict a dearth of pulsars in the inner Galaxy. We show that this model can also bias the inferred pulsar scaleheight with respect to the Galactic plane. Combining our results with other Parkes multibeam surveys we find that the population is best described by an exponential distribution with a scaleheight of 330 pc. Surveys underway at Parkes and Arecibo are expected to improve the knowledge of the radial distribution outside the solar circle, and to discover several hundred new pulsars in the inner Galaxy.     

Current Flows in Pulsar Magnetospheres

The global structure of current flows in pulsar magnetosphere is investigated, with rough calculations of the circuit elements. It is emphasized that the potential of the critical field lines (the field lines that intersect the null surface at the light cylinder radius) should be the same as that of interstellar medium, and that pulsars whose rotation axes and magnetic dipole axes are parallel should be positively charged, in order to close the pulsar's current flows. The statistical relation between the radio luminosity and pulsar's electric charge (or the spindown power) may hint that the millisecond pulsars could be low-mass bare strange stars.     

大规模脉冲星候选体信号的无监督聚类分析研究

随着500 m口径球面射电望远镜(Five-hundred-meter Aperture Spherical radio Telescope, FAST)等大型射电望远镜的建设和使用,脉冲星巡天数据进入PB时代.为解决如此大量高速采样的标量数据挖掘问题,促进新天文现象的发现,提出一种基于无监督聚类的脉冲星候选体筛选方案.该方案采用基于密度层次、划分方法的混合聚类算法,结合MapReduce/Spark并行计算模型和基于滑动窗口的分组策略,进而提高大量候选体信号筛选的效率.通过在脉冲星数据集HTRU2 (High Time Resolution Universe)上的对比实验,结果表明该算法能取得较高的精确度和召回率,分别是0.946和0.905,并且当并行节点足够时,该算法的时间复杂度相比串行执行明显下降.可见,该方法为脉冲星观测大数据的分析挖掘提供一种可行思路.