EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 586 (February 2016) A&A, 586 (2016) A141 Full HTML
Free Access
Issue
A&A
Volume 586, February 2016
Article Number A141
Number of page(s) 17
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201526506
Published online 09 February 2016

© ESO, 2016

1. Introduction

Recently, Planck1 has reported an asymmetry in power between the dust E- and B-modes in its 353 GHz observations (Planck Collaboration Int. XXX 2016; Planck Collaboration X 2015). This power asymmetry has been observed outside masks covering 20 to 70% of the sky, excluding the Galactic plane. The ratio of the dust B- to E-mode power amplitudes is about a half over the multipole range 40 <ℓ< 600 (Planck Collaboration Int. XXX 2016). The source of this power asymmetry in the dust polarization data is currently unknown. Models of the Galactic magnetic field (GMF) used in the Planck sky model (PSM, Delabrouille et al. 2013) and the FGPol model (O’Dea et al. 2012) produce an equal amount of power in E- and B-modes outside the regions covered by the sky masks. These models, which were used to estimate the dust polarization foreground (BICEP2 Collaboration 2014; BICEP2/Keck Array and Planck Collaborations 2015), include an analytical model of the large-scale GMF (BGMF) and a statistical description of the turbulent component of the magnetic field.

The Planck maps of thermal dust emission display filaments distributed over the whole sky (Planck Collaboration XI 2014). The filamentary structure of the diffuse interstellar matter is also a striking feature of dust observations at higher angular resolution, performed by Herschel, and of spectroscopic Hi observations (e.g. Miville-Deschênes et al. 2010; André et al. 2014; Clark et al. 2014). The analysis of Planck dust polarization data in the diffuse interstellar medium (ISM), at low and intermediate Galactic latitudes, indicates that the structures of interstellar matter tend to be aligned with the plane of the sky (POS) projection of the magnetic field (BPOS, Planck Collaboration Int. XXXII 2016). This preferential relative orientation is also observed in simulations of magneto-hydrodynamic (MHD) turbulence of the diffuse ISM (Hennebelle 2013; Soler et al. 2013). Such a coupling between the structure of interstellar matter and BPOS is not included in the PSM or FGPol models of the dust polarization sky (Delabrouille et al. 2013; O’Dea et al. 2012).

The goal of this paper is to test whether the correlation between the filamentary structures of the intensity map and BPOS in the diffuse ISM accounts for the observed E−B asymmetry. Zaldarriaga (2001) describes the E- and B-modes decomposition of simple patterns of polarized emission, including filaments with a homogeneous polarization degree and orientation. The presence of E-modes is related to invariance by parity of the polarization pattern. There is E-only power if BPOS is either parallel or perpendicular to the filaments. If BPOS is oriented at + 45° or − 45° with respect to the filaments, there is B-only power.

In this paper, we filter the Planck intensity and polarization maps to isolate filaments in the range of angular scales where the E−B asymmetry is observed. We identify coherent elongated filaments within regions of low column density at high Galactic latitude using a filament-finding algorithm. We evaluate the mean polarization angle in each of these filaments and compare it to the mean orientation of each filament. In doing so, we extend the analysis presented in Planck Collaboration Int. XXXII (2016) to the relevant region of the sky for CMB polarization observations at high Galactic latitude. In order to achieve a high signal-to-noise ratio and enhance the contrast with respect to the local background dust emission, we stack the Stokes I, Q, U, and also E and B maps, for the filaments we select in the Planck dust intensity map. We use the stacked images to quantify the power asymmetry in E- and B-modes associated with the filaments.

This paper is organized as follows. In Sect. 2, we introduce the Planck353 GHz data used in this study. The filament-finding algorithm is presented in Sect. 3. Section 4 presents the study of relative orientation between the filaments, BPOS, and the POS component of the large-scale GMF (Bm,POS) at high Galactic latitude. In Sect. 5, we present the stacking of both intensity and polarization maps and derive the mean polarization fraction of the filaments. In Sect. 6, we discuss the relation between the relative orientation of the filaments and BPOS and the E−B asymmetry. Section 7 presents our results in the context of earlier studies and its relation to Galactic astrophysics. Finally, we present our conclusions in Sect. 8. This paper has four appendices. Appendix A details the Hessian analysis implemented to identify the filaments in the dust intensity map. The application of the filament-finding algorithm to a simulated Gaussian dust sky is detailed in Appendix B. In Appendix C, we study the impact of our selection of the filaments on the main results of the paper. The computation of all the angle uncertainties that we use in our analysis is presented in Appendix D.

2. Planck data

The Planck satellite has observed the sky polarization in seven frequency bands from 30 to 353 GHz (Planck Collaboration I 2015). In this paper, we only use the 2015 (“DX11d”) data from the High Frequency Instrument (HFI, Lamarre et al. 2010) at 353 GHz, since they are best suited to study the structure of dust polarization (Planck Collaboration Int. XIX 2015; Planck Collaboration Int. XX 2015; Planck Collaboration Int. XXI 2015; Planck Collaboration Int. XXII 2015). The data processing, map-making, and calibration of the HFI data are described in Planck Collaboration VII (2016) and Planck Collaboration VIII (2016). In our analysis, we ignore the dust and CO spectral mismatch leakage from intensity to polarization (Planck Collaboration VIII 2016). Planck Collaboration Int. XXX (2016) has shown that the amplitude of the dust spectral mismatch leakage at high latitude (fsky = 0.5) is small compared to the total polarization signal in E and B modes. No CO emission is detected at 353 GHz away from the Galactic plane and the brightest molecular clouds (Planck Collaboration XIII 2014) and so we do not consider it in our analysis.

To quantify the statistical noise and systematic effects on the results presented in this paper, we use the two HalfMission (HM), two DetSet (DS), and two HalfRing (HR) Planck353 GHz polarization maps (Planck Collaboration VII 2016). The two HM maps are made from the two halves of the full-mission Planck data, whereas the two HR maps are produced by splitting each ring (also called stable pointing period) into two equal duration parts. The two DS maps are constructed using two subsets of polarization-sensitive bolometers at a given frequency. The noise is uncorrelated between the two HM, HR, and DS maps. We only use them to compute the error bars on the relevant quantities that we measure in this paper.

The total polarization intensity (P353) and the polarization angle (ψ) are derived from the full-mission Stokes Q353 and U353 maps at 353 GHz using the relations where the two-argument function atan2(−U353,Q353) is used to compute atan(−U353/Q353) avoiding the π ambiguity. To recover the correct full range of polarization angles ([− π/ 2,π/ 2] as used for ψ here), attention must be paid to the signs of both U353 and Q353, not just their ratio. We use the IAU convention for ψ, which is measured from the Galactic North (GN) and positive to the East. The minus sign in Eq. (2)converts the convention provided in the Planck data to that of the IAU (see Planck Collaboration Int. XIX 2015). The orientation angle (χ) of BPOS is defined within the π ambiguity by adding π/ 2 to the polarization angle (3)For the dust intensity at 353 GHz, we use the model map D353, computed from a modified blackbody fit to the Planck data at ν ≥ 353 GHz and the IRAS 100  μm map (Planck Collaboration XI 2014). This map has lower noise than the 353 GHz Stokes I map and is corrected for zodiacal light emission, CMB anisotropies, and the cosmic infrared background monopole. We neglect the contribution of the CMB polarization at 353 GHz for this study.

