Free Access
Issue
A&A
Volume 494, Number 2, February I 2009
Page(s) 611 - 622
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361:20078912
Published online 27 October 2008

WSRT Faraday tomography of the Galactic ISM at ${\lambda \sim 0.86}$ m

I. The GEMINI data set at (l, b) = (181$^\circ $, 20$^\circ $)

D. H. F. M. Schnitzeler1,2 - P. Katgert1 - A. G. de Bruyn3,4


1 - Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
2 - Australia Telescope National Facility, CSIRO, Marsfield, NSW 2122, Australia
3 - ASTRON, PO Box 2, 7990 AA Dwingeloo, The Netherlands
4 - Kapteyn Institute, PO Box 800, 9700 AV Groningen, The Netherlands

Received 24 October 2007 / Accepted 3 October 2008

Abstract
Aims. We investigate the properties of the Galactic ISM by applying Faraday tomography to a radio polarization data set in the direction of the Galactic anti-centre.
Methods. We address the problem of missing large-scale structure in our data, and show that this does not play an important role for the results we present.
Results. The main peak of the Faraday depth spectra in our data set is not measurably resolved for about 8% of the lines of sight. An unresolved peak indicates a separation between the regions with Faraday rotation and synchrotron emission. However, cosmic rays pervade the ISM, and synchrotron emission would therefore also be produced where there is Faraday rotation. We suggest that the orientation of the magnetic field can separate the two effects. By modelling the thermal electron contribution to the Faraday depth, we map the strength of the magnetic field component along the line of sight. Polarized point sources in our data set have rotation measures that are comparable to the Faraday depths of the diffuse emission in our data. Our Faraday depth maps show narrow canals of low polarized intensity. We conclude that depolarization over the telescope beam produces at least some of these canals. Finally, we investigate the properties of one conspicuous region in this data set and argue that it is created by a decrease in line-of-sight depolarization compared to its surroundings.

Key words: magnetic fields - radio continuum: ISM - ISM: magnetic fields - techniques: polarimetric - polarization

1 Introduction

Faraday rotation, the rotation of the plane of linear polarization due to the birefringence of a magneto-ionic medium, provides a powerful tool for exploring the magnetic universe, complementing effects such as Zeeman splitting and the alignment of ellipsoidal dust grains perpendicular to magnetic field lines. From external galaxies we have learnt about the global structure of the magnetic field, about its alignment with the galactic spiral arms, and of the importance of the magnetic pressure in ISM dynamics, which is comparable to the thermal pressure in the cold and warm neutral phases of the ISM (see Beck 2007a,b).

In the Milky Way, rotation measures of pulsars, combined with dispersion measures, yield a rather complicated picture of the large-scale Galactic magnetic field, with evidence for several field reversals (see e.g. Han et al. 2006). Rotation measures of extragalactic sources from the International Galactic Plane Survey (IGPS) indicate a simpler structure of the large-scale field, without much evidence for many field reversals (Brown et al. 2007).

Studies of Faraday rotation of the diffuse Galactic ISM offer a big advantage over other investigations based on pulsars and extragalactic sources, which is that the diffuse emission is visible in all directions. This has been spectacularly demonstrated by the 1.41 GHz survey carried out with the DRAO 26 m telescope by Wolleben et al. (2006). A milestone in this type of work is the polarization survey by Brouw & Spoelstra in the late 1970s (1976; Spoelstra 1984), which has the largest sky coverage in rotation measures of the diffuse emission to date. Since the 1990s interferometric studies at low frequencies have led to new insights in e.g. the properties of the ISM on small scales, and on turbulence in the Galactic ISM (Wieringa et al. 1993; Haverkorn et al. 2004, 2006).

Recently, a method called Faraday tomography has been developed with which one can study the relative distribution of synchrotron-emitting and Faraday-rotating regions along the line of sight (see Brentjens & De Bruyn 2005). Application of this technique for studying the Galactic ISM is a very new field of research, and not much literature is available on this subject. In Schnitzeler et al. (2007b) we discussed a low-frequency data set that covers $\approx$ $50^{\Box}$ in the direction of (l,b) = (181$^\circ $, 20$^\circ $), that we obtained with the WSRT. We concluded that many lines of sight contain at least one Faraday screen, a region in which Faraday rotation and synchrotron emission are not mixed. This result appears counterintuitive, since one would expect synchrotron-emitting cosmic rays to be present everywhere there are magnetic fields. In this paper we present the full analysis of this data set.

Faraday tomography promises to become an interesting and exciting way to investigate the Galactic ISM that will become more and more available in the near future, with the advent of new surveys and radio telescopes that combine high spatial resolution with high frequency resolution, like the GALFACTS survey that is carried out with the Arecibo telescope (Taylor 2004), the LOFAR array (Röttgering et al. 2006), and the MWA SKA-precursor (Bowman et al. 2006). A better understanding of the Galactic foreground will also be useful for studies of the cosmic microwave background.

This paper is the first in a series of articles in which we discuss different regions in the second Galactic quadrant. Here we discuss some of the observational features that we encounter, and we extend some of the techniques that have been used to study the diffuse Galactic ISM with a relatively small number of frequency channels to the many-channel regime. First, we give a short overview of Faraday tomography in Sect. 2, and we illustrate its potential with a couple of examples. In Sect. 3 we discuss the observational characteristics of the GEMINI data set, and in Sect. 4 we present the GEMINI data. We show that Faraday modulation in the foreground can convert polarized emission on large angular scales to smaller angular scales that can be picked up by an interferometer that has no short baselines. We derive the magnetic field strengths implied by the observed Faraday depths of the diffuse emission in Sect. 5, and in Sect. 6 we compare the Faraday depths we derive for the GEMINI data set with the rotation measures that we find for polarized point sources. One conspicuous feature in the GEMINI data are the dark canals in the maps at constant Faraday depth, and we investigate their origin in Sect. 7. In Sect. 8 we present an explanation for a large and bright polarized region in our data set.

2 Notes on Faraday tomography

2.1 Faraday tomography

The wide coverage in $\lambda ^2$ space of our data (from 0.6 m $^2\la \lambda^2 \la 0.9$ m2), combined with the relatively small channelwidth ( $\delta\lambda^2 \approx 1.5$ $\times $ 10-3 m2), allows us to do Faraday tomography, also known as Rotation Measure Synthesis (see e.g. Brentjens & De Bruyn 2005). In this section we introduce the concept of Faraday tomography, and we illustrate how it can be used to study the properties and the relative distribution of regions with Faraday rotation and synchrotron emission along the line of sight.

The observed Stokes parameters $Q(\lambda)$ and $U(\lambda)$ are the 2 orthogonal components of the polarization vector $\vec{P}(\lambda)$ = $Q(\lambda)+{\rm i}U(\lambda)$. $\vec{P}(\lambda)$ is the vector sum of all polarization vectors that are emitted along the line of sight and that are Faraday rotated in the ISM between the point of emission and the observer:

\begin{displaymath}%
\vec{P}(\lambda)\ =\ \int\limits_{{0}}^{\infty} \vec{P}({x}...
...{R})\ {\rm e}^{2{\rm i}\mathcal{R}\lambda^2}{\rm d}\mathcal{R}
\end{displaymath} (1)

where the first integral is over physical distance ``x'', and the Faraday depth of point ``x''

\begin{displaymath}%
\mathcal{R}(x)\ [{\rm rad/m}^2] = 0.81\int_{{\rm source\ at...
...}^{-3}]\ \vec{B}\ [\mu{\rm G}]\cdot {\rm d}\vec{l}\ [{\rm pc}]
\end{displaymath} (2)

measures the total amount of Faraday rotation between the point of emission ``x'' and the observer. In the second integral of Eq. (1) we have replaced the integral over physical distance by an integral over Faraday depth. If there is a one-to-one correspondence between ``x'' and  $\mathcal{R}$, this change in coordinate is trivial. However, a change in the direction of the magnetic field component along the line of sight $B_{\Vert}\ =\ \vec{B}\cdot {\rm d}\vec{l}$ assigns the same $\mathcal{R}$ to different ``x''. To go back from $\mathcal{R}$ to ``x'' thus requires assumptions on the magnetic field geometry. In a future article we will return to this issue. Faraday tomography is based on the inversion of the second integral in Eq. (1):


\begin{displaymath}%
\vec{P}(\mathcal{R}) = K\int\limits_{-\infty}^{\infty} W(\l...
...a^2)\ {\rm e}^{-2{\rm i}\mathcal{R}\lambda^2} {\rm d}\lambda^2
\end{displaymath} (3)

where

\begin{displaymath}K = \left(\int\limits_{-\infty}^{\infty} W(\lambda^2)\ {\rm d}\lambda^2\right)^{-1}
\end{displaymath} (4)

