基于分布式执行框架的低频射电干涉阵列成像管线优化

平方公里阵列(Square Kilometre Array,SKA)项目是建设全球最大射电望远镜的国际合作项目,其灵敏度和测量速度将比当前所有的射电望远镜都要高出一个数量级.连续谱巡天是SKA的主要观测模式之一,基于连续谱成像建立巡天区域的标准星图,将能为后续天文科学研究奠定重要基础.银河系与河外星系全天默奇森宽场阵列拓展巡天(GaLactic and Extragalactic All-sky Murchison Widefield Array survey eXtended,GLEAM-X)是2018—2020年利用SKA先导望远镜默奇森宽场阵列(Murchison Wide-field Array,MWA)二期拓展阵列开展的新的射电连续谱巡天项目,观测期间积累了大量的低频巡天观测数据.海量观测数据的自动化、大批量处理是SKA望远镜项目所面临的的最大挑战和难题之一,基于分布式执行框架的成像管线优化经验将有助于解决海量数据处理问题.详细介绍了GLEAM-X成像管线并对其进行整合和改进,在中国SKA区域中心原型机(China SKA Regional Centre Prototype,...     

Research on Parallel Algorithms for uv-faceting Imaging

The uv-faceting imaging is one of the widely used large field of view imaging technologies, and will be adopted for the data processing of the low-frequency array in the first stage of the Square Kilometre Array (SKA1). Due to the scale of the raw data of SKA1 is unprecedentedly large, the efficiency of data processing directly using the original uv-faceting imaging will be very low. Therefore, a uv-faceting imaging algorithm based on the MPI (Message Passing Interface)+OpenMP (Open Multi-Processing) and a uv-faceting imaging algorithm based on the MPI+CUDA (Compute Unified Device Architecture) are proposed. The most time-consuming data reading and gridding in the algorithm are optimized in parallel. The verification results show that the results of the proposed two algorithms are basically consistent with that obtained by the current mainstream data processing software CASA (Common Astronomy Software Applications), which indicates that the proposed two algorithms are basically correct. Further analysis of the accuracy and total running time shows that the MPI+CUDA method is better than the MPI+OpenMP method in both the correctness rate and running speed. The performance test results show that the proposed algorithms are effective and have certain extensibility.     

Radio source calibration for the Very Small Array and other cosmic microwave background instruments at around 30 GHz

Accurate calibration of data is essential for the current generation of cosmic microwave background (CMB) experiments. Using data from the Very Small Array (VSA), we describe procedures which will lead to an accuracy of 1 per cent or better for experiments such as the VSA and CBI. Particular attention is paid to the stability of the receiver systems, the quality of the site and frequent observations of reference sources. At 30 GHz the careful correction for atmospheric emission and absorption is shown to be essential for achieving 1 per cent precision.
The sources for which a 1 per cent relative flux density calibration was achieved included Cas A, Cyg A, Tau A and NGC 7027 and the planets Venus, Jupiter and Saturn. A flux density, or brightness temperature in the case of the planets, was derived at 33 GHz relative to Jupiter which was adopted as the fundamental calibrator. A spectral index at ∼30 GHz is given for each.
Cas A, Tau A, NGC 7027 and Venus were examined for variability. Cas A was found to be decreasing at 0.394 ± 0.019 per cent yr⁻¹ over the period 2001 March to 2004 August. In the same period Tau A was decreasing at 0.22 ± 0.07 per cent yr⁻¹. A survey of the published data showed that the planetary nebula NGC 7027 decreased at 0.16 ± 0.04 per cent yr⁻¹ over the period 1967–2003. Venus showed an insignificant (1.5 ± 1.3 per cent) variation with Venusian illumination. The integrated polarization of Tau A at 33 GHz was found to be 7.8 ± 0.6 per cent at position angle =148°± 3°.     

Interferometer phase calibration sources — II. The region 0° ≤ δB1950≤ +20°

We present a catalogue of 781 compact radio sources in the declination range 0° ≤ δ_B1950≤ +20° whose positions have been measured to an rms accuracy of about 14 mas with the Very Large Array (VLA). These sources are primarily intended for use as phase calibration sources for the Jodrell Bank MERLIN. However, they will also be suitable as phase calibrators for the VLA and can be considered as candidate phase calibrators for very long-baseline interferometry (VLBI) networks.     

Wide-field VLBI imaging

We discuss the technique of Wide-field imaging as it applies to Very Long Baseline Interferometry (VLBI). In the past VLBI data sets were usually averaged so severely that the field-of-view was typically restricted to regions extending a few hundred milliarcseconds from the phase centre of the field. Recent advances in data analysis techniques, together with increasing data storage capabilities, and enhanced computer processing power, now permit VLBI images to be made whose angular size represents a significant fraction of an individual antenna's primary beam. This technique has recently been successfully applied to several large separation gravitational lens systems, compact Supernova Remnants in the starburst galaxy M82, and two faint radio sources located within the same VLA FIRST field. It seems likely that other VLBI observing programmes might benefit from this wide-field approach to VLBI data analysis. With the raw sensitivity of global VLBI set to improve by a factor 4–5 over the coming few years, the number of sources that can be detected in a given field will rise considerably. In addition, a continued progression in VLBI's ability to image relatively faint and extended low brightness temperature features (such as hot-spots in large-scale astrophysical jets) is also to be expected. As VLBI sensitivity approaches the μJy level, a wide-field approach to data analysis becomes inevitable.     

Numerical simulations of pinhole and single-mode fibre spatial filters for optical interferometers

We use a numerical simulation to investigate the effectiveness of pinhole spatial filters for optical/IR interferometers and to compare them with single-mode optical fibre spatial filters and interferometers without spatial filters. We show that fringe visibility measurements in interferometers containing spatial filters are much less affected by changing seeing conditions than equivalent measurements without spatial filters. This reduces visibility calibration uncertainties, and hence can reduce the need for frequent observations of separate astronomical sources for calibration of visibility measurements. We also show that spatial filters can increase the signal-to-noise ratios (SNRs) of visibility measurements and that pinhole filters give SNRs within 17 per cent of the values obtained with single-mode fibres for aperture diameters up to 3 r ₀. Given the simplicity of the use of pinhole filters we suggest that it represents a competitive, if not optimal, technique for spatial filtering in many current and next generation interferometers.     

Research on the Application of Fast-steering Mirror in Stellar Interferometers

In stellar interferometers, the fast-steering mirror (FSM) is widely utilized to correct the wavefront tilt caused by the atmospheric turbulence and internal instrumental vibration, because of its high resolution and fast response frequency. In this study, the non-coplanar error between the FSM and the actuator deflection axis introduced by the manufacturing, assembly, and adjustment is analyzed systematically. Via a numerical method, the additional optical path difference (OPD) caused by the above factors is studied, and its effect on the fringe tracking accuracy of a stellar interferometer is also discussed. On the other hand, the starlight parallelism between the beams of two arms is one of the main factors for the loss of fringe visibility. By analyzing the influence of wavefront tilt caused by the atmospheric turbulence on fringe visibility, a simple and efficient real-time correction scheme of starlight parallelism is proposed based on a single array detector. The feasibility of this scheme is demonstrated by a laboratory experiment. The results show that after the correction of fast-steering mirror, the starlight parallelism meets preliminarily the requirement of a stellar interferometer on the wavefront tilt.     

Overview and current status of the VSOP mission

An overview and status report of the VLBI Space Observatory Programme is presented.     

An Australia Telescope survey for CMB anisotropies

We have surveyed six distinct 'empty fields' using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration with the aim of imaging, with a high brightness sensitivity, any arcminute-scale brightness-temperature anisotropies in the background radio sky. The six well-separated regions were observed at a frequency of 8.7 GHz, and the survey regions were limited by the ATCA primary beams which have a full width at half-maximum of 6 arcmin at this frequency; all fields were observed with a resolution of 2 arcmin and an rms thermal noise of 24 μJy beam⁻¹. After subtracting foreground confusion detected in higher resolution images of the fields, residual fluctuations in Stokes I images are consistent with the expectations from thermal noise and weaker (unidentified) foreground sources; the Stokes Q and U images are consistent with expectations from thermal noise.
Within the sensitivity of our observations, we have no reason to believe that there are any Sunyaev–Zeldovich holes in the microwave sky surveyed. Assuming Gaussian-form CMB anisotropy with a 'flat' spectrum, we derive 95 per cent confidence upper limits of Q _flat<10–11 μK in polarized intensity and Q _flat<25 μK in total intensity. The ATCA filter function peaks at l =4700 and has half-maximum values at l =3350 and 6050.     

The Australia Telescope search for cosmic microwave background anisotropy

In an attempt to detect cosmic microwave background (CMB) anisotropy on arcmin scales, we have made an 8.7-GHz image of a sky region with a resolution of 2 arcmin and high surface brightness sensitivity using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration. The foreground discrete-source confusion was estimated from observations with higher resolution at the same frequency and in a scaled array at a lower frequency. Following the subtraction of the foreground confusion, the field shows no features in excess of the instrument noise. This limits the CMB anisotropy flat-band power to Q _flat < 23.6 μ K with 95 per cent confidence; the ATCA filter function (which is available at the website www.atnf.csiro.au/Research/cmbr/cmbr_atca.html) F _l in multipole l -space peaks at l _eff = 4700 and has half-maximum values at l  = 3350 and 6050.     

Maximum-likelihood astrometric geometry calibration of interferometric telescopes: application to the Very Small Array

Interferometers require accurate determination of the array configuration in order to produce reliable observations. A method is presented for finding the maximum-likelihood estimate of the telescope geometry, and of other instrumental parameters, astrometrically from the visibility timelines obtained from observations of celestial calibrator sources. The method copes systematically with complicated and unconventional antenna and array geometries, with electronic bandpasses that are different for each antenna radiometer, and with low signal-to-noise ratios for the calibrators. The technique automatically focuses on the geometry errors that are most significant for astronomical observation. We apply this method to observations made with the Very Small Array and constrain some 450 telescope parameters, such as the antenna positions, effective observing frequencies and correlator amplitudes and phaseshifts; this requires only ∼1 h of CPU time on a typical workstation.     

A Composite Model for GHz Peaked Spectra of Radio Sources

GHz-Peaked Spectrum (GPS) radio sources are powerful and compact sources with convex spectra. With increasing observational findings, it has been realized that either Synchrotron Self-Absorption (SSA) alone or Free-Free Absorption (FFA) alone is not enough to account for all the spectral features. We present a model consisting of an SSA region partially covered by FFA plasma, and derive a composite spectral formula. By applying the model to a sample of 19 GPS sources having strong absorption, it is found that the external FFA process makes the SSA peak frequency linearly shift to a higher (observed) peak frequency. The shift indicates that the FFA does play a role at the frequency close to the observed peak frequency.     

An active fiber sensor for mirror vibration metrology in astronomical interferometers

We present a fiber sensor based on an active integrated component which could be effectively used to measure the longitudinal vibration modes of telescope mirrors in an interferometric array. We demonstrate the possibility to measure vibrations with frequencies up to ⋍100 Hz with a precision better than 10 nm (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Smear fitting: a new image-deconvolution method for interferometric data

A new technique is presented for producing images from interferometric data. The method, 'smear fitting', makes the constraints necessary for interferometric imaging double as a model, with uncertainties, of the sky brightness distribution. It does this by modelling the sky with a set of functions and then convolving each component with its own elliptical Gaussian to account for the uncertainty in its shape and location that arises from noise. This yields much sharper resolution than clean for significantly detected features, without sacrificing any sensitivity. Using appropriate functional forms for the components provides both a scientifically interesting model and imaging constraints that tend to be better than those used by traditional deconvolution methods. This allows it to avoid the most serious problems that limit the imaging quality of those methods. Comparisons of smear fitting to clean and maximum entropy are given, using both real and simulated observations. It is also shown that the famous Rayleigh criterion (resolution = wavelength/baseline) is inappropriate for interferometers as it does not consider the reliability of the measurements.     

基于天文开源软件的卫星干涉测量数据处理系统

VLBI (Very Long Baseline Interferometry)技术观测卫星需要对干涉测量数据进行相关和后处理,通过相关、时延校准、条纹搜索,最终得到卫星的基线几何时延.基于天文开源软件建立起一套卫星干涉测量数据处理系统.该系统可工作在实时和事后两种状态,实现相关、中性大气、电离层、钟模型以及仪器硬件的时延校准、条纹搜索、生成基线时延和时延率序列.使用该系统处理北斗GEO (Geosynchronous Earth Orbit)卫星的干涉测量试验数据,得到了精度在1–2 ns量级的卫星基线时延序列.     

Global very long-baseline interferometry observations of compact radio sources in M82

Observations of the starburst galaxy, M82, have been made with a 20-station global very long-baseline interferometry (VLBI) array at λ 18 cm. Maps are presented of the brightest young supernova remnants (SNR) in M82 and the wide-field mapping techniques used in making images over a field of view of ∼1 arcmin with 3-milliarcsecond resolution are discussed. A limit has been placed on the power-law deceleration of the young supernova remnant (SNR) 43.31+592 with an index greater than 0.73±0.11 from observations with the European VLBI Network. Using the global array we have resolved compact knots of radio emission in the source which, with future global observations, will enable better constraints to be placed on the expansion parameters of this SNR.
The latest global observations have also provided high-resolution images of the most compact radio source in M82, 41.95+575. We determine an upper limit to the radial expansion rate along the major axis of 2000 km s⁻¹. However, the new images also show structure resembling that of collimated ejection which brings into question the previous explanation of the source as being a result of the confinement of a supernova by a high-density circumstellar medium.
It is apparent that we are now able to image the brightest supernova remnants in M82 with a linear scale which allows direct comparison with galactic SNR such as Cassiopeia A.     

Neural networks and the separation of cosmic microwave background and astrophysical signals in sky maps

We implement an independent component analysis (ICA) algorithm to separate signals of different origin in sky maps at several frequencies. Owing to its self-organizing capability, it works without prior assumptions on either the frequency dependence or the angular power spectrum of the various signals; rather, it learns directly from the input data how to identify the statistically independent components, on the assumption that all but, at most, one of the components have non-Gaussian distributions.
We have applied the ICA algorithm to simulated patches of the sky at the four frequencies (30, 44, 70 and 100 GHz) used by the Low Frequency Instrument of the European Space Agency's Planck satellite. Simulations include the cosmic microwave background (CMB), the synchrotron and thermal dust emissions, and extragalactic radio sources. The effects of the angular response functions of the detectors and of instrumental noise have been ignored in this first exploratory study. The ICA algorithm reconstructs the spatial distribution of each component with rms errors of about 1 per cent for the CMB, and 10 per cent for the much weaker Galactic components. Radio sources are almost completely recovered down to a flux limit corresponding to ≃0.7 σ _CMB, where σ _CMB is the rms level of the CMB fluctuations. The signal recovered has equal quality on all scales larger than the pixel size. In addition, we show that for the strongest components (CMB and radio sources) the frequency scaling is recovered with per cent precision. Thus, algorithms of the type presented here appear to be very promising tools for component separation. On the other hand, we have been dealing here with a highly idealized situation. Work to include instrumental noise, the effect of different resolving powers at different frequencies and a more complete and realistic characterization of astrophysical foregrounds is in progress.