The full-mission Planck Stokes Q353 and U353 maps are provided in HEALPix2 format (Górski et al. 2005) at 4.́8 resolution and D353 at 5′. To increase the signal-to-noise ratio, we smooth the three maps to a common resolution of 15′, taking into account the effective beam response of each map, and reduce to a HEALPix resolution of Nside = 512. For the polarization data, we decompose the Stokes Q353 and U353 maps into E353 and B353aℓms (Eℓm and Bℓm) using the “ianafast” routine of HEALPix, apply the Gaussian smoothing in harmonic space (after deconvolving the effective azimuthally symmetric beam response of each map), and transform the smoothed E353 and B353aℓms back to Q353 and U353 maps using the “isynfast” routine at Nside = 512. We also transform the E353 and B353aℓms to E353 and B353 maps at Nside = 512 using the relations (4)All the maps that we use are in thermodynamic units (μKCMB).

In this paper, we work with the bandpass-filtered dust intensity map, , to identify and isolate filaments over the filtering scale using a filament-finding algorithm. By filtering out large-scale and small-scale modes, we enhance the contrast of the filaments with respect to the diffuse background and reduce the instrumental noise, which is critical for accurately measuring the polarization orientations of the filaments within regions of low column density at high Galactic latitude.

thumbnail Fig. 1

The representative bandpass filter, retaining only the scales between ℓ = 30 and 300.

For filtering, we apply the three-dimensional spline wavelet decomposition based on the undecimated wavelet transform, as described by Starck et al. (2006). We use the publicly available package Interactive Sparse Astronomical Data Analysis Packages (ISAP3) to compute the map at Nside = 512 resolution. The spline wavelet used in this analysis provides less oscillation in position space compared to Meyers or needlet ones (Lanusse et al. 2012). The filtering is done in pixel space; the corresponding bandpass filter in harmonic space varies a little over the sky. Figure 1 presents the typical shape of this bandpass filter, which selects the scales between ℓ = 30 and 300. The filtering scale is chosen in such a way that it highlights all the bright filaments present in the PlanckD353 map. We also compute the bandpass-filtered polarization maps, , , , and .

thumbnail Fig. 2

Data processing steps implemented to identify filaments from the Planck data. We start with the PlanckD353 map (upper left panel) smoothed at 15′ resolution. The bandpass-filtered map (upper right panel) is produced using the spline wavelet decomposition, retaining only the scales between ℓ = 30 and 300. The lower eigenvalue map of the Hessian matrix, λ−, is shown in the lower left panel. Structures identified in the high-latitude sky λ− map are shown in the lower right panel. The superimposed graticule is plotted in each image and labelled only on the lower right panel. It shows lines of constant longitude separated by 60° and lines of constant latitude separated by 30°. The same graticule is used in all plots of the paper.

3. Filament-finding algorithm

3.1. Methodology

Identification of filaments as coherent structures is a crucial part of this analysis. Previous studies used algorithms such as DisPerSE (Sousbie 2011; Arzoumanian et al. 2011), getfilaments (Men’shchikov 2013), and the rolling Hough transform (Clark et al. 2014). Hennebelle (2013) and Soler et al. (2013) used the inertia matrix and the gradient of the density and column density fields to identify filaments in numerical simulations of MHD turbulence.

In this paper, we employ the Smoothed Hessian Major Axis Filament Finder (SMAFF, Bond et al. 2010a) algorithm, which has been used to identify filaments in the three-dimensional galaxy distribution (Bond et al. 2010b). SMAFF is primarily based on the Hessian analysis. The Hessian analysis has also been used to analyse the Planck dust total intensity map in Planck Collaboration Int. XXXII (2016), Herschel images of the L1641 cloud in Orion A (Polychroni et al. 2013), and large-scale structure in simulations of the cosmic web (Colombi et al. 2000; Forero-Romero et al. 2009). Planck Collaboration Int. XXXII (2016) has reported good agreement between the filament orientations derived from the Hessian and inertia matrix algorithms.

3.2. Implementation

In this study, we apply the two-dimensional version of SMAFF to the map, which is shown in the upper right panel of Fig. 2. From the Hessian matrix, we compute an all-sky map of the lower eigenvalue λ− and the orientation angle θ− of the perpendicular to the corresponding eigenvector, measured with respect to the GN. The details of the Hessian analysis are provided in Appendix A. The map of λ− is presented in the lower left panel of Fig. 2. For the subsequent analysis, we consider only the high-latitude sky, defined as | b | > 30°, with the Large Magellanic Cloud and Small Magellanic Cloud regions masked out.

thumbnail Fig. 3

Distribution of the eigenvalues λ− over the unmasked pixels in the high-latitude sky. The grey region represents the pixels that were used in the SMAFF algorithm to find strong filaments.

The map of λ− highlights filaments in the map with an orientation angle θ−, which we refer to as the Hessian angle hereafter. The distribution of λ− over the unmasked pixels is shown in Fig. 3. This distribution of λ− is non-Gaussian with an extended tail. We use the median absolute deviation (MAD, Hampel 1974; Komm et al. 1999) to measure the width, σλ−, of the distribution, as given by (5)where mλ− is the median of the λ− distribution. We select filaments using an upper threshold (K−) on λ− given by (6)Hereafter, we refer to the filaments satisfying λ−<K− as “strong”. This threshold K− separates the strong filaments from the weak ones, as detailed in Appendices B and C. By construction, the threshold K− rejects pixels where λ− is positive, since those pixels do not correspond to local maxima.

We seek coherent elongated structures in the map. In SMAFF, this is achieved by placing an upper limit C on the difference between Hessian angles within a given structure. For our purpose, we set the value of C = 15° to identify relatively straight filaments.

We start with the pixel having the most negative λ− and denote the corresponding Hessian orientation angle by . We identify its neighbouring pixels using the “neighbours_ring” routine of HEALPix and look for pixels with λ−<K− and orientation angle such that . If both conditions are satisfied, we count that neighbouring pixel as a part of the filament and move on to that neighbouring pixel. The neighbouring pixel becomes the new reference point and we search for its neighbours that satisfy both conditions. In our algorithm, is fixed by the starting pixel, which has the most negative λ−. We continue this friend-of-friend algorithm to connect pixels until one of the conditions is no longer satisfied. We limit our selection to filaments with a length (L, defined as the maximum angular distance between pixels within a given structure) larger than or equal to the threshold length L0, which we choose to be 2°. This process yields a set of 259 elongated filaments, as shown in the lower left panel of Fig. 2. Hereafter, we refer to this set as our filament sample. Selected sky pixels represent 2.2% of the high-latitude sky considered in our analysis. There is no overlap between the filaments in our sample.

The column density is computed from the map using the conversion factor, 0.039 MJy sr-1 per 1020H cm-2; this was derived in Planck Collaboration Int. XVII (2014) by correlating the Planck353 GHz dust total emission map with an Hi column density map over the southern Galactic polar cap. We average the column density along each filament and assign one mean column density, , to each. This column density is computed on the filtered intensity map. The histogram of for the filament sample is presented in Fig. 4. The number of filaments per is represented by NF.

thumbnail Fig. 4

Histogram of the mean column density of the filament sample. The column density is computed from the map using the conversion factor derived in Planck Collaboration Int. XVII (2014).

