脉冲星VLBI观测软件相关处理进展

VLBI是进行高分辨率脉冲星观测研究的重要手段.脉冲星信号是非常微弱的脉冲序列,其VLBI观测面临多种挑战.在数据相关过程中采用有效技术提取脉冲星信号可提高观测成功率和精度.DiFX(Distributed FX-style Software Correlator)是目前国际上流行的开源软件相关处理机,它采用非相干消色散技术和“脉冲星数据分箱技术”(Pulsar Binning),在脉冲星VLBI观测数据的相关处理方面具有优异的性能.介绍了DiFX的构架,安装、调试方法,并对利用DiFX处理脉冲星VLBI观测数据的进展情况进行论述.利用单机环境下的DiFX,在普通模式和Pulsar Binning模式下对中国VLBI网(CVN)第一次脉冲星观测数据进行相关处理;利用德国马普射电天文研究所集群计算环境下的Bonn-DiFX,在Pulsar Binning模式下对流量仅有2.6 mJy的毫秒脉冲星PSRJ1022+1001的欧洲VLBI网(EVN)观测数据成功进行相关处理.最后,对使用DiFX处理脉冲星观测数据情况进行总结,并展望了今后CVN开展脉冲星观测研究的前景.     

Spatial modulation search applied to the search and confirmation of highly scintillated pulsars at FAST with a pulsar discovered in M3

We present a pulsar candidate identification and confirmation procedure based on a position-switch mode during the pulsar search observations. This method enables the simultaneous search and confirmation of a pulsar in a single observation, by utilizing the different spatial features of a pulsar signal and radio frequency interference(RFI). Based on this method, we performed test pulsar search observations in globular clusters M3, M15 and M92. We discovered and confirmed a new pulsar, M3 F, and detected the known pulsars M3 B, M15 A to G(except C) and M92 A.     

快速射电暴观测数据干扰缓解方法研究

从海量的天文观测数据中快速搜寻罕见的快速射电暴(Fast Radio Burst, FRB)事件, 干扰缓解是其中一项关键而具有挑战的工作. 射频干扰(Radio Frequency Interference, RFI)会淹没真实的天文事件, 还会导致搜寻管线输出大量的假阳性候选体. 由于干扰来源及其种类的复杂性, 目前并没有一种通用的方法可以解决这个问题. 为了降低干扰对FRB观测搜寻的影响, 分析和研究了南山26m射电望远镜L波段观测数据中的干扰情况, 针对主要的窄带干扰和宽带干扰建立了3层次的干扰缓解处理流程, 从而有效缓解了观测数据的干扰污染情况. 将该流程嵌入到FRB色散动态谱搜寻(Dispersed Dynamic Spectra Search, DDSS)管线中, 实验结果表明, 搜寻管线的检测率和检测精度得到了进一步的提高. 该方法为FRB观测数据干扰缓解处理提供了有价值的参考.     

Radio frequency interference mitigation using pseudoinverse learning autoencoders

Radio frequency interference(RFI) is an important challenge in radio astronomy. RFI comes from various sources and increasingly impacts astronomical observation as telescopes become more sensitive. In this study, we propose a fast and effective method for removing RFI in pulsar data. We use pseudo-inverse learning to train a single hidden layer auto-encoder(AE). We demonstrate that the AE can quickly learn the RFI signatures and then remove them from fast-sampled spectra, leaving real pulsar signals. This method has the advantage over traditional threshold-based filter method in that it does not completely remove contaminated channels, which could also contain useful astronomical information.     

基于卷积神经网络的快速射电暴候选体分类

针对目前从海量的快速射电暴(Fast Radio Burst, FRB)候选体中人工筛选FRB事件难以为继的现状,提出了一种基于卷积神经网络(Convolutional Neural Networks, CNN)的FRB候选体分类方法.首先,通过真实的观测数据和仿真FRB组成训练和测试样本集.其次,建立了二输入的深度卷积神经网络模型,并对其进行训练、测试和优化,获取FRB候选体分类器.最后,利用来自脉冲星的单脉冲数据对该分类器的有效性和性能进行了验证.实验结果表明,该方法可以快速从大量候选体中准确识别出单脉冲事件,极大地提高了FRB候选体的处理速率和效率.     

Repeating fast radio bursts reveal the secret of pulsar magnetospheric activity

