Free Access
Issue
A&A
Volume 600, April 2017
Article Number A63
Number of page(s) 8
Section Astronomical instrumentation
DOI https://doi.org/10.1051/0004-6361/201629571
Published online 31 March 2017

© ESO, 2017

1. Introduction

Linearly polarised emission is a powerful tool in astrophysics. Scattering of light generates optical polarisation that can constrain the geometry of reflection nebulae (e.g. Scarrott et al. 1986). Photons can also be polarised by scattering or extinction at elongated dust grains aligned in interstellar magnetic fields (e.g. Hoang & Lazarian 2014), which allows mapping these magnetic fields in the Milky Way (Fosalba et al. 2002) and in the Small Magellanic Cloud (Gomes et al. 2015). Elongated dust grains emit polarised emission at sub-mm wavelengths, which is useful to study the magnetic fields in molecular clouds (e.g. Tang et al. 2009; Pillai et al. 2015) or the halos of galaxies (e.g. Greaves et al. 2000). Recently, the Planck mission provided all-sky dust polarisation maps revealing large-scale magnetic fields in the Milky Way (Planck Collaboration Int. XIX 2015).

Synchrotron emission is up to 75% linearly polarised, with its B-vector intrinsically parallel to the magnetic field. Optical synchrotron emission allows the investigation of magnetic fields of jets emerging from galactic nuclei (e.g. Perlman et al. 2011) and radio synchrotron emission the investigation of magnetic fields in the Milky Way (Wolleben et al. 2006), other spiral galaxies (e.g. Beck 2016), and radio galaxies (e.g. Laing & Bridle 2014).

Linearly polarised emission is usually described in terms of the Stokes parameters U and Q, defined as U = Psin(2 χ) and Q = Pcos(2 χ), where P is the polarised intensity and χ is the polarisation angle. If the receiving system delivers orthogonally polarised signals with amplitudes X and Y, U and Q are computed as U = 2 XYcos δ (where δ is the phase between X and Y) and Q = X2Y2. If circularly polarised signals with amplitudes L and R are delivered by the receiving system, U = 2 LRsin δ and Q = 2 LRcos δ, where δ is the phase between L and R. Low-frequency radio telescopes consisting of dipoles, like the LOw Frequency ARray (LOFAR), deliver projections of X and Y signals, whereas most higher-frequency radio telescopes use correlating devices to transform L and R signals into U and Q.

Measuring polarised signals suffers from a fundamental problem. In addition to the true polarised signal PT of a source, root mean squared (rms) noise obtained from the receiving system is also detected. If we use the standard formula to calculate \hbox{$\hat{P}$} from the measured data Û and \hbox{$\hat{Q}$}, \hbox{$\hat{P}$} can be expressed by the polarised components of the source UT and QT and their noise contributions NU and NQ: =2+2=(UT+NU)2+(QT+NQ)2.\begin{equation} \begin{array}{lcl} \hat{P} &=& \sqrt{\hat{U}^2 + \hat{Q}^2} \\ &=& \sqrt{(U_T+N_U)^2+(Q_T+N_Q)^2} . \end{array} \end{equation}(1)The noise always delivers a positive bias to the true polarised intensity PT that cannot be separated out for small signal-to-noise ratios s (\hbox{$s = \hat{P}/\sigma$}, where σ is the rms noise in the maps of Û and \hbox{$\hat{Q}$}). Though this could be overcome by clipping maps of polarised intensity below a certain value, the bias accumulates by integration of polarised intensities and, if it cannot be effectively suppressed, prevents the determination of reliable flux densities in polarisation. Further, due to the noise bias, the distribution of \hbox{$\hat{P}$} is Ricean. As a result, further data processing, for example smoothing directly to a larger beam size (this is usually done instead by smoothing the Û and \hbox{$\hat{Q}$} maps, which is subject to depolarisation).

The distribution of \hbox{$\hat{P}$} is Ricean for small s and becomes Gaussian for large s. Vinokur (1965, Eqs. (81) and (82)) showed that a separation of signal PT and standard deviation of the noise distribution σ is possible only for large s and the expectation value is: P=PT+σ2/(2PT),\begin{equation} \langle P \rangle = P_T + \sigma^2/(2 \,P_T) , \end{equation}(2)while for small sP=π/2σ(1+s2/4).\begin{equation} \langle P \rangle = \sqrt{\pi/2} \, \sigma \, (1 + s^2/4) . \end{equation}(3)Several other methods have been developed to correct for the noise bias. The most widely used method is that of Wardle & Kronberg (1974) who proposed the following bias correction for \hbox{$\hat{P}$} maps: P=2σ2\begin{equation} P^* = \sqrt{\hat{P}^2 - \sigma^2} \end{equation}(4)where σ is the rms standard deviation of the noise distributions in Û and \hbox{$\hat{Q}$}.

