Calibration Challenges for Future Radio Telescopes

Instruments for radio astronomical observations have come a long way. While the first telescopes were based on very large dishes and 2-antenna interferometers, current instruments consist of dozens of steerable dishes, whereas future instruments will be even larger distributed sensor arrays with a hierarchy of phased array elements. For such arrays to provide meaningful output (images), accurate calibration is of critical importance. Calibration must solve for the unknown antenna gains and phases, as well as the unknown atmospheric and ionospheric disturbances. Future telescopes will have a large number of elements and a large field of view. In this case the parameters are strongly direction dependent, resulting in a large number of unknown parameters even if appropriately constrained physical or phenomenological descriptions are used. This makes calibration a daunting parameter estimation task, that is reviewed from a signal processing perspective in this article.Comment: 12 pages, 7 figures, 20 subfigures The title quoted in the meta-data is the title after release / final editing

Redundancy Calibration of Phased Array Stations

Our aim is to assess the benefits and limitations of using the redundant visibility information in regular phased array systems for improving the calibration. Regular arrays offer the possibility to use redundant visibility information to constrain the calibration of the array independent of a sky model and a beam models of the station elements. It requires a regular arrangement in the configuration of array elements and identical beam patterns. We revised a calibration method for phased array stations using the redundant visibility information in the system and applied it successfully to a LOFAR station. The performance and limitations of the method were demonstrated by comparing its use on real and simulated data. The main limitation is the mutual coupling between the station elements, which leads to non-identical beams and stronger baseline dependent noise. Comparing the variance of the estimated complex gains with the Cramer-Rao Bound (CRB) indicates that redundancy is a stable and optimum method for calibrating the complex gains of the system. Our study shows that the use of the redundant visibility does improve the quality of the calibration in phased array systems. In addition it provides a powerful tool for system diagnostics. Our results demonstrate that designing redundancy in both the station layout and the array configuration of future aperture arrays is strongly recommended. In particular in the case of the Square Kilometre Array with its dynamic range requirement which surpasses any existing array by an order of magnitude.Comment: 16 pages, 15 figures, accepted for publication in the A&A in Section 13, acceptance date: 1st May 2012. NOTE: Please contact the first author for high resolution figure

Non-linear Redundancy Calibration

For radio interferometric arrays with a sufficient number of redundant spacings the multiplicity of measurements of the same sky visibility can be used to determine both the antenna gains as well as the true visibilities. Many of the earlier approaches to this problem focused on linearized versions of the relation between the measured and the true visibilities. Here we propose to use a standard non-linear minimization algorithm to solve for both the antenna gains as well as the true visibilities. We show through simulations done in the context of the ongoing upgrade to the Ooty Radio Telescope that the non-linear minimization algorithm is fast compared to the earlier approaches. Further, unlike the most straightforward linearized approach, which works with the logarithms of the visibilities and the gains, the non-linear minimization algorithm leads to unbiased solutions. Finally we present error estimates for the estimated gains and visibilities. Monte-Carlo simulations establish that the estimator is indeed statistically efficient, achieving the Cramer-Rao bound.Comment: 9 pages, 5 figures. Accepted for publication in MNRAS. The definitive version will be available at http://mnras.oxfordjournals.or

Understanding synthesis imaging dynamic range

We develop a general framework for quantifying the many different contributions to the noise budget of an image made with an array of dishes or aperture array stations. Each noise contribution to the visibility data is associated with a relevant correlation timescale and frequency bandwidth so that the net impact on a complete observation can be assessed. All quantities are parameterised as function of observing frequency and the visibility baseline length. We apply the resulting noise budget analysis to a wide range of existing and planned telescope systems that will operate between about 100 MHz and 5 GHz to ascertain the magnitude of the calibration challenges that they must overcome to achieve thermal noise limited performance. We conclude that calibration challenges are increased in several respects by small dimensions of the dishes or aperture array stations. It will be more challenging to achieve thermal noise limited performance using 15 m class dishes rather than the 25 m dishes of current arrays. Some of the performance risks are mitigated by the deployment of phased array feeds and more with the choice of an (alt,az,pol) mount, although a larger dish diameter offers the best prospects for risk mitigation. Many improvements to imaging performance can be anticipated at the expense of greater complexity in calibration algorithms. However, a fundamental limitation is ultimately imposed by an insufficient number of data constraints relative to calibration variables. The upcoming aperture array systems will be operating in a regime that has never previously been addressed, where a wide range of effects are expected to exceed the thermal noise by two to three orders of magnitude. Achieving routine thermal noise limited imaging performance with these systems presents an extreme challenge. The magnitude of that challenge is inversely related to the aperture array station diameter.Comment: 27 pages, 24 figures, accepted in A&A, final versio

Precision Calibration of Radio Interferometers Using Redundant Baselines

Growing interest in 21 cm tomography has led to the design and construction of broadband radio interferometers with low noise, moderate angular resolution, high spectral resolution, and wide fields of view. With characteristics somewhat different from traditional radio instruments, these interferometers may require new calibration techniques in order to reach their design sensitivities. Self-calibration or redundant calibration techniques that allow an instrument to be calibrated off complicated sky emission structures are ideal. In particular, the large number of redundant baselines possessed by these new instruments makes redundant calibration an especially attractive option. In this paper, we explore the errors and biases in existing redundant calibration schemes through simulations, and show how statistical biases can be eliminated. We also develop a general calibration formalism that includes both redundant baseline methods and basic point source calibration methods as special cases, and show how slight deviations from perfect redundancy and coplanarity can be taken into account.Comment: 18 pages, 13 figures; Replaced to match accepted MNRAS versio

Phased Array Feed Calibration, Beamforming and Imaging

Phased array feeds (PAFs) for reflector antennas offer the potential for increased reflector field of view and faster survey speeds. To address some of the development challenges that remain for scientifically useful PAFs, including calibration and beamforming algorithms, sensitivity optimization, and demonstration of wide field of view imaging, we report experimental results from a 19 element room temperature L-band PAF mounted on the Green Bank 20-Meter Telescope. Formed beams achieved an aperture efficiency of 69% and system noise temperature of 66 K. Radio camera images of several sky regions are presented. We investigate the noise performance and sensitivity of the system as a function of elevation angle with statistically optimal beamforming and demonstrate cancelation of radio frequency interference sources with adaptive spatial filtering.Comment: 19 pages, 13 figure

Multisource Self-calibration for Sensor Arrays

Calibration of a sensor array is more involved if the antennas have direction dependent gains and multiple calibrator sources are simultaneously present. We study this case for a sensor array with arbitrary geometry but identical elements, i.e. elements with the same direction dependent gain pattern. A weighted alternating least squares (WALS) algorithm is derived that iteratively solves for the direction independent complex gains of the array elements, their noise powers and their gains in the direction of the calibrator sources. An extension of the problem is the case where the apparent calibrator source locations are unknown, e.g., due to refractive propagation paths. For this case, the WALS method is supplemented with weighted subspace fitting (WSF) direction finding techniques. Using Monte Carlo simulations we demonstrate that both methods are asymptotically statistically efficient and converge within two iterations even in cases of low SNR.Comment: 11 pages, 8 figure