No Arabic abstract
Heterodyne receivers register the sky signal on either a circular polarization basis (where it is split into left-hand and right-hand circular polarization) or a linear polarization basis (where it is split into horizontal and vertical linear polarization). We study the problem of interferometric observations performed with telescopes that observe on different polarization bases, hence producing visibilities that we call mixed basis (i.e., linear in one telescope and circular in the other). We present novel algorithms for the proper calibration and treatment of such interferometric observations and test our algorithms with both simulations and real data. The use of our algorithms will be important for the optimum calibration of forthcoming observations with the Atacama Large mm/submm Array (ALMA) in very-long-baseline interferometry (VLBI) mode. Our algorithms will also allow us to optimally calibrate future VLBI observations at very high data rates (i.e., wide bandwidths), where linear-polarization feeds will be preferable at some stations, to overcome the polarimetric limitations due to the use of quarter-wave plates.
We report the development of a semi-automatic pipeline for the calibration of 86 GHz full-polarization observations performed with the Global Millimeter-VLBI array (GMVA) and describe the calibration strategy followed in the data reduction. Our calibration pipeline involves non-standard procedures, since VLBI polarimetry at frequencies above 43 GHz is not yet well established. We also present, for the first time, a full-polarization global-VLBI image at 86 GHz (source 3C 345), as an example of the final product of our calibration pipeline, and discuss the effect of instrumental limitations on the fidelity of the polarization images. Our calibration strategy is not exclusive for the GMVA, and could be applied on other VLBI arrays at millimeter wavelengths. The use of this pipeline will allow GMVA observers to get fully-calibrated datasets shortly after the data correlation.
New and upcoming radio interferometers will produce unprecedented amounts of data that demand extremely powerful computers for processing. This is a limiting factor due to the large computational power and energy costs involved. Such limitations restrict several key data processing steps in radio interferometry. One such step is calibration where systematic errors in the data are determined and corrected. Accurate calibration is an essential component in reaching many scientific goals in radio astronomy and the use of consensus optimization that exploits the continuity of systematic errors across frequency significantly improves calibration accuracy. In order to reach full consensus, data at all frequencies need to be calibrated simultaneously. In the SKA regime, this can become intractable if the available compute agents do not have the resources to process data from all frequency channels simultaneously. In this paper, we propose a multiplexing scheme that is based on the alternating direction method of multipliers (ADMM) with cyclic updates. With this scheme, it is possible to simultaneously calibrate the full dataset using far fewer compute agents than the number of frequencies at which data are available. We give simulation results to show the feasibility of the proposed multiplexing scheme in simultaneously calibrating a full dataset when a limited number of compute agents are available.
The detection of the primordial $B$-mode spectrum of the polarized cosmic microwave background (CMB) signal may provide a probe of inflation. However, observation of such a faint signal requires excellent control of systematic errors. Interferometry proves to be a promising approach for overcoming such a challenge. In this paper we present a complete simulation pipeline of interferometric observations of CMB polarization, including systematic errors. We employ two different methods for obtaining the power spectra from mock data produced by simulated observations: the maximum likelihood method and the method of Gibbs sampling. We show that the results from both methods are consistent with each other, as well as, within a factor of 6, with analytical estimates. Several categories of systematic errors are considered: instrumental errors, consisting of antenna gain and antenna coupling errors, and beam errors, consisting of antenna pointing errors, beam cross-polarization and beam shape (and size) errors. In order to recover the tensor-to-scalar ratio, $r$, within a 10% tolerance level, which ensures the experiment is sensitive enough to detect the $B$-signal at $r=0.01$ in the multipole range $28 < ell < 384$, we find that, for a QUBIC-like experiment, Gaussian-distributed systematic errors must be controlled with precisions of $|g_{rms}| = 0.1$ for antenna gain, $|epsilon_{rms}| = 5 times 10^{-4}$ for antenna coupling, $delta_{rms} approx 0.7^circ$ for pointing, $zeta_{rms} approx 0.7^circ$ for beam shape, and $mu_{rms} = 5 times 10^{-4}$ for beam cross-polarization.
Calibration of radio interferometric observations becomes increasingly difficult towards lower frequencies. Below ~300 MHz, spatially variant refractions and propagation delays of radio waves traveling through the ionosphere cause phase rotations that can vary significantly with time, viewing direction and antenna location. In this article we present a description and first results of SPAM (Source Peeling and Atmospheric Modeling), a new calibration method that attempts to iteratively solve and correct for ionospheric phase errors. To model the ionosphere, we construct a time-variant, 2-dimensional phase screen at fixed height above the Earths surface. Spatial variations are described by a truncated set of discrete Karhunen-Loeve base functions, optimized for an assumed power-law spectral density of free electrons density fluctuations, and a given configuration of calibrator sources and antenna locations. The model is constrained using antenna-based gain phases from individual self-calibrations on the available bright sources in the field-of-view. Application of SPAM on three test cases, a simulated visibility data set and two selected 74 MHz VLA data sets, yields significant improvements in image background noise (5-75 percent reduction) and source peak fluxes (up to 25 percent increase) as compared to the existing self-calibration and field-based calibration methods, which indicates a significant improvement in ionospheric phase calibration accuracy.
The daytime sky has been recently demonstrated as a useful calibration tool for deriving polarization cross-talk properties of large astronomical telescopes. The Daniel K Inouye Solar Telescope (DKIST) and other large telescopes under construction can benefit from precise polarimetric calibration of large mirrors. Several atmospheric phenomena and instrumental errors potentially limit the techniques accuracy. At the 3.67m AEOS telescope on Haleakala, we have performed a large observing campaign with the HiVIS spectropolarimeter to identify limitations and develop algorithms for extracting consistent calibrations. Effective sampling of the telescope optical configurations and filtering of data for several derived parameters provide robustness to the derived Mueller matrix calibrations. Second-order scattering models of the sky show that this method is relatively insensitive to multiple-scattering in the sky provided calibration observations are done in regions of high polarization degree. The technique is also insensitive to assumptions about telescope induced polarization provided the mirror coatings are highly reflective. Zemax-derived polarization models show agreement between the functional dependence of polarization predictions and the corresponding on-sky calibrations.