produces the correct normalization for $\vec{P}(\mathcal{R})$. $W(\lambda^2)$ = 1 if the $\lambda ^2$ has been observed, otherwise $W(\lambda^2)$ = 0. This inversion can be interpreted either as a ``coherent addition'' of the observed  $\vec{P}(\lambda)$ by using an assumed Faraday depth to derotate the polarization vectors, or, similarly, as a Fourier transform from $\lambda ^2$ space to $\mathcal{R}$ space. Note that the integral is also calculated over negative values of $\lambda ^2$, which is physically impossible. If the measured  $P(\lambda ^2)$ are in K (or mJy/beam), then also the $P(\mathcal{R})$ will be in K (or mJy/beam), or, more precisely, in K/RMSF width (or mJy/beam/RMSF width). If we want to compare different data sets, with different RMSFs, it is necessary to convert this quantity to K/rad/m2 instead of K/RMSF (equivalently for mJy/rad/m2). Here we only consider one data set, and we did not convert the measured  $P(\mathcal{R})$ to K/rad/m2.

The complex polarization vector $\vec{P}(\lambda)$ can be written in terms of the observed polarized intensity  $P(\lambda)$ and polarization angle  $\Phi(\lambda)$ as $\vec{P}(\lambda) = P(\lambda)\
{\rm e}^{2{\rm i}\Phi(\lambda)}$. Similarly, $\vec{P}(\mathcal{R}) = P(\mathcal{R})\
{\rm e}^{2{\rm i}\Phi(\mathcal{R})}$, where $P(\mathcal{R})$ is the intensity of the polarized emission at Faraday depth  $\mathcal{R}$, and $\Phi(\mathcal{R})$ is the orientation of the electric field vector of the synchrotron radiation emitted at Faraday depth  $\mathcal{R}$, which is perpendicular to the local direction of the magnetic field at that Faraday depth. A Faraday depth spectrum (or $\mathcal{R}$ spectrum) can be constructed by calculating $P(\mathcal{R})$ and $\Phi(\mathcal{R})$ for many $\mathcal{R}$ values.

The Fourier-transform nature of Eq. (3) implies that it shares some characteristics with other Fourier-transform based methods like radio interferometry. Equations (61) to (63) from Brentjens & De Bruyn (2005) describe this behaviour quantitatively, and for convenience we reproduce these equations here. We assume here that the observations uniformly span the $\lambda ^2$ range from $\lambda^2_{{\rm min}}$ to $\lambda^2_{{\rm max}} = \lambda^2_{{\rm min}}+\Delta\lambda^2$ with weight 1, and with weight 0 outside this range.

The finite extent in the $\Delta\lambda^2$ coverage of the data introduces an instrumental response along the $\mathcal{R}$ axis, known as the RMSF (Rotation Measure Spread Function), with which the $P(\mathcal{R})$ spectrum is convolved. The shape of the (normalised) RMSF is given by

\begin{displaymath}%
{\rm RMSF}\ (\mathcal{R})\ =\ K \int\limits_{-\infty}^{\inf...
...^2)\ {\rm e}^{-2{\rm i}\mathcal{R}\lambda^2} {\rm d}\lambda^2.
\end{displaymath} (5)

The FWHM of the RMSF depends on the total range in $\lambda ^2$ space covered by the observations, $\Delta\lambda^2$, according to

\begin{displaymath}%
{\rm FWHM}\ [{\rm radians/m^2}]\ =\ \frac{3.8}{\Delta\lambda^2}\cdot
\end{displaymath} (6)

We replaced the factor of $2\sqrt{3}$ from the original paper by Brentjens & De Bruyn by 3.8, because 3.8 is closer to the FWHM of the sinc RMSF response than $2\sqrt{3} \approx 3.5$.

In general the observations will not cover all wavelengths down to 0 m2, which means that $\mathcal{R}$-extended structures will be missing from the $P(\mathcal{R})$ spectra. This is similar to the missing large-scale structure problem interferometers suffer from when observing extended emission. If the source has a Gaussian $\vec{P}(\mathcal{R})$ distribution along the line of sight, Brentjens & De Bruyn calculate that in the reconstructed $\vec{P}(\mathcal{R})$ spectrum $\mathcal{R}$ scales larger than

\begin{displaymath}%
{\rm maximum\ scale\ [radians/m^2]} \approx \frac{\pi}{\lambda^2_{{\rm min}}}
\end{displaymath} (7)

are suppressed by more than 50% by the missing short wavelengths.

Analogous to how in an interferometer the size of the individual dishes sets the field of view, the channelwidth expressed in $\lambda ^2$, $\delta\lambda^2$, determines the $\mathcal{R}$ for which the sensitivity has dropped to 50% due to smearing of the polarization angles over individual channels, $\mathcal{R}_{{\rm max}}$:

\begin{displaymath}%
\vert\mathcal{R}_{{\rm max}}\vert\ {\rm [radians/m^2]} = \frac{1.9}{\delta\lambda^2}\cdot
\end{displaymath} (8)

Also in this equation we replaced the $\sqrt{3}$ from the paper by Brentjens & De Bruyn by 1.9 because this is closer to the correct value for the HWHM of the sinc RMSF response. Brentjens & De Bruyn discuss the (minor) effect of using (possibly non-contiguous) channels with constant width in frequency instead of $\lambda ^2$.

The addition of polarization vectors originating at different distances along the line of sight in general leads to depolarization, either because the vectors are emitted with different orientations and/or because different amounts of Faraday rotation induce misalignment of the polarization vectors. Faraday tomography can separate the emission coming from subregions with different  $\mathcal{R}$ along the line of sight, thereby reducing the influence of these sources of depolarization.

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig1.ps}
\end{figure} Figure 1:

$P(\lambda ^2)$ and $\Phi (\lambda ^2)$ distributions ( top 2 panels), and $P(\mathcal{R})$ and $\Phi(\mathcal{R})$ spectra ( bottom 2 panels) for a line of sight with a barely resolved main peak that is similar to a Faraday screen. In the text we describe how we derotated the polarization angles in the bottom panel of this figure. The $\chi ^2_{{\rm red}}$ of the RM fit is 0.99 (determined by using the RM with the maximum P(RM)). The $\mathcal{R}$ spectra are sampled from -100 rad/m2 to +100 rad/m2 in steps of 1 rad/m2. In the bottom 2 panels, solid lines indicate the calculated $P(\mathcal{R})$ and $\Phi(\mathcal{R})$ spectra, and the dashed lines indicate the fitted RMSF.

Open with DEXTER

Brentjens & De Bruyn (2005) note that one can also derotate the polarization angles of $\vec{P}(\mathcal{R})$ to $\lambda^2 \ne 0$. In that case the $\lambda ^2$ in the complex exponentials of Eqs. (3) and (5) are replaced by ( $\lambda^2-\lambda_0^2$). Brentjens & De Bruyn showed that derotating to the weighted average of the observed $\lambda ^2$,

\begin{displaymath}%
\lambda_0^2\ =\ \frac{\int\limits_{-\infty}^{\infty}W(\lamb...
...2}{\int\limits_{-\infty}^{\infty}W(\lambda^2){\rm d}\lambda^2}
\end{displaymath} (9)

minimises the polarization angle variation over the main peak of the RMSF. However, now the polarization angle at the peak of the RMSF at Faraday depth  $\mathcal{R}_{{\rm peak}}$ is no longer perpendicular to the magnetic field component in the plane of the sky. To determine the orientation of this magnetic field component, one needs to correct the polarization angles $\Phi $ by

\begin{displaymath}%
\Phi_0(\mathcal{R})\ =\ \Phi(\mathcal{R}) - \mathcal{R}_{{\rm peak}}\lambda_0^2.
\end{displaymath} (10)

We derotated the polarization angles in the bottom panels of Figs. 1 and 2 in this way.

2.2 Examples

We illustrate for a number of geometries of synchrotron-emitting and Faraday-rotating regions the resulting  $P(\lambda)$ and $\Phi(\lambda)$, as well as the $\vec{P}(\mathcal{R})$ that is reconstructed from the observations.

First, consider the case where there is only one infinitely thin source of emission along the line of sight, at Faraday depth  $\mathcal{R_{{\rm 1}}}$: $\vec{P}(\mathcal{R})$ = $P_{{\rm 1}}\ \delta(\mathcal{R}-\mathcal{R_{{\rm 1}}})$, where $\delta(x)$ is the Dirac delta function. Then $\vec{P}(\lambda)$ = $P_{{\rm 1}}\ {\rm e}^{2{\rm i}\mathcal{R_{{\rm 1}}}\lambda^2}$: the $P(\lambda)$ are identical for all observing wavelengths $\lambda$, where we assumed that the synchrotron emission has a spectral index of 0. As Brentjens & De Bruyn (2005) have discussed, a non-zero spectral index introduces only a distortion in the RMSF away from the main peak. $\Phi $ depends linearly on $\lambda ^2$, and since the rotation measure RM $\equiv\ \partial\ \Phi /\partial\ \lambda^2$, RM = $\mathcal{R}$ in this case. The $P(\mathcal{R})$ spectrum that is reconstructed from the observed  $P(\lambda)$ and $\Phi(\lambda)$ is identical to the original  $P(\mathcal{R})$ spectrum, convolved with the RMSF. The situation described in this paragraph occurs in our data: in Fig. 1 we show the measured  $\vec{P}(\lambda^2)$ and reconstructed $\vec{P}(\mathcal{R})$ spectra for a line of sight with an essentially unresolved $P(\mathcal{R})$ spectrum.

Next, consider the case where there are two peaks of height $P_{{\rm 1}}$ and $P_{{\rm 2}}$ at Faraday depths $\mathcal{R_{{\rm 1}}}$ and $\mathcal{R}_{\ {\rm 2}}$. In this case $\vec{P}(\lambda)$ = $P_{{\rm 1}}{\rm e}^{2{\rm i}\mathcal{R}_{\ {\rm 1}}\lambda^2} + P_{{\rm 2}}{\rm e}^{2{\rm i}\mathcal{R}_{\ {\rm 2}}\lambda^2}$. This configuration produces a beat in $P(\lambda)$, because the 2 polarization vectors of length  $P_{{\rm 1}}$ and $P_{{\rm 2}}$ rotate at different speeds $\mathcal{R_{{\rm 1}}}$ and $\mathcal{R}_{\ {\rm 2}}$ in the Stokes (QU) plane. Note that $\Phi $ no longer depends linearly on $\lambda ^2$, which means that RM will not be constant. The $P(\mathcal{R})$ spectrum that is reconstructed from the observed $\vec{P}(\lambda)$ will show the sum of 2 $\delta$-functions convolved with the RMSF. The peaks in the reconstructed $P(\mathcal{R})$ spectrum will not lie at exactly $\mathcal{R_{{\rm 1}}}$ and $\mathcal{R}_{\ {\rm 2}}$, because the presence of the other peak will influence the shape of both peaks. Therefore, since $\Phi(\mathcal{R})$ is derived for the Faraday depth of the peak in the reconstructed $P(\mathcal{R})$ spectrum, the derived  $\Phi(\mathcal{R})$ will not be the intrinsic position angle of the electric field of the emitted synchrotron radiation. When a polarized extragalactic point source is producing one of the 2 peaks, and the diffuse Galactic emission produces the second peak, the $\chi ^2_{{\rm red}}$ of the linear fit of $\Phi $ versus $\lambda ^2$ will be very high. In Fig. 2 we illustrate this case. The peak at $\mathcal{R}$ = -28 rad/m2 is produced by an extragalactic source, and the region around $\mathcal{R}$ = +5 rad/m2 is produced by Galactic emission in the foreground, as we established by comparing this line of sight to adjacent lines of sight, where the extragalactic source is not present.

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig2.ps}
\end{figure} Figure 2:

Identical to Fig. 1, but for a line of sight towards a polarized extragalactic source that also contains diffuse Galactic emission. The  $\chi ^2_{{\rm red}}$ of the RM fit is 8.72.

Open with DEXTER

The previous two examples dealt with infinitely thin emission regions. If the emission region has a finite depth, but is not mixed with Faraday-rotating ISM, the results we discussed in the previous sections still hold. However, if the emitting and Faraday-rotating regions are mixed, the $P(\mathcal{R})$ distribution along the line of sight will no longer be a $\delta$-function, and the reconstructed $P(\mathcal{R})$ spectrum will be wider than the width of the RMSF. If the synchrotron-emitting and Faraday-rotating regions are fully mixed, and if the emissivity and the amount of Faraday rotation per parsec are independent of the position along the line of sight, then

\begin{displaymath}%
\vec{P}(\lambda)\ =\ P_{{\rm i}}\ \frac{\sin(\mathcal{R}_{{...
...\rm e}^{2{\rm i}(\frac{1}{2}\mathcal{R}_{{\rm max}})\lambda^2}
\end{displaymath} (11)

where $\mathcal{R}_{{\rm max}}$ is the Faraday depth of the far side of the slab, and assuming that all radiation along the line of sight is emitted with a position angle of 0$^\circ $. This result was first obtained by Burn (1966), and the distribution of synchrotron-emitting and Faraday-rotating regions along the line of sight that produces it has become known as a ``Burn slab''. Note that $\Phi $ again depends linearly on $\lambda ^2$, but that RM = 0.5 $\mathcal{R}_{{\rm max}}$. Sokoloff et al. (1998) have shown that a linear dependence of $\Phi $ on $\lambda ^2$ occurs not only for a Burn slab, but for every distribution of synchrotron-emitting and Faraday-rotating regions that is symmetric along the line of sight between the observer and $\mathcal{R}_{{\rm max}}$. The reconstructed $P(\mathcal{R})$ spectrum of a Burn slab has a finite extent, and may even be wider than the width of a single RMSF.

It is important to realise that for a general distribution of emitting and Faraday-rotating regions, only those structures show up in the $P(\mathcal{R})$ spectrum that illuminate a column of Faraday-rotating ISM from the back. If Faraday rotation occurs also at larger distances, but is not illuminated from the back by synchrotron emission, this will not appear in the $\vec{P}(\mathcal{R})$ spectrum.

3 The data

Table 1:   Characteristics of the GEMINI data set. Observing dates and times are given for each of the 12 h observing runs, which have been indicated by their shortest baseline length. We calculated the conversion factor between mJy/beam and K at 345 MHz, the average of the $\lambda ^2$ sampling of the (usable) frequency channels in our data set.

The present data set was obtained with the WSRT, a 14-element E-W interferometer of which 4 elements are moveable to improve (u,v) coverage. Each of the telescope dishes is 25 m in diameter. The central coordinates of the region we study (in the constellation Gemini) are $\alpha $ = $\rm 7^h 18^m$ and $\delta$ = 36$^\circ $24 $\hbox{$^\prime$ }$ (J2000.0), which is $l \approx 181^{\circ}$ and $b \approx 20^{\circ}$ in Galactic coordinates. The GEMINI region was observed in 6 12 h observing runs in December 2002 and January 2003 (see Table 1). This yielded visibilities at baselines from 36 to 2760 m, with an increment of 12 m. We tapered the individual frequency channel maps in such a way that the synthesized beamsize for all maps is that of the 385 MHz beam of 2.76$^\prime $ $\times $ 4.70$^\prime $ (RA $\times $ DEC). Combining the 6 12 h observing runs puts the first grating ring at 4.1$^\circ $ (at 350 MHz) from the pointing centre, outside the 3$^\circ $ $\times $ 3$^\circ $ area that we mapped for each individual pointing.

We mapped an area of about 9$^\circ $ $\times $ 9$^\circ $ with a 7 $\times $ 7 pointing mosaic. In each night the same field was observed about 16 times, resulting in visibilities on 16 ``spokes'' in the (u,v) plane, each time integrating for 40 s before moving to the next field. The 1.2$^\circ $ distance between pointing centres suppresses off-axis instrumental polarization to less than 1% (Wieringa et al. 1993). In our analysis we leave out the edges of the mosaic where instrumental polarization effects are not suppressed by mosaicking.

The observations cover the frequency range between 324 and 387 MHz. 202 independent spectral channels are usable from a maximum of 224, each channel being 0.4 MHz wide (where we used a Hamming taper). 10 channels were flagged because Stokes V was contaminated by radio frequency interference, and one channel was manually flagged in all 6 12 h runs.

The data were reduced using the NEWSTAR data reduction package. Dipole gains and phases and leakage corrections were determined using the unpolarized calibrators 3C48, 3C147 and 3C295. The flux scales of both unpolarized and polarized calibrators are set by the calibrated flux of 3C286 (26.93 Jy at 325 MHz - Baars et al. 1977). Due to an a priori unknown phase offset between the horizontal and vertical dipoles, signal can leak from Stokes U into Stokes V. We corrected for this by rotating the polarization vector in the Stokes (UV) plane back to the U axis, assuming that there is no signal in V. The polarized calibrator sources 3C345 and DA240 defined the sense of derotation (i.e. to the positive or negative U-axis). Special care was taken to avoid automatic flagging of real signal on the shortest baselines. From Stokes V, which we assume to be empty, we estimate that the average noise level in the mosaics of the individual channels is 6.2 mJy (2.0 K). Instrumental polarization increases towards the edges of the maps, therefore we excluded these regions in determining the noise level.

These observations were carried out in the evening and at night to limit solar interference and to reduce the importance of ionospheric RM variations. From 2 lines of sight with a strong polarized signal we estimate that the amounts of ionospheric Faraday rotation in the 6 nights are identical to within $\sim$10$^\circ $, so we did not correct for this. Ionosphere models indicate that the ionospheric contribution to RM during the observing nights was only 0.6 rad/m2(Johnston-Hollitt, private communication).

An interferometer will not cover all baseline lengths down to 0 m, which means that maps of the sky that were made using an interferometer will miss structure on large angular scales. We will return to this point in Sect. 4.2.

4 Analysis

4.1 The $\mathcal{R}$ datacube

We calculated $\vec{P}(\mathcal{R})$ maps of the sky for Faraday depths from -1000 rad/m2 to +998 rad/m2 in steps of 6 rad/m2. As the width of the RMSF along the $\mathcal{R}$ axis of the datacube is about 12 radians/m2 for our data set, this gives Nyquist sampling in Faraday depth. Due to the finite channelwidth of the data, the sensitivity of the $P(\mathcal{R})$ spectra will drop at large  $\mathcal{R}$, as described by Eq. (8). This becomes important for $\vert\mathcal{R}\vert$ $\ga$  $1.9/{\rm\delta\lambda^2}$ $\approx$ 1250 rad/m2.

 \begin{figure}
\par\includegraphics[width=17cm,clip]{8912fig3.ps}
\end{figure} Figure 3:

Slices through the $P(\mathcal{R})$ cube that show the strongest diffuse emission. The $\mathcal{R}$ of each slice is indicated in rad/m2 in the upper right corner of each image. The RMSF has a FWHM of 12 rad/m2, and we sampled each RMSF with 2 images. The linear intensity scale is the same for all panels, and it saturates at 6.4 K.

Open with DEXTER

In Fig. 3 we show slices through our $\mathcal{R}$ datacube that show strong Galactic emission. All of the images saturate at 6.4 K. In Figs. 4 and 5 we summarise the information in the $\mathcal{R}$ datacube by plotting for each line of sight the maximum  $P(\mathcal{R})$ along that line of sight, and the $\mathcal{R}$ at which this maximum occurs. We only plotted lines of sight where the main peak is more than twice as high as the second highest peak. From these figures it is clear that the majority of the lines of sight satisfies this condition. Because we treated the data in a slightly different way from what we described in Schnitzeler et al. (2007b), there are some minor differences with the corresponding figures in that paper. From the $P(\mathcal{R})$ maps at $\vert\mathcal{R}\vert > 200$ rad/m2 we estimate that the noise level in the $P(\mathcal{R})$ slices is 0.5 mJy (0.14 K).

The complexity of a line of sight depends on whether the main peak is resolved, and on how strong second and higher-order peaks are. The $\Delta $ criterion that we introduced in Schnitzeler et al. (2007b) can be used to address the issue of resolution. We defined $\Delta $ as the root-mean-square vertical separation in the $\mathcal{R}$ spectrum between the main peak in the $P(\mathcal{R})$ spectrum and the best-fitting RMSF, and we compare the observed $\Delta $ to the distribution of $\Delta $ that we simulated for an input signal of the appropriate strength + noise. In Fig. 6 we plot the lines of sight that have an unresolved main peak as green pixels. ``Unresolved'' means that the $\Delta $ value of the main peak is $\le$ the $\Delta $ that we found for 99% of our simulations of signal + noise ($\approx$ the 3$\sigma $ level of the Rayleigh distribution).

Lines of sight with linear $\Phi (\lambda ^2)$ relations have special geometries of Faraday-rotating and synchrotron-emitting regions, as discussed in Sect. 2. We therefore fitted a straight line to our $\Phi (\lambda ^2)$ datapoints, where we included the periodicity of the data by using the procedure described in Schnitzeler et al. (2007b). Lines of sight with $\chi ^2_{{\rm red}}$ < 2 are indicated in Fig. 6 as yellow pixels, and lines of sight that have both a not measurably resolved main peak and a low $\chi ^2_{{\rm red}}$ are indicated as red pixels. Since many lines of sight have a low  $\chi ^2_{{\rm red}}$, but a high $\Delta $, there are many yellow pixels, but relatively few red pixels.

We simulated the $P(\mathcal{R})$ response from the main peak plus noise by calculating the Fourier transform from Eq. (3) for an input signal of a given strength + the noise levels that we derived from the Stokes V maps of the individual frequency channels. We assume that the Stokes V maps contain no signal. We ran our simulations for a grid of input signal strengths, that was matched to the strengths of the main peak in the $P(\mathcal{R})$ spectra that we encountered in the data. In this way we calculated for each $\mathcal{R}$ in the spectrum the level below which 99% of the RMSF plus noise simulations lie.

The main peak in the $P(\mathcal{R})$ spectrum is not measurably resolved for about 8% of the lines of sight in the current data set, in the sense that its value of $\Delta $ is lower than the $\Delta $ that we find for 99% of our simulations of a RMSF + noise. This value is lower than the 14% quoted in Schnitzeler et al. (2007b). This difference might be produced because we now use frequency channel maps that all have been tapered to the same beamsize. A not measurably resolved peak means that a Faraday-rotating region is physically separated from a region with synchrotron emission; mixing of the 2 regions produces broad peaks in the $P(\mathcal{R})$ spectrum. Since we only tested the main peak in our $P(\mathcal{R})$ spectra, and no other features in the $P(\mathcal{R})$ spectra, emission regions with a Faraday screen in front of them will be even more numerous than in 8% of the lines of sight. RFI and calibration errors can only increase $\Delta $, which is another reason why there are probably more Faraday screens present in the ISM.

Intuitively one would not expect such a separation between the emitting and Faraday-rotating regions along the line of sight: cosmic rays that produce synchrotron radiation are expected to be present wherever there is a magnetic field. Therefore one possibility could be that we simply do not see the synchrotron emission because the full magnetic field vector is pointing along the line of sight. The magnetic field component perpendicular to the line of sight, which determines the synchrotron emissivity, is in that case 0. Such a region would from our point of view then only produce Faraday rotation. Even though the strengths of the magnetic field component along the line of sight that we determine in Sect. 5 are low, these are averages over the line of sight, and they are biased towards regions where the thermal (i.e. 104 K) electron density is higher.

The majority of lines of sight in our data have a secondary peak that is higher than what would be expected for one in 104 noise realisations. However, the main peak is typically more than twice as high as the secondary peak, and thus dominates most of the lines of sight in the GEMINI region. This makes these lines of sight relatively simple, and we did not deconvolve the $P(\mathcal{R})$ spectra. The reason why these lines of sight are not very complex could be because we are looking in the direction of the Galactic anti-centre, at an intermediate Galactic latitude.

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig4.ps}
\end{figure} Figure 4:

$P(\mathcal{R})$ of the highest peak in the $\mathcal{R}$ spectrum. Each pixel represents one telescope beam, and the pixel size is 2.8$^\prime $ $\times $ 4.9$^\prime $. The lines of sight that are shown in this plot have a second peak in the $P(\mathcal{R})$ spectrum that is at most half the strength of the main peak. Out of a total of 22.800 lines of sight, 2720 do not obey this selection criterion, and they are shown in black, as are the 769 lines of sight that have $P(\mathcal{R})$ < 2 K. The 1$\sigma $ noise level in the $P(\mathcal{R})$ maps is 0.14 K.

Open with DEXTER

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig5.ps}
\end{figure} Figure 5:

$\mathcal{R}$ corresponding to the $P(\mathcal{R})$ plotted in Fig. 4. The lines of sight plotted in these figures are selected by using the same criterion. Lines of sight indicated in white have $\mathcal{R}$ > 30 rad/m2. Lines of sight indicated in black either did not pass the selection criterion or have $\mathcal{R}$ < -6 rad/m2 (2720 and 727 lines of sight respectively).

Open with DEXTER

 \begin{figure}
\par\includegraphics[width=8cm,clip]{8912fig6.ps}
\end{figure} Figure 6:

Distribution of independent lines of sight with a main peak that satisfies the $\Delta $ criterion (green pixels), or with $\chi ^2_{{\rm red}}$ < 2 and max( $P(\mathcal{R})$) > 2 K (yellow pixels), or that satisfy all three criteria (red pixels). Lines of sight that did not pass these criteria are shown as black pixels. The outer edges of the mosaic show an increase in  $\chi ^2_{{\rm red}}$, which explains the absence of yellow pixels.

Open with DEXTER

4.2 Missing large-scale structure

Structure on large angular scales is missing from interferometric observations if information on short baselines is not available. Differences in Faraday rotation in between the source of the polarized emission and the observer can for neighbouring lines of sight ``transfer'' polarized signal from large angular scales to small enough angular scales so that it can be picked up by the interferometer. Since this effect does not operate on Stokes I, this is the canonical explanation why structure in polarized intensity does not appear to have a counterpart in total intensity. In this section we expand on this idea, and we show that modulation by a linear gradient in Faraday depth  $\mathcal{R}$ shifts the angular frequency spectrum as a whole.

For a linear variation of $\mathcal{R}$ with position l (i.e. a fixed gradient  $\nabla{\mathcal{R}}$), and wavelength $\lambda$, this modulation of the polarization vector $\vec{P}$ can be expressed as

\begin{displaymath}%
\vec{P}_{{\rm obs}}(l)\ =\ \vec{P}_{{\rm em}}(l)\ {\rm e}^{2{\rm i}\nabla{\mathcal{R}}\lambda^2 l}
\end{displaymath} (12)

where subscripts ``em'' and ``obs'' refer to the emission before and after it gets modulated by the foreground gradient in $\mathcal{R}$. This linear modulation with position shifts the angular frequency spectrum as a whole:
           $\displaystyle %
{\rm FT}(\vec{P}_{{\rm obs}})(u)$ = $\displaystyle \int_{-\infty}^{\infty} {\rm d}l\ \vec {P}_{{\rm obs}}(l)\ {\rm e}^{-2\pi {\rm i}ul}$ (13)
  = $\displaystyle \int_{-\infty}^{\infty} {\rm d}l\ \vec {P}_{{\rm em}}(l)\ {\rm e}^{2{\rm i}\nabla{\mathcal{R}}\lambda^2l} {\rm e}^{-2\pi {\rm i}ul}$  
  = $\displaystyle \int_{-\infty}^{\infty} {\rm d}l\ \vec {P}_{{\rm em}}(l)\ {\rm e}^{-2\pi {\rm i}\left(u-\frac{1}{\pi}\nabla{\mathcal{R}}\lambda^2\right)\ l}$  
  = $\displaystyle {\rm FT}(\vec{P}_{{\rm em}})\left(u-\frac{1}{\pi}\nabla{\mathcal{R}}\lambda^2\right).$ (14)

One consequence of this is that foreground Faraday modulation makes the 0-angular frequency from $\vec{P}_{{\rm em}}(l)$ visible if $\vert\nabla\mathcal{R}\vert$ is large enough. To make it detectable at the 25 m (projected) baseline length, which can be obtained with the WSRT, requires a $\vert\nabla\mathcal{R}\vert$ $\approx$ 6 rad/m2/field of view at 350 MHz (where 1 field of view is 3$^\circ $ wide for our observations). Since we observe gradients of this magnitude in the individual pointings of the GEMINI data set (Fig. 5), Figs. 4 and 5 contain at least the 0-angular frequency = total polarized intensity component of  $\vec{P}_{{\rm em}}$.

The direction of the shift in Eq. (14) is set by the sign of $\nabla{\mathcal{R}}$. Since the WSRT only measures one half of the (u,v)-plane, this means that the angular frequency spectrum could move away from the measurement points in the (u,v)-plane. However, $\vec{P}_{{\rm obs}}(l)=Q_{{\rm obs}}(l)+{\rm i}U_{{\rm obs}}(l)$, and since $Q_{{\rm obs}}(l)$ and $U_{{\rm obs}}(l)$ are real quantities, their Fourier transforms are hermitian. The gradient from Eq. (14) always moves both the positive and the negative angular frequencies of $Q_{{\rm em}}(l)$ and $U_{{\rm em}}(l)$ towards or away from the origin in the (u,v)-plane. Therefore, $\vec{P}_{{\rm obs}}(l)$ can be reconstructed even if only one half of the (u,v)-plane is observed.

Note that foreground Faraday modulation shifts some angular frequencies from $\vec{P}_{{\rm em}}(l)$ towards smaller (in an absolute sense) angular frequencies. At some point these angular frequencies from $\vec{P}_{{\rm em}}(l)$ will no longer be visible for an interferometer, because it can only measure visibilities from a certain shortest baseline. For the analysis that we present in the remainder of this article it is sufficient that our data contain the total polarized intensity component of $\vec{P}_{{\rm em}}$. Single-dish observations are therefore still required to produce maps that contain all angular frequencies.

If Faraday rotation in the foreground produces strong variations in $\mathcal{R}$ over the telescope beam, then this will also lead to beam depolarization. For a Gaussian beam, a gradient of more than 12 radians/m2/field of view is required to give more than 10% beam depolarization (Sokoloff et al. 1998), and such steep gradients are not observed in the GEMINI data.

As a final note in this section, we showed in Schnitzeler et al. (2007a) that the polarization angle gradients in the WENSS data in the region 130$^\circ $ < l < 170$^\circ $, -5$^\circ $ < b < 30$^\circ $ are on average about 2 radians/m2/$^\circ $ along Galactic latitude, which translates into the 6 radians/m2/field of view required to shift the angular frequency spectrum towards observable frequencies. These WENSS gradients therefore also show at least the total polarized intensity component of the angular frequency spectrum.

5 The line-of-sight component of the large-scale magnetic field

In this section we derive the strength of the magnetic field component along the line of sight, $B_{\Vert}$. We use the Faraday depths from Fig. 5, and we determine the contribution from thermal electrons to $\mathcal{R}$ (Eq. (2)) by converting H$\alpha $ intensities from the WHAM survey (Haffner et al. 2003) first to emission measures EM = $\int_{0}^{\infty} n^2_{{\rm e}}{\rm d}l$, and then to dispersion measures DM = $\int_{0}^{\infty} n_{{\rm e}}{\rm d}l$, using an empirical relation between EM and DM that Berkhuijsen et al. (2006) derived for Galactic pulsars. The integrals for EM and DM are over the entire line of sight, [EM] = cm-6 pc, [DM] = cm-3 pc, $[n_{{\rm e}}]$ = cm-3 and [dl] = pc.

To estimate DM, we could also have used the electron density models by Reynolds (1991) or by Cordes & Lazio (2003) as alternatives to the relation between EM and DM that Berkhuijsen et al. derived. In the Reynolds model, 40% of the line of sight contains free electrons, at a density of 0.08 cm-3, and there are no free electrons in the remainder of the line of sight. Cordes & Lazio fitted a smooth electron density model consisting of a thin and thick disc + spiral arms to the DM measured for about 1200 Galactic pulsars and 100 extragalactic sources. Small-scale structure in the Galactic ISM is included in their model as a list of local clumps and voids. Berkhuijsen et al. used DM measured for a sample of 157 pulsars, and accurately determined the EM between us and the pulsar by correcting for the EM that is produced beyond the pulsar, and for reddening occurring in front of the pulsar. By combining Berkhuijsen et al.'s statistical description of the ISM with the WHAM H$\alpha $ intensities that are measured on a $\approx$$1^{\circ}$ $\times $ $1^{\circ}$ grid, we can more accurately determine DM in specific directions than is possible with the models by Reynolds or by Cordes & Lazio.

Significant numbers of free electrons are present in both the warm ionized medium (WIM) and hot intercloud medium (HIM) phases of the ISM. Snowden et al. (1997) modelled the Galactic X-ray halo from ROSAT diffuse X-ray background maps. Their model consists of a cylindrical bulge of radius 5.6 kpc, with a HIM electron density that decays exponentially away from the Galactic midplane. The HIM electron density in their model is in the Galactic midplane a factor of 10 smaller than the average WIM electron density (3.5 $\times $ 10 $^{-4}~{\rm cm}^{-3}$ compared to $\sim$$\times $ $10^{-3}~{\rm cm}^{-3}$). If we assume that the HIM electron densities in the outer Galaxy (in the direction of GEMINI) are comparable to the densities in the inner Galaxy, then the DM produced in the HIM is small compared to the DM of the WIM, and we therefore neglected the former.

In the top panel of Fig. 7 we plot the H$\alpha $ intensities (in Rayleigh) for WHAM lines of sight in the GEMINI region. To convert the H$\alpha $ intensities to EM, we use Eq. (1) from Haffner et al. (1998):

\begin{displaymath}%
{\rm EM}=2.75 T_4^{0.9}\ I_{{\rm H}\alpha}\ {\rm e}^{{2.2\ E(B-V)}}.
\end{displaymath} (15)

T4 is the temperature of the WIM gas in units of 104 K, which is typically 0.8 (Reynolds 1985). $I_{{\rm H}\alpha}$ is the H$\alpha $ intensity in Rayleigh, and E(B-V) is the interstellar B-V reddening. We determined EM for each of the WHAM lines of sight by using reddening values from Schlegel et al. (1998) and the programs available on Marc Davis' website[*].

The 157 pulsars with $\vert b\vert> 5^{\circ}$ in the sample by Berkhuijsen et al. follow the relation

\begin{displaymath}%
{\rm EM}= (0.042 \pm 0.014)\ {\rm DM}^{1.47 \pm 0.09}.
\end{displaymath} (16)

With the $\mathcal{R}$ from Fig. 5 and the DM that we calculated from the WHAM H$\alpha $ intensities, we can then calculate the electron-density weighted average $B_{\Vert}$ along the line of sight, $\langle B_{\Vert}\rangle$:

\begin{displaymath}%
\langle B_{\Vert}\rangle \equiv \frac{\int\limits_{{\rm lin...
...ight}} n_{\rm e} {\rm d}l} = \frac{\mathcal{R}}{0.81~{\rm DM}}
\end{displaymath} (17)

where $\langle B_{\Vert}\rangle$ is in $\mu $G. In the bottom panel of Fig. 7 we show the $\langle B_{\Vert}\rangle$ that we calculated in this way. Since the WHAM beam (FWHM $\approx$ 1$^\circ $) is much larger than the synthesized WSRT beam, we used the $P(\mathcal{R})$-weighted $\mathcal{R}$ average over the WHAM beam to calculate $\langle B_{\Vert}\rangle$. Also note that since Schlegel et al. give reddening values for the entire line of sight through the Galaxy, the EM and DM that we calculate are also values for the entire line of sight. Since we do not know where the polarized emission that we observe originates, the $\vert\langle B_{\Vert}\rangle\vert$ that we derive from Eq. (17) are lower limits to the actual  $\vert\langle B_{\Vert}\rangle\vert$.

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig7.ps}
\end{figure} Figure 7:

WHAM H$\alpha $ intensities in Rayleigh ( top panel) and the electron-density weighted strengths of the magnetic field component parallel to the line of sight $\langle B_{\Vert}\rangle$ [$\mu $G] that we derived ( bottom panel). Open circles indicate positive  $\langle B_{\Vert}\rangle$; there are no negative $\langle B_{\Vert}\rangle$ < -0.1 $\mu $G in this figure. $\vert\langle B_{\Vert}\rangle\vert$ < 0.1 $\mu $G are shown as ``+''.

Open with DEXTER

An upper limit for $B_{\Vert}$ can be estimated by using a pitch angle of 8$^\circ $ and a strength of 4 $\mu $G for the large-scale magnetic field (both values from Beck 2007a), which gives $B_{{\rm reg},\Vert,{\rm max}} \approx 0.6$ $\mu $G. In the presence of a random magnetic field, $B_{\Vert}$ approaches the $B_{{\rm reg},\Vert}$ component of the regular field if either the reversals in the direction of the random field are unimportant, or if the line of sight is long enough so that the contributions to $\mathcal{R}$ from such small-scale reversals cancel each other along the line of sight. As discussed in Haverkorn et al. (2004), the random magnetic field has a correlation length of order 10 parsec. Shortly we will estimate that the length of the line of sight through ionized regions in the direction of GEMINI is of the order of 130 parsec. The number of ``draws'' of the random field is therefore small, which means that the random magnetic field will have a non-negligible contribution (compared to the large-scale field) to the $\langle B_{\Vert}\rangle$ that we derive. The $\langle B_{\Vert}\rangle$ in the bottom panel of Fig. 7 already show that the $\langle B_{\Vert}\rangle$ in GEMINI vary much more than what we would expect on the basis of the large-scale magnetic field alone.

