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

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

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

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

Real-time signal processor for pulsar studies

This paper describes the design, tests and preliminary results of a real-time parallel signal processor built to aid a wide variety of pulsar observations. The signal processor reduces the distortions caused by the effects of dispersion, Faraday rotation, doppler acceleration and parallactic angle variations, at a sustained data rate of 32 Msamples/sec. It also folds the pulses coherently over the period and integrates adjacent samples in time and frequency to enhance the signal-to-noise ratio. The resulting data are recorded for further off-line analysis of the characteristics of pulsars and the intervening medium. The signal processing for analysis of pulsar signals is quite complex, imposing the need for a high computational throughput, typically of the order of a Giga operations per second (GOPS). Conventionally, the high computational demand restricts the flexibility to handle only a few types of pulsar observations. This instrument is designed to handle a wide variety of Pulsar observations with the Giant Metre Wave Radio Telescope (GMRT), and is flexible enough to be used in many other high-speed, signal processing applications. The technology used includes field-programmable-gate-array(FPGA) based data/code routing interfaces, PC-AT based control, diagnostics and data acquisition, digital signal processor (DSP) chip based parallel processing nodes and C language based control software and DSP-assembly programs for signal processing. The architecture and the software implementation of the parallel processor are fine-tuned to realize about 60 MOPS per DSP node and a multiple-instruction-multiple-data (MIMD) capability.     

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.     

Radio Phase-Resolved Spectra of the Conal-Double Pulsar B1133 16

Based on dividing the profile into a number of absolute phase intervals,the phase-resolved spectra (PHRS) are derived from published time-aligned average profiles at radio frequencies over two decades for the classic conal-double pulsar B1133 16. The relative spectral index,defined as the difference between the spectral indices of a reference and the given arbitrary phase interval,is obtained by power-law fit at each phase interval. The derived phase-resolved spectra show an "M-like" shape,of which the leading part and trailing part are approximately symmetrical. The basic feature of the PHRS is that the spectrum first flat-tens then steepens as the pulse phase sweeps from the profile center to the profile edges. The PHRS provide a coherent explanation of the major features of profile evolution of B1133 16,namely,the pulse width shrinkage with increasing frequency and the frequency evolution of the relative intensity between the leading and trailing conal components,and the bridge emission. The PHRS may be an indicator for emission spectral variation across the pulsar magnetosphere. Possible mapping from PHRS to emission-location-dependent spectral vari-ation is presented,and some intrinsic mechanisms are discussed.     

The clustering of galaxy clusters in cosmological models with non-Gaussian initial conditions: predictions for future surveys

We predict the biasing and clustering properties of galaxy clusters that are expected to be observed in the catalogues produced by two forthcoming X-ray and Sunyaev–Zel'dovich effect surveys. We study a set of flat cosmological models where the primordial density probability distribution shows deviations from Gaussianity in agreement with current observational bounds form the background radiation. We consider both local and equilateral shapes for the primordial bispectrum in non-Gaussian models. The two catalogues investigated are those produced by the e ROSITA wide survey and from a survey based on South Pole Telescope observations. It turns out that both the bias and observed power spectrum of galaxy clusters are severely affected in non-Gaussian models with local shape of the primordial bispectrum, especially at large scales. On the other hand, models with equilateral shape of the primordial bispectrum show only a mild effect at all scales, that is difficult to be detected with clustering observations. Between the two catalogues, the one performing better is the e ROSITA one, since it contains only the largest masses that are more sensitive to primordial non-Gaussianity.     

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.     

Processing and Analysis of X-ray Timing Data for Crab Pulsartwo

The X-ray timing data for the Crab pulsar obtained by the Chinese X-ray pulsar navigation test satellite are processed and analyzed. The method to build the integrated and standard X-ray pulse profiles of the Crab pulsar by using the X-ray pulsar observation data and the satellite orbit data is described. The principle and algorithm for determining the pulsar's pulse time of arrival (toa) in the frequency domain are briefly introduced. The pulsar's pulse time of arrival is calculated by using the timing data of 50 min integration for each set of observational data. By the comparison between the observed Crab pulsar's pulse time of arrival at the solar system barycenter and that predicted with the Crab pulsar ephemeris, it is found that the timing accuracy is about 14 μs after the systematic error is removed by a quadratic polynomial fitting.     

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

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

Mechanism of generation of the emission bands in the dynamic spectrum of the Crab pulsar