4. Interplay between the filament orientation and the magnetic field

In this section, we study the orientations of matter structures and BPOS in our filament sample (Sect. 3). The orientation angle of BPOS is derived from the observed Stokes Q353 and U353 maps using Eqs. (2)and (3). We also consider the orientation angle () of Bm,POS, as estimated from starlight polarization observations (Heiles 1996) and pulsar rotation measures (Rand & Lyne 1994; Han et al. 1999). We compare these three orientations, as represented in Fig. 5. Our analysis follows Planck Collaboration Int. XXXII (2016), which used a set of pixels representing approximately 4% of the sky at low and intermediate Galactic latitudes. Only 25% of the pixels in our current filament sample were considered in this earlier study.

thumbnail Fig. 5

Sketch of the mean orientation angle of the filament (), the magnetic field (), the polarization angle (), and the large-scale GMF () along the filament. All the angles are defined with respect to the GN and follow the IAU convention.

4.1. Relative orientation of the filaments and the magnetic field

We study the angle difference between the orientations of the filaments and BPOS in our sample. First we associate one POS orientation angle with each of the filaments with respect to the GN. By construction, due to our selection criteria on the angles, the filaments are fairly straight and, hence, they may be described with a single orientation angle. Given one filament, we measure the mean orientation angle, , over the n pixels that belongs to it. We make use of the pseudo-vector field with unit length computed from the values of θ− for each pixel. This pseudo-vector has components Q− = cos 2θ− and U− = −sin 2θ− (following the HEALPix convention for the Q− and U− components). The mean POS orientation angle of the filament is obtained by first averaging Q− and U− over all n pixels and then calculating the position angle of this averaged pseudo-vector. It is given by (7)If we rotate the Stokes Q353 and U353 maps by , i.e., into the frame where the axis of the filament is in the North-South direction, the rotated Q353 and U353 can be written as Combining Eqs. (8)and (9)with Eqs. (1)and (2), we get where the orientation angle χ is defined in Eq. (3). Similar to the computation of , we average and over all n pixels and then calculate the position angle of this averaged pseudo-vector where

The angle difference measures the weighted mean of the angle difference per pixel between the orientations of the given filament and BPOS. The index refers to the mean orientation angle of BPOS along the filament. Note that we directly measure the angle difference between the filament and BPOS, without computing for each filament.

The histogram of relative orientation (HRO) between the filament and BPOS for our filament sample is presented in the upper panel of Fig. 6. The mean value of the histogram is computed using the equivalent of Eq. (7). Our histogram agrees with the pixel-by-pixel analysis at intermediate and low Galactic latitudes presented in Planck Collaboration Int. XXXII (2016). Like in this earlier study, we find that the filaments are statistically aligned with BPOS. A similar alignment between the filaments in the intensity map and BPOS has been reported for synchrotron emission observed by WMAP at 23 GHz (Vidal et al. 2015).

To quantify the shape of the histogram of , we fit it with a Gaussian plus a constant. The Gaussian has a 1σ dispersion of 19°. The constant may be accounted for by the projection of the magnetic field and filament orientations on the POS as discussed in Planck Collaboration Int. XXXII (2016).

thumbnail Fig. 6

Upper panel: HRO between the filaments and BPOS. Middle panel: HRO between BPOS and Bm,POS. Lower panel: HRO between the filaments and Bm,POS.

thumbnail Fig. 7

Mean images of the Planck⟨ D353 ⟩, , and maps over the λ− filaments at 15′ resolution.

4.2. Relative orientation of the magnetic field and the large-scale Galactic magnetic field

Here, we compare the orientation of BPOS on the filaments with that of Bm,POS. Heiles (1996) derived the orientation of Bm,POS pointing towards and from the polarization pseudo-vectors of stars more distant than 500 parsecs. Slightly different l0 values have been reported in other studies. From the rotation measures of nearby pulsars within a few hundred parsecs of the Sun, Rand & Lyne (1994) found the direction of Bm,POS pointing towards l0 = 88° ± 5°. In another study of pulsar rotation measures, Han et al. (1999) derived the direction of Bm,POS as l0 ≃ 82°. These two studies do not report values for b0, which is assumed to be zero. Based on these observations, we assume that the mean orientation of Bm,POS in the solar neighbourhood is l0 = 84° ± 10° and b0 = 0° ± 10°, with the same uncertainty on l0 and b0.

We construct a pseudo-vector field with unit length based on the uniform orientation of Bm,POS. This pseudo-vector has components: Qm = cos 2ψm = cos 2(χm−π/ 2) and Um = −sin 2ψm = −sin 2(χm−π/ 2) (following the HEALPix convention for the Qm and Um maps), where ψm is the polarization angle of Bm,POS. The procedure to go from the uniform Bm,POS pointing towards (l0,b0) to ψm is detailed by Heiles (1996). The mean orientation angle () of Bm,POS for each filament is obtained by first averaging Qm and Um over all n pixels within a filament and then calculating the position angle of this averaged pseudo-vector. We compute the angle difference, , between the orientations of BPOS and Bm,POS on the filament in a similar manner to the method described in Sect. 4.1.

The HRO between BPOS and Bm,POS for our filament sample is presented in the middle panel of Fig. 6. The mean value of the histogram is where the uncertainty is computed by changing the mean orientation of Bm,POS within its quoted uncertainties. This HRO has a larger dispersion than that between the filaments and BPOS shown in the upper panel of Fig. 6. To quantify the shape of the histogram of , we fit it with a Gaussian plus a constant. The Gaussian has a 1σ dispersion of 36°.

We conclude that BPOS of our filament sample is statistically aligned with Bm,POS. Planck Collaboration Int. XXXII (2016) reports a similar correlation, for the low and intermediate Galactic latitudes, when comparing the polarization measured on the filaments with their background polarization maps. The scatter measured by the HRO may be interpreted considering both the turbulent component of the magnetic field and projection effects.

4.3. Relative orientation of the large-scale Galactic magnetic field and the filaments

We combine the results obtained in Sects. 4.1 and 4.2 to assess statistically the orientation of Bm,POS in the solar neighbourhood with respect to the filaments. BPOS is statistically aligned with the filaments in our sample and with Bm,POS. From both results, one would intuitively expect Bm,POS to be statistically aligned with the filaments. To test this expectation, for each filament, we compute the angle difference, , between the orientations of Bm,POS and the filament. The angle difference is computed in a similar manner to the method described in Sect. 4.1.

The HRO between the filament and Bm,POS for our filament sample is presented in the lower panel of Fig. 6. A correlation between the orientation angles of the filament and Bm,POS is present, but the HRO shows more scatter than the HRO between the filaments and BPOS and that between BPOS and Bm,POS. The histogram has a mean value of . To quantify the shape of the histogram of , we fit with a Gaussian plus a constant. The Gaussian has a 1σ dispersion of 54°. Planck Collaboration Int. XXXII (2016) reported a similar loss of correlation when comparing the orientations of the filaments with that of BPOS derived from their local background polarization maps.

5. Mean polarization properties of filaments

In this section, we present stacked images of the filaments in Stokes I, Q, and U, after rotation to align the filaments and to compute Q and U with respect to their orientation. The images are used to compute the average polarization fraction of our sample of filaments.

5.1. Stacking filaments

