[Paper Review] Calibration of Quasi-Redundant Interferometers

Jonathan Sievers|arXiv (Cornell University)|Jan 7, 2017

Advanced Measurement and Metrology TechniquesEngineering16 citations

TL;DR

This paper introduces correlation calibration, a new method for calibrating quasi-redundant radio interferometers that relaxes the strict redundancy assumption by modeling baseline covariance and incorporating known source positions. It improves amplitude and phase reconstruction by a factor of ~2 and ~5 over traditional redundant calibration, respectively, especially under realistic array imperfections and noise levels.

ABSTRACT

High precison calibration is essential for a new generation of radio interferometers looking for Epoch of Reionization and Baryon Acoustic Oscillation signatures in neutral hydrogen. These arrays have so far been calibrated by redundant calibration, which usually assumes baselines intended to be identical are perfectly so. We present a new calibration scheme that relaxes the assumption of explicit redundancy by calculating the expected covariance of baselines. The technique also allows one to take advantage of partial knowledge of the sky, such as point sources with known positions but unknown fluxes. We describe a 2-level sparse matrix inverse to make the calibration tractable for 1,000-element class interferometers. We provide a reference implementation and use it to test the calibration of simulations of an array with imperfectly located antennas observing Euclidean-distributed point sources. Including position information for a handful of the brightest sources, we find the amplitude/phase reconstruction improves by a factor of $\sim$2/5 over redundant calibration for the noise levels/position errors adopted in the simulations. Inclusion of source positions also allows us to measure the overall phase gradient across the array, information which is lost in traditional redundant calibration.

Motivation & Objective

• To address the limitations of traditional redundant calibration, which assumes perfect baseline redundancy and ignores partial sky knowledge.
• To develop a calibration method robust to real-world array imperfections such as imperfect antenna placement and beam variations.
• To improve calibration accuracy by incorporating known positions of bright sources while maintaining sensitivity to diffuse foregrounds.
• To enable phase gradient estimation across the array, which is lost in standard redundant calibration.
• To provide a scalable, sparse matrix-based solution for large arrays (e.g., 1,000-element systems).

Proposed method

• The method formulates calibration as a likelihood problem using the expected covariance between baselines, derived from the sky power spectrum and baseline geometry.
• It introduces a 2-level sparse matrix inversion technique to make calibration computationally tractable for large arrays.
• The formalism incorporates known source positions by modeling their contribution to baseline covariance, improving calibration precision.
• It generalizes beyond strictly redundant arrays by treating baseline redundancy as approximate, allowing for small deviations due to positioning and beam errors.
• The approach uses a Gaussian random field assumption for the sky, enabling efficient likelihood-based inference without requiring a full sky model.
• A reference implementation is provided, supporting single-frequency calibration and extensible to polarization and bandpass calibration.

Experimental results

Research questions

• RQ1 How can calibration be improved in quasi-redundant interferometers when baseline redundancy is imperfect due to real-world array imperfections?
• RQ2 To what extent does including known source positions enhance calibration accuracy for diffuse signal detection?
• RQ3 Can a calibration method be developed that is robust to small baseline deviations and does not assume perfect redundancy?
• RQ4 How does the inclusion of phase information from known sources affect the reconstruction of overall array phase gradients?
• RQ5 What is the impact of sky model non-Gaussianity on the calibration process, particularly in the presence of phase correlations?

Key findings

• Including known positions for a handful of the brightest sources improves amplitude reconstruction by a factor of approximately 2 compared to traditional redundant calibration.
• Phase reconstruction accuracy improves by a factor of approximately 5 when source positions are included, reducing phase errors in the calibration solution.
• The method successfully recovers the overall phase gradient across the array, a quantity lost in standard redundant calibration.
• The calibration remains robust to non-Gaussian sky features, as phase correlations do not bias the solution, and amplitude non-Gaussianity has minimal impact under typical conditions.
• Errors in the assumed sky power spectrum have only mild effects on calibration, with no significant bias introduced under smooth, spatially varying power spectrum errors.
• The 2-level sparse matrix inversion enables tractable calibration for 1,000-element interferometers, making the method scalable to next-generation arrays.