The $\langle B_{\Vert}\rangle$ in the bottom panel of Fig. 7 are smaller than or comparable to the $B_{{\rm reg},\Vert,{\rm max}}$ that we estimated. The $\langle B_{\Vert}\rangle$ we derive are small, which can be expected since we are looking in the direction of the Galactic anti-centre, where the regular magnetic field is mostly perpendicular to the line of sight. It is clear from the bottom panel in Fig. 7 that variations in EM alone are not enough to explain the measured variation in $\mathcal{R}$, but that variations in $B_{\Vert}$ are required as well. Further support for this is given by the fact that $\mathcal{R}$ in Fig. 5 changes sign in the area covered by the GEMINI data, which can only be explained by a change in the direction of $B_{\Vert}$. The $B_{\Vert}$ structure in Fig. 5 is not associated with any of the 3 largest non-thermal loops in the Galaxy, the North Polar Spur, Cetus Arc and Loop III. The locations of these structures are described in Elliott (1970).

Berkhuijsen et al. also constructed a model of the electron density distribution along the line of sight, that we can use to estimate the length of the line of sight that produces the observed EM. Similar to the model by Reynolds (1991), their model consists of cells that are empty or filled with a constant electron density. Contrary to the Reynolds model, the electron density of the filled cells (``clouds''), $n_{{\rm e,cloud}}$, is allowed to vary between lines of sight:

