(3) |
(4) |
By solving for antenna-based errors, the closure phase on a given triangle ijk is preserved, which implies the constraint:
< E_{i} -E_{j}>+ < E_{j} -E_{k}>- < E_{i} -E_{k}> =0 . | (5) |
We can see that by substituting Eq. (3) and Eq. (4) in Eq. (5), we have:
(6) | ||
(where hereafter we omit the baseline subscript for the visibility phases.)
By defining and the observed and the model closure phases (Eq. 2), we have:
(7) |
Where we omit the subscript "ijk'' for the closure phases.
In other words the estimation of the 's depend on
the closure phase. How the difference between the observed and model
closure phases is distributed among the individual errors it depends upon the weights attributed to the various baselines.
Since we want here to illustrate in general the role played by the closure
phase in the solution, we assume for simplicity equal weights, and
estimations of the errors like:
(8) |
(9) |
In other words a small closure phase gives rise to an erroneously low value for the estimated error in the model. This wrong estimation of the error in the model propagates itself into a wrong estimation of the antenna error,
(10) |
which becomes since and both vanish. As a result the corrected observed phase is equal to
(11) |
Let us now assume that the model phase determined by the set of clean components, produces, due to the use of proper CLEAN boxes, an smaller than before and more comparable to the small closure phase. In this case self calibration should converge toward the proper solution, i.e., each iteration should remove the bias of the preceding one. Why, then, does this not happen? In a linear process such as described in Eq. (10), the solution derived from a given iteration of self calibration (sc) is simply added to the solution from the previous iteration, that is
(12) | ||
If the self calibration process were indeed linear, Eq. (12) would produce the same result as from a single iteration of Loop "B'':
(13) |
But we have shown empirically that in the case of Test-2
(14) |
This implies a non-linear equation for .
As we will show in Sect. 4.2 self calibration is in fact a non linear process. This explains why Loop "B'' works better than Loop "A''. At each iteration Loop "B'' discards the previous (E_{i}-E_{k}) estimation and simply starts again, approaching the proper solution for the antenna corrections as .On the contrary Loop "A'' accumulates previous, incorrect estimations and being unable to correct for them, biases the total antenna corrections even though .We have checked this surmise by setting CLEAN boxes to exclude the spurious features, forcing the CLEAN model to be the "true" one. Self calibration nevertheless determines corrections which still reproduce the artificial counter-jet.
However, if we accept that in general Loop "B'' works better than Loop "A'' since self calibration is non-linear, why does, in some cases, Loop "A'' converge to the same solution as Loop "B''? There must be some linearizing quantity in Test-1 that is lacking in Test-2. We remind the reader that Test-2 has the same u-v coverage as Test-1, except for the single intermediate baseline.
The presence of
large errors in the model violates the condition of linearity;
we will see below quantitatively what "large" means.
Let us consider in Eq. (1) only phase errors, ,
that the gain terms g_{i} = g_{k} = 1 and neglect the amplitude and
weight terms.
Then, minimizing Eq. (1) with respect to the
phase error terms gives
(15) |
Equation (15) then becomes
, or
(16) |
This implies that for the linear case, that is for Eq. (16) to be equal to the linear case expressed in Eq. (8), it must hold true that .
Assuming an initial point-like source, the condition
for linearity is further reduced to
(17) |
Since the errors for a point-like model are exactly the phase of our simulated source from FAKE, we can illustrate this result.
Figure 9a shows the closure phase available in Test-1 along with the combination of the errors. The closure phase is in this case more than one order of magnitude higher than the combination of the errors; these last appear in fact quite unsignificant in the plot. This implies that Eq. (17) is for Test-1 satisfied and as a consequence self calibration should reduce to the linear case, that is Eq. (12) should be satisfied as well. And in fact we have shown as Loop "A'' give same results as Loop "B'' for Test-1.
On the contrary we see in Fig. 9b that when the closure phase drops to a few degrees, as is the case for unbalanced triangles, the errors in the model, previous negligible, become a significant fraction of the closure phase. In other words the condition of Eq. (17) no longer applies and the non linear terms in self calibration cannot be neglected; we are therefore in the condition expressed in Eq. (14). A conclusion verified in Sect. 2 with test2 which gives different results for Test-2A and Test-2B.
Figure 9c illustrates that when the source has a complicated structure and a point source has been used as an initial model, the errors are so large to become an important fraction of the closure phase even for large closure phase values. In the case of Fig. 9d finally we show the worst combination of small closure phase values and large errors giving rise to a solution dominated by the non linear terms.
Figure 9: Closure phase (continuous line) and non-linear error terms (crosses) in the self calibration solution (see Sects. 4 and 5) |
Copyright The European Southern Observatory (ESO)