Over the high-latitude sky, the signal-to-noise ratio of the 353 GHz Planck polarization maps is low and it is not possible to measure the polarization fraction of individual dust intensity filaments in our sample. In order to increase the signal-to-noise ratio, we therefore stack images of the 259 filaments and their surroundings.

thumbnail Fig. 8

Radial profiles of the mean stacked PlanckD353, , and images as functions of distance from the centre of the filament (blue line). The grey shaded region shows the 1σ dispersion from the data values at a given radial distance from the filament axis that we average. The dashed line is the Gaussian fit to the filament profile plus a constant background emission.

For each filament in the sample, using the gnomview routine of HEALPix, we extract from the Planck maps a local, flat-sky, image (7° × 5° patch) centred on the filament centre and rotated by in the clockwise direction to align the filament in the North-South direction. We stack the images of the filaments in D353, , and (as defined in Eqs. (8) and (9)) after aligning all the maps in the North-South direction. We produce mean stacked images, denoted with angle brackets ⟨ ... ⟩, by dividing the sum of the individual images by the total number of filaments in our sample; they are presented in Fig. 7. The 1σ errorbar both on the and images is 1.3  μKCMB, as computed from the difference of two polarization HM maps. All the features presented in Fig. 7 are significant compared to the data systematics and statistical noise. The average filament appears as a negative feature with respect to the background in the image and is not seen in the image. This result is a direct consequence of the alignment between the filaments and BPOS (Sect. 4.1 and Zaldarriaga 2001). The background in both the and images is rather homogeneous. This reflects the smoothness of BPOS within the 7° × 5° patches.

We perform a null test to assess the significance of the stacking of filaments. This test is made by stacking 259 randomly chosen 7° × 5° patches in the high-latitude sky. Each patch is rotated in the clockwise direction, with the orientation angle θ− of the central pixel. The images of ⟨ D353 ⟩, , and for random patches are consistent with noise. The amplitude of and images is comparable to that of the difference between stacked images obtained when applying the same analysis to each of the two polarization HM maps. This confirms the hypothesis that the filaments detected in Fig. 7 are indeed real and are rotated with a well-determined angle .

5.2. Polarization fraction

Instead of using individual pixels, we collapse the mean stacked images in the filament direction to draw the radial profiles (R) of the D353, , and images, which are presented in Fig. 8. The shaded area in Fig. 8 represents the 1σ dispersion from the data values at a given radial distance from the filament axis that we average. We clearly identify the profile of the filament on top of the constant background emission in the D353 and radial profiles, while the radial profile of is consistent with a constant background emission. The radial profiles of D353, , and can be decomposed into the filament (F) and the background (B) contributions as We fit the radial profiles of D353 and in Fig. 8 with a Gaussian profile for the filament emission plus a constant for the background emission. We find that the centre of the Gaussian profile is zero and that their 1σ dispersion is 27′ ± 1′ for both the D353 and radial profiles.

Following Eqs. (12)and (13), we can express the average Stokes and for one given filament as where and are the average specific intensity and polarization intensity of the filament. The superscript F represents the contribution from the filament only. The polarization fraction () of a filament is defined by (21)From Fig. 4, we know that most filaments in our sample have comparable column densities and hence roughly the same . The mean stacked Stokes and values for all the filaments can be approximated as where is the mean polarization fraction of our filament sample. For the radial profiles of the filament emission, we have The angle difference between the filaments and BPOS (Eq. (14)) is used for . The histogram of is roughly symmetric around 0°, implying .

Similarly, the radial profiles from the background emission can be written as where is the average polarization fraction of the background emission and is the angle difference between a given filament and BPOS from its local background polarization.

The observed radial profile of the stacked image in Fig. 8 can be written as (28)where is the scaling parameter and is the offset of the linear fit between the and RD353 profiles. The linear fit between and RD353 is shown in Fig. 9. The best-fit parameter values from the linear fit are

In our estimate of , the noise bias is negligible, because the noise averages out in the stacking of the and maps. The HRO between the filament and BPOS (upper panel in Fig. 6) is used to compute , which is 0.48. Taking this factor into account, the mean polarization fraction of our filament sample is (31)We have computed using the two independent subsets of the Planck data (HM maps). For both data sets, we find the same mean value of 11%. This is expected, because the HRO between the filaments and BPOS is determined with high accuracy (Appendix D). It shows that the measurement error on is small. Obviously may be different for another set of filaments, because it depends on the angles of the filaments with respect to the plane of the sky.

thumbnail Fig. 9

Linear fit (dashed line) to the correlation between the mean radial profiles of the stacked and D353 images. The grey shaded region shows the 1σ dispersion from the data values at a given radial distance from the filament axis that we average.

thumbnail Fig. 10

Mean images of the Planck, , , , and maps over the λ− filaments.

thumbnail Fig. 11

Mean images of the Planck, , , , and maps over the λ− filaments.

The mean polarization fraction of the filaments determined from the filament sample is smaller than the maximum degree of polarization reported in Planck Collaboration Int. XIX (2015), which is pmax = 19.6%. This result suggests that there is some depolarization due to changes in the magnetic field orientation within the filaments (Planck Collaboration Int. XX 2015; Planck Collaboration Int. XXXIII 2016).

The offset b from the linear fit to the radial profiles of the and D353 images is positive. This means that (32)The distribution of for our filament sample is needed to compute . We make a simple approximation using the HRO between the filament and Bm,POS (lower panel in Fig. 6) as a proxy for . We then find for our filament sample. Combining Eqs. (29)and (32), we put a upper limit on the mean polarization fraction of the background emission as (33)This upper limit on depends on the mean orientation of BGMF and it ranges from 48% to 14% for the 10° uncertainty on l0 and b0 (Sect. 4.2). The profile of in the lower panel of Fig. 8 is roughly constant and positive. This result indicates that BPOS is smooth within the 7° × 5° patches and is positive (Eq. (27)), which follows if the distribution of is not symmetric with respect to 0°. We point out that the histogram in the lower panel in Fig. 6 may not be used to estimate the sign of because the uncertainty on l0 and b0 is too large.

6. E–B asymmetry

In this section, we quantify the E−B asymmetry of dust polarization, over the ℓ range 30 to 300, for our filament sample, and relate it to the relative orientation between filaments and BPOS. First, we present the stacked images of filtered maps that we use to quantify the E−B asymmetry. Second, we compute the contribution of the pixels used in the stacking to the variance in the high-latitude sky. Last, we present an analytical approximation that relates the HRO between the filaments and BPOS to the E−B asymmetry.

6.1. Stacked images of filtered maps

Here we use the Planck, , , , and maps to analyse the ℓ range, between 30 and 300, over which we identify the filaments in the dust intensity map. Our choice of ℓ range has a large overlap with the angular scales, 40 <ℓ< 600, where the E−B asymmetry has been measured with the power spectra of dust polarisation. We stack all the bandpass-filtered maps, rotated by the mean orientation angle of the filament, as described in Sect. 5.1. The and maps are the filtered Stokes Q353 and U353 maps computed with respect to the axis of the filament. The mean stacked images of the bandpass-filtered maps are presented in Fig. 10. The sidelobes that appear in the images of , , and of Fig. 10 on both sides of the filament centre are coming from the filtering. The 1σ errorbar on the , , , and images is 0.2  μKCMB, as computed from the difference of two polarization HM maps. The filaments appear as a negative feature in the image and a positive feature in that of . The and images are consistent with a mean value zero.

