PSR J0034-0721强单个脉冲的观测研究

利用南山基地25 m射电望远镜在1.54 GHz频段对脉冲星PSRJ0034-0721强单个脉冲进行了观测.使用单脉冲探测方法,从1 h观测数据中探测到116个信噪比R_(SN)N≥5的单个脉冲信号.在1.54 GHz频段探测到的单个脉冲R_(SN)从5到10.5,峰值流量约是平均脉冲峰值流量的14～29倍,远小于典型巨脉冲强度与其平均脉冲强度的比.这些脉冲的强度累积分布基本符合幂律谱,拟合得到谱指数α=-4.3±0.4.本次观测对R_(SN)≥5的单个脉冲探测率为3%,对R_(SN)≥10的单个脉冲探测率约为0.08%.这些脉冲的半峰线宽(W_(50))从1.6到8 ms,平均为3.9 ms.探测到的绝大多数强单脉冲发射相位集中分布在平均脉冲轮廓的峰值相位处,但也探测到两个R_(SN)8.5的强脉冲相位提前平均脉冲轮廓峰值相位33 ms左右,表明强脉冲发射区可能有两个,符合前人在40 MHz和111 MHz频段的观测结果,但1.54 GHz频段的平均脉冲轮廓只显示一个成份.     

PSR J0835-4510单脉冲的观测研究