The puzzling mechanism of coherent radio emission remains unknown, but fortunately, repeating fast radio bursts (FRBs) provide a precious opportunity, with extremely bright subpulses created in a clear and vacuum-like pulsar magnetosphere. FRBs are millisecond-duration signals that are highly dispersed at distant galaxies but with uncertain physical origin(s). Coherent curvature radiation by bunches has already been proposed for repeating FRBs. The charged particles are created during central star's quakes, which can form bunches streaming out along curved magnetic field lines, so as to trigger FRBs. The nature of narrow-band radiation with time-frequency drifting can be a natural consequence that bunches could be observed at different times with different curvatures. Additionally, high linear-polarization can be seen if the line of sight is confined to the beam angle, whereas the emission could be highly circular-polarized if off-beam. It is also discussed that pulsar surface may be full of small hills (i.e., zits) which would help producing bulk of energetic bunches for repeating FRBs as well as for rotation-powered pulsars.     

The radio environment of the 21 Centimeter Array:RFI detection and mitigation

Detection and mitigation of radio frequency interference(RFI) is the first and also the key step for data processing in radio observations, especially for ongoing low frequency radio experiments towards the detection of the cosmic dawn and epoch of reionization(Eo R). In this paper we demonstrate the technique and efficiency of RFI identification and mitigation for the 21 Centimeter Array(21CMA), a radio interferometer dedicated to the statistical measurement of Eo R. For terrestrial, man-made RFI, we concentrate mainly on a statistical approach by identifying and then excising non-Gaussian signatures, in the sense that the extremely weak cosmic signal is actually buried under thermal and therefore Gaussian noise. We also introduce the so-called visibility correlation coefficient instead of conventional visibility, which allows a further suppression of rapidly time-varying RFI. Finally, we briefly discuss removals of the sky RFI, the leakage of sidelobes from off-field strong radio sources with time-invariant power and a featureless spectrum. It turns out that state of the art technique should allow us to detect and mitigate RFI to a satisfactory level in present low frequency interferometer observations such as those acquired with the 21 CMA, and the accuracy and efficiency can be greatly improved with the employment of low-cost, high-speed computing facilities for data acquisition and processing.     

用于雷达天文研究的高速基带数据记录系统

为了开展雷达天文科学研究, 将射电望远镜接收的雷达回波信号进行采集和记录, 研究基于SNAP (Smart Network ADC (Analog to Digital Converter) Processor)硬件实验板和快速存储服务器设计并开发了雷达天文基带数据采集与记录系统. 该系统采用CASPER (Collaboration for Astronomy Signal Processing and Electronics Research)提供的图形化FPGA (Field Programmable Gate Array)开发工具流, 设计了双通道、256MHz带宽信号采集和VDIF (VLBI (Very Long Baseline Interferometry) Data Interchange Format)基带数据输出固件程序; 基于HASHPIPE (High Availability SHared PIPeline Engine)多线程管理引擎开发了双万兆以太网口实时基带数据存储程序, 存储带宽达到1GB/s; 最后编写了VDIF格式到雷达天文格式的转换程序. 经过脉冲星信号观测实验检测, 该系统准确、可靠.     

EVN observations of low-luminosity flat-spectrum active galactic nuclei

We present and discuss the results of very-long baseline interferometry (VLBI, EVN) observations of three low-luminosity     broad emission line active galactic nuclei (AGNs) carefully selected from a sample of flat-spectrum radio sources (CLASS). Based on the total and the extended radio power at 5 and at 1.4 GHz respectively, these objects should be technically classified as radio-quiet AGN and thus the origin of their radio emission is not clearly understood. The VLBI observations presented in this paper have revealed compact radio cores which imply a lower limit on the brightness temperature of about     . This result rules out a thermal origin for the radio emission and strongly suggests an emission mechanism similar to that observed in more powerful radio-loud AGNs. Since, by definition, the three objects show a flat (or inverted) radio spectrum between 1.4 and 8.4 GHz, the observed radio emission could be relativistically beamed. Multi-epoch VLBI observations can confirm this possibility in two years' time.     

A model of fast radio bursts: collisions between episodic magnetic blobs

Fast radio bursts(FRBs) are bright radio pulses from the sky with millisecond durations and Jansky-level flux densities. Their origins are still largely uncertain. Here we suggest a new model for FRBs. We argue that the collision of a white dwarf with a black hole can generate a transient accretion disk, from which powerful episodic magnetic blobs will be launched. The collision between two consecutive magnetic blobs can result in a catastrophic magnetic reconnection, which releases a large amount of free magnetic energy and forms a forward shock. The shock propagates through the cold magnetized plasma within the blob in the collision region, radiating through the synchrotron maser mechanism,which is responsible for a non-repeating FRB signal. Our calculations show that the theoretical energetics, radiation frequency, duration timescale and event rate can be very consistent with the observational characteristics of FRBs.     

乌鲁木齐25m射电望远镜的 单天线观测研究

