A&A 375, 344-350 (2001)
DOI: 10.1051/0004-6361:20010799
S. Bhatnagar - R. Nityananda
National Centre for Radio Astrophysics (TIFR), Post Bag 3, Ganeshkhind, Pune 411 007, India
Received 10 April 2001 / Accepted 30 May 2001
Abstract
Mechanical and electronic imperfections can result into
polarization leakage in individual antennas of a radio interferometer.
Such leakages manifest themselves as closure errors even in co-polar visibility measurements of unpolarized sources. This paper describes
and tests a method for the study of polarization leakage for
radio interferometric telescopes using only the nominally
co-polar visibilities for unpolarized calibrators.
Interpretation of the resulting closure phases on the Poincaré
sphere is presented. Since unpolarized sources are used, the actual
solutions for leakage parameters is subject to a degeneracy which is
discussed. This however, does not affect the correction of closure
errors in our scheme.
Key words: methods: data analysis - techniques: polarimetric - techniques: interferometric
The mutual coherence function (also called the visibility function) for an unresolved and unpolarized source, measured by an interferometer array can be modeled as a product of antenna based complex gains. These complex gains can be derived from the measured visibility function using the standard algorithm, which we call antsol. antsol forms the central engine of most amplitude and phase calibration schemes used for radio interferometric data. (The earliest published reference for an algorithm for antsol of which we are aware is Thompson & D'Addario 1982.)
Usually one measures the components of the incident radiation along two orthogonal polarization states by using two separate feeds. The signals from the two feeds travel through essentially independent paths till the correlator. However, due to mechanical imperfections in the feed or imperfections in the electronics, the two signals can leak into each other at various points in the signal chain.
At the correlator, signals from all the antennas are multiplied with each other and the results averaged to produce the visibilities. The signals of same polarization are multiplied to produce the co-polar visibilities while the signals of orthogonal polarizations are multiplied to produce the cross-polar visibilities. The co-polar and cross-polar visibilities can be used to compute the full Stokes visibility function. Antenna based instrumental polarization and polarization leakage can be derived from the full Stokes coherence function for a source of known structure (usually an unresolved source) (Hamaker et al. 1996; Sault et al. 1996, henceforth HBS).
The correlator used for the Giant Metrewave Radio Telescope (GMRT) by default computes the co-polar visibilities using the Indian mode of the VLBA Multiplier and Accumulator (MAC) chip. Here we describe a method, which we call leaky antsol, for the computation of the leakages using only the co-polar visibility function for an unpolarized source. Following the notation used by HBS, we label the two orthogonal polarizations by and to remind us that the formulation is independent of the precise orthogonal pair of polarization states chosen.
Section 2 describes the motivation which led to this analysis. For orientation, Sect. 3 starts with the problem of solving for the usual complex antenna based gains and sets up an iterative scheme for the solution. The problem of simultaneously solving for the complex antenna gains and leakages is then posed in Sect. 3.1 and a similar iterative scheme is set up. Section 3.2 presents the results of the simulations done to test the scheme. Section 4.1 presents some results using the GMRT at 150 MHz. Also, we were fortunate to have the L-band feeds of one of the GMRT antennas converted from linear to circular polarization. We observed 3C147 in this mode where all baselines with this special antenna measured the correlation between nominally linear and circular polarization. Results of this experiment demonstrate that the leakage solutions are indeed giving information about the polarization properties of the feeds. These results and their interpretation on the Poincaré sphere are presented in Sect. 4.2. Section 5 gives the interpretation of the leakage solutions and discusses closure errors due to polarization leakage using the Poincaré sphere.
Rogers (1983) pointed out in the context of the VLBA, that the polarization leakage cause closure errors even in nominally co-polar visibilities. Massi et al. (1997) have carried out a detailed study of this effect for the telescopes of the European VLBI Network (EVN). Our motivation in this paper is that the current single sideband GMRT correlator uses the so called Indian mode of the VLBA MAC chips to produce only the co-polar visibilities. Also, the planned Walsh switching has not yet been implemented at the GMRT and in any case, would not eliminate leakage generated before the switching point. Tests done using strong point source dominated fields show unaccounted closure errors at a few percent level. The motivation behind developing an algorithm to solve for gains and leakages simultaneously, using only the co-polar visibilities was to determine if the measured closure errors could be due to polarization leakage in the system. Estimates of leakage can then be used in the primary calibration to remove the effects of polarization leakage. This is where this paper differs from the earlier work of HBS which is about the calibration using the full Stokes visibility function, needed for observations of polarized sources. The polarization leakage in some of the EVN antennas corrupts the co-polar visibilities at a level visible as a reduction in the dynamic range of the maps (Massi & Aaron 1997a,b; Massi et al. 1998). Thus such a method can also be used in imaging data from the EVN and other telescopes affected by such closure errors.
Let
represent the complex gain for the -polarization
channel of the ith antenna and
represent the leakage of
the q-polarization signal into the p-polarization channel. The
electric field measured by antenna i can then be written as
For an unpolarized point source
where I is the total flux
density. Writing
we get
Assuming s to be negligible, the usual antsol algorithm estimates s such that is minimized (see Sect. 3). Normally, Walsh switching (Thompson et al. 1986) is used to eliminate the polarization leakage due to cross-talk between the signal paths, such that . However, s can also be finite due to mechanical imperfections in the feed or the cross-polar primary beam, which cannot be eliminated by Walsh switching.
In the case of significant antenna based polarization leakage (compared to ), the second term in Eq. (3) involving s will combine with the closure noise . The polarization leakage therefore manifests itself as increased closure errors (see Sect. 5 for a geometric explanation on the Poincaré sphere).
In the absence of any polarization leakage, g_{i}s can be estimated by
minimizing
In Eq. (2), if accurately represents the source structure, X_{ij}^{pp} will have no source structure dependent terms and is purely a product of two antenna dependent complex gains. For a resolved source, can be estimated from the image of the source.
Evaluating
and equating it to
zero^{} (Bhatnagar 1999), we get
Since the antenna dependent complex gains also appear on the
right-hand side of Eq. (5), it has to be solved iteratively
starting with some initial guess for g_{j}s or initializing them all
to 1. Equation (5) can be written in the iterative form as:
In the presence of significant polarization leakage,
Eq. (3) can be used to re-write Eq. (4) as
(7) |
Equating the partial derivatives
,
to
zero, we get
Equation (3), which expresses the observed visibilities on a point source unpolarized calibrator in terms of the gains and leakage coefficients of the antennas, would take the same form if written in an arbitrary orthogonal basis. It is clear that the g's and the 's will change when we change the basis, so this means that the equations cannot have a unique solution. This situation is familiar from ordinary self-calibration, when only relative phases of antennas are determinate, with one antenna acting as an arbitrary reference. For observations of unpolarized sources, we can similarly say that any feed can be chosen as a reference polarization, with zero leakage, and other feeds have gains and leakages in the basis defined by this reference. Other conventions may be more convenient, as discussed in the appendix which discusses degeneracy in detail.
We simulated visibilities with varying fraction of polarization leakage in the antennas to test the algorithm as follows. The antenna based signal and leakage were constructed as and where and were drawn from the same Gaussian random population. The visibility from two antennas i and j was then constructed as for . This is equivalent to a visibility of an unpolarized point source of unit strength with a complex antenna based gain g_{i} and leakage of strength proportional to f. Equation (6) was then used to compute g_{i} and residual computed as . The computed values of were then used to compute improved estimates for by simultaneously solving for and using the iterative forms of Eqs. (8) and (9). The derived values of and matched the true values to within the tolerance limit. A new was computed as . The values of and as a function of f are plotted in Fig. 1. The two curves become distinguishable when the leakage is significantly greater than (for f greater than 1%). After that, the value of is consistently lower than , where the contribution of antenna based leakage has not been removed. Also notice that remains constant while quadratically increases as a function of f. This is due to the fact that antsol treats the antenna based polarization leakage as closure errors resulting in an increased with increasing fractional leakage.
Figure 1: Figure showing the results of the simulations. The top curve is the value of using the classical antsol ( ). The bottom curve is the value of using the leaky antsol ( ) as a function of the percentage polarization leakage. | |
Open with DEXTER |
Engineering measurements for polarization isolation at 150 MHz for the GMRT show significant polarization leakage in the system. We therefore used leaky antsol to calibrate the data from the Galactic plane phase calibrator 1830-36 which is known to be less than polarized at 1.4 GHz. The percentage polarization at 150 MHz is not known, but it is expected to decrease further and it was taken to be an unpolarized point source.
Fractional polarization leakage ( ) of up to 100% was measured for most of the antennas, which is consistent with the estimated leakage measured from system engineering tests. Again, and were computed and the results are shown in Fig. 2. The 150-MHz GMRT band suffers from severe radio frequency interference (RFI). The sharp rise in the value of around sample number 10 is due to one such RFI spike. This spike is present in the total power data from all antennas at this time. On an average, the reduces by 60% when leakage calibration is applied ( ). This is consistent with polarization leakage being a major source of non-closure at this frequency.
Figure 2: Figure showing the results using GMRT data at 150 MHz for the compact Galactic plane source 1830-36. The top curve is the value of using the classical antsol ( ). The bottom curve is the value of using the leaky antsol ( ) as a function of time. | |
Open with DEXTER |
The GMRT L-band feeds are linearly polarized. For the purpose of a VLBI experiment conducted in December 2000, the L-band feed of one of the antennas was converted to a circularly polarized feed. The rest of the L-band feeds were linearly polarized and we took this opportunity to measure correlations between the circularly polarized antenna with other linearly polarized antennas using the source 3C147. Two scans of approximately one hour long observations were done using the single side band GMRT correlator. This correlator computes only co-polar visibilities. With this configuration of feeds, visibilities between the circularly polarized antenna and all other linearly polarized antennas corresponds to correlation between the nominal X- and R-polarizations, labeled by RX, were recorded in the first scan. The polarization of the circularly polarized antenna was then flipped for the second scan to record the correlation between the nominal X- and L-polarization states, labeled by LX.
The VLA Calibrator Manual^{} lists the percentage polarization ( ) for 3C147 at L-band <0.1%. The cross-polar terms in Eq. (2), which are assumed to be zero, will therefore contribute an error of the order of . These cross-polar terms will be, however, multiplied by gains of type . Since and are both assumed to be uncorrelated between antennas, this error will manifest as random noise in Eq. (3). Within the limits of other sources of errors, the source 3C147 can therefore be considered to be a completely unpolarized source.
Figure 3: Figure showing the results using visibilities with one circularly polarized antenna and all other linearly polarized antennas at L-band. The x- and y-axis denote the real and imaginary parts of respectively. and were solved for every integration time (17 s). All linearly polarized antennas are close to the origin in this plot. The solutions for the circularly polarized antenna (C03) are the set of points away from the origin (shown by open circles and triangles). The two sets of points for this antenna, separated from each other by are solutions for the right- and left-circular polarization channel. The points denoted by open circles are from correlation between the right- and left-circular polarization of C03 with nominal linear X-polarization of the other antennas (labeled as C03^{RX}and C03^{LX} respectively). The points denoted by triangles are from correlation of C03 with nominal Y-polarization of the other antennas (labeled as C03^{RY} and C03^{LY}). | |
Open with DEXTER |
The response of an ideal circularly polarized antenna to unpolarized
incident radiation can be expressed as a superposition of two linear
polarization states as
where, the superscripts R, Xand Y denote the right circular and the two linear polarization
states respectively.
is half the phase difference between the two linear
polarization states and is equal to
for right-circular
polarization and
for left-circular polarization. Writing the
general Eq. (1) for right-circularly polarized antenna
as
and
substituting for
and
we get
(10) |
X^{RX}_{ij} | = | (11) | |
Let ( for the circularly polarized antenna). Then, the amplitude of P_{i}is a measure of the fractional polarization leakage in the antenna while the phase of P_{i} gives the phase difference between the signal from one of the feeds and the leaked signal from the other feed. For an ideal circularly polarized antenna, . A plot of the real and imaginary parts of this quantity for all antennas should therefore clearly show P_{i} for the circularly polarized antenna with an amplitude of 1 and at an angle of with respect to the nominal X-axis.
The real and imaginary parts of P_{i} for all antennas from this experiment are shown in Fig. 3. The solutions were computed for every integration cycle of 17 s and the points on this plot represent the tip of phasor P_{i}. The collection of points near the origin are for all the linearly polarized antennas while the collection of two sets of points away from the origin, approximately at an angle of from each other, are for the circularly polarized antenna. The solutions found by leaky antsol match the expected results quite well. This therefore constitutes a reasonably controlled test with real data showing that the solutions indeed provide information about the polarization leakage in the system.
This experiment however provides much more information about the polarization properties of the antenna feeds used. The collection of points in the first quadrant denoted by open circles are the values of P_{i} derived from the correlation between the nominal right-circularly polarized signal and the linearly polarized signals along the nominal X-axis from all other antennas. Points in the third quadrant are similarly derived using the left-circular signals. The set of points denoted by triangles in the second and fourth quadrant are derived using correlations of right- and left-circularly polarized signals with the linearly polarized signals along the nominal Y-axis from all other antennas.
A larger spread in the solutions using the left-circularly polarized signals indicates that the closure noise (from other unknown sources) in these signals is higher. The fact that the amplitude of P_{i}derived using the right-circularly polarized signals is 0.5 indicates that the nominal circularly polarized feed is in fact elliptically polarized with this axial ratio. The spread of about the origin is indicative of polarization leakage at the level of few percent in the linearly polarized antennas as well. The leakage in one of the linearly polarized antennas is significantly larger (4%). Since this kind of data is routinely taken on primary calibrators during GMRT observations for synthesis imaging, leaky antsol provides a useful diagnostic of system health, polarization performance and numbers needed to correct the data in high accuracy work.
The following test was also carried out to check that the closure phase on a triangle involving the circular feed was indeed mainly due to polarization effects. The three baselines making up this triangle were flagged as bad baselines from the input data and a new solution found for the gains and leakages of all antennas. This solution predicted the same closure phase (to within errors) as actually observed.
In this section we use right- and left-circular polarization states as the basis. A general elliptically polarized state can be written as a superposition of two states represented by the vector . Clearly, corresponds to linear polarization and to elliptical polarization. Increasing by rotates the direction of the linear state or the major axis of the ellipse by . We can chose the phase of the basis so that corresponding to linear polarization along the x-axis. The Poincaré sphere representation of the state of polarization maps the general elliptic state to the point ( ) on the sphere. The properties of this representation are reviewed by Ramachandran & Ramaseshan (1961). We are concerned here with one remarkable property, discovered by Pancharatnam (1956, 1975). Whenever there is constructive interference between two sources of radiation, it is natural to regard them as in phase. An unexpected property of this simple definition manifests itself when we consider 3 sources of radiation of different polarization - that if a source A is in phase with B and B in phase with C, C in general need not be in phase with A. The phase difference between A and C is known in the optics literature as the geometric or Pancharatanam phase (see also Ramaseshan & Nityananda 1986; Berry 1987). We show that this naturally occurs in radio interferometry of an unpolarized source with three antennas of different polarizations.
Let the polarization states of the three antennas be represented by , , and in a circular basis. Denoting the vector by , one clearly see that the visibility on the 1-2 baseline is proportional to . Hence the closure phase around a triangle made by antennas 1, 2, and 3 is the phase of the complex number (also called the triple product) . In the quantum mechanical literature, this type of quantity goes by the name of Bargmann's invariant and its connection to the geometric phase was made clear by Samuel & Bhandari (1988). With some work, one can give a general proof that the closure phase (phase of V_{123}) is equal to half the solid angle subtended at the centre of the Poincaré sphere by the points represented by , , and on the surface of the sphere. For the case where the polarization state of the three antennas are same, this phase is zero in general. However, when the polarization states of the antennas are different, this phase is non-zero.
The well known result that an arbitrary polarization state can be represented as a superposition of two orthogonal polarization states translates to representing any point on the Poincaré sphere by the superposition of two diametrically opposite states on a great circle passing through that point. For example, circular polarization can be expressed by two linear polarizations, each with intensity . In the context of the present work, the nominally circularly polarized antenna maps to a point away from the equator on the Poincaré sphere (it would be exactly on the pole if it is purely circular) while the rest of the antennas map close to the equator (they would be exactly on the equator if they are purely linear and map to a single point if they were also identical). The visibility phase due to the extra baseline based term in Eq. (3) due to polarization mis-match is a consequence of the Pancharatanam phase mentioned above. This phase, on a triangle involving the circularly polarized antenna, will be close to the angle between the two linear antennas. For example, if , , and , the phase of V_{123} will be . This picture can be depicted by plotting the real and imaginary parts of , which is done in Fig. 3. The circularly polarized antenna can be clearly located in this figure as the set of point away from the origin while the linearly polarized antennas as the set of points close to the origin. The collection of points located away but almost symmetrically about the origin represents the nominal right- and left-circularly polarized feeds. Points on the equator, but significantly away from the origin represents an imperfect linearly polarized antenna. Note that the average closure phase between the nominally linear antennas is close to zero, which defines the mean reference frame in Fig. 3.
Rogers (1983) pointed out that non ideal feed polarizations of the individual antennas of a radio interferometer can result into closure errors in the co-polar visibilities. In this paper we describe and demonstrate a method to measure the polarization leakage of individual antennas using the nominally co-polar visibilities for an unpolarized calibrator. This method can therefore be used as a useful tool for studying the polarization purity of the antennas of radio interferometers from the observations of unpolarized calibrators. However, since only unpolarized calibrators are used, the actual solution for the leakage parameters is subject to a degeneracy. This degeneracy does not affect the correction of the visibilities and can be used to remove the closure errors due to polarization leakage. Massi et al. (1997) have shown that such polarization leakage induced closure errors in the data from the EVN is the dominant effect of instrumental polarization. For the EVN, this effect can be seen as a reduction in the dynamic range of the images. Our method can be used for such data to remove these closure errors for unpolarized sources.
The general elliptic state of the polarization of radiation can be represented by a point on the Poincaré sphere. The phase difference between three coherent sources of radiation but with different states of polarization goes by the name of Pancharatanam or geometric phase in the optics literature. We interpret the co-polar visibilities with polarization leakages on the Poincaré sphere and show that the polarization induced closure phase errors in radio interferometers is same as the Pancharatanam phase of optics. The antenna based leakages also map to points on the Poincaré sphere and the ambiguity in the solution can be understood as a rigid rotation of the Poincaré sphere, which leaves the leakage solutions unchanged relative to each other.
Acknowledgements
We thank A. Omar and S. Roy who suggested that there could be closure errors in the data due to polarization leakage. Many useful discussions on this with A. P. Rao and Divya Oberoi are thankfully acknowledged. We also thank the telescope operators and other GMRT staff for their help and cooperation in making possible the observations used for this work. The GMRT is run by the National Centre for Radio Astrophysics of the Tata Institute of Fundamental Research.
This research has made use of NASA's Astrophysics Data System Abstract Service. All of this work was done using computers running the GNU/Linux operating system and it is a pleasure to thank all the numerous contributors to this software.
We discuss the non-uniqueness of the solutions of
Eq. (3), and possible convenient conventions for choosing
a specific solution. One obvious degeneracy is that multiplication of
all the 's by one common phase factor independent of antenna,
and all the g's, by another, in general different, common factor,
does not affect the right hand side of Eq. (3). Also, the
equation was written in a specific basis, say right and left
circular. But it would have had the same form when using any other
orthogonal pair as basis since the source is unpolarized. Hence
we are free to apply this change of basis to one solution to get
another solution of Eq. (3). Under such a change, the
coefficients transform according to
(A.1) |
A more geometric view of this degeneracy is obtained when we think in terms of the Poincaré sphere representation of the states of polarization of all the feeds. The cross correlation between the outputs of two feeds, both of which receive unpolarized radiation, has a magnitude equal to the cosine of half the angle between the representative points on the sphere. Measurements of all such cross correlations with unpolarized radiation fixes the relative geometry of the points on the sphere, while leaving a two parameter degeneracy corresponding to overall rigid rotations of the sphere. This degeneracy can be lifted by the measurement of one polarized source at many parallactic angles.
Finally, we note that for the purpose of correcting the observations of unpolarized sources for the effects of non-identical feed polarization, the degeneracy is unimportant, because the correction factor is precisely the right hand side of Eq. (3) which is unaffected by all the transformations we have discussed.