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

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

毫秒脉冲星定时研究进展

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

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

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

基于实测数据的综合脉冲星时研究

介绍了脉冲星计时在几个重要方面的应用,给出了脉冲星时间尺度的简要定义.基于澳大利亚PARKES天文台的毫秒脉冲星实测计时数据,建立了基于4颗毫秒脉冲星的综合脉冲星时,并将其稳定度σz与原子时进行比较,最后分析了影响综合脉冲星时的几个重要因素及其对脉冲星计时精度的影响与应用.     

脉冲星时间特征和多星计时观测

在评介脉冲星时研究发展趋势的基础上,讨论了基于毫秒脉冲星自转的脉冲星时间的短期特征和长期稳定度,并对影响其稳定度的因素作了初步分析.50m射电望远镜计时观测的目的之一就是综合脉冲星时研究.对于50m望远镜计时观测的选星和由多星观测进行综合脉冲星时研究中的数据分析问题进行了讨论.     

毫秒脉冲星计时观测进展

在评述国际上毫秒脉冲星长期计时观测成果的基础上，比较了由毫秒脉冲星自转定义的脉冲星时PT与原子时AT的长期稳定度。毫秒脉冲星PSR B1855 09 15yr的长期稳定度优于原子时TAI—USNO，而接近原子时TAI—PTB；其短期稳定度主要受到计时观测误差的限制。毫秒脉冲星PSR J0437-47 1yr以上的稳定度明显优于上述两原子时系统。毫秒脉冲星PSR B193T 21和PSR J1713 07的长期稳定度受到低频噪声影响，但仍能对综合脉冲星时PTens的建立与保持作出贡献。另外，介绍和评述了国际上脉冲星计时阵的现状及发展。随着脉冲星计时阵的实旅，在脉冲星时间标准研究、原子时误差、太阳系历表误差研究和引力波探测方面将会取得重要成果；展望了国际上平方公里望远镜(SKA)在毫秒脉冲星计时观测方面的发展前景。     

毫秒(ms)脉冲星计时观

简介了国际上脉冲星计时研究的成果。主要评述中国在脉冲星计时研究方面的进展和观点:某些ms脉冲星可以用作天然的高精度的频率标准源。ms脉冲星时(PTens)的长期稳定度随观测时间的加长而提高,同时,它具有寿命长,稳定、可靠、不需要维修等优点。尽早适时地建立专门用于测量ms脉冲星时(PTens)的脉冲星钟系统网络是有益的。我国正在建设用于脉冲星脉冲到达时间(TOA)测量的50m射电望远镜,和正在研究(FAST)500m射电望远镜。论述了脉冲星时间尺度的建立是可行的,脉冲星时和原子时可以"交叉比对,相互检验,并肩发展"。研究表明综合脉冲星时的长期稳定度明显优于原子时。阐述了把小波分析方法应用于时间尺度,研究结果比较传统方法有明显的优势。对国际间的合作提出建议。并展望了中国的脉冲星计时规划。     

毫秒（ms）脉冲星计时观

简介了国际上脉冲星计时研究的成果。主要评述中国在脉冲星计时研究方面的进展和观点：某些ms脉冲星可以用作天然的高精度的频率标准源。ms脉冲星时(PTens)的长期稳定度随观测时间的加长而提高，同时，它具有寿命长，稳定、可靠、不需要维修等优点。尽早适时地建立专门用于测量ms脉冲星时(PFens)的脉冲星钟系统网络是有益的。我国正在建设用于脉冲星脉冲到达时间(TOA)测量的50m射电望远镜，和正在研究(FAST)500m射电望远镜。论述了脉冲星时间尺度的建立是可行的，脉冲星时和原子时可以“交叉比对，相互检验，并肩发展”。研究表明综合脉冲星时的长期稳定度明显优于原子时。阐述了把小波分析方法应用于时间尺度，研究结果比较传统方法有明显的优势。对国际间的合作提出建议。并展望了中国的脉冲星计时规划。     

脉冲星时稳定度及可能应用

给出了毫秒脉冲星长期计时观测最新结果和脉冲星时与原子时频率稳定度的比较。指出脉冲星时频率稳定度的提高受到计时观测误差的限制,讨论了提高毫秒脉冲星计时观测精度的方法。随着脉冲星计时阵的实施,脉冲星时的应用已为期不远。可能的应用包括毫秒脉冲星与原子钟结合守时、建立综合脉冲星时和用脉冲星计时阵检测原子时误差等方面。     

脉冲星对氢原子钟的频率驾驭算法研究

氢原子钟具有较高的短期稳定度,将其作为主钟可在短期内产生高精度的本地时间信号.但氢钟存在频率漂移现象,导致其长期稳定度较差,从而影响本地时间的准确性.毫秒脉冲星自转高度稳定,借助于其长期稳定度高的特性,可定期实现对氢原子钟的频率驾驭,并对实时信号加以控制.首先分析了国际脉冲星计时阵(International Pulsar Timing Array,IPTA)第二批发布数据中四颗毫秒脉冲星的稳定度随时间的变化,同时采用哈达玛方差分析了中国科学院国家授时中心(National Time Service Center,NTSC)一台氢钟的频率稳定性能,最终给出了利用脉冲星驾驭氢原子钟频率的方法.     

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.     

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⁻⁹  .     

Wavelet Analysis Algorithm for Synthetic Pulsar Time

The pulsar time defined by a single pulsar is affected by several noises and in order to weaken their effects and acquire a much more stable time scale we define a synthetic pulsar time from many single pulsar times. Synthesis of the two pulsars, PSR B1855 + 09 and PSR B1937 + 21 is implemented by two methods: the classical weighting algorithm and the wavelet decomposition algorithm. The results are compared. The classical weighting algorithm is unable to take into consideration the different degrees of stability at different frequencies while the wavelet algorithm can, and thereby get better results.     

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.     

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.     

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.