自1993年乌鲁木齐天文站25m射电望远镜建成以来,除了不断完善VLBI观测系统外,还选择发展了具有科学价值的单天线天体物理观测课题。其中基于常温接收机的脉冲星脉冲到达时间观测系统已经于1999年5～6月间建成。该系统建立在25m射电天线的18cm波段上,消色散采用了2×128×2．5MHz多通道滤波器和数字化器,并由PC机完成数据采集。同时进行的脉冲星工作还有92cm及其它波段的脉冲轮廓监测,对0329+54的多波段观测得到了它的频谱。在25m天线的1．3cm波段上建立了基于声表面波频谱仪和频率综合器的分子谱线观测系统,对水脉泽的观测已经发现了十几个可能的水脉泽源,观测结果正在认证当中。     

Design of a multi-function high-speed digital baseband data acquisition system

A multi-function digital baseband data acquisition system is designed for the sampling,distribution and recording of wide-band multi-channel astronomical signals.The system hires a SNAP2board as a digital baseband converter to digitize,channelize and packetize the received signal.It can be configured dynamically from a single channel to eight channels with a maximum bandwidth of 4096 MHz.Eight parallel HASHPIPE instances run on four servers,each carrying two NVMe SSD cards,achieving a total continuous write rate of 8 GB s~(-1).Data are recorded in the standard VDIF file format.The system is deployed on a 25-meter radio telescope to verify its functionality based on pulsar observations.Our results indicate that during the 30-minute observation period,the system achieved zero data loss at a data recording rate of 1 GB s~(-1)on a single server.The system will serve as a verification platform for testing the functions of the QTT (Qi Tai radio Telescope) digital backend system.In addition,it can be used as a baseband/VLBI(Very Long Baseline Interferometry) recorder or D-F-engine of correlator/beamformer as well.     

乌鲁木齐天文站建立脉冲星相干消色散观测系统

脉冲星发射的辐射信号经过星际介质到达观测天线的过程中,存在色散效应。该效应导致有一定带宽脉冲信号的不同频率成份到达天线时间有延迟,影响对脉冲星的观测。消色散技术是脉冲星观测的关键技术,它对脉冲星观测系统的灵敏度和观测精度至关重要。脉冲星相干消色散过程是:通过对观测信号进行Nyquist采样,对采样数据做傅立叶变换,变换后的频域信号与星际介质Chirp函数乘积,然后再做逆傅立叶变换回时域,得到消色散后的时域信号。乌鲁木齐天文站依托现有的南山25m射电望远镜和VLBI记录终端MK5A系统,自行开发的相干消色散处理软件(Linux操作系统下C语言调用MPI库)和4节点机群系统,建立了脉冲星相干消色散观测系统。     

Peculiarities of processing VLBI observations of pulsars

We determined the features of pulsars that were disregarded in standard amplitude-calibration procedures for VLBI observations. We suggest additional amplitude-calibration methods. These methods were used to process the VLBI observations of the pulsar PSR B0329+54 carried out with the HALCA ground—space interferometer. Data from the space radio telescope are corrected for a nonuniform reception band. We estimated the diameter of the scattering disk for this pulsar at a frequency of 1600 MHz: \( < 1\mathop .\limits^{''} 8 \times 10^{ - 3}\).     

上海天文台25 m射电望远镜首次单天线脉冲星观测

介绍上海天文台25 m射电望远镜首次单天线脉冲星观测。2010年4月23日,使用上海天文台位于佘山观测基地的25 m射电望远镜对脉冲星J0332+5434在L波段进行了观测,此次观测使用VLBI终端进行数据采集记录,通过对观测数据进行非相干消色散和周期折叠,成功获得目标源的平均轮廓。此次观测的成功,表明该天线具备开展单天线脉冲星观测的条件,并为上海天文台建设中的65 m天线的天文观测提供了参考,为将来自主研发脉冲星终端进行了技术储备。     

基于小波变换的射频干扰消除方法的研究

在射电天文观测中,射频干扰(Radio Frequency Interference, RFI)会以多种形式混入望远镜接收系统,给观测带来误判或者降低观测信噪比.近年来国内国际射电天文快速发展,国内国际大型射电望远镜和阵列先后建设,观测灵敏度大为提高,射频干扰的影响尤为突出.随着科技发展和人类活动的加剧,射频干扰日益严重且不可逆转.提出利用2维离散小波变换的方法分析射电天文观测的数据,对望远镜系统输出的时间频率序列进行小波变换,根据小波系数分离出原始信号中各分量,每个分量统计得到相应的阈值,将各分量与阈值相比较识别干扰成分并标记去除.利用该方法对实际观测数据进行了处理,结果表明该方法能够很好地标记并消减干扰信号,且提高了观测的信噪比.     

The central point source in G76.9+1.0