Next, we stack the products of two quantities, i.e., , , , , and . We follow the same methodology as discussed in Sect. 5.1 with oriented stacking at the centre of the filament. To avoid a noise bias in the square quantities, and , we compute the cross-product of the two HalfMission maps (HM1 and HM2). For other quantities, , , and , we use the full-mission maps. The images produced by stacking the products are presented in Fig. 11. Note that the maps shown in Figs. 10 and 11 are not used for data analysis, but still we comment on these images. We computed the mean stacked images from each of the two polarization DS or HR maps and all the features presented in Fig. 11 appear for both independent subset of Planck data. In the idealized description of the filaments in Zaldarriaga (2001), we will expect the average filament to appear similarly in the EE, BB, and TE maps. The differences between these three images shows that the reality of the dust sky is more complex than the idealized model.

6.2. Measured E–B asymmetry

We measure the E−B asymmetry at the filtering scale over the high-latitude (HL) region. We compute the variances (V) of the and maps using the relations where NHL is the total number of pixels in the HL region. The ratio of the filtered B353 and E353 variances is (36)The uncertainty on VEE is computed by repeating the calculation of Eq. (34)using the different cross-products, i.e., the two HalfRing (HR1 and HR2) and the two DetSet (DS1 and DS2) maps. We use the cross-HalfMissions as a reference for mean VEE. The 1σ uncertainty on VEE comes from the differences of these cross-products (DetSets minus HalfMissions and HalfRings minus HalfMissions). This 1σ uncertainty is dominated by data systematics rather than statistical noise. This is in agreement with the uncertainties on power spectra over the same ℓ range, as shown for fsky = 0.5 in Fig. 2 of Planck Collaboration Int. XXX (2016). The statistical noise on VEE is estimated from the cross-HalfMissions between their two HalfRing half-differences, as (37)A similar procedure is applied to compute the uncertainty on VBB. The ratio of the B353 and E353 variances differs slightly from the measurement at the power spectra level (Planck Collaboration Int. XXX 2016), probably because of the multipole range over which the ratio of the B353 and E353 variances are computed.

thumbnail Fig. 12

Map of the selected pixels (grey colour) used in the stacking analysis. It covers 28% of the high-latitude sky. Each tile in the image is a 7° × 5° patch around the filament centre and rotated by .

We compute the covariances of the three maps over the high-latitude sky using the relations The above covariances divided by VEE are listed in Table 1. The ratio VTE/VEE that we find is consistent with the measurements of Planck Collaboration Int. XXX (2016) at the power spectrum level.

6.3. Contribution of filaments to the variance of the E and B maps

In this section, we compute the variance of the and maps over the sky pixels used to produce the stacked images in Fig. 10, and compare the values with those measured over the HL region. These pixels are within the 7° × 5° patches, with an orientation angle , centred on the filaments (Sect. 5.1). These pixels define the grey regions in Fig. 12. We label them as SP and the rest of the high-latitude sky as O. The stacking procedure includes the filaments along with their surrounding background emission and, hence, effectively increases the selected fraction of the high-latitude sky, f1, from 2.2% (filament pixels as described in Sect. 3.2) to 28%.

Table 1

Ratios of the variances computed from the selected pixels (SP) used in the stacking analysis and the rest (O) of the high-latitude sky (HL).

We compute the variance from the SP pixels using the relation given in Eq. (34), (41)The sky variance of the map in the high-latitude sky can be written as the sum of contributions from SP (f1 = 0.28) and O (1−f1 = 0.72) regions. It is given by (42)The ratio (RSP) of the variance from the stacked pixels to the total sky variance is given by (43)The value of RSP is expected to be high, since the filaments are bright structures on the sky. The pixels we used for stacking contribute 83% of the total sky variance in the high-latitude sky. A similar result has been reported for the synchrotron emission, where bright filaments/shells are also measured to contribute most of the sky variance in polarization (Vidal et al. 2015).

It has been noted that the structure in the Planck353 GHz dust polarization maps is not fully accounted for by the filaments seen in the total dust intensity map. In particular, the local dispersion of the polarization angle shows structures in the polarization maps that have no counterpart in total intensity (Planck Collaboration Int. XIX 2015). These structures are thought to trace morphology of BPOS uncorrelated with matter structures (Planck Collaboration Int. XX 2015). However, as our RSP value shows, these polarization structures do not contribute much to the variance of the dust polarization.

In the same way as in Sect. 6.2, we compute the variances VEE, VBB, VTE, VTB, and VEB over the SP and O regions. Table 1 presents the ratios of the variances computed over different sky regions.

6.4. Analytical approximation

Following Zaldarriaga (2001) and the description of the filaments in Eqs. (19)and (20), we can express the E- and B-modes of a given filament as where is the mean polarization intensity of the filament. The direct relation between and , and and only holds for an idealised filament. For N idealized filaments oriented arbitrarily on the sky with respect to the GN, the ratio between the variances of the B and E maps is given by (46)which is expanded in terms of power spectra and under the assumption of statistical isotropy and homogeneity. The beam window function (wℓ) is the filter function. Both the observed and dust power spectra follow a power-law model with the same slope α (Planck Collaboration Int. XXX 2016). This reduces Eq. (46)to (47)From the histogram of the upper panel of Fig. 6, the distribution of the angle is known for our filament sample. Similar to the assumption made in Sect. 5.2, we assume that all the filaments have the same polarized intensity and therefore (48)We have computed the ratio ABB/AEE using the two independent subsets of the Planck data (HM maps) and find the same mean value of 0.66. This value of the ABB/AEE ratio based on this analytical model matches the observed mean value of 0.62 ± 0.03 (Sect. 6.2). We note that the model value is directly inferred from the distribution of for our filament sample. If the HRO of was flat with uniform probability between − 90° and + 90°, we would have found equal variances in both the E353 and B353 maps.

In summary, we propose that the alignment between BPOS and the filament orientations accounts for the E−B asymmetry in the range of angular scales 30 <ℓ< 300. The mean value of VBB/VEE for O region is consistent with the HL value within 1.9σ (Table 1). This shows that a similar alignment between the matter structures and BPOS can be inferred over the rest of the high-latitude sky. Some high Galactic latitude sky areas, such as the BICEP2 field (BICEP2 Collaboration 2014), do not include any of the strong filaments from our study. The Planck 353 GHz polarization maps do not have the required signal-to-noise ratio to measure the ABB/AEE ratio for individual BICEP2-like fields (Planck Collaboration Int. XXX 2016). Therefore more sensitive observations will be needed to test whether our interpretation is relevant there.

7. Relation to Galactic astrophysics

In this section we place the paper results in the context of earlier studies about the filamentary structure of interstellar matter and its correlation with the Galactic magnetic field.

Over the last decades, observations of interstellar gas and dust have been revealing the filamentary structure of the ISM in increasing details. Before Planck and Herschel, the discovery of the infrared cirrus with the IRAS and Hi all-sky surveys was a main milestone in our perception of the structure of the diffuse ISM (Boulanger 1994; Kalberla & Kerp 2009). At high Galactic latitude, the Planck dust emission is tightly correlated with Hi emission at local velocities (Planck Collaboration Int. XVII 2014). In particular, all of the filaments in our sample have an Hi counterpart. They are selected on the 353 GHz map but are seen at far-infrared wavelengths, in particular the IRAS 100  μm map.