We show that the proportionately spaced emission bands in the dynamic spectrum of the Crab pulsar fit the oscillations of the square of a Bessel function whose argument exceeds its order. This function has already been encountered in the analysis of the emission from a polarization current with a superluminal distribution pattern: a current whose distribution pattern rotates (with an angular frequency ω) and oscillates (with a frequency  Ω > ω  differing from an integral multiple of ω) at the same time. Using the results of our earlier analysis, we find that the dependence on frequency of the spacing and width of the observed emission bands can be quantitatively accounted for by an appropriate choice of the value of the single free parameter  Ω/ω  . In addition, the value of this parameter, thus implied by Hankins & Eilek's data, places the last peak in the amplitude of the oscillating Bessel function in question at a frequency  (∼Ω³/ω²)  that agrees with the position of the observed ultraviolet peak in the spectrum of the Crab pulsar. We also show how the suppression of the emission bands by the interference of the contributions from differing polarizations can account for the differences in the time and frequency signatures of the interpulse and the main pulse in the Crab pulsar. Finally, we put the emission bands in the context of the observed continuum spectrum of the Crab pulsar by fitting this broad-band spectrum (over 16 orders of magnitude of frequency) with that generated by an electric current with a superluminally rotating distribution pattern.     

The possible role of r-modes in post-glitch relaxation of the Crab pulsar

The loss of angular momentum through gravitational radiation, driven by the excitation of r-modes, is considered for neutron stars that have rotation frequencies lower than the associated critical frequency. We find that for reasonable values of the initial amplitudes of such pulsation modes of the star, being excited at the event of a glitch in a pulsar, the total post-glitch losses correspond to a negligible fraction of the initial rise of the spin frequency in the case of Vela and older pulsars. However, for the Crab pulsar the same effect would result, within a few months, in a decrease in its spin frequency by an amount larger than its glitch-induced frequency increase. This could provide an explanation for the peculiar behaviour observed in the post-glitch relaxations of the Crab pulsar.     

面向单脉冲信号分类的集成特征选择与评价

受大量射频干扰信号影响, 快速从海量观测数据中准确识别出单脉冲信号已成为天文数据处理的一项重要任务, 而设计和提取有效数据特征, 是利用机器学习进行单脉冲信号高效识别的决定因素. 针对如何选择最优特征, 进而提升单脉冲信号的分类精度这一关键问题, 设计了面向单脉冲信号分类的集成特征选择方法. 方法首先混合单脉冲信号的参数特征、统计特征和抽象特征, 然后分别利用5种单一特征选择方法选出各自的最优特征集, 最后利用贪心策略对5种单一方法获取的最优特征集进行集成筛选, 获取最优集成特征集. 实验表明, 最优特征集合既包含统计特征也包含抽象特征. 在相同特征数量下, 利用集成特征选择比单一特征选择能获得更高的模型精度, 可使F1值最高提升1.8%. 在海量数据背景下, 集成特征选择对减少特征数量、提升分类性能和加快数据处理速度具有重要作用.     

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.     

A Research on the Pulsar Timing Based on Kalman Filtering

The pulsar is a high-speed rotationing neutron star with a stable rotational period, being not disturbed and destroyed artificially, and can be taken as the reference quantity of the absolute time. In this article a kind of pulsar time service method based on the Kalman filtering algorithm is proposed, and the simulation analysis of the clock error control based on the Kalman filtering and of the effect of the pulsar catalogue error and the measuring accuracy of the pulsar time of arrival (TOA) on the accuracy of time service is made by taking a certain solar synchronous orbit as an example. The result shows that by utilizing this method the clock error of the satellite-borne clock can be effectively eliminated and its time-dependent increase is restrained, thereby solving the problem that the accuracy of the spacecraft-borne low cost clock can not meet the needs.     

Large-scale non-Gaussianities in the 21-cm background anisotropies from the era of reionization

The brightness temperature fluctuations in the 21-cm background related to the neutral hydrogen distribution provide a probe of the physics related to the era of reionization, when the intergalactic medium changed from being completely neutral to partially ionized. We formulate statistics of 21-cm brightness temperature anisotropies in terms of the angular power spectrum, the bispectrum, and the trispectrum. Using the trispectrum, we estimate the covariance related to the power spectrum measurements and show that correlations resulting from non-Gaussianities are below a per cent, at most. While all-sky observations of the 21-cm background at arcminute-scale resolution can be used to measure the bispectrum with a cumulative signal-to-noise ratio of the order of a few tens, in the presence of foregrounds and instrumental noise related to first-generation interferometers, the measurement is unlikely to be feasible. For most purposes, non-Gaussianities in 21-cm fluctuations can be ignored and the distribution can be described with Gaussian statistics. Because 21-cm fluctuations are significantly contaminated by foregrounds, such as galactic synchrotron or low-frequency radio point sources, the lack of significant non-Gaussianity in the signal suggests that any significant detection of non-Gaussianity could be the result of foregrounds. Similarly, in addition to the frequency information that is now proposed to separate 21-cm fluctuations from foregrounds, if the non-Gaussian structure of foregrounds is known a priori, this additional information could potentially be used to reduce the confusion further.     