An astrophysical interpretation requires that a map of bias-corrected P has a base level of about zero in regions containing only noise. This requires introduction of negative values of P at locations where the true signal P is very weak: P=2(Cσ)2(if2(Cσ)2),P=2+(Cσ)2(if2<(Cσ)2).\begin{equation} \begin{array}{lcl} P^* &=& \sqrt{\hat{P}^2 - (C\, \sigma)^2} \ \ (\mathrm{if} \,\, \hat{P}^2 \ge (C\, \sigma)^2) \, ,\\ P^* &=& - \sqrt{- \hat{P}^2 + (C\, \sigma)^2} \ \ (\mathrm{if} \,\, \hat{P}^2 < (C\, \sigma)^2) \, . \end{array} \end{equation}(5)The smallest possible value of P = −Cσ occurs for \hbox{$\hat{U} + \hat{Q} = 0$}. The factor C was introduced later to adjust the bias correction: C = 1.2 is used in the NOD2 software package (Haslam 1974; Andernach 1985) and AIPS (Greisen 2003) uses C = 1.253 in the option POLC of the task COMB. If \hbox{$\hat{P}$} is determined from multi-channel data with help of rotation measure (RM) synthesis, the additional uncertainty in RM requires that C is increased to 1.5 (George et al. 2012).

The AIPS task POLCO uses the maximum likelihood correction introduced by Killeen et al. (1986) for large signal-to-noise ratios (s> 2) and the Wardle & Kronberg method for smaller s.

Simmons & Stewart (1985) discussed various estimators of the true signal P, including a maximum likelihood, a median estimator, and the Wardle & Kronberg estimator. All estimators agree asymptotically for large s. All estimators yield a positive residual bias (relative to P) for s< 1, but a negative for 1 <s< 4 and hence are not a good representation of the true signal for s< 4. Montier et al. (2015) introduced a Bayesian estimator that provides a very low relative bias for s< 1.2, allowing reliable estimates of P in regions with low s, while the maximum likelihood estimator performs best for s> 2. However, none of the methods works well at small s and large s simultaneously.

If the polarisation angle is perfectly well known, the estimator proposed by Vidal et al. (2016) is able to completely correct the polarisation bias. This method can be applied to regions where the polarisation angle is expected to be constant, such as for large-scale magnetic fields as observed in sub-mm dust emission with WMAP and Planck.

The new method proposed here works for all values of s and does not need prior knowledge of the polarisation angle. All methods to compute P maps require that the base levels in the maps of Û and \hbox{$\hat{Q}$} are about zero (i.e. their mean values are smaller than about 20% of the rms noise values) in regions without sources. Baseline shifts may remain even after processing and combining single dish maps.

2. The method

When detecting astronomical signals at the telescope, statistical receiver noise is added to the U and Q signals. Any linear and non-linear combination of the noisy Û and \hbox{$\hat{Q}$} signals should take into account the probability density of the noise distribution in the resulting data. If we could mathematically express noise and the signal as a linear combination of both, it is possible to deal with the noise separately. This is not the case if we compute the square root of U2 plus Q2.

thumbnail Fig. 1

Geometrical sketch showing the noise components of the polarisation vectors. The noise vector \hbox{$N_{\hat{P}}$} can be split in two different orthogonal components. Statistically, all NU and NQ contain the measured noise distributions, so the other two projected noise components calculated by the projection of \hbox{$\hat{P}$} onto PT, respectively Pm, contain the same distribution.