\begin{displaymath}%
n_{{\rm e,cloud}} = \frac{{\rm EM}}{{\rm DM}} = 0.12\ {\rm EM}^{\ 0.32}
\end{displaymath} (18)

where we used Eq. (16). Berkhuijsen et al. found that the filling factor f of electron clouds depends inversely on the electron density in the clouds: $f= (0.0184\pm 0.0011)\ n_{{\rm e,cloud}}^{-1.07\pm0.03}$. We can use this relation to determine the length of the line of sight D (in parsec) that produces the observed EM; note that D includes both the fraction of the line of sight that passes through the electron clouds and also the electron-free regions between the clouds:

\begin{displaymath}%
D = \frac{1}{f}\left(\frac{{\rm DM}^2}{{\rm EM}}\right) = 404\ {\rm EM}^{\ 0.70}.
\end{displaymath} (19)

E(B-V) = 0.071 mag for ( $\alpha ,\delta $) = (110$^\circ $, 37$^\circ $), which means that the average H$\alpha $ intensity in GEMINI of 1.78 R translates into an EM = 4.68 cm-6 pc. This gives $n_{{\rm e,cloud}}$ = 0.19 cm-3, f = 0.11, and D = 1.2 kpc. This filling factor and D imply that the combined pathlength through the electron clouds is f $\times $ D $\approx$ 130 parsec. We calculated EM and DM for the entire line of sight, but it is unclear whether we can observe the polarized emission coming from the furthest parts along the line of sight. D is therefore an upper limit to the length of the line of sight that produces the $\mathcal{R}$ in Fig. 5.

6 Polarized point sources

Since the line of sight towards extragalactic sources passes through the entire Milky Way, whereas different depolarization effects can limit the line of sight of the diffuse emission, the information we obtain from the $\mathcal{R}$ of extragalactic sources can complement what we learn from the diffuse emission. However, extragalactic sources will also have an intrinsic  $\mathcal{R}$.

To find polarized point sources, we first made maps that only contain baselines >250 m, to filter out the strong diffuse emission from our $P(\mathcal{R})$ datacube. These maps use a Gaussian taper that reaches a value of 0.25 at a baseline length of 2500 m. We then Fourier transformed these channel maps from $\lambda ^2$ space to $\mathcal{R}$ space, and constructed $P(\mathcal{R})$ maps from $\mathcal{R}$ = -1000 rad/m2 to $\mathcal{R}$ = +998 rad/m2, in a way identical to what we described in Sect. 4.1.

The polarized intensity of an extragalactic source with rotation measure RM will be attenuated by bandwidth depolarization according to sinc(|RM $\vert\delta\lambda^2$), if the polarization angle varies linearly with $\lambda ^2$. Here $\delta\lambda^2$ is the $\lambda ^2$ width of the frequency channels. The channels that we use have on average a $\delta\lambda^2$ $\approx$ 1.5 $\times $ 10-3 m2, therefore bandwidth depolarization of a point source with $\vert\mathcal{R}\vert$ = 1000 rad/m2 will reduce the signal to 66%, which would still be detectable.

We looked for each line of sight in the $P(\mathcal{R})$ datacube if 1) there was a Stokes I counterpart brighter than 150 mJy; 2) if the maximum $P(\mathcal{R})$ along the line of sight was more than a 5$\sigma $ detection; and 3) if the ratio of the maximum $P(\mathcal{R})$ along the line of sight/ total intensity I was larger than 0.02, meaning that this source is more than 2% polarized (the instrumental polarization level is about 1%). We also exclude sources with $\vert\mathcal{R}\vert$ $\le$ 4 rad/m2 because instrumentally polarized point sources will show up in the $P(\mathcal{R})$ datacube at 0 rad/m2, and the finite width of the RMSF along the $P(\mathcal{R})$ axis will produce a strong signal also in the vicinity of 0 rad/m2. Furthermore we excluded a region of about 0.9$^\circ $ from the edge of the mosaic, where there is no overlap between pointings, and instrumental polarization levels are therefore much higher there. We then fitted a RM to the $\Phi (\lambda ^2)$ distribution of the sources that satisfy these criteria in the way described in Schnitzeler et al. (2007b).

In Table 2 we list the properties of the sources that we found, and we plot the RM of these sources in Fig. 8. The low $\chi ^2_{{\rm red}}$ of the RM fits indicate that there is only one $P(\mathcal{R})$ component along the $\mathcal{R}$ axis, and that we have effectively filtered out the diffuse Galactic contribution to the $P(\mathcal{R})$ spectrum. This, in combination with the fact that these sources are detected in maps where we left out low angular frequencies on the sky, means that the sources that we detected are compact in all 3 dimensions of the $P(\mathcal{R})$ cube, and are therefore likely to be either polarized pulsars or polarized extragalactic sources. Since the surface density of the latter is so much higher than that of the former, the sources that we detected are probably extragalactic in origin.

Most of the sources in the upper left quadrant of Fig. 8 show positive RMs of $\approx$ 5-15 rad/m2, which is similar to the Faraday depths of the diffuse emission from Fig. 5. In other regions the source density is insufficient to reach definite conclusions. The RM that we observe for an extragalactic source is the combination of a Galactic RM and a RM that is intrinsic to the source. However, since a number of sources show the same RM, and since the intrinsic RM of these sources are uncorrelated, what we see is probably mostly the Galactic RM contribution for the entire line of sight through the Galaxy.

Table 2:   Properties of the polarized point sources we found in our data. Shown are the equatorial coordinates (J2000.0) of the source in decimal degrees, the RM that we fitted and its error (both in rad/m2), the reduced $\chi ^2$ ( $\chi ^2_{{\rm red}}$) of the fit, the maximum $P(\mathcal{R})$ along the line of sight and the total intensity I (both in mJy), and the polarization fraction (in %).

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig8.ps}
\end{figure} Figure 8:

RM distribution for the polarized point sources from Table 2. The size of the circles is proportional to RM, shown in the scale on top of the figure in units of radians/m2. Black circles indicate negative RM, and white circles indicate positive RM. |RM|  $\protect\la$ 4 rad/m2 are missing from this figure because of our selection criteria. We excluded the grey region when looking for polarized point sources, because instrumental polarization levels are higher there.

Open with DEXTER

7 Depolarization canals in P ${(\mathcal{R})}$ maps

Dark, narrow channels (``canals'') can be clearly seen in the images of Fig. 3. Canals have been known to exist in $P(\lambda)$ images (e.g. Haverkorn et al. 2000), but, as far as we know, canals in $P(\mathcal{R})$ maps have only been mentioned by De Bruyn et al. (2006).

In Fig. 9 we show slices through the $P(\mathcal{R})$ cube that have particularly strong depolarization canals. It appears from this figure that the canals essentially are not changing position with changing $\mathcal{R}$. This is demonstrated by the reference crosses that we positioned on three canals in this figure. In the first and second rows of Fig. 10 we plot the $P(\lambda ^2)$ (left column) and $\Phi (\lambda ^2)$ (right column) spectra for the two canal pixels from the top row of Fig. 9 (``CH 1'' and ``CH 2'' respectively). In the bottom row of Fig. 10 we plot the spectra for the canal pixel in the bottom row of Fig. 9.

In Fig. 10 we also plot the $P(\lambda ^2)$ and $\Phi (\lambda ^2)$ spectra for 2 lines of sight that lie on either side of the central canal pixel, at half a (synthesized) beamwidth away. These adjacent lines of sight share pixels with the line of sight that is centred on the canal pixel, and with these lines of sight we can investigate the nature of the canal pixels.

The $\Phi (\lambda ^2)$ spectra of the neighbouring lines of sight were unwrapped by fitting a straight line to the polarization angles by hand. As a first guess for the slope of this line, we used the $\mathcal{R}$ at which the $P(\mathcal{R})$ of the neighbouring lines of sight is high. Furthermore we required that for each line of sight the difference in polarization angle between the successive (plotted) wavelength2 channels must be smaller than 90$^\circ $.

The difference in $\Phi (\lambda ^2)$ for lines of sight adjacent to the canal pixels is large, often even almost 90$^\circ $, averaged over the available $\lambda ^2$. Also, the polarized brightness temperature of the 2 neighbouring lines of sight of each canal pixel are on average comparable. These 2 facts combined indicates that beam depolarization, i.e. summing polarization vectors with similar lengths but different orientations over the (synthesized) telescope beam, is an important factor in creating canals in the $P(\mathcal{R})$ maps. In Haverkorn et al. (2000) the same mechanism was identified to be producing canals in $P(\lambda ^2)$ maps. Line-of-sight depolarization, occurring for example as the nulls in Eq. (11), produces canals only at certain wavelength2. Shukurov & Berkhuijsen (2003) suggested that this mechanism is producing canals that are seen towards M 31. This type of canal is characterised as extended along the $\mathcal{R}$ axis of the $P(\mathcal{R})$ cube, contrary to the properties of the canal pixels we present here.

 \begin{figure}
\par {\hspace*{5mm}\includegraphics[width=15.5cm,clip]{8912fig9.p...
...mm}}\vspace*{3mm}
\includegraphics[width=16.5cm,clip]{8912fig10.ps}
\end{figure} Figure 9:

Depolarization canals in the $P(\mathcal{R})$ datacube. The $\mathcal{R}$ [rad/m2] of each of the $P(\mathcal{R})$ slices is indicated in the lower left. All maps have the same linear greyscale that saturates at 4.8 K. The area shown in the top row of this figure lies in the upper left of Fig. 3. For the 2 pixels that are indicated by crosses, and that lie in depolarization canals, we plot the $P(\lambda ^2)$ and $\Phi (\lambda ^2)$ spectra in the first row of Fig. 10 (left cross, ``CH 1'' in Fig. 10) and middle row of Fig. 10 (right cross, ``CH 2''). The bottom row shows a canal in a different part of the mosaic, and the $P(\lambda ^2)$ and $\Phi (\lambda ^2)$ spectra of the canal pixel that is indicated by a cross are shown in the bottom row of Fig. 10 (``CH 3'').

Open with DEXTER
 \begin{figure}
\includegraphics[width=16.3cm,clip]{8912fig11.ps}
\end{figure} Figure 10:

Polarized brightness temperature  $T_{{\rm b,pol}}$ [K] and polarization angle $\Phi $ [degrees] spectra for 3 canal pixels, and two adjacent pixels half a beamwidth from each canal pixel. Top row: canal pixel at ( $\alpha ,\delta $) = (111.24$^\circ $, 39.72$^\circ $), and two adjacent lines of sight at $\alpha $ = 111.27$^\circ $ (dotted line), and $\alpha $ = 111.21$^\circ $ (dashed line). Middle row: canal pixel at ( $\alpha ,\delta $) = (110.98$^\circ $, 39.57$^\circ $), and adjacent lines of sight at $\alpha $ = 111.01$^\circ $ dotted line, and $\alpha $ = 110.95$^\circ $, dashed line. Bottom row: canal pixel at ( $\alpha ,\delta $) = (111.50$^\circ $, 34.94$^\circ $), and two adjacent lines of sight at $\alpha $ = 111.53$^\circ $, dotted line, and $\alpha $ = 111.47$^\circ $, dashed line. The 1$\sigma $ errors for polarized intensity and polarization angle are 2 K and 15$^\circ $ (for a S/N = 2 detection) respectively Every 8 ${{\rm th}}$ channel is plotted. Bad frequency channels were selected using the same selection criteria as in Sect. 3, and discarded.

Open with DEXTER

8 The area around (l, b) = (109$^\circ $, 34.5$^\circ $)

The 2$^\circ $ $\times $ 2$^\circ $ area centred on (l,b) = (109.5$^\circ $, 34.5$^\circ $) is bright in polarized intensity (Fig. 4), shows a very uniform distribution in $\mathcal{R}$ (Fig. 5), and has a relatively small surface density of strong depolarization canals. Furthermore the $\chi ^2_{{\rm red}}$ fitted to the $\Phi (\lambda ^2)$ spectra are low (Fig. 6): 30% of the lines of sight in this area have a $\chi ^2_{{\rm red}}$ < 2, and nearly 75% have a $\chi ^2_{{\rm red}}$ < 3. There is no enhanced H$\alpha $ emission in this region (see the top panel of Fig. 7), and $B_{\Vert}$ varies only little (bottom panel of Fig. 7). The single-dish 408 MHz total intensity map by Haslam et al. (1982) shows hardly an increase in brightness temperature in this region. There are not enough lines of sight in the Brouw & Spoelstra (1976) data set to discern the area around (l,b) = (109.5$^\circ $, 34.5$^\circ $) from its surroundings.

The increased polarized brightness temperature with respect to its surroundings could be due to enhanced synchrotron emission, or a local decrease in depolarization, or to a combination of these. The polarized brightness temperature of this region is more than 3 K larger than that of its surroundings, and this would translate into a brightness temperature difference of at least 3 K/0.7 = 4.3 K, assuming that the radiation is 70% polarized. Lower polarization percentages increase this brightness temperature contrast even further. An enhancement in synchrotron emission can be ruled out as an explanation for the increased  $T_{{\rm b,pol}}$ of this area, since this is not seen in the Haslam data, and since beam averaging is not strong enough to attenuate the total intensity signal (the Haslam beam measures 0.85$^\circ $ $\times $ 0.85$^\circ $).

The small variation in $\mathcal{R}$ over this region, combined with an underabundance of strong depolarization canals, furthermore indicates that depolarization across the synthesized telescope beam is not a strong effect. This leaves a decrease in the amount of depolarization along the line of sight as an explanation for the observed increase in $T_{{\rm b,pol}}$. The low  $\chi ^2_{{\rm red}}$ of lines of sight in this area, in combination with secondary peaks that are not stronger than half the strength of the main peak, also indicates that the lines of sight in this area are not very complex, and line-of-sight depolarization effects are therefore not very important. We thus think that the origin of this particular region lies in a decreased amount of depolarization compared to its surroundings.

9 Conclusions

We applied Faraday tomography to high spectral-resolution radio polarization data to study the properties of the magnetised Galactic ISM. We showed that differential Faraday modulation in the foreground can shift the emitted angular frequency spectrum as a whole. In particular the 0-angular frequency from the emitted radiation is modulated in this way towards smaller angular scales, that are observable with an interferometer. We quantified how strong a linear gradient in $\mathcal{R}$ should be to accomplish this, and we showed that the gradients in $\mathcal{R}$ that we find for the main peak in the $P(\mathcal{R})$ spectra are indeed strong enough so that we are sensitive to the 0-angular frequency that is emitted at these Faraday depths.