The interstellar filaments seen in the Planck353 GHz dust intensity map are not all straight. With our filament-finding algorithm, we have identified the straight segments with lengths L ≥ 2°. Some of these segments are pieces of longer non-straight filaments. The 259 filaments in our sample make most of the sky variance in polarization as measured in our analysis. This is not a complete sample but other filaments at high Galactic latitude do not contribute much to the dust power in E-modes.

A number of studies, starting with the pioneering work of Goodman et al. (1990), have used the polarization of background starlight to investigate the relative orientation between the interstellar filaments and BPOS. While most studies have targeted filaments identified in extinction maps of molecular clouds in the solar neighbourhood, a few studies have focussed on filaments seen in Hi emission (McClure-Griffiths et al. 2006; Clark et al. 2014). These last two papers report a preferred alignment between the filaments and BPOS in the diffuse ISM. The analysis of Planck data has complemented earlier studies providing greater statistics and sensitivity. Planck Collaboration Int. XXXII (2016) compare the orientations of matter structures identified in the Planck353 GHz map with that of the Galactic magnetic field at intermediate latitudes. The alignment between the filaments and BPOS reported in this paper becomes weaker for increasing column density (see Fig. 15 of Planck Collaboration Int. XXXII 2016). Towards molecular clouds the relative orientation is observed to change progressively from preferentially parallel in areas with the lowest column density to preferentially perpendicular in the areas with the highest column density (Planck Collaboration Int. XXXV 2016). The transition occurs at a column density of 1021.7 cm-2. All the filaments considered in our analysis are much below this transition limit and observed, as expected from earlier studies, to be statistically aligned with BPOS.

Planck Collaboration Int. XXXII (2016) and Planck Collaboration Int. XXXV (2016) discuss these observational results in light of MHD simulations, which quantify the respective roles of the magnetic field, turbulence, and gas self-gravity in the formation of structures in the magnetized ISM. The alignment between the magnetic field and matter structures in the diffuse ISM is thought to be a signature of turbulence. Simulations show that turbulent flows will tend to stretch gas condensations into sheets and filaments, which appear elongated in column density maps (Hennebelle 2013). These structures will tend to be aligned with the magnetic field where the gas velocity is dynamically aligned with the field (Brandenburg & Lazarian 2013). Alignment also results from the fact that matter and the magnetic field are stretched in the same direction because the field is frozen into matter. The change in relative orientation observed within molecular clouds might be a signature of the formation of gravitationally bound structures in the presence of a dynamically important magnetic field. Indeed, Soler et al. (2013) report a change in the relative orientation between matter structures and the magnetic field, from parallel to perpendicular, for gravitationally bound structures in MHD simulations. This change is most significant for their simulation with the highest magnetization.

8. Conclusion

We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarization foreground to the CMB. The main results of our work are summarized as follows.

We filter the intensity and polarization maps to isolate structures over the range of angular scales where the E−B power asymmetry is observed. From a Hessian analysis of the Planck total dust intensity map at 353 GHz, we identify a sample of 259 filaments in the high-latitude sky with lengths L ≥ 2°. We measure the mean orientation angle of each filament in this sample and find that the filaments are statistically aligned with the plane of the sky component of the magnetic field, BPOS, inferred from the polarization angles measured by Planck. We also find that the orientation of BPOS is correlated with that of Bm,POS in the solar neighbourhood. Our results show that the correlation between the structures of interstellar matter and BPOS in the diffuse ISM reported in Planck Collaboration Int. XXXII (2016) for intermediate Galactic latitudes also applies to the lower column density filaments (a few 1019 cm-2) observed at high Galactic latitude.

We present mean images of our filament sample in dust intensity and Stokes Q353 and U353 with respect to the filament orientation ( and ), computed by stacking individual 7° × 5° patches centred on each filament. The stacked images show that the contribution of the filaments is a negative feature with respect to the background in the image and is not seen in the image. This result directly follows from the fact that the histogram of relative orientation between the filaments and BPOS peaks and is symmetric around 0°. Combining the stacked images and the histogram, we estimate the mean polarization fraction of the filaments to be 11%.

We relate the E−B asymmetry discovered in the power spectrum analysis of Planck 353 GHz polarization maps (Planck Collaboration Int. XXX 2016) to the alignment between the filaments and BPOS in the diffuse ISM. The set of 7° × 5° patches we stack represents 28% of the sky area at high Galactic latitude. The power of the E-mode dust polarization computed over this area amounts to 83% of the total dust polarization power in the high-latitude sky. We show with an analytical approximation of the filaments (based on the work of Zaldarriaga 2001), that the HRO between the filaments and BPOS may account for the ratio measured over the high-latitude sky. Our interpretation could also apply to the E−B asymmetry reported for the synchrotron emission (Planck Collaboration X 2015), since there is also a correlation between the orientation angle of BPOS and the filamentary structures of the synchrotron intensity map (Vidal et al. 2015; Planck Collaboration XXV 2015).

Present models of the dust polarization sky (e.g., O’Dea et al. 2012; Delabrouille et al. 2013) produce an equal amount of power in E- and B-modes for masks excluding the Galactic plane, because they ignore the correlation between the structure of the magnetic field and that of matter. Our work should motivate a quantitative modelling of the polarized sky, which will take into account the observed correlations between the Galactic magnetic field and the structure of interstellar matter.

1

Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states and led by Principal Investigators from France and Italy, telescope reflectors provided through a collaboration between ESA and a scientific consortium led and funded by Denmark, and additional contributions from NASA (USA).

Acknowledgments

The Planck Collaboration acknowledges the support of: ESA; CNES, and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, JA and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); ERC and PRACE (EU). A description of the Planck Collaboration and a list of its members, indicating which technical or scientific activities they have been involved in, can be found at http://www.cosmos.esa.int/web/planck/planck-collaboration. The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007−2013)/ERC grant agreement No. 267934. Some of the results in this paper have been derived using the HEALPix package.

