4. Receiver Response

4.5. Calibration

The true visibility of a target source V_\text{ij}^\text{target}(u,v) observed by two telescopes i and j of an array is related to the observed visibility V_\text{ij}^\text{obs}(u,v) by

\displaystyle V_\text{ij}^\text{obs}(u,v)=G_\text{ij}(t)\cdot V_\text{ij}^\text{target}(u,v)\text{,}

in which G_\text{ij}(t) is the time-variable and complex gain factor of the correlation product of the two telescopes i and j. This gain factor includes variations caused by different effects, e.g., weather effects, effects of the atmospheric paths to the telescopes, and instrumental effects. Therefore, to determine the true visibility, the gain factor has to be calculated. This can be done via frequent observations of calibrator sources with known flux densities S_\text{cal} and accurately known positions. Furthermore, these sources have to be point sources, meaning that their angular size \Theta_\text{cal} is

\displaystyle \Theta_\text{cal} \ll \frac{\lambda}{D_\text{max}}\text{.}

Then, the complex gain G_\text{ij}(t) of each interferometer can be determined by

\displaystyle V_\text{ij}^\text{cal}(u,v)=G_\text{ij}(t)\cdot S_\text{cal}\text{,}

in which V_\text{ij}^\text{cal}(u,v) is the observed visibility of a calibrator source. This complex gain factor can then be used to calibrate the observed visibility of the target source. The true visibility of a target source V_\text{ij}^\text{target}(u,v) is then given by

\displaystyle V_\text{ij}^\text{target}(u,v)=\frac{V_\text{ij}^\text{obs}(u,v)}{G_\text{ij}(t)}=V_\text{ij}^\text{obs}(u,v)\cdot\frac{S_\text{cal}}{V_\text{ij}^\text{cal}(u,v)}\text{.}

Note that the gains G_\text{ij} can also be expressed by the voltage gain factors g_\text{i} and g_\text{j} of the individual i-th and j-th antennas by

\displaystyle G_\text{ij} = g_\text{i}\cdot g_\text{j}^\star \,\text{.}

This reduces the amount of calibration data, since there are many more correlated antenna pairs than antennas in large arrays. Furthermore, this makes the calibration procedure more flexible, since, for example, a source that is resolved at the longest baselines of an array can be used to compute the gain factors of the individual antennas from measurements made only at shorter baselines, at which the source appears unresolved.

Gibb's phenomenon

Fig. 2.34 Real and imaginary part of a frequency power spectrum for an idealized rectangular bandpass.

In continuum observations the correlation of the signals is usually measured for zero time lags, meaning that the correlation of the measured signal is averaged over the bandwidth, which is defined via a bandpass. In spectral line observations (see Chapt. 3 Sect. 4), measurements at different frequencies across the bandpass are required, which can be implemented by correlating the signals as a function of time lags. After Fourier transforming this correlation output, this yields the so-called cross power spectrum. This cross power spectrum is obtained for K frequency channels (for more detailed information see Chapt. 3 Sect. 4).

To calibrate an interferometer over the bandpass, which represents the filter characteristics of the last IF stage, a continuum point source is used to determine the complex gains of the interferometer, which yields the true cross-correlation spectrum of the observed target source. For calibrator sources with a constant continuum spectrum over the bandpass, the behavior of the complex gains g_\text{i} and g_\text{j} across the bandpass \Delta\nu are given by the correlated power as a function of frequency for each interferometer of the observing array. The real and imaginary parts of this correlated power, measured for a continuum source with a flat spectrum, are shown in Fig. 2.31 for an idealized rectangular bandpass. 

The observed visibility of an unknown target source is given by 

 \displaystyle V_\text{ij}^\text{obs}(\nu) = f(\nu)\star\left[g_\text{i}(\nu)\cdot g_\text{j}^\star(\nu)\cdot V(\nu)\right]\,\text{,}

in which f(\nu) is the Fourier transform of the truncating function in the frequency domain, which is a sinc function given by

 \displaystyle f(\nu) = \frac{\sin(\pi\cdot\nu\cdot K/\Delta\nu)}{\pi\cdot\nu}\,\text{.}

Here, V_\text{ij}^\text{obs} is given by a convolution with f(\nu), since the cross-correlation has only a finite maximum lag of 2K\cdot\Delta\nu. The observed visibility V_\text{ij}^\text{cal} of a calibrator source with |V| = 1 and \varphi_\text{V} = 0 is then given by

 \displaystyle V_\text{ij}^\text{cal}(\nu)= f(\nu)\star\left[g_\text{i}(\nu)\cdot g_\text{j}^\star(\nu)\right]\,\text{.} 

Finally, the true visibility V_\text{ij}^\text{true} of the observed target source is given by 

 \displaystyle  V_\text{ij}^\text{true}(\nu) = \frac{V_\text{ij}^\text{obs}(\nu)}{V_\text{ij}^\text{cal}(\nu)} = \frac{f(\nu)\star\left[g_\text{i}(\nu)\cdot g_\text{j}^\star(\nu)\cdot V(\nu)\right]}{f(\nu)\star\left[g_\text{i}(\nu)\cdot g_\text{j}^\star(\nu)\right]}\,\text{,}

in which the gains do not cancel out, due to the convolution. Therefore, the original frequency power spectrum, shown in Fig. 2.31, is convolved with the sinc function f(\nu), due to the finite lag, which will lead to large errors, if the target source has a spectral line near \nu=0. This is known as Gibbs' phenomenon.

Since the cross-correlation function has a complex power spectrum, Gibbs' phenomenon has strong influence on the imaginary part of the true visibility, because of the ratio in the equation for V_\text{ij}^\text{true}, which leads to non-negligible errors, especially for strong line emission near the step at \nu=0. Since the step in the imaginary part is large, these errors have a strong effect on the phases, producing strong ripples.