The main peak in the $P(\mathcal{R})$ spectrum is not measurably resolved for 8% of the lines of sight in our data set. An unresolved peak means that there is no synchrotron emission occurring in the Faraday-rotating region. This is unexpected, since cosmic rays that produce synchrotron emission pervade the ISM. We propose that the magnetic field orientation plays an important role in this, in that a magnetic field that is oriented along the line of sight produces only synchrotron emission perpendicular to the line of sight, hiding it from our view. By using the observed emission measures from the WHAM survey as input, and by simulating the thermal electron contribution to $\mathcal{R}$, we could determine the magnetic field strength along the line of sight, and we mapped this quantity over the area covered by our GEMINI data. The polarized point sources we found in our data have RMs that are comparable to the $\mathcal{R}$ we find for the diffuse emission. In our $P(\mathcal{R})$ maps we found depolarization canals, narrow structures with only a small fraction of the polarized intensity of adjacent lines of sight. We established that for a number of deep canals the polarization angle changes in many frequency channels by a large amount over the canal, in some cases 90$^\circ $. We therefore argue that depolarization over the synthesized telescope beam is producing (at least) these canals. Finally, we investigated the properties of a large, conspicuous area in our data, and we argue that a decrease in depolarization along the line of sight as compared to its surroundings is probably responsible for the observational features of this region.

Acknowledgements
DS would like to thank the referee, Rob Reid, for his useful comments on the original manuscript, and Melanie Johnston-Hollitt for supplying the ionospheric RM data. The Westerbork Synthesis Radio Telescope is operated by The Netherlands Foundation for Research in Astronomy (NFRA) with financial support from the Netherlands Organization for Scientific Research (NWO). The Wisconsin H-Alpha Mapper is funded by the National Science Foundation.

References

Footnotes

... website[*]
http://astro.berkeley.edu/~marc/dust/data/data.html

All Tables

Table 1:   Characteristics of the GEMINI data set. Observing dates and times are given for each of the 12 h observing runs, which have been indicated by their shortest baseline length. We calculated the conversion factor between mJy/beam and K at 345 MHz, the average of the $\lambda ^2$ sampling of the (usable) frequency channels in our data set.

Table 2:   Properties of the polarized point sources we found in our data. Shown are the equatorial coordinates (J2000.0) of the source in decimal degrees, the RM that we fitted and its error (both in rad/m2), the reduced $\chi ^2$ ( $\chi ^2_{{\rm red}}$) of the fit, the maximum $P(\mathcal{R})$ along the line of sight and the total intensity I (both in mJy), and the polarization fraction (in %).

All Figures

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig1.ps}
\end{figure} Figure 1:

$P(\lambda ^2)$ and $\Phi (\lambda ^2)$ distributions ( top 2 panels), and $P(\mathcal{R})$ and $\Phi(\mathcal{R})$ spectra ( bottom 2 panels) for a line of sight with a barely resolved main peak that is similar to a Faraday screen. In the text we describe how we derotated the polarization angles in the bottom panel of this figure. The $\chi ^2_{{\rm red}}$ of the RM fit is 0.99 (determined by using the RM with the maximum P(RM)). The $\mathcal{R}$ spectra are sampled from -100 rad/m2 to +100 rad/m2 in steps of 1 rad/m2. In the bottom 2 panels, solid lines indicate the calculated $P(\mathcal{R})$ and $\Phi(\mathcal{R})$ spectra, and the dashed lines indicate the fitted RMSF.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig2.ps}
\end{figure} Figure 2:

Identical to Fig. 1, but for a line of sight towards a polarized extragalactic source that also contains diffuse Galactic emission. The  $\chi ^2_{{\rm red}}$ of the RM fit is 8.72.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=17cm,clip]{8912fig3.ps}
\end{figure} Figure 3:

Slices through the $P(\mathcal{R})$ cube that show the strongest diffuse emission. The $\mathcal{R}$ of each slice is indicated in rad/m2 in the upper right corner of each image. The RMSF has a FWHM of 12 rad/m2, and we sampled each RMSF with 2 images. The linear intensity scale is the same for all panels, and it saturates at 6.4 K.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig4.ps}
\end{figure} Figure 4:

$P(\mathcal{R})$ of the highest peak in the $\mathcal{R}$ spectrum. Each pixel represents one telescope beam, and the pixel size is 2.8$^\prime $ $\times $ 4.9$^\prime $. The lines of sight that are shown in this plot have a second peak in the $P(\mathcal{R})$ spectrum that is at most half the strength of the main peak. Out of a total of 22.800 lines of sight, 2720 do not obey this selection criterion, and they are shown in black, as are the 769 lines of sight that have $P(\mathcal{R})$ < 2 K. The 1$\sigma $ noise level in the $P(\mathcal{R})$ maps is 0.14 K.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig5.ps}
\end{figure} Figure 5:

$\mathcal{R}$ corresponding to the $P(\mathcal{R})$ plotted in Fig. 4. The lines of sight plotted in these figures are selected by using the same criterion. Lines of sight indicated in white have $\mathcal{R}$ > 30 rad/m2. Lines of sight indicated in black either did not pass the selection criterion or have $\mathcal{R}$ < -6 rad/m2 (2720 and 727 lines of sight respectively).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{8912fig6.ps}
\end{figure} Figure 6:

Distribution of independent lines of sight with a main peak that satisfies the $\Delta $ criterion (green pixels), or with $\chi ^2_{{\rm red}}$ < 2 and max( $P(\mathcal{R})$) > 2 K (yellow pixels), or that satisfy all three criteria (red pixels). Lines of sight that did not pass these criteria are shown as black pixels. The outer edges of the mosaic show an increase in  $\chi ^2_{{\rm red}}$, which explains the absence of yellow pixels.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig7.ps}
\end{figure} Figure 7:

WHAM H$\alpha $ intensities in Rayleigh ( top panel) and the electron-density weighted strengths of the magnetic field component parallel to the line of sight $\langle B_{\Vert}\rangle$ [$\mu $G] that we derived ( bottom panel). Open circles indicate positive  $\langle B_{\Vert}\rangle$; there are no negative $\langle B_{\Vert}\rangle$ < -0.1 $\mu $G in this figure. $\vert\langle B_{\Vert}\rangle\vert$ < 0.1 $\mu $G are shown as ``+''.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{8912fig8.ps}
\end{figure} Figure 8:

RM distribution for the polarized point sources from Table 2. The size of the circles is proportional to RM, shown in the scale on top of the figure in units of radians/m2. Black circles indicate negative RM, and white circles indicate positive RM. |RM|  $\protect\la$ 4 rad/m2 are missing from this figure because of our selection criteria. We excluded the grey region when looking for polarized point sources, because instrumental polarization levels are higher there.

Open with DEXTER
In the text

  \begin{figure}
\par {\hspace*{5mm}\includegraphics[width=15.5cm,clip]{8912fig9.p...
...mm}}\vspace*{3mm}
\includegraphics[width=16.5cm,clip]{8912fig10.ps}
\end{figure} Figure 9:

Depolarization canals in the $P(\mathcal{R})$ datacube. The $\mathcal{R}$ [rad/m2] of each of the $P(\mathcal{R})$ slices is indicated in the lower left. All maps have the same linear greyscale that saturates at 4.8 K. The area shown in the top row of this figure lies in the upper left of Fig. 3. For the 2 pixels that are indicated by crosses, and that lie in depolarization canals, we plot the $P(\lambda ^2)$ and $\Phi (\lambda ^2)$ spectra in the first row of Fig. 10 (left cross, ``CH 1'' in Fig. 10) and middle row of Fig. 10 (right cross, ``CH 2''). The bottom row shows a canal in a different part of the mosaic, and the $P(\lambda ^2)$ and $\Phi (\lambda ^2)$ spectra of the canal pixel that is indicated by a cross are shown in the bottom row of Fig. 10 (``CH 3'').

Open with DEXTER
In the text

  \begin{figure}
\includegraphics[width=16.3cm,clip]{8912fig11.ps}
\end{figure} Figure 10:

Polarized brightness temperature  $T_{{\rm b,pol}}$ [K] and polarization angle $\Phi $ [degrees] spectra for 3 canal pixels, and two adjacent pixels half a beamwidth from each canal pixel. Top row: canal pixel at ( $\alpha ,\delta $) = (111.24$^\circ $, 39.72$^\circ $), and two adjacent lines of sight at $\alpha $ = 111.27$^\circ $ (dotted line), and $\alpha $ = 111.21$^\circ $ (dashed line). Middle row: canal pixel at ( $\alpha ,\delta $) = (110.98$^\circ $, 39.57$^\circ $), and adjacent lines of sight at $\alpha $ = 111.01$^\circ $ dotted line, and $\alpha $ = 110.95$^\circ $, dashed line. Bottom row: canal pixel at ( $\alpha ,\delta $) = (111.50$^\circ $, 34.94$^\circ $), and two adjacent lines of sight at $\alpha $ = 111.53$^\circ $, dotted line, and $\alpha $ = 111.47$^\circ $, dashed line. The 1$\sigma $ errors for polarized intensity and polarization angle are 2 K and 15$^\circ $ (for a S/N = 2 detection) respectively Every 8 ${{\rm th}}$ channel is plotted. Bad frequency channels were selected using the same selection criteria as in Sect. 3, and discarded.

Open with DEXTER
In the text


Copyright ESO 2009

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.