Narrow depolarized canals are common in maps of the polarized synchrotron emission of the Milky Way. Two physical effects that can produce these canals have been identified: the presence of Faraday rotation measure (RM) gradients in a foreground screen, and the cumulative cancellation of polarization known as differential Faraday rotation. We show that the behaviour of the Stokes parameters Q and U in the vicinity of a canal can be used to identify its origin. In the case of canals produced by a Faraday screen we demonstrate that, if the polarization angle changes by 90° across the canal, as is observed in all fields to-date, the gradients in RM must be discontinuous. Shocks are an obvious source of such discontinuities and we derive a relation of the expected mean separation of canals to the abundance and Mach number of supernova-driven shocks, and compare this with recent observations by Haverkorn, Katgert & de Bruyn. We also predict the existence of less common canals with polarization angle changes other than 90°. Differential Faraday rotation can produce canals in a uniform magneto-ionic medium, but as the emitting layer becomes less uniform the canals will disappear. We show that for moderate differences in emissivity in a two-layer medium, of up to 1/2, and for Faraday depth fluctuations of standard deviation ≲1 rad, canals produced by differential rotation will still be visible.
Polarized (synchrotron) radio emission is a rich source of information about the relativistic and thermal plasmas and magnetic fields in the interstellar medium (ISM). Recent observations have revealed an abundance of unexpected features that arise from the propagation of the emission through the turbulent ISM (Wieringa et al. 1993; Uyaniker et al. 1998; Duncan et al. 1999; Gray et al. 1999; Haverkorn, Katgert & de Bruyn 2000; Gaensler et al. 2001; Wolleben et al. 2006). Arguably the most eye-catching is the pattern of depolarized canals: a random network of dark narrow regions, clearly visible against a bright polarized background. These canals evidently carry information about the ISM, but it is still not quite clear how this information can be extracted. Two theories for the origin of the canals have been proposed; both attribute the canals to the effects of Faraday rotation, but one invokes steep gradients of the Faraday rotation measure (RM) across the telescope beam, in a Faraday screen (Haverkorn et al. 2000; Haverkorn, Katgert & de Bruyn 2004), whereas the other relies on the line-of-sight effects producing differential Faraday rotation (Beck 1999; Shukurov & Berkhuijsen 2003). In order to use the canal properties to derive parameters of the ISM, one must correctly identify their origin.
In this Letter we briefly discuss a few important aspects of the two theories describing the origin of canals. Detailed discussion of relevant depolarization mechanisms is presented in Fletcher & Shukurov (2006, and in preparation). Canals produced by Faraday rotation measure gradients are discussed in Section 2 where we show that these canals require a discontinuous distribution of free electrons and/or magnetic field; we further suggest an interpretation of the discontinuities in terms of interstellar shocks. The case of differential Faraday rotation is discussed in Section 3; in Section 4 we suggest an observational test that can be used to identify the specific mechanism that produces a given canal.
The defining features of the canals are as follows:
the observed polarized emission P approaches the polarization noise level σP, P≲σP;
the canal is about one beam wide;
the canal passes through a region of significant polarized intensity, say P≳ 3σP;
the canal is not related to any structure in total intensity, and so cannot be readily explained by e.g. an intervening gas or magnetic filament.
Polarized intensity vanishes when emission within the telescope beam consists of two equal parts with mutually orthogonal polarization planes. Thus feature (i) most often arises because the polarization angle Ψ changes by 90° across the canal. However, in Section 2 we predict a new type of canal across which the observed polarization angle does not change.
We only consider canals occurring in properly calibrated maps: Haverkorn et al. (2004) argue that the observations we discuss in Section 2.1 do not suffer from missing large-scale structure; Reich (2006) discusses the calibration of radio polarization maps in depth.1
The value of RM is related to F, but often in a complicated manner (see e.g. Burn 1966; Sokoloff et al. 1998). Simplest is the case of a Faraday screen, where the source of synchrotron emission is located behind a magneto-ionic region (e.g. because relativistic and thermal electrons occupy disjoint regions): then RM =Fλ−2. In a homogeneous region, where relativistic and thermal electrons are uniformly mixed, RM = 0.5Fλ−2.
Canals Produced by a Faraday Screen
Consider polarized emission from a uniform background source passing through a Faraday screen – a layer where no further emission occurs, but which rotates the polarized plane. If F varies with r⊥, adjacent lines of sight within the beam are subject to different amounts of Faraday rotation and the observed degree of polarization decreases. If the variation of F within the beam, ΔFD (Fig. 1), produces a 90° difference in Ψ, one might expect that the depolarization will be complete as illustrated in Fig. 2(a), and a canal will be observed along contours defined by ΔFD= (n+ 1/2)π with n= 0, 1, 2, … .
The situation is more subtle, though. Fig. 2(b) shows Ψ changing smoothly by 90° across a beam, as a result of a monotonic change in Faraday depth by the same amount, ΔFD=π/2, as in panel (a). It is obvious that polarized emission does not cancel within the beam; the polarized intensity will be lower than for a uniform arrangement of angles but there will still be a polarized signal detected with a polarization angle of about −45°.
Thus gradients in F across the beam can produce complete depolarization in a Faraday screen if F changes by ΔFD≃ (n+ 1/2)π within a small fraction of the beamwidth. As shown in Fig. 3 (solid line), obtained using equations (2) and (4), the depth of a canal will only be less than 10 per cent of the surrounding polarized emission if the gradient in F occurs over one-fifth or less of the beam: Δx/D≲ 0.2, where Δx is the extent of the region where the gradient occurs and D is the beamwidth. Haverkorn et al. (2004) reached a similar conclusion in an analysis of the depth and profile of canals observed at λ84 cm (see their section 4.4).
As illustrated in Figs 2(c) and 3 (dotted line), a variation in F of a magnitude ΔFD=π across the FWHM of the Gaussian beam (i.e. ΔF=π and Δx/D≃ 1) can also produce strong depolarization but now the change in polarization angle across the beam is ΔΨ≃ 0°. Similarly, other gradients in F can significantly reduce the degree of polarization: for example, ΔF= 0.7π (dashed line in Fig. 3) across a region that is about 0.7× the beamwidth will produce complete depolarization and a change in angle of 54° across the beam.
Larger gradients in F also result in strong depolarization and when ΔF > π the degree of depolarization is less sensitive to the resolution. This is illustrated in Fig. 3 for the case where ΔF= 1.5π (dash–dotted line); complete depolarization occurs when Δx/D= 0, 1.5 but in the range 0 < Δx/D≲ 2 the degree of polarization remains below the 10 per cent level. Thus an increment of ΔF= 1.5π, and more generally ΔF= (n+ 1/2)π, will produce one-beam-wide canals when the gradient occurs in a region roughly equal to or less than the beamwidth, i.e. where Δx/D≲ 1. Other large gradients in F, such as ΔF=nπ, depolarize to around 10 per cent of the background p and therefore can also produce canals; however, in this case one-beam-wide canals will only be formed when Δx/D≃ 1.
So, if F is a continuous, random function of position (with sufficiently strong fluctuations on the scale of the beam) we predict the existence of canals with every change in polarization angle 0° < ΔΨ < 90° across the beam. On the other hand a discontinuous change in F by (n+ 1/2)π will produce canals across which the polarization angle changes by ΔΨ= 90°. This insight is important, as to-date all observed canals have ΔΨ≃ 90°. The only other way in which Faraday screens can consistently produce canals with this change in angle is if all the gradients in Faraday depth on the approximate scale of the beamwidth are of the magnitude ΔF≃ (n+ 1/2)π, n= 1, 2, … (n= 0 in this case only produces weak depolarization, as noted above); this is extremely unlikely in a random medium.
If a system of observed canals has ΔΨ= 90° and is caused by a foreground Faraday screen – see Section 4 for a diagnostic test for the origin of canals – then the distribution of rotation measure in the screen must be discontinuous on the scale of the beam. Haverkorn & Heitsch (2004) used magnetohydrodynamic simulations of Mach 10 turbulence to show that RM gradients can be steep enough to produce depolarized canals on smoothing with a sufficiently large simulated beam. Note that a finite beamsize is essential for the production of canals by this mechanism; with perfect resolution the width of the shock front will be resolved and no canal will appear. Furthermore, these canals will not fill up if the data are smoothed (solid and dash–dotted curves for Δx/D→ 0 in Fig. 3) whereas canals produced by gradients other than ΔF= (n+ 1/2)π will disappear under smoothing at a rate shown with dashed and dotted curves in Fig. 3.
The mean separation of canals produced by shocks
An obvious source of Faraday depth discontinuities is shocks in the ISM. Then the mean distance between the canals produced by a Faraday screen will be related to the distance between the shock fronts. The expected separation of suitable shock fronts can be derived using the model of interstellar shock-wave turbulence of Bykov & Toptygin (1987) and then compared with the separation of canals observed by Haverkorn et al. (2003), claimed to arise in a Faraday screen.
For supernova-driven shocks, Bykov & Toptygin (1987) derived
If the magnetic field is frozen into the gas and the gas density increases by a factor of ερ at the shock then so will the field strength (we neglect a factor arising from the relative alignment of the shock and field). A canal forms where F increases discontinuously by ΔFD= (n+ 1/2)π. So the required density compression ratio isLandau & Lifshitz 1960)
In the field of canals observed by Haverkorn et al. (2003) at λ84 cm, 〈RM〉≃−3.4 rad m−2 and a by-eye estimate of their mean separation gives arcmin. The most abundant canals are produced by the shocks which can generate ΔFD=π/2 and from equation (12) this requires . Using fcl= 0.25, taking α= 4.5 and the parameter values used to normalize equation (9), we obtain L≃ 10 pc for shocks with . Our estimate of L includes shocks with that will not produce canals, but the strong dependence of in equation (8) on Mach number means that L will underestimate the separation of canal-generating shocks insignificantly. Using equation (7), with L≃ 10 pc and arcmin, we find that the canals observed by Haverkorn et al. (2003) are compatible with a system of shocks occurring in a Faraday screen with a depth of d≃ 100 pc. The maximum distance in this field, beyond which emission is completely depolarized, is estimated to be ∼600 pc by Haverkorn et al. (2003). Thus this model for the origin of the canals requires that the nearest 100 pc in the direction (l, b) = (161°, 16°) is effectively devoid of cosmic ray electrons; if the cosmic rays normally spread over a distance vAt∼ 1 kpc, where vA≃ 20 km s−1 is the Alfvén speed and t≃ 3 × 107 yr their lifetime, such a condition is difficult to explain.
Canals Produced by Differential Faraday Rotation
In Section 2 we discussed the depolarization effects of Faraday rotation measure gradients transverse to the line of sight acting on smooth polarized background emission. Now we will consider a uniform layer in which both emission and Faraday rotation occur. This gives rise to the well-known effect of differential Faraday rotation (Burn 1966; Sokoloff et al. 1998), where polarized emission from two positions along the line of sight separated by a rotation of Δφ=π/2 exactly cancels, thus reducing the degree of polarization. When a line of sight has a total rotation of F=nπ there is total cancellation of all polarized emission in the layer and the degree of polarization is p= 0. Canals produced by this mechanism are discussed by Shukurov & Berkhuijsen (2003).
Differential Faraday rotation in a non-uniform medium
We now investigate how canals produced by differential Faraday rotation are affected by deviations from uniformity in the magneto-ionic medium. First we will discuss the case of a two-layer medium in which the synchrotron emissivity is different in each layer. Then we allow for random fluctuations of F (the latter case produces what is known as Faraday dispersion).Sokoloff et al. (1998), but with typos corrected.] We assume that F=F1+F2=π, the condition for canals to form as a result of differential Faraday rotation, and investigate what will happen to the canals if I1≠I2. We parametrize the difference in the synchrotron emissivity in the two layers as ε= (I2−I1)/I1, choose equal Faraday rotation through each layer for simplicity so that F1=F2=π/2, and substitute these into equation (13). We then obtain the dependence of the minimum degree of polarization on ε: Shukurov & Berkhuijsen 2003). Where the difference in emissivity is greater, say ε∼ 2 or more as one might expect viewing distant bright spiral arm emission near the Galactic plane, the canals will more readily fill up, and become much less distinct and more difficult to interpret confidently.
Now let us consider the effect of random fluctuations in Faraday rotation, sometimes called internal Faraday dispersion. We start with equation (34) of Sokoloff et al. (1998):Sokoloff et al. 1998): for large enough σF depolarization is complete. However, as long as σF < 1, canals produced by differential Faraday rotation will only fill up by one-third or less. For example, at λ20 cm the canals will not be destroyed as long as the dispersion in Faraday rotation measure is less than 30 rad m−2. This is why canals are still visible in the λ20-cm map studied by Shukurov & Berkhuijsen (2003) which has a dispersion in rotation measure of about σRM≃ 10 rad m−2, i.e. σF=σRMλ2≃ 0.4.
Observational Diagnostics for the Origin of a Canal
The complex polarization emitted by a layer producing differential Faraday rotation can be written as (Burn 1966)
The main points of this paper can be summarized as follows.
The behaviour of the Stokes parameters Q and U in the vicinity of a canal allows one to identify whether a foreground Faraday screen or differential Faraday rotation is the cause of the canal (Section 4).
A foreground Faraday screen can produce canals with any polarization angle change across the canal, 0 < ΔΨ < 90°. However, discontinuous jumps in the Faraday depth will only produce canals with ΔΨ≃ 90° (Section 2).
If shocks produce the discontinuities in a foreground Faraday screen that generate canals, the mean separation of the canals can provide information about the Mach number and separation of shocks in the screen (Section 2.1).
Canals produced by differential Faraday rotation are sensitive to non-uniformity in the medium along the line of sight, systematic or random. However, they will remain recognizable if the synchrotron emissivity varies by less than a factor of about 2 or if the standard deviation of the Faraday depth is σF < 1 (Section 3.1).
This work was supported by the Leverhulme Trust under research grant F/00 125/N.