References

  1. André, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, Protostars and Planets VI, 27 [Google Scholar]
  2. Arzoumanian, D., André, P., Didelon, P., et al. 2011, A&A, 529, L6 [Google Scholar]
  3. BICEP2 Collaboration 2014, Phys. Rev. Lett., 112, 241101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. BICEP2/Keck Array and Planck Collaborations 2015, Phys. Rev. Lett., 114, 101301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  5. Bond, N. A., Strauss, M. A., & Cen, R. 2010a, MNRAS, 406, 1609 [NASA ADS] [Google Scholar]
  6. Bond, N. A., Strauss, M. A., & Cen, R. 2010b, MNRAS, 409, 156 [NASA ADS] [CrossRef] [Google Scholar]
  7. Boulanger, F. 1994, in The First Symposium on the Infrared Cirrus and Diffuse Interstellar Clouds, eds. R. M. Cutri, & W. B. Latter, ASP Conf. Ser., 58, 101 [Google Scholar]
  8. Brandenburg, A., & Lazarian, A. 2013, Space Sci. Rev., 178, 163 [NASA ADS] [CrossRef] [Google Scholar]
  9. Clark, S. E., Peek, J. E. G., & Putman, M. E. 2014, ApJ, 789, 82 [NASA ADS] [CrossRef] [Google Scholar]
  10. Colombi, S., Pogosyan, D., & Souradeep, T. 2000, Phys. Rev. Lett., 85, 5515 [NASA ADS] [CrossRef] [Google Scholar]
  11. Delabrouille, J., Betoule, M., Melin, J.-B., et al. 2013, A&A, 553, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Forero-Romero, J. E., Hoffman, Y., Gottlöber, S., Klypin, A., & Yepes, G. 2009, MNRAS, 396, 1815 [NASA ADS] [CrossRef] [Google Scholar]
  13. Goodman, A. A., Bastien, P., Menard, F., & Myers, P. C. 1990, ApJ, 359, 363 [NASA ADS] [CrossRef] [Google Scholar]
  14. Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
  15. Hampel, F. R. 1974, J. Am. Stat. Assoc., 69, 383 [CrossRef] [Google Scholar]
  16. Han, J. L., Manchester, R. N., & Qiao, G. J. 1999, MNRAS, 306, 371 [NASA ADS] [CrossRef] [Google Scholar]
  17. Heiles, C. 1996, ApJ, 462, 316 [NASA ADS] [CrossRef] [Google Scholar]
  18. Hennebelle, P. 2013, A&A, 556, A153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Kalberla, P. M. W., & Kerp, J. 2009, ARA&A, 47, 27 [NASA ADS] [CrossRef] [Google Scholar]
  20. Komm, R. W., Gu, Y., Hill, F., Stark, P. B., & Fodor, I. K. 1999, ApJ, 519, 407 [NASA ADS] [CrossRef] [Google Scholar]
  21. Lamarre, J., Puget, J., Ade, P. A. R., et al. 2010, A&A, 520, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Lanusse, F., Rassat, A., & Starck, J.-L. 2012, A&A, 540, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. McClure-Griffiths, N. M., Dickey, J. M., Gaensler, B. M., Green, A. J., & Haverkorn, M. 2006, ApJ, 652, 1339 [NASA ADS] [CrossRef] [Google Scholar]
  24. Men’shchikov, A. 2013, A&A, 560, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Miville-Deschênes, M.-A., Martin, P. G., Abergel, A., et al. 2010, A&A, 518, L104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Monteserín, C., Barreiro, R. B., Sanz, J. L., & Martínez-González, E. 2005, MNRAS, 360, 9 [NASA ADS] [CrossRef] [Google Scholar]
  27. O’Dea, D. T., Clark, C. N., Contaldi, C. R., & MacTavish, C. J. 2012, MNRAS, 419, 1795 [NASA ADS] [CrossRef] [Google Scholar]
  28. Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Planck Collaboration XIII. 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Planck Collaboration I. 2015, A&A, submitted [arXiv:1502.01582] [Google Scholar]
  31. Planck Collaboration VII. 2016, A&A, in press, DOI: 10.1051/0004-6361/201525844 [Google Scholar]
  32. Planck Collaboration VIII. 2016, A&A, in press, DOI: 10.1051/0004-6361/201525820 [Google Scholar]
  33. Planck Collaboration X. 2015, A&A, submitted [arXiv:1502.01588] [Google Scholar]
  34. Planck Collaboration XXV. 2015, A&A, submitted [arXiv:1506.06660] [Google Scholar]
  35. Planck Collaboration Int. XVII. 2014, A&A, 566, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Planck Collaboration Int. XIX. 2015, A&A, 576, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Planck Collaboration Int. XX. 2015, A&A, 576, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Planck Collaboration Int. XXI. 2015, A&A, 576, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Planck Collaboration Int. XXII. 2015, A&A, 576, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Planck Collaboration Int. XXX. 2016, A&A, 586, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Planck Collaboration Int. XXXII. 2016, A&A, 586, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Planck Collaboration Int. XXXIII. 2016, A&A, 586, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Planck Collaboration Int. XXXV. 2016, A&A, 586, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Polychroni, D., Schisano, E., Elia, D., et al. 2013, ApJ, 777, L33 [NASA ADS] [CrossRef] [Google Scholar]
  45. Rand, R. J., & Lyne, A. G. 1994, MNRAS, 268, 497 [NASA ADS] [CrossRef] [Google Scholar]
  46. Soler, J. D., Hennebelle, P., Martin, P. G., et al. 2013, ApJ, 774, 128 [NASA ADS] [CrossRef] [Google Scholar]
  47. Sousbie, T. 2011, MNRAS, 414, 350 [NASA ADS] [CrossRef] [Google Scholar]
  48. Starck, J.-L., Moudden, Y., Abrial, P., & Nguyen, M. 2006, A&A, 446, 1191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Vidal, M., Dickinson, C., Davies, R. D., & Leahy, J. P. 2015, MNRAS, 452, 656 [NASA ADS] [CrossRef] [Google Scholar]
  50. Zaldarriaga, M. 2001, Phys. Rev. D, 64, 103001 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Hessian analysis

In this Appendix, we detail the implementation of the Hessian analysis to determine the orientations of the filaments. We start with the map at HEALPix resolution Nside = 512. For each pixel on the sky, we estimate the first and second derivatives of with respect to the Galactic longitude l and latitude b (as described in Monteserín et al. 2005). The Hessian matrix of is defined as (A.1)where

thumbnail Fig. A.1

Filtered realization of the Gaussian dust sky (upper panel) and the corresponding eigenvalue λ− (middle panel). The distribution of the eigenvalue λ− (lower panel) over the unmasked pixels is Gaussian.

We decompose the map into aℓm coefficients using the ianafast routine of HEALPix and then use the isynfast routine in HEALPix to compute the second partial derivatives at each pixel. This method is computationally faster than the one used in Planck Collaboration Int. XXXII (2016). The smallest eigenvalue, λ−, of the Hessian matrix is calculated as (A.5)where . The map of λ−, displayed in the lower left panel of Fig. 2, highlights the filaments in the map. The filament orientation angle, θ−, is calculated using the relation (A.6)This follows the IAU convention, since it is measured from the GN and positive to the East direction. The formula for the filament orientation angle θ− is equivalent to equation (9) of Planck Collaboration Int. XXXII (2016).

Appendix B: Gaussian realization of the dust sky

Here we apply the SMAFF algorithm on a Gaussian realization of the dust total intensity map using the same methodology as described in Sect. 3.2. The power spectrum is modelled as , where α is the slope of the power law. We choose α = −2.4 in the multipole range from ℓ = 30 to 300 (Planck Collaboration Int. XXII 2015).

We compute a Gaussian realization of a sky map on a HEALPix grid with an ℓ-2.4 power spectrum, and use the spline wavelet transform to filter it. Hereafter, we refer to the filtered simulated Gaussian map as the “Sb map”. We compute the smallest eigenvalue λ− of the Sb map using the Hessian analysis described in Appendix A. The Sb map and its corresponding λ− map are shown in the upper and middle panels of Fig. A.1.

We only consider the high-latitude sky. With the upper threshold K− = mλ−−3σλ−, we remove pixels as shown in the lower panel of Fig. A.1. Running our friend-of-friend algorithm on the map, we do not detect any filament with L ≥ 2°. For higher values of the threshold K−, we detect such filaments in the Sb map. We call them “weak” filaments. The threshold factor K−, used in our study, is a key factor to separate the strong filaments from the weak ones. The main results of this paper follow from the statistical properties of strong filaments. However, in Appendix C we demonstrate that they still hold for the weak filaments.

