No Arabic abstract
Since more than a century astronomers measure the position angle of the major axis of the polarization ellipse starting from the North direction and increasing counter-clockwise, when looking at the source. This convention has been enforced by the IAU with a Resolution in 1973. Much later the WMAP satellite, which has observed the polarization of the cosmic microwave background, has unfortunately adopted the opposite convention: the polarization position angle is measured starting from the South and increasing clockwise, when looking at the source. This opposite convention has been followed by most cosmic microwave background polarization experiments and is causing obvious problems and misunderstandings. The attempts and prospects to enforce the official IAU convention are described.
With the forthcoming release of high precision polarization measurements, such as from the Planck satellite, the metrology of polarization needs to improve. In particular, it is crucial to take into account full knowledge of the noise properties when estimating polarization fraction and angle, which suffer from well-known biases. While strong simplifying assumptions have usually been made in polarization analysis, we present a method for including the full covariance matrix of the Stokes parameters in estimates for the distributions of the polarization fraction and angle. We thereby quantify the impact of the noise properties on the biases in the observational quantities. We derive analytical expressions for the pdf of these quantities, taking into account the full complexity of the covariance matrix, including the Stokes I intensity components. We perform simulations to explore the impact of the noise properties on the statistical variance and bias of the polarization fraction and angle. We show that for low variations of the effective ellipticity between the Q and U components around the symmetrical case the covariance matrix may be simplified as is usually done, with negligible impact on the bias. For S/N on intensity lower than 10 the uncertainty on the total intensity is shown to drastically increase the uncertainty of the polarization fraction but not the relative bias, while a 10% correlation between the intensity and the polarized components does not significantly affect the bias of the polarization fraction. We compare estimates of the uncertainties affecting polarization measurements, addressing limitations of estimates of the S/N, and we show how to build conservative confidence intervals for polarization fraction and angle simultaneously. This study is the first of a set of papers dedicated to the analysis of polarization measurements.
With the forthcoming release of high precision polarization measurements, such as from the Planck satellite, it becomes critical to evaluate the performance of estimators for the polarization fraction and angle. These two physical quantities suffer from a well-known bias in the presence of measurement noise, as has been described in part I of this series. In this paper, part II of the series, we explore the extent to which various estimators may correct the bias. Traditional frequentist estimators of the polarization fraction are compared with two recent estimators: one inspired by a Bayesian analysis and a second following an asymptotic method. We investigate the sensitivity of these estimators to the asymmetry of the covariance matrix which may vary over large datasets. We present for the first time a comparison among polarization angle estimators, and evaluate the statistical bias on the angle that appears when the covariance matrix exhibits effective ellipticity. We also address the question of the accuracy of the polarization fraction and angle uncertainty estimators. The methods linked to the credible intervals and to the variance estimates are tested against the robust confidence interval method. From this pool of estimators, we build recipes adapted to different use-cases: build a mask, compute large maps, and deal with low S/N data. More generally, we show that the traditional estimators suffer from discontinuous distributions at low S/N, while the asymptotic and Bayesian methods do not. Attention is given to the shape of the output distribution of the estimators, and is compared with a Gaussian. In this regard, the new asymptotic method presents the best performance, while the Bayesian output distribution is shown to be strongly asymmetric with a sharp cut at low S/N.Finally, we present an optimization of the estimator derived from the Bayesian analysis using adapted priors.
This paper defines the mathematical convention adopted to describe an electromagnetic wave and its polarisation state, as implemented in the PSRCHIVE software and represented in the PSRFITS definition. Contrast is made between the convention that has been widely accepted by pulsar astronomers and the IAU/IEEE definitions of the Stokes parameters. The former is adopted as the PSR/IEEE convention, and a set of useful parameters are presented for describing the differences between the PSR/IEEE standard and the conventions (either implicit or explicit) that form part of the design of observatory instrumentation. To aid in the empirical determination of instrumental convention parameters, well-calibrated average polarisation profiles of PSR J0304+1932 and PSR J0742-2822 are presented at radio wavelengths of approximately 10, 20, and 40 cm.
Rotations of the electric vector position angle (EVPA) in blazars are often close to an integral multiple of 180$^circ$. There are multiple examples of this in the literature, and our analysis here, of the optical polarization data from the RoboPol monitoring program, strengthens the evidence by showing that $npi$ rotations occur more frequently than expected by chance. We explain this with a model consisting of two polarized emission components: a jet that is constant in time, and a burst that is variable. The EVPA of the combination is $rm EVPA_{jet}$ at both the beginning and the end of the burst, so the net rotation across the burst must be $npi$. Examples are analyzed on the Stokes plane, where the winding number for the Stokes vector of the combination gives the value of $n$. The main conclusion is that the EVPA rotation can be much larger than the physical rotation of the emission region around the axis of the jet, but this requires the EVPAs of the jet and the burst to be nearly orthogonal. A shock-in-jet calculation by Zhang et al. can provide a physical model for our toy model, and in addition automatically gives the needed orthogonality. The model is illustrated with data on OJ287 published by Myserlis et al., and we suggest that the large rapid EVPA rotation seen there might be a phase effect and not representative of a physical rotation.
We present a method of cross-calibrating the polarization angle of a polarimeter using BICEP Galactic observations. bicep was a ground based experiment using an array of 49 pairs of polarization sensitive bolometers observing from the geographic South Pole at 100 and 150 GHz. The BICEP polarimeter is calibrated to +/-0.01 in cross-polarization and less than +/-0.7 degrees in absolute polarization orientation. BICEP observed the temperature and polarization of the Galactic plane (R.A= 100 degrees ~ 270 degrees and Dec. = -67 degrees ~ -48 degrees). We show that the statistical error in the 100 GHz BICEP Galaxy map can constrain the polarization angle offset of WMAP Wband to 0.6 degrees +- 1.4 degrees. The expected 1 sigma errors on the polarization angle cross-calibration for Planck or EPIC are 1.3 degrees and 0.3 degrees at 100 and 150 GHz, respectively. We also discuss the expected improvement of the BICEP Galactic field observations with forthcoming BICEP2 and Keck observations.