We express the polarised intensity P in a different way. For each pixel, we convert U and Q to polar coordinates and define the geometrical angle θ = tan-1(U/Q) (twice the polarisation angle), so that: U=Psin(θ),Q=Pcos(θ),P=Usin(θ)+Qcos(θ).\begin{equation} \begin{array}{lcl} U &=& P\sin(\theta) \, , \\ Q &=& P\cos(\theta) \, , \\ P &=& U\sin(\theta)+ Q\cos(\theta) \, . \end{array} \end{equation}(6)It is then possible to separate the true signals UT, QT and PT from their noise contributions NU and NQ. The noise contribution of NP depends on the method applied and is given by the projection of NU and NQ if the true angle θT is known: PT+NP=(UT+NU)sin(θT)+(QT+NQ)cos(θT),θT=tan-1(UT/QT),\begin{equation} \begin{array}{rcl} P_T+N_P &=& (U_T+N_U)\sin(\theta_T) + (Q_T+N_Q)\cos(\theta_T) \, ,\\ \theta_T &=& \tan^{-1}(U_T/Q_T) \, , \label{pinew} \end{array} \end{equation}(7)where θT is the true angle in the polar representation of the noise-free signals UT and QT, as illustrated in Fig. 1. The observed quantities are defined as Û = UT + NU, \hbox{$\hat{Q} = Q_T + N_Q$}, so that \hbox{$\hat{\theta} = \tan^{-1}(\hat{U}/\hat{Q})$} and =2+2\hbox{$\hat{P} = \sqrt{\hat{U}^2 + \hat{Q}^2\ }$}. Equation (7) can also be expressed as the projection of the vector \hbox{$\hat{P}$} onto the direction of the true PT as follows: PT+NP=cos(θTθ̂)=cos(Δθ),\begin{equation} P_T + N_P = \hat{P}\cos(\theta_T - \hat{\theta}) = \hat{P}\cos(\Delta\theta) , \label{piproj} \end{equation}(8)which can easily be proven by expanding \hbox{$\cos(\theta_T - \hat{\theta})$} and comparing with \hbox{$\hat{U}\sin(\theta_T) + \hat{Q}\cos(\theta_T)$}.

Application of this equation requires the knowledge of the true angle θT = tan-1(UT/QT) of the signal. This is the crucial step for this method. For an individual pair of measured Û and \hbox{$\hat{Q}$} it is impossible to know their noise contributions. For a number n of single observations on a source, or a map of n pixels around the source, we reduce the angle error θnθT significantly by applying an averaging filter to adjacent pixels. We choose the median filter, because it reduces the bias approximately the square root of the number of the averaging pixels, and gives a reliable value for sudden changes of the polarisation angles at the adjacent pixels. However, a problem occurs with the discontinuity of the angle θ from –π to π. Therefore the components x = cos(θ) and y = sin(θ) are taken for the median filtering and the arctan2(ym,xm) is used to calculate the angle θ from the medians of x and y. The adjacent pixels are not independent from each other due to the smoothing of the telescope beam. The pixel size is usually chosen to be one third of the beam size, for which a median filter size of 5 × 5 pixels is a good choice. In this case the angle error Δθm of the median value θm is reduced by a factor of approximately five.

We have to take into account that the angle \hbox{$\hat{\theta} = \tan^{-1}(\hat{U}/\hat{Q})$} represents the bias at the centre position of the filter box, which leads to an undercorrection of the bias. Therefore we remove this angle at the centre position of the filter box to suppress this bias again at noise level (“modified median filter”). However, disregarding this pixel leads to an overcorrection. In Fig. 2 we used the weighted average (1:2) of the median and the modified median of the filter box, which is an empirical choice to optimise the bias correction.

In case of overlapping sources with different polarisation angles, there is a depolarisation effect due to the beam smoothing. The application of the modified median filter may increase the error in P. We note that the averaging of the angle does not affect the angular resolution of the P map, other than smoothing the Û and \hbox{$\hat{Q}$} maps.

From the observations, we get Û = UT + NU and \hbox{$\hat{Q}=Q_T+N_Q$} and the modified median filtered angle θm, and we can separate Pm and its noise contribution NPm. The term Pm + NPm is the resulting polarised intensity P determined with our new method, calculated with the angle θm that should be very close to the true angle θT: P=Pm+NPm=(UT+NU)sin(θm)+(QT+NQ)cos(θm)=\begin{eqnarray} P^* &=& P_{\rm m} + N_{P_{\rm m}} = (U_T+N_U)\sin(\theta_{\rm m}) + (Q_T+N_Q)\cos(\theta_{\rm m}) \nonumber\\ &=& \hat{U}\sin(\theta_{\rm m}) + \hat{Q}\cos(\theta_{\rm m}) , \label{piobs} \end{eqnarray}(9)where θmθT. According to Eq. (8) the noise-free PT is somewhat different from Pm: PT+δP=Pmcos(θTθm).\begin{equation} P_T + \delta{P} = P_{\rm m}\cos(\theta_T - \theta_{\rm m}) . \end{equation}(10)The assumption is that the median angle θm is very close to the true angle θT, so |δθ| = |θTθm| ≪ 1, and sin(δθ) <δθ.