利用中国科学院国家授时中心昊平观测站40m射电望远镜, 在L波段对Vela脉冲星(PSR J0835-4510)进行了单个脉冲观测研究. 在56min的观测数据中, 共观测到38040个单脉冲. 探测到观测时间内辐射的所有单脉冲信号, 其中单脉冲的半峰线宽(half-maximum line width, $W_{50     

基于泊松模型的X射线脉冲星信号的最大似然TOA估计

讨论了X射线脉冲星光子到达的周期平稳Poisson模型和到达时刻(Time of arrival,TOA)估计问题.在此基础上,得到了脉冲到达时刻的最大似然估计和克拉美-罗限(Cramer-Rao boundary,CRB),并推导了低信噪比情况下的似然函数近似表达和克拉美-罗限.利用解析的脉冲轮廓,对PSR B1821-24的TOA估计进行了Monte Carlo仿真,讨论了不同观测时间和信噪比下的估计误差,给出了相应的信噪比门限.研究表明:该分析方法能够有效估计X射线脉冲星的TOA定时精度,有利于评估其在各种应用中的性能.     

基于小波模极大值相关信息的脉冲星标准脉冲轮廓累积方法

讨论了产生脉冲星标准脉冲轮廓的重要性.在分析经典脉冲轮廓生成方法的基础上,提出了基于小波模极大值相关处理的脉冲轮廓累加方法.首先将观测数据按周期分组,分别进行小波分解,生成粗尺度上以脉峰为特征点的小波系数;然后,选择某组高信噪比信号的系数为参考数据,分别与其它各组变换系数进行相关,确定脉峰的相对时延.最后,对分组的原始观测数据按此相对时延量进行累加,生成脉冲轮廓.观测信号可先通过背景噪声估计和平滑处理而得到信噪比较高的归一化标准轮廓.实例分析以及与模板相关法的比较结果表明:该方法简单而有效,不需要设计近似模板,且对噪声类型不敏感,具有较高的工程应用价值.     

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

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

PSR B1237+25的强脉冲辐射特性的观测与研究

强脉冲是一种特殊的单脉冲辐射,表现为较强的射电爆发.利用新疆天文台南山25 m射电望远镜研究了PSR B1237+25的强脉冲辐射特性.发现探测到的793个强脉冲出现在积分轮廓的所有辐射成分中,峰值流量密度是平均脉冲的10.2至82.5倍.用对数正态分布对强脉冲的峰值流量密度比、信噪比和脉冲宽度的分布进行了拟合.在15...     

脉冲星核成份和锥成份的宽度关系

来自脉冲星辐射束的中心和外围的成份一般分别称之为核成份和锥成份.前人研究发现脉冲星的核成份与锥成份的宽度之间有很好的0.7倍比例关系,并将此写入教科书.收集整理了33颗核、锥成份非常明显的脉冲星的观测数据,对这些脉冲星的积分脉冲轮廓进行高斯拟合和成份分离,从而得到核、锥成份的宽度及误差.分析了三峰和五峰脉冲星的核锥成份,发现核锥脉冲星其核和锥成份的宽度在统计上没有明显差别,并不存在前人发现的0.7倍的比例关系.     

脉冲星脉冲宽度的测量方法

脉冲宽度的测量方法有几种.最常用的脉冲宽度有W_pp（两峰值点的间距）、W₅₀（最高峰值强度50%处的宽度）和W₁₀（最高峰值强度10%处的宽度）.另外,还有测量两个最外成分各自强度50%处的宽度W_c50和各自强度10%处的宽度W_c10.利用16颗双锥峰脉冲星和7颗核锥三峰脉冲星样本在1.4 GHz的脉冲轮廓数据,测量了5种脉冲宽度.利用有磁倾角α和撞击角β数据的脉冲星,导出了对应的辐射束半径ρ,并验证其与脉冲星周期P之间的关系ρ∝P^-0.5.通过比较5种脉冲宽度导出的ρ与P关系的好坏,发现最外成分之间的宽度比前面3种宽度更好,其中W_c50宽度得出的相关关系最好.由此认为W_c50是最能反映辐射束宽度的测量值.还讨论了脉冲轮廓的对称性,发现后随的锥辐射成分与先导的锥辐射成分相比,更靠近核心辐射成分,两个锥辐射成分的宽度在统计上基本相同.     

地面稳相频标传输系统设计和测试

设计了一款1.5 GHz的稳相频标传输装置,具有实时和事后补偿能力,并且兼容了同轴电缆和光纤两种传输介质.对装置的稳相原理和性能进行了分析和评估,然后对同轴电缆稳相系统的性能进行了测试.结果表明,基于目前的测试环境,与未稳相时相比,实时补偿模式可以改善相位变化112倍以上,积分时间长于7 s时,频率稳定度明显提升,长于60s的频率稳定度有1个数量级的改善;事后补偿模式可以改善40倍相位波动,积分时间长于2.5 s频率稳定度明显提升,长于40 s的频率稳定度有1个数量级的改善.对于更长时间积分的频率稳定度,实时和事后补偿模式均可以改善1.5个数量级以上,装置能有效抑制因传输介质的温度系数等因素导致的缓慢伸缩效应.     

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

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

On the subpulse modulation, polarization and sub-beam carousel configuration of pulsar B1857−26

New Giant Metre-Wave Radio Telescope (GMRT) observations of the five-component pulsar B1857−26 provide detailed insight into its pulse-sequence modulation phenomena for the first time. The outer conal components exhibit a 7.4-rotation period, longitude-stationary modulation. Several lines of evidence indicate a carousel circulation time     of about 147 stellar rotations, characteristic of a pattern with 20 beamlets. The pulsar nulls some 20 per cent of the time, usually for only a single pulse, and these nulls show no discernible order or periodicity. Finally, the pulsar's polarization-angle traverse raises interesting issues: if most of its emission comprises a single polarization mode, the full traverse exceeds 180°; or if both polarization modes are present, then the leading and the trailing halves of the profiles exhibit two different modes. In either case, the rotating-vector model fails to fit the polarization-angle traverse of the core component.     

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

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

On the measurement precision of solar p-mode eigenfrequencies

We make use of 3456 d of observations of the low-ℓ p-mode oscillations of the Sun in order to study the evolution over time of the measurement precision of the radial eigenfrequencies. These data were collected by the ground-based Birmingham Solar-Oscillations Network (BiSON) between 1991 January and 2000 June. When the power spectrum of the complete time series is fitted, the analysis yields frequency uncertainties that are close to those expected from the returned coherence times of the modes. The slightly elevated levels compared with the prediction appear to be consistent with a degradation of the signal-to-noise ratio in the spectrum that is the result of the influence of the window function of the observations (duty cycle 71 per cent). The fractional frequency precision reaches levels of a several parts in 10⁶ for many of the modes. The corresponding errors reported from observations made by the GOLF instrument on board the ESA/NASA SOHO satellite, when extrapolated to the length of the BiSON data set, are shown to be (on average) about ∼25 per cent smaller than their BiSON counterparts owing to the uninterrupted nature of the data from which they were derived.
An analysis of the BiSON data in contiguous segments of different lengths, T , demonstrates that the frequency uncertainties scale as T ^−1/2. This is to be expected in the regime where the coherence (life) times of the modes, τ _{n ℓ}, are smaller than the observing time T (the 'oversampled' regime). We show that mode detections are only now beginning to encroach on the 'undersampled' regime (where   T < τ _{n ℓ})  .     

Giant radio pulses from the millisecond pulsar PSR B1937+21 at 327 MHz

Seven giant radio pulses were recorded from the millisecond pulsar PSR B1937+21 during ≈8.1 min observation by the Ooty Radio Telescope (ORT) at 326.5 MHz. Although sparse, these observations support most of the giant pulse behaviour reported at higher radio frequencies (430 to 2380 MHz). Within the main component of the integrated profile, they are emitted only in a narrow (≲47 μs) window of pulse phase, close to its peak. This has important implications for doing super-high precision timing of PSR B1937+21 at low radio frequencies.     

Pulsar PSR B2303+30: a single system of drifting subpulses, moding and nulling

Analyses of multiple pulse sequences of the pulsar PSR B2303+30 reveal two distinct emission modes. One mode (B) follows a steady even–odd pattern and is more intense. The second mode (Q) is characteristically weak, but has intermittent drift bands with a periodicity of approximately 3 P ₁/cycle, and nulls much more frequently than the B mode. Both modes occur with roughly equal frequency, and their profiles have a similar single-humped form with a slight asymmetry. Our observations and analyses strongly suggest that the subpulse drift rates in both modes are linked in a series of cycles, which can be modelled as relaxing oscillations in the underlying circulation rate.     

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.     

The topology and polarization of subbeams associated with the 'drifting' subpulse emission of pulsar B0943+10 – V. A new look at the low-frequency burst-mode emission

This paper reports new observations of pulsar B0943+10 carried out at the Pushchino Radio Astronomy Observatory (PRAO) at the low radio frequencies of 42, 62 and 112 MHz. B0943+10 is well known for its exquisitely regular burst-mode (B-mode) drifting subpulses as well as its weaker and chaotic quiescent mode. Earlier Arecibo investigations at 327 MHz have identified remarkable, continuous changes in its B-mode subpulse drift rate and integrated-profile shape with durations of several hours. These PRAO observations reveal that the changes in profile shape during the B-mode lifetime are strongly frequency dependent – namely the measured changes in the component amplitude ratio are more dramatic at 327 and 112 MHz as compared with those at 62 and 42 MHz. The differences, however, are most marked during the first several tens of minutes after B-mode onset; after an hour or so the profile shape changes tend to be more similar at all four frequencies. We also have found that the linear polarization of the integrated profile increases continuously throughout the lifetime of the B mode, going from hardly 10 per cent just after onset to some 40–50 per cent after several hours. Pulsar B0943+10's B mode thus provides a unique new opportunity to investigate continuous systematic changes in the plasma flow within the polar flux tube. While refraction in the pulsar's magnetosphere may well play some role, we find that the various frequency-dependent effects, both between and within the two modes, can largely be understood geometrically. If the modes and B-mode decay reflect systematic variations in the carousel-'spark' radius and emission height then a specific set of profile and linear polarization changes would be expected.     

Interaction between nulls and emission in the pulsar B0834+06

We present a detailed study of the single pulses of the bright radio pulsar B0834+06, and offer evidence that the dominant periodic modulation in this pulsar's emission governs the occurrence of nulls. The nulls of B0834+06 constitute approximately 9 per cent of the total pulses and we demonstrate that they do not occur at random in the pulse sequence. On the contrary, they are found to occur preferentially close to the minimum of the pulsar's emission cycle, whose period jitters around a central value of P ₃≈ 2.17 rotation periods. It is likely that the intrinsic duration of the nulls averages about 0.2 times the pulsar rotation period. Surprisingly, the clearly distinct population of nulls and partial nulls of B0834+06 exhibit a two-peak profile slightly broader than that of the normal emission. This is in contrast to the profile of extremely weak normal pulses, which is narrower than the overall profile. A flow/counterflow model for the pulsar's two components can reproduce the essential observed features of the emission in its dominant mode, with nulls occurring at the point where the minima of the two systems are aligned. This suggests that the observed nulling rate is determined by the chance positioning of our sightline with respect to the system. If the flow is interpreted as part of a circulating carousel, a fit yields a best estimate of 14 'sparks'.     

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.