Modelling interaction of relativistic and non-relativistic winds in binary system PSR B1259−63/SS2883 – I. Hydrodynamical limit

In this paper, we present a detailed hydrodynamical study of the properties of the flow produced by the collision of a pulsar wind with the surrounding in a binary system. This work is the first attempt to simulate interaction of the ultrarelativistic flow (pulsar wind) with the non-relativistic stellar wind. Obtained results show that the wind collision could result in the formation of an 'unclosed' (at spatial scales comparable to the binary system size) pulsar wind termination shock even when the stellar wind ram pressure exceeds significantly the pulsar wind kinetic pressure. Moreover, the post-shock flow propagates in a rather narrow region, with very high bulk Lorentz factor (γ∼ 100). This flow acceleration is related to adiabatic losses which are purely hydrodynamical effects. Interestingly, in this particular case, no magnetic field is required for formation of the ultrarelativistic bulk outflow. The obtained results provide a new interpretation for the orbital variability of radio, X-ray and gamma-ray signals detected from binary pulsar system PSR B1259−63/SS2883.     

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

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

一种基于短时傅立叶变换的脉冲星消色散算法

叙述了一种基于短时傅立叶变换的脉冲星消色散算法。首先介绍了脉冲星的色散形成原因以及几种当今流行的消色散算法。接着介绍了短时傅立叶变换,并在此基础上提出了一种基于短时傅立叶变换的消色散算法。然后详细讲述了这种消色散算法的具体实现步骤,并且比较了在选取不同长度的时间窗函数的情况下,消色散的处理过程和残余色散量。最后,通过将这种算法的计算量以及结果和其他几种算法做横向比较,得出了结论:这种消色散算法的计算量小,实现简单,可以有效的进行消色散处理。     

The cross-correlation between the microwave and X-ray backgrounds: foregrounds and systematics

In this paper, we describe a novel experiment for the accurate estimation of pulsar dispersion measures (DMs) using the Giant Metrewave Radio Telescope. This experiment was carried out for a sample of 12 pulsars, over a period of more than one year (2001 January to 2002 May) with observations about once every fortnight. At each epoch, the pulsar DMs were obtained from simultaneous dual-frequency observations, without requiring any absolute timing information. The DM estimates were obtained from both the single-pulse data streams and from the average profiles. The accuracy of the DM estimates at each epoch is ∼1 part in 10⁴ or better, making the data set useful for many different kinds of studies.
The time-series of DMs shows significant variations on time-scales of weeks to months for most of the pulsars. An analysis of the mean DM values from these data shows significant deviations from catalogue values (as well as from other estimates in the literature) for some of the pulsars, with PSR B1642−03 showing the most notable differences. From our analysis results it appears that the constancy of pulsar DMs (at the level of 1 in 10³ or better) cannot be taken for granted. For PSR B2217+47, we see evidence of a large-scale DM gradient over a 1-yr period, which is modelled as being due to a blob of enhanced electron density sampled by the line of sight. For some pulsars, including pulsars with fairly simple profiles such as PSR B1642−03, we find evidence for small changes in DM values for different frequency pairs of measurement, a result that needs to be investigated in detail. Another interesting result is that we find significant differences in DM values obtained from average profiles and single-pulse data.     

Interstellar holography

The dynamic spectrum of a radio pulsar is an in-line digital hologram of the ionized interstellar medium. It has previously been demonstrated that such holograms permit image reconstruction, in the sense that one can determine an approximation to the complex electric field values as a function of Doppler shift and delay, but to date the quality of the reconstructions has been poor. Here we report a substantial improvement in the method which we have achieved by simultaneous optimization of the thousands of coefficients that describe the electric field. For our test spectrum of PSR B0834+06 we find that the model provides an accurate representation of the data over the full 63 dB dynamic range of the observations: residual differences between model and data are noise like. The advent of interstellar holography enables detailed quantitative investigation of the interstellar radio-wave propagation paths for a given pulsar at each epoch of observation. We illustrate this using our test data which show the scattering material to be structured and highly anisotropic. The temporal response of the medium exhibits a scattering tail which extends to beyond  100 μs  , and the centroid of the pulse at this frequency and this epoch of observation is delayed by approximately  15 μs  as a result of multipath propagation in the interstellar medium.