No Arabic abstract
In this paper we present DMC, a model and associated tool for polarimetric imaging of very long baseline interferometry datasets that simultaneously reconstructs the full-Stokes emission structure along with the station-based gain and leakage calibration terms. DMC formulates the imaging problem in terms of posterior exploration, which is achieved using Hamiltonian Monte Carlo sampling. The resulting posterior distribution provides a natural quantification of uncertainty in both the image structure and in the data calibration. We run DMC on both synthetic and real datasets, the results of which demonstrate its ability to accurately recover both the image structure and calibration quantities as well as to assess their corresponding uncertainties. The framework underpinning DMC is flexible, and its specific implementation is under continued development.
The Murchison Widefield Array (MWA), located in Western Australia, is one of the low-frequency precursors of the international Square Kilometre Array (SKA) project. In addition to pursuing its own ambitious science program, it is also a testbed for wide range of future SKA activities ranging from hardware, software to data analysis. The key science programs for the MWA and SKA require very high dynamic ranges, which challenges calibration and imaging systems. Correct calibration of the instrument and accurate measurements of source flux densities and polarisations require precise characterisation of the telescopes primary beam. Recent results from the MWA GaLactic Extragalactic All-sky MWA (GLEAM) survey show that the previously implemented Average Embedded Element (AEE) model still leaves residual polarisations errors of up to 10-20 % in Stokes Q. We present a new simulation-based Full Embedded Element (FEE) model which is the most rigorous realisation yet of the MWAs primary beam model. It enables efficient calculation of the MWA beam response in arbitrary directions without necessity of spatial interpolation. In the new model, every dipole in the MWA tile (4 x 4 bow-tie dipoles) is simulated separately, taking into account all mutual coupling, ground screen and soil effects, and therefore accounts for the different properties of the individual dipoles within a tile. We have applied the FEE beam model to GLEAM observations at 200 - 231 MHz and used false Stokes parameter leakage as a metric to compare the models. We have determined that the FEE model reduced the magnitude and declination-dependent behaviour of false polarisation in Stokes Q and V while retaining low levels of false polarisation in Stokes U.
Extension of very long baseline interferometry (VLBI) to observing wavelengths shorter than 1.3mm provides exceptional angular resolution (~20 micro arcsec) and access to new spectral regimes for the study of astrophysical phenomena. To maintain phase coherence across a global VLBI array at these wavelengths requires that ultrastable frequency references be used for the heterodyne receivers at all participating telescopes. Hydrogen masers have traditionally been used as VLBI references, but atmospheric turbulence typically limits (sub) millimeter VLBI coherence times to ~1-30 s. Cryogenic Sapphire Oscillators (CSO) have better stability than Hydrogen masers on these time scale and are potential alternatives to masers as VLBI references. Here, We describe the design, implementation and tests of a system to produce a 10 MHz VLBI frequency standard from the microwave (11.2 GHz) output of a CSO. To improve long-term stability of the new reference, the CSO was locked to the timing signal from the Global Positioning System satellites and corrected for the oscillator aging. The long-term performance of the CSO was measured by comparison against a hydrogen maser in the same laboratory. The superb short-term performance, along with the improved long-term performance achieved by conditioning, makes the CSO a suitable reference for VLBI at wavelengths less than 1.3mm.
Space very long baseline interferometry (VLBI) has unique applications in high-resolution imaging of fine structure of astronomical objects and high-precision astrometry due to the key long space-Earth or space-space baselines beyond the Earths diameter. China has been actively involved in the development of space VLBI in recent years. This review briefly summarizes Chinas research progress in space VLBI and the future development plan.
Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.
New and upcoming radio interferometers will produce unprecedented amounts of data that demand extremely powerful computers for processing. This is a limiting factor due to the large computational power and energy costs involved. Such limitations restrict several key data processing steps in radio interferometry. One such step is calibration where systematic errors in the data are determined and corrected. Accurate calibration is an essential component in reaching many scientific goals in radio astronomy and the use of consensus optimization that exploits the continuity of systematic errors across frequency significantly improves calibration accuracy. In order to reach full consensus, data at all frequencies need to be calibrated simultaneously. In the SKA regime, this can become intractable if the available compute agents do not have the resources to process data from all frequency channels simultaneously. In this paper, we propose a multiplexing scheme that is based on the alternating direction method of multipliers (ADMM) with cyclic updates. With this scheme, it is possible to simultaneously calibrate the full dataset using far fewer compute agents than the number of frequencies at which data are available. We give simulation results to show the feasibility of the proposed multiplexing scheme in simultaneously calibrating a full dataset when a limited number of compute agents are available.