In this paper we present the mathematical framework for
calibration of direction-dependent effects (DDE) not known
a priori. We also include a brief overview of the direction-
independent (DI) calibration and the mathematical formula-
tion of the existing DI SelfCal algorithm.
As an example of a DD calibration algorithm, we present
the Pointing Selfcal (PSC) algorithm. As in DI SelfCal, given
a model of the sky brightness distribution, the pointing off-
set vector per antenna is solved by iteratively minimizing
the residual vector with respect to the antenna-based point
offsets. We verified the PSC solver, first by applying it to
simulated wide-band data for the VLA at L-band and show
that the pointing offset vector is recovered correctly. We then
also apply the PSC algorithm to on-sky data from the VLA
at L-band. The VLA antenna optics has a polarization squint
which results in an angular separation between the right- and
left-circular polarization beams. To verify the PSC algo-
rithm, the solver was set-up to solve for the pointing offset
vector separately for the two polarizations. In the antenna
Az-El plane, the difference between the pointing offset vec-
tors for the two polarizations is a measure of the squint vec-
tor and the separation of the center of squint vector from the
origin gives a measure of the mechanical antenna pointing
offset. We verified that the PSC solver indeed recovers the
average squint vector, though antenna-to-antenna variations
were also large and significant. Some of the antennas had
large mechanical pointing offsets, which were subsequently
verified via independent measurement of the expected im-
provement in the antenna gain after correcting for them in
the telescope software.
We also discuss the noise budget for the PSC algorithm.
Analysis of the signal-to-noise ratio (SNR) available for the
PSC solver as a function of the wide-band sky brightness dis-
tribution and telescope parameters leads to the following con-
clusions:
1. The PSC algorithm is optimal in utilizing the SNR due
to the sky brightness distribution in the entire antenna
field of view, rather than, for example, a few bright
sources.
2. Antenna pointing offsets can be solved-for at high sig-
nificance with the instantaneous sensitivity of most
modern radio interferometric telescopes with wide-
band receivers and typical sky brightness distributions.
3. While the sky brightness at the center of the PB does
not contribute signal for the PSC solver, the contribu-
tion from around the center increases with the mag-
nitude of the pointing offsets. This signal also peaks
around the half-power points of the PB. The PSC al-
gorithm therefore works well for observations where
the sky brightness is distributed across the FoV. The
degradation in the imaging performance due to the an-
tenna pointing errors is also more significant for such
observations. Emission spread across the FoV is typi-
cal at frequencies below a few GHz and at much higher
frequencies where mosaic imaging of emission much
larger than the antenna PB is often necessary.
4. We expect PSC to scale well in a parallel comput-
ing environment. Simple parallelization by data-
partitioning is possible and efficient. Reduction in
the run time by large factors using multi-threaded re-
samplers (Golap 2015) deployed on massively parallel
hardware remains a possibility, though more work is
needed in this area to arrive at an optimal computing
architecture.
5. The PSC is typically set-up for relatively short solu-
tion interval in time. For low-frequency observations,
the increase in the run time due to the w-term may be
mitigated by treating the array as co-planar for each
The mathematical framework for DD calibration presented
here can be extended for a Shape SelfCal (SSC) algorithm
to account for the change in the shape of the aperture us-
ing the low-order A-Solver approach. In the A-Solver ap-
proach, the parameters describing the physical structure of
the antenna are determined using a geometric optics predictor
for the antenna aperture illumination pattern (AIP) and holo-
Based on our estimate of the SNR typically available, it may
be possible to develop an SSC algorithm using a model of
the sky brightness distribution. More work is required, and
in progress, in this area.