The accuracy of the determination of P depends directly on the angle error δθ. Therefore we can express the error in Pm as a function of δθ: δP(δθ)=Pmθmcos(θTθm)δθ=\begin{eqnarray} \delta{P(\delta{\theta)}} &=& P_{\rm m}\ \frac{\partial{}}{\partial{\theta_{\rm m}}} \cos(\theta_T-\theta_{\rm m}) \ \delta{\theta}\nonumber \\ &=& P_{\rm m}\sin(\theta_T-\theta_{\rm m}) \ \delta{\theta}< P_{\rm m}\,(\delta{\theta})^2 \, . \label{eqnoise} \end{eqnarray}(11)On the other hand, the angle error Δθ generally depends on the errors of U and Q: (Δθ(U,Q))2=[∂Utan-1(U/Q)ΔU]2+[∂Qtan-1(U/Q)ΔQ]2=\begin{eqnarray} (\Delta{\theta}(U,Q))^2 &=& \left[\frac{\partial{}}{\partial{U}} {\tan^{-1}(U/Q)}\ \Delta{U}\right]^2 + \left[\frac{\partial{}}{\partial{Q}} {\tan^{-1}(U/Q)}\ \Delta{Q}\right]^2\nonumber \\ &=& [(Q\,\Delta{U})^2 + (U\,\Delta{Q})^2]\,/\,(U^2 + Q^2)^2 \, . \label{dP_theta} \end{eqnarray}(12)In addition, we have to take into account the contribution of the modified median filter to the angle error. If we apply an n × n filter box we get δθm(U,Q)Δθ(U,Q)/n.\begin{equation} \delta{\theta}_{\rm m}(U,Q) \approx \Delta{\theta}(U, Q)/n \, . \end{equation}(13)Using sin(δθm) <δθm, the error of Pm with respect to the angle error depending on ΔU and ΔQ is: δP(δθm)=Pmsin(δθm)δθm<[(QΔU)2+(UΔQ)2]/[n2Pm3].\begin{equation} \begin{array}{lcl} \delta{P}(\delta{\theta_{\rm m}}) &=& P_{\rm m}\sin(\delta{\theta_{\rm m}})\ \delta{\theta_{\rm m}}\\ &&< [(Q\,\Delta{U})^2 + (U\,\Delta{Q})^2]\,/\,[n^2\,P_{\rm m}^3] \, . \end{array} \end{equation}(14)The higher the signal-to-noise ratio, the smaller the angle error δθm. In the case where ΔU ≈ ΔQ, δP(δθm)<(ΔUΔQ)/(n2Pm),\begin{equation} \delta{P}(\delta{\theta_{\rm m}}) < (\Delta{U}\,\Delta{Q})/(n^2 \ P_{\rm m}) \, , \end{equation}(15)and we can neglect the error of P with respect to the angle error δθm in all cases.

The above section explained that the contribution of Δθ to the error of ΔP is negligible. Thus the error depends only on ΔU and ΔQ, following Eq. (9): (ΔP(U,Q))2=sin2(θm)(ΔU)2+cos2(θm)(ΔQ)2.\begin{equation} (\Delta{P}(U, Q))^2 = \sin^2(\theta_{\rm m})\,(\Delta{U})^2 + \cos^2(\theta_{\rm m})\,(\Delta{Q})^2 \, . \end{equation}(16)If U and Q are obtained from correlated signals, then ΔU = ΔQ, and ΔP = ΔU = ΔQ.

According to Fig. 2 the corrected polarised intensity P is always overestimated when no bias correction is applied, whereas the Wardle & Kronberg method always underestimates P. Our new method recovers the true signal P more precisely, except for very small signal-to-noise ratios.

thumbnail Fig. 2

Averaged signal-to-noise ratio (P/σ) as a function of the true ratio PT/σ after applying: the bias suppression method; the Wardle & Kronberg method (using C = 1.2); our new method; and without any bias correction. To show the improvement of the bias suppression with our new method, we created a 100 × 100 pixel Û map and \hbox{$\hat{Q}$} map that contained only Gaussian noise. We added a constant positive value A to the entire Û map and a negative value A to the \hbox{$\hat{Q}$} map. The values vary between 0 < | A | < 3σ. For each A we calculated the corresponding P as the average value over the entire map for all methods given in the plot.

thumbnail Fig. 3

Polarised intensity reference map. P is calculated from the noise-free U and Q maps using P=U2+Q2\hbox{$P = \sqrt{U^2 + Q^2\ }$}. The lines indicate the orientation of the polarised emission.

thumbnail Fig. 4

Polarised intensity map calculated with our new method.

Our method is fundamentally different from that of Wardle & Kronberg. We are able to determine the polarised intensity, including the noise, by projecting the observed values Û and \hbox{$\hat{Q}$} onto the direction of the median filtered polarisation angle θm (Eq. (9)). The noise distribution of P is the same as that of U und Q (Fig. 7) and hence averages to zero over large areas. In contrast to our method, Wardle & Kronberg use the mean noise (rms) of the entire map to correct P for the bias. We do not need to know the mean errors of Û and \hbox{$\hat{Q}$}. Furthermore, the distributions of these errors do not necessarily need to be Gaussian.