Appendix C: Effect of λ− thresholding on the filament count

We apply the threshold K− = mλ−−3σλ− in Sect. 3.2 to find strong filaments in the map. In this section, we change the threshold K− to quantify its effect on the HRO between the filaments and BPOS. We choose different thresholds and divide the selected λ− coherent structures into two categories, namely strong and weak filaments. The strong filaments are selected from the sky pixels that have λ−<mλ−−3σλ−. Whereas, the weak filaments are selected from the sky pixels that have mλ−−3σλ− ≤ λ−<mλ−−pσλ−, where p is a factor to defined the threshold K− = mλ−−pσλ−. In Table C.1, we list the number of weak filaments for different values of p.

By choosing a threshold K− = mλ−−1σλ− and running our friend-of-friend algorithm on the map, we double the number of filaments with L ≥ 2°. For this larger set, we compute the HRO between the filaments and BPOS, and present it in Fig. C.1. The 1σ dispersion of the angle difference remains the same as for our nominal set. This shows that the alignment between the structures of interstellar matter and BPOS holds even when including the weak filaments.

Including the weak filaments and their surrounding background emission in the analysis of Sect. 6.3 effectively increases the sky fraction from 28% to 50% of the high-latitude sky. With more sky coverage, the ratio of the variances, i.e., VBB/VEE, VTE/VEE, VTB/VEE, and VEB/VEE is close to the HL region values quoted in Table 1.

thumbnail Fig. C.1

Same as the upper panel of Fig. 6, when including angle differences from the weak filaments. The width of the histogram does not change.

Table C.1

Total filament count, including the strong and weak filaments, as a function of the threshold K− on the λ− map derived from the map.

Appendix D: Uncertainties on the angles and

In Sect. 4 we present our statistical analysis of the relative orientation between the filaments and BPOS. The filament orientation angle, , is a single number computed from the map. To estimate the error bar on , we propagate the noise in the D353 map (ND353) through the Hessian analysis. We simulate a Gaussian realization of the ND353 map and add it to the D353 map.

We filter the noise added to the D353 map in the same manner as the D353 map. For each filament in our sample, we compute the change of its orientation angle from the added noise using the formula (D.1)The distribution of the angle difference is presented with red colour in Fig. D.1. The mean of the distribution is 0° and the 1σ dispersion is 0.2°. This shows that the noise in the D353 map has little impact on our estimate of the filament orientation angle.

Next, we compute the uncertainty on the orientation angle of BPOS. We use the two HM Planck polarization maps, and compute the orientation angle of BPOS using Eq. (2). For each filament in our sample, we compute the difference between the two values obtained from the two HM maps, (D.2)The histogram of the angle difference is presented with blue colour in Fig. D.1. The mean of the distribution is 0° and the 1σ dispersion is 2°. The error on the orientation angle of BPOS is small compared to the width of the HROs in Fig. 6.

Last, we assess the impact of the map filtering on our analysis, computing the mean orientations of the filaments, over the same set of selected pixels, on the D353 map. We compute the angle difference of the filament orientation derived from the filtered and unfiltered Planck 353 GHz maps (D.3)The histogram of ΔΘ is plotted with green colour in Fig. D.1. The mean of the distribution is 0° and the 1σ dispersion is 7.0°. All the uncertainties measured in and are small compared to the dispersion of the HRO measured in Fig. 6. This means that the data noise (for the intensity and polarization maps) and the filtering of the data are not critical for our study.

thumbnail Fig. D.1

Histogram of changes in (red) and (blue) due to noise present in the Planck data. The histogram of changes in from the difference of the filtered and unfiltered D353 maps is shown in green.

All Tables

Table 1

Ratios of the variances computed from the selected pixels (SP) used in the stacking analysis and the rest (O) of the high-latitude sky (HL).

In the text
Table C.1

Total filament count, including the strong and weak filaments, as a function of the threshold K− on the λ− map derived from the map.

In the text

All Figures

thumbnail Fig. 1

The representative bandpass filter, retaining only the scales between ℓ = 30 and 300.
In the text
thumbnail Fig. 2

Data processing steps implemented to identify filaments from the Planck data. We start with the PlanckD353 map (upper left panel) smoothed at 15′ resolution. The bandpass-filtered map (upper right panel) is produced using the spline wavelet decomposition, retaining only the scales between ℓ = 30 and 300. The lower eigenvalue map of the Hessian matrix, λ−, is shown in the lower left panel. Structures identified in the high-latitude sky λ− map are shown in the lower right panel. The superimposed graticule is plotted in each image and labelled only on the lower right panel. It shows lines of constant longitude separated by 60° and lines of constant latitude separated by 30°. The same graticule is used in all plots of the paper.

In the text
thumbnail Fig. 3

Distribution of the eigenvalues λ− over the unmasked pixels in the high-latitude sky. The grey region represents the pixels that were used in the SMAFF algorithm to find strong filaments.
In the text
thumbnail Fig. 4

Histogram of the mean column density of the filament sample. The column density is computed from the map using the conversion factor derived in Planck Collaboration Int. XVII (2014).
In the text
thumbnail Fig. 5

Sketch of the mean orientation angle of the filament (), the magnetic field (), the polarization angle (), and the large-scale GMF () along the filament. All the angles are defined with respect to the GN and follow the IAU convention.
In the text
thumbnail Fig. 6

Upper panel: HRO between the filaments and BPOS. Middle panel: HRO between BPOS and Bm,POS. Lower panel: HRO between the filaments and Bm,POS.

In the text
thumbnail Fig. 7

Mean images of the Planck⟨ D353 ⟩, , and maps over the λ− filaments at 15′ resolution.

In the text
thumbnail Fig. 8

Radial profiles of the mean stacked PlanckD353, , and images as functions of distance from the centre of the filament (blue line). The grey shaded region shows the 1σ dispersion from the data values at a given radial distance from the filament axis that we average. The dashed line is the Gaussian fit to the filament profile plus a constant background emission.
In the text
thumbnail Fig. 9

Linear fit (dashed line) to the correlation between the mean radial profiles of the stacked and D353 images. The grey shaded region shows the 1σ dispersion from the data values at a given radial distance from the filament axis that we average.

In the text
thumbnail Fig. 10

Mean images of the Planck, , , , and maps over the λ− filaments.

In the text
thumbnail Fig. 11

Mean images of the Planck, , , , and maps over the λ− filaments.

In the text
thumbnail Fig. 12

Map of the selected pixels (grey colour) used in the stacking analysis. It covers 28% of the high-latitude sky. Each tile in the image is a 7° × 5° patch around the filament centre and rotated by .
In the text
thumbnail Fig. A.1

Filtered realization of the Gaussian dust sky (upper panel) and the corresponding eigenvalue λ− (middle panel). The distribution of the eigenvalue λ− (lower panel) over the unmasked pixels is Gaussian.

In the text
thumbnail Fig. C.1

Same as the upper panel of Fig. 6, when including angle differences from the weak filaments. The width of the histogram does not change.
In the text
thumbnail Fig. D.1

Histogram of changes in (red) and (blue) due to noise present in the Planck data. The histogram of changes in from the difference of the filtered and unfiltered D353 maps is shown in green.
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.