We describe the serendipitous discovery of a very steep-spectrum radio point source in low-frequency Giant Metrewave Radio Telescope (GMRT) images of the supernova remnant (SNR) G76.9+1.0. The steep spectrum, as well as the location of the point source near the centre of this SNR confirm that this indeed is the pulsar J2022+3842. Archival Chandra X-ray data shows a point source coincident with the radio point source. However, no pulsed radio emission was detected despite deep searches at 610 MHz and 1160 MHz – which can be understood to be due to temporal broadening of the pulses. Weak pulsed emission has indeed been seen at 2 GHz with the Green Bank Telescope (GBT), establishing the fact that scattering is responsible for its non-detection at low radio frequencies. We underline the usefulness of low-frequency radio imaging as a good technique to prospect for pulsar candidates.     

大口径射电望远镜电磁兼容评估技术研究

电磁兼容性是设备或系统的重要性能指标, 也是保障系统的工作效能和提高系统可靠性的重要因素. 大口径射电望远镜运行阶段, 台址周围无线电业务及内部潜在的电磁干扰会降低观测系统灵敏度、影响天文观测的质量. 本论文针对拟建的新疆110 m全向可动射电望远镜(Qi Tai raido Telescope, QTT)开展了系统电磁兼容评估技术及控制方法研究, 具有重要的工程应用价值. 首先, 依据现有电波环境测量方法的不足, 深入分析了仪器设备的关键参数配置方法及测量时间计算方法, 采用Y因子法校准测量数据, 提出一种准实时电波环境测量方法. 面向高重复性宽带频谱, 分析了宽带频谱信号和噪声特征, 结合标准差理论, 提出一种基于邻值比较的信噪分离方法, 并采用邻值统计方法优化关键参数, 提高信噪分离精度. 针对QTT台址, 开发了自动化电波环境监测系统, 该系统6 GHz以下频段系统增益大于40 dB, 系统噪声系数小于2 dB, 测量不确定度小于1.49 dB, 具有极高的系统灵敏度和测量精度; 分析了频谱监测数据流, 设计了基于HDF5 (Hierarchical Data Format version 5)的数据存储格式, 开发了自动化电波环境测量和监控软件及数据处理软件. 依据QTT台址长期监测数据, 评估分析了台址电磁环境、主要干扰源特征及其影响. 其次, 提出大口径射电望远镜馈源口面干扰电平限值量化方法, 建立了基于台址地形的电波传播模型, 分析了现有电波传播模型的优缺点及适应性, 结合QTT台址实际地形及地质特征, 采用Longley-Rice和Two-Ray电波传播模型, 预测分析了QTT台址潜在干扰区域电磁干扰达到射电望远镜的电波路径衰减, 结合大口径射电望远镜天线增益量化方法, 提出设备所在位置干扰电平限值量化方法, 运用该方法对QTT台址潜在干扰区域的干扰电平限值进行量化. 依据设备所在位置干扰电平限值, 调研分析了国内外军用、民用电磁兼容测量标准, 结合电磁干扰对射电天文观测的影响, 提出一种大口径射电望远镜电磁兼容控制方法, 解决了现有电波暗室测量系统无法直接测量评估电子设备电磁兼容的问题, 该电磁兼容控制方法计划应用于QTT建设及运行阶段, 确保系统拥有良好的电磁兼容性. 最后, 依据QTT台址潜在干扰区域干扰电平限值, 结合典型电子设备电磁辐射频谱, 分析了QTT电磁兼容设计需求, 提出电磁兼容设计初步方案. 另外, 针对台址建筑设施内的中低电磁辐射干扰源, 提出一种低成本建筑屏蔽方法, 应用于QTT台址现有建筑.     

Study of Radio Frequency Interference Mitigation Method Based on Wavelet Transform

In radio astronomy observations, radio frequency interference (RFI) is mixed into the telescope receiving system in various forms. The existence of RFI brings misjudgment to the observation or reduces the observational signal-to-noise ratio. In recent years, the domestic and overseas radio astronomy has developed rapidly. Large-scale radio telescopes and telescope arrays at home and abroad have been constructed successively. Observation sensitivity is greatly improved, and the influence of RFI is particularly prominent. With the development of science and technology and the intensification of human activities, the RFI has become increasingly serious and irreversible. We propose to use the two-dimensional discrete wavelet transform method to analyze the data of radio astronomy observations, and to perform wavelet transform on the time-frequency sequence output by the telescope system. According to the wavelet coefficients, the every component in the original signal is separated, and the corresponding threshold value is obtained by the statistics on each component. Each component is compared with the threshold to identify the interference component, and to mark it for removal. This method is used to process the actual observation data. The results show that this method can well mark and eliminate the interference signal, and improve the observational signal-to-noise ratio.