Systematic errors (in addition to pure noise), such as fluctuations of the background emission level in Stokes U and Q on the scales of the beam or instrumental polarisation (usually below 1% of the total emission), lead to additional signals in polarised intensity that cannot be corrected by the methods of bias suppression. Instrumental polarisation should be reduced by applying the Mueller matrix to U and Q before the determination of P.

3. Numerical simulations

To test our new method, simulated images are created that contain an artificial source (e.g. a galaxy) and four box shaped structures of size 40 × 33 pixels. These simulate abrupt changes at the edges and constant polarised intensity inside the boxes. The artificial galaxy is composed of two two-dimensional Gaussians, a circular one at the central region and a large elliptical one simulating the halo. The simulated Stokes U contains only the halo, whereas the Q map contains both a halo of the artificial galaxy component and the centre source. The boxes have different amplitudes of ± 0.5·rms, ± 1.5·rms, ± 3.0·rms, and ± 5.0·rms, in which the component of Q is always positive and U always negative. The amplitudes of the boxes increases from the upper left anticlockwise to the upper right.

To test the effect of different polarisation angles on the scale of the beam we add three overlapping Gaussian sources with amplitudes of + 5·rms, 2·rms, −5·rms to Q and + 2·rms, −5·rms, + 2·rms to U. We note that in this case we get depolarisation in the overlapping regions of the sources. From the two maps of U and Q the polarised intensity is constructed as P=U2+Q2\hbox{$P = \sqrt{U^2 + Q^2}$}, which is the noise-free reference map (Fig. 3).

Gaussian noise NU and NQ with rms = 5 is added to both the U and Q maps. Subsequently the bias correction performed with our method by applying the modified median filter (MMF) to the angle map θm = MMF(tan-1 [ (Ui + NU) / (Qi + NQ) ]). Then the bias corrected polarised intensity P is computed using Eq. (9) which is shown in Fig. 4.

After subtracting the reference map, the residual map should contain only noise, if the method has worked properly. This is the case for our method (Fig. 5) but not for the method of Wardle & Kronberg (Fig. 6). We note that even in the case of three overlapping sources with different polarisation angles our method works well, due to the behaviour of the median filter. The histogram (Fig. 7) of the residual map shows Gaussian noise distribution with rms = 5 and centred close to zero. The rms is the same that was added to U and Q, as expected. Compared with the Wardle & Kronberg method (Fig. 8) there is a significant improvement. At low intensities the Wardle & Kronberg method shows artefacts which are not seen in our new method.

thumbnail Fig. 5

Residual map from our new method. Only Gaussian noise is showing after subtracting the reference map.

thumbnail Fig. 6

Residual map from the Wardle & Kronberg method. Artefacts are showing at low intensities.

thumbnail Fig. 7

Histogram of the residual map from our new method (Fig. 5). Intensities are given in arbitrary numbers.

thumbnail Fig. 8

Histogram of the double-peaked residual map from the Wardle & Kronberg method (Fig. 6). Intensities are given in arbitrary numbers.

Table 1 shows the results of the two methods, comparing the total integral of the entire maps. The relative error of our new method is significantly smaller than the relative error of the Wardle & Kronberg method. Also the relative errors of the boxes of our new method are smaller than those of the Wardle & Kronberg method (Table 2).

Table 1

Sum of polarised intensities of the simulated images over all pixels.

Table 2

Sum of polarised intensities of the simulated images over each box.

4. Application to real data: M 31

The radio data of the Andromeda galaxy M 31 (Fig. 9), observed with the Effelsberg 100-m telescope at 4.85 GHz (Berkhuijsen et al. 2003), can be used to effectively demonstrate the advantages of our new method. The maps have 3′ resolution and a pixel size of 1′. The polarisation map determined with our new method (Fig. 10) agrees with that obtained with the Wardle & Kronberg method (Fig. 11) at medium and high intensities, but the maps differ significantly at low intensities. The histogram for the new map is smooth (Fig. 12), whereas the histogram of the map determined with the Wardle & Kronberg method (Fig. 12) shows a sharp minimum at P = 0 and an increase towards negative values of P, which is the reason for the deep minima in Fig. 11. The lowest value in the P map created with our new method (Fig. 11) is –740 mJy/beam, while the map using the Wardle & Kronberg method is limited to about –293 mJy/beam. Using our new method, the integrated polarised flux density is 14% larger than that of the Wardle & Kronberg method.

thumbnail Fig. 9

Total intensity of M 31 at 4.85 GHz (in mJy/beam), observed with the Effelsberg 100-m telescope at an angular resolution of 3′ and 1′ pixel size. Bright background sources have been subtracted (from Berkhuijsen et al. 2003).

thumbnail Fig. 10

Polarised intensity of M 31 at 4.85 GHz (in mJy/beam) determined with our method and an angular resolution of 3′.

thumbnail Fig. 11

Polarised intensity of M 31 at 4.85 GHz (in mJy/beam) determined with the Wardle & Kronberg method, assuming a 0.2 mJy/beam rms noise for U and Q.

thumbnail Fig. 12

Histograms of polarised intensity of M 31 at 4.85 GHz: our new method (red) and the Wardle & Kronberg method (green).

The radio polarisation data of the Andromeda galaxy M 31 observed with the Very Large Array (VLA) at 1.4 GHz (Beck et al. 1998) reveals “depolarisation canals” at locations where the values of Û and \hbox{$\hat{Q}$} both cross the zero level and the polarisation angle jumps by 90° across the canals. The width is observed to be smaller than the telescope beamwidth, which needs steep gradients in Û and \hbox{$\hat{Q}$}. Such canals could be produced either by differential Faraday rotation (Shukurov & Berkhuijsen 2003) or by beam depolarisation in a turbulent magneto-ionic medium (Haverkorn & Heitsch 2004). The length, depth, and separation of the canals provide tools for studies of interstellar turbulence. The widths of the canals at half intensity, 60″−100″ in the map obtained with the Wardle & Kronberg method (Fig. 14), are 20%−30% larger in the new map (Fig. 13). This relaxes the need for steep gradients in Faraday rotation measures and hence reduces the amplitude needed for studies of interstellar turbulence. The difference between the maps of Figs. 13 and 14 is particularly significant in the canals and other regions of low intensity (Fig. 15). The canals are less deep when applying our method.

Equation (7) describes the projection of NU and NQ to NP. We then show that the noise distribution of NP has the same characteristics as NU and NQ. Therefore any other projections of two orthogonal noise components that result in NP also have the same characteristics. The projection described by Eq. (8) can be used to estimate the noise contribution NPm. The noise biased \hbox{$\hat{P}$} is the vector addition of the true polarised intensity PT and the noise contribution \hbox{$N_{\hat{P}}$} of Û and \hbox{$\hat{Q}$}. \hbox{$N_{\hat{P}}$} is also projected onto the direction of PTPm. Although the cosine component P = Pm + NPm does not allow us to separate the noise contribution, the sine component NP\hbox{$N_P'$} is due to pure noise: NP=sin(ΔΘ),whereΔΘ=θmθ̂,\begin{equation} N_P' = \hat{P}\sin(\Delta{\Theta})\, , \mbox{\ \ where \ \ } \Delta{\Theta} = \theta_{\rm m} - \hat{\theta} \, , \\ \label{error} \end{equation}(17)which is shown in Fig. 16. This equation can be used as an estimate of the noise distribution of the P map created with our method. It does, however, not include systematic errors arising from deviations between θm and θT. These become relevant for very small signal-to-noise ratios (the remaining bias is shown in Fig. 2) and in cases of strong small-scale variations of θ, in which the modified median cannot be a good representation of θm. In the former case, the statistics according to Eq. (17) still represent the noise well but do not represent residual systematic errors.

thumbnail Fig. 13

Polarised intensity of M 31 at 1.4 GHz (in mJy/beam) determined with our method and an angular resolution of 3′.

thumbnail Fig. 14

Polarised intensity of M 31 at 1.4 GHz (in mJy/beam) determined using the Wardle & Kronberg method and assuming an rms noise of 0.2 mJy/beam.

thumbnail Fig. 15

Difference map of our method (Fig. 13) and the Wardle & Kronberg method (Fig. 14).

The noise components \hbox{$\hat{P}_{\rm m}\cos(\Delta{\Theta}) - P_{\rm m}$} and \hbox{$\hat{P}\sin(\Delta{\Theta})$} of \hbox{$N_{\hat{P}}$} (as shown in Fig. 16) have the same statistical characteristics as the noise components NU and NQ of the same vector \hbox{$N_{\hat{P}}$}. Statistically, both noise components are equally distributed.

thumbnail Fig. 16

Histogram of the probability density of the noise distribution taken from Eq. (17).

This enables us to determine the probability density of the noise from the error of the angles between the observed angle \hbox{$\hat{\theta} = \tan^{-1}(\hat{U}/\hat{Q})$} and the mean angle θm, which should be close to the true angle θT. Figure 16 shows that the noise distribution is almost Gaussian. The half-power width of about 157 mJy/beam is a reliable measure of the rms noise σ in the P map.

5. Conclusions

Our simulated data show that our new method to suppress the polarisation bias almost perfectly reproduces the true polarised intensity even in regions of very low signal-to noise ratios, which was not possible before with any of the previously used methods. In particular, the achievements of our new method are:

  • It estimates the polarisation angle of the source signal in a noisy environment with help of a modified median filter. The corrected polarised intensities P do not suffer from a general residual bias as is the case for all other methods.

  • It works best for smooth variations of the polarisation angle in the source. In the case of sharp jumps in polarisation angle, the modified median filter slightly increases the statistical angle error.

  • It delivers a reliable value for the rms noise. If the noise distributions in the maps of Stokes Û and \hbox{$\hat{Q}$} are Gaussian, the noise distribution in the corrected P map is also Gaussian.

  • The signal-to-noise ratios are measured directly from the P map, without using the maps of Û and \hbox{$\hat{Q}$}.

  • It can also be applied if the noise distributions in the maps of Stokes Û and \hbox{$\hat{Q}$} are different and/or if the distributions are not Gaussian.

  • The maps of corrected polarised intensities P and polarisation angles are reliable even in regions with weak signals and hence allow us to analyse the distribution of polarised intensities and polarisation angles in faint sources. The maps are free of the artefacts produced by the Wardle & Kronberg method (1974).

  • The corrected intensities of the polarised emission P provide reliable integrated flux densities and degrees of polarisation without a cumulative effect of the bias, especially for faint sources.

  • The P map has the same noise distribution as U and Q, which allows us to directly convolve P signals to a larger beamsize and hence to increase the signal-to-noise ratio for diffuse extended emission. However, caution is needed for very small signal-to-noise ratios s, for example when s< 0.5.

  • Features at low intensity levels, like “depolarisation canals”, are smoother in our maps than in those produced by using the previous methods.

Acknowledgments

We thank Axel Jessner and Aritra Basu for useful discussions and Olaf Wucknitz for critical comments on the manuscript.

References

  1. Andernach, H. 1985, Hitch-Hiker’s Guide to NOD2, MPfR, Bonn [Google Scholar]
  2. Beck, R. 2016, Astron. Astrophys. Rev., 24, 4 [NASA ADS] [CrossRef] [Google Scholar]
  3. Beck, R., Berkhuijsen, E. M., & Hoernes, P. 1998, A&AS, 129, 329 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Berkhuijsen, E. M., Beck, R., & Hoernes, P. 2003, A&A, 398, 937 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Fosalba, P., Lazarian, A., Prunet, S., & Tauber, J. A. 2002, ApJ, 546, 762 [NASA ADS] [CrossRef] [Google Scholar]
  6. George, S. J., Stil, J. M., & Keller, B. W. 2012, PASA, 29, 214 [NASA ADS] [CrossRef] [Google Scholar]
  7. Gomes, A. L., Magalhães, A. M., Pereyra, A., & Rodrigues, C. V. 2015, ApJ, 806, 94 [NASA ADS] [CrossRef] [Google Scholar]
  8. Greaves, J. S., Holland, W. S., Jenness, T., & Hawarden, T. G. 2000, Nature, 404, 732 [NASA ADS] [CrossRef] [Google Scholar]
  9. Greisen, E. W. 2003, in Information Handling in Astronomy – Historical Vistas, ed. A. Heck (Dordrecht: Kluwer Academic Publishers), Astrophys. Space Sci. Lib., 285, 109 [Google Scholar]
  10. Haslam, C. G. T. 1974, A&AS, 15, 133 [Google Scholar]
  11. Haverkorn, M., & Heitsch, F. 2004, A&A, 421, 1011 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Hoang, T., & Lazarian, A. 2014, MNRAS, 438, 680 [NASA ADS] [CrossRef] [Google Scholar]
  13. Killeen, N. E. B., Bicknell, G. V., & Ekers, R. D. 1986, ApJ, 302, 306 [NASA ADS] [CrossRef] [Google Scholar]
  14. Laing, R. A., & Bridle, A. H. 2014, MNRAS, 437, 3405 [NASA ADS] [CrossRef] [Google Scholar]
  15. Montier, L., Plaszczynski, S., Levrier, F., et al. 2015, A&A, 574, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Perlman, E. S., Adams, S. C., Cara, M., et al. 2011, ApJ, 743, 119 [NASA ADS] [CrossRef] [Google Scholar]
  17. Pillai, T., Kauffmann, J., Tan, J. C., et al. 2015, ApJ, 799, 74 [NASA ADS] [CrossRef] [Google Scholar]
  18. Planck Collaboration Int. XIX. 2015, A&A, 576, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Scarrott, S. M., Brosch, N., Ward-Thompson, D., & Warren-Smith, R. F. 1986, MNRAS, 223, 505 [NASA ADS] [CrossRef] [Google Scholar]
  20. Shukurov, A., & Berkhuijsen, E. M. 2003, MNRAS, 342, 496 [NASA ADS] [CrossRef] [Google Scholar]
  21. Simmons, J. F. L., & Stewart, B. G. 1985, A&A, 142, 100 [NASA ADS] [Google Scholar]
  22. Tang, Y.-W., Ho, P. T. P., Koch, P. M., et al. 2009, ApJ, 700, 251 [NASA ADS] [CrossRef] [Google Scholar]
  23. Vidal, M., Leahy, J. P., & Dickinson, C. 2016, MNRAS, 461, 698 [NASA ADS] [CrossRef] [Google Scholar]
  24. Vinokur, M. 1965, Ann. Astrophys., 28, 412 [NASA ADS] [Google Scholar]
  25. Wardle, J. F. C., & Kronberg, P. P. 1974, ApJ, 194, 249 [NASA ADS] [CrossRef] [Google Scholar]
  26. Wolleben, M., Landecker, T. L., Reich, W., & Wielebinski, R. 2006, A&A, 448, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Sum of polarised intensities of the simulated images over all pixels.

Table 2

Sum of polarised intensities of the simulated images over each box.

All Figures

thumbnail Fig. 1

Geometrical sketch showing the noise components of the polarisation vectors. The noise vector \hbox{$N_{\hat{P}}$} can be split in two different orthogonal components. Statistically, all NU and NQ contain the measured noise distributions, so the other two projected noise components calculated by the projection of \hbox{$\hat{P}$} onto PT, respectively Pm, contain the same distribution.

In the text
thumbnail Fig. 2

Averaged signal-to-noise ratio (P/σ) as a function of the true ratio PT/σ after applying: the bias suppression method; the Wardle & Kronberg method (using C = 1.2); our new method; and without any bias correction. To show the improvement of the bias suppression with our new method, we created a 100 × 100 pixel Û map and \hbox{$\hat{Q}$} map that contained only Gaussian noise. We added a constant positive value A to the entire Û map and a negative value A to the \hbox{$\hat{Q}$} map. The values vary between 0 < | A | < 3σ. For each A we calculated the corresponding P as the average value over the entire map for all methods given in the plot.

In the text
thumbnail Fig. 3

Polarised intensity reference map. P is calculated from the noise-free U and Q maps using P=U2+Q2\hbox{$P = \sqrt{U^2 + Q^2\ }$}. The lines indicate the orientation of the polarised emission.

In the text
thumbnail Fig. 4

Polarised intensity map calculated with our new method.

In the text
thumbnail Fig. 5

Residual map from our new method. Only Gaussian noise is showing after subtracting the reference map.

In the text
thumbnail Fig. 6

Residual map from the Wardle & Kronberg method. Artefacts are showing at low intensities.

In the text
thumbnail Fig. 7

Histogram of the residual map from our new method (Fig. 5). Intensities are given in arbitrary numbers.

In the text
thumbnail Fig. 8

Histogram of the double-peaked residual map from the Wardle & Kronberg method (Fig. 6). Intensities are given in arbitrary numbers.

In the text
thumbnail Fig. 9

Total intensity of M 31 at 4.85 GHz (in mJy/beam), observed with the Effelsberg 100-m telescope at an angular resolution of 3′ and 1′ pixel size. Bright background sources have been subtracted (from Berkhuijsen et al. 2003).

In the text
thumbnail Fig. 10

Polarised intensity of M 31 at 4.85 GHz (in mJy/beam) determined with our method and an angular resolution of 3′.

In the text
thumbnail Fig. 11

Polarised intensity of M 31 at 4.85 GHz (in mJy/beam) determined with the Wardle & Kronberg method, assuming a 0.2 mJy/beam rms noise for U and Q.

In the text
thumbnail Fig. 12

Histograms of polarised intensity of M 31 at 4.85 GHz: our new method (red) and the Wardle & Kronberg method (green).

In the text
thumbnail Fig. 13

Polarised intensity of M 31 at 1.4 GHz (in mJy/beam) determined with our method and an angular resolution of 3′.

In the text
thumbnail Fig. 14

Polarised intensity of M 31 at 1.4 GHz (in mJy/beam) determined using the Wardle & Kronberg method and assuming an rms noise of 0.2 mJy/beam.

In the text
thumbnail Fig. 15

Difference map of our method (Fig. 13) and the Wardle & Kronberg method (Fig. 14).

In the text
thumbnail Fig. 16

Histogram of the probability density of the noise distribution